A masculinizing supergene underlies an exaggerated male reproductive morph in a spider

Karyotype and genome size analyses for two spiders of the lycosidae family

Background

Karyotype and genome size are critical genetic characteristics with significant value for cytogenetics, taxonomy, phylogenetics, evolution, and molecular biology. The Lycosidae family, known for its diverse spiders with varying ecological habits and behavioral traits, has seen limited exploration of its karyotype and genome size.

Methods

We utilized an improved tissue drop technique to prepare chromosome slides and compare the features of male and female karyotypes for two wolf spiders with different habits of Lycosidae. Furthermore, we predicted their genome sizes using flow cytometry (FCM) and K-mer analysis.

Results

The karyotypes of female and male Hippasa lycosina were 2n♀ = 26 = 14 m + 12 sm and 2n♂ = 24 = 10 m + 14 sm, respectively, and were composed of metacentric (m) and submetacentric (sm) chromosomes. In contrast, the karyotypes of Lycosa grahami consisted of telocentric (t) and subtelocentric (st) chromosomes (2n♀ = 20 = 20th and 2n♂ = 18 = 12th + 6t, for females and males). The sex chromosomes were both X1X2O. The estimated sizes of the H. lycosina and L. grahami genomes were 1966.54-2099.89 Mb and 3692.81-4012.56 Mb, respectively. Flow cytometry yielded slightly smaller estimates for genome size compared to k-mer analysis. K-mer analysis revealed a genome heterozygosity of 0.42% for H. lycosina and 0.80% for L. grahami, along with duplication ratios of 21.39% and 54.91%, respectively.

Conclusion

This study describes the first analysis of the genome sizes and karyotypes of two spiders from the Lycosidae that exhibit differential habits and provides essential data for future phylogenetic, cytogenetic, and genomic studies.

Genomic and transcriptomic analyses of Heteropoda venatoria reveal the expansion of P450 family for starvation resistance in spiders

BACKGROUND

Research on the mechanism of starvation resistance can help reveal how animals adjust their physiology and behavior to adapt to the uncertainty of food resources. A low metabolic rate is a significant characteristic of spider physiological activity and can increase spider starvation resistance and adapt to complex ecological environments.

RESULTS

We sequenced the genome of Heteropoda venatoria and discovered significant expansions in gene families related to lipid metabolism, such as cytochrome P450 and steroid hormone biosynthesis genes, through comparative genomic analysis. We also systematically analyzed the gene expression characteristics of H. venatoria at different starvation resistance stages and reported that the fat body plays a crucial role during starvation in spiders. This study indicates that during the early stages of starvation, H. venatoria relies on glucose metabolism to meet its energy demands. In the middle stage, gene expression stabilizes, whereas in the late stage of starvation, pathways for fatty acid metabolism and protein degradation are significantly activated, and autophagy is increased, serving as a survival strategy under extreme starvation. Notably, analysis of expanded P450 gene families revealed that H. venatoria has many duplicated CYP3 clan genes that are highly expressed in the fat body, which may help maintain a low-energy metabolic state, allowing H. venatoria to endure longer periods of starvation. We also observed that the motifs of P450 families in H. venatoria are less conserved than those in insects are, which may be related to the greater polymorphism of spider genomes.

CONCLUSIONS

This research not only provides important genetic and transcriptomic evidence for understanding the starvation mechanisms of spiders but also offers new insights into the adaptive evolution of arthropods.

Three Novel Spider Genomes Unveil Spidroin Diversification and Hox Cluster Architecture: Ryuthela nishihirai (Liphistiidae), Uloborus plumipes (Uloboridae) and Cheiracanthium punctorium (Cheiracanthiidae)

Spiders are a hyperdiverse taxon and among the most abundant predators in nearly all terrestrial habitats. Their success is often attributed to key developments in their evolution such as silk and venom production and major apomorphies such as a whole-genome duplication. Resolving deep relationships within the spider tree of life has been historically challenging, making it difficult to measure the relative importance of these novelties for spider evolution. Whole-genome data offer an essential resource in these efforts, but also for functional genomic studies. Here, we present de novo assemblies for three spider species: Ryuthela nishihirai (Liphistiidae), a representative of the ancient Mesothelae, the suborder that is sister to all other extant spiders; Uloborus plumipes (Uloboridae), a cribellate orbweaver whose phylogenetic placement is especially challenging; and Cheiracanthium punctorium (Cheiracanthiidae), which represents only the second family to be sequenced in the hyperdiverse Dionycha clade. These genomes fill critical gaps in the spider tree of life. Using these novel genomes along with 25 previously published ones, we examine the evolutionary history of spidroin gene and structural hox cluster diversity. Our assemblies provide critical genomic resources to facilitate deeper investigations into spider evolution. The near chromosome-level genome of the 'living fossil' R. nishihirai represents an especially important step forward, offering new insights into the origins of spider traits.

A supergene affects androgen concentrations during early development in a bird with alternative reproductive morphs

Species with alternative reproductive tactics typically show pronounced phenotypic variation between and within sexes. In some species, this variation culminates in discrete reproductive morphs that are genetically determined, facilitating studies on how genetic variation translates into phenotypic variation. In ruffs (Calidris pugnax), an autosomal inversion polymorphism underlies three reproductive morphs (Independents, Satellites and Faeders), which differ in circulating steroid concentrations in adults. Yet, it remains unknown whether morph differences in steroid concentrations already arise before adulthood. We examined variation in circulating testosterone, androstenedione and progesterone concentrations between morphs and sexes in ruff chicks and juveniles and compared the differences to those in adults. Since measured hormone concentrations only provide momentary states and show high within- and between-individual variation, we took repeated measurements to compare means, variances and skewness between groups. We found clear differences between morphs but not the sexes in early life. Between morphs, androgen concentrations in young ruffs differed in variance and skewness, but not in their means. For testosterone, Independents had a higher variance than Satellites/Faeders, whereas for androstenedione, we observed the opposite pattern. For progesterone, we did not detect clear differences between groups. Skewness values mirrored differences in morph variances. Compared to adults, premature ruffs had lower androgen concentrations. In both life stages, we detected morph-specific associations between androgen concentrations: androstenedione concentrations increased with testosterone concentrations more in Satellites/Faeders than in Independents. These observed morph differences during early life are consistent with a supergene-mediated regulation of androgen variation that underlies the diversification of adult behavioural phenotypes.

Supergenes are not necessary to explain the maintenance of complex alternative phenotypes

Evolutionary biology aims to explain the diversity seen in nature. Evolutionary theory provides frameworks to understand how simple polymorphisms or continuous variation are maintained, but phenotypes inherited as discrete suites of quantitative traits are difficult to fit into this framework. Supergenes have been proposed as a solution to this problem-if causal genes are co-located, they can be inherited as if a single gene, thus bridging the gap between simple polymorphisms and continuous traits. We develop models to ask: how are critical supergenes for maintaining phenotypic diversity? In our simplest model, without explicit genetic architectures, three alternative reproductive morphs are maintained in many of the parameter combinations we evaluated. For these same parameter values, models with demographic stochasticity, recombination and mutation (but without explicit genetic architecture) maintained only two of these three morphs, with stochasticity determining which morphs persisted. With explicit genetic architectures, regardless of whether causal loci were co-located in a supergene or distributed randomly, this stochasticity in which morphs are maintained was reduced. Even when phenotypic variation was lost, genetic diversity was maintained. Altogether, categorical traits with polygenic bases exhibited similar evolutionary dynamics to those determined by supergenes. Our work suggests that supergenes are not the only answer to the puzzle of how discrete polygenic phenotypic variation is maintained.

Identification and Evolutionary Analysis of the Widely Distributed CAP Superfamily in Spider Venom

Venom plays a crucial role in the defense and predation of venomous animals. Spiders (Araneae) are among the most successful predators and have a fascinating venom composition. Their venom mainly contains disulfide-rich peptides and large proteins. Here, we analyzed spider venom protein families, utilizing transcriptomic and genomic data, and highlighted their similarities and differences. We show that spiders have specific combinations of toxins for better predation and defense, typically comprising a core toxin expressed alongside several auxiliary toxins. Among them, the CAP superfamily is widely distributed and highly expressed in web-building Araneoidea spiders. Our analysis of evolutionary relationships revealed four subfamilies (subA-subD) of the CAP superfamily that differ in structure and potential functions. CAP proteins are composed of a conserved CAP domain and diverse C-terminal domains. CAP subC shares similar domains with the snake ion channel regulator svCRISP proteins, while CAP subD possesses a sequence similar to that of insect venom allergen 5 (Ag5). Furthermore, we show that gene duplication and selective expression lead to increased expression of CAP subD, making it a core member of the CAP superfamily. This study sheds light on the functional diversity of CAP subfamilies and their evolutionary history, which has important implications for fully understanding the composition of spider venom proteins and the core toxin components of web-building spiders.

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 genomics and evolution of inter-sexual mimicry and female-limited polymorphisms in damselflies

Sex-limited morphs can provide profound insights into the evolution and genomic architecture of complex phenotypes. Inter-sexual mimicry is one particular type of sex-limited polymorphism in which a novel morph resembles the opposite sex. While inter-sexual mimics are known in both sexes and a diverse range of animals, their evolutionary origin is poorly understood. Here, we investigated the genomic basis of female-limited morphs and male mimicry in the common bluetail damselfly. Differential gene expression between morphs has been documented in damselflies, but no causal locus has been previously identified. We found that male mimicry originated in an ancestrally sexually dimorphic lineage in association with multiple structural changes, probably driven by transposable element activity. These changes resulted in ~900 kb of novel genomic content that is partly shared by male mimics in a close relative, indicating that male mimicry is a trans-species polymorphism. More recently, a third morph originated following the translocation of part of the male-mimicry sequence into a genomic position ~3.5 mb apart. We provide evidence of balancing selection maintaining male mimicry, in line with previous field population studies. Our results underscore how structural variants affecting a handful of potentially regulatory genes and morph-specific genes can give rise to novel and complex phenotypic polymorphisms.

Lampshade web spider Ectatosticta davidi chromosome-level genome assembly provides evidence for its phylogenetic position

The spider of Ectatosticta davidi, belonging to the lamp-shade web spider family, Hypochilidae, which is closely related to Hypochilidae and Filistatidae and recovered as sister of the rest Araneomorphs spiders. Here we show the final assembled genome of E. davidi with 2.16 Gb in 15 chromosomes. Then we confirm the evolutionary position of Hypochilidae. Moreover, we find that the GMC gene family exhibit high conservation throughout the evolution of true spiders. We also find that the MaSp genes of E. davidi may represent an early stage of MaSp and MiSp genes in other true spiders, while CrSp shares a common origin with AgSp and PySp but differ from MaSp. Altogether, this study contributes to addressing the limited availability of genomic sequences from Hypochilidae spiders, and provides a valuable resource for investigating the genomic evolution of spiders.

A sexually selected male weapon characterized by strong additive genetic variance and no evidence for sexually antagonistic polyphenic maintenance

Sexual selection and sexual antagonism are important drivers of eco-evolutionary processes. The evolution of traits shaped by these processes depends on their genetic architecture, which remains poorly studied. Here, implementing a quantitative genetics approach using diallel crosses of the bulb mite, Rhizoglyphus robini, we investigated the genetic variance that underlies a sexually selected weapon that is dimorphic among males and female fecundity. Previous studies indicated that a negative genetic correlation between these two traits likely exists. We found male morph showed considerable additive genetic variance, which is unlikely to be explained solely by mutation-selection balance, indicating the likely presence of large-effect loci. However, a significant magnitude of inbreeding depression also indicates that morph expression is likely to be condition-dependent to some degree and that deleterious recessives can simultaneously contribute to morph expression. Female fecundity also showed a high degree of inbreeding depression, but the variance in female fecundity was mostly explained by epistatic effects, with very little contribution from additive effects. We found no significant genetic correlation, nor any evidence for dominance reversal, between male morph and female fecundity. The complex genetic architecture underlying male morph and female fecundity in this system has important implications for our understanding of the evolutionary interplay between purifying selection and sexually antagonistic selection.

Recent progress in understanding the genomic architecture of sexual conflict

Genomic conflict between the sexes over shared traits is widely assumed to be resolved through the evolution of sex-biased expression and the subsequent emergence of sexually dimorphic phenotypes. However, while there is support for a broad relationship between genome-wide patterns of expression level and sexual conflict, recent studies suggest that sex differences in the nature and strength of interactions between loci are instead key to conflict resolution. Furthermore, the advent of new technologies for measuring and perturbing expression means we now have much more power to detect genomic signatures of sexual conflict. Here, we review our current understanding of the genomic architecture of sexual conflict in the light of these new studies and highlight the potential for novel approaches to address outstanding knowledge gaps.

One to host them all: genomics of the diverse bacterial endosymbionts of the spider Oedothorax gibbosus

Bacterial endosymbionts of the groups Wolbachia, Cardinium and Rickettsiaceae are well known for their diverse effects on their arthropod hosts, ranging from mutualistic relationships to reproductive phenotypes. Here, we analysed a unique system in which the dwarf spider Oedothorax gibbosus is co-infected with up to five different endosymbionts affiliated with Wolbachia, 'Candidatus Tisiphia' (formerly Torix group Rickettsia), Cardinium and Rhabdochlamydia. Using short-read genome sequencing data, we show that the endosymbionts are heterogeneously distributed among O. gibbosus populations and are frequently found co-infecting spider individuals. To study this intricate host-endosymbiont system on a genome-resolved level, we used long-read sequencing to reconstruct closed genomes of the Wolbachia, 'Ca. Tisiphia' and Cardinium endosymbionts. We provide insights into the ecology and evolution of the endosymbionts and shed light on the interactions with their spider host. We detected high quantities of transposable elements in all endosymbiont genomes and provide evidence that ancestors of the Cardinium, 'Ca. Tisiphia' and Wolbachia endosymbionts have co-infected the same hosts in the past. Our findings contribute to broadening our knowledge about endosymbionts infecting one of the largest animal phyla on Earth and show the usefulness of transposable elements as an evolutionary 'contact-tracing' tool.

Sex-specific morphs: the genetics and evolution of intra-sexual variation

Sex-specific morphs exhibit discrete phenotypes, often including many disparate traits, that are observed in only one sex. These morphs have evolved independently in many different animals and are often associated with alternative mating strategies. The remarkable diversity of sex-specific morphs offers unique opportunities to understand the genetic basis of complex phenotypes, as the distinct nature of many morphs makes it easier to both categorize and compare genomes than for continuous traits. Sex-specific morphs also expand the study of sexual dimorphism beyond traditional bimodal comparisons of male and female averages, as they allow for a more expansive range of sexualization. Although ecological and endocrinological studies of sex-specific morphs have been advancing for some time, genomic and transcriptomic studies of morphs are far more recent. These studies reveal not only many different paths to the evolution of sex-specific morphs but also many commonalities, such as the role of sex-determining genes and hormone signalling in morph development, and the mixing of male and female traits within some morphs.

Chromosome-level genome and the identification of sex chromosomes in Uloborus diversus

The orb web is a remarkable example of animal architecture that is observed in families of spiders that diverged over 200 million years ago. While several genomes exist for araneid orb-weavers, none exist for other orb-weaving families, hampering efforts to investigate the genetic basis of this complex behavior. Here we present a chromosome-level genome assembly for the cribellate orb-weaving spider Uloborus diversus. The assembly reinforces evidence of an ancient arachnid genome duplication and identifies complete open reading frames for every class of spidroin gene, which encode the proteins that are the key structural components of spider silks. We identified the 2 X chromosomes for U. diversus and identify candidate sex-determining loci. This chromosome-level assembly will be a valuable resource for evolutionary research into the origins of orb-weaving, spidroin evolution, chromosomal rearrangement, and chromosomal sex determination in spiders.

A Novel Widespread MITE Element in the Repeat-Rich Genome of the Cardinium Endosymbiont of the Spider Oedothorax gibbosus

Free-living bacteria have evolved multiple times to become host-restricted endosymbionts. The transition from a free-living to a host-restricted lifestyle comes with a number of different genomic changes, including a massive loss of genes. In host-restricted endosymbionts, gene inactivation and genome reduction are facilitated by mobile genetic elements, mainly insertion sequences (ISs). ISs are small autonomous mobile elements, and one of, if not the most, abundant transposable elements in bacteria. Proliferation of ISs is common in some facultative endosymbionts, and is likely driven by the transmission bottlenecks, which increase the level of genetic drift. In this study, we present a manually curated genome annotation for a Cardinium endosymbiont of the dwarf spider Oedothorax gibbosus. Cardinium species are host-restricted endosymbionts that, similarly to ColbachiaWolbachia spp., include strains capable of manipulating host reproduction. Through the focus on mobile elements, the annotation revealed a rampant spread of ISs, extending earlier observations in other Cardinium genomes. We found that a large proportion of IS elements are pseudogenized, with many displaying evidence of recent inactivation. Most notably, we describe the lineage-specific emergence and spread of a novel IS-derived Miniature Inverted repeat Transposable Element (MITE), likely being actively maintained by intact copies of its parental IS982-family element. This study highlights the relevance of manual curation of these repeat-rich endosymbiont genomes for the discovery of novel MITEs, as well as the possible role these understudied elements might play in genome streamlining. IMPORTANCE Cardinium bacteria, a widespread symbiont lineage found across insects and nematodes, have been linked to reproductive manipulation of their hosts. However, the study of Cardinium has been hampered by the lack of comprehensive genomic resources. The high content of mobile genetic elements, namely, insertion sequences (ISs), has long complicated the analyses and proper annotations of these genomes. In this study, we present a manually curated annotation of the Cardinium symbiont of the spider Oedothorax gibbosus. Most notably, we describe a novel IS-like element found exclusively in this strain. We show that this mobile element likely evolved from a defective copy of its parental IS and then spread throughout the genome, contributing to the pseudogenization of several other mobile elements. We propose this element is likely being maintained by the intact copies of its parental IS element and that other similar elements in the genome could potentially follow this route.

A butterfly pan-genome reveals that a large amount of structural variation underlies the evolution of chromatin accessibility

Despite insertions and deletions being the most common structural variants (SVs) found across genomes, not much is known about how much these SVs vary within populations and between closely related species, nor their significance in evolution. To address these questions, we characterized the evolution of indel SVs using genome assemblies of three closely related Heliconius butterfly species. Over the relatively short evolutionary timescales investigated, up to 18.0% of the genome was composed of indels between two haplotypes of an individual Heliconius charithonia butterfly and up to 62.7% included lineage-specific SVs between the genomes of the most distant species (11 Mya). Lineage-specific sequences were mostly characterized as transposable elements (TEs) inserted at random throughout the genome and their overall distribution was similarly affected by linked selection as single nucleotide substitutions. Using chromatin accessibility profiles (i.e., ATAC-seq) of head tissue in caterpillars to identify sequences with potential cis-regulatory function, we found that out of the 31,066 identified differences in chromatin accessibility between species, 30.4% were within lineage-specific SVs and 9.4% were characterized as TE insertions. These TE insertions were localized closer to gene transcription start sites than expected at random and were enriched for sites with significant resemblance to several transcription factor binding sites with known function in neuron development in Drosophila We also identified 24 TE insertions with head-specific chromatin accessibility. Our results show high rates of structural genome evolution that were previously overlooked in comparative genomic studies and suggest a high potential for structural variation to serve as raw material for adaptive evolution.

Ecology and evolution of chlamydial symbionts of arthropods

The phylum Chlamydiae consists of obligate intracellular bacteria including major human pathogens and diverse environmental representatives. Here we investigated the Rhabdochlamydiaceae, which is predicted to be the largest and most diverse chlamydial family, with the few described members known to infect arthropod hosts. Using published 16 S rRNA gene sequence data we identified at least 388 genus-level lineages containing about 14 051 putative species within this family. We show that rhabdochlamydiae are mainly found in freshwater and soil environments, suggesting the existence of diverse, yet unknown hosts. Next, we used a comprehensive genome dataset including metagenome assembled genomes classified as members of the family Rhabdochlamydiaceae, and we added novel complete genome sequences of Rhabdochlamydia porcellionis infecting the woodlouse Porcellio scaber, and of 'Candidatus R. oedothoracis' associated with the linyphiid dwarf spider Oedothorax gibbosus. Comparative analysis of basic genome features and gene content with reference genomes of well-studied chlamydial families with known host ranges, namely Parachlamydiaceae (protist hosts) and Chlamydiaceae (human and other vertebrate hosts) suggested distinct niches for members of the Rhabdochlamydiaceae. We propose that members of the family represent intermediate stages of adaptation of chlamydiae from protists to vertebrate hosts. Within the genus Rhabdochlamydia, pronounced genome size reduction could be observed (1.49-1.93 Mb). The abundance and genomic distribution of transposases suggests transposable element expansion and subsequent gene inactivation as a mechanism of genome streamlining during adaptation to new hosts. This type of genome reduction has never been described before for any member of the phylum Chlamydiae. This study provides new insights into the molecular ecology, genomic diversity, and evolution of representatives of one of the most divergent chlamydial families.