Complex gene expression in the dragline silk producing glands of the Western black widow (Latrodectus hesperus)

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 evolutionary history of cribellate orb-weaver capture thread spidroins

Background

Spiders have evolved two types of sticky capture threads: one with wet adhesive spun by ecribellate orb-weavers and another with dry adhesive spun by cribellate spiders. The evolutionary history of cribellate capture threads is especially poorly understood. Here, we use genomic approaches to catalog the spider-specific silk gene family (spidroins) for the cribellate orb-weaver Uloborus diversus.

Results

We show that the cribellar spidroin, which forms the puffy fibrils of cribellate threads, has three distinct repeat units, one of which is conserved across cribellate taxa separated by ~ 250 Mya. We also propose candidates for a new silk type, paracribellar spidroins, which connect the puffy fibrils to pseudoflagelliform support lines. Moreover, we describe the complete repeat architecture for the pseudoflagelliform spidroin (Pflag), which contributes to extensibility of pseudoflagelliform axial fibers.

Conclusions

Our finding that Pflag is closely related to Flag, supports homology of the support lines of cribellate and ecribellate capture threads. It further suggests an evolutionary phase following gene duplication, in which both Flag and Pflag were incorporated into the axial lines, with subsequent loss of Flag in uloborids, and increase in expression of Flag in ecribellate orb-weavers, explaining the distinct mechanical properties of the axial lines of these two groups.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12862-022-02042-5.

Egg Case Protein 3: A Constituent of Black Widow Spider Tubuliform Silk

Spider silk has outstanding mechanical properties, rivaling some of the best materials on the planet. Biochemical analyses of tubuliform silk have led to the identification of TuSp1, egg case protein 1, and egg case protein 2. TuSp1 belongs to the spidroin superfamily, containing a non-repetitive N- and C-terminal domain and internal block repeats. ECP1 and ECP2, which lack internal block repeats and sequence similarities to the highly conserved N- and C-terminal domains of spidroins, have cysteine-rich N-terminal domains. In this study, we performed an in-depth proteomic analysis of tubuliform glands, spinning dope, and egg sacs, which led to the identification of a novel molecular constituent of black widow tubuliform silk, referred to as egg case protein 3 or ECP3. Analysis of the translated ECP3 cDNA predicts a low molecular weight protein of 11.8 kDa. Real-time reverse transcription–quantitative PCR analysis performed with different silk-producing glands revealed ECP3 mRNA is predominantly expressed within tubuliform glands of spiders. Taken together, these findings reveal a novel protein that is secreted into black widow spider tubuliform silk.

Gene expression profiling reveals candidate genes for defining spider silk gland types

Molecular studies of the secretory glands involved in spider silk production have revealed candidate genes for silk synthesis and a complicated history of spider silk gene evolution. However, differential gene expression profiles of the multiple silk gland types within an individual orb-web weaving spider are lacking. Each of these gland types produces a functionally distinct silk type. Comparison of gene expression among spider silk gland types would provide insight into the genes that define silk glands generally from non-silk gland tissues, and the genes that define silk glands from each other. Here, we perform 3' tag digital gene expression profiling of the seven silk gland types of the silver garden orb weaver Argiope argentata. Five of these gland types produce silks that are non-adhesive fibers, one silk includes both fibers and glue-like adhesives, and one silk is exclusively glue-like. We identify 1275 highly expressed, significantly upregulated, and tissue specific silk gland specific transcripts (SSTs). These SSTs include seven types of spider silk protein encoding genes known as spidroin genes. We find that the fiber-producing major ampullate and minor ampullate silk glands have more similar expression profiles than any other pair of glands. We also find that a subset of the SSTs is enriched for transmembrane transport and oxidoreductases, and that these transcripts highlight differences and similarities among the major ampullate, minor ampullate, and aggregate silk glands. Furthermore, we show that the wet glue-producing aggregate glands have the most unique SSTs, but still share some SSTs with fiber producing glands. Aciniform glands were the only gland type to share a majority of SSTs with other silk gland types, supporting previous hypotheses that duplication of aciniform glands and subsequent divergence of the duplicates gave rise to the multiple silk gland types within an individual spider.

Shifts in morphology, gene expression, and selection underlie web loss in Hawaiian Tetragnatha spiders

BACKGROUND

A striking aspect of evolution is that it often converges on similar trajectories. Evolutionary convergence can occur in deep time or over short time scales, and is associated with the imposition of similar selective pressures. Repeated convergent events provide a framework to infer the genetic basis of adaptive traits. The current study examines the genetic basis of secondary web loss within web-building spiders (Araneoidea). Specifically, we use a lineage of spiders in the genus Tetragnatha (Tetragnathidae) that has diverged into two clades associated with the relatively recent (5 mya) colonization of, and subsequent adaptive radiation within, the Hawaiian Islands. One clade has adopted a cursorial lifestyle, and the other has retained the ancestral behavior of capturing prey with sticky orb webs. We explore how these behavioral phenotypes are reflected in the morphology of the spinning apparatus and internal silk glands, and the expression of silk genes. Several sister families to the Tetragnathidae have undergone similar web loss, so we also ask whether convergent patterns of selection can be detected in these lineages.

RESULTS

The cursorial clade has lost spigots associated with the sticky spiral of the orb web. This appears to have been accompanied by loss of silk glands themselves. We generated phylogenies of silk proteins (spidroins), which showed that the transcriptomes of cursorial Tetragnatha contain all major spidroins except for flagelliform. We also found an uncharacterized spidroin that has higher expression in cursorial species. We found evidence for convergent selection acting on this spidroin, as well as genes involved in protein metabolism, in the cursorial Tetragnatha and divergent cursorial lineages in the families Malkaridae and Mimetidae.

CONCLUSIONS

Our results provide strong evidence that independent web loss events and the associated adoption of a cursorial lifestyle are based on similar genetic mechanisms. Many genes we identified as having evolved convergently are associated with protein synthesis, degradation, and processing, which are processes that play important roles in silk production. This study demonstrates, in the case of independent evolution of web loss, that similar selective pressures act on many of the same genes to produce the same phenotypes and behaviors.

Structural Characterization of Black Widow Spider Dragline Silk Proteins CRP1 and CRP4

Spider dragline silk represents a biomaterial with outstanding mechanical properties, possessing high-tensile strength and toughness. In black widows at least eight different proteins have been identified as constituents of dragline silk. These represent major ampullate spidroins MaSp1, MaSp2, MaSp’, and several low-molecular weight cysteine-rich protein (CRP) family members, including CRP1, CRP2, and CRP4. Molecular modeling predicts that CRPs contain a cystine slipknot motif, but experimental evidence to support this assertion remains to be reported. To advance scientific knowledge regarding CRP function, we recombinantly expressed and purified CRP1 and CRP4 from bacteria and investigated their secondary structure using circular dichroism (CD) under different chemical and physical conditions. We demonstrate by far-UV CD spectroscopy that these proteins contain similar secondary structure, having substantial amounts of random coil conformation, followed by lower levels of beta sheet, alpha helical and beta turn structures. CRPs are thermally and pH stable; however, treatment with reagents that disrupt disulfide bonds impact their structural conformations. Cross-linking mass spectrometry (XL-MS) data also support computational models of CRP1. Taken together, the chemical and thermal stability of CRPs, the cross-linking data, coupled with the structural sensitivity to reducing agents, are experimentally consistent with the supposition CRPs are cystine slipknot proteins.

Towards engineering and production of artificial spider silk using tools of synthetic biology

Spider silk is one of the strongest biomaterials available in nature. Its mechanical properties make it a good candidate for applications in various fields ranging from protective armour to bandages for wound dressing to coatings for medical implants. Spider silk is formed by an intricate arrangement of spidroins, which are extremely large proteins containing long stretches of repeating segments rich in alanine and glycine. A large amount of research has been directed towards harnessing the spectacular potential of spider silks and using them for different applications. The interdisciplinary approach of synthetic biology is an ideal tool to study these spider silk proteins and work towards the engineering and production of synthetic spider silk. This review aims to highlight the recent progress that has been made in the study of spider silk proteins using different branches of synthetic biology. Here, the authors discuss the different computational approaches, directed evolution techniques and various expression platforms that have been tested for the successful production of spider silk. Future challenges facing the field and possible solutions offered by synthetic biology are also discussed.

Major ampullate silk gland transcriptomes and fibre proteomes of the golden orb-weavers, Nephila plumipes and Nephila pilipes (Araneae: Nephilidae)

Natural spider silk is one of the world’s toughest proteinaceous materials, yet a truly biomimetic spider silk is elusive even after several decades of intense focus. In this study, Next-Generation Sequencing was utilised to produce transcriptomes of the major ampullate gland of two Australian golden orb-weavers, Nephila plumipes and Nephila pilipes, in order to identify highly expressed predicted proteins that may co-factor in the construction of the final polymer. Furthermore, proteomics was performed by liquid chromatography tandem-mass spectroscopy to analyse the natural solid silk fibre of each species to confirm highly expressed predicted proteins within the silk gland are present in the final silk product. We assembled the silk gland transcriptomes of N. plumipes and N. pilipes into 69,812 and 70,123 contigs, respectively. Gene expression analysis revealed that silk gene sequences were among the most highly expressed and we were able to procure silk sequences from both species in excess of 1,300 amino acids. However, some of the genes with the highest expression values were not able to be identified from our proteomic analysis. Proteome analysis of “reeled” silk fibres of N. plumipes and N. pilipes revealed 29 and 18 proteins, respectively, most of which were identified as silk fibre proteins. This study is the first silk gland specific transcriptome and proteome analysis for these species and will assist in the future development of a biomimetic spider silk.

Multiscale mechanisms of nutritionally induced property variation in spider silks

Variability in spider major ampullate (MA) silk properties at different scales has proven difficult to determine and remains an obstacle to the development of synthetic fibers mimicking MA silk performance. A multitude of techniques may be used to measure multiscale aspects of silk properties. Here we fed five species of Araneoid spider solutions that either contained protein or were protein deprived and performed silk tensile tests, small and wide-angle X-ray scattering (SAXS/WAXS), amino acid composition analyses, and silk gene expression analyses, to resolve persistent questions about how nutrient deprivation induces variations in MA silk mechanical properties across scales. Our analyses found that the properties of each spider's silk varied differently in response to variations in their protein intake. We found changes in the crystalline and non-crystalline nanostructures to play specific roles in inducing the property variations we found. Across treatment MaSp expression patterns differed in each of the five species. We found that in most species MaSp expression and amino acid composition variations did not conform with our predictions based on a traditional MaSp expression model. In general, changes to the silk's alanine and proline compositions influenced the alignment of the proteins within the silk's amorphous region, which influenced silk extensibility and toughness. Variations in structural alignment in the crystalline and non-crystalline regions influenced ultimate strength independent of genetic expression. Our study provides the deepest insights thus far into the mechanisms of how MA silk properties vary from gene expression to nanostructure formations to fiber mechanics. Such knowledge is imperative for promoting the production of synthetic silk fibers.

Evolutionary shifts in gene expression decoupled from gene duplication across functionally distinct spider silk glands

Spider silk synthesis is an emerging model for the evolution of tissue-specific gene expression and the role of gene duplication in functional novelty, but its potential has not been fully realized. Accordingly, we quantified transcript (mRNA) abundance in seven silk gland types and three non-silk gland tissues for three cobweb-weaving spider species. Evolutionary analyses based on expression levels of thousands of homologous transcripts and phylogenetic reconstruction of 605 gene families demonstrated conservation of expression for each gland type among species. Despite serial homology of all silk glands, the expression profiles of the glue-forming aggregate glands were divergent from fiber-forming glands. Also surprising was our finding that shifts in gene expression among silk gland types were not necessarily coupled with gene duplication, even though silk-specific genes belong to multi-paralog gene families. Our results challenge widely accepted models of tissue specialization and significantly advance efforts to replicate silk-based high-performance biomaterials.

The Nephila clavipes genome highlights the diversity of spider silk genes and their complex expression

Spider silks are the toughest known biological materials, yet are lightweight and virtually invisible to the human immune system, and they thus have revolutionary potential for medicine and industry. Spider silks are largely composed of spidroins, a unique family of structural proteins. To investigate spidroin genes systematically, we constructed the first genome of an orb-weaving spider: the golden orb-weaver (Nephila clavipes), which builds large webs using an extensive repertoire of silks with diverse physical properties. We cataloged 28 Nephila spidroins, representing all known orb-weaver spidroin types, and identified 394 repeated coding motif variants and higher-order repetitive cassette structures unique to specific spidroins. Characterization of spidroin expression in distinct silk gland types indicates that glands can express multiple spidroin types. We find evidence of an alternatively spliced spidroin, a spidroin expressed only in venom glands, evolutionary mechanisms for spidroin diversification, and non-spidroin genes with expression patterns that suggest roles in silk production.

Silk gene expression of theridiid spiders: implications for male-specific silk use

Spiders (order Araneae) rely on their silks for essential tasks, such as dispersal, prey capture, and reproduction. Spider silks are largely composed of spidroins, members of a protein family that are synthesized in silk glands. As needed, silk stored in silk glands is extruded through spigots on the spinnerets. Nearly all studies of spider silks have been conducted on females; thus, little is known about male silk biology. To shed light on silk use by males, we compared silk gene expression profiles of mature males to those of females from three cob-web weaving species (Theridiidae). We de novo assembled species-specific male transcriptomes from Latrodectus hesperus, Latrodectus geometricus, and Steatoda grossa followed by differential gene expression analyses. Consistent with their complement of silk spigots, male theridiid spiders express appreciable amounts of aciniform, major ampullate, minor ampullate, and pyriform spidroin genes but not tubuliform spidroin genes. The relative expression levels of particular spidroin genes varied between sexes and species. Because mature males desert their prey-capture webs and become cursorial in their search for mates, we anticipated that major ampullate (dragline) spidroin genes would be the silk genes most highly expressed by males. Indeed, major ampullate spidroin genes had the highest expression in S. grossa males. However, minor ampullate spidroin genes were the most highly expressed spidroin genes in L. geometricus and L. hesperus males. Our expression profiling results suggest species-specific adaptive divergence of silk use by male theridiids.

Duplication and concerted evolution of MiSp-encoding genes underlie the material properties of minor ampullate silks of cobweb weaving spiders

BACKGROUND

Orb-web weaving spiders and their relatives use multiple types of task-specific silks. The majority of spider silk studies have focused on the ultra-tough dragline silk synthesized in major ampullate glands, but other silk types have impressive material properties. For instance, minor ampullate silks of orb-web weaving spiders are as tough as draglines, due to their higher extensibility despite lower strength. Differences in material properties between silk types result from differences in their component proteins, particularly members of the spidroin (spider fibroin) gene family. However, the extent to which variation in material properties within a single silk type can be explained by variation in spidroin sequences is unknown. Here, we compare the minor ampullate spidroins (MiSp) of orb-weavers and cobweb weavers. Orb-web weavers use minor ampullate silk to form the auxiliary spiral of the orb-web while cobweb weavers use it to wrap prey, suggesting that selection pressures on minor ampullate spidroins (MiSp) may differ between the two groups.

RESULTS

We report complete or nearly complete MiSp sequences from five cobweb weaving spider species and measure material properties of minor ampullate silks in a subset of these species. We also compare MiSp sequences and silk properties of our cobweb weavers to published data for orb-web weavers. We demonstrate that all our cobweb weavers possess multiple MiSp loci and that one locus is more highly expressed in at least two species. We also find that the proportion of β-spiral-forming amino acid motifs in MiSp positively correlates with minor ampullate silk extensibility across orb-web and cobweb weavers.

CONCLUSIONS

MiSp sequences vary dramatically within and among spider species, and have likely been subject to multiple rounds of gene duplication and concerted evolution, which have contributed to the diverse material properties of minor ampullate silks. Our sequences also provide templates for recombinant silk proteins with tailored properties.

Candidate egg case silk genes for the spider Argiope argentata from differential gene expression analyses

Orb-web weaving spiders produce a variety of task-specific silks from specialized silk glands. The genetics underlying the synthesis of specific silk types are largely unknown, and transcriptome analysis could be a powerful approach for identifying candidate genes. However, de novo assembly and expression profiling of silk glands with RNA-sequencing (RNAseq) are problematic because the few known gene transcripts for silk proteins are extremely long and highly repetitive. To identify candidate genes for tubuliform (egg case) silk synthesis by the orb-weaver Argiope argentata (Araneidae), we estimated transcript abundance using two sequencing methods: RNAseq reads from throughout the length of mRNA molecules, and 3' digital gene expression reads from the 3' region of mRNA molecules. Both analyses identified similar sets of genes as differentially expressed when comparing tubuliform and nonsilk gland tissue. However, incompletely assembled silk gene transcripts were identified as differentially expressed because of RNAseq read alignments to highly repetitive regions, confounding interpretation of RNAseq results. Homologues of egg case silk protein (ECP) genes were upregulated in tubuliform glands. This discovery is the first description of ECP homologues in an araneid. We also propose additional candidate genes involved in synthesis of tubuliform or other silk types.

Comprehensive Proteomic Analysis of Spider Dragline Silk from Black Widows: A Recipe to Build Synthetic Silk Fibers

The outstanding material properties of spider dragline silk fibers have been attributed to two spidroins, major ampullate spidroins 1 and 2 (MaSp1 and MaSp2). Although dragline silk fibers have been treated with different chemical solvents to elucidate the relationship between protein structure and fiber mechanics, there has not been a comprehensive proteomic analysis of the major ampullate (MA) gland, its spinning dope, and dragline silk using a wide range of chaotropic agents, inorganic salts, and fluorinated alcohols to elucidate their complete molecular constituents. In these studies, we perform in-solution tryptic digestions of solubilized MA glands, spinning dope and dragline silk fibers using five different solvents, followed by nano liquid chromatography coupled to tandem mass spectrometry (LC-MS/MS) analysis with an Orbitrap Fusion™ Tribrid™. To improve protein identification, we employed three different tryptic peptide fragmentation modes, which included collision-induced dissociation (CID), electron transfer dissociation (ETD), and high energy collision dissociation (HCD) to discover proteins involved in the silk assembly pathway and silk fiber. In addition to MaSp1 and MaSp2, we confirmed the presence of a third spidroin, aciniform spidroin 1 (AcSp1), widely recognized as the major constituent of wrapping silk, as a product of dragline silk. Our findings also reveal that MA glands, spinning dope, and dragline silk contain at least seven common proteins: three members of the Cysteine-Rich Protein Family (CRP1, CRP2 and CRP4), cysteine-rich secretory protein 3 (CRISP3), fasciclin and two uncharacterized proteins. In summary, this study provides a proteomic blueprint to construct synthetic silk fibers that most closely mimic natural fibers.

Proteomic Evidence for Components of Spider Silk Synthesis from Black Widow Silk Glands and Fibers

Spider silk research has largely focused on spidroins, proteins that are the primary components of spider silk fibers. Although a number of spidroins have been characterized, other types of proteins associated with silk synthesis are virtually unknown. Previous analyses of tissue-specific RNA-seq libraries identified 647 predicted genes that were differentially expressed in silk glands of the Western black widow, Latrodectus hesperus. Only ∼5% of these silk-gland specific transcripts (SSTs) encode spidroins; although the remaining predicted genes presumably encode other proteins associated with silk production, this is mostly unverified. Here, we used proteomic analysis of multiple silk glands and dragline silk fiber to investigate the translation of the differentially expressed genes. We find 48 proteins encoded by the differentially expressed transcripts in L. hesperus major ampullate, minor ampullate, and tubuliform silk glands and detect 17 SST encoded proteins in major ampullate silk fibers. The observed proteins include known silk-related proteins, but most are uncharacterized, with no annotation. These unannotated proteins likely include novel silk-associated proteins. Major and minor ampullate glands have the highest overlap of identified proteins, consistent with their shared, distinctive ampullate shape and the overlapping functions of major and minor ampullate silks. Our study substantiates and prioritizes predictions from differential expression analysis of spider silk gland transcriptomes.

Spider Transcriptomes Identify Ancient Large-Scale Gene Duplication Event Potentially Important in Silk Gland Evolution

The evolution of specialized tissues with novel functions, such as the silk synthesizing glands in spiders, is likely an influential driver of adaptive success. Large-scale gene duplication events and subsequent paralog divergence are thought to be required for generating evolutionary novelty. Such an event has been proposed for spiders, but not tested. We de novo assembled transcriptomes from three cobweb weaving spider species. Based on phylogenetic analyses of gene families with representatives from each of the three species, we found numerous duplication events indicative of a whole genome or segmental duplication. We estimated the age of the gene duplications relative to several speciation events within spiders and arachnids and found that the duplications likely occurred after the divergence of scorpions (order Scorpionida) and spiders (order Araneae), but before the divergence of the spider suborders Mygalomorphae and Araneomorphae, near the evolutionary origin of spider silk glands. Transcripts that are expressed exclusively or primarily within black widow silk glands are more likely to have a paralog descended from the ancient duplication event and have elevated amino acid replacement rates compared with other transcripts. Thus, an ancient large-scale gene duplication event within the spider lineage was likely an important source of molecular novelty during the evolution of silk gland-specific expression. This duplication event may have provided genetic material for subsequent silk gland diversification in the true spiders (Araneomorphae).

Dragline silk: a fiber assembled with low-molecular-weight cysteine-rich proteins

Dragline silk has been proposed to contain two main protein constituents, MaSp1 and MaSp2. However, the mechanical properties of synthetic spider silks spun from recombinant MaSp1 and MaSp2 proteins have yet to approach natural fibers, implying the natural spinning dope is missing critical factors. Here we report the discovery of novel molecular constituents within the spinning dope that are extruded into dragline silk. Protein studies of the liquid spinning dope from the major ampullate gland, coupled with the analysis of dragline silk fibers using mass spectrometry, demonstrate the presence of a new family of low-molecular-weight cysteine-rich proteins (CRPs) that colocalize with the MA fibroins. Expression of the CRP family members is linked to dragline silk production, specifically MaSp1 and MaSp2 mRNA synthesis. Biochemical data support that CRP molecules are secreted into the spinning dope and assembled into macromolecular complexes via disulfide bond linkages. Sequence analysis supports that CRP molecules share similarities to members that belong to the cystine slipknot superfamily, suggesting that these factors may have evolved to increase fiber toughness by serving as molecular hubs that dissipate large amounts of energy under stress. Collectively, our findings provide molecular details about the components of dragline silk, providing new insight that will advance materials development of synthetic spider silk for industrial applications.

Spider genomes provide insight into composition and evolution of venom and silk

Spiders are ecologically important predators with complex venom and extraordinarily tough silk that enables capture of large prey. Here we present the assembled genome of the social velvet spider and a draft assembly of the tarantula genome that represent two major taxonomic groups of spiders. The spider genomes are large with short exons and long introns, reminiscent of mammalian genomes. Phylogenetic analyses place spiders and ticks as sister groups supporting polyphyly of the Acari. Complex sets of venom and silk genes/proteins are identified. We find that venom genes evolved by sequential duplication, and that the toxic effect of venom is most likely activated by proteases present in the venom. The set of silk genes reveals a highly dynamic gene evolution, new types of silk genes and proteins, and a novel use of aciniform silk. These insights create new opportunities for pharmacological applications of venom and biomaterial applications of silk.

Spiders use self-produced venom and silk for their daily survival. Here, the authors report the assembled genome of the social velvet spider and a draft assembly of the tarantula genome and, together with proteomic data, provide insights into the evolution of genes that affect venom and silk production.