The first myriapod genome sequence reveals conservative arthropod gene content and genome organisation in the centipede Strigamia maritima

Hepatincolaceae (Alphaproteobacteria) are Distinct From Holosporales and Independently Evolved to Associate With Ecdysozoa

The Hepatincolaceae (Alphaproteobacteria) are a group of bacteria that inhabit the gut of arthropods and other ecdysozoans, associating extracellularly with microvilli. Previous phylogenetic studies, primarily single-gene analyses, suggested their relationship to the Holosporales, which includes intracellular bacteria in protist hosts. However, the genomics of Hepatincolaceae is still in its early stages. In this study, the number of available Hepatincolaceae genomes was increased to examine their evolutionary and functional characteristics. It was found that the previous phylogenetic grouping with Holosporales was incorrect due to sequence compositional biases and that Hepatincolaceae form an independent branch within the Hepatincolaceae. This led to a reinterpretation of their features, proposing a new evolutionary scenario that involves an independent adaptation to host association compared to the Holosporales, with distinct specificities. The Hepatincolaceae exhibit greater nutritional flexibility, utilising various molecules available in the host gut and thriving in anaerobic conditions. However, they have a less complex mechanism for modulating host interactions, which are likely less direct than those of intracellular bacteria. In addition, representatives of Hepatincolaceae show several lineage-specific traits related to differences in host species and life conditions.

Functional characterization of CCHamides and deorphanization of their receptors in the yellow fever mosquito, Aedes aegypti

As a widely distributed anthropophilic mosquito species and vector of various arboviruses, Aedes aegypti poses a significant threat to human health on a global scale. Investigating mosquito neuropeptides allows us to better understand their physiology. The neuropeptides CCHamide1 (CCHa1) and CCHamide2 (CCHa2) along with their associated G protein-coupled receptors (CCHa1R and CCHa2R) were recently identified and studied across insects. However, expression profiles and physiological roles of CCHamides and their receptors in many other insects, including A. aegypti, remain unclear. This research aimed to quantify and localize the expression of CCHamides along with their receptors and gain insight on their physiological function in the yellow fever mosquito. RT-qPCR analysis revealed transcript abundance of CCHamides and receptors changes over development. Differential expression was also observed in tissues/organs of adult mosquitoes indicating CCHa1 and CCHa2 transcripts are enriched in the midgut, while receptors are expressed across various tissues. CCHamide immunoreactivity was observed in neurons in the brain and ventral nerve cord along with enteroendocrine cells in the posterior midgut adjacent to the midgut-hindgut junction, corroborating their transcript expression profiles. Using different mass spectrometrical approaches, presence of CCHamides were confirmed in the brain of both sexes, including the pars intercerebralis of female mosquitoes, as well as in the gut of adult mosquitoes. For chemical identification of predicted CCHamides, we analyzed brain and gut extracts by ESI-Q Exactive Orbitrap MS and resulting fragmentations confirmed CCHa1 and CCHa2 in brain and midgut samples of both male and female mosquitoes. A heterologous functional assay was used to confirm the specificity and sensitivity of the two CCHamide receptors by assessing their activation in response to diverse mosquito peptides, which confirmed CCHa1 and CCHa2 as natural ligands. Finally, using a capillary feeder (CAFE) bioassay, our results suggest that CCHa2 modulates feeding behaviour in female mosquitoes.

The moulting arthropod: a complete genetic toolkit review

Exoskeletons are a defining character of all arthropods that provide physical support for their segmented bodies and appendages as well as protection from the environment and predation. This ubiquitous yet evolutionarily variable feature has been instrumental in facilitating the adoption of a variety of lifestyles and the exploitation of ecological niches across all environments. Throughout the radiation that produced the more than one million described modern species, adaptability afforded by segmentation and exoskeletons has led to a diversity that is unrivalled amongst animals. However, because of the limited extensibility of exoskeleton chitin and cuticle components, they must be periodically shed and replaced with new larger ones, notably to accommodate the growing individuals encased within. Therefore, arthropods grow discontinuously by undergoing periodic moulting events, which follow a series of steps from the preparatory pre-moult phase to ecdysis itself and post-moult maturation of new exoskeletons. Each event represents a particularly vulnerable period in an arthropod's life cycle, so processes must be tightly regulated and meticulously executed to ensure successful transitions for normal growth and development. Decades of research in representative arthropods provide a foundation of understanding of the mechanisms involved. Building on this, studies continue to develop and test hypotheses on the presence and function of molecular components, including neuropeptides, hormones, and receptors, as well as the so-called early, late, and fate genes, across arthropod diversity. Here, we review the literature to develop a comprehensive overview of the status of accumulated knowledge of the genetic toolkit governing arthropod moulting. From biosynthesis and regulation of ecdysteroid and sesquiterpenoid hormones, to factors involved in hormonal stimulation responses and exoskeleton remodelling, we identify commonalities and differences, as well as highlighting major knowledge gaps, across arthropod groups. We examine the available evidence supporting current models of how components operate together to prepare for, execute, and recover from ecdysis, comparing reports from Chelicerata, Myriapoda, Crustacea, and Hexapoda. Evidence is generally highly taxonomically imbalanced, with most reports based on insect study systems. Biases are also evident in research on different moulting phases and processes, with the early triggers and late effectors generally being the least well explored. Our synthesis contrasts knowledge based on reported observations with reasonably plausible assumptions given current taxonomic sampling, and exposes weak assumptions or major gaps that need addressing. Encouragingly, advances in genomics are driving a diversification of tractable study systems by facilitating the cataloguing of putative genetic toolkits in previously under-explored taxa. Analysis of genome and transcriptome data supported by experimental investigations have validated the presence of an "ultra-conserved" core of arthropod genes involved in moulting processes. The molecular machinery has likely evolved with elaborations on this conserved pathway backbone, but more taxonomic exploration is needed to characterise lineage-specific changes and novelties. Furthermore, linking these to transformative innovations in moulting processes across Arthropoda remains hampered by knowledge gaps and hypotheses based on untested assumptions. Promisingly however, emerging from the synthesis is a framework that highlights research avenues from the underlying genetics to the dynamic molecular biology through to the complex physiology of moulting.

Chromosome-level genome of Scolopendra mutilans provides insights into its evolution

Report the first chromosome level genome of myriapod Scolopendra mutilans. Reveal gene expansions for importance to adapt. Annotate nine Hox cluster genes in this genome.

Irreducible Complexity of Hox Gene: Path to the Canonical Function of the Hox Cluster

The evolution of major taxa is often associated with the emergence of new gene families. In all multicellular animals except sponges and comb jellies, the genomes contain Hox genes, which are crucial regulators of development. The canonical function of Hox genes involves colinear patterning of body parts in bilateral animals. This general function is implemented through complex, precisely coordinated mechanisms, not all of which are evolutionarily conserved and fully understood. We suggest that the emergence of this regulatory complexity was preceded by a stage of cooperation between more ancient morphogenetic programs or their individual elements. Footprints of these programs may be present in modern animals to execute non-canonical Hox functions. Non-canonical functions of Hox genes are involved in maintaining terminal nerve cell specificity, autophagy, oogenesis, pre-gastrulation embryogenesis, vertical signaling, and a number of general biological processes. These functions are realized by the basic properties of homeodomain protein and could have triggered the evolution of ParaHoxozoa and Nephrozoa subsequently. Some of these non-canonical Hox functions are discussed in our review.

New genes helped acorn barnacles adapt to a sessile lifestyle

Barnacles are the only sessile lineages among crustaceans, and their sessile life begins with the settlement of swimming larvae (cyprids) and the formation of protective shells. These processes are crucial for adaptation to a sessile lifestyle, but the underlying molecular mechanisms remain poorly understood. While investigating these mechanisms in the acorn barnacle, Amphibalanus amphitrite, we discovered a new gene, bcs-6, which is involved in the energy metabolism of cyprid settlement and originated from a transposon by acquiring the promoter and cis-regulatory element. Unlike mollusks, the barnacle shell comprises alternate layers of chitin and calcite and requires another new gene, bsf, which generates silk-like fibers that efficiently bind chitin and aggregate calcite in the aquatic environment. Our findings highlight the importance of exploring new genes in unique adaptative scenarios, and the results will provide important insights into gene origin and material development.

Molecular Characterization, Expression and In Situ Hybridization Analysis of a Pedal Peptide/Orcokinin-type Neuropeptide in Cuttlefish Sepiella japonica

BACKGROUND

Neuropeptide pedal peptide (PP) and orcokinin (OK), which are structurally related active peptides, have been widely discovered in invertebrates and constitute the PP/OK neuropeptide family. They have complex structures and play myriad roles in physiological processes. To date, there have been no related reports of PP/OK-type neuropeptide in cephalopods, which possess a highly differentiated multi-lobular brain.

METHODS

Rapid Amplification of cDNA Ends (RACE) was employed to obtain the open reading frame (ORF) of PP/OK-type neuropeptide in Sepiella japonica (termed as Sj-PP/OK). Various software were used for sequence analysis. Semi-quantitative PCR was applied to analyze the tissue distribution profile, quantitative real-time PCR (qRT-PCR) was used to study spatio-temporal expression throughout the entire growth and development period, and in situ hybridization (ISH) was employed to observe the tissue location of Sj-PP/OK.

RESULTS

in the present study, we identified the ORF of Sj-PP/OK. The putative precursor of Sj-PP/ OK encodes 22 mature peptides, of which only tridecapeptides could undergo post-translationally amidated at C-terminus. Each of these tridecapeptides possesses the most conserved and frequent N-terminus Asp-Ser-Ile (DSI). Sequence analysis revealed that Sj-PP/OK shared comparatively low identity with other invertebrates PP or OK. The tissue distribution profile showed differences in the expression level of Sj-PP/OK between male and female. qRT-PCR data demonstrated that Sj-PP/OK was widely distributed in various tissues, with its expression level increasing continuously in the brain, optic lobe, liver, and nidamental gland throughout the entire growth and development stages until gonad maturation. ISH detected that Sj-PP/OK positive signals existed in almost all regions of the optic lobe except the plexiform zone, the outer edge of all functional lobes in the brain, epithelial cells and the outer membrane layer of the accessory nidamental gland. These findings suggest that Sj-PP/OK might play a role in the regulation of reproduction, such as vitellogenin synthesis, restoration, and ova encapsulation.

CONCLUSION

The study indicated that Sj-PP/OK may be involved in the neuroendocrine regulation in cephalopods, providing primary theoretical basis for further studies of its regulation role in reproduction.

Fast cycling culture of the annelid model Platynereis dumerilii

Platynereis dumerilii, a marine annelid, is a model animal that has gained popularity in various fields such as developmental biology, biological rhythms, nervous system organization and physiology, behaviour, reproductive biology, and epigenetic regulation. The transparency of P. dumerilii tissues at all developmental stages makes it easy to perform live microscopic imaging of all cell types. In addition, the slow-evolving genome of P. dumerilii and its phylogenetic position as a representative of the vast branch of Lophotrochozoans add to its evolutionary significance. Although P. dumerilii is amenable to transgenesis and CRISPR-Cas9 knockouts, its relatively long and indefinite life cycle, as well as its semelparous reproduction have been hindrances to its adoption as a reverse genetics model. To overcome this limitation, an adapted culturing method has been developed allowing much faster life cycling, with median reproductive age at 13-14 weeks instead of 25-35 weeks using the traditional protocol. A low worm density in boxes and a strictly controlled feeding regime are important factors for the rapid growth and health of the worms. This culture method has several advantages, such as being much more compact, not requiring air bubbling or an artificial moonlight regime for synchronized sexual maturation and necessitating only limited water change. A full protocol for worm care and handling is provided.

Evolution of the Spider Homeobox Gene Repertoire by Tandem and Whole Genome Duplication

Gene duplication generates new genetic material that can contribute to the evolution of gene regulatory networks and phenotypes. Duplicated genes can undergo subfunctionalization to partition ancestral functions and/or neofunctionalization to assume a new function. We previously found there had been a whole genome duplication (WGD) in an ancestor of arachnopulmonates, the lineage including spiders and scorpions but excluding other arachnids like mites, ticks, and harvestmen. This WGD was evidenced by many duplicated homeobox genes, including two Hox clusters, in spiders. However, it was unclear which homeobox paralogues originated by WGD versus smaller-scale events such as tandem duplications. Understanding this is a key to determining the contribution of the WGD to arachnopulmonate genome evolution. Here we characterized the distribution of duplicated homeobox genes across eight chromosome-level spider genomes. We found that most duplicated homeobox genes in spiders are consistent with an origin by WGD. We also found two copies of conserved homeobox gene clusters, including the Hox, NK, HRO, Irx, and SINE clusters, in all eight species. Consistently, we observed one copy of each cluster was degenerated in terms of gene content and organization while the other remained more intact. Focussing on the NK cluster, we found evidence for regulatory subfunctionalization between the duplicated NK genes in the spider Parasteatoda tepidariorum compared to their single-copy orthologues in the harvestman Phalangium opilio. Our study provides new insights into the relative contributions of multiple modes of duplication to the homeobox gene repertoire during the evolution of spiders and the function of NK genes.

Confluence and convergence of Dscam and Pcdh cell-recognition codes

The ability of neurites of the same neuron to avoid each other (self-avoidance) is a conserved feature in both invertebrates and vertebrates. The key to self-avoidance is the generation of a unique subset of cell-surface proteins in individual neurons engaging in isoform-specific homophilic interactions that drive neurite repulsion rather than adhesion. Among these cell-surface proteins are fly Dscam1 and vertebrate clustered protocadherins (cPcdhs), as well as the recently characterized shortened Dscam (sDscam) in the Chelicerata. Herein, we review recent advances in our understanding of how cPcdh, Dscam, and sDscam cell-surface recognition codes are expressed and translated into cellular functions essential for neural wiring.

The genome of the deep-sea anemone Actinernus sp. contains a mega-array of ANTP-class homeobox genes

Members of the phylum Cnidaria include sea anemones, corals and jellyfish, and have successfully colonized both marine and freshwater habitats throughout the world. The understanding of how cnidarians adapt to extreme environments such as the dark, high-pressure deep-sea habitat has been hindered by the lack of genomic information. Here, we report the first chromosome-level deep-sea cnidarian genome, of the anemone Actinernus sp., which was 1.39 Gbp in length and contained 44 970 gene models including 14 806 tRNA genes and 30 164 protein-coding genes. Analyses of homeobox genes revealed the longest chromosome hosts a mega-array of Hox cluster, HoxL, NK cluster and NKL homeobox genes; until now, such an array has only been hypothesized to have existed in ancient ancestral genomes. In addition to this striking arrangement of homeobox genes, analyses of microRNAs revealed cnidarian-specific complements that are distinctive for nested clades of these animals, presumably reflecting the progressive evolution of the gene regulatory networks in which they are embedded. Also, compared with other sea anemones, circadian rhythm genes were lost in Actinernus sp., which likely reflects adaptation to living in the dark. This high-quality genome of a deep-sea cnidarian thus reveals some of the likely molecular adaptations of this ecologically important group of metazoans to the extreme deep-sea environment. It also deepens our understanding of the evolution of genome content and organization of animals in general and cnidarians in particular, specifically from the viewpoint of key developmental control genes like the homeobox-encoding genes, where we find an array of genes that until now has only been hypothesized to have existed in the ancient ancestor that pre-dated both the cnidarians and bilaterians.

The genome sequence of the centipede Strigamia acuminata (Leach, 1816)

We present a genome assembly from an individual male Strigamia acuminata (centipede; Arthropoda; Chilopoda; Geophilomorpha; Geophilidae; Linotaeniinae). The genome sequence is 237.5 megabases in span. Most of the assembly is scaffolded into 11 chromosomal pseudomolecules, including the X and Y sex chromosomes. The mitochondrial genome has also been assembled and is 15.07 kilobases in length.

Genomic and transcriptomic analyses support a silk gland origin of spider venom glands

BACKGROUND

Spiders comprise a hyperdiverse lineage of predators with venom systems, yet the origin of functionally novel spider venom glands remains unclear. Previous studies have hypothesized that spider venom glands originated from salivary glands or evolved from silk-producing glands present in early chelicerates. However, there is insufficient molecular evidence to indicate similarity among them. Here, we provide comparative analyses of genome and transcriptome data from various lineages of spiders and other arthropods to advance our understanding of spider venom gland evolution.

RESULTS

We generated a chromosome-level genome assembly of a model spider species, the common house spider (Parasteatoda tepidariorum). Module preservation, GO semantic similarity, and differentially upregulated gene similarity analyses demonstrated a lower similarity in gene expressions between the venom glands and salivary glands compared to the silk glands, which questions the validity of the salivary gland origin hypothesis but unexpectedly prefers to support the ancestral silk gland origin hypothesis. The conserved core network in the venom and silk glands was mainly correlated with transcription regulation, protein modification, transport, and signal transduction pathways. At the genetic level, we found that many genes in the venom gland-specific transcription modules show positive selection and upregulated expressions, suggesting that genetic variation plays an important role in the evolution of venom glands.

CONCLUSIONS

This research implies the unique origin and evolutionary path of spider venom glands and provides a basis for understanding the diverse molecular characteristics of venom systems.

The enormous repetitive Antarctic krill genome reveals environmental adaptations and population insights

Antarctic krill (Euphausia superba) is Earth's most abundant wild animal, and its enormous biomass is vital to the Southern Ocean ecosystem. Here, we report a 48.01-Gb chromosome-level Antarctic krill genome, whose large genome size appears to have resulted from inter-genic transposable element expansions. Our assembly reveals the molecular architecture of the Antarctic krill circadian clock and uncovers expanded gene families associated with molting and energy metabolism, providing insights into adaptations to the cold and highly seasonal Antarctic environment. Population-level genome re-sequencing from four geographical sites around the Antarctic continent reveals no clear population structure but highlights natural selection associated with environmental variables. An apparent drastic reduction in krill population size 10 mya and a subsequent rebound 100 thousand years ago coincides with climate change events. Our findings uncover the genomic basis of Antarctic krill adaptations to the Southern Ocean and provide valuable resources for future Antarctic research.

Comparative genomics of Hox and ParaHox genes among major lineages of Branchiopoda with emphasis on tadpole shrimps

Hox and ParaHox genes (HPHGs) are key developmental genes that pattern regional identity along the anterior–posterior body axis of most animals. Here, we identified HPHGs in tadpole shrimps (Pancrustacea, Branchiopoda, Notostraca), an iconic example of the so-called “living fossils” and performed a comparative genomics analysis of HPHGs and the Hox cluster among major branchiopod lineages. Notostraca possess the entire Hox complement, and the Hox cluster seems to be split into two different subclusters, although we were not able to support this finding with chromosome-level assemblies. However, the genomic structure of Hox genes in Notostraca appears more derived than that of Daphnia spp., which instead retains the plesiomorphic condition of a single compact cluster. Spinicaudata and Artemia franciscana show instead a Hox cluster subdivided across two or more genomic scaffolds with some orthologs either duplicated or missing. Yet, branchiopod HPHGs are similar among the various clades in terms of both intron length and number, as well as in their pattern of molecular evolution. Sequence substitution rates are in fact generally similar for most of the branchiopod Hox genes and the few differences we found cannot be traced back to natural selection, as they are not associated with any signals of diversifying selection or substantial switches in selective modes. Altogether, these findings do not support a significant stasis in the Notostraca Hox cluster and further confirm how morphological evolution is not tightly associated with genome dynamics.

Diversity, duplication, and genomic organization of homeobox genes in Lepidoptera

Homeobox genes encode transcription factors with essential roles in patterning and cell fate in developing animal embryos. Many homeobox genes, including Hox and NK genes, are arranged in gene clusters, a feature likely related to transcriptional control. Sparse taxon sampling and fragmentary genome assemblies mean that little is known about the dynamics of homeobox gene evolution across Lepidoptera or about how changes in homeobox gene number and organization relate to diversity in this large order of insects. Here we analyze an extensive data set of high-quality genomes to characterize the number and organization of all homeobox genes in 123 species of Lepidoptera from 23 taxonomic families. We find most Lepidoptera have around 100 homeobox loci, including an unusual Hox gene cluster in which the lab gene is repositioned and the ro gene is next to pb A topologically associating domain spans much of the gene cluster, suggesting deep regulatory conservation of the Hox cluster arrangement in this insect order. Most Lepidoptera have four Shx genes, divergent zen-derived loci, but these loci underwent dramatic duplication in several lineages, with some moths having over 165 homeobox loci in the Hox gene cluster; this expansion is associated with local LINE element density. In contrast, the NK gene cluster content is more stable, although there are differences in organization compared with other insects, as well as major rearrangements within butterflies. Our analysis represents the first description of homeobox gene content across the order Lepidoptera, exemplifying the potential of newly generated genome assemblies for understanding genome and gene family evolution.

The diversity of invertebrate visual opsins spanning Protostomia, Deuterostomia, and Cnidaria

Across eumetazoans, the ability to perceive and respond to visual stimuli is largely mediated by opsins, a family of proteins belonging to the G protein-coupled receptor (GPCR) superclass. Lineage-specific gains and losses led to a striking diversity in the numbers, types, and spectral sensitivities conferred by visual opsin gene expression. Here, we review the diversity of visual opsins and differences in opsin gene expression from well-studied protostome, invertebrate deuterostome, and cnidarian groups. We discuss the functional significance of opsin expression differences and spectral tuning among lineages. In some cases, opsin evolution has been linked to the detection of relevant visual signals, including sexually selected color traits and host plant features. In other instances, variation in opsins has not been directly linked to functional or ecological differences. Overall, the array of opsin expression patterns and sensitivities across invertebrate lineages highlight the diversity of opsins in the eumetazoan ancestor and the labile nature of opsins over evolutionary time.

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.

Comparative transcriptomics analyses of chemosensory genes of antenna in male red swamp crayfish Procambarus clarkii

The red swamp crayfish, Procambarus clarkii, is a globally invasive species and has caused huge damage to aquaculture, biodiversity, and ecology worldwide. Antenna-expressed receptors are important for P. clarkii to detect chemosensory cues for mate attraction. In this study, we tested the behavior of male P. clarkii to the conditioned water from female P. clarkii during the mating and non-mating periods, and performed RNA sequencing to investigate the chemosensory-related genes of the antenna of male P. clarkii. The results of the behavioral assay have shown that for the female-conditioned water, male P. clarkii within the mating period can be significantly attracted, but not during the non-mating period. This suggested that the expressions of chemosensory-related genes in the antenna of male P. clarkii may change significantly with mating seasonal variation. Antenna transcriptomes found that a total of 59,218 unigenes with an average length of 1,056.41 bp, and 4,889 differentially expressed unigenes (DEGs), among which 2,128 were upregulated, while 2,761 were downregulated were obtained. A total of 12 upregulated and nine downregulated DEGs were associated with chemical reception, including four ionotropic receptors (IRs) or ionotropic glutamate receptors (iGluRs), eight G-protein-coupled receptors (GPCRs), five transient receptor potential channels (TRP channels), one sodium–calcium exchanger, one isomerase, and two uncharacterized proteins (chemosensory proteins-like, CSPs). CSPs were preliminarily classified as pheromone receptors in the antenna of male P. clarkii. Furthermore, the calcium transduction-related pathways may play an important role in the sex pheromone reception of the male P. clarkii’s antenna. The results of quantitative real-time reverse transcriptase PCR (RT-qPCR) showed that the trends of expression of eight selected unigenes were consistent with RNA-Seq results. Our results provide more comprehensive data for chemical communication mechanisms after P. clarkii enter the mating period and eventually would develop better control strategies in further.

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.

Rethinking Sesquiterpenoids: A Widespread Hormone in Animals

The sesquiterpenoid hormone juvenile hormone (JH) controls development, reproduction, and metamorphosis in insects, and has long been thought to be confined to the Insecta. While it remains true that juvenile hormone is specifically synthesized in insects, other types or forms of sesquiterpenoids have also been discovered in distantly related animals, such as the jellyfish. Here, we combine the latest literature and annotate the sesquiterpenoid biosynthetic pathway genes in different animal genomes. We hypothesize that the sesquiterpenoid hormonal system is an ancestral system established in an animal ancestor and remains widespread in many animals. Different animal lineages have adapted different enzymatic routes from a common pathway, with cnidarians producing farnesoic acid (FA); non-insect protostomes and non-vertebrate deuterostomes such as cephalochordate and echinoderm synthesizing FA and methyl farnesoate (MF); and insects producing FA, MF, and JH. Our hypothesis revolutionizes the current view on the sesquiterpenoids in the metazoans, and forms a foundation for a re-investigation of the roles of this important and yet neglected type of hormone in different animals.

Chromosome-level genome of the globe skimmer dragonfly (Pantala flavescens)

BACKGROUND

The globe skimmer dragonfly (Pantala flavescens) is a notable Odonata insect distributed in nature fields and farmlands worldwide, and it is commonly recognized as a natural enemy because it preys on agricultural pests and health pests. As one of the sister groups of winged insects, odonatan species are key to understanding the evolution of insect wings.

FINDINGS

We present a high-quality reference genome of P. flavescens, which is the first chromosome-level genome in the Palaeoptera (Odonata and Ephemeroptera). The assembled genome size was 662 Mb, with a contig N50 of 16.2 Mb. Via Hi-C scaffolding, 648 Mb (97.9%) of contig sequences were clustered, ordered, and assembled into 12 large scaffolds, each corresponding to a natural chromosome. The X chromosome was identified by sequence coverage depth. The repetitive sequences and gene density of the X chromosome are similar to those of autosomal sequences, but the X chromosome shows a much lower degree of heterozygosity. Our analysis shows that the effective population size experienced 3 declining events, which may have been caused by climate change and environmental pollution.

CONCLUSIONS

The genome of P. flavescens provides more information on the biology and evolution of insects and will help for the use of this species in pest control.

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.

Genome of the ramshorn snail Biomphalaria straminea-an obligate intermediate host of schistosomiasis

BACKGROUND

Schistosomiasis, or bilharzia, is a parasitic disease caused by trematode flatworms of the genus Schistosoma. Infection by Schistosoma mansoni in humans results when cercariae emerge into water from freshwater snails in the genus Biomphalaria and seek out and penetrate human skin. The snail Biomphalaria straminea is native to South America and is now also present in Central America and China, and represents a potential vector host for spreading schistosomiasis. To date, genomic information for the genus is restricted to the neotropical species Biomphalaria glabrata. This limits understanding of the biology and management of other schistosomiasis vectors, such as B. straminea.

FINDINGS

Using a combination of Illumina short-read, 10X Genomics linked-read, and Hi-C sequencing data, our 1.005 Gb B. straminea genome assembly is of high contiguity, with a scaffold N50 of 25.3 Mb. Transcriptomes from adults were also obtained. Developmental homeobox genes, hormonal genes, and stress-response genes were identified, and repeat content was annotated (40.68% of genomic content). Comparisons with other mollusc genomes (including Gastropoda, Bivalvia, and Cephalopoda) revealed syntenic conservation, patterns of homeobox gene linkage indicative of evolutionary changes to gene clusters, expansion of heat shock protein genes, and the presence of sesquiterpenoid and cholesterol metabolic pathway genes in Gastropoda. In addition, hormone treatment together with RT-qPCR assay reveal a sesquiterpenoid hormone responsive system in B. straminea, illustrating that this renowned insect hormonal system is also present in the lophotrochozoan lineage.

CONCLUSION

This study provides the first genome assembly for the snail B. straminea and offers an unprecedented opportunity to address a variety of phenomena related to snail vectors of schistosomiasis, as well as evolutionary and genomics questions related to molluscs more widely.

Deeply conserved synteny and the evolution of metazoan chromosomes

Animal genomes show networks of deeply conserved gene linkages whose phylogenetic scope and chromosomal context remain unclear. Here, we report chromosome-scale conservation of synteny among bilaterians, cnidarians, and sponges and use comparative analysis to reconstruct ancestral chromosomes across major animal groups. Comparisons among diverse metazoans reveal the processes of chromosome evolution that produced contemporary karyotypes from their Precambrian progenitors. On the basis of these findings, we introduce a simple algebraic representation of chromosomal change and use it to establish a unified systematic framework for metazoan chromosome evolution. We find that fusion-with-mixing, a previously unappreciated mode of chromosome change, has played a central role. We find that relicts of several metazoan chromosomal units are preserved in unicellular eukaryotes. These conserved pre-metazoan linkages include the chromosomal unit that encodes the most diverse set of metazoan homeobox genes, suggesting a candidate genomic context for the early diversification of this key gene family.

Multiscale analysis of the randomization limits of the chromosomal gene organization between Lepidoptera and Diptera

How chromosome gene organization and gene content evolve among distantly related and structurally malleable genomes remains unresolved. This is particularly the case when considering different insect orders. We have compared the highly contiguous genome assemblies of the lepidopteran Danaus plexippus and the dipteran Drosophila melanogaster, which shared a common ancestor around 290 Ma. The gene content of 23 out of 30 D. plexippus chromosomes was significantly associated with one or two of the six chromosomal elements of the Drosophila genome, denoting common ancestry. Despite the phylogenetic distance, 9.6% of the 1-to-1 orthologues still reside within the same ancestral genome neighbourhood. Furthermore, the comparison D. plexippus-Bombyx mori indicated that the rates of chromosome repatterning are lower in Lepidoptera than in Diptera, although still within the same order of magnitude. Concordantly, 14 developmental gene clusters showed a higher tendency to retain full or partial clustering in D. plexippus, further supporting that the physical association between the SuperHox and NK clusters existed in the ancestral bilaterian. Our results illuminate the scope and limits of the evolution of the gene organization and content of the ancestral chromosomes to the Lepidoptera and Diptera while helping reconstruct portions of the genome in their most recent common ancestor.

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.

Transcriptomic analysis of brine shrimp Artemia franciscana across a wide range of salinities

Brine shrimp Artemia franciscana, a commercially important species, can thrive in a wide range of salinities and is commonly found in hypersaline lakes and solar salterns. Transcriptome analysis can enhance the understanding of the adaptative mechanisms of brine shrimp in aquaculture. RNA sequencing (RNAseq) data was generated from A. franciscana adults that were salt-adapted for 2-4 weeks at five salinities: 35, 50, 100, 150, and 230 psu. Long-read isoform sequencing (IsoSeq) data was used to construct a high-quality transcriptome assembly. Also, the gene expression patterns in A. franciscana adults were examined. Notably, the transcriptional response of A. franciscana's acclimation to intermediate salinities (50-150 psu) displayed frequently and differentially U-shaped or inverted U-shaped expression patterns. In addition, the types of genes showing two nonmonotonic expression patterns were distinct from each other. The coordinated shifts in gene expression suggest different homeostatic strategies of A. franciscana at specific salinities; such strategies may enhance population fitness at extreme salinities. Our study should promote a scientific concept for the gene expression patterns of A. franciscana along a broad salinity gradient, and a variety of salinity and prey should be monitored for testing the gene expression pattern of this important aquaculture species.

Widespread retention of ohnologs in key developmental gene families following whole-genome duplication in arachnopulmonates

Whole-genome duplications (WGDs) have occurred multiple times during animal evolution, including in lineages leading to vertebrates, teleosts, horseshoe crabs, and arachnopulmonates. These dramatic events initially produce a wealth of new genetic material, generally followed by extensive gene loss. It appears, however, that developmental genes such as homeobox genes, signaling pathway components and microRNAs are frequently retained as duplicates (so-called ohnologs) following WGD. These not only provide the best evidence for WGD, but an opportunity to study its evolutionary consequences. Although these genes are well studied in the context of vertebrate WGD, similar comparisons across the extant arachnopulmonate orders are patchy. We sequenced embryonic transcriptomes from two spider species and two amblypygid species and surveyed three important gene families, Hox, Wnt, and frizzled, across these and 12 existing transcriptomic and genomic resources for chelicerates. We report extensive retention of putative ohnologs, further supporting the ancestral arachnopulmonate WGD. We also found evidence of consistent evolutionary trajectories in Hox and Wnt gene repertoires across three of the six arachnopulmonate orders, with interorder variation in the retention of specific paralogs. We identified variation between major clades in spiders and are better able to reconstruct the chronology of gene duplications and losses in spiders, amblypygids, and scorpions. These insights shed light on the evolution of the developmental toolkit in arachnopulmonates, highlight the importance of the comparative approach within lineages, and provide substantial new transcriptomic data for future study.

PhylomeDB V5: an expanding repository for genome-wide catalogues of annotated gene phylogenies

PhylomeDB is a unique knowledge base providing public access to minable and browsable catalogues of pre-computed genome-wide collections of annotated sequences, alignments and phylogenies (i.e. phylomes) of homologous genes, as well as to their corresponding phylogeny-based orthology and paralogy relationships. In addition, PhylomeDB trees and alignments can be downloaded for further processing to detect and date gene duplication events, infer past events of inter-species hybridization and horizontal gene transfer, as well as to uncover footprints of selection, introgression, gene conversion, or other relevant evolutionary processes in the genes and organisms of interest. Here, we describe the latest evolution of PhylomeDB (version 5). This new version includes a newly implemented web interface and several new functionalities such as optimized searching procedures, the possibility to create user-defined phylome collections, and a fully redesigned data structure. This release also represents a significant core data expansion, with the database providing access to 534 phylomes, comprising over 8 million trees, and homology relationships for genes in over 6000 species. This makes PhylomeDB the largest and most comprehensive public repository of gene phylogenies. PhylomeDB is available at http://www.phylomedb.org.

The Nereid on the rise: Platynereis as a model system

The Nereid Platynereis dumerilii (Audouin and Milne Edwards (Annales des Sciences Naturelles 1:195-269, 1833) is a marine annelid that belongs to the Nereididae, a family of errant polychaete worms. The Nereid shows a pelago-benthic life cycle: as a general characteristic for the superphylum of Lophotrochozoa/Spiralia, it has spirally cleaving embryos developing into swimming trochophore larvae. The larvae then metamorphose into benthic worms living in self-spun tubes on macroalgae. Platynereis is used as a model for genetics, regeneration, reproduction biology, development, evolution, chronobiology, neurobiology, ecology, ecotoxicology, and most recently also for connectomics and single-cell genomics. Research on the Nereid started with studies on eye development and spiralian embryogenesis in the nineteenth and early twentieth centuries. Transitioning into the molecular era, Platynereis research focused on posterior growth and regeneration, neuroendocrinology, circadian and lunar cycles, fertilization, and oocyte maturation. Other work covered segmentation, photoreceptors and other sensory cells, nephridia, and population dynamics. Most recently, the unique advantages of the Nereid young worm for whole-body volume electron microscopy and single-cell sequencing became apparent, enabling the tracing of all neurons in its rope-ladder-like central nervous system, and the construction of multimodal cellular atlases. Here, we provide an overview of current topics and methodologies for P. dumerilii, with the aim of stimulating further interest into our unique model and expanding the active and vibrant Platynereis community.

Maternal Transcripts of Hox Genes Are Found in Oocytes of Platynereis dumerilii (Annelida, Nereididae)

Hox genes are some of the best studied developmental control genes. In the overwhelming majority of bilateral animals, these genes are sequentially activated along the main body axis during the establishment of the ground plane, i.e., at the moment of gastrulation. Their activation is necessary for the correct differentiation of cell lines, but at the same time it reduces the level of stemness. That is why the chromatin of Hox loci in the pre-gastrulating embryo is in a bivalent state. It carries both repressive and permissive epigenetic markers at H3 histone residues, leading to transcriptional repression. There is a paradox that maternal RNAs, and in some cases the proteins of the Hox genes, are present in oocytes and preimplantation embryos in mammals. Their functions should be different from the zygotic ones and have not been studied to date. Our object is the errant annelid Platynereis dumerilii. This model is convenient for studying new functions and mechanisms of regulation of Hox genes, because it is incomparably simpler than laboratory vertebrates. Using a standard RT-PCR on cDNA template which was obtained by reverse transcription using random primers, we found that maternal transcripts of almost all Hox genes are present in unfertilized oocytes of worm. We assessed the localization of these transcripts using WMISH.

The genome of a daddy-long-legs (Opiliones) illuminates the evolution of arachnid appendages

Chelicerate arthropods exhibit dynamic genome evolution, with ancient whole-genome duplication (WGD) events affecting several orders. Yet, genomes remain unavailable for a number of poorly studied orders, such as Opiliones (daddy-long-legs), which has hindered comparative study. We assembled the first harvestman draft genome for the species Phalangium opilio, which bears elongate, prehensile appendages, made possible by numerous distal articles called tarsomeres. Here, we show that the genome of P. opilio exhibits a single Hox cluster and no evidence of WGD. To investigate the developmental genetic basis for the quintessential trait of this group-the elongate legs-we interrogated the function of the Hox genes Deformed (Dfd) and Sex combs reduced (Scr), and a homologue of Epidermal growth factor receptor (Egfr). Knockdown of Dfd incurred homeotic transformation of two pairs of legs into pedipalps, with dramatic shortening of leg segments in the longest leg pair, whereas homeosis in L3 is only achieved upon double Dfd + Scr knockdown. Knockdown of Egfr incurred shortened appendages and the loss of tarsomeres. The similarity of Egfr loss-of-function phenotypic spectra in insects and this arachnid suggest that repeated cooption of EGFR signalling underlies the independent gains of supernumerary tarsomeres across the arthropod tree of life.

There and Back Again: Hox Clusters Use Both DNA Strands

Bilaterian animals operate the clusters of Hox genes through a rich repertoire of diverse mechanisms. In this review, we will summarize and analyze the accumulated data concerning long non-coding RNAs (lncRNAs) that are transcribed from sense (coding) DNA strands of Hox clusters. It was shown that antisense regulatory RNAs control the work of Hox genes in cis and trans, participate in the establishment and maintenance of the epigenetic code of Hox loci, and can even serve as a source of regulatory peptides that switch cellular energetic metabolism. Moreover, these molecules can be considered as a force that consolidates the cluster into a single whole. We will discuss the examples of antisense transcription of Hox genes in well-studied systems (cell cultures, morphogenesis of vertebrates) and bear upon some interesting examples of antisense Hox RNAs in non-model Protostomia.

Sawfly Genomes Reveal Evolutionary Acquisitions That Fostered the Mega-Radiation of Parasitoid and Eusocial Hymenoptera

The tremendous diversity of Hymenoptera is commonly attributed to the evolution of parasitoidism in the last common ancestor of parasitoid sawflies (Orussidae) and wasp-waisted Hymenoptera (Apocrita). However, Apocrita and Orussidae differ dramatically in their species richness, indicating that the diversification of Apocrita was promoted by additional traits. These traits have remained elusive due to a paucity of sawfly genome sequences, in particular those of parasitoid sawflies. Here, we present comparative analyses of draft genomes of the primarily phytophagous sawfly Athalia rosae and the parasitoid sawfly Orussus abietinus. Our analyses revealed that the ancestral hymenopteran genome exhibited traits that were previously considered unique to eusocial Apocrita (e.g., low transposable element content and activity) and a wider gene repertoire than previously thought (e.g., genes for CO₂ detection). Moreover, we discovered that Apocrita evolved a significantly larger array of odorant receptors than sawflies, which could be relevant to the remarkable diversification of Apocrita by enabling efficient detection and reliable identification of hosts.

Cell-specific expression and individual function of prohormone convertase PC1/3 in Tribolium larval growth highlights major evolutionary changes between beetle and fly neuroendocrine systems

Background

The insect neuroendocrine system acts in the regulation of physiology, development and growth. Molecular evolution of this system hence has the potential to allow for major biological differences between insect groups. Two prohormone convertases, PC1/3 and PC2, are found in animals and both function in the processing of neuropeptide precursors in the vertebrate neurosecretory pathway. Whereas PC2-function is conserved between the fly Drosophila and vertebrates, ancestral PC1/3 was lost in the fly lineage and has not been functionally studied in any protostome.

Results

In order to understand its original functions and the changes accompanying the gene loss in the fly, we investigated PC1/3 and PC2 expression and function in the beetle Tribolium castaneum. We found that PC2 is broadly expressed in the nervous system, whereas surprisingly, PC1/3 expression is restricted to specific cell groups in the posterior brain and suboesophageal ganglion. Both proteases have parallel but non-redundant functions in adult beetles’ viability and fertility. Female infertility following RNAi is caused by a failure to deposit sufficient yolk to the developing oocytes. Larval RNAi against PC2 produced moulting defects where the larvae were not able to shed their old cuticle. This ecdysis phenotype was also observed in a small subset of PC1/3 knockdown larvae and was strongest in a double knockdown. Unexpectedly, most PC1/3-RNAi larvae showed strongly reduced growth, but went through larval moults despite minimal to zero weight gain.

Conclusions

The cell type-specific expression of PC1/3 and its essential requirement for larval growth highlight the important role of this gene within the insect neuroendocrine system. Genomic conservation in most insect groups suggests that it has a comparable individual function in other insects as well, which has been replaced by alternative mechanisms in flies.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13227-021-00179-w.

A Review of Pedal Peptide/Orcokinin-type Neuropeptides

Neuropeptides are endogenous active substances that play important roles in a number of physiological processes and are ubiquitous in the nervous tissue in vivo. The gene encoding pedal peptide/orcokinin-type (PP/OK-type) neuropeptide is an important member of the neuropeptide gene family and is ubiquitous in invertebrates of Bilateria; orcokinin (OK) is mainly found in Arthropoda, while pedal peptide (PP) is mainly found in Mollusca. OK and PP are also present in other animals. PP/OK-type neuropeptides are a kind of multifunctional neuropeptides predominantly expressed in the nervous tissue and play important roles in the nerve regulation of movement. Moreover, OK has a number of other physiological functions. This review describes the distribution, expression, function and maturation of PP/OK-type neuropeptides to facilitate investigations of new functions and receptors of PP/OK-type neuropeptides, providing the theoretical foundation for the potential use of PP/OK-type neuropeptides in the prevention and control of agricultural and forestry pests, as an additive for skin care products and in the screening of drugs for the treatment of diabetes.

Taxonomic Sampling and Rare Genomic Changes Overcome Long-Branch Attraction in the Phylogenetic Placement of Pseudoscorpions

Long-branch attraction is a systematic artifact that results in erroneous groupings of fast-evolving taxa. The combination of short, deep internodes in tandem with long-branch attraction artifacts has produced empirically intractable parts of the Tree of Life. One such group is the arthropod subphylum Chelicerata, whose backbone phylogeny has remained unstable despite improvements in phylogenetic methods and genome-scale data sets. Pseudoscorpion placement is particularly variable across data sets and analytical frameworks, with this group either clustering with other long-branch orders or with Arachnopulmonata (scorpions and tetrapulmonates). To surmount long-branch attraction, we investigated the effect of taxonomic sampling via sequential deletion of basally branching pseudoscorpion superfamilies, as well as varying gene occupancy thresholds in supermatrices. We show that concatenated supermatrices and coalescent-based summary species tree approaches support a sister group relationship of pseudoscorpions and scorpions, when more of the basally branching taxa are sampled. Matrix completeness had demonstrably less influence on tree topology. As an external arbiter of phylogenetic placement, we leveraged the recent discovery of an ancient genome duplication in the common ancestor of Arachnopulmonata as a litmus test for competing hypotheses of pseudoscorpion relationships. We generated a high-quality developmental transcriptome and the first genome for pseudoscorpions to assess the incidence of arachnopulmonate-specific duplications (e.g., homeobox genes and miRNAs). Our results support the inclusion of pseudoscorpions in Arachnopulmonata (new definition), as the sister group of scorpions. Panscorpiones (new name) is proposed for the clade uniting Scorpiones and Pseudoscorpiones.

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.

Unconventional Actions of Glycoprotein Hormone Subunits: A Comprehensive Review

The glycoprotein hormones (GPH) are heterodimers composed of a common α subunit and a specific β subunit. They act by activating specific leucine-rich repeat G protein-coupled receptors. However, individual subunits have been shown to elicit responses in cells devoid of the receptor for the dimeric hormones. The α subunit is involved in prolactin production from different tissues. The human chorionic gonadotropin β subunit (βhCG) plays determinant roles in placentation and in cancer development and metastasis. A truncated form of the thyrotropin (TSH) β subunit is also reported to have biological effects. The GPH α- and β subunits are derived from precursor genes (gpa and gpb, respectively), which are expressed in most invertebrate species and are still represented in vertebrates as GPH subunit paralogs (gpa2 and gpb5, respectively). No specific receptor has been found for the vertebrate GPA2 and GPB5 even if their heterodimeric form is able to activate the TSH receptor in mammals. Interestingly, GPA and GPB are phylogenetically and structurally related to cysteine-knot growth factors (CKGF) and particularly to a group of antagonists that act independently on any receptor. This review article summarizes the observed actions of individual GPH subunits and presents the current hypotheses of how these actions might be induced. New approaches are also proposed in light of the evolutionary relatedness with antagonists of the CKGF family of proteins.

Four myriapod relatives - but who are sisters? No end to debates on relationships among the four major myriapod subgroups

Background

Phylogenetic relationships among the myriapod subgroups Chilopoda, Diplopoda, Symphyla and Pauropoda are still not robustly resolved. The first phylogenomic study covering all subgroups resolved phylogenetic relationships congruently to morphological evidence but is in conflict with most previously published phylogenetic trees based on diverse molecular data. Outgroup choice and long-branch attraction effects were stated as possible explanations for these incongruencies. In this study, we addressed these issues by extending the myriapod and outgroup taxon sampling using transcriptome data.

Results

We generated new transcriptome data of 42 panarthropod species, including all four myriapod subgroups and additional outgroup taxa. Our taxon sampling was complemented by published transcriptome and genome data resulting in a supermatrix covering 59 species. We compiled two data sets, the first with a full coverage of genes per species (292 single-copy protein-coding genes), the second with a less stringent coverage (988 genes). We inferred phylogenetic relationships among myriapods using different data types, tree inference, and quartet computation approaches. Our results unambiguously support monophyletic Mandibulata and Myriapoda. Our analyses clearly showed that there is strong signal for a single unrooted topology, but a sensitivity of the position of the internal root on the choice of outgroups. However, we observe strong evidence for a clade Pauropoda+Symphyla, as well as for a clade Chilopoda+Diplopoda.

Conclusions

Our best quartet topology is incongruent with current morphological phylogenies which were supported in another phylogenomic study. AU tests and quartet mapping reject the quartet topology congruent to trees inferred with morphological characters. Moreover, quartet mapping shows that confounding signal present in the data set is sufficient to explain the weak signal for the quartet topology derived from morphological characters. Although outgroup choice affects results, our study could narrow possible trees to derivatives of a single quartet topology. For highly disputed relationships, we propose to apply a series of tests (AU and quartet mapping), since results of such tests allow to narrow down possible relationships and to rule out confounding signal.

Ecdysis triggering hormone modulates molt behaviour in the redclaw crayfish Cherax quadricarinatus, providing a mechanistic evidence for conserved function in molt regulation across Pancrustacea

Molting enables growth and development across ecdysozoa. The molting process is strictly controlled by hormones - ecdysteroids. Ecdysteroidogenesis occurs in theprothoracic glands and stimulated by prothoracicotropic hormone in insects, while it ensues in the Y-organ and regulated by the molt inhibiting hormone in crustaceans. A peak in ecdysteroids in the hemolymph induces a cascade of multiple neuropeptides including Ecdysis Triggering Hormone (ETH) and Corazonin. The role of ETH is well defined in controlling the molt process in insects, but it is yet to be defined in crustaceans. In this study, we investigated the behavioral response of intermolt crayfish to ETH and Corazonin injections as well as the impact of ETH on the molt period using in vivo assays. Injection of Corazonin and ETH resulted in a clear and immediate eye twitching response to these two neuropeptides. The Corazonin injection induced eye twitching in slow and asynchronous manner, while ETH injection caused eye twitching in a relatively fast and synchronous way. A single injection of ETH to crayfish resulted in a remarkable prolong molt period, at twice the normal molting cycle, suggesting that ETH plays a key role in controlling the molt cycle in decapod crustaceans. Given the key significance of ETH in molt regulation and its plausible application in pest control, we characterized ETH across the pancrustacean orders. Bioinformatic analysis shows the mature ETH sequence is identical in all studied decapod species. ETH can be classified into specific groups based on the associated motif in each insect order and shows an insect motif -KxxPRx to be conserved in crustaceans.

Genome streamlining in a minute herbivore that manipulates its host plant

The tomato russet mite, Aculops lycopersici, is among the smallest animals on earth. It is a worldwide pest on tomato and can potently suppress the host's natural resistance. We sequenced its genome, the first of an eriophyoid, and explored whether there are genomic features associated with the mite's minute size and lifestyle. At only 32.5 Mb, the genome is the smallest yet reported for any arthropod and, reminiscent of microbial eukaryotes, exceptionally streamlined. It has few transposable elements, tiny intergenic regions, and is remarkably intron-poor, as more than 80% of coding genes are intronless. Furthermore, in accordance with ecological specialization theory, this defense-suppressing herbivore has extremely reduced environmental response gene families such as those involved in chemoreception and detoxification. Other losses associate with this species' highly derived body plan. Our findings accelerate the understanding of evolutionary forces underpinning metazoan life at the limits of small physical and genome size.

Reconstruction of ancient homeobox gene linkages inferred from a new high-quality assembly of the Hong Kong oyster (Magallana hongkongensis) genome

BACKGROUND

Homeobox-containing genes encode crucial transcription factors involved in animal, plant and fungal development, and changes to homeobox genes have been linked to the evolution of novel body plans and morphologies. In animals, some homeobox genes are clustered together in the genome, either as remnants from ancestral genomic arrangements, or due to coordinated gene regulation. Consequently, analyses of homeobox gene organization across animal phylogeny provide important insights into the evolution of genome organization and developmental gene control, and their interaction. However, homeobox gene organization remains to be fully elucidated in several key animal ancestors, including those of molluscs, lophotrochozoans and bilaterians.

RESULTS

Here, we present a high-quality chromosome-level genome assembly of the Hong Kong oyster, Magallana hongkongensis (2n = 20), for which 93.2% of the genomic sequences are contained on 10 pseudomolecules (~ 758 Mb, scaffold N50 = 72.3 Mb). Our genome assembly was scaffolded using Hi-C reads, facilitating a larger scaffold size compared to the recently published M. hongkongensis genome of Peng et al. (Mol Ecol Resources, 2020), which was scaffolded using the Crassostrea gigas assembly. A total of 46,963 predicted gene models (45,308 protein coding genes) were incorporated in our genome, and genome completeness estimated by BUSCO was 94.6%. Homeobox gene linkages were analysed in detail relative to available data for other mollusc lineages.

CONCLUSIONS

The analyses performed in this study and the accompanying genome sequence provide important genetic resources for this economically and culturally valuable oyster species, and offer a platform to improve understanding of animal biology and evolution more generally. Transposable element content is comparable to that found in other mollusc species, contrary to the conclusion of another recent analysis. Also, our chromosome-level assembly allows the inference of ancient gene linkages (synteny) for the homeobox-containing genes, even though a number of the homeobox gene clusters, like the Hox/ParaHox clusters, are undergoing dispersal in molluscs such as this oyster.

Millipede genomes reveal unique adaptations during myriapod evolution

The Myriapoda, composed of millipedes and centipedes, is a fascinating but poorly understood branch of life, including species with a highly unusual body plan and a range of unique adaptations to their environment. Here, we sequenced and assembled 2 chromosomal-level genomes of the millipedes Helicorthomorpha holstii (assembly size = 182 Mb; shortest scaffold/contig length needed to cover 50% of the genome [N50] = 18.11 Mb mainly on 8 pseudomolecules) and Trigoniulus corallinus (assembly size = 449 Mb, N50 = 26.78 Mb mainly on 17 pseudomolecules). Unique genomic features, patterns of gene regulation, and defence systems in millipedes, not observed in other arthropods, are revealed. Both repeat content and intron size are major contributors to the observed differences in millipede genome size. Tight Hox and the first loose ecdysozoan ParaHox homeobox clusters are identified, and a myriapod-specific genomic rearrangement including Hox3 is also observed. The Argonaute (AGO) proteins for loading small RNAs are duplicated in both millipedes, but unlike in insects, an AGO duplicate has become a pseudogene. Evidence of post-transcriptional modification in small RNAs—including species-specific microRNA arm switching—providing differential gene regulation is also obtained. Millipedes possesses a unique ozadene defensive gland unlike the venomous forcipules found in centipedes. We identify sets of genes associated with the ozadene that play roles in chemical defence as well as antimicrobial activity. Macro-synteny analyses revealed highly conserved genomic blocks between the 2 millipedes and deuterostomes. Collectively, our analyses of millipede genomes reveal that a series of unique adaptations have occurred in this major lineage of arthropod diversity. The 2 high-quality millipede genomes provided here shed new light on the conserved and lineage-specific features of millipedes and centipedes. These findings demonstrate the importance of the consideration of both centipede and millipede genomes—and in particular the reconstruction of the myriapod ancestral situation—for future research to improve understanding of arthropod evolution, and animal evolutionary genomics more widely.

Myriapods were among the first arthropods to invade the land over 400 million years ago, and survive today as the herbivorous millipedes and venomous centipedes. This study describes the genome sequences of two millipedes, Helicorthomorpha holstii and Trigoniulus corallinus, revealing unique adaptations not found in other arthropods.

Neuropeptide and microRNA regulators of juvenile hormone production

The sesquiterpenoid juvenile hormone(s) (JHs) of insects are the primary regulators of growth, metamorphosis, and reproduction in most insect species. As a consequence, it is essential that JH production be precisely regulated so that it is present only during appropriate periods necessary for the control of these processes. The presence of JH at inappropriate times results in disruption to metamorphosis and development and, in some cases, to disturbances in female reproduction. Neuropeptides regulate the timing and production of JH by the corpora allata. Allatostatin and allatotropin were the names coined for neuropeptides that serve as inhibitors or stimulators of JH biosynthesis, respectively. Three different allatostatin neuropeptide families are capable of inhibiting juvenile hormone but only one family is utilized for that purpose dependent on the insect studied. The function of allatotropin also varies in different insects. These neuropeptides are pleiotropic in function acting on diverse physiological processes in different insects such as muscle contraction, sleep and neuromodulation. Genome projects and expression studies have assigned individual neuropeptide families to their respective receptors. An understanding of the localization of these receptors is providing clues as to how numerous peptide families might be integrated in regulating physiological functions. In recent years microRNAs have been identified that down-regulate enzymes and transcription factors that are involved in the biosynthesis and action of juvenile hormone.

Genomic resources and toolkits for developmental study of whip spiders (Amblypygi) provide insights into arachnid genome evolution and antenniform leg patterning

BACKGROUND

The resurgence of interest in the comparative developmental study of chelicerates has led to important insights, such as the discovery of a genome duplication shared by spiders and scorpions, inferred to have occurred in the most recent common ancestor of Arachnopulmonata (a clade comprising the five arachnid orders that bear book lungs). Nonetheless, several arachnid groups remain understudied in the context of development and genomics, such as the order Amblypygi (whip spiders). The phylogenetic position of Amblypygi in Arachnopulmonata posits them as an interesting group to test the incidence of the proposed genome duplication in the common ancestor of Arachnopulmonata, as well as the degree of retention of duplicates over 450 Myr. Moreover, whip spiders have their first pair of walking legs elongated and modified into sensory appendages (a convergence with the antennae of mandibulates), but the genetic patterning of these antenniform legs has never been investigated.

RESULTS

We established genomic resources and protocols for cultivation of embryos and gene expression assays by in situ hybridization to study the development of the whip spider Phrynus marginemaculatus. Using embryonic transcriptomes from three species of Amblypygi, we show that the ancestral whip spider exhibited duplications of all ten Hox genes. We deploy these resources to show that paralogs of the leg gap genes dachshund and homothorax retain arachnopulmonate-specific expression patterns in P. marginemaculatus. We characterize the expression of leg gap genes Distal-less, dachshund-1/2 and homothorax-1/2 in the embryonic antenniform leg and other appendages, and provide evidence that allometry, and by extension the antenniform leg fate, is specified early in embryogenesis.

CONCLUSION

This study is the first step in establishing P. marginemaculatus as a chelicerate model for modern evolutionary developmental study, and provides the first resources sampling whip spiders for comparative genomics. Our results suggest that Amblypygi share a genome duplication with spiders and scorpions, and set up a framework to study the genetic specification of antenniform legs. Future efforts to study whip spider development must emphasize the development of tools for functional experiments in P. marginemaculatus.

Evolutionary History of Major Chemosensory Gene Families across Panarthropoda

Chemosensory perception is a fundamental biological process of particular relevance in basic and applied arthropod research. However, apart from insects, there is little knowledge of specific molecules involved in this system, which is restricted to a few taxa with uneven phylogenetic sampling across lineages. From an evolutionary perspective, onychophorans (velvet worms) and tardigrades (water bears) are of special interest since they represent the closest living relatives of arthropods, altogether comprising the Panarthropoda. To get insights into the evolutionary origin and diversification of the chemosensory gene repertoire in panarthropods, we sequenced the antenna- and head-specific transcriptomes of the velvet worm Euperipatoides rowelli and analyzed members of all major chemosensory families in representative genomes of onychophorans, tardigrades, and arthropods. Our results suggest that the NPC2 gene family was the only family encoding soluble proteins in the panarthropod ancestor and that onychophorans might have lost many arthropod-like chemoreceptors, including the highly conserved IR25a receptor of protostomes. On the other hand, the eutardigrade genomes lack genes encoding the DEG-ENaC and CD36-sensory neuron membrane proteins, the chemosensory members of which have been retained in arthropods; these losses might be related to lineage-specific adaptive strategies of tardigrades to survive extreme environmental conditions. Although the results of this study need to be further substantiated by an increased taxon sampling, our findings shed light on the diversification of chemosensory gene families in Panarthropoda and contribute to a better understanding of the evolution of animal chemical senses.

Transcriptomic characterisation of neuropeptides and their putative cognate G protein-coupled receptors during late embryo and stage-1 juvenile development of the Aotearoa-New Zealand crayfish, Paranephrops zealandicus

We de novo assembled a transcriptome for early life-stages of the Aotearoa-New Zealand crayfish, Paranephrops zealandicus, establishing the first genetic resource for this under-developed aquaculture species and for the Paranephrops genus. Mining of this transcriptome for neuropeptides and their putative cognate G protein-coupled receptors (GPCRs) yielded a comprehensive catalogue of neuropeptides, but few putative neuropeptide GPCRs. Of the neuropeptides commonly identified from decapod transcriptomes, only crustacean female sex hormone and insulin-like peptide were absent from our trinity de novo transcriptome assembly, and also RNA-sequence reads. We identified 63 putative neuropeptide precursors from 43 families, predicted to yield 122 active peptides. Transcripts encoding 26 putative neuropeptide GPCRs were identified but were often incomplete. Putative GPCRs for 15 of the neuropeptides identified here were absent from our transcriptome and RNAseq reads. These data highlight the diverse neuropeptide systems already present at the early development life stages sampled here for P. zealandicus.

Disparate expression specificities coded by a shared Hox-C enhancer

Can a single regulatory sequence be shared by two genes undergoing functional divergence? Here we describe a single promiscuous enhancer within the Drosophila Antennapedia Complex, EO053, that directs aspects of the expression of two adjacent genes, pb (a Hox2 ortholog) and zen2 (a divergent Hox3 paralog), with disparate spatial and temporal expression patterns. We were unable to separate the pb-like and zen2-like specificities within EO053, and we identify sequences affecting both expression patterns. Importantly, genomic deletion experiments demonstrate that EO053 cooperates with additional pb- and zen2-specific enhancers to regulate the mRNA expression of both genes. We examine sequence conservation of EO053 within the Schizophora, and show that patterns of synteny between the Hox2 and Hox3 orthologs in Arthropods are consistent with a shared regulatory relationship extending prior to the Hox3/zen divergence. Thus, EO053 represents an example of two genes having evolved disparate outputs while utilizing this shared regulatory region.

Editorial note: This article has been through an editorial process in which the authors decide how to respond to the issues raised during peer review. The Reviewing Editor's assessment is that all the issues have been addressed (see decision letter).