Allelic resolution of insect and spider silk genes reveals hidden genetic diversity.

Arthropod silk is vital to the evolutionary success of hundreds of thousands of species. The primary proteins in silks are often encoded by long, repetitive gene sequences. Until recently, sequencing and assembling these complex gene sequences has proven intractable given their repetitive structure. Here, using high-quality long-read sequencing, we show that there is extensive variation-both in terms of length and repeat motif order-between alleles of silk genes within individual arthropods. Further, this variation exists across two deep, independent origins of silk which diverged more than 500 Mya: the insect clade containing caddisflies and butterflies and spiders. This remarkable convergence in previously overlooked patterns of allelic variation across multiple origins of silk suggests common mechanisms for the generation and maintenance of structural protein-coding genes. Future genomic efforts to connect genotypes to phenotypes should account for such allelic variation.

Rapid molecular diversification and homogenization of clustered major ampullate silk genes in Argiope garden spiders.

The evolutionary diversification of orb-web weaving spiders is closely tied to the mechanical performance of dragline silk. This proteinaceous fiber provides the primary structural framework of orb web architecture, and its extraordinary toughness allows these structures to absorb the high energy of aerial prey impact. The dominant model of dragline silk molecular structure involves the combined function of two highly repetitive, spider-specific, silk genes (spidroins)-MaSp1 and MaSp2. Recent genomic studies, however, have suggested this framework is overly simplistic, and our understanding of how MaSp genes evolve is limited. Here we present a comprehensive analysis of MaSp structural and evolutionary diversity across species of Argiope (garden spiders). This genomic analysis reveals the largest catalog of MaSp genes found in any spider, driven largely by an expansion of MaSp2 genes. The rapid diversification of Argiope MaSp genes, located primarily in a single genomic cluster, is associated with profound changes in silk gene structure. MaSp2 genes, in particular, have evolved complex hierarchically organized repeat units (ensemble repeats) delineated by novel introns that exhibit remarkable evolutionary dynamics. These repetitive introns have arisen independently within the genus, are highly homogenized within a gene, but diverge rapidly between genes. In some cases, these iterated introns are organized in an alternating structure in which every other intron is nearly identical in sequence. We hypothesize that this intron structure has evolved to facilitate homogenization of the coding sequence. We also find evidence of intergenic gene conversion and identify a more diverse array of stereotypical amino acid repeats than previously recognized. Overall, the extreme diversification found among MaSp genes requires changes in the structure-function model of dragline silk performance that focuses on the differential use and interaction among various MaSp paralogs as well as the impact of ensemble repeat structure and different amino acid motifs on mechanical behavior.

Expression of spidroin proteins in the silk glands of golden orb-weaver spiders.

The expression of spidroins in the major ampullate, minor ampullate, flagelliform, and tubuliform silk glands of Trichonephila clavipes spiders was analyzed using proteomics analysis techniques. Spidroin peptides were identified and assigned to different gene products based on sequence concurrence when compared with the whole genome of the spider. It was found that only a relatively low proportion of the spidroin genes are expressed as proteins in any of the studied glands. In addition, the expression of spidroin genes in different glands presents a wide range of patterns, with some spidroins being found in a single gland exclusively, while others appear in the content of several glands. The combination of precise genomics, proteomics, microstructural, and mechanical data provides new insights both on the design principles of these materials and how these principles might be translated for the production of high-performance bioinspired artificial fibers.

Connecting materials, performance and evolution: a case study of the glue of moth-catching spiders (Cyrtarachninae).

Morphological structures and extended phenotypes are made possible by materials that are encoded by the genome. Nearly all biomaterials are viscoelastic, which means that to understand performance, one must understand the strain rate-dependent properties of these materials in relevant ecological interactions, as the behavior of a material can vary dramatically and rapidly. Spider silks are an example of materials whose properties vary substantially intra- and inter-specifically. Here, we focus on aggregate silk, which functions as a biological adhesive. As a case study to understand how a material manifests from genome through organism to ecology, we highlight moth-specialist spiders, the Cyrtarachninae, and their glues as an ideal experimental system to investigate the relationship between genomics and ecologically variable performance of a biological material. There is a clear eco-evolutionary innovation that Cyrtarachne akirai and related species have evolved, a unique trait not found in other spiders, a glue which overcomes the scales of moths. By examining traditional orb-weavers, C. akirai and other subfamily members using biomechanical testing and genomic analysis, we argue that we can track the evolution of this novel bioadhesive and comment on the selection pressures influencing prey specialization. The importance of the ecological context of materials testing is exemplified by the poor performance of C. akirai glue on glass and the exceptional spreading ability and adhesive strength on moths. The genetic basis for these performance properties is experimentally tractable because spider silk genes are minimally pleiotropic and advances in genomic technologies now make possible the discovery of complete silk gene sequences.

Ultrastructures and Mechanics of Annealed Nephila clavipes Major Ampullate Silk.

The semicrystalline protein structure and impressive mechanical properties of major ampullate (MA) spider silk make it a promising natural alternative to polyacrylonitrile (PAN) fibers for carbon fiber manufacture. However, when annealed using a similar procedure to carbon fiber production, the tensile strength and Young's modulus of MA silk decrease. Despite this, MA silk fibers annealed at 600 °C remain stronger and tougher than similarly annealed PAN but have a lower Young's modulus. Although MA silk and PAN graphitize to similar extents, annealing disrupts the hydrogen bonding that controls crystal alignment within MA silk. Consequently, unaligned graphite crystals form in annealed MA silk, causing it to weaken, while graphite crystals in PAN maintain alignment along the fiber axis, strengthening the fibers. These shortcomings of spider silk when annealed provide insights into the selection and design of future alternative carbon fiber precursors.

Phylogeny of the orb-weaving spider family Araneidae (Araneae: Araneoidea).

We present a new phylogeny of the spider family Araneidae based on five genes (28S, 18S, COI, H3 and 16S) for 158 taxa, identified and mainly sequenced by us. This includes 25 outgroups and 133 araneid ingroups representing the subfamilies Zygiellinae Simon, 1929, Nephilinae Simon, 1894, and the typical araneids, here informally named the "ARA Clade". The araneid genera analysed here include roughly 90% of all currently named araneid species. The ARA Clade is the primary focus of this analysis. In taxonomic terms, outgroups comprise 22 genera and 11 families, and the ingroup comprises three Zygiellinae and four Nephilinae genera, and 85 ARA Clade genera (ten new). Within the ARA Clade, we recognize ten informal groups that contain at least three genera each and are supported under Bayesian posterior probabilities (≥ 0.95): "Caerostrines" (Caerostris, Gnolus and Testudinaria), "Micrathenines" (Acacesia, Micrathena, Ocrepeira, Scoloderus and Verrucosa), "Eriophorines" (Acanthepeira, Alpaida, Eriophora, Parawixia and Wagneriana), "Backobourkiines" (Acroaspis, Backobourkia, Carepalxis, Novakiella, Parawixia, Plebs, Singa and three new genera), "Argiopines" (Arachnura, Acusilas, Argiope, Cyrtophora, Gea, Lariniaria and Mecynogea), "Cyrtarachnines" (Aranoethra, Cyrtarachne, Paraplectana, Pasilobus and Poecilopachys), "Mastophorines" (Celaenia, Exechocentrus and Mastophora,), "Nuctenines" (Larinia, Larinioides and Nuctenea), "Zealaraneines" (Colaranea, Cryptaranea, Paralarinia, Zealaranea and two new genera) and "Gasteracanthines" (Augusta, Acrosomoides, Austracantha, Gasteracantha, Isoxya, Macracantha, Madacantha, Parmatergus and Thelacantha). Few of these groups are currently corroborated by morphology, behaviour, natural history or biogeography. We also include the large genus Araneus, along with Aculepeira, Agalenatea, Anepsion, Araniella, Cercidia, Chorizopes, Cyclosa, Dolophones, Eriovixia, Eustala, Gibbaranea, Hingstepeira, Hypognatha, Kaira, Larinia, Mangora, Metazygia, Metepeira, Neoscona, Paraplectanoides, Perilla, Poltys, Pycnacantha, Spilasma and Telaprocera, but the placement of these genera was generally ambiguous, except for Paraplectanoides, which is strongly supported as sister to traditional Nephilinae. Araneus, Argiope, Eriophora and Larinia are polyphyletic, Araneus implying nine new taxa of genus rank, and Eriophora and Larinia two each. In Araneus and Eriophora, polyphyly was usually due to north temperate generic concepts being used as dumping grounds for species from southern hemisphere regions, e.g. South-East Asia, Australia or New Zealand. Although Araneidae is one of the better studied spider families, too little natural history and/or morphological data are available across these terminals to draw any strong evolutionary conclusions. However, the classical orb web is reconstructed as plesiomorphic for Araneidae, with a single loss in "cyrtarachnines"-"mastophorines". Web decorations (collectively known as stabilimenta) evolved perhaps five times. Sexual dimorphism generally results from female body size increase with few exceptions; dimorphic taxa are not monophyletic and revert to monomorphism in a few cases.

The developmental and genetic trajectory of coloration in the guppy (Poecilia reticulata).

Examining the association between trait variation and development is crucial for understanding the evolution of phenotypic differences. Male guppy ornamental caudal fin coloration is one trait that shows a striking degree of variation within and between guppy populations. Males initially have no caudal fin coloration, then gradually develop it as they reach sexual maturity. For males, there is a trade-off between female preference for caudal fin coloration and increased visibility to predators. This trade-off may reach unique endpoints in males from different predation regimes. Caudal fin coloration includes black melanin, orange/yellow pteridines or carotenoids, and shimmering iridescence. This study examined the phenotypic trajectory and genetics associated with color development. We found that black coloration always developed first, followed by orange/yellow, then iridescence. The ordering and timing of color appearance was the same regardless of predation regime. The increased expression of melanin synthesis genes correlated well with the visual appearance of black coloration, but there was no correlation between carotenoids or pteridine synthesis gene expression and the appearance of orange/yellow. The lack of orange/yellow coloration in earlier male caudal fin developmental stages may be due to reduced expression of genes underlying the development of orange/yellow xanthophores.

Cryptic elevational zonation in trapdoor spiders (Araneae, Antrodiaetidae, Aliatypus janus complex) from the California southern Sierra Nevada.

The relative roles of ecological niche conservatism versus niche divergence in promoting montane speciation remains an important topic in biogeography. Here, our aim was to test whether lineage diversification in a species complex of trapdoor spiders corresponds with riverine barriers or with an ecological gradient associated with elevational tiering. Aliatypus janus was sampled from throughout its range, with emphasis on populations in the southern Sierra Nevada Mountains of California. We collected multi-locus genetic data to generate a species tree for A. janus and its close relatives. Coalescent based hypothesis tests were conducted to determine if genetic breaks within A. janus conform to riverine barriers. Ecological niche models (ENM) under current and Last Glacial Maximum (LGM) conditions were generated and hypothesis tests of niche conservatism and divergence were performed. Coalescent analyses reveal deeply divergent genetic lineages within A. janus, likely corresponding to cryptic species. Two primary lineages meet along an elevational gradient on the western slopes of the southern Sierra Nevada Mountains. ENMs under both current and LGM conditions indicate that these groups occupy largely non-overlapping niches. ENM hypothesis testing rejected niche identity between the two groups, and supported a sharp ecological gradient occurring where the groups meet. However, the niche similarity test indicated that the two groups may not inhabit different background niches. The Sierra Nevada Mountains provide a natural laboratory for simultaneously testing ecological niche divergence and conservatism and their role in speciation across a diverse range of taxa. Aliatypus janus represents a species complex with cryptic lineages that may have diverged due to parapatric speciation along an ecological gradient, or been maintained by the evolution of ecological niche differences following allopatric speciation.

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.

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.

Phylogenomic analyses resolve an ancient trichotomy at the base of Ischyropsalidoidea (Arachnida, Opiliones) despite high levels of gene tree conflict and unequal minority resolution frequencies.

Phylogenetic resolution of ancient rapid radiations has remained problematic despite major advances in statistical approaches and DNA sequencing technologies. Here we report on a combined phylogenetic approach utilizing transcriptome data in conjunction with Sanger sequence data to investigate a tandem of ancient divergences in the harvestmen superfamily Ischyropsalidoidea (Arachnida, Opiliones, Dyspnoi). We rely on Sanger sequences to resolve nodes within and between closely related genera, and use RNA-seq data from a subset of taxa to resolve a short and ancient internal branch. We use several analytical approaches to explore this succession of ancient diversification events, including concatenated and coalescent-based analyses and maximum likelihood gene trees for each locus. We evaluate the robustness of phylogenetic inferences using a randomized locus sub-sampling approach, and find congruence across these methods despite considerable incongruence across gene trees. Incongruent gene trees are not recovered in frequencies expected from a simple multispecies coalescent model, and we reject incomplete lineage sorting as the sole contributor to gene tree conflict. Using these approaches we attain robust support for higher-level phylogenetic relationships within Ischyropsalidoidea.

Inactivation of C4orf26 in toothless placental mammals.

Previous studies have reported inactivated copies of six enamel-related genes (AMBN, AMEL, AMTN, ENAM, KLK4, MMP20) and one dentin-related gene (DSPP) in one or more toothless vertebrates and/or vertebrates with enamelless teeth, thereby providing evidence that these genes are enamel or tooth-specific with respect to their critical functions that are maintained by natural selection. Here, we employ available genome sequences for edentulous and enamelless mammals to evaluate the enamel specificity of four genes (WDR72, SLC24A4, FAM83H, C4orf26) that have been implicated in amelogenesis imperfecta, a condition in which proper enamel formation is abrogated during tooth development. Coding sequences for WDR72, SCL24A4, and FAM83H are intact in four edentulous taxa (Chinese pangolin, three baleen whales) and three taxa (aardvark, nine-banded armadillo, Hoffmann's two-toed sloth) with enamelless teeth, suggesting that these genes have critical functions beyond their involvement in tooth development. By contrast, genomic data for C4orf26 reveal inactivating mutations in pangolin and bowhead whale as well as evidence for deletion of this gene in two minke whale species. Hybridization capture of exonic regions and PCR screens provide evidence for inactivation of C4orf26 in eight additional baleen whale species. However, C4orf26 is intact in all three species with enamelless teeth that were surveyed, as well as in 95 additional mammalian species with enamel-capped teeth. Estimates of selection intensity suggest that dN/dS ratios on branches leading to taxa with enamelless teeth are similar to the dN/dS ratio on branches leading to taxa with enamel-capped teeth. Based on these results, we conclude that C4orf26 is tooth-specific, but not enamel-specific, with respect to its essential functions that are maintained by natural selection. A caveat is that an alternative splice site variant, which translates exon 3 in a different reading frame, is putatively functional in Catarrhini and may have evolved an additional role in this primate clade.

Efficacy of Several Pesticide Products on Brown Widow Spider (Araneae: Theridiidae) Egg Sacs and Their Penetration Through the Egg Sac Silk.

Information on pesticide effects on spiders is less common than for insects; similar information for spider egg sacs is scarcer in the open literature. Spider egg sacs are typically covered with a protective silk layer. When pesticides are directly applied to egg sacs, the silk might prevent active ingredients from reaching the eggs, blocking their insecticidal effect. We investigated the impact of six water-based pesticide sprays and four oil-based aerosol products against egg sacs of brown widow spiders, Latrodectus geometricus C. L. Koch. All water-based spray products except one failed to provide significant mortality to egg sacs, resulting in successful spiderling emergence from treated egg sacs at a similar rate to untreated egg sacs. In contrast to water-based sprays, oil-based aerosols provided almost complete control, with 94-100% prevention of spiderling emergence. Penetration studies using colored pesticide products indicated that oil-based aerosols were significantly more effective in penetrating egg sac silk than were the water-based sprays, delivering the active ingredients on most (>99%) of the eggs inside the sac. The ability of pesticides to penetrate spider egg sac silk and deliver lethal doses of active ingredients to the eggs is discussed in relation to the chemical nature of egg sac silk proteins. Our study suggests that pest management procedures primarily relying on perimeter application of water-based sprays might not provide satisfactory control of brown widow spider eggs. Determination of the most effective active ingredients and carrier characteristics warrant further research to provide more effective control options for spider egg sacs.

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.

Intragenic homogenization and multiple copies of prey-wrapping silk genes in Argiope garden spiders.

BACKGROUND: Spider silks are spectacular examples of phenotypic diversity arising from adaptive molecular evolution. An individual spider can produce an array of specialized silks, with the majority of constituent silk proteins encoded by members of the spidroin gene family. Spidroins are dominated by tandem repeats flanked by short, non-repetitive N- and C-terminal coding regions. The remarkable mechanical properties of spider silks have been largely attributed to the repeat sequences. However, the molecular evolutionary processes acting on spidroin terminal and repetitive regions remain unclear due to a paucity of complete gene sequences and sampling of genetic variation among individuals. To better understand spider silk evolution, we characterize a complete aciniform spidroin gene from an Argiope orb-weaving spider and survey aciniform gene fragments from congeneric individuals. RESULTS: We present the complete aciniform spidroin (AcSp1) gene from the silver garden spider Argiope argentata (Aar_AcSp1), and document multiple AcSp1 loci in individual genomes of A. argentata and the congeneric A. trifasciata and A. aurantia. We find that Aar_AcSp1 repeats have >98% pairwise nucleotide identity. By comparing AcSp1 repeat amino acid sequences between Argiope species and with other genera, we identify regions of conservation over vast amounts of evolutionary time. Through a PCR survey of individual A. argentata, A. trifasciata, and A. aurantia genomes, we ascertain that AcSp1 repeats show limited variation between species whereas terminal regions are more divergent. We also find that average dN/dS across codons in the N-terminal, repetitive, and C-terminal encoding regions indicate purifying selection that is strongest in the N-terminal region. CONCLUSIONS: Using the complete A. argentata AcSp1 gene and spidroin genetic variation between individuals, this study clarifies some of the molecular evolutionary processes underlying the spectacular mechanical attributes of aciniform silk. It is likely that intragenic concerted evolution and functional constraints on A. argentata AcSp1 repeats result in extreme repeat homogeneity. The maintenance of multiple AcSp1 encoding loci in Argiope genomes supports the hypothesis that Argiope spiders require rapid and efficient protein production to support their prolific use of aciniform silk for prey-wrapping and web-decorating. In addition, multiple gene copies may represent the early stages of spidroin diversification.

Dramatic expansion of the black widow toxin arsenal uncovered by multi-tissue transcriptomics and venom proteomics.

BACKGROUND: Animal venoms attract enormous interest given their potential for pharmacological discovery and understanding the evolution of natural chemistries. Next-generation transcriptomics and proteomics provide unparalleled, but underexploited, capabilities for venom characterization. We combined multi-tissue RNA-Seq with mass spectrometry and bioinformatic analyses to determine venom gland specific transcripts and venom proteins from the Western black widow spider (Latrodectus hesperus) and investigated their evolution. RESULTS: We estimated expression of 97,217 L. hesperus transcripts in venom glands relative to silk and cephalothorax tissues. We identified 695 venom gland specific transcripts (VSTs), many of which BLAST and GO term analyses indicate may function as toxins or their delivery agents. ~38% of VSTs had BLAST hits, including latrotoxins, inhibitor cystine knot toxins, CRISPs, hyaluronidases, chitinase, and proteases, and 59% of VSTs had predicted protein domains. Latrotoxins are venom toxins that cause massive neurotransmitter release from vertebrate or invertebrate neurons. We discovered ≥ 20 divergent latrotoxin paralogs expressed in L. hesperus venom glands, significantly increasing this biomedically important family. Mass spectrometry of L. hesperus venom identified 49 proteins from VSTs, 24 of which BLAST to toxins. Phylogenetic analyses showed venom gland specific gene family expansions and shifts in tissue expression. CONCLUSIONS: Quantitative expression analyses comparing multiple tissues are necessary to identify venom gland specific transcripts. We present a black widow venom specific exome that uncovers a trove of diverse toxins and associated proteins, suggesting a dynamic evolutionary history. This justifies a reevaluation of the functional activities of black widow venom in light of its emerging complexity.

Multi-tissue transcriptomics of the black widow spider reveals expansions, co-options, and functional processes of the silk gland gene toolkit.

BACKGROUND: Spiders (Order Araneae) are essential predators in every terrestrial ecosystem largely because they have evolved potent arsenals of silk and venom. Spider silks are high performance materials made almost entirely of proteins, and thus represent an ideal system for investigating genome level evolution of novel protein functions. However, genomic level resources remain limited for spiders. RESULTS: We de novo assembled a transcriptome for the Western black widow (Latrodectus hesperus) from deeply sequenced cDNAs of three tissue types. Our multi-tissue assembly contained ~100,000 unique transcripts, of which > 27,000 were annotated by homology. Comparing transcript abundance among the different tissues, we identified 647 silk gland-specific transcripts, including the few known silk fiber components (e.g. six spider fibroins, spidroins). Silk gland specific transcripts are enriched compared to the entire transcriptome in several functions, including protein degradation, inhibition of protein degradation, and oxidation-reduction. Phylogenetic analyses of 37 gene families containing silk gland specific transcripts demonstrated novel gene expansions within silk glands, and multiple co-options of silk specific expression from paralogs expressed in other tissues. CONCLUSIONS: We propose a transcriptional program for the silk glands that involves regulating gland specific synthesis of silk fiber and glue components followed by protecting and processing these components into functional fibers and glues. Our black widow silk gland gene repertoire provides extensive expansion of resources for biomimetic applications of silk in industry and medicine. Furthermore, our multi-tissue transcriptome facilitates evolutionary analysis of arachnid genomes and adaptive protein systems.

Molecular evolution of α-latrotoxin, the exceptionally potent vertebrate neurotoxin in black widow spider venom.

Black widow spiders (members of the genus Latrodectus) are widely feared because of their potent neurotoxic venom. α-Latrotoxin is the vertebrate-specific toxin responsible for the dramatic effects of black widow envenomation. The evolution of this toxin is enigmatic because only two α-latrotoxin sequences are known. In this study, ~4 kb α-latrotoxin sequences and their homologs were characterized from a diversity of Latrodectus species, and representatives of Steatoda and Parasteatoda, establishing the wide distribution of latrotoxins across the mega-diverse spider family Theridiidae. Across black widow species, α-latrotoxin shows ≥ 94% nucleotide identity and variability consistent with purifying selection. Multiple codon and branch-specific estimates of the nonsynonymous/synonymous substitution rate ratio also suggest a long history of purifying selection has acted on α-latrotoxin across Latrodectus and Steatoda. However, α-latrotoxin is highly divergent in amino acid sequence between these genera, with 68.7% of protein differences involving non-conservative substitutions, evidence for positive selection on its physiochemical properties and particular codons, and an elevated rate of nonsynonymous substitutions along α-latrotoxin's Latrodectus branch. Such variation likely explains the efficacy of red-back spider, L. hasselti, antivenom in treating bites from other Latrodectus species, and the weaker neurotoxic symptoms associated with Steatoda and Parasteatoda bites. Long-term purifying selection on α-latrotoxin indicates its functional importance in black widow venom, even though vertebrates are a small fraction of their diet. The greater differences between Latrodectus and Steatoda α-latrotoxin, and their relationships to invertebrate-specific latrotoxins, suggest a shift in α-latrotoxin toward increased vertebrate toxicity coincident with the evolution of widow spiders.

Ancient properties of spider silks revealed by the complete gene sequence of the prey-wrapping silk protein (AcSp1).

Spider silk fibers have impressive mechanical properties and are primarily composed of highly repetitive structural proteins (termed spidroins) encoded by a single gene family. Most characterized spidroin genes are incompletely known because of their extreme size (typically >9 kb) and repetitiveness, limiting understanding of the evolutionary processes that gave rise to their unusual gene architectures. The only complete spidroin genes characterized thus far form the dragline in the Western black widow, Latrodectus hesperus. Here, we describe the first complete gene sequence encoding the aciniform spidroin AcSp1, the primary component of spider prey-wrapping fibers. L. hesperus AcSp1 contains a single enormous (∼19 kb) exon. The AcSp1 repeat sequence is exceptionally conserved between two widow species (∼94% identity) and between widows and distantly related orb-weavers (∼30% identity), consistent with a history of strong purifying selection on its amino acid sequence. Furthermore, the 16 repeats (each 371-375 amino acids long) found in black widow AcSp1 are, on average, >99% identical at the nucleotide level. A combination of stabilizing selection on amino acid sequence, selection on silent sites, and intragenic recombination likely explains the extreme homogenization of AcSp1 repeats. In addition, phylogenetic analyses of spidroin paralogs support a gene duplication event occurring concomitantly with specialization of the aciniform glands and the tubuliform glands, which synthesize egg-case silk. With repeats that are dramatically different in length and amino acid composition from dragline spidroins, our L. hesperus AcSp1 expands the knowledge base for developing silk-based biomimetic technologies.

Special Prey, Special Glue: NMR Spectroscopy on Aggregate Glue Components of Moth-Specialist Spiders, Cyrtarachninae.

Orb-weaver spiders produce upwards of seven different types of silk, each with unique material properties. We focus on the adhesive within orb-weaving spider webs, aggregate glue silk. These droplets are composed of three main components: water, glycoproteins, and a wide range of low molecular mass compounds (LMMCs). These LMMCs are known to play a crucial role in maintaining the material properties of the glycoproteins, aid in water absorption from the environment, and increase surface adhesion. Orb-weavers within the Cyrtarachninae subfamily are moth specialists and have evolved glue droplets with novel material properties. This study investigated the biochemical composition and diversity of the LMMCs present in the aggregate glue of eight moth-specialist species and compared them with five generalist orb-weavers using nuclear magnetic resonance (NMR) spectroscopy. We hypothesized that the novel drying ability of moth-specialist glue was accompanied by novel LMMCs and lower overall percentages by silk weight of LMMCs. We measured no difference in LMMC weight by the type of prey specialization, but observed novel compositions in the glue of all eight moth-catching species. Further, we quantified the presence of a previously reported but unidentified compound that appears in the glue of all moth specialists. These silks can provide insight into the functions of bioadhesives and inform our own synthetic adhesives.