Genome-wide analysis of homeobox genes from Mesobuthus martensii reveals Hox gene duplication in scorpions

Insight into the adaptive role of arachnid genome-wide duplication through chromosome-level genome assembly of the Western black widow spider

Although spiders are one of the most diverse groups of arthropods, the genetic architecture of their evolutionary adaptations is largely unknown. Specifically, ancient genome-wide duplication occurring during arachnid evolution ~450 mya resulted in a vast assembly of gene families, yet the extent to which selection has shaped this variation is understudied. To aid in comparative genome sequence analyses, we provide a chromosome-level genome of the Western black widow spider (Latrodectus hesperus)-a focus due to its silk properties, venom applications, and as a model for urban adaptation. We used long-read and Hi-C sequencing data, combined with transcriptomes, to assemble 14 chromosomes in a 1.46 Gb genome, with 38,393 genes annotated, and a BUSCO score of 95.3%. Our analyses identified high repetitive gene content and heterozygosity, consistent with other spider genomes, which has led to challenges in genome characterization. Our comparative evolutionary analyses of eight genomes available for species within the Araneoidea group (orb weavers and their descendants) identified 1,827 single-copy orthologs. Of these, 155 exhibit significant positive selection primarily associated with developmental genes, and with traits linked to sensory perception. These results support the hypothesis that several traits unique to spiders emerged from the adaptive evolution of ohnologs-or retained ancestrally duplicated genes-from ancient genome-wide duplication. These comparative spider genome analyses can serve as a model to understand how positive selection continually shapes ancestral duplications in generating novel traits today within and between diverse taxonomic groups.

The genetic adaptations of Toxoptera aurantii facilitated its rapid multiple plant hosts dispersal and invasion

Toxoptera aurantii Boyer de Fonscolombe (Hemiptera: Aphididae) can attack many plant hosts, including tea (Camellia sinensis L.), citrus (Citrus spp.), lychee (Litchi chinensis Sonn.), banana (Musa spp.), and pineapple (Ananas comasus L.) among others. It is a widely distributed hexapod and one of the most destructive pests in tea plantations, causing enormous economic losses in tea production each year. A high-quality reference genome is important to study the phylogenetics and evolution of T. aurantii because its genome is highly heterozygous and repetitive. We obtained a de novo genome assembly of T. aurantii at the chromosome level using a combination of long Nanopore reads from sequencing with high-throughput chromosome conformation capture technology. When finally assembled, the genome was 318.95 Mb on four chromosomes with a 15.19 Mb scaffold N50. A total of 12,162 genes encoded proteins, while there were 22.01% repetitive sequences that totaled 67.73 Mb. Phylogenetic analyses revealed that T. aurantii and Aphis gossypii parted ways approximately 7.6 million years ago (Mya). We used a combination of long-read single-molecule sequencing with Hi-C-based chromatin interaction maps that resulted in a reference chromosomal level reference genome of T. aurantii that was high quality. Our results will enable the exploration of the genetics behind the special biological features of T. aurantii and also provide a source of data that should be useful to compare the compare genome among the Hemiptera.

Mutations in Growth-Related Genes Induced by EMS Treatment in Scallops

Background: The goal of genetic breeding is to select variants with mutations that are related to expected traits, such as fast growth. Artificial induction has been widely used to obtain strains with more mutations for further selection. Ethylmethylsulfone (EMS) is one of the most commonly used chemical mutagens in plant and microorganism breeding. However, the application of EMS mutagenesis in shellfish has not been reported. The aim of this study is to evaluate the potential use of EMS as a mutagen in scallop breeding, especially in characterization of mutations in growth-related genes. Results: Our results indicated that hatching of about 50% of fertilized eggs was blocked by treatment with 20 mM EMS for 3 h and the resulted larvae developed normally into adult stages. We then evaluated the mutagenic effects of EMS by sequencing the genomes of 4 adult scallops from the control group and 12 from the treatment group at 8 months after fertilization. On average, after removing shared types of mutations, there were 1,151,380 ± 258,188 SNPs (Single Nucleotide Polymorphisms) and 229,256 ± 51,714 InDels (insertion-deletion) in each animal in the EMS treatment group, while there were only134841 ± 10,115 SNPs and 42,605 ± 5,136 InDels in the control group. The average mutation rate in the genome of the EMS treatment group (0.0137 ± 0.0013%) was about 9 times that of the control group (0.0015 ± 0.0002%). GO (Gene Ontology) annotation and KEGG (Kyoto Encyclopedia of Genes and Genomes) enrichment analyses revealed that mutations induced by EMS occurred evenly in most biological processes, cellular components and functions, as well in most pathways. However, significant lower percentage of mutations were found in the exonic region, in non-synonymous or Stopgain/Stoploss SNPs and in coding domains, suggesting apparent DNA repair or selection during grow-out stage. Analyses of the growth-related genes with mutations indicated that mutations in MFS (Major Facilitator Superfamily) and Tubulin were only found in the large-sized group (Five largest scallops: Treated-1, Treated-2, Treated-3, Treated-4, and Treated-5) and Homeobox and Socs (Suppressor of cytokine signaling) only in the small group (Two smallest scallops: Treated-11 and Treated-12). These results suggested that these genes may be involved in the regulation of growth in these animals, although further verification is certainly warranted. Conclusion: Treatment of fertilized eggs with 20 mM EMS for 3 h induced 9 times more mutations in scallop genomes. We found that mutations in MFS and Tubulin may be related to fast growth in the large-sized group and those mutations in Homeobox and SOCs may be involved in the slow growth in the small-sized scallops. EMS can be used to accelerate selection of economically important traits in molluscs.

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.

Homeotic Genes: Clustering, Modularity, and Diversity

Hox genes code for transcription factors and are evolutionarily conserved. They regulate a plethora of downstream targets to define the anterior-posterior (AP) body axis of a developing bilaterian embryo. Early work suggested a possible role of clustering and ordering of Hox to regulate their expression in a spatially restricted manner along the AP axis. However, the recent availability of many genome assemblies for different organisms uncovered several examples that defy this constraint. With recent advancements in genomics, the current review discusses the arrangement of Hox in various organisms. Further, we revisit their discovery and regulation in Drosophila melanogaster. We also review their regulation in different arthropods and vertebrates, with a significant focus on Hox expression in the crustacean Parahyale hawaiensis. It is noteworthy that subtle changes in the levels of Hox gene expression can contribute to the development of novel features in an organism. We, therefore, delve into the distinct regulation of these genes during primary axis formation, segment identity, and extra-embryonic roles such as in the formation of hair follicles or misregulation leading to cancer. Toward the end of each section, we emphasize the possibilities of several experiments involving various organisms, owing to the advancements in the field of genomics and CRISPR-based genome engineering. Overall, we present a holistic view of the functioning of Hox in the animal world.

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.

Homeobox Gene Duplication and Divergence in Arachnids

Homeobox genes are key toolkit genes that regulate the development of metazoans and changes in their regulation and copy number have contributed to the evolution of phenotypic diversity. We recently identified a whole genome duplication (WGD) event that occurred in an ancestor of spiders and scorpions (Arachnopulmonata), and that many homeobox genes, including two Hox clusters, appear to have been retained in arachnopulmonates. To better understand the consequences of this ancient WGD and the evolution of arachnid homeobox genes, we have characterized and compared the homeobox repertoires in a range of arachnids. We found that many families and clusters of these genes are duplicated in all studied arachnopulmonates (Parasteatoda tepidariorum, Pholcus phalangioides, Centruroides sculpturatus, and Mesobuthus martensii) compared with nonarachnopulmonate arachnids (Phalangium opilio, Neobisium carcinoides, Hesperochernes sp., and Ixodes scapularis). To assess divergence in the roles of homeobox ohnologs, we analyzed the expression of P. tepidariorum homeobox genes during embryogenesis and found pervasive changes in the level and timing of their expression. Furthermore, we compared the spatial expression of a subset of P. tepidariorum ohnologs with their single copy orthologs in P. opilio embryos. We found evidence for likely subfunctionlization and neofunctionalization of these genes in the spider. Overall our results show a high level of retention of homeobox genes in spiders and scorpions post-WGD, which is likely to have made a major contribution to their developmental evolution and diversification through pervasive subfunctionlization and neofunctionalization, and paralleling the outcomes of WGD in vertebrates.

In silico identification of the neuropeptidome of the pond wolf spider Pardosa pseudoannulata

Neuropeptides have been successfully documented in numerous arthropod species via in silico prediction from transcriptomic and genomic data. We recently sequenced the genome and nine transcriptomes of a chelicerate species, the pond wolf spider, Pardosa pseudoannulata. Here 43 neuropeptide families encoded by 87 neuropeptide genes were identified, among which 84 genes were presented with complete open reading frames. The neuropeptide genes often had paralogs and paralogous genes showed different expression profiles in nine transcriptomes. Six crustacean hyperglycemic hormone/ion transport peptide-like (CHH/ITP) genes were predicted and CHH/ITP6 was expressed much higher than the others. Orcokinin 1 and orcokinin 2 genes were both expressed in brain at a similar level. But, interestingly, orcokinin 1 gene was ubiquitously expressed in appendages while orcokinin 2 gene was enriched in venom gland to an extreme extent. The expression profiling of neuropeptide genes offers clues for further functional investigation. Paralogous genes were also found to be clustered at scaffolds such as nine insulin-like peptide genes at three scaffolds and six pyrokinin genes at two scaffolds, indicating a result of local gene duplication. In contrast, the four C-type allatostatin family members were scattered at five scaffolds, different from their closely associated locations reported in many arthropod species including several spiders. The comprehensive inventory of P. pseudoannulata neuropeptides here expands our repository of chelicerate neuropeptides and further promotes our understanding of neuropeptide evolution and functions.

Expression of NK cluster genes in the onychophoran Euperipatoides rowelli: implications for the evolution of NK family genes in nephrozoans

Background

Understanding the evolution and development of morphological traits of the last common bilaterian ancestor is a major goal of the evo-devo discipline. The reconstruction of this "urbilaterian" is mainly based on comparative studies of common molecular patterning mechanisms in recent model organisms. The NK homeobox genes are key players in many of these molecular pathways, including processes regulating mesoderm, heart and neural development. Shared features seen in the expression patterns of NK genes have been used to determine the ancestral bilaterian characters. However, the commonly used model organisms provide only a limited view on the evolution of these molecular pathways. To further investigate the ancestral roles of NK cluster genes, we analyzed their expression patterns in the onychophoran Euperipatoides rowelli.

Results

We identified nine transcripts of NK cluster genes in E. rowelli, including single copies of NK1, NK3, NK4, NK5, Msx, Lbx and Tlx, and two copies of NK6. All of these genes except for NK6.1 and NK6.2 are expressed in different mesodermal organs and tissues in embryos of E. rowelli, including the anlagen of somatic musculature and the heart. Furthermore, we found distinct expression patterns of NK3, NK5, NK6, Lbx and Msx in the developing nervous system. The same holds true for the NKL gene NK2.2, which does not belong to the NK cluster but is a related gene playing a role in neural patterning. Surprisingly, NK1, Msx and Lbx are additionally expressed in a segment polarity-like pattern early in development-a feature that has been otherwise reported only from annelids.

Conclusion

Our results indicate that the NK cluster genes were involved in mesoderm and neural development in the last common ancestor of bilaterians or at least nephrozoans (i.e., bilaterians to the exclusion of xenacoelomorphs). By comparing our data from an onychophoran to those from other bilaterians, we critically review the hypothesis of a complex "urbilaterian" with a segmented body, a pulsatile organ or heart, and a condensed mediolaterally patterned nerve cord.

Recent progress and prospects for advancing arachnid genomics

Arachnids exhibit tremendous species richness and adaptations of biomedical, industrial, and agricultural importance. Yet genomic resources for arachnids are limited, with the first few spider and scorpion genomes becoming accessible in the last four years. We review key insights from these genome projects, and recommend additional genomes for sequencing, emphasizing taxa of greatest value to the scientific community. We suggest greater sampling of spiders whose genomes are understudied but hold important protein recipes for silk and venom production. We further recommend arachnid genomes to address significant evolutionary topics, including the phenotypic impact of genome duplications. A barrier to high-quality arachnid genomes are assemblies based solely on short-read data, which may be overcome by long-range sequencing and other emerging methods.

The house spider genome reveals an ancient whole-genome duplication during arachnid evolution

BACKGROUND

The duplication of genes can occur through various mechanisms and is thought to make a major contribution to the evolutionary diversification of organisms. There is increasing evidence for a large-scale duplication of genes in some chelicerate lineages including two rounds of whole genome duplication (WGD) in horseshoe crabs. To investigate this further, we sequenced and analyzed the genome of the common house spider Parasteatoda tepidariorum.

RESULTS

We found pervasive duplication of both coding and non-coding genes in this spider, including two clusters of Hox genes. Analysis of synteny conservation across the P. tepidariorum genome suggests that there has been an ancient WGD in spiders. Comparison with the genomes of other chelicerates, including that of the newly sequenced bark scorpion Centruroides sculpturatus, suggests that this event occurred in the common ancestor of spiders and scorpions, and is probably independent of the WGDs in horseshoe crabs. Furthermore, characterization of the sequence and expression of the Hox paralogs in P. tepidariorum suggests that many have been subject to neo-functionalization and/or sub-functionalization since their duplication.

CONCLUSIONS

Our results reveal that spiders and scorpions are likely the descendants of a polyploid ancestor that lived more than 450 MYA. Given the extensive morphological diversity and ecological adaptations found among these animals, rivaling those of vertebrates, our study of the ancient WGD event in Arachnopulmonata provides a new comparative platform to explore common and divergent evolutionary outcomes of polyploidization events across eukaryotes.

The house spider genome reveals an ancient whole-genome duplication during arachnid evolution

Background

The duplication of genes can occur through various mechanisms and is thought to make a major contribution to the evolutionary diversification of organisms. There is increasing evidence for a large-scale duplication of genes in some chelicerate lineages including two rounds of whole genome duplication (WGD) in horseshoe crabs. To investigate this further, we sequenced and analyzed the genome of the common house spider Parasteatoda tepidariorum.

Results

We found pervasive duplication of both coding and non-coding genes in this spider, including two clusters of Hox genes. Analysis of synteny conservation across the P. tepidariorum genome suggests that there has been an ancient WGD in spiders. Comparison with the genomes of other chelicerates, including that of the newly sequenced bark scorpion Centruroides sculpturatus, suggests that this event occurred in the common ancestor of spiders and scorpions, and is probably independent of the WGDs in horseshoe crabs. Furthermore, characterization of the sequence and expression of the Hox paralogs in P. tepidariorum suggests that many have been subject to neo-functionalization and/or sub-functionalization since their duplication.

Conclusions

Our results reveal that spiders and scorpions are likely the descendants of a polyploid ancestor that lived more than 450 MYA. Given the extensive morphological diversity and ecological adaptations found among these animals, rivaling those of vertebrates, our study of the ancient WGD event in Arachnopulmonata provides a new comparative platform to explore common and divergent evolutionary outcomes of polyploidization events across eukaryotes.

Electronic supplementary material

The online version of this article (doi:10.1186/s12915-017-0399-x) contains supplementary material, which is available to authorized users.

Expression and function of spineless orthologs correlate with distal deutocerebral appendage morphology across Arthropoda

The deutocerebral (second) head segment is putatively homologous across Arthropoda, in spite of remarkable disparity of form and function of deutocerebral appendages. In Mandibulata this segment bears a pair of sensory antennae, whereas in Chelicerata the same segment bears a pair of feeding appendages called chelicerae. Part of the evidence for the homology of deutocerebral appendages is the conserved function of homothorax (hth), which has been shown to specify antennal or cheliceral fate in the absence of Hox signaling, in both mandibulate and chelicerate exemplars. However, the genetic basis for the morphological disparity of antenna and chelicera is not understood. To test whether downstream targets of hth have diverged in a lineage-specific manner, we examined the evolution of the function and expression of spineless (ss), which in two holometabolous insects is known to act as a hth target and distal antennal determinant. Toward expanding phylogenetic representation of gene expression data, here we show that strong expression of ss is observed in developing antennae of a hemimetabolous insect, a centipede, and an amphipod crustacean. By contrast, ss orthologs are not expressed throughout the cheliceral limb buds of spiders or harvestmen during developmental stages when appendage fate is specified. RNA interference-mediated knockdown of ss in Oncopeltus fasciatus, which bears a simple plesiomorphic antenna, resulted in homeotic distal antenna-to-leg transformation, comparable to data from holometabolous insect counterparts. Knockdown of hth in Oncopeltus fasciatus abrogated ss expression, suggesting conservation of upstream regulation. These data suggest that ss may be a flagellar (distal antennal) determinant more broadly, and that this function was acquired at the base of Mandibulata.

Adaptive evolution of insect selective excitatory β-type sodium channel neurotoxins from scorpion venom

Insect selective excitatory β-type sodium channel neurotoxins from scorpion venom (β-NaScTxs) are composed of about 70-76 amino acid residues and share a common scaffold stabilized by four unique disulfide bonds. The phylogenetic analysis of these toxins was hindered by limited sequence data. In our recent study, two new insect selective excitatory β-NaScTxs, LmIT and ImIT, were isolated from Lychas mucronatus and Isometrus maculatus, respectively. With the sequences previously reported, we examined the adaptive molecular evolution of insect selective excitatory β-NaScTxs by estimating the nonsynonymous-to-synonymous rate ratio (ω=d_N/d_S). The results revealed 12 positively selected sites in the genes of insect selective excitatory β-NaScTxs. Moreover, these positively selected sites match well with the sites important for interacting with sodium channels, as demonstrated in previous mutagenesis study. These results reveal that adaptive evolution after gene duplication is one of the most important genetic mechanisms of scorpion neurotoxin diversification.

Segmentation and tagmosis in Chelicerata

Patterns of segmentation and tagmosis are reviewed for Chelicerata. Depending on the outgroup, chelicerate origins are either among taxa with an anterior tagma of six somites, or taxa in which the appendages of somite I became increasingly raptorial. All Chelicerata have appendage I as a chelate or clasp-knife chelicera. The basic trend has obviously been to consolidate food-gathering and walking limbs as a prosoma and respiratory appendages on the opisthosoma. However, the boundary of the prosoma is debatable in that some taxa have functionally incorporated somite VII and/or its appendages into the prosoma. Euchelicerata can be defined on having plate-like opisthosomal appendages, further modified within Arachnida. Total somite counts for Chelicerata range from a maximum of nineteen in groups like Scorpiones and the extinct Eurypterida down to seven in modern Pycnogonida. Mites may also show reduced somite counts, but reconstructing segmentation in these animals remains challenging. Several innovations relating to tagmosis or the appendages borne on particular somites are summarised here as putative apomorphies of individual higher taxa. We also present our observations within the concept of pseudotagma, whereby the true tagmata - the prosoma and opisthosoma - can be defined on a fundamental change in the limb series while pseudotagmata, such as the cephalosoma/proterosoma, are expressed as divisions in sclerites covering the body without an accompanying change in the appendages.

Transcriptome analysis on the exoskeleton formation in early developmetal stages and reconstruction scenario in growth-moulting in Litopenaeus vannamei

Exoskeleton construction is an important issue in shrimp. To better understand the molecular mechanism of exoskeleton formation, development and reconstruction, the transcriptome of the entire developmental process in Litopenaeus vannamei, including nine early developmental stages and eight adult-moulting stages, was sequenced and analysed using Illumina RNA-seq technology. A total of 117,539 unigenes were obtained, with 41.2% unigenes predicting the full-length coding sequence. Gene Ontology, Clusters of Orthologous Group (COG), the Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis and functional annotation of all unigenes gave a better understanding of the exoskeleton developmental process in L. vannamei. As a result, more than six hundred unigenes related to exoskeleton development were identified both in the early developmental stages and adult-moulting. A cascade of sequential expression events of exoskeleton-related genes were summarized, including exoskeleton formation, regulation, synthesis, degradation, mineral absorption/reabsorption, calcification and hardening. This new insight on major transcriptional events provide a deep understanding for exoskeleton formation and reconstruction in L. vannamei. In conclusion, this is the first study that characterized the integrated transcriptomic profiles cover the entire exoskeleton development from zygote to adult-moulting in a crustacean, and these findings will serve as significant references for exoskeleton developmental biology and aquaculture research.

Neuropeptide Evolution: Chelicerate Neurohormone and Neuropeptide Genes may reflect one or more whole genome duplications

Four genomes and two transcriptomes from six Chelicerate species were analyzed for the presence of neuropeptide and neurohormone precursors and their GPCRs. The genome from the spider Stegodyphus mimosarum yielded 87 neuropeptide precursors and 101 neuropeptide GPCRs. High neuropeptide transcripts were also found in the trancriptomes of three other spiders, Latrodectus hesperus, Parasteatoda tepidariorum and Acanthoscurria geniculata. For the scorpion Mesobuthus martensii the numbers are 79 and 74 respectively. The very small genome of the house dust mite, Dermatophagoides farinae, on the other hand contains a much smaller number of such genes. A few new putative Arthropod neuropeptide genes were discovered. Thus, both spiders and the scorpion have an achatin gene and in spiders there are two different genes encoding myosuppressin-like peptides while spiders also have two genes encoding novel LGamides. Another finding is the presence of trissin in spiders and scorpions, while neuropeptide genes that seem to be orthologs of Lottia LFRYamide and Platynereis CCRFamide were also found. Such genes were also found in various insect species, but seem to be lacking from the Holometabola. The Chelicerate neuropeptide and neuropeptide GPCR genes often have paralogs. As the large majority of these are probably not due to local gene duplications, is not impossible that they reflect the effects of one or more ancient whole genome duplications.

Pervasive microRNA Duplication in Chelicerates: Insights from the Embryonic microRNA Repertoire of the Spider Parasteatoda tepidariorum

MicroRNAs are small (∼22 nt) noncoding RNAs that repress translation and therefore regulate the production of proteins from specific target mRNAs. microRNAs have been found to function in diverse aspects of gene regulation within animal development and many other processes. Among invertebrates, both conserved and novel, lineage specific, microRNAs have been extensively studied predominantly in holometabolous insects such as Drosophila melanogaster However little is known about microRNA repertoires in other arthropod lineages such as the chelicerates. To understand the evolution of microRNAs in this poorly sampled subphylum, we characterized the microRNA repertoire expressed during embryogenesis of the common house spider Parasteatoda tepidariorum We identified a total of 148 microRNAs in P. tepidariorum representing 66 families. Approximately half of these microRNA families are conserved in other metazoans, while the remainder are specific to this spider. Of the 35 conserved microRNAs families 15 had at least two copies in the P. tepidariorum genome. A BLAST-based approach revealed a similar pattern of duplication in other spiders and a scorpion, but not among other chelicerates and arthropods, with the exception of a horseshoe crab. Among the duplicated microRNAs we found examples of lineage-specific tandem duplications, and the duplication of entire microRNA clusters in three spiders, a scorpion, and in a horseshoe crab. Furthermore, we found that paralogs of many P. tepidariorum microRNA families exhibit arm switching, which suggests that duplication was often followed by sub- or neofunctionalization. Our work shows that understanding the evolution of microRNAs in the chelicerates has great potential to provide insights into the process of microRNA duplication and divergence and the evolution of animal development.

Evidence of duplicated Hox genes in the most recent common ancestor of extant scorpions

Scorpions (order Scorpiones) are unusual among arthropods, both for the extreme heteronomy of their bauplan and for the high gene family turnover exhibited in their genomes. These phenomena appear to be correlated, as two scorpion species have been shown to possess nearly twice the number of Hox genes present in most arthropods. Segmentally offset anterior expression boundaries of a subset of Hox paralogs have been shown to correspond to transitions in segmental identities in the scorpion posterior tagmata, suggesting that posterior heteronomy in scorpions may have been achieved by neofunctionalization of Hox paralogs. However, both the first scorpion genome sequenced and the developmental genetic data are based on exemplars of Buthidae, one of 19 families of scorpions. It is therefore not known whether Hox paralogy is limited to Buthidae or widespread among scorpions. We surveyed 24 high throughput transcriptomes and the single whole genome available for scorpions, in order to test the prediction that Hox gene duplications are common to the order. We used gene tree parsimony to infer whether the paralogy was consistent with a duplication event in the scorpion common ancestor. Here we show that duplicated Hox genes in non-buthid scorpions occur in six of the ten Hox classes. Gene tree topologies and parsimony-based reconciliation of the gene trees are consistent with a duplication event in the most recent common ancestor of scorpions. These results suggest that a Hox paralogy, and by extension the model of posterior patterning established in a buthid, can be extended to non-Buthidae scorpions.

Genome Sequencing of the Phytoseiid Predatory Mite Metaseiulus occidentalis Reveals Completely Atomized Hox Genes and Superdynamic Intron Evolution

Metaseiulus occidentalis is an eyeless phytoseiid predatory mite employed for the biological control of agricultural pests including spider mites. Despite appearances, these predator and prey mites are separated by some 400 Myr of evolution and radically different lifestyles. We present a 152-Mb draft assembly of the M. occidentalis genome: Larger than that of its favored prey, Tetranychus urticae, but considerably smaller than those of many other chelicerates, enabling an extremely contiguous and complete assembly to be built-the best arachnid to date. Aided by transcriptome data, genome annotation cataloged 18,338 protein-coding genes and identified large numbers of Helitron transposable elements. Comparisons with other arthropods revealed a particularly dynamic and turbulent genomic evolutionary history. Its genes exhibit elevated molecular evolution, with strikingly high numbers of intron gains and losses, in stark contrast to the deer tick Ixodes scapularis Uniquely among examined arthropods, this predatory mite's Hox genes are completely atomized, dispersed across the genome, and it encodes five copies of the normally single-copy RNA processing Dicer-2 gene. Examining gene families linked to characteristic biological traits of this tiny predator provides initial insights into processes of sex determination, development, immune defense, and how it detects, disables, and digests its prey. As the first reference genome for the Phytoseiidae, and for any species with the rare sex determination system of parahaploidy, the genome of the western orchard predatory mite improves genomic sampling of chelicerates and provides invaluable new resources for functional genomic analyses of this family of agriculturally important mites.

Composition and genomic organization of arthropod Hox clusters

Background

The ancestral arthropod is believed to have had a clustered arrangement of ten Hox genes. Within arthropods, Hox gene mutations result in transformation of segment identities. Despite the fact that variation in segment number/character was common in the diversification of arthropods, few examples of Hox gene gains/losses have been correlated with morphological evolution. Furthermore, a full appreciation of the variation in the genomic arrangement of Hox genes in extant arthropods has not been recognized, as genome sequences from each major arthropod clade have not been reported until recently. Initial genomic analysis of the chelicerate Tetranychusurticae suggested that loss of Hox genes and Hox gene clustering might be more common than previously assumed. To further characterize the genomic evolution of arthropod Hox genes, we compared the genomic arrangement and general characteristics of Hox genes from representative taxa from each arthropod subphylum.

Results

In agreement with others, we find arthropods generally contain ten Hox genes arranged in a common orientation in the genome, with an increasing number of sampled species missing either Hox3 or abdominal-A orthologs. The genomic clustering of Hox genes in species we surveyed varies significantly, ranging from 0.3 to 13.6 Mb. In all species sampled, arthropod Hox genes are dispersed in the genome relative to the vertebrate Mus musculus. Differences in Hox cluster size arise from variation in the number of intervening genes, intergenic spacing, and the size of introns and UTRs. In the arthropods surveyed, Hox gene duplications are rare and four microRNAs are, in general, conserved in similar genomic positions relative to the Hox genes.

Conclusions

The tightly clustered Hox complexes found in the vertebrates are not evident within arthropods, and differential patterns of Hox gene dispersion are found throughout the arthropods. The comparative genomic data continue to support an ancestral arthropod Hox cluster of ten genes with a shared orientation, with four Hox gene-associated miRNAs, although the degree of dispersion between genes in an ancestral cluster remains uncertain. Hox3 and abdominal-A orthologs have been lost in multiple, independent lineages, and current data support a model in which inversions of the Abdominal-B locus that result in the loss of abdominal-A correlate with reduced trunk segmentation.

Electronic supplementary material

The online version of this article (doi:10.1186/s13227-016-0048-4) contains supplementary material, which is available to authorized users.

Allatostatins C, double C and triple C, the result of a local gene triplication in an ancestral arthropod

Allatostatin C is the arthropod homolog of vertebrate somatostatin. The gene went through a local gene triplification leading to the existence of three genes coding such peptides, allatostatins C, CC and CCC. All three genes are still present in several chelicerates, such as the horseshoe crab Limulus polyphemus, several spiders and the scorpion Mesobuthus martensii, the myriapod Strigamia maritima, as well as at least two insect species, Locusta migratoria and Athalia rosae, a sawfly. All three peptides have well conserved primary structures and peptides can easily be classified as either allatostatin C, CC or CCC. In most insect species only two of the genes have been preserved. In many species, these are CC and CCC, but in Diptera, Lepidoptera and Coleoptera it are allatostatins C and CC that are still present. In some arthropod species two or even all three genes can still be found closely associated in the genome and are present on the same scaffold showing that a local amplification was at the origin of these genes.

Neuropeptide evolution: Chelicerate neurohormone and neuropeptide genes may reflect one or more whole genome duplications

Four genomes and two transcriptomes from six Chelicerate species were analyzed for the presence of neuropeptide and neurohormone precursors and their GPCRs. The genome from the spider Stegodyphus mimosarum yielded 87 neuropeptide precursors and 120 neuropeptide GPCRs. Many neuropeptide transcripts were also found in the transcriptomes of three other spiders, Latrodectus hesperus, Parasteatoda tepidariorum and Acanthoscurria geniculata. For the scorpion Mesobuthus martensii the numbers are 79 and 93 respectively. The very small genome of the house dust mite, Dermatophagoides farinae, on the other hand contains a much smaller number of such genes. A few new putative Arthropod neuropeptide genes were discovered. Thus, both spiders and the scorpion have an achatin gene and in spiders there are two different genes encoding myosuppressin-like peptides while spiders also have two genes encoding novel LGamides. Another finding is the presence of trissin in spiders and scorpions, while neuropeptide genes that seem to be orthologs of Lottia LFRYamide and Platynereis CCRFamide were also found. Such genes were also found in various insect species, but seem to be lacking from the Holometabola. The Chelicerate neuropeptide and neuropeptide GPCR genes often have paralogs. As the large majority of these are probably not due to local gene duplications, is plausible that they reflect the effects of one or more ancient whole genome duplications.