A phylogenomic approach to resolve the arthropod tree of life

Revisiting the four Hexapoda classes: Protura as the sister group to all other hexapods

Insects represent the most diverse animal group, yet previous phylogenetic analyses based on morphological and molecular data have failed to agree on the evolutionary relationships of early insects and their six-legged relatives (together constituting the clade Hexapoda). In particular, the phylogenetic positions of the three early-diverging hexapod lineages-the coneheads (Protura), springtails (Collembola), and two-pronged bristletails (Diplura)-have been debated for over a century, with alternative topologies implying drastically different scenarios of the evolution of the insect body plan and hexapod terrestrialization. We addressed this issue by sampling all hexapod orders and experimenting with a broad range of across-site compositional heterogeneous models designed to tackle ancient divergences. Our analyses support Protura as the earliest-diverging hexapod lineage ("Protura-sister") and Collembola as a sister group to Diplura, a clade corresponding to the original composition of Entognatha, and characterized by the shared possession of internal muscles in the antennal flagellum. The previously recognized 'Elliplura' hypothesis is recovered only under the site-homogeneous substitution models with partial supermatrices. Our cross-validation analysis shows that the site-heterogeneous CAT-GTR model, which recovers "Protura-sister," fits significantly better than homogeneous models. Furthermore, the morphologically unusual Protura are also supported as the earliest-diverging hexapod lineage by other lines of evidence, such as mitogenomes, comparative embryology, and sperm morphology, which produced results similar to those in this study. Our backbone phylogeny of hexapods will facilitate the exploration of the underpinnings of hexapod terrestrialization and megadiversity.

Orthoptera-specific target enrichment (OR-TE) probes resolve relationships over broad phylogenetic scales

Phylogenomic data are revolutionizing the field of insect phylogenetics. One of the most tenable and cost-effective methods of generating phylogenomic data is target enrichment, which has resulted in novel phylogenetic hypotheses and revealed new insights into insect evolution. Orthoptera is the most diverse insect order within polyneoptera and includes many evolutionarily and ecologically interesting species. Still, the order as a whole has lagged behind other major insect orders in terms of transitioning to phylogenomics. In this study, we developed an Orthoptera-specific target enrichment (OR-TE) probe set from 80 transcriptomes across Orthoptera. The probe set targets 1828 loci from genes exhibiting a wide range of evolutionary rates. The utility of this new probe set was validated by generating phylogenomic data from 36 orthopteran species that had not previously been subjected to phylogenomic studies. The OR-TE probe set captured an average of 1037 loci across the tested taxa, resolving relationships across broad phylogenetic scales. Our detailed documentation of the probe design and bioinformatics process is intended to facilitate the widespread adoption of this tool.

An updated world checklist of velvet worms (Onychophora) with notes on nomenclature and status of names

More than a decade has passed since the publication of the only world checklist available for Onychophora. During this period, numerous nomenclatural acts and taxonomic changes have been suggested within the group and a wealth of novel data has been published on many taxa. Herein, the up-to-date taxonomic scenario within Onychophora is presented, with appraisal of name status. This checklist covers both extant (Peripatidae and Peripatopsidae) and fossil taxa, and each species is accompanied by information on synonyms, type designation, holotype location, type locality, and language of original description. Additional remarks include nomenclatural inconsistencies, synonymizations, name misspellings, conflicting collecting event data, availability of taxonomically informative molecular data, etc. According to the data, 237 species are currently assigned to Onychophora: 140 of Peripatopsidae, 92 of Peripatidae, and five fossil species with unclear relationship to extant taxa. Since the previous checklist, 37 species have been added to Onychophora, representing an increase of 18.5% in the diversity described for the group. Yet, taxonomic descriptions seem slow-paced, with an average of 3.6 onychophoran species being described annually. From the taxonomic standpoint, 216 species are valid, although many of them require morphological revision and molecular characterization; 21 species exhibit major taxonomic ambiguities and have been regarded as nomina dubia. Recurrent taxonomic issues identified in the literature include inaccurate collecting event data, doubtful taxonomic assignment of molecular sequences, and non-observance of nomenclatural rules. These and other taxonomic aspects are addressed herein in the light of the directives established by the International Code of Zoological Nomenclature.

Exploring objective feature sets in constructing the evolution relationship of animal genome sequences

BACKGROUND

Exploring evolution regularities of genome sequences and constructing more objective species evolution relationships at the genomic level are high-profile topics. Based on the evolution mechanism of genome sequences proposed in our previous research, we found that only the 8-mers containing CG or TA dinucleotides correlate directly with the evolution of genome sequences, and the relative frequency rather than the actual frequency of these 8-mers is more suitable to characterize the evolution of genome sequences.

RESULT

Therefore, two types of feature sets were obtained, they are the relative frequency sets of CG1 + CG2 8-mers and TA1 + TA2 8-mers. The evolution relationships of mammals and reptiles were constructed by the relative frequency set of CG1 + CG2 8-mers, and two types of evolution relationships of insects were constructed by the relative frequency sets of CG1 + CG2 8-mers and TA1 + TA2 8-mers respectively. Through comparison and analysis, we found that evolution relationships are consistent with the known conclusions. According to the evolution mechanism, we considered that the evolution relationship constructed by CG1 + CG2 8-mers reflects the evolution state of genome sequences in current time, and the evolution relationship constructed by TA1 + TA2 8-mers reflects the evolution state in the early stage.

CONCLUSION

Our study provides objective feature sets in constructing evolution relationships at the genomic level.

Whole-Genome Phylogenetic Reconstruction as a Powerful Tool to Reveal Homoplasy and Ancient Rapid Radiation in Waterflea Evolution

Although phylogeny estimation is notoriously difficult in radiations that occurred several hundred million years ago, phylogenomic approaches offer new ways to examine relationships among ancient lineages and evaluate hypotheses that are key to evolutionary biology. Here, we reconstruct the deep-rooted relationships of one of the oldest living arthropod clades, the branchiopod crustaceans, using a kaleidoscopic approach. We use concatenation and coalescent tree-building methods to analyze a large multigene data set at the nucleotide and amino acid level and examine gene tree versus species tree discordance. We unequivocally resolve long-debated relationships among extant orders of the Cladocera, the waterfleas, an ecologically relevant zooplankton group in global aquatic and marine ecosystems that is famous for its model systems in ecology and evolution. To build the data set, we assembled eight de novo genomes of key taxa including representatives of all extant cladoceran orders and suborders. Our phylogenetic analysis focused on a BUSCO-based set of 823 conserved single-copy orthologs shared among 23 representative taxa spanning all living branchiopod orders, including 11 cladoceran families. Our analysis supports the monophyly of the Cladocera and reveals remarkable homoplasy in their body plans. We found large phylogenetic distances between lineages with similar ecological specializations, indicating independent evolution in major body plans, such as in the pelagic predatory orders Haplopoda and Onychopoda (the “Gymnomera”). In addition, we assessed rapid cladogenesis by estimating relative timings of divergence in major lineages using reliable fossil-calibrated priors on eight nodes in the branchiopod tree, suggesting a Paleozoic origin around 325 Ma for the cladoceran ancestor and an ancient rapid radiation around 252 Ma at the Perm/Triassic boundary. These findings raise new questions about the roles of homoplasy and rapid radiation in the diversification of the cladocerans and help examine trait evolution from a genomic perspective in a functionally well understood, ancient arthropod group. [Cladocera; Daphnia; evolution; homoplasy; molecular clock; phylogenomics; systematics; waterfleas.]

Evolutionary history and divergence times of Odonata (dragonflies and damselflies) revealed through transcriptomics

Dragonflies and damselflies are among the earliest flying insects with extant representatives. However, unraveling details of their long evolutionary history, such as egg laying (oviposition) strategies, is impeded by unresolved phylogenetic relationships, particularly in damselflies. Here we present a transcriptome-based phylogenetic reconstruction of Odonata, analyzing 2,980 protein-coding genes in 105 species representing nearly all the order's families. All damselfly and most dragonfly families are recovered as monophyletic. Our data suggest a sister relationship between dragonfly families of Gomphidae and Petaluridae. According to our divergence time estimates, both crown-Zygoptera and -Anisoptera arose during the late Triassic. Egg-laying with a reduced ovipositor apparently evolved in dragonflies during the late Jurassic/early Cretaceous. Lastly, we also test the impact of fossil choice and placement, particularly, of the extinct fossil species, †Triassolestodes asiaticus, and †Proterogomphus renateae on divergence time estimates. We find placement of †Proterogomphus renateae to be much more impactful than †Triassolestodes asiaticus.

Monophyletic blowflies revealed by phylogenomics

BACKGROUND

Blowflies are ubiquitous insects, often shiny and metallic, and the larvae of many species provide important ecosystem services (e.g., recycling carrion) and are used in forensics and debridement therapy. Yet, the taxon has repeatedly been recovered to be para- or polyphyletic, and the lack of a well-corroborated phylogeny has prevented a robust classification.

RESULTS

We here resolve the relationships between the different blowfly subclades by including all recognized subfamilies in a phylogenomic analysis using 2221 single-copy nuclear protein-coding genes of Diptera. Maximum likelihood (ML), maximum parsimony (MP), and coalescent-based phylogeny reconstructions all support the same relationships for the full data set. Based on this backbone phylogeny, blowflies are redefined as the most inclusive monophylum within the superfamily Oestroidea not containing Mesembrinellidae, Mystacinobiidae, Oestridae, Polleniidae, Sarcophagidae, Tachinidae, and Ulurumyiidae. The constituent subfamilies are re-classified as Ameniinae (including the Helicoboscinae, syn. nov.), Bengaliinae, Calliphorinae (including Aphyssurinae, syn. nov., Melanomyinae, syn. nov., and Toxotarsinae, syn. nov.), Chrysomyinae, Luciliinae, Phumosiinae, Rhiniinae stat. rev., and Rhinophorinae stat. rev. Metallic coloration in the adult is shown to be widespread but does not emerge as the most likely ground plan feature.

CONCLUSIONS

Our study provides the first phylogeny of oestroid calyptrates including all blowfly subfamilies. This allows settling a long-lasting controversy in Diptera by redefining blowflies as a well-supported monophylum, and blowfly classification is adjusted accordingly. The archetypical blowfly trait of carrion-feeding maggots most likely evolved twice, and the metallic color may not belong to the blowfly ground plan.

The evolution of insect biodiversity

Insects comprise over half of all described animal species. Together with the Protura (coneheads), Collembola (springtails) and Diplura (two-pronged bristletails), insects form the Hexapoda, a terrestrial arthropod lineage characterised by possessing six legs. Exponential growth of genome-scale data for the hexapods has substantially altered our understanding of the origin and evolution of insect biodiversity. Phylogenomics has provided a new framework for reconstructing insect evolutionary history, resolving their position among the arthropods and some long-standing internal controversies such as the placement of the termites, twisted-winged insects, lice and fleas. However, despite the greatly increased size of phylogenomic datasets, contentious relationships among key insect clades remain unresolved. Further advances in insect phylogeny cannot rely on increased depth and breadth of genome and taxon sequencing. Improved modelling of the substitution process is fundamental to countering tree-reconstruction artefacts, while gene content, modelling of duplications and deletions, and comparative morphology all provide complementary lines of evidence to test hypotheses emerging from the analysis of sequence data. Finally, the integration of molecular and morphological data is key to the incorporation of fossil species within insect phylogeny. The emerging integrated framework of insect evolution will help explain the origins of insect megadiversity in terms of the evolution of their body plan, species diversity and ecology. Future studies of insect phylogeny should build upon an experimental, hypothesis-driven approach where the robustness of hypotheses generated is tested against increasingly realistic evolutionary models as well as complementary sources of phylogenetic evidence.

A phylotranscriptomic framework for flesh fly evolution (Diptera, Calyptratae, Sarcophagidae)

The Sarcophagidae (flesh flies) comprise a large and widely distributed radiation within the Calyptratae (Diptera). Larval feeding habits are ecologically diverse and include sarcosaprophagy, coprophagy, herbivory, invertebrate and vertebrate predation, and kleptoparasitism. To elucidate the geographic origin and evolution of flesh fly life-history, we inferred a backbone phylogeny based on transcriptomic data from 26 sarcophagid species covering all three subfamilies plus 15 outgroups. The phylogeny was inferred using maximum parsimony and maximum likelihood methods based on a series of supermatrices, one set with overall information content improved by MARE (2290 loci), one set with 100% gene coverage for all included species (587 loci), and the last set including mitochondrial and nuclear genes (589 loci) and additional taxa. In order to obtain a more detailed hypothesis, we utilized the supertree approach to combine results from the present study with previously published hypotheses. This resulted supertree covers 84 of the one hundred currently recognized sarcophagid genera and formed the basis for the ancestral state reconstructions. The monophyletic Sarcophagidae is well-supported as sister to {Mystacinobiidae + Oestridae}, and relationships at the subfamily level are inferred as {Sarcophaginae, (Paramacronychiinae + Miltogramminae)}. The Sarcophagidae and each subfamily originated in the Americas, with Sarcophaginae diversifying mainly in the Neotropics, whereas the major radiation of both Miltogramminae and Paramacronychiinae occurred in the Palaearctic. Sarcosaprophagy is reconstructed as the ancestral larval feeding habit of the family Sarcophagidae and each subfamily. The ancestral sarcophagid larva probably utilized dead invertebrates as food, and the food spectrum expanded together with the diversification of breeding strategies. Particularly, kleptoparasitism in Miltogramminae is derived from sarcosaprophagy and may be seen as having derived from the breeding biology of 'lower' miltogrammines, the larvae of which feed on buried vertebrate carrion.

The Rapid Rise of Next-Generation Natural History

Many ecologists have lamented the demise of natural history and have attributed this decline to a misguided view that natural history is outdated and unscientific. Although there is a perception that the focus in ecology and conservation have shifted away from descriptive natural history research and training toward hypothetico-deductive research, we argue that natural history has entered a new phase that we call “next-generation natural history.” This renaissance of natural history is characterized by technological and statistical advances that aid in collecting detailed observations systematically over broad spatial and temporal extents. The technological advances that have increased exponentially in the last decade include electronic sensors such as camera-traps and acoustic recorders, aircraft- and satellite-based remote sensing, animal-borne biologgers, genetics and genomics methods, and community science programs. Advances in statistics and computation have aided in analyzing a growing quantity of observations to reveal patterns in nature. These robust next-generation natural history datasets have transformed the anecdotal perception of natural history observations into systematically collected observations that collectively constitute the foundation for hypothetico-deductive research and can be leveraged and applied to conservation and management. These advances are encouraging scientists to conduct and embrace detailed descriptions of nature that remain a critically important component of the scientific endeavor. Finally, these next-generation natural history observations are engaging scientists and non-scientists alike with new documentations of the wonders of nature. Thus, we celebrate next-generation natural history for encouraging people to experience nature directly.

A review of the hexapod tracheal system with a focus on the apterygote groups

Respiratory systems are key innovations for the radiation of terrestrial arthropods. It is therefore surprising that there is still a considerable lack of knowledge. In this review of the available information on tracheal systems of hexapods (with a focus on the apterygote lineages Protura, Collembola, Diplura, Archaeognatha and Zygentoma), we summarize available data on the spiracles (number, position and morphology), the shape and variability of tracheal branching patterns including anastomoses, the tracheal fine structure and the respiratory proteins. The available data are strongly fragmented, and information for most subgroups is missing. In various cases, individual observations for one species account for the knowledge of the entire order. The available data show that there are strong differences between but also within apterygote orders. We conclude that the available data are insufficient to derive detailed conclusions on the hexapod ground plan and outline the possible evolutionary scenarios for the tracheal system in this group.

Increasing 28 mitogenomes of Ephemeroptera, Odonata and Plecoptera support the Chiastomyaria hypothesis with three different outgroup combinations

BACKGROUND

The phylogenetic relationships of Odonata (dragonflies and damselflies) and Ephemeroptera (mayflies) remain unresolved. Different researchers have supported one of three hypotheses (Palaeoptera, Chiastomyaria or Metapterygota) based on data from different morphological characters and molecular markers, sometimes even re-assessing the same transcriptomes or mitochondrial genomes. The appropriate choice of outgroups and more taxon sampling is thought to eliminate artificial phylogenetic relationships and obtain an accurate phylogeny. Hence, in the current study, we sequenced 28 mt genomes from Ephemeroptera, Odonata and Plecoptera to further investigate phylogenetic relationships, the probability of each of the three hypotheses, and to examine mt gene arrangements in these species. We selected three different combinations of outgroups to analyze how outgroup choice affected the phylogenetic relationships of Odonata and Ephemeroptera.

METHODS

Mitochondrial genomes from 28 species of mayflies, dragonflies, damselflies and stoneflies were sequenced. We used Bayesian inference (BI) and Maximum likelihood (ML) analyses for each dataset to reconstruct an accurate phylogeny of these winged insect orders. The effect of outgroup choice was assessed by separate analyses using three outgroups combinations: (a) four bristletails and three silverfish as outgroups, (b) five bristletails and three silverfish as outgroups, or (c) five diplurans as outgroups.

RESULTS

Among these sequenced mitogenomes we found the gene arrangement IMQM in Heptageniidae (Ephemeroptera), and an inverted and translocated tRNA-Ile between the 12S RNA gene and the control region in Ephemerellidae (Ephemeroptera). The IMQM gene arrangement in Heptageniidae (Ephemeroptera) can be explained via the tandem-duplication and random loss model, and the transposition and inversion of tRNA-Ile genes in Ephemerellidae can be explained through the recombination and tandem duplication-random loss (TDRL) model. Our phylogenetic analysis strongly supported the Chiastomyaria hypothesis in three different outgroup combinations in BI analyses. The results also show that suitable outgroups are very important to determining phylogenetic relationships in the rapid evolution of insects especially among Ephemeroptera and Odonata. The mt genome is a suitable marker to investigate the phylogeny of inter-order and inter-family relationships of insects but outgroup choice is very important for deriving these relationships among winged insects. Hence, we must carefully choose the correct outgroup in order to discuss the relationships of Ephemeroptera and Odonata.

Host Phylogeny and Diet Shape Gut Microbial Communities Within Bamboo-Feeding Insects

The gut microbiome plays an important role in a host’s development and adaption to its dietary niche. In this study, a group of bamboo-feeding insects are used to explore the potential role of the gut microbiota in the convergent adaptation to extreme diet specialization. Specifically, using a 16S rRNA marker and an Illumina sequencing platform, we profiled the microbial communities of 76 gut samples collected from nine bamboo-feeding insects, including both hemimetabolous (Orthoptera and Hemiptera) and holometabolous (Coleoptera and Lepidoptera) species, which are specialized in three distinct dietary niches: bamboo leaf, shoot, and sap. The gut microbiota of these insects were dominated by Proteobacteria, Firmicutes, and Bacteroidetes and were clustered into solid (leaf and shoot) and liquid (sap) dietary niches. The gut bacterial communities of insects feeding on solid diet overlapped significantly, even though these insects belong to phylogenetically distant lineages representing different orders. In addition, the presence of cellulolytic bacterial communities within the gut microbiota allows bamboo-feeding insects to adapt to a highly specialized, fiber-rich diet. Although both phylogeny and diet can impact the structure and composition of gut microbiomes, phylogeny is the primary driving force underlying the convergent adaptation to a highly specialized diet, especially when the related insect species harbor similar gut microbiomes and share the same dietary niche over evolutionary timescales. These combined findings lay the foundation for future research on how convergent feeding strategies impact the interplays between hosts and their gut microbiomes and how the gut microbiota may facilitate convergent evolution in phylogenetically distant species in adaptation to the shared diet.

Resolving the Early Divergence Pattern of Teleost Fish Using Genome-Scale Data

Regarding the phylogenetic relationship of the three primary groups of teleost fishes, Osteoglossomorpha (bonytongues and others), Elopomorpha (eels and relatives), Clupeocephala (the remaining teleost fish), early morphological studies hypothesized the first divergence of Osteoglossomorpha, whereas the recent prevailing view is the first divergence of Elopomorpha. Molecular studies supported all the possible relationships of the three primary groups. This study analyzed genome-scale data from four previous studies: 1) 412 genes from 12 species, 2) 772 genes from 15 species, 3) 1,062 genes from 30 species, and 4) 491 UCE loci from 27 species. The effects of the species, loci, and models used on the constructed tree topologies were investigated. In the analyses of the data sets (1)–(3), although the first divergence of Clupeocephala that left the other two groups in a sister relationship was supported by concatenated sequences and gene trees of all the species and genes, the first divergence of Elopomorpha among the three groups was supported using species and/or genes with low divergence of sequence and amino-acid frequencies. This result corresponded to that of the UCE data set (4), whose sequence divergence was low, which supported the first divergence of Elopomorpha with high statistical significance. The increase in accuracy of the phylogenetic construction by using species and genes with low sequence divergence was predicted by a phylogenetic informativeness approach and confirmed by computer simulation. These results supported that Elopomorpha was the first basal group of teleost fish to have diverged, consistent with the prevailing view of recent morphological studies.

More than one way to smell ashore - Evolution of the olfactory pathway in terrestrial malacostracan crustaceans

Crustaceans provide a fascinating opportunity for studying adaptations to a terrestrial lifestyle because within this group, the conquest of land has occurred at least ten times convergently. The evolutionary transition from water to land demands various morphological and physiological adaptations of tissues and organs including the sensory and nervous system. In this review, we aim to compare the brain architecture between selected terrestrial and closely related marine representatives of the crustacean taxa Amphipoda, Isopoda, Brachyura, and Anomala with an emphasis on the elements of the olfactory pathway including receptor molecules. Our comparison of neuroanatomical structures between terrestrial members and their close aquatic relatives suggests that during the convergent evolution of terrestrial life-styles, the elements of the olfactory pathway were subject to different morphological transformations. In terrestrial anomalans (Coenobitidae), the elements of the primary olfactory pathway (antennules and olfactory lobes) are in general considerably enlarged whereas they are smaller in terrestrial brachyurans compared to their aquatic relatives. Studies on the repertoire of receptor molecules in Coenobitidae do not point to specific terrestrial adaptations but suggest that perireceptor events - processes in the receptor environment before the stimuli bind - may play an important role for aerial olfaction in this group. In terrestrial members of amphipods (Amphipoda: Talitridae) as well as of isopods (Isopoda: Oniscidea), however, the antennules and olfactory sensilla (aesthetascs) are largely reduced and miniaturized. Consequently, their primary olfactory processing centers are suggested to have been lost during the evolution of a life on land. Nevertheless, in terrestrial Peracarida, the (second) antennae as well as their associated tritocerebral processing structures are presumed to compensate for this loss or rather considerable reduction of the (deutocerebral) primary olfactory pathway. We conclude that after the evolutionary transition from water to land, it is not trivial for arthropods to establish aerial olfaction. If we consider insects as an ingroup of Crustacea, then the Coenobitidae and Insecta may be seen as the most successful crustacean representatives in this respect.

Global Analysis of Gene Expression Profiles Provides Novel Insights into the Development and Evolution of the Large Crustacean Eriocheir sinensis

Chinese mitten crab (Eriocheir sinensis) is an important aquaculture species in Crustacea. Functional analysis, although essential, has been hindered due to the lack of sufficient genomic or transcriptomic resources. In this study, transcriptome sequencing was conducted on 59 samples representing diverse developmental stages (fertilized eggs, zoea, megalopa, three sub-stages of larvae, juvenile crabs, and adult crabs) and different tissues (eyestalk, hepatopancreas, and muscle from juvenile crabs, and eyestalk, hepatopancreas, muscle, heart, stomach, gill, thoracic ganglia, intestine, ovary, and testis from adult crabs) of E. sinensis. A comprehensive reference transcriptome was assembled, including 19,023 protein-coding genes. Hierarchical clustering based on 128 differentially expressed cuticle-related genes revealed two distinct expression patterns during the early larval developmental stages, demonstrating the distinct roles of these genes in “crab-like” cuticle formation during metamorphosis and cuticle calcification after molting. Phylogenetic analysis of 1406 one-to-one orthologous gene families identified from seven arthropod species and Caenorhabditis elegans strongly supported the hypothesis that Malacostraca and Branchiopoda do not form a monophyletic group. Furthermore, Branchiopoda is more phylogenetically closely related to Hexapoda, and the clade of Hexapoda and Branchiopoda and the clade of Malacostraca belong to the Pancrustacea. This study offers a high-quality transcriptome resource for E. sinensis and demonstrates the evolutionary relationships of major arthropod groups. The differentially expressed genes identified in this study facilitate further investigation of the cuticle-related gene expression networks which are likely associated with “crab-like” cuticle formation during metamorphosis and cuticle calcification after molting.

Inferring Evolutionary Timescales without Independent Timing Information: An Assessment of "Universal" Insect Rates to Calibrate a Collembola (Hexapoda) Molecular Clock

Previous estimates of nucleotide substitution rates are routinely applied as secondary or “universal” molecular clock calibrations for estimating evolutionary timescales in groups that lack independent timing information. A major limitation of this approach is that rates can vary considerably among taxonomic groups, but the assumption of rate constancy is rarely evaluated prior to using secondary rate calibrations. Here I evaluate whether an insect mitochondrial DNA clock is appropriate for estimating timescales in Collembola—a group of insect-like arthropods characterized by high levels of cryptic diversity. Relative rates of substitution in cytochrome oxidase subunit 1 (COI) were inferred via Bayesian analysis across a topologically constrained Hexapod phylogeny using a relaxed molecular clock model. Rates for Collembola did not differ significantly from the average rate or from the rates estimated for most other groups (25 of 30), suggesting that (1) their apparent cryptic diversity cannot be explained by accelerated rates of molecular evolution and (2) clocks calibrated using “universal” insect rates may be appropriate for estimating evolutionary timescales in this group. However, of the 31 groups investigated, 10 had rates that deviated significantly from the average (6 higher, 4 lower), underscoring the need for caution and careful consideration when applying secondary insect rate calibrations. Lastly, this study exemplifies a relatively simple approach for evaluating rate constancy within a taxonomic group to determine whether the use of secondary rates are appropriate for molecular clock calibrations.

Optimizing the widely used nuclear protein-coding gene primers in beetle phylogenies and their application in the genus Sasajiscymnus Vandenberg (Coleoptera: Coccinellidae)

Advances in genomic biology and the increasing availability of genomic resources allow developing hundreds of nuclear protein-coding (NPC) markers, which can be used in phylogenetic research. However, for low taxonomic levels, it may be more practical to select a handful of suitable molecular loci for phylogenetic inference. Unfortunately, the presence of degenerate primers of NPC markers can be a major impediment, as the amplification success rate is low and they tend to amplify nontargeted regions. In this study, we optimized five NPC fragments widely used in beetle phylogenetics (i.e., two parts of carbamoyl-phosphate synthetase: CADXM and CADMC, Topoisomerase, Wingless and Pepck) by reducing the degenerate site of primers and the length of target genes slightly. These five NPC fragments and 6 other molecular loci were amplified to test the monophyly of the coccinellid genus Sasajiscymnus Vandenberg. The analysis of our molecular data set clearly supported the genus Sasajiscymnus may be monophyletic but confirmation with an extended sampling is required. A fossil-calibrated chronogram was generated by BEAST, indicating an origin of the genus at the end of the Cretaceous (77.87 Myr). Furthermore, a phylogenetic informativeness profile was generated to compare the phylogenetic properties of each gene more explicitly. The results showed that COI provides the strongest phylogenetic signal among all the genes, but Pepck, Topoisomerase, CADXM and CADMC are also relatively informative. Our results provide insight into the evolution of the genus Sasajiscymnus, and also enrich the molecular data resources for further study.

Evolutionary bedfellows: Reconstructing the ancestral state of autotomy and regeneration

Some form of regeneration occurs in all lifeforms and extends from single-cell organisms to humans. The degree to which regenerative ability is distributed across different taxa, however, is harder to ascertain given the potential for phylogenetic constraint or inertia, and adaptive processes to shape this pattern. Here, we examine the phylogenetic history of regeneration in two groups where the trait has been well-studied: arthropods and reptiles. Because autotomy is often present alongside regeneration in these groups, we performed ancestral state reconstructions for both traits to more precisely assess the timing of their origins and the degree to which these traits coevolve. Using an ancestral trait reconstruction, we find that autotomy and regeneration were present at the base of the arthropod and reptile trees. We also find that when autotomy is lost it does not re-evolve easily. Lastly, we find that the distribution of regeneration is intimately connected to autotomy with the association being stronger in reptiles than in arthropods. Although these patterns suggest that decoupling autotomy and regeneration at a broad phylogenetic scale may be difficult, the available data provides useful insight into their entanglement. Ultimately, our reconstructions provide the important groundwork to explore how selection may have played a role during the loss of regeneration in specific lineages.

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.

Comparison of transcriptomes from two chemosensory organs in four decapod crustaceans reveals hundreds of candidate chemoreceptor proteins

Crustaceans express genes for at least three classes of putative chemosensory proteins. These are: Ionotropic Receptors (IRs), derived from the heterotetrameric ionotropic glutamate receptors (iGluRs); Transient Receptor Potential (TRP) channels, a diverse set of sensor-channels that include several families of chemoreceptor channels; and Gustatory Receptor Like receptors (GRLs), ionotropic receptors that are homologues of Gustatory Receptors (GRs) of insects and are expressed sparingly in most crustaceans so far studied. IRs are typically numerically the most dominant of these receptor proteins in crustaceans and include two classes: co-receptor IRs, which are necessary for making a functional receptor-channel; and tuning IRs, whose specific combination in the IR subunits in the heterotetramer confers chemical specificity. Previous work showed that the transcriptomes from two major chemosensory organs-the lateral flagellum of the antennule (LF) and the tips of the legs (dactyls)-of the Caribbean spiny lobster Panulirus argus express four co-receptor IRs and over 100 tuning IRs. In this paper, we examined and compared the transcriptomes from the LF and dactyls of P. argus and three other decapod crustaceans-the clawed lobster Homarus americanus, red swamp crayfish Procambarus clarkii, and the blue crab Callinectes sapidus. Each species has at least ca. 100 to 250 IRs, 1 to 4 GRLs, and ca. 15 TRP channels including those shown to be involved in chemoreception in other species. The IRs show different degrees of phylogenetic conservation: some are arthropod-conserved, others are pancrustacean-conserved, others appear to be crustacean-conserved, and some appear to be species-specific. Many IRs appear to be more highly expressed in the LF than dactyl. Our results show that decapod crustaceans express an abundance of genes for chemoreceptor proteins of different types, phylogenetic conservation, and expression patterns. An understanding of their functional roles awaits determining their expression patterns in individual chemosensory neurons and the central projections of those neurons.

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.

Adult neurogenesis in the mushroom bodies of red flour beetles (Tribolium castaneum, HERBST) is influenced by the olfactory environment

Several studies showed adult persisting neurogenesis in insects, including the red flour beetle Tribolium castaneum, while it is absent in honeybees, carpenter ants, and vinegar flies. In our study, we focus on cell proliferation in the adult mushroom bodies of T. castaneum. We reliably labelled the progenies of the adult persisting mushroom body neuroblasts and determined the proliferation rate under several olfactory conditions within the first week after adult eclosion. We found at least two phases of Kenyon cell proliferation in the early adult beetle. Our results suggest that the generation of Kenyon cells during the first three days after adult eclosion is mainly genetically predetermined and a continuation of the developmental processes (nature), whereas from day four on proliferation seems to be mainly dependent on the odour environment (nurture). Considering that the mushroom bodies are linked to learning and memory, neurogenesis in the mushroom bodies is part of the remodelling of neuronal circuits leading to the adaption to the environment and optimization of behaviour.

Developmental gene expression as a phylogenetic data class: support for the monophyly of Arachnopulmonata

Despite application of genome-scale datasets, the phylogenetic placement of scorpions within arachnids remains contentious between two different phylogenetic data classes. Paleontologists continue to recover scorpions in a basally branching position, partly owing to their morphological similarity to extinct marine orders like Eurypterida (sea scorpions). Phylogenomic datasets consistently recover scorpions in a derived position, as the sister group of Tetrapulmonata (a clade of arachnids that includes spiders). To adjudicate between these hypotheses using a rare genomic change (RGC), we leveraged the recent discovery of ancient paralogy in spiders and scorpions to assess phylogenetic placement. We identified homologs of four transcription factors required for appendage patterning (dachshund, homothorax, extradenticle, and optomotor blind) in arthropods that are known to be duplicated in spiders. Using genomic resources for a spider, a scorpion, and a harvestman, we conducted gene tree analyses and assayed expression patterns of scorpion gene duplicates. Here we show that scorpions, like spiders, retain two copies of all four transcription factors, whereas arachnid orders like mites and harvestmen bear a single copy. A survey of embryonic expression patterns of the scorpion paralogs closely matches those of their spider counterparts, with one paralog consistently retaining the putatively ancestral pattern found in the harvestman, as well as the mite, and/or other outgroups. These data comprise a rare genomic change in chelicerate phylogeny supporting the inference of a distal placement of scorpions. Beyond demonstrating the diagnostic power of developmental genetic data as a phylogenetic data class, a derived placement of scorpions within the arachnids, together with an array of stem-group Paleozoic scorpions that occupied marine habitats, effectively rules out a scenario of a single colonization of terrestrial habitat within Chelicerata, even in tree topologies contrived to recover the monophyly of Arachnida.

The Genome of the Blind Soil-Dwelling and Ancestrally Wingless Dipluran Campodea augens: A Key Reference Hexapod for Studying the Emergence of Insect Innovations

The dipluran two-pronged bristletail Campodea augens is a blind ancestrally wingless hexapod with the remarkable capacity to regenerate lost body appendages such as its long antennae. As sister group to Insecta (sensu stricto), Diplura are key to understanding the early evolution of hexapods and the origin and evolution of insects. Here we report the 1.2-Gb draft genome of C. augens and results from comparative genomic analyses with other arthropods. In C. augens, we uncovered the largest chemosensory gene repertoire of ionotropic receptors in the animal kingdom, a massive expansion that might compensate for the loss of vision. We found a paucity of photoreceptor genes mirroring at the genomic level the secondary loss of an ancestral external photoreceptor organ. Expansions of detoxification and carbohydrate metabolism gene families might reflect adaptations for foraging behavior, and duplicated apoptotic genes might underlie its high regenerative potential. The C. augens genome represents one of the key references for studying the emergence of genomic innovations in insects, the most diverse animal group, and opens up novel opportunities to study the under-explored biology of diplurans.

Ordered phylogenomic subsampling enables diagnosis of systematic errors in the placement of the enigmatic arachnid order Palpigradi

The miniaturized arachnid order Palpigradi has ambiguous phylogenetic affinities owing to its odd combination of plesiomorphic and derived morphological traits. This lineage has never been sampled in phylogenomic datasets because of the small body size and fragility of most species, a sampling gap of immediate concern to recent disputes over arachnid monophyly. To redress this gap, we sampled a population of the cave-inhabiting species Eukoenenia spelaea from Slovakia and inferred its placement in the phylogeny of Chelicerata using dense phylogenomic matrices of up to 1450 loci, drawn from high-quality transcriptomic libraries and complete genomes. The complete matrix included exemplars of all extant orders of Chelicerata. Analyses of the complete matrix recovered palpigrades as the sister group of the long-branch order Parasitiformes (ticks) with high support. However, sequential deletion of long-branch taxa revealed that the position of palpigrades is prone to topological instability. Phylogenomic subsampling approaches that maximized taxon or dataset completeness recovered palpigrades as the sister group of camel spiders (Solifugae), with modest support. While this relationship is congruent with the location and architecture of the coxal glands, a long-forgotten character system that opens in the pedipalpal segments only in palpigrades and solifuges, we show that nodal support values in concatenated supermatrices can mask high levels of underlying topological conflict in the placement of the enigmatic Palpigradi.

Inside the head of a cybertype – three-dimensional reconstruction of the head muscles of Ommatoiulus avatar (Diplopoda: Juliformia: Julidae) reveals insights into the feeding movements of Juliformia

The origin and diversification of the arthropod head is one of the major topics in the field of evolutionary morphology of Arthropoda. Among the major arthropod groups, Myriapoda and, more precisely Diplopoda, are generally poorly studied regarding their head anatomy. However, this group is of pivotal importance to understand the evolutionary functional morphology of the arthropod head. In this study, we investigate the complete musculoskeletal system of the diplopod head with a detailed description of the cephalic anatomy of the recently described species Ommatoiulus avatar. The comparison of our data with the literature on the few other species available show that the morphology of the musculoskeletal system within Juliformia, a subgroup of the Diplopoda, is relatively conservative. Using video recordings of the feeding movements in addition to the anatomical data, we revise the mechanism of the mandibular movements in Juliformia. There was a controversy whether mandibular abduction is an active process, facilitated by contraction of an abductor muscle, or if it is a passive process, mediated by tentorial and gnathochilarial movements not involving a direct abduction by muscular contraction. We show that mandibular abduction in Ommatoiulus is an active movement involving the contraction of an abductor muscle. This is similar to the mandibular abduction in other arthropod groups.

Phylogenomics resolves major relationships and reveals significant diversification rate shifts in the evolution of silk moths and relatives

BACKGROUND

Silkmoths and their relatives constitute the ecologically and taxonomically diverse superfamily Bombycoidea, which includes some of the most charismatic species of Lepidoptera. Despite displaying spectacular forms and diverse ecological traits, relatively little attention has been given to understanding their evolution and drivers of their diversity. To begin to address this problem, we created a new Bombycoidea-specific Anchored Hybrid Enrichment (AHE) probe set and sampled up to 571 loci for 117 taxa across all major lineages of the Bombycoidea, with a newly developed DNA extraction protocol that allows Lepidoptera specimens to be readily sequenced from pinned natural history collections.

RESULTS

The well-supported tree was overall consistent with prior morphological and molecular studies, although some taxa were misplaced. The bombycid Arotros Schaus was formally transferred to Apatelodidae. We identified important evolutionary patterns (e.g., morphology, biogeography, and differences in speciation and extinction), and our analysis of diversification rates highlights the stark increases that exist within the Sphingidae (hawkmoths) and Saturniidae (wild silkmoths).

CONCLUSIONS

Our study establishes a backbone for future evolutionary, comparative, and taxonomic studies of Bombycoidea. We postulate that the rate shifts identified are due to the well-documented bat-moth "arms race". Our research highlights the flexibility of AHE to generate genomic data from a wide range of museum specimens, both age and preservation method, and will allow researchers to tap into the wealth of biological data residing in natural history collections around the globe.

The expression and characterization of recombinant cp19k barnacle cement protein from Pollicipes pollicipes

Adhesive proteins of barnacle cement have potential as environmentally friendly adhesives owing to their ability to adhere to various substrates in aqueous environments. By understanding the taxonomic breath of barnacles with different lifestyles, we may uncover commonalities in adhesives produced by these specialized organisms. The 19 kDa cement protein (cp19k) of the stalked barnacle Pollicipes pollicipes was expressed in Escherichia coli BL21 to investigate its adhesive properties. Initial expression of hexahistidine-tagged protein (rPpolcp19k-his) yielded low levels of insoluble protein. Co-overproduction of E. coli molecular chaperones GroEL-GroES and trigger factor (TF) increased soluble protein yields, although TF co-purified with the target protein (TF-rPpolcp19k-his). Surface coat analysis revealed high levels of adsorption of the TF-rPpolcp19k-his complex and of purified E. coli TF on both hydrophobic and hydrophilic surfaces, while low levels of adsorption were observed for rPpolcp19k-his. Tag-free rPpolcp19k protein also exhibited low adsorption compared to fibrinogen and Cell-Tak controls on hydrophobic, neutral hydrophilic and charged self-assembled monolayers under surface plasmon resonance assay conditions designed to mimic the barnacle cement gland or seawater. Because rPpolcp19k protein displays low adhesive capability, this protein is suggested to confer the ability to self-assemble into a plaque within the barnacle cement complex. This article is part of the theme issue 'Transdisciplinary approaches to the study of adhesion and adhesives in biological systems'.

Phylotranscriptomics confirms Alveopora is sister to Montipora within the family Acroporidae

The genus Alveopora is a scleractinian coral taxon whose phylogenetic classification has recently changed from the family Poritidae to Acroporidae. This change, which was made based on single-locus genetic data, has led to uncertainty about the placement of Alveopora and the ability for deep evolutionary relationships in these groups to be accurately recovered and represented by limited genetic datasets. We sought to characterize the higher-level position of Alveopora using newly available transcriptome data to confirm its placement within Acroporidae and resolve its closest ancestor. Here we present an analysis of a new 2031 gene dataset that confirms the placement of Alveopora within Acroporidae corroborating other single-locus (COI, 16S and ITS) analyses and a mitogenome dataset. We also resolve the position of Alveopora as sister to the genus Montipora. This has allowed the re-interpretation of morphology, and a rediagnosis of the family Acroporidae and the genus Alveopora.

Phylogenomics of the longitarsal Colossendeidae: The evolutionary history of an Antarctic sea spider radiation

Sea spiders (Pycnogonida) constitute a group of marine benthic arthropods that has a particularly high species diversity in the Southern Ocean. The "longitarsal" group of the sea spider family Colossendeidae is especially abundant in this region. However, this group also includes some representatives from other oceans, which raises the question where the group originates from. Therefore, we here investigated the phylogeny of the group with a hybrid enrichment approach that yielded a dataset of 1607 genes and over one million base pairs. We obtained a well-resolved phylogeny of the group, which is mostly consistent with morphological data. The data support an Antarctic origin of the longitarsal Colossendeidae and multiple dispersal events to other regions, which occurred at different timescales. This scenario is consistent with evidence found in other groups of marine invertebrates and highlights the role of the Southern Ocean as a source for non-Antarctic biota, especially of the deep sea. Our results suggest an initially slow rate of diversification followed by a more rapid radiation possibly correlated with the mid-Miocene cooling of Antarctica, similar to what is found in other taxa.

Going Deeper into High and Low Phylogenetic Relationships of Protura

Proturans are small, wingless, soil-dwelling arthropods, generally associated with the early diversification of Hexapoda. Their bizarre morphology, together with conflicting results of molecular studies, has nevertheless made their classification ambiguous. Furthermore, their limited dispersal capability (due to the primarily absence of wings) and their euedaphic lifestyle have greatly complicated species-level identification. Mitochondrial and nuclear markers have been applied herein to investigate and summarize proturan systematics at different hierarchical levels. Two new mitochondrial genomes are described and included in a phylum-level phylogenetic analysis, but the position of Protura could not be resolved with confidence due to an accelerated rate of substitution and extensive gene rearrangements. Mitochondrial and nuclear loci were also applied in order to revise the intra-class systematics, recovering three proturan orders and most of the families/subfamilies included as monophyletic, with the exception of the subfamily Acerentominae. At the species level, most morphologically described species were confirmed using molecular markers, with some exceptions, and the advantages of including nuclear, as well as mitochondrial, markers and morphology are discussed. At all levels, an enlarged taxon sampling and the integration of data from different sources may be of significant help in solving open questions that still persist on the evolutionary history of Protura.

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.

Comparative Analysis of the Adhesive Proteins of the Adult Stalked Goose Barnacle Pollicipes pollicipes (Cirripedia: Pedunculata)

Adhesion in barnacles is still poorly understood. The cement gland secretes an insoluble multi-protein complex, which adheres very strongly to a variety of substrates in the presence of water. This adhesion mechanism is bioinspiring for the engineering of new adhesive materials, but to replicate this adhesive system, the genes coding for the cement constitutive proteins must be identified and elucidated, and their products characterised. Here, the complete sequences of three cement protein (CP) genes (CP-100K, CP-52K, and CP-19K) isolated from the cement gland of the stalked barnacle Pollicipes pollicipes (order Scalpelliformes) were obtained using RACE PCR. The three genes were compared to the 23 other acorn barnacle CP genes so far sequenced (order Sessilia) to determine common and differential patterns and molecular properties, since the adhesives of both orders have visibly different characteristics. A shotgun proteomic analysis was performed on the cement, excreted at the membranous base of specimens, where the products of the three genes sequenced in the gland were identified, validating their function as CPs. A principal component analysis (PCA) was performed, to cluster CPs into groups with similar amino acid composition. This analysis uncovered three CP groups, each characterised by similar residue composition, features in secondary structure, and some biochemical properties, including isoelectric point and residue accessibility to solvents. The similarity among proteins in each defined group was low despite comparable amino acid composition. PCA can identify putative adhesive proteins from NGS transcriptomic data regardless of their low homology. This analysis did not highlight significant differences in residue composition between homologous acorn and stalked barnacle CPs. The characteristics responsible for the structural differences between the cement of stalked and acorn barnacles are described, and the presence of nanostructures, such as repetitive homologous domains and low complexity regions, and repetitive β-sheets are discussed relatively to self-assembly and adhesion.

Integration of phylogenomics and molecular modeling reveals lineage-specific diversification of toxins in scorpions

Scorpions have evolved a variety of toxins with a plethora of biological targets, but characterizing their evolution has been limited by the lack of a comprehensive phylogenetic hypothesis of scorpion relationships grounded in modern, genome-scale datasets. Disagreements over scorpion higher-level systematics have also incurred challenges to previous interpretations of venom families as ancestral or derived. To redress these gaps, we assessed the phylogenomic relationships of scorpions using the most comprehensive taxonomic sampling to date. We surveyed genomic resources for the incidence of calcins (a type of calcium channel toxin), which were previously known only from 16 scorpion species. Here, we show that calcins are diverse, but phylogenetically restricted only to parvorder Iurida, one of the two basal branches of scorpions. The other branch of scorpions, Buthida, bear the related LKTx toxins (absent in Iurida), but lack calcins entirely. Analysis of sequences and molecular models demonstrates remarkable phylogenetic inertia within both calcins and LKTx genes. These results provide the first synapomorphies (shared derived traits) for the recently redefined clades Buthida and Iurida, constituting the only known case of such traits defined from the morphology of molecules.

Tetraconatan phylogeny with special focus on Malacostraca and Branchiopoda: highlighting the strength of taxon-specific matrices in phylogenomics

Understanding the evolution of Tetraconata or Pancrustacea-the clade that includes crustaceans and insects-requires a well-resolved hypothesis regarding the relationships within and among its constituent taxa. Here, we assembled a taxon-rich phylogenomic dataset focusing on crustacean lineages based solely on genomes and new-generation Illumina-generated transcriptomes, including 89 representatives of Tetraconata. This constitutes, to our knowledge, the first phylogenomic study specifically addressing internal relationships of Malacostraca (with 26 species included) and Branchiopoda (36 species). Seven matrices comprising 81-684 orthogroups and 17 690-242 530 amino acid positions were assembled and analysed under five different analytical approaches. To maximize gene occupancy and to improve resolution, taxon-specific matrices were designed for Malacostraca and Branchiopoda. Key tetraconatan taxa (i.e. Oligostraca, Multicrustacea, Branchiopoda, Malacostraca, Thecostraca, Copepoda and Hexapoda) were monophyletic and well supported. Within Branchiopoda, Phyllopoda, Diplostraca, Cladoceromorpha and Cladocera were monophyletic. Within Malacostraca, the clades Eumalacostraca, Decapoda and Reptantia were well supported. Recovery of Caridoida or Peracarida was highly dependent on the analysis for the complete matrix, but it was consistently monophyletic in the malacostracan-specific matrices. From such examples, we demonstrate that taxon-specific matrices and particular evolutionary models and analytical methods, namely CAT-GTR and Dayhoff recoding, outperform other approaches in resolving certain recalcitrant nodes in phylogenomic analyses.

Multispecies coalescent analysis confirms standing phylogenetic instability in Hexapoda

The multispecies coalescent (MSC) has been increasingly used in phylogenomic analyses due to the accommodation of gene tree topological heterogeneity by taking into account population-level processes, such as incomplete lineage sorting. In this sense, the phylogeny of insect species, which are characterized by their large effective population sizes, is suitable for a coalescent-based analysis. Furthermore, studies so far recovered short internal branches at early divergences of the insect tree of life, indicating fast evolutionary radiations that increase the probability of incomplete lineage sorting in deep time. Here, we investigated the performance of the MSC for a phylogenomic data set of hexapods compiled by Misof et al. (2014, Science 346:763). Our analysis recovered the monophyly of most insect orders, and major phylogenetic relationships were in agreement with current insect systematics. We identified, however, some evolutionary associations that were consistently problematic. Most noticeable, Hexapod monophyly was disrupted by the sister group relationship between the remiped crustacean and Insecta. Additionally, the interordinal relationships within Polyneoptera and Neuropteroidea were found to be phylogenetically unstable. We show that these controversial phylogenetic arrangements were also poorly supported by previous analyses, and therefore, we evaluated their robustness to stochastic errors from sampling sites and terminals, confirming standing problems in hexapod phylogeny in the genomics age.

Fragmentary Gene Sequences Negatively Impact Gene Tree and Species Tree Reconstruction

Species tree reconstruction from genome-wide data is increasingly being attempted, in most cases using a two-step approach of first estimating individual gene trees and then summarizing them to obtain a species tree. The accuracy of this approach, which promises to account for gene tree discordance, depends on the quality of the inferred gene trees. At the same time, phylogenomic and phylotranscriptomic analyses typically use involved bioinformatics pipelines for data preparation. Errors and shortcomings resulting from these preprocessing steps may impact the species tree analyses at the other end of the pipeline. In this article, we first show that the presence of fragmentary data for some species in a gene alignment, as often seen on real data, can result in substantial deterioration of gene trees, and as a result, the species tree. We then investigate a simple filtering strategy where individual fragmentary sequences are removed from individual genes but the rest of the gene is retained. Both in simulations and by reanalyzing a large insect phylotranscriptomic data set, we show the effectiveness of this simple filtering strategy.

Reanalyzing the Palaeoptera problem - The origin of insect flight remains obscure

The phylogenetic relationships of the winged insect lineages - mayflies (Ephemeroptera), damselflies and dragonflies (Odonata), and all other winged insects (Neoptera) - are still controversial with three hypotheses supported by different datasets: Palaeoptera, Metapterygota and Chiastomyaria. Here, we reanalyze available phylogenomic data with a focus on detecting confounding and alternative signal. In this context, we provide a framework to quantitatively evaluate and assess incongruent molecular phylogenetic signal inherent in phylogenomic datasets. Despite overall support for the Palaeoptera hypothesis, we also found considerable signal for Chiastomyaria, which is not easily detectable by standardized tree inference approaches. Analyses of the accumulation of signal across gene partitions showed that signal accumulates gradually. However, even in case signal only slightly supported one over the other hypothesis, topologies inferred from large datasets switch from statistically strongly supported Palaeoptera to strongly supported Chiastomyaria. From a morphological point of view, Palaeoptera currently appears to be the best-supported hypothesis; however, recent analyses were restricted to head characters. Phylogenetic approaches covering all organ systems including analyses of potential functional or developmental convergence are still pending so that the Palaeoptera problem has to be considered an open question in insect systematics.

The mitochondrial genome of Caenis sp. (Ephemeroptera: Caenidae) and the phylogeny of Ephemeroptera in Pterygota

The phylogenetic relationship between Ephemeroptera (mayflies) and Odonata (dragonflies and damselflies) remains hotly debated in the insect evolution community. We sequenced the complete mitochondrial genome of Caenis sp. (Ephemeroptera: Caenidae) to discuss the phylogenetic relationship of Palaeoptera. The mitochondrial genome of Caenis sp. is a circular molecule of 15,254 bp in length containing 37 genes (13 protein-coding genes, 22 tRNAs, and 2 rRNAs), which showed the typical insect mitochondrial gene arrangement. In BI and ML phylogenetic trees using 71 species of 12 orders, our results support the Ephemeroptera as the basal group of winged insects.

Trajectories and Drivers of Genome Evolution in Surface-Associated Marine Phaeobacter

The extent of genome divergence and the evolutionary events leading to speciation of marine bacteria have mostly been studied for (locally) abundant, free-living groups. The genus Phaeobacter is found on different marine surfaces, seems to occupy geographically disjunct habitats, and is involved in different biotic interactions, and was therefore targeted in the present study. The analysis of the chromosomes of 32 closely related but geographically spread Phaeobacter strains revealed an exceptionally large, highly syntenic core genome. The flexible gene pool is constantly but slightly expanding across all Phaeobacter lineages. The horizontally transferred genes mostly originated from bacteria of the Roseobacter group and horizontal transfer most likely was mediated by gene transfer agents. No evidence for geographic isolation and habitat specificity of the different phylogenomic Phaeobacter clades was detected based on the sources of isolation. In contrast, the functional gene repertoire and physiological traits of different phylogenomic Phaeobacter clades were sufficiently distinct to suggest an adaptation to an associated lifestyle with algae, to additional nutrient sources, or toxic heavy metals. Our study reveals that the evolutionary trajectories of surface-associated marine bacteria can differ significantly from free-living marine bacteria or marine generalists.

Phylogeny of a cosmopolitan family of morphologically conserved trapdoor spiders (Mygalomorphae, Ctenizidae) using Anchored Hybrid Enrichment, with a description of the family, Halonoproctidae Pocock 1901

The mygalomorph family Ctenizidae has a world-wide distribution and currently contains nine genera and 135 species. However, the monophyly of this group has long been questioned on both morphological and molecular grounds. Here, we use Anchored Hybrid Enrichment (AHE) to gather hundreds of loci from across the genome for reconstructing the phylogenetic relationships among the nine genera and test the monophyly of the family. We also reconstruct the possible ancestral ranges of the most inclusive clade recovered. Using AHE, we generate a supermatrix of 565 loci and 115,209 bp for 27 individuals. For the first time, analyses using all nine genera produce results definitively establishing the non-monophyly of Ctenizidae. A lineage formed exclusively by representatives of South African Stasimopus was placed as the sister group to the remaining taxa in the tree, and the Mediterranean Cteniza and Cyrtocarenum were recovered with high support as sister to exemplars of Euctenizidae, Migidae, and Idiopidae. All the remaining genera-Bothriocyrtum, Conothele, Cyclocosmia, Hebestatis, Latouchia, and Ummidia-share a common ancestor. Based on these results, we formally elevate this clade to the level of family. Our results definitively establish both the non-monophyly of the Ctenizidae and non-validity of the subfamilies Ummidiinae and Ctenizinae. In order to establish the placement of the remaining three ctenizid genera, Cteniza, Cyrtocarenum, and Stasimopus, thorough analyses within the context of a complete mygalomorph phylogenetic framework are needed. We formally describe the family Halonoproctidae Pocock 1901 and infer that the family's most recent common ancestor was likely distributed in western North America and Asia.