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.

Spider phylosymbiosis: divergence of widow spider species and their tissues' microbiomes.

BACKGROUND: Microbiomes can have profound impacts on host biology and evolution, but to date, remain vastly understudied in spiders despite their unique and diverse predatory adaptations. This study evaluates closely related species of spiders and their host-microbe relationships in the context of phylosymbiosis, an eco-evolutionary pattern where the microbial community profile parallels the phylogeny of closely related host species. Using 16S rRNA gene amplicon sequencing, we characterized the microbiomes of five species with known phylogenetic relationships from the family Theridiidae, including multiple closely related widow spiders (L. hesperus, L. mactans, L. geometricus, S. grossa, and P. tepidariorum). RESULTS: We compared whole animal and tissue-specific microbiomes (cephalothorax, fat bodies, venom glands, silk glands, and ovary) in the five species to better understand the relationship between spiders and their microbial symbionts. This showed a strong congruence of the microbiome beta-diversity of the whole spiders, cephalothorax, venom glands, and silk glands when compared to their host phylogeny. Our results support phylosymbiosis in these species and across their specialized tissues. The ovary tissue microbial dendrograms also parallel the widow phylogeny, suggesting vertical transfer of species-specific bacterial symbionts. By cross-validating with RNA sequencing data obtained from the venom glands, silk glands and ovaries of L. hesperus, L. geometricus, S. grossa, and P. tepidariorum we confirmed that several microbial symbionts of interest are viably active in the host. CONCLUSION: Together these results provide evidence that supports the importance of host-microbe interactions and the significant role microbial communities may play in the evolution and adaptation of their hosts.

Structural and functional characterization of sticholysin III: A newly discovered actinoporin within the venom of the sea anemone Stichodactyla helianthus.

Actinoporins are a family of pore-forming toxins produced by sea anemones as part of their venomous cocktail. These proteins remain soluble and stably folded in aqueous solution, but when interacting with sphingomyelin-containing lipid membranes, they become integral oligomeric membrane structures that form a pore permeable to cations, which leads to cell death by osmotic shock. Actinoporins appear as multigenic families within the genome of sea anemones: several genes encoding very similar actinoporins are detected within the same species. The Caribbean Sea anemone Stichodactyla helianthus produces three actinoporins (sticholysins I, II and III; StnI, StnII and StnIII) that differ in their toxic potency. For example, StnII is about four-fold more effective than StnI against sheep erythrocytes in causing hemolysis, and both show synergy. However, StnIII, recently discovered in the S. helianthus transcriptome, has not been characterized so far. Here we describe StnIII's spectroscopic and functional properties and show its potential to interact with the other Stns. StnIII seems to maintain the well-preserved fold of all actinoporins, characterized by a high content of ß-sheet, but it is significantly less thermostable. Its functional characterization shows that the critical concentration needed to form active pores is higher than for either StnI or StnII, suggesting differences in behavior when oligomerizing on membrane surfaces. Our results show that StnIII is an interesting and unexpected piece in the puzzle of how this Caribbean Sea anemone species modulates its venomous activity.

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.

House spider genome uncovers evolutionary shifts in the diversity and expression of black widow venom proteins associated with extreme toxicity.

BACKGROUND: Black widow spiders are infamous for their neurotoxic venom, which can cause extreme and long-lasting pain. This unusual venom is dominated by latrotoxins and latrodectins, two protein families virtually unknown outside of the black widow genus Latrodectus, that are difficult to study given the paucity of spider genomes. Using tissue-, sex- and stage-specific expression data, we analyzed the recently sequenced genome of the house spider (Parasteatoda tepidariorum), a close relative of black widows, to investigate latrotoxin and latrodectin diversity, expression and evolution. RESULTS: We discovered at least 47 latrotoxin genes in the house spider genome, many of which are tandem-arrayed. Latrotoxins vary extensively in predicted structural domains and expression, implying their significant functional diversification. Phylogenetic analyses show latrotoxins have substantially duplicated after the Latrodectus/Parasteatoda split and that they are also related to proteins found in endosymbiotic bacteria. Latrodectin genes are less numerous than latrotoxins, but analyses show their recruitment for venom function from neuropeptide hormone genes following duplication, inversion and domain truncation. While latrodectins and other peptides are highly expressed in house spider and black widow venom glands, latrotoxins account for a far smaller percentage of house spider venom gland expression. CONCLUSIONS: The house spider genome sequence provides novel insights into the evolution of venom toxins once considered unique to black widows. Our results greatly expand the size of the latrotoxin gene family, reinforce its narrow phylogenetic distribution, and provide additional evidence for the lateral transfer of latrotoxins between spiders and bacterial endosymbionts. Moreover, we strengthen the evidence for the evolution of latrodectin venom genes from the ecdysozoan Ion Transport Peptide (ITP)/Crustacean Hyperglycemic Hormone (CHH) neuropeptide superfamily. The lower expression of latrotoxins in house spiders relative to black widows, along with the absence of a vertebrate-targeting α-latrotoxin gene in the house spider genome, may account for the extreme potency of black widow venom.

Spit and venom from scytodes spiders: a diverse and distinct cocktail.

Spiders from the family Scytodidae have a unique prey capturing technique: they spit a zig-zagged silken glue to tether prey to a surface. Effectiveness of this sticky mixture is based on a combination of contraction and adhesion, trapping prey until the spider immobilizes it by envenomation and then feeds. We identify components expressed in Scytodes thoracica venom glands using combined transcriptomic and proteomic analyses. These include homologues of toxic proteins astacin metalloproteases and potentially toxic proteins including venom allergen, longistatin, and translationally controlled tumor protein (TCTP). We classify 19 distinct groups of candidate peptide toxins; 13 of these were detected in the venom, making up 35% of the proteome. Six have significant similarity to toxins from spider species spanning mygalomorph and nonhaplogyne araneomorph lineages, suggesting their expression in venom is phylogenetically widespread. Twelve peptide toxin groups have homologues in venom gland transcriptomes of other haplogynes. Of the transcripts, approximately 50% encode glycine-rich peptides that may contribute to sticky fibers in Scytodes spit. Fifty-one percent of the identified venom proteome is a family of proteins that is homologous to sequences from Drosophila sp. and Latrodectus hesperus with uncharacterized function. Characterization of these components holds promise for discovering new functional activity.

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.

Diverse formulas for spider dragline fibers demonstrated by molecular and mechanical characterization of spitting spider silk.

Spider silks have outstanding mechanical properties. Most research has focused on dragline silk proteins (major ampullate spidroins, MaSps) from orb-weaving spiders. Using silk gland expression libraries from the haplogyne spider Scytodes thoracica, we discovered two novel spidroins (S. thoracica fibroin 1 and 2). The amino acid composition of S. thoracica silk glands and dragline fibers suggest that fibroin 1 is the major component of S. thoracica dragline silk. Fibroin 1 is dominated by glycine-alanine motifs, and lacks sequence motifs associated with orb-weaver MaSps. We hypothesize fibroin 2 is a piriform or aciniform silk protein, based on amino acid composition, spigot morphology, and phylogenetic analyses. S. thoracica's dragline silk is less tough than previously reported, but is still comparable to other dragline silks. Our analyses suggest that dragline silk proteins evolved multiple times. This demonstrates that spider dragline silk is more diverse than previously understood, providing alternative high performance silk designs.

Correlation between protein secondary structure and mechanical performance for the ultra-tough dragline silk of Darwin's bark spider.

The spider major ampullate (MA) silk exhibits high tensile strength and extensibility and is typically a blend of MaSp1 and MaSp2 proteins with the latter comprising glycine-proline-glycine-glycine-X repeating motifs that promote extensibility and supercontraction. The MA silk from Darwin's bark spider (Caerostris darwini) is estimated to be two to three times tougher than the MA silk from other spider species. Previous research suggests that a unique MaSp4 protein incorporates proline into a novel glycine-proline-glycine-proline motif and may explain C. darwini MA silk's extraordinary toughness. However, no direct correlation has been made between the silk's molecular structure and its mechanical properties for C. darwini. Here, we correlate the relative protein secondary structure composition of MA silk from C. darwini and four other spider species with mechanical properties before and after supercontraction to understand the effect of the additional MaSp4 protein. Our results demonstrate that C. darwini MA silk possesses a unique protein composition with a lower ratio of helices (31%) and ß-sheets (20%) than other species. Before supercontraction, toughness, modulus and tensile strength correlate with percentages of ß-sheets, unordered or random coiled regions and ß-turns. However, after supercontraction, only modulus and strain at break correlate with percentages of ß-sheets and ß-turns. Our study highlights that additional information including crystal size and crystal and chain orientation is necessary to build a complete structure-property correlation model.

Untangling spider silk evolution with spidroin terminal domains.

BACKGROUND: Spidroins are a unique family of large, structural proteins that make up the bulk of spider silk fibers. Due to the highly variable nature of their repetitive sequences, spidroin evolutionary relationships have principally been determined from their non-repetitive carboxy (C)-terminal domains, though they offer limited character data. The few known spidroin amino (N)-terminal domains have been difficult to obtain, but potentially contain critical phylogenetic information for reconstructing the diversification of spider silks. Here we used silk gland expression data (ESTs) from highly divergent species to evaluate the functional significance and phylogenetic utility of spidroin N-terminal domains. RESULTS: We report 11 additional spidroin N-termini found by sequencing approximately 1,900 silk gland cDNAs from nine spider species that shared a common ancestor > 240 million years ago. In contrast to their hyper-variable repetitive regions, spidroin N-terminal domains have retained striking similarities in sequence identity, predicted secondary structure, and hydrophobicity. Through separate and combined phylogenetic analyses of N-terminal domains and their corresponding C-termini, we find that combined analysis produces the most resolved trees and that N-termini contribute more support and less conflict than the C-termini. These analyses show that paralogs largely group by silk gland type, except for the major ampullate spidroins. Moreover, spidroin structural motifs associated with superior tensile strength arose early in the history of this gene family, whereas a motif conferring greater extensibility convergently evolved in two distantly related paralogs. CONCLUSIONS: A non-repetitive N-terminal domain appears to be a universal attribute of spidroin proteins, likely retained from the origin of spider silk production. Since this time, spidroin N-termini have maintained several features, consistent with this domain playing a key role in silk assembly. Phylogenetic analyses of the conserved N- and C-terminal domains illustrate dramatic radiation of the spidroin gene family, involving extensive duplications, shifts in expression patterns and extreme diversification of repetitive structural sequences that endow spider silks with an unparalleled range of mechanical properties.

Ovarian Transcriptomic Analyses in the Urban Human Health Pest, the Western Black Widow Spider.

Due to their abundance and ability to invade diverse environments, many arthropods have become pests of economic and health concern, especially in urban areas. Transcriptomic analyses of arthropod ovaries have provided insight into life history variation and fecundity, yet there are few studies in spiders despite their diversity within arthropods. Here, we generated a de novo ovarian transcriptome from 10 individuals of the western black widow spider (Latrodectus hesperus), a human health pest of high abundance in urban areas, to conduct comparative ovarian transcriptomic analyses. Biological processes enriched for metabolism-specifically purine, and thiamine metabolic pathways linked to oocyte development-were significantly abundant in L. hesperus. Functional and pathway annotations revealed overlap among diverse arachnid ovarian transcriptomes for highly-conserved genes and those linked to fecundity, such as oocyte maturation in vitellogenin and vitelline membrane outer layer proteins, hormones, and hormone receptors required for ovary development, and regulation of fertility-related genes. Comparative studies across arachnids are greatly needed to understand the evolutionary similarities of the spider ovary, and here, the identification of ovarian proteins in L. hesperus provides potential for understanding how increased fecundity is linked to the success of this urban pest.

Expanding Canonical Spider Silk Properties through a DNA Combinatorial Approach.

The properties of native spider silk vary within and across species due to the presence of different genes containing conserved repetitive core domains encoding a variety of silk proteins. Previous studies seeking to understand the function and material properties of these domains focused primarily on the analysis of dragline silk proteins, MaSp1 and MaSp2. Our work seeks to broaden the mechanical properties of silk-based biomaterials by establishing two libraries containing genes from the repetitive core region of the native Latrodectus hesperus silk genome (Library A: genes masp1, masp2, tusp1, acsp1; Library B: genes acsp1, pysp1, misp1, flag). The expressed and purified proteins were analyzed through Fourier Transform Infrared Spectrometry (FTIR). Some of these new proteins revealed a higher portion of ß-sheet content in recombinant proteins produced from gene constructs containing a combination of masp1/masp2 and acsp1/tusp1 genes than recombinant proteins which consisted solely of dragline silk genes (Library A). A higher portion of ß-turn and random coil content was identified in recombinant proteins from pysp1 and flag genes (Library B). Mechanical characterization of selected proteins purified from Library A and Library B formed into films was assessed by Atomic Force Microscopy (AFM) and suggested Library A recombinant proteins had higher elastic moduli when compared to Library B recombinant proteins. Both libraries had higher elastic moduli when compared to native spider silk proteins. The preliminary approach demonstrated here suggests that repetitive core regions of the aforementioned genes can be used as building blocks for new silk-based biomaterials with varying mechanical properties.

Diversity despite dispersal: colonization history and phylogeography of Hawaiian crab spiders inferred from multilocus genetic data.

The Hawaiian archipelago is often cited as the premier setting to study biological diversification, yet the evolution and phylogeography of much of its biota remain poorly understood. We investigated crab spiders (Thomisidae, Mecaphesa) that demonstrate contradictory tendencies: (i) dramatic ecological diversity within the Hawaiian Islands, and (ii) accompanying widespread distribution of many species across the archipelago. We used mitochondrial and nuclear genetic data sampled across six islands to generate phylogenetic hypotheses for Mecaphesa species and populations, and included penalized likelihood molecular clock analyses to estimate arrival times on the different islands. We found that 17 of 18 Hawaiian Mecaphesa species were monophyletic and most closely related to thomisids from the Marquesas and Society Islands. Our results indicate that the Hawaiian species evolved from either one or two colonization events to the archipelago. Estimated divergence dates suggested that thomisids may have colonized the Hawaiian Islands as early as ~10 million years ago, but biogeographic analyses implied that the initial diversification of this group was restricted to the younger island of Oahu, followed by back-colonizations to older islands. Within the Hawaiian radiation, our data revealed several well-supported genetically distinct terminal clades corresponding to species previously delimited by morphological taxonomy. Many of these species are codistributed across multiple Hawaiian Islands and some exhibit genetic structure consistent with stepwise colonization of islands following their formation. These results indicate that dispersal has been sufficiently limited to allow extensive ecological diversification, yet frequent enough that interisland migration is more common than speciation.

Alternative Transcription at Venom Genes and Its Role as a Complementary Mechanism for the Generation of Venom Complexity in the Common House Spider.

The complex composition of venom, a proteinaceous secretion used by diverse animal groups for predation or defense, is typically viewed as being driven by gene duplication in conjunction with positive selection, leading to large families of diversified toxins with selective venom gland expression. Yet, the production of alternative transcripts at venom genes is often overlooked as another potentially important process that could contribute proteins to venom, and requires comprehensive datasets integrating genome and transcriptome sequences together with proteomic characterization of venom to be fully documented. In the common house spider, Parasteatoda tepidariorum, we used RNA sequencing of four tissue types in conjunction with the sequenced genome to provide a comprehensive transcriptome annotation. We also used mass spectrometry to identify a minimum of 99 distinct proteins in P tepidariorum venom, including at least 33 latrotoxins, pore-forming neurotoxins shared with the confamilial black widow. We found that venom proteins are much more likely to come from multiple transcript genes, whose transcripts produced distinct protein sequences. The presence of multiple distinct proteins in venom from transcripts at individual genes was confirmed for eight loci by mass spectrometry, and is possible at 21 others. Alternative transcripts from the same gene, whether encoding or not encoding a protein found in venom, showed a range of expression patterns, but were not necessarily restricted to the venom gland. However, approximately half of venom protein encoding transcripts were found among the 1,318 transcripts with strongly venom gland biased expression. Our findings revealed an important role for alternative transcription in generating venom protein complexity and expanded the traditional model of venom evolution.

The transcriptome of Darwin's bark spider silk glands predicts proteins contributing to dragline silk toughness.

Darwin's bark spider (Caerostris darwini) produces giant orb webs from dragline silk that can be twice as tough as other silks, making it the toughest biological material. This extreme toughness comes from increased extensibility relative to other draglines. We show C. darwini dragline-producing major ampullate (MA) glands highly express a novel silk gene transcript (MaSp4) encoding a protein that diverges markedly from closely related proteins and contains abundant proline, known to confer silk extensibility, in a unique GPGPQ amino acid motif. This suggests C. darwini evolved distinct proteins that may have increased its dragline's toughness, enabling giant webs. Caerostris darwini's MA spinning ducts also appear unusually long, potentially facilitating alignment of silk proteins into extremely tough fibers. Thus, a suite of novel traits from the level of genes to spinning physiology to silk biomechanics are associated with the unique ecology of Darwin's bark spider, presenting innovative designs for engineering biomaterials.

Recent progress and prospects for advancing arachnid genomics.

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

Stichodactyla helianthus' de novo transcriptome assembly: Discovery of a new actinoporin isoform.

Transcriptomic profiling of venom producing tissues from different animals is an effective approach for discovering new toxins useful in biotechnological and pharmaceutical applications, as well in evolutionary comparative studies of venomous animals. Stichodactyla helianthus is a Caribbean sea anemone which produces actinoporins as part of its toxic venom. This family of pore forming toxins is multigenic and at least two different isoforms, encoded by separate genes, are produced by S. helianthus. These isoforms, sticholysins I and II, share 93% amino acid identity but differ in their pore forming activity and act synergistically. This observation suggests that other actinoporin isoforms, if present in the venomous mixture, could offer an advantageous strategy to modulate whole venom activity. Using high-throughput sequencing we generated a de novo transcriptome of S. helianthus and determined the relative expression of assembled transcripts using RNA-Seq to better characterize components of this species' venom, focusing on actinoporin diversity. Applying this approach, we have discovered at least one new actinoporin variant from S. helianthus in addition to several other putative venom components.