The evolutionary history of cribellate orb-weaver capture thread spidroins
BMC Ecology and Evolution volume 22, Article number: 89 (2022)
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.
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.
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.
Spiders are known for making a variety of task-specific silk fibers with diverse material properties . The evolutionary success of spiders can be attributed to their wide ranging use of silk. For example, orb-web weavers use multiple silks in web construction, with each type synthesized in a morphologically specialized silk gland. The combination of different silk types makes the orb-web a strong yet extensible prey-catching system [2, 3]. The web frame and radii are primarily made of major ampullate silk, the attachment disks are made of pyriform silk, and the capture spiral is a composite of a stretchy filament and sticky silks. There are two types of capture spirals. One type has a pair of proteinaceous fibers produced by flagelliform glands that are coated with moist, adhesive droplets produced by aggregate glands . This type of capture spiral is spun by over 4,600 species of ecribellate spiders in seven spider families (superfamily Araneoidea), the largest and most commonly encountered being the family Araneidae . The protein adhesive in araneid capture lines is also found in the gumfoot lines of cobwebs which are constructed by the ~ 2,800 species of Theridiidae and Nesticidae [6, 7]. In contrast, cribellate capture threads are supported by a pair of pseudoflagelliform fibers, which are covered with tangles of dry, mechanically sticky cribellate nanofibrils [8,9,10,11]. Cribellate threads also are spun by various basal members of the infraorder Araneomorphae, which encompasses the largest number of living spider species. However, members of only 21 of the 96 families in Araneomorphae continue to spin cribellate threads . Among orb-web building spiders, the transition from dry cribellate to viscous sticky capture threads is not well understood and whether cribellate and ecribellate webs evolved in parallel [13,14,15,16,17,18] or cribellate silk was replaced by sticky lines [19, 20] is still uncertain.
Among the cribellate spiders are members of the family Uloboridae. These small to medium-sized spiders are one of only two spider families to lack venom glands [21, 22]. The family Uloboridae currently encompasses 19 genera and 287 species distributed worldwide . Uloborids are known to have a laborious prey wrapping technique that results in tight swaddling, thought to be a way of compensating for the lack of venom to immobilize prey [10, 23,24,25]. Most uloborid species construct fully circular orb-webs while other species construct reduced orb-webs, such as the triangle-webs of Hyptiotes or the single-line capture threads of Miagrammopes [10, 26].
Cribellate silk is produced by a very complex and specific spinning process; thousands of ultrafine fibrils are spun simultaneously from a plate-like spinning organ with numerous spigots, called the cribellum, onto an extensible core fiber of pseudoflagelliform silk. The fibrils are vigorously combed into woolly snarls by a row of specialized leg bristles—the calamistrum [27,28,29]. In the uloborid capture spiral, the cribellar fibrils surround the supporting axial fibers [8, 9, 30]. The capture spiral also possesses connecting fibers produced in paracribellate glands [31, 32]. Non-uloborid cribellate species can have a more complex cribellate silk that includes highly folded reserve warp fibers (also called undulating fibers) that are likely produced in minor ampullate glands [4, 30, 33,34,35]. Cribellate capture threads can stick to many different surfaces and they achieve this adhesion via a combination of hygroscopic forces, van der Waals’ forces, entanglement of cuticular structures, and the absorption of cuticular waxes from prey [36,37,38,39]. The mechanical properties of capture thread are influenced by the use of different fibers, making cribellate silk both stiff and extensible. While the pseudoflagelliform axial fibers provide stiffness, strength, and initial elasticity, the irreversible unfolding of cribellar fibrils allow the capture thread to extend as much as 500% relative to its original length .
The mechanical behavior of spider silk has also been attributed to its protein composition. Spider silks are mainly composed of silk proteins called spidroins (a contraction of spider fibroin; ). Spidroins are encoded by a single gene family and are highly expressed in silk glands [42,43,44,45,46,47,48,49]. While silk proteins for most silk types have been described, the molecular composition of cribellate threads has only recently gained attention. Sequences from several cribellate spiders show that cribellar fibrils and pseudoflagelliform fibers are composed of distinct spidroins specific to cribellate and pseudoflagelliform glands, respectively [49, 50]. Yet, detailed examination of the genetics of cribellate silk and how particular genes relate to the evolution of cribellate silk components is not well understood. For instance, candidates for the paracribellar spidroins, the main components of paracribellar connecting fibers, have not been identified. Here, we describe a comprehensive spidroin set for the uloborid orb-weaver Uloborus diversus (Walckenaer, 1841) using genomic and transcriptomic data and examine, in detail, the spidroins associated with the three types of silk in cribellate capture threads: cribellar fibrils, pseudoflagelliform fibers, and paracribellar fibers. We present the complete sequence of pseudoflagelliform spidroin and identify potential spidroin candidates for paracribellar silk. We also show that pseudoflagelliform spidroins are closely related to the flagelliform spidroins, suggesting a shared origin of the two fiber types.
Uloborus diversus genome and transcriptome assembly
We obtained a total of 1571 billion reads from the genome sequencing of U. diversus. These reads were assembled using Supernova 2.0.1 (10X Genomics) and the resulting genome assembly was 1.49 Gb and included 47,680 scaffolds greater than 10 kb (scaffold N50: 42.36 kb). Based on published flow cytometry data, one Uloborus species has a genome size of approximately 2.5 Gb , suggesting our assembly is incomplete. Nevertheless, the assembled genome was found to have 90% of the Benchmarking Universal Single-Copy Orthologs (BUSCO ) genes as either complete (65%) or fragmented (25%). In total, 17,878 high quality Iso-Seq transcripts were obtained with a mean length of 3980.6 bp. For the RNA-Seq assembly, an average of 23.5 and 23.2 million cleaned reads were obtained from the two total silk libraries and two cephalothorax libraries, respectively.
The de novo assembly of the RNA-Seq reads had 542,127 transcripts with an N50 of 794. Since Trinity transcriptomes can contain both multiple alleles per gene and potentially low-quality transcripts, we developed a novel trimming and annotation pipeline called TrTAP to reduce the de novo transcriptome to a set of high confidence genes (Additional file 3: Figure S1). The pipeline begins by comparing the transcripts to multiple gene-specific databases (here spidroins) or to high-quality complete gene sets from related species using BLAST. The pipeline then parses the BLAST results to select one representative allele per Trinity cluster as well as exclude chimeric and RNA transcripts. For the clusters where there are no BLAST matches, the longest ORF is used in the selection, with such ORFs greater than 50 amino acids. The pipeline annotates all allele clusters using comparisons to PFAM and SwissProt and calculates their expression abundance. In the final step, the pipeline removes transcripts that are not best hits to the gene databases and have low coverage of matching genes or low expression using RSEM Transcripts per Million (TPM). The resulting pruned transcriptome is composed of 115,961 proteins, of which 26,088 had homologous mapping to known genes. This pruned set contains a majority of the expression and annotated genes as well as differentially expressed genes (Additional file 2: Table S1). Single copy orthologs analysis, BUSCO , indicates that the trimmed transcriptome is 91% complete. The pruned set has matches to 87.4% of the BUSCO genes. The missing 3.6% BUSCO genes not found in the pruned set matched to genes removed for having low expression values.
Spidroin complement and expression in U. diversus
We identified 18 spidroin genes expressed in U. diversus silk glands. For some spidroin genes, full length nucleotide sequences were obtained from single 10X contigs, for those that were not, reconstruction was done using a combination of 10X contigs and Iso-Seq reads. Full or partial length spidroin genes were assembled for all expected spidroin types. Aciniform (AcSp), tubuliform (TuSp), minor ampullate (MiSp), major ampullate (MaSp), pseudoflagelliform (Pflag) spidroins, and two spidroins that could not be assigned to a known spidroin type, annotated with the general name Spidroin and a version letter (Sp_vA and Sp_vB), all appear to be complete, full-length genes. The U. diversus ampullate spidroin (U. div_AmSp) also appears to be full-length but could not be categorized as either MaSp or MiSp. The region bounded by the terminal domains of U. div_AmSp is not repetitive and lacks the characteristic amino acid motifs of MaSp or MiSp, but the terminal regions share conserved spidroin amino acid signatures (Additional file 3: Figure S2). The N- terminal region has the conserved amino acid motifs found to be restricted to MaSp spidroins (e.g. amino acid motifs PW, AxxFxxxxF from ) and the C-terminal region has the charged amino acid residues that form salt bridges and aid in monomer stabilization on MaSp and MiSp spidroins [53,54,55]. We recovered all of the previously published partial spidroin sequences for U. diversus by Garb and colleagues [56, 57]. Spidroin annotation was based on spidroin-specific gene characteristics, such as the presence of spidroin-specific motifs in the terminal domains (see Collin et al. ) and repeat composition. Since spidroins are members of a gene family, annotation was also supported by evolutionary relationship information from spidroin gene trees (Fig. 1).
The set of spidroin genes in U. diversus is consistent with the set recently described for six species in the uloborid genus Octonoba by Kono et al. . The amino acid composition, approximate length, and organization of U. diversus spidroins are consistent not only with those of Octonoba but also with the spidroins of non-uloborid spiders [49, 57,58,59]. For example, the complete sequence of AcSp1 is composed of ten tandem repeats that are on average 357 amino acids long and account for ~ 91% of the protein (Additional file 3: Figure S2). Like previously described AcSp sequences, U. diversus AcSp1 is enriched in the amino acids serine (25%), alanine (14%), and glycine (8%). Phylogenetic analyses of N- and C- terminal coding regions show Uloborus and Octonoba spidroin sequences cluster together with high bootstrap support (e.g. ≥ 99%) within clades of the same type of spidroin from other species (e.g., U. div and O. syb AcSp1 in an AcSp1 clade; Fig. 1 orange box).
Confirmation of U. diversus spidroin functionality was obtained from expression analyses using transcripts derived from total silk gland and cephalothorax tissues (Additional file 3: Figure S3). We found that all identified spidroin genes are expressed in silk glands (relative to the cephalothorax; DESeq analysis, padj < 0.001), most at high levels (Additional file 3: Figure S3). Moreover, the type, number, and expression of spidroins in U. diversus is similar to the silk genes expressed in Octonoba species . Aside from spidroin genes, there are 1,357 genes that are up-upregulated in uloborid silk glands compared to cephalothorax. The functions of these proteins suggest they are involved in oxidoreductase and transmembrane transport, similar to what has been described in theridiid and araneid species [48, 60, 61]. Given the functional similarities of the dry adhesive used by cribellate spiders and the wet glue in ecribellate orb-web and cobweb weavers, we searched for genes in U. diversus genome with sequence similarly to Theridiidae aggregate gland specific OESTs (Over Expressed Sequence Transcripts ). These aggregate OESTs are thought to be important in the production of aggregate glue droplets. We found 21 U. diversus genes with a match to 11 OESTs (e-values ≤ e−50). However, mapping RNA-seq reads to these 21 genes shows that they are lowly expressed in both total silk and cephalothorax tissues (Additional file 3: Table S2). These results suggest that the U. diversus genome has homologs to only a fraction of the of cobweb weaver aggregate-specific OESTs (> 300 per species), and that these do not have similar roles to cobweb weaver proteins.
Cribellate capture thread spidroins in Uloboridae
Uloborid capture threads are composed of three different silk types, thus we expected to find genes corresponding to those silks: cribellar spidroin (CrSp), pseudoflagelliform spidroin (Pflag), and at least one spidroin candidate for paracribellar silk, which has yet to be identified in any species [49, 50]. Indeed, we found gene sequences that match the cribellar spidroin CrSp (Fig. 2). The repetitive region of U. diversus CrSp was found to have three distinct types of repeat modules (Fig. 2B). Each module type has a distinctive amino acid combination and is repeated multiple times within the repetitive region of U. diversus CrSp (Fig. 2A, B). These repeat modules are not unique to U. diversus, they are also present in CrSp orthologs from other cribellate spiders such as Octonoba, Tengella, Stegodyphus, and Badumna [49, 50, 59]. For example, repeat module 1 has a 79% amino acid (aa) identity within Uloboridae and a range of 28–34% aa identity between uloborids and other families, suggesting it has been conserved for at least ~ 250 Mya (Additional file 3: Figure S4) [17, 62]. Repeat modules 2 and 3 do not align well outside of Uloboridae while having 75% and 69% identity within Uloboridae, respectively (Additional file 3: Figure S4). Comparison of CrSp repeat structure across species shows that unlike previously reported by Kono et al. , Octonoba CrSp is also composed of three repeat units that are homologous to those of U. diversus CrSp but differ significantly at the higher level repeat arrangement (Additional file 3: Figure S4A). The terminal regions of cribellar spidroins show conservation across species, with N- and C-terminal regions sharing 56% and 76% identity at the aa level, respectively (Fig. 2C, D). Phylogenetic analysis further grouped CrSp sequences into a well-supported monophyletic clade (Fig. 1 and 99% BT).
The full-length U. diversus Pflag gene is approximately 8 kbp in size and codes for an ~ 2600 aa long protein (Fig. 3). The repetitive sequence of Pflag has 48 ensemble repeats that cluster into four distinct repeat types (Additional file 3: Figure S5). Each repeat has an 11–13 aa spacer motif, which is a highly conserved glycine-poor region that interrupts glycine-rich motifs , followed by variations of the aa motif XPSSGGXGGXEK, where X can be the amino acids: G, S, A, Q, E, or D (Fig. 3). These repeats are further organized in recurring 4-mer repeats (Fig. 3). Comparative analysis of the repeat structure of Pflag led us to define multiple ensemble repeats in Octonoba Pflag, which is different than the single unit described by Kono et al.  (Additional file 3: Figure S6). The sequence similarity among the multiple ensemble repeats in Octonoba Pflag is not as strongly correlated as those of CrSp, yet it shows a 4-mer repeat in the higher repeat structure (Additional file 3: Figure S6). Unlike the flagelliform spidroin gene (Flag) from ecribellate orb-web weavers [63,64,65], uloborid Pflag does not appear to have introns. Furthermore, we found Pflag sequences to have a lower percentage of glycine and proline amino acids (G: ~38%, P: ~11% for U. diversus and O. sybotides) compared to Flag sequences (G: ~50%, P: ~14% for A. ventricosus and T. clavipes). Moreover, Pflag amino acid motifs are not the same as in Flag; instead Pflag has the proline-containing motifs GPS(X), KPS(X), and QPS(X), with KPS(X) being the dominant motif in U. diversus Pflag (Fig. 3).
The Pflag N-terminal domain is most closely related to Sp5803 (Fig. 1), a spidroin that is expressed in flagelliform silk glands (but is not Flag) and that also lacks introns . The spidroin Sp5803, described in the orb-web weaving Trichonephila clavipes, has a unique repetitive region relative to all other described spidroin types and lacks the conserved spidroin C-terminal region [64, 66]. BLAST searches to published spider genomes, revealed the presence of Sp5803 in Araneus ventricosus (Accession GBL75419.1) based on the presence of the conserved spidroin N-terminal region, similar repetitive region to T. clavipes Sp5803, and absence of conserved spidroin C-terminal region. Phylogenetic analyses show that U. diversus Pflag is in a clade with Pflag spidroins from other cribellate spiders (Figs. 1, 4). Moreover, we found phylogenetic evidence for ecribellate orb-web weaving Flag spidroins to share a most recent common ancestor with Pflag/Sp5803 (69% BT), which together are closely related to CrSp (Figs. 1, 4).
The molecular components of paracribellar silk are unknown. We propose two candidate paracribellar spidroins: U. div Sp_vA and U. div Sp_vB. These spidroins appear to be full-length and have divergent repetitive sequence from each other and other spidroin types. While U. div_Sp_vA and U. div_Sp_vB possess the conserved spidroin terminal domains that include positionally conserved amino acid motifs , they did not group with well-established silk gene types in our phylogenetic analysis, suggesting they may represent new spidroin types. The repetitive regions of U. div_Sp_vA and U. div_Sp_vB, like known spidroins, are organized into tandemly arranged repeat units. U. div_Sp_vA has eight repeat units that are highly similar to each other with 86% identity at the aa level, and the repeats are enriched for the amino acids alanine (21%), serine (16%), glycine (11%) and valine (10%) (Additional file 3: Figure S2). U. div_Sp_vB repeat region is rich in serine (17%) and alanine (12%) amino acids; it has a 28 aa repeat module that is repeated five times, followed by a region with no recognizable modules (Additional file 3: Figure S2). Considering overall amino acid composition, U. div_Sp_vA and U. div_Sp_vB are similar to tubuliform and aciniform spidroins, but the repeat modules do not align well to these other paralog groups.
We examined spigot morphological data within a phylogenetic framework to investigate how capture thread spidroins are associated with evolutionary changes in spigot types. Gene family analyses of concatenated N- and C- terminal nucleotide sequences of spidroins genes involved or thought to be involved in capture thread production, show that U. div_Sp_vA and U. div_Sp_vB cluster within a clade with other “Sp”. Other spidroins in that clade also have divergent repetitive regions and have not been assigned a spidroin type (Fig. 4). Unexpectedly, D. spinosa spidroin grouped with U. div_Sp_vB, likely a result of D. spi_Sp sequence lacking an N-terminal region and having a divergent C-terminal region that has similarities to that of U. div_Sp_vB. Moreover, unlike U. div_Sp_vB, the repetitive region of D. spi_Sp is like other Pflag and Flag sequences . Two transcripts containing the N-terminal region only, B. lon_Sp_NvA and T. per_Sp_N grouped together as a sister clade of Flag (Figs. 4 and 57% BT), a grouping consistent with previous spidroin gene tree analysis . The available repetitive regions of T. per_Sp_N and B. lon_Sp_NvA are short, enriched in the amino acids glycine (22% and 23% for T. per_Sp_N and B. lon_Sp_NvA, respectively), glutamic acid (10% and 11%), and serine (9% and 10%), and lack amino acid motifs corresponding to Flag or Pflag spidroins.
Cribellate silk is spun by thousands of araneomorph spider species with a wide ecological and taxonomic distribution [68, 69]. Our understanding of spider silks and silk genetics is largely based on ecribellate orb-web weavers, leaving cribellate silk genetics relatively unstudied. Our identification and annotation of the silk gene catalog of U. diversus greatly increases our knowledge of spidroin diversity in cribellate species and sheds light into the evolution of cribellate silk genes. Sequence similarity and phylogenetic affinity of U. diversus spidroins to those of Octonoba species suggest that the ancestor of these uloborids had a similar complement of silk proteins likely dating back to ~ 145 Mya . Spidroin sequences not related to the capture threads (e.g. AcSp, TuSp, PySp) show conservation in their terminal regions, as well as in repetitive regions, with previously described sequences [45, 49, 50, 64, 71,72,73,74]. This sequence conservation suggests similar selective pressures have acted on these spidroins in diverse spider lineages. Uloborid ampullate sequences form monophyletic clades with moderate support (U. div + O. syb MaSp1, MaSp2, and MiSp; Fig. 1), within a large, diverse MaSp and MiSp spidroin clade (Fig. 1). The grouping of uloborid MaSp1 and MaSp2 to the exclusion of MaSp1 and MaSp2 of araneoid spiders suggests gene conversion of the terminal domains in the ancestor of uloborids, similar to what has been proposed for black widows . It is also possible that extensive convergent evolution of MaSp1-like and MaSp2-like repetitive regions has occurred in divergent orb-web weaver clades.
Since cribellate silk production is ancestral for the hyper-diverse infraorder Araneomorphae, a fundamental knowledge of the molecular constituents of cribellate capture thread is pivotal to our understanding of spider silk evolution. Uloborid cribellate silk is a composite of multiple fiber types that are produced in distinct silk glands: the core axial fibers are produced in a pair of pseudoflagelliform silk glands, cribellar nanofibers that form the “puffs” of cribellate silk are produced in several thousand cribellar silk glands, and paracribellar silk, which connects nanofibers to the axial fiber, come from ~ 20 paracribellate silk glands [8, 10, 11, 76]. The combination of these fibers gives cribellate silk its woolly yarn appearance, adhesiveness, and mechanical properties [37,38,39,40, 77,78,79].
Flagelliform silk is the main component of the core axial fibers of capture threads in ecribellate spiders and is produced in flagelliform silk glands . Pseudoflagelliform silk of cribellate orb-web weaving spiders is thought to be homologous to flagelliform silk of ecribellate orb-web and cob-web weaving spiders because the spigots that produce these fibers are located on the same pair of spinnerets in both groups. That is, they have positional homology [19, 30, 80, 81]. It follows that the glands associated with those spigots are also homologous. Furthermore, under the glandular affiliation hypothesis, we would expect proteins to have evolved in association with the silk gland where they are mainly expressed . Gene tree analyses show that spidroins associated with core axial fibers are within a clade of Flag and Pflag plus Sp5803 spidroins (Figs. 1 and 4). The well-supported grouping of uloborid Pflag spidroins with Sp5803 and S. mimosarum Sp2a suggests these latter spidroins are, in fact Pflag orthologs. We also propose that the D. spinosa spidroin (ABD61590) described by Garb et al.  is likely to be expressed in pseudoflagelliform silk glands, but is neither a Pflag ortholog nor a Flag spidroin as previously thought. In terms of Flag and Pflag evolution, we posit that these two spidroin types are the results of a duplication event in the common ancestor of araneoids, uloborids, the RTA-clade, and eresids (Entelegynae), approximately 250 Mya [17, 62]. Flag appears to have been lost in uloborids and eresids, but retained in at least two RTA-clade families, Zoropsidae and Desidae (Fig. 4). The phylogenetic relationship of Pflag and Sp5803 suggests that this spidroin type was maintained in the evolution of araneoid spiders but was subsequently lost in cobweb weavers (Theridiidae); at least searches of the Parasteatoda tepidariorum genome and L. hesperus transcriptomes have failed to recover orthologs of Pflag/Sp5803 [48, 60, 82]. In the orb-web weaving spider T. clavipes, Sp5803 is highly expressed in flagelliform silk glands and thought to be associated with capture webs . Theridiid cobwebs are highly modified relative to orb-webs, and while they have flagelliform glands, in most species glue droplets are deposited on major ampullate silk fibers, suggesting less selective pressures on genes associated with flagelliform fibers [83, 84].
Capture threads of cribellate and ecribellate spiders serve the same function, prey retention, but they achieve it using different mechanisms. Cribellate capture thread is a dry adhesive that uses a combination of hygroscopic forces, van der Waals’ forces, entanglement of cuticular structurers (e.g., setae and spines), and the absorption of cuticular waxes of prey [11, 36, 38]. By contrast, ecribellate capture thread is a wet adhesive that uses viscoelastic glue droplets composed of glycoproteins and phosphoproteins surrounded by aqueous coats of salts [84,85,86,87,88]. Unlike the spigots of the core axial fiber (see above), spigots associated with the production of cribellar fibrils and aqueous glue droplets in cribellate and ecribellate capture threads are not homologous . In terms of amino acid repeat structure, we found no similarity in the silk genes associated with cribellar fibrils or aggregate droplets, CrSp and AgSp respectively. Moreover, phylogenetic analyses show that while CrSp and AgSp each form monophyletic clades, they are not each other’s closest relatives (Figs. 1 and 4). AgSp appears to result from a duplication event in araneoid ecribellate spiders that is independent from the origin of CrSp.
In Uloboridae, paracribellar silk connects the core axial fibers to the surrounding cribellar nanofibers [4, 8, 89]. Paracribellar spigots have been identified in some cribellate taxa and are usually located on the posterior median spinneret  (Fig. 4B, C). The main proteins associated with this silk type have not been identified, in part due to lack of silk-oriented studies in cribellate spider taxa. From the U. diversus spidroin catalog, we propose U. div Sp_vA and U. div Sp_vB as candidate genes of paracribellar silk based on their phylogenetic relationships. Spidroin types almost always form well supported monophyletic groups (e.g. all TuSp together, all AcSp together, etc.; Fig. 1; [44, 49, 57]). In contrast, U. div Sp_vA and U. div Sp_vB, failed to group with any known spidroin type, and instead both U. div Sp_vA and U. div Sp_vB clustered weakly outside of the Flag, Pflag/Sp5803, and CrSp clade (Fig. 1). To better understand the origins of U. div Sp_vA and U. div Sp_vB sequences, spidroins exclusively involved in capture thread production, including “Sp” spidroins from other cribellate species, were included in a gene tree analyses (Fig. 4). The two uloborid sequences, U. div Sp_vA and U. div Sp_vB, are positioned within a clade that includes other spidroins with unknown glandular origin (Fig. 4). Because these sequences do not share the characteristics of CrSp or Pflag spidroins and share a monophyletic origin, we propose that spidroins within U. div Sp_vA and U. div Sp_vB containing clade are likely to be expressed in paracribellate silk glands. While repetitive sequences within this clade are divergent, this could be related to a mechanism that connects the axial fiber and the cribellar fibrils to each other or a strategy for dealing with different environments. A wider sampling of cribellate spider genomes and expression studies of paracribellar glands are needed to obtain a more specific annotation of “Sp” spidroins.
Differences in the mechanical performance of spider silks have been attributed to fiber number and arrangement [33, 34], spinning rate , humidity [79, 91], temperature , and diet . The molecular components of silk, especially spidroins, are also known to influence mechanical performance. The mechanical properties of cribellate capture thread have been recorded for a few species, showing it can stretch up to 500% its original length [40, 77, 93]. This extensibility is achieved because of the composite design of cribellate capture thread, with pseudoflagelliform axial fibers contributing the initial strength and extension (~ 150%) and after its breakage, the multiple cribellar fibrils unfold to allow exceptional extensibility. Ecribellate capture silk can stretch as much as 1000% when covered in aqueous glue [94, 95]. This extensibility is attributed to the presence of the tandemly arrayed amino acid motif GPG(X) in Flag that forms spring-like helices , with a greater number of GPG(X) motifs associated with more extensible fibers . Compared to native wet araneoid silks, dry cribellate capture silks tend to have lower extensibilities and higher strengths [77, 96]. Based on these performance comparisons, it was suggested that the core fiber of cribellate capture thread spidroin (Pflag) would have similar amino acid motifs to those of ecribellate Flag spidroins, but fewer instances of proline-containing motifs, or would express different spidroin paralogs to achieve extensibility intermediate between other fiber types and araneoid wet capture silk . Consistent with both hypotheses, Flag sequences indeed have a higher percentage of glycine and proline amino acids compared to Pflag sequences. Uloborid Pflag spidroin does not have the same amino acid motifs as Flag, instead, it has the motifs GPS(X), KPS(X), and QPS(X), with KPS(X) being the dominant motif in U. diversus Pflag (Fig. 3). These proline-containing amino acid motifs are also different from the Deinopis spinosa spidroin motif GPQ(X) .
The differences in amino acid motifs in Flag and Pflag may also partly explain differences in how axial fibers of ecribellate and cribellate capture spirals respond to water. Specifically, ecribellate flagelliform fibers are more than twice as extensible (> 100% higher) when covered in aqueous glue than when cleaned of aqueous glue [96,97,98]. In contrast, uloborid (cribellate) pseudoflagelliform fibers only increase extensibility by 20–30% when covered in water , which is more similar to proline-poor silks (i.e., minor ampullate and major ampullate silk of some species) than to proline-rich flagelliform, even though Pflag has a much higher percentage of proline than most minor ampullate spidroins . The presence of glycine in the GPGX motif of Flag versus the higher prevalence of glycine-lacking motifs in uloborid Pflag may contribute to their differences in response to water and maximum extensibility. Consistent with that idea, D. spinosa capture threads stretch significantly more than those of Uloboridae , and have a higher percentage of glycine containing GPS(X) motifs in their spidroin (ABD61590) than uloborids. However, the additional undulating fibers present in D. spinosa cribellate capture thread  may be more important for its greater extensibility than abundance of GPS(X) motifs. Extensibility differences across cribellate and ecribellate capture threads are likely due to variation and abundance of specific amino acid motifs. Within cribellate species, it seems that different lineages have developed different strategies to achieve stretchiness, such as the addition of undulating fibers and the use of multiple diverse motifs in Pflag spidroins.
The overall structure of CrSp and Pflag spidroins in U. diversus highlights the complexity of repeat organization within these genes. Both genes are comprised of distinct repeats arranged in a stereotypical pattern. The occurrence of hierarchically organized repetitive units (i.e., ensemble repeats) has been found in other spidroins such as MaSp1  and AgSp [99, 100] suggesting they may be critical to the mechanical performance of these proteins. Similarly, comparison of the U. diversus sequences with Octonoba reveals shared repetitive features that may have functional significance. For CrSp, there is clear homology between species in each of the primary repeat modules but striking differences in the overall arrangement of these modules. For Pflag, there is higher repetitive unit sequence divergence among species, but similar ensemble repeat structure. Collecting full-length CrSp and Pflag sequence for additional species within the Uloboridae will be critical to mapping the pattern of evolutionary changes in repeat organization in these genes and to understanding how these changes impact capture thread performance.
Cribellate capture silk is a complex mixture of many fibers that work as a unit to adhere to prey. How these fibers are woven together affects the mechanical properties as each fiber has a distinct biophysical behavior. We described in detail the molecular components of each silk type for U. diversus, including a newly discovered spidroin type. For each silk type, we analyzed associated spidroins and elucidated their role in the mechanical properties of cribellate silk. Furthermore, we showed that the molecular composition of the capture thread of cribellate spiders has a complex evolutionary history involving at least six duplication events dating back to ~ 250 mya [17, 62]. While Pflag and Flag spidroins resulted from a more recent duplication and likely were co-expressed ancestrally, CrSp spidroins do not appear homologous to AgSp (do not share a recent gene duplication event) suggesting both punctuated and gradual evolution of the capture spiral. Given that cribellate threads are mainly spun by the basal members of the infraorder Araneomorphae, further knowledge of the underlying molecular composition of more primitive cribellate capture thread will surely shed light into the evolution of capture thread fibers and adhesives.
Materials and methods
Sample preparation and nucleic acid extractions
Mature female U. diversus were collected in Riverside, California, USA by T. Dugger. High molecular weight genomic DNA was extracted from a single whole individual using the Gentra PureGene tissue kit (Qiagen, Valencia, CA, USA). For RNA, cephalothorax and total silk gland tissue (a combination of all silk gland types attached to the spinnerets) were isolated from individual spiders, flash frozen in liquid nitrogen, and stored at − 80 °C. Tissues were homogenized in TRIzol reagent using a Fisherbrand bead mill 24 homogenizer (Fisher Scientific, Waltham, MA, USA). RNA was then purified using the PureLink RNA Mini kit with on-column DNAse treatment (Ambion, ThermoFisher Scientific, Wilmington, MA, USA). Because of the small size of U. diversus, for each RNA extraction, tissues of the same type were combined from four individuals prior to homogenization. Nucleic acid quantification was done using a Qubit Fluorometer (Thermo Fisher Scientific, Wilmington, MA, USA) and RNA integrity was assessed with a Bioanalyzer (Agilent Technologies, Santa Clara CA, USA).
Genomic library preparation and sequencing
A total of 950 ng of high molecular weight genomic DNA extracted from a single individual was sent to the New York Genome Center (New York, NY, USA) for 10× Chromium library preparation (10X Genomics, Pleasanton, CA, USA), which included size selection (with PippinHT at 30 kb) and quality control. The library was paired-end sequenced (150 base pairs) on one lane of an Illumina HiSeqX.
Iso-Seq library preparation and sequencing
Iso-Seq libraries were constructed from a total silk gland tissue RNA extraction following the Pacific Bioscience’s Iso-Seq protocol using the SMARTer PCR cDNA Synthesis kit (Clontech, Mountain View, CA, USA). Amplification of cDNA was done with 16 cycles for 4 reactions and 18 cycles for another four reactions. Amplification products were combined and purified using normal-phase 2× AMPure XP (Beckman Coulter, Brea, CA, USA). Cleaned cDNA was sent to the Genomics Resource Center, Institute for Genome Sciences at the University of Maryland School of Medicine (Baltimore, MD, USA), for library preparation and sequencing. In total, two SRMT cells were prepared and sequenced on a PacBio Sequel System.
RNA-Seq library preparation and sequencing
RNA-Seq libraries were made from four RNA extractions, two from total silk gland tissue and two from cephalothorax tissue. RNA samples were sent to Novogene (Sacramento, CA, USA) for library preparation and sequenced on an Illumina HiSeq System, paired-end 150 bp.
Transcriptome assembly and expression analyses
Raw sequencing reads were trimmed of low quality bases and adaptors using Trimmomatic . Quality of filtered reads was evaluated with FastQC. All reads from cephalothorax and total silk gland tissues were combined for a de novo assembly transcriptome using Trinity v2.8.5 with default parameters . We developed a novel Transcriptome Trimming and Annotation Pipeline (TrTAP), to select a single allele per Trinity defined gene and to remove low-confidence transcripts (Additional file 3: Figure S1). Briefly, TrTAP uses BLAST comparisons to high quality spider and arthropod genomes to identify transcripts with well-supported alignments to homologs and RSEM  to identify transcripts meeting a threshold of expression. Transcripts identified as chimeric using scripts from Clarke et al. 2014 or as rRNAs and tRNAs using either SILVA 18 S rRNA v123  or tRNAscan  are excluded. The expression of representative transcripts (one per Trinity-defined gene) was obtained using RSEM v1.3.1 and each transcript was then compared to taxonomically-diverse databases such as PFAM v34.0  and SwissProt v2018_03  for annotation and GO term assignment. The retained transcripts are those that: (1) are the reciprocal best BLAST hits (BEST) or are the best match to a database gene, even if another allele in the Trinity-defined gene is the reciprocal best BLAST hit (2BEST); (2) are not the reciprocal best hit, but have sufficient expression and a significant BLAST hit with sufficient coverage (GOOD); or (3) have no BLASTP hit, but an ORF with sufficient length and expression (LONGORF). The rest of the transcripts are discarded. We used cutoffs of BLAST coverage of 20% of the full database gene, amino acids length of 50, and expression of 1 TPM in at least one library. Scripts and documentation for TrTAP are available at http://www.github.com/thclarke/TrTAP. Transcriptome quality was approximated with N50 and completeness evaluated by comparison to the arthropod v9 set of Universal Single-Copy Orthologs (BUSCO v 3.0.2).
Cleaned RNA-Seq reads were mapped back to the TrTAP single-copy transcriptome using Bowtie2 v2.4.2 with default parameters . Differentially expressed genes were obtained from Bowtie counts of unambiguous matches using the R package DESeq2 with default parameters . Genes with a Benjamini-Hochberg adjusted P-value < 0.001 were considered differentially expressed . GOSlim term enrichment of the differentially expressed genes were calculated using the goseq module  with GO Slim terms obtained from GO Slim viewer  with the TrTAP assigned GO terms.
Silk gene expression was evaluated using GMAP  to map and align transcripts to the genome with default parameters, followed by STAR  which mapped RNA-seq reads to the genome, limiting the spidroin read counts to the terminal regions only. TPM values were calculated for all transcripts.
Gene annotation and phylogenetic analyses
In addition to the TrTAP based-annotations, the genomic and Iso-Seq transcriptomic assemblies were subject to BLASTX searches (e-value < e-10) against published spider genomes: Araneus ventricosus (BGPR00000000.1), Trichonephila clavipes (MWRG00000000.1), and Stegodyphus mimosarum (AZAQ00000000.1), in addition to a protein database of published silk proteins downloaded from UniProtKB/SwissProt. Contigs were visually inspected and error corrected manually using Sequencher v5.4.6 and Geneious . Spidroin contigs were assigned a spidroin type based on their repetitive sequence and phylogenetic relationships in the spidroin gene trees based on N- and C-terminal regions.
Selected spidroin contigs were translated, and the amino (N-) and carboxyl (C-) terminal regions were combined with spidroin sequences from published cribellate and ecribellate species (Supplementary Table S3). Alignments of the N- and C- terminal regions were done separately and then concatenated using MUSCLE implemented in Geneious (Additional file 1). Nucleotide sequences of selected AgSp, Flag, Pflag, Sp, and CrSp terminal region sequences (Additional file 3: Table S4) as well as an alignment of tandemly arrayed U. diversus Pflag repeat unit sequences were obtained separately using MUSCLE. Maximum likelihood analyses were constructed with 10,000 bootstrap replicates by RAxML v8 . GAMMA amino acid model was used for the N- and C- concatenated and the Pflag/Flag analyses. For the Pflag repeat nucleotide alignment, sequences were clustered with the Unweighted Pair Group Method with Arithmetic Mean (UPGMA) implemented in Geneious. FigTree v1.4.4 was used to visualize the resulting trees from all analyses.
All sequencing data used in this study are publicly available in the National Center for Biotechnology Information (NCBI) accession number PRJNA747999.
Blackledge TA, Hayashi CY. Silken toolkits: biomechanics of silk fibers spun by the orb web spider Argiope argentata (Fabricius 1775). J Exp Biol. 2006;209:2452–61. https://doi.org/10.1242/jeb.02275.
Sensenig A, Agnarsson I, Blackledge TA. Behavioural and biomaterial coevolution in spider orb webs. J Evol Biol. 2010;23:1839–56. https://doi.org/10.1111/j.1420-9101.2010.02048.x.
Sensenig AT, Kelly SP, Lorentz KA, Lesher B, Blackledge TA. Mechanical performance of spider orb webs is tuned for high-speed prey. J Exp Biol. 2013;216:3388–94. https://doi.org/10.1242/jeb.085571.
Peters HM. Fine structure and function of capture threads. Berlin: Spiders; 1987. p. 187–202. https://doi.org/10.1007/978-3-642-71552-5_13.
World Spider Catalog. World Spider Catalog. World Spider Cat. 2021. http://www.wsc.nmbe.ch. Accessed 7 Feb 2018.
Eberhard WG, Agnarsson I, Levi HW. Web forms and the phylogeny of theridiid spiders (Araneae: Theridiidae): chaos from order. Syst Biodivers. 2007;6:415–75.
Liu J, May-Collado LJ, Pekár S, Agnarsson I. A revised and dated phylogeny of cobweb spiders (Araneae, Araneoidea, Theridiidae): A predatory Cretaceous lineage diversifying in the era of the ants (Hymenoptera, Formicidae). Mol Phylogenet Evol 2015. https://doi.org/10.1016/j.ympev.2015.09.023.
Peters HM. The spinning apparatus of Uloboridae in relation to the structure and construction of capture threads (Arachnida, Araneida). Zoomorphology. 1984;104:96–104.
Peters HM. On the spinning apparatus and the structure of the capture threads of Deinopis subrufus (Araneae, Deinopidae). Zoomorphology. 1992;112:27–37. https://doi.org/10.1007/BF01632992.
Lubin YD. Web buiding and prey capture in the Uloboridae. Spiders Webs Behav Evol. 1986;8:132–71.
Opell B. Cribellar thread. Spider Ecophysiol. Wolfgang Nentwig. New York: Springer; 2013. p. 303–18.
Griswold CE, Coddington JA, Platnick NI, Forster RR. Towards a phylogeny of entelegyne spiders (Araneae, Araneomorphae, Entelegynae). J Arachnol. 1999;34:53–63.
Bond JE, Garrison NL, Hamilton CA, Godwin RL, Hedin M, Agnarsson I. Phylogenomics resolves a spider backbone phylogeny and rejects a prevailing paradigm for orb web evolution. Curr Biol. 2014;24:1765–71. https://doi.org/10.1016/j.cub.2014.06.034.
Kullmann EJ. The convergent development of orb-webs in cribellate and ecribellate spiders. Am Zool. 1972;12:395–405. https://doi.org/10.1093/icb/12.3.395.
Kaston BJ. The evolution of spider webs. Am Zool. 1964;4:191–207.
Fernández R, Kallal RJ, Dimitrov D, Ballesteros JA, Arnedo MA, Giribet G, et al. Phylogenomics, diversification dynamics, and comparative transcriptomics across the spider tree of life. Curr Biol. 2018;28:1489-1497.e5. https://doi.org/10.1016/j.cub.2018.03.064.
Wheeler WC, Coddington JA, Crowley LM, Dimitrov D, Goloboff PA, Griswold CE, et al. The spider tree of life: phylogeny of Araneae based on target-gene analyses from an extensive taxon sampling. Cladistics. 2017;33:574–616. https://doi.org/10.1111/cla.12182.
Dimitrov D, Lopardo L, Giribet G, Arnedo MA, Álvarez-Padilla F, Hormiga G. Tangled in a sparse spider web: single origin of orb weavers and their spinning work unravelled by denser taxonomic sampling. Proc R Soc B Biol Sci. 2012;279:1341–50. https://doi.org/10.1098/rspb.2011.2011.
Coddington JA. The monophyletic origin of the orb web. Spiders Webs Behav Evol Shear Ed Stanf Univ Press Stanf Calif. 1986;45:319–63.
Hormiga G, Griswold CE. Systematics, phylogeny, and evolution of orb-weaving spiders. Annu Rev Entomol. 2014;59:487–512. https://doi.org/10.1146/annurev-ento-011613-162046.
Opell B. Revision of the genera and tropical American species of the spider family Uloboridae Revisión de los géneros de las especies americanas tropicales de arañas de la familia Uloboridae. Bull Mus Comp Zool. 1979;148:443–549.
Forster RR, Platnick NI. A review of the archaeid spiders and their relatives, with notes on the limits of the superfamily Palpimanoidea (Arachnida, Araneae). Bull AMNH. 1984;178:1.
Opell BD. Prey handling and food extraction by the triangle-web spider Hyptiotes cavatus (Uloboridae). J Arachnol. 1988;16:272–4.
Eberhard WG, Barrantes G, Weng J-L. Tie them up tight: wrapping by Philoponella vicina spiders breaks, compresses and sometimes kills their prey. Naturwissenschaften. 2006;93:251–4. https://doi.org/10.1007/s00114-006-0094-1.
Weng J-L, Barrantes G, Eberhard WG. Feeding by Philoponella vicina (Araneae, Uloboridae) and how uloborid spiders lost their venom glands. Can J Zool. 2006;84:1752–62. https://doi.org/10.1139/z06-149.
Opell BD. Functional similarities of spider webs with diverse architectures. Am Nat. 1996;148:630–48.
Bertkau P. Über das cribellum und calamistrum: ein beitrag zur histiologie, biologie und systematik der spinnen. 1882.
Eberhard WG. Combing and sticky silk attachment behaviour by cribellate spiders and its taxonomic implications. Bull Br Arachnol Soc. 1988;67:45.
Joel A-C, Kappel P, Adamova H, Baumgartner W, Scholz I. Cribellate thread production in spiders: Complex processing of nano-fibres into a functional capture thread. Arthropod Struct Dev. 2015;44:568–73. https://doi.org/10.1016/j.asd.2015.07.003.
Eberhard W, Pereira F. Ultrastructure of cribellate silk of nine species in eight families and possible taxonomic implications (Araneae: Amaurobiidae, Deinopidae, Desidae, Dictynidae, Filistatidae, Hypochilidae, Stiphidiidae, Tengellidae). J Arachnol. 1993;21:161–74.
Peters HM, Kovoor J. A complement to the spinning apparatus in Uloboridae (Araneae): the paracribellum and its glands. Zoomorphologie. 1980;45:56.
Peters HM. Structure and Formation of the Tentacles of cribellate spiders (Arachnida: Araneae). Verhandlungen: Naturwissenschaftlichen Ver Hambg; 1983. p. 241–54.
Michalik P, Piorkowski D, Blackledge TA, Ramírez MJ. Functional trade-offs in cribellate silk mediated by spinning behavior. Sci Rep. 2019. https://doi.org/10.1038/s41598-019-45552-x.
Grannemann CCF, Meyer M, Reinhardt M, Ramírez MJ, Herberstein ME, Joel A-C. Small behavioral adaptations enable more effective prey capture by producing 3D-structured spider threads. Sci Rep. 2019;9:17273. https://doi.org/10.1038/s41598-019-53764-4.
Griswold CE, Ramírez MJ, Coddington JA, Platnick NI. Atlas of phylogenetic data for Entelegyne spiders (Araneae: Araneomorphae: Entelegynae), with comments on their phylogeny. Proc-Calif Acad Sci. 2005;56:1.
Opell BD. The ability of spider cribellar prey capture thread to hold insects with different surface features. Funct Ecol. 1994;8:145–50. https://doi.org/10.2307/2389897.
Hawthorn AC, Opell BD. Evolution of adhesive mechanisms in cribellar spider prey capture thread: evidence for van der waals and hygroscopic forces: evolution of adhesive mechanisms. Biol J Linn Soc. 2002;77:1–8. https://doi.org/10.1046/j.1095-8312.2002.00099.x.
Hawthorn AC, Opell BD. van der Waals and hygroscopic forces of adhesion generated by spider capture threads. J Exp Biol. 2003;206:3905–11. https://doi.org/10.1242/jeb.00618.
Bott RA, Baumgartner W, Bräunig P, Menzel F, Joel A-C. Adhesion enhancement of cribellate capture threads by epicuticular waxes of the insect prey sheds new light on spider web evolution. Proc R Soc B. 2017;284:20170363. https://doi.org/10.1098/rspb.2017.0363.
Blackledge TA, Hayashi CY. Unraveling the mechanical properties of composite silk threads spun by cribellate orb-weaving spiders. J Exp Biol. 2006;209:3131–40. https://doi.org/10.1242/jeb.02327.
Hinman MB, Lewis RV. Isolation of a clone encoding a second dragline silk fibroin. Nephila clavipes dragline silk is a two-protein fiber. J Biol Chem. 1992;267:19320–4.
Guerette PA, Ginzinger DG, Weber BHF, Gosline JM. Silk properties determined by gland-specific expression of a spider fibroin gene family. Science. 1996;272:112–5. https://doi.org/10.1126/science.272.5258.112.
Hayashi CY, Lewis RV. Evidence from flagelliform silk cDNA for the structural basis of elasticity and modular nature of spider silks. J Mol Biol. 1998;275:773–84. https://doi.org/10.1006/jmbi.1997.1478.
Ayoub NA, Hayashi CY. Multiple recombining loci encode MaSp1, the primary constituent of dragline silk, in widow spiders (Latrodectus: Theridiidae). Mol Biol Evol. 2008;25:277–86. https://doi.org/10.1093/molbev/msm246.
Blasingame E, Tuton-Blasingame T, Larkin L, Falick AM, Zhao L, Fong J, et al. Pyriform spidroin 1, a novel member of the silk gene family that anchors dragline silk fibers in attachment discs of the black widow spider Latrodectus hesperus. J Biol Chem. 2009;284:29097–108. https://doi.org/10.1074/jbc.M109.021378.
Lane AK, Hayashi CY, Whitworth GB, Ayoub NA. Complex gene expression in the dragline silk producing glands of the Western black widow (Latrodectus hesperus). BMC Genomics. 2013;14:846. https://doi.org/10.1186/1471-2164-14-846.
Collin MA, Clarke TH, Ayoub NA, Hayashi CY. Evidence from multiple species that spider silk glue component ASG2 is a spidroin. Sci Rep. 2016;6:21589. https://doi.org/10.1038/srep21589.
Clarke TH, Garb JE, Haney RA, Chaw RC, Hayashi CY, Ayoub NA. Evolutionary shifts in gene expression decoupled from gene duplication across functionally distinct spider silk glands. Sci Rep. 2017. https://doi.org/10.1038/s41598-017-07388-1.
Correa-Garhwal SM, Chaw RC, Iii THC, Alaniz LG, Chan FS, Alfaro RE, et al. Silk genes and silk gene expression in the spider Tengella perfuga (Zoropsidae), including a potential cribellar spidroin (CrSp). PLoS ONE. 2018;13: e0203563. https://doi.org/10.1371/journal.pone.0203563.
Kono N, Nakamura H, Mori M, Tomita M, Arakawa K. Spidroin profiling of cribellate spiders provides insight into the evolution of spider prey capture strategies. Sci Rep. 2020;10:15721. https://doi.org/10.1038/s41598-020-72888-6.
Gregory TR, Shorthouse DP. Genome Sizes of Spiders. J Hered. 2003;94:285–90. https://doi.org/10.1093/jhered/esg070.
Simão FA, Waterhouse RM, Ioannidis P, Kriventseva EV, Zdobnov EM. BUSCO: assessing genome assembly and annotation completeness with single-copy orthologs. Bioinformatics. 2015;31:3210–2. https://doi.org/10.1093/bioinformatics/btv351.
Collin MA, Clarke TH, Ayoub NA, Hayashi CY. Genomic perspectives of spider silk genes through target capture sequencing: Conservation of stabilization mechanisms and homology-based structural models of spidroin terminal regions. Int J Biol Macromol. 2018;113:829–40. https://doi.org/10.1016/j.ijbiomac.2018.02.032.
Ittah S, Michaeli A, Goldblum A, Gat U. A model for the structure of the C-terminal domain of dragline spider silk and the role of its conserved cysteine. Biomacromol. 2007;8:2768–73. https://doi.org/10.1021/bm7004559.
Gauthier M, Leclerc J, Lefèvre T, Gagné SM, Auger M. Effect of pH on the structure of the recombinant C-terminal domain of Nephila clavipes dragline silk protein. Biomacromol. 2014;15:4447–54. https://doi.org/10.1021/bm501241n.
Garb JE, Hayashi CY. Modular evolution of egg case silk genes across orb-weaving spider superfamilies. Proc Natl Acad Sci U S A. 2005;102:11379–84. https://doi.org/10.1073/pnas.0502473102.
Garb JE, Dimauro T, Vo V, Hayashi CY. Silk genes support the single origin of orb webs. Science. 2006;312:1762. https://doi.org/10.1126/science.1127946.
Sanggaard KW, Bechsgaard JS, Fang X, Duan J, Dyrlund TF, Gupta V, et al. Spider genomes provide insight into composition and evolution of venom and silk. Nat Commun. 2014;5:3765. https://doi.org/10.1038/ncomms4765.
Correa-Garhwal SM, Clarke TH, Janssen M, Crevecoeur L, McQuillan BN, Simpson AH, et al. Spidroins and Silk Fibers of Aquatic Spiders. Sci Rep. 2019;9:13656. https://doi.org/10.1038/s41598-019-49587-y.
Clarke TH, Garb JE, Hayashi CY, Haney RA, Lancaster AK, Corbett S, et al. Multi-tissue transcriptomics of the black widow spider reveals expansions, co-options, and functional processes of the silk gland gene toolkit. BMC Genomics. 2014. https://doi.org/10.1186/1471-2164-15-365.
Chaw RC, Clarke TH, Arensburger P, Ayoub NA, Hayashi CY. Gene expression profiling reveals candidate genes for defining spider silk gland types. Insect Biochem Mol Biol. 2021;135: 103594. https://doi.org/10.1016/j.ibmb.2021.103594.
Kallal RJ, Kulkarni SS, Dimitrov D, Benavides LR, Arnedo MA, Giribet G, et al. Converging on the orb: denser taxon sampling elucidates spider phylogeny and new analytical methods support repeated evolution of the orb web. Cladistics. 2020. https://doi.org/10.1111/cla.12439.
Hayashi CY, Lewis RV. Spider flagelliform silk: lessons in protein design, gene structure, and molecular evolution. BioEssays. 2001;23:750–6.
Babb PL, Lahens NF, Correa-Garhwal SM, Nicholson DN, Kim EJ, Hogenesch JB, et al. The Nephila clavipes genome highlights the diversity of spider silk genes and their complex expression. Nat Genet. 2017;49:895–903. https://doi.org/10.1038/ng.3852.
Kono N, Nakamura H, Ohtoshi R, Moran DAP, Shinohara A, Yoshida Y, et al. Orb-weaving spider Araneus ventricosus genome elucidates the spidroin gene catalogue. Sci Rep. 2019;9:1–13. https://doi.org/10.1038/s41598-019-44775-2.
Correa-Garhwal SM, Babb PL, Voight BF, Hayashi CY. Golden orb-weaving spider (Trichonephila clavipes) silk genes with sex-biased expression and atypical architectures. G3 Genes Genomes Genetics. 2021;11:1–10. https://doi.org/10.1093/g3journal/jkaa039.
Kallal RJ, Dimitrov D, Arnedo MA, Giribet G, Hormiga G. Monophyly, Taxon Sampling, and the Nature of Ranks in the Classification of Orb-Weaving Spiders (Araneae: Araneoidea). Syst Biol. 2020;69:401–11. https://doi.org/10.1093/sysbio/syz043.
Foelix R. Biology of Spiders. Oxford: Oxford University Press; 2011.
Sahni V, Blackledge TA, Dhinojwala A. A Review on Spider Silk Adhesion. J Adhes. 2011. https://doi.org/10.1080/00218464.2011.583588.
Vollrath F, Selden P. The Role of Behavior in the Evolution of Spiders, Silks, and Webs. Annu Rev Ecol Evol Syst. 2007;38:819–46. https://doi.org/10.1146/annurev.ecolsys.37.091305.110221.
Hayashi CY, Blackledge TA, Lewis RV. Molecular and mechanical characterization of aciniform silk: uniformity of iterated sequence modules in a novel member of the spider silk fibroin gene family. Mol Biol Evol. 2004;21:1950–9. https://doi.org/10.1093/molbev/msh204.
Tian M, Lewis RV. Molecular characterization and evolutionary study of spider tubuliform (eggcase) silk protein. Biochemistry. 2005;44:8006–12.
Ayoub NA, Garb JE, Kuelbs A, Hayashi CY. Ancient properties of spider silks revealed by the complete gene sequence of the prey-wrapping silk protein (AcSp1). Mol Biol Evol. 2013;30:589–601. https://doi.org/10.1093/molbev/mss254.
Chaw RC, Saski CA, Hayashi CY. Complete gene sequence of spider attachment silk protein (PySp1) reveals novel linker regions and extreme repeat homogenization. Insect Biochem Mol Biol. 2017;81:80–90. https://doi.org/10.1016/j.ibmb.2017.01.002.
Ayoub NA, Garb JE, Tinghitella RM, Collin MA, Hayashi CY. Blueprint for a high-performance biomaterial: full-length spider dragline silk genes. PLoS ONE. 2007;2: e514. https://doi.org/10.1371/journal.pone.0000514.
Opell BD. Factors governing the stickiness of cribellar prey capture threads in the spider family Uloboridae. J Morphol. 1994;221:111–9. https://doi.org/10.1002/jmor.1052210109.
Piorkowski D. Punctuated evolution of viscid silk in spider orb webs supported by mechanical behavior of wet cribellate silk 2017:12.
Piorkowski D, Blackledge TA, Liao C-P, Joel A-C, Weissbach M, Wu C-L, et al. Uncoiling springs promote mechanical functionality of spider cribellate silk. J Exp Biol. 2020. https://doi.org/10.1242/jeb.215269.
Piorkowski D, Liao C-P, Joel A-C, Wu C-L, Doran N, Blamires S, et al. Adhesion of spider cribellate silk enhanced in high humidity by mechanical plasticization of the underlying fiber. J Mech Behav Biomed Mater. 2021;114: 104200. https://doi.org/10.1016/j.jmbbm.2020.104200.
Hajer J. Notes on the spinning of the spiders Hyptiotes paradoxus CLK, 1834, and Uloborus wakkenaerius Latr, 1806 (Araneae: Uloboridae). Bulletin de la Société des Science Naturelles de Neuchâtel 1991;116:99–103. https://doi.org/10.5169/seals-89371.
Alfaro RE, Griswold CE, Miller KB. Comparative spigot ontogeny across the spider tree of life. PeerJ. 2018;6: e4233. https://doi.org/10.7717/peerj.4233.
Schwager EE, Sharma PP, Clarke T, Leite DJ, Wierschin T, Pechmann M, et al. The house spider genome reveals an ancient whole-genome duplication during arachnid evolution. BMC Biol. 2017;15:62. https://doi.org/10.1186/s12915-017-0399-x.
Blackledge TA, Summers AP, Hayashi CY. Gumfooted lines in black widow cobwebs and the mechanical properties of spider capture silk. Zool Jena Ger. 2005;108:41–6. https://doi.org/10.1016/j.zool.2004.11.001.
Ayoub NA, Friend K, Clarke T, Baker R, Correa-Garhwal S, Crean A, et al. Protein composition and associated material properties of cobweb spiders’ gumfoot glue droplets. Integr Comp Biol. 2021. https://doi.org/10.1093/icb/icab086.
Vollrath F, Fairbrother WJ, Williams RJP, Tillinghast EK, Bernstein DT, Gallagher KS, et al. Compounds in the droplets of the orb spider’s viscid spiral. Nature. 1990;345:526–8. https://doi.org/10.1038/345526a0.
Vollrath F, Tillinghast EK. Glycoprotein glue beneath a spider web’s aqueous coat. Naturwissenschaften. 1991;78:557–9. https://doi.org/10.1007/BF01134447.
Opell BD, Hendricks ML. Adhesive recruitment by the viscous capture threads of araneoid orb-weaving spiders. J Exp Biol. 2007;210:553–60. https://doi.org/10.1242/jeb.02682.
Opell BD, Jain D, Dhinojwala A, Blackledge TA. Tuning orb spider glycoprotein glue performance to habitat humidity. J Exp Biol. 2018. https://doi.org/10.1242/jeb.161539.
Opell BD. Functional associations between the cribellum spinning plate and capture threads of Miagrammopes animotus (Araneida, Uloboridae). Zoomorphology. 1989;108:263–7. https://doi.org/10.1007/BF00312159.
Vollrath F, Madsen B, Shao Z. The effect of spinning conditions on the mechanics of a spider’s dragline silk. Proc R Soc B Biol Sci. 2001;268:2339–46. https://doi.org/10.1098/rspb.2001.1590.
Blackledge TA, Boutry C, Wong S-C, Baji A, Dhinojwala A, Sahni V, et al. How super is supercontraction? Persistent versus cyclic responses to humidity in spider dragline silk. J Exp Biol. 2009;212:1981–9. https://doi.org/10.1242/jeb.028944.
Craig CL, Riekel C, Herberstein ME, Weber RS, Kaplan D, Pierce NE. Evidence for diet effects on the composition of silk proteins produced by spiders. Mol Biol Evol. 2000;17:1904–13.
Liao X, Yin G, Huang Z, Yao Y, Gu J, Han D. Supercontraction on cribellate spider spiral silk with wet-rebuilt micro-structure. Mater Sci Eng C. 2011;31:128–33. https://doi.org/10.1016/j.msec.2010.08.002.
Opell BD, Bond JE. Capture thread extensibility of orb-weaving spiders: testing punctuated and associative explanations of character evolution. Biol J Linn Soc. 2000;70:107–20. https://doi.org/10.1111/j.1095-8312.2000.tb00203.x.
Becker N, Oroudjev E, Mutz S, Cleveland JP, Hansma PK, Hayashi CY, et al. Molecular nanosprings in spider capture-silk threads. Nat Mater. 2003;2:278–83. https://doi.org/10.1038/nmat858.
Swanson BO, Blackledge TA, Hayashi CY. Spider capture silk: performance implications of variation in an exceptional biomaterial. J Exp Zool Part Ecol Genet Physiol. 2007;307:654–66. https://doi.org/10.1002/jez.420.
Piorkowski D, Blackledge TA. Punctuated evolution of viscid silk in spider orb webs supported by mechanical behavior of wet cribellate silk. Sci Nat. 2017;104:67. https://doi.org/10.1007/s00114-017-1489-x.
Gosline JM, Denny MW, DeMont ME. Spider silk as rubber. Nature. 1984;309:551–2. https://doi.org/10.1038/309551a0.
Stellwagen SD, Renberg RL. Toward spider glue: long read scaffolding for extreme length and repetitious silk family genes AgSp1 and AgSp2 with insights into functional adaptation. GGenesGenomesGenetics. 2019;9:1909–19. https://doi.org/10.1534/g3.119.400065.
Diaz C, Baker RH, Long JH, Hayashi CY. Connecting materials, performance and evolution: a case study of the glue of moth-catching spiders (Cyrtarachninae). J Exp Biol. 2022;225:243271. https://doi.org/10.1242/jeb.243271.
Bolger AM, Lohse M, Usadel B. Trimmomatic: a flexible trimmer for Illumina sequence data. Bioinformatics. 2014;30:2114–20. https://doi.org/10.1093/bioinformatics/btu170.
Haas BJ, Papanicolaou A, Yassour M, Grabherr M, Blood PD, Bowden J, et al. De novo transcript sequence reconstruction from RNA-seq using the Trinity platform for reference generation and analysis. Nat Protoc. 2013;8:1494–512. https://doi.org/10.1038/nprot.2013.084.
Li B, Dewey CN. RSEM: accurate transcript quantification from RNA-Seq data with or without a reference genome. BMC Bioinformatics. 2011;12:323. https://doi.org/10.1186/1471-2105-12-323.
Quast C, Pruesse E, Yilmaz P, Gerken J, Schweer T, Yarza P, et al. The SILVA ribosomal RNA gene database project: improved data processing and web-based tools. Nucleic Acids Res. 2013;41:D590–6. https://doi.org/10.1093/nar/gks1219.
Chan PP, Lowe TM. tRNAscan-SE: Searching for tRNA Genes in Genomic Sequences. Methods Mol Biol Clifton NJ. 2019;1962:1–14. https://doi.org/10.1007/978-1-4939-9173-0_1.
Mistry J, Chuguransky S, Williams L, Qureshi M, Salazar GA, Sonnhammer ELL, et al. Pfam: The protein families database in 2021. Nucleic Acids Res. 2021;49:D412–9. https://doi.org/10.1093/nar/gkaa913.
The UniProt Consortium. UniProt: the universal protein knowledgebase in 2021. Nucleic Acids Res. 2021;49:D480–9. https://doi.org/10.1093/nar/gkaa1100.
Langmead B, Salzberg SL. Fast gapped-read alignment with Bowtie 2. Nat Methods. 2012;9:357–9. https://doi.org/10.1038/nmeth.1923.
Love MI, Huber W, Anders S. Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Biol. 2014;15:550. https://doi.org/10.1186/s13059-014-0550-8.
Benjamini Y, Hochberg Y. Controlling the false discovery rate: a practical and powerful approach to multiple testing. J R Stat Soc Ser B Methodol. 1995;57:289–300.
Young MD, Wakefield MJ, Smyth GK, Oshlack A. Gene ontology analysis for RNA-seq: accounting for selection bias. Genome Biol. 2010;11:R14. https://doi.org/10.1186/gb-2010-11-2-r14.
McCarthy FM, Bridges SM, Wang N, Magee GB, Williams WP, Luthe DS, et al. AgBase: a unified resource for functional analysis in agriculture. Nucleic Acids Res. 2007;35:D599-603. https://doi.org/10.1093/nar/gkl936.
Wu TD, Reeder J, Lawrence M, Becker G, Brauer MJ. GMAP and GSNAP for genomic sequence alignment: enhancements to speed, accuracy, and functionality. Methods Mol Biol Clifton NJ. 2016;1418:283–334. https://doi.org/10.1007/978-1-4939-3578-9_15.
Dobin A, Davis CA, Schlesinger F, Drenkow J, Zaleski C, Jha S, et al. STAR: ultrafast universal RNA-seq aligner. Bioinforma Oxf Engl. 2013;29:15–21. https://doi.org/10.1093/bioinformatics/bts635.
Kearse M, Moir R, Wilson A, Stones-Havas S, Cheung M, Sturrock S, et al. Geneious Basic: An integrated and extendable desktop software platform for the organization and analysis of sequence data. Bioinformatics. 2012;28:1647–9. https://doi.org/10.1093/bioinformatics/bts199.
Stamatakis A. RAxML version 8: A tool for phylogenetic analysis and post-analysis of large phylogenies. Bioinformatics. 2014;9:1312. https://doi.org/10.1093/bioinformatics/btu033.
We thank Thomas Dugger for his assistance in collecting spiders. We thank Brent Opell, Gustavo Hormiga, Kyle Friend, and John Gatesy for improving the manuscript.
This research was supported by the National Science Foundation (Award Numbers: IOS-1755142 to NAA, IOS-1754979 to CYH and RB, and IntBIO-2128028 to SCG, CYH, and RB).
Ethics approval and consent to participate
Consent for publication
The authors declare no competing interest.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Protein sequence alignment of sequences used in phylogenetic analyses for Figure 1.
Counts and percentages of the alleles U. diversus transcriptome in each of the final TrTAP category both by total and by annotation.
Supplementary Figures S1–S6 and Supplementary Tables S2–S4.
About this article
Cite this article
Correa-Garhwal, S.M., Baker, R.H., Clarke, T.H. et al. The evolutionary history of cribellate orb-weaver capture thread spidroins. BMC Ecol Evo 22, 89 (2022). https://doi.org/10.1186/s12862-022-02042-5