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.

Rapid specialization of counter defenses enables two-spotted spider mite to adapt to novel plant hosts.

Genetic adaptation, occurring over a long evolutionary time, enables host-specialized herbivores to develop novel resistance traits and to efficiently counteract the defenses of a narrow range of host plants. In contrast, physiological acclimation, leading to the suppression and/or detoxification of host defenses, is hypothesized to enable broad generalists to shift between plant hosts. However, the host adaptation mechanisms used by generalists composed of host-adapted populations are not known. Two-spotted spider mite (TSSM; Tetranychus urticae) is an extreme generalist herbivore whose individual populations perform well only on a subset of potential hosts. We combined experimental evolution, Arabidopsis thaliana genetics, mite reverse genetics, and pharmacological approaches to examine mite host adaptation upon the shift of a bean (Phaseolus vulgaris)-adapted population to Arabidopsis. We showed that cytochrome P450 monooxygenases are required for mite adaptation to Arabidopsis. We identified activities of two tiers of P450s: general xenobiotic-responsive P450s that have a limited contribution to mite adaptation to Arabidopsis and adaptation-associated P450s that efficiently counteract Arabidopsis defenses. In approximately 25 generations of mite selection on Arabidopsis plants, mites evolved highly efficient detoxification-based adaptation, characteristic of specialist herbivores. This demonstrates that specialization to plant resistance traits can occur within the ecological timescale, enabling the TSSM to shift to novel plant hosts.

Past climate cooling promoted global dispersal of amphipods from Tian Shan montane lakes to circumboreal lakes.

Climate changes have substantial impacts on the geographic distribution of montane lakes and evolutionary dynamics of cold-adapted species. Past climate cooling is hypothesized to have promoted the dispersal of cold-adapted species via montane lakes, while future climate warming is thought to constrain their distributions. We test this hypothesis by using phylogeographic analysis and niche modeling of the Holarctic crustacean Gammarus lacustris with global sampling comprised of 567 sequenced individuals and 3180 occurrence records. We found that the species arose in Tian Shan in Central Asia and dispersed into montane lakes along the Alps, Himalayas, Tibet, East Asia, and the North American Rocky Mountain ranges, with accelerated diversification rates outside Tian Shan. Climatically suitable regions for geographic lineages of G. lacustris were larger during cooling periods (LGM), but smaller during warming periods (Mid-Holocene). In the future (2070) scenario, potential distributions in the Himalayas, North Tibet, South Tibet and North America are predicted to expand, whereas ranges in East Asia, Europe and Tian Shan will decline. Our results suggest that Mid-Miocene-to-Pleistocene continuous cooling promoted multiple independent dispersal events out of Tian Shan due to increased availability of montane lakes via "budding" of lineages. Montane lakes are conduits through which cold-adapted amphipods globally dispersed, dominating circumboreal lakes. However, future climate warming is likely to force organisms to shift upward in altitude and northward in latitude, leading to a future change in local populations. These findings highlight the importance of conservation of montane lakes, especially in the context of climate change.

Geography alone cannot explain Tetranychus truncatus (Acari: Tetranychidae) population abundance and genetic diversity in the context of the center-periphery hypothesis.

The center-periphery hypothesis (CPH) states that the genetic diversity, genetic flow, and population abundance of a species are highest at the center of the species' geographic distribution. However, most CPH studies have focused on the geographic distance and have ignored ecological and historical effects. Studies using niche models to define the center and periphery of a distribution and the interactions among geographical, ecological, and historical gradients have rarely been done in the framework of the CPH, especially in biogeographical studies of animal species. Here, we examined the CPH for a widely distributed arthropod, Tetranychus truncatus (Acari: Tetranychidae), in eastern China using three measurements: geographic distance to the center of the distribution (geography), ecological suitability based on current climate data (ecology), and historical climate data from the last glacial maximum (history). We found that the relative abundances of different populations were more strongly related to ecology than to geography and history. Genetic diversity within populations and genetic differentiation among populations based on mitochondrial marker were only significantly related to history. However, the genetic diversity and population differentiation based on microsatellites were significantly related to all three CPH measurements. Overall, population abundance and genetic pattern cannot be explained very well by geography alone. Our results show that ecological gradients explain the variation in population abundance better than geographic gradients and historical factors, and that current and historical factors strongly influence the spatial patterns of genetic variation. This study highlights the importance of examining more than just geography when assessing the CPH.

Phylogenetic signals in pest abundance and distribution range of spider mites.

BACKGROUND: Attributes of pest species like host range are frequently reported as being evolutionarily constrained and showing phylogenetic signal. Because these attributes in turn could influence the abundance and impact of species, phylogenetic information could be useful in predicting the likely status of pests. In this study, we used regional (China) and global datasets to investigate phylogenetic patterns in occurrence patterns and host ranges of spider mites, which constitute a pest group of many cropping systems worldwide. RESULTS: We found significant phylogenetic signal in relative abundance and distribution range both at the regional and global scales. Relative abundance and range size of spider mites were positively correlated with host range, although these correlations became weaker after controlling for phylogeny. CONCLUSIONS: The results suggest that pest impacts are evolutionarily constrained. Information that is easily obtainable - including the number of known hosts and phylogenetic position of the mites - could therefore be useful in predicting future pest risk of species.

High genetic diversity in a 'recent outbreak' spider mite, Tetranychus pueraricola, in mainland China.

Tetranychus pueraricola is a newly reported spider mite that occurs frequently in mainland China. It is possible that this species was introduced from elsewhere and became a serious pest recently. However, the correct identification of red-pigmented spider mites has repeatedly proven problematic. There is also the possibility that T. pueraricola in China was long misidentified as its sibling species, Tetranychus urticae (red form). To test which of these two scenarios is the more likely, individuals from 14 populations of T. pueraricola and five populations of T. urticae (red form) in China were sampled and genotyped using mitochondrial COI and microsatellite loci. Unlike a recent invasive species, the genetic diversity of T. pueraricola was very high with high mitochondrial genetic diversity (16 haplotypes), high effective alleles (Ne = 2.038 ± 0.081) and expected heterozygosity (He = 0.395 ± 0.016). Surprisingly, we found that all T. urticae (red form) populations shared only one mitochondrial haplotype and showed quite low genetic diversity (Ne = 1.443 ± 0.055; He = 0.234 ± 0.025) which was even lower than that of the green form of T. urticae in mainland China from a previous study. We did not detect significant signals of recent bottlenecks for most populations from both species. These results suggest T. pueraricola is unlikely to be a recent invasive pest but a species that has existed in China for a long time. It is probable that T. pueraricola in China has long been misidentified as T. urticae (red form).

Comparative transcriptomes and reciprocal best hit analysis revealed potential pigment genes in two color forms of Tetranychus urticae.

Tetranychus urticae Koch is a worldwide agricultural pest. There are two color forms: red and green. The molecular mechanism underlying this color variation is unknown. To elucidate the mechanism, we characterized differentially expressed pigment pathway genes shared in the transcriptomes of these two forms using RNA sequencing and reciprocal best hit analysis. Differentially expressed pigment pathway genes were determined by qRT-PCR to confirm the accuracy of RNA-Seq. The transcriptomes revealed 963 differentially expressed genes (DEGs), of which 687 DEGs were higher in the green form. KEGG enrichment analysis revealed carotenoid biosynthesis genes in T. urticae. Reciprocal best hit analysis revealed 817 putative pigment pathway genes, 38 of which were differentially expressed and mainly classified into four categories: heme, melanin, ommochrome and rhodopsin. Phylogenetic analysis of homologous ommochrome genes showed that tetur09g01950 is closely related to Ok. This study revealed putative pigment pathway genes in the two forms of T. urticae, and might provide a new resource for understanding the mechanism of color variation.

Development of microsatellite markers for six Tetranychus species by transfer from Tetranychus urticae genome.

Microsatellite markers are frequently used to explore the population genetic structure of organisms. Spider mites (genus Tetranychus) are important agricultural pests. Several markers have been developed for T. urticae, but for other spider mites, few such markers are available, hampering studies of their population genetics. In this study, we developed and characterized microsatellite markers for six non-model spider mite species (T. truncatus, T. kanzawai, T. ludeni, T. piercei, T. phaselus and T. pueraricola) by cross-species amplification of markers in the T. urticae genome, in order to better understand the population structure of Tetranychus species. Among 228 screened loci, many were polymorphic, including 13 loci in T. urticae, 11 loci in T. truncatus, 15 loci in T. pueraricola, 23 loci in T. kanzawai, 19 loci in T. piercei, 11 loci in T. phaselus and 9 loci in T. ludeni. Sequence analysis determined that the fragment length variations of the transferred microsatellites were mainly due to the variations of the numbers of repeats. These new microsatellite markers should be useful for studying the population genetics of the seven Tetranychus species.

Single-cell transcriptomics reveals the brain evolution of web-building spiders.

Spiders are renowned for their efficient capture of flying insects using intricate aerial webs. How the spider nervous systems evolved to cope with this specialized hunting strategy and various environmental clues in an aerial space remains unknown. Here we report a brain-cell atlas of >30,000 single-cell transcriptomes from a web-building spider (Hylyphantes graminicola). Our analysis revealed the preservation of ancestral neuron types in spiders, including the potential coexistence of noradrenergic and octopaminergic neurons, and many peptidergic neuronal types that are lost in insects. By comparing the genome of two newly sequenced plesiomorphic burrowing spiders with three aerial web-building spiders, we found that the positively selected genes in the ancestral branch of web-building spiders were preferentially expressed (42%) in the brain, especially in the three mushroom body-like neuronal types. By gene enrichment analysis and RNAi experiments, these genes were suggested to be involved in the learning and memory pathway and may influence the spiders' web-building and hunting behaviour. Our results provide key sources for understanding the evolution of behaviour in spiders and reveal how molecular evolution drives neuron innovation and the diversification of associated complex behaviours.

Marine-montane transitions coupled with gill and genetic convergence in extant crustacean.

Marine-terrestrial transition represents an important aspect of organismal evolution that requires numerous morphological and genetic innovations and has been hypothesized to be caused by geological changes. We used talitrid crustaceans with marine-coastal-montane extant species at a global scale to investigate the marine origination and terrestrial adaptation. Using genomic data, we demonstrated that marine ancestors repeatedly colonized montane terrestrial habitats during the Oligocene to Miocene. Biological transitions were well correlated with plate collisions or volcanic island formation, and top-down cladogenesis was observed on the basis of a positive relationship between ancestral habitat elevation and divergence time for montane lineages. We detected convergent variations of convoluted gills and convergent evolution of SMC3 associated with montane transitions. Moreover, using CRISPR-Cas9 mutagenesis, we proposed that SMC3 potentially regulates the development of exites, such as talitrid gills. Our results provide a living model for understanding biological innovations and related genetic regulatory mechanisms associated with marine-terrestrial transitions.

Chromosomal-level genome of a sheet-web spider provides insight into the composition and evolution of venom.

Spiders are the most abundant venomous predators in the world. Previous research related to spider venom has mostly relied on transcriptomes and proteomes, with only a few high-quality genomes available. This is far from consistent with the species diversity of spiders. In this study, we constructed a high-quality chromosome-level genome assembly of Hylyphantes graminicola, which contained 13 chromosomes, with a genome length of 931.68 Mb and scaffold N50 of 77.07 Mb. Integrating genome, transcriptome, and proteome profiling, we identified a total of 59 coding genes among nine toxin gene families. Among them, Group 7 allergen (ALL7) protein was reported in spider venom for the first time. Its coding genes had a predicted signal peptide and maintained high expression levels in the venom, suggesting that ALL7 plays an important role in venom and maybe is a type of newly discovered venom toxin in the spider. By implementing comparative genomics, we found a similar gene number of main toxin gene families in spiders and the scorpion genome with conservative evolutionary rates, indicating that these toxin genes could be an ancient (~400 million years) and a conserved "basic toolkit" for spiders and scorpions to perform primary defence functions. Obtaining high-quality chromosome-level genomes from spiders not only facilitates venom research and toxin resource application, but also can improve comparative genomic analysis in other important traits, like the evolution of silk or behaviour.

Population genomic data in spider mites point to a role for local adaptation in shaping range shifts.

Local adaptation is particularly likely in invertebrate pests that typically have short generation times and large population sizes, but there are few studies on pest species investigating local adaptation and separating this process from contemporaneous and historical gene flow. Here, we use a population genomic approach to investigate evolutionary processes in the two most dominant spider mites in China, Tetranychus truncatus Ehara and Tetranychus pueraricola Ehara et Gotoh, which have wide distributions, short generation times, and large population sizes. We generated genome resequencing of 246 spider mites mostly from China, as well as Japan and Canada at a combined total depth of 3,133×. Based on demographic reconstruction, we found that both mite species likely originated from refugia in southwestern China and then spread to other regions, with the dominant T. truncatus spreading ~3,000 years later than T. pueraricola. Estimated changes in population sizes of the pests matched known periods of glaciation and reinforce the recent expansion of the dominant spider mites. T. truncatus showed a greater extent of local adaptation with more genes (76 vs. 17) associated with precipitation, including candidates involved in regulation of homeostasis of water and ions, signal transduction, and motor skills. In both species, many genes (135 in T. truncatus and 95 in T. pueraricola) also showed signatures of selection related to elevation, including G-protein-coupled receptors, cytochrome P450s, and ABC-transporters. Our results point to historical expansion processes and climatic adaptation in these pests which could have contributed to their growing importance, particularly in the case of T. truncatus.

Environmental RNA interference in two-spotted spider mite, Tetranychus urticae, reveals dsRNA processing requirements for efficient RNAi response.

Comprehensive understanding of pleiotropic roles of RNAi machinery highlighted the conserved chromosomal functions of RNA interference. The consequences of the evolutionary variation in the core RNAi pathway genes are mostly unknown, but may lead to the species-specific functions associated with gene silencing. The two-spotted spider mite, Tetranychus urticae, is a major polyphagous chelicerate pest capable of feeding on over 1100 plant species and developing resistance to pesticides used for its control. A well annotated genome, susceptibility to RNAi and economic importance, make T. urticae an excellent candidate for development of an RNAi protocol that enables high-throughput genetic screens and RNAi-based pest control. Here, we show that the length of the exogenous dsRNA critically determines its processivity and ability to induce RNAi in vivo. A combination of the long dsRNAs and the use of dye to trace the ingestion of dsRNA enabled the identification of genes involved in membrane transport and 26S proteasome degradation as sensitive RNAi targets. Our data demonstrate that environmental RNAi can be an efficient reverse genetics and pest control tool in T. urticae. In addition, the species-specific properties together with the variation in the components of the RNAi machinery make T. urticae a potent experimental system to study the evolution of RNAi pathways.

A change in the bacterial community of spider mites decreases fecundity on multiple host plants.

Bacterial symbionts may influence the fitness of their herbivore hosts, but such effects have been poorly studied across most invertebrate groups. The spider mite, Tetranychus truncatus, is a polyphagous agricultural pest harboring various bacterial symbionts whose function is largely unknown. Here, by using a high-throughput 16S rRNA amplicon sequencing approach, we characterized the bacterial diversity and community composition of spider mites fed on five host plants after communities were modified following tetracycline exposure. We demonstrated that spider mite bacterial diversity and community composition were significantly affected by host plants and antibiotics. In particular, the abundance of the maternally inherited endosymbionts Wolbachia and Spiroplasma significantly differed among spider mites that were reared on different plant species and were completely removed by antibiotics. There was an overall tendency for daily fecundity to be lower in the mites with reduced bacterial diversity following the antibiotic treatment. Our data suggest that host plants and antibiotics can shape spider mite bacterial communities and that bacterial symbionts improve mite performance.

RNAi-based reverse genetics in the chelicerate model Tetranychus urticae: A comparative analysis of five methods for gene silencing.

RNA interference (RNAi) can be used for the protection against agricultural pests through the silencing of genes required for pest fitness. To assess the potential of RNAi approaches in the two-spotted spider mite, Tetranychus urticae, we compared 5 methods for the delivery of double-stranded RNA (dsRNA). These methods include mite feeding on either (i) leaves floating on a dsRNA solution, (ii) dsRNA-expressing plants, (iii) artificial diet supplemented with dsRNA, or (iv) dsRNA-coated leaves, and (v) mite soaking in a dsRNA solution. In all cases, the gene targeted for method validation was the Vacuolar-type H+-ATPase (TuVATPase), encoding a constitutively expressed ATP-driven proton pump located in the membrane. Down-regulation of TuVATPase increased mortality and/or reduced fecundity in all methods, but with variable efficiency. The most efficient methods for dsRNA delivery were direct soaking of mites in the dsRNA solution and mite feeding on dsRNA-coated leaves that mimics dsRNA application as a sprayable pesticide. Both resulted in a dark-body phenotype not observed in mites treated with a control dsRNA. Although with lower efficiency, dsRNA designed for TuVATPase silencing and expressed in transgenic Arabidopsis plants impacted the fitness of mites feeding on these plants. RNAi may thus be a valuable strategy to control spider mite populations, either as a sprayable pesticide or through transgenic crops. This comparative methodological study focusing on the induction of RNAi-based gene silencing in T. urticae paves the way for reverse genetics approaches in this model chelicerate system and prepares large-scale systematic RNAi screens as a first step towards the development of specific RNA-based pesticides. Such alternative molecules may help control spider mites that cause significant damages to crops and ornamental plant species, as well as other chelicerates detrimental to agriculture and health.

The complete mitochondrial genome of Tetranychus truncatus Ehara (Acari: Tetranychidae).

The complete mitochondrial genome of Tetranychus truncatus Ehara (Acari: Tetranychidae) is a typical circular DNA with length of 13,089 bp (GenBank accession number: KM111296). The genome contains all 13 protein-coding genes (PCGs), an A + T-rich region, two rRNA genes and 22 tRNA genes. The A + T content of the mitochondrial genome is 84.5%. The AT-skew is positive (0.032) while the GC-skew is negative (-0.058). The gene arrangement is conserved in T. urticae, Panonychus citri and P. ulmi which are in the same family (Tetranychidae). The A + T-rich region is only 43 bp in length with high A + T content (97.7%). All the PCGs start with typical ATD codons. Eight PCGs have complete TAA stop codons, while five PCGs have an incomplete stop codon (T).

The complete mitochondrial genomes of six species of Tetranychus provide insights into the phylogeny and evolution of spider mites.

Many spider mites belonging to the genus Tetranychus are of agronomical importance. With limited morphological characters, Tetranychus mites are usually identified by a combination of morphological characteristics and molecular diagnostics. To clarify their molecular evolution and phylogeny, the mitochondrial genomes of the green and red forms of Tetranychus urticae as well as T. kanzawai, T. ludeni, T. malaysiensis, T. phaselus, T. pueraricola were sequenced and compared. The seven mitochondrial genomes are typical circular molecules of about 13,000 bp encoding and they are composed of the complete set of 37 genes that are usually found in metazoans. The order of the mitochondrial (mt) genes is the same as that in the mt genomes of Panonychus citri and P. ulmi, but very different from that in other Acari. The J-strands of the mitochondrial genomes have high (∼ 84%) A+T contents, negative GC-skews and positive AT-skews. The nucleotide sequence of the cox1 gene, which is commonly used as a taxon barcode and molecular marker, is more highly conserved than the nucleotide sequences of other mitochondrial genes in these seven species. Most tRNA genes in the seven genomes lose the D-arm and/or the T-arm. The functions of these tRNAs need to be evaluated. The mitochondrial genome of T. malaysiensis differs from the other six genomes in having a slightly smaller genome size, a slight difference in codon usage, and a variable loop in place of the T-arm of some tRNAs by a variable loop. A phylogenic analysis shows that T. malaysiensis first split from other Tetranychus species and that the clade of the family Tetranychoidea occupies a basal position in the Trombidiformes. The mt genomes of the green and red forms of T. urticae have limited divergence and short evolutionary distance.