Skip to main content

Are personalities genetically determined? Inferences from subsocial spiders



Recent research has revealed that polymorphic behavioral strategies shape intra-and interspecific interactions and contribute to fitness in many animal species. A better understanding of the proximate mechanisms underlying these behavioral syndromes will enhance our grasp this phenomenon. Spiders in the genus Anelosimus exhibit inter-individual behavioral variation on several axes: individuals have consistent responses to stimuli (e.g. bold vs. shy individuals) and they are subsocial (exhibiting extended maternal care and sibling cooperation) across most of their range, but they sometimes form permanent social groups in northern temperate regions. Here, we seek genetic variants associated with boldness and with social structure in a socially polymorphic population of the spider Anelosimus studiosus. We also develop preliminary genomic resources, including a genome assembly and linkage map, that support this and future genomic research on this group.


Remarkably, we identify a small genomic scaffold (~ 1200 bp) that harbors seven single nucleotide polymorphisms (SNPs) associated with boldness. Moreover, heterozygotes are less common than expected based on Hardy-Weinberg equilibrium, suggesting that either assortative mating or selection against heterozygotes may be occurring in this system. We find no loci significantly associated with social organization. Our draft genome assembly allows us to localize SNPs of interest in this study and to carry out genetic comparisons with other published genomes, although it remains highly fragmented.


By identifying a locus associated with a well-studied animal personality trait, this study opens up avenues for future research to link behavioral studies of animal personality with genotype and fitness.


Much of classical behavioral ecology has been led by optimality models, which tend to view inter-individual variation as statistical noise [17]. Nowadays, however, researchers recognize that many species harbor individual variants with highly consistent responses to different stimuli [50]. These inter-individual differences, sometimes referred to as ‘personalities’, modulate within- and between-species interactions, including competition, predator defense, and mating strategy [53]. Although repeatable personality differences have been documented in animals ranging from arthropods to mammals and are known to be heritable in some systems [19], we still know little about the specific molecular mechanisms underlying the maintenance of this variation in most systems [32]. While some behavioral traits are generated through epigenetic mechanisms, the observation that personality can be heritable suggests that genetic variation in these traits is present [19].

Recent advances in sequencing technology provide unprecedented opportunities to develop genomic resources for non-model systems and allow researchers to seek genetic associations with a broad suite of phenotypic traits, including behaviors [9]. The application of genome-wide association studies (GWAS) to animal personality traits has lagged behind other phenotypic traits, due in part to the general assumption that behaviors are either entirely phenotypically plastic at one extreme or highly polygenic at the other [51]. So far, much of our understanding of the heritability of behavioral traits comes from quantitative genetic studies of captive rodent populations (reviewed by [51]), and, in these and other systems, bold versus shy or aggressive versus docile strategies appear to have a substantial genetic component [19]. However, studies of rodents and other organisms in natural populations reveal that gene-by-environment interactions can profoundly influence behavioral phenotypes (reviewed by [34]), supporting the importance of studying behavioral traits in natural populations as well. These quantitative genetic approaches are valuable in revealing the role of genetic inheritance in shaping behavioral variation; however, we are now able to take these studies a step further and ask what genetic architecture underlies personality traits.

The spider genus Anelosimus harbors a broad range of social strategies, from highly “social” species that live in large groups throughout their life cycle, to “subsocial” species that exhibit extended maternal care and cooperation among subadult siblings, to solitary species that disperse early in development [3]. In recent years, this group has been the subject of much research on inter-individual behavioral variation, with studies carried out on both subsocial and highly social species (e.g. [39, 40]). Anelosimus studiosus is among the most broadly distributed species in the genus, with its range extending from the northeastern United States through the Andes to Patagonia [1]. Across most of its range, this species is subsocial, but some larger multifemale social groups form in northern populations [20, 26, 27]. Colony social structure and personality are independent traits, such that bold and shy individuals may live in either social or solitary colonies [38]. Pruitt and Goodnight [38] showed that the composition of spider personalities in social colonies influenced their survival and reproduction, and that the optimal composition depended upon their environment. Thus, while individual personality is not causally linked with sociality in this species, the ratio of types of individuals within colonies is a major determinant of a group’s success.

In the present study, we develop genomic resources to support the integration of research on social evolution and animal personality with genetics. We pursue four goals. We first assemble a preliminary genome (goal 1) and high-density linkage map (goal 2) for A. studiosus. Using double digest Restriction-site Associated DNA sequencing (ddRADseq), we then carry out a genome-wide association study to identify genomic regions associated with spider personality (goal 3) and with colony social structure (goal 4). Identifying the genetic basis of well-studied behaviors will open new avenues for research on the mechanisms underlying selection for alternative behavioral strategies in A. studiosus and perhaps the broader Anelosimus genus.


Genome assembly

Based on our draft assembly, the A. studiosus genome is 2.22 Gpb with a GC content of 27.95%. The genome assembly N50 is 4622 bp; the maximum scaffold size is 78,878 bp. Using the function in BBMap [14], we estimate that the total haploid genome size is approximately 1.58 Gbp. The difference in estimated genome size and assembled genome size could indicate that our genome assembly includes some duplicated variant regions; alternatively, the ancient genome duplication present in spiders [48] could affect the precision of kmer-based estimates of genome size.

Linkage map assembly

Through maternal linkage maps for two families, we were able to place 284 and 541 markers in linkage groups containing five or more loci for families P1 and P2, respectively (Additional file 1: Tables S2, S3). We excluded small linkage groups (fewer than five loci) from further analyses (P1: 399 additional markers were on 301 small linkage groups, most of which contained a single locus; P2: 175 additional markers were on 135 small linkage groups). We obtained 24 and 19 linkage groups for families P1 and P2, respectively. Based on karyotypes, Anelosimus species have a haploid chromosome number of 12; females have 24 chromosomes (sex chromosomes are X1X1X2X2), while males have 22 chromosomes (sex chromosomes are X1X2O [5];). Thus, our linkage groups are oversplit, which is a common outcome for RADseq-based linkage maps produced with fewer than 50 offspring (Pers. Obs.). In total, 16 scaffolds harboring 28 SNPs overlapped between the two linkage maps on informative linkage groups, allowing us to determine several syntenic regions of the linkage maps (Additional file 1: Figure S1). During the linkage map assembly process, we noticed that the offspring groups appeared to have unexpected genetic structure. By carrying out principal component analyses (PCA) on offspring and the mother from both families, we were able to infer that both mothers probably mated with two different males (Fig. 1). We confirmed this by investigating relatedness patterns across a combined dataset containing spiders from our linkage mapping and association studies (Additional file 1: Figure S2). The average relatedness between individuals in putative half-sib groups was approximately half the relatedness between putative full siblings and mother/offspring pairs, which supports this inference.

Fig. 1

Principal component analyses of genetic structure between mother and offspring from family P1 (a) and family P2 (b) reveal a pattern consistent with multiple paternity of offspring. In both cases, the mother (enlarged black dot) is positioned between two offspring clusters (open dots) that are separated along principal component axis 1. For a follow-up assessment of relatedness (Additional file 1: Figure S2), we consider individuals with PC axis 1 values less than the mother to be putative full sib group 1, and those with values greater than the mother to be full sib group 2

Association study

We identified seven tightly linked markers that are strongly associated with boldness (p < = 1.06 × 10− 5; Fig. 2). These are all found on a small scaffold (1239 bp) in the A. studiosus genome assembly, and explain approximately 11% of the variance in boldness. The presence of both homozygotes and of heterozygotes and confirmation that the depth of coverage did not differ systematically across each genotype suggest that these genetic variants are not the result of copy number variation in our sample. We also note that both the reference and the alternate allele for the associated SNPs were present in all four populations, despite the rarity of the alternate allele (Fig. 3). We detect a total of five distinct genotypes at the markers associated with boldness on this scaffold: the majority of individuals were homozygous for the reference allele at all seven markers, the six heterozygous individuals were heterozygous at all markers, and eight out of nine individuals homozygous for the alternative allele (Fig. 2) were consistent at all markers. We detected one individual from Melton Hill that was heterozygous at three markers, and homozygous for the reference allele for the remaining four, but we note that the depth of coverage was relatively low (less than 10x) for the heterozygous markers in this individual. We also detected one individual from Melton Lake that was homozygous for the alternative allele at one marker and homozygous for the reference allele at the remaining markers, with a depth of coverage greater than 100 for all loci in this scaffold. Taken together, these results suggest that there are two major haplotypes on this scaffold, but that recombination may occur, breaking apart the multi-SNP haplotypes occasionally.

Fig. 2

Comparison of boldness score by genotype (a) in individuals at the most highly associated locus of seven highly significant SNPs on scaffold 422,272 (enlarged red dot, B). The outlier loci on the boxplots are shown with light grey open circles, while all data points (including outliers) are shown as jittered, black filled circles. The QQ plot (b) compares the expected and observed p-values with the grey bar indicating 95% confidence intervals, revealing that the observed associations exhibit more significant p-values than expected by chance

Fig. 3

The four sampling localities in eastern Tennessee are shown with orange stars (a), and the total number of nests collected that were either solitary (blue) or social (green) are shown with pie charts. The distribution of latency to move after perturbation (a measure of boldness) is shown for individuals from each of the four populations (b: Melton Lake, c: Melton Hill, d: Boone Lake, and e: Chilhowee Dam). Individuals that moved quickly (low latency score) are bold, whereas individuals with a slow response (high latency score) are shy. There was no association between nest type (social or solitary) and boldness; populations differed in the frequency of the alternate allele at the locus that is most strongly associated with boldness

Alignment to the house spider Parasteatoda tepidariorum revealed a partial hit of this scaffold to unplaced genomic scaffold with an e-value of 0.001 (Table 1). We then used blastn to align the focal A. studiosus scaffold to the broader NCBI database. Here, the top hit was to the Argiope argentata pyriform spidroin 1 gene (which encodes a silk protein). Since the sequence partially hit a coding region of a spider gene, we used blastx to examine hits to proteins (Table 1). Given the discovery that the top hit matches a transposable element in the ant Lasius niger, we suggest that the genomic region associated with boldness may be in a transposable element that is present in different fragments in the genomes of related spiders and other arthropods. We did not find any RADtags associated with the social structure of the colonies.

Table 1 NCBI blast query results based on search


Research on social spiders has given rise to broad insights about social evolution, sex ratio theory, and animal personality [2, 4, 6,7,8, 22,23,24, 39], but they have lagged behind other organisms in the development and application of genomic resources (but see [37, 49]). We identify multiple variants in a small genomic region that are strongly associated with boldness in the spider Anelosimus studiosus (Fig. 2), and we develop a preliminary genome assembly and linkage map that will facilitate future genomic studies of this genus (Additional file 1: Figure S1, Tables S2, S3). Previous studies demonstrated that highly repeatable inter-individual behavioral variants, including boldness, were strongly heritable in A. studiosus [38, 41]. This study takes our understanding of this trait a step further by identifying a putative genetic basis.

We discovered that SNPs associated with boldness all occurred on a single scaffold in our dataset. This scaffold also contained four SNPs that were not associated with boldness. Using one representative SNP with the most significant association to boldness, we showed that homozygotes for the alternative allele are the most likely to have a long duration of the death feign response, indicative of a shy individual, heterozygotes are intermediate, and homozygotes for the reference allele have the lowest score (indicative of bold individuals), although there is variation in behavioral response in all categories. We note that even traits that are strongly affected by genotype may also be influenced by the environment (e.g. via gene-by-environment interactions mediated by epigenetic mechanisms) or by other loci of minor effect. The alternative allele was present in all four sampled populations. Interestingly, we observed a deviation from Hardy-Weinberg equilibrium in this marker, with fewer heterozygous individuals than expected given the allele frequency (Chi-squared = 16.58, p = 0.00025). This pattern could be explored by seeking fitness differences between homozygous and heterozygous individuals at this locus. We speculate that having an intermediate boldness level may be disadvantageous in this system. Alternatively, assortative mating based upon behavioral syndrome could also account for the relative scarcity of heterozygotes.

The scaffold harboring the boldness-associated loci had marginal alignments to the house spider Parasteatoda tepidariorum genome and to Argiope argentata pyriform spidroin 1 gene, which encodes a silk protein (Table 1). Given its position in an A. argentata coding region, we translated our DNA sequence and reran a protein search. This yielded a well-supported hit to a transposase from the ant Lasius niger (as well as a secondary hit to a hypothetical protein in the spider Stegodyphus mimosarum). Given this finding, we speculate that our scaffold may overlap with a transposable element in the A. studiosus genome. Transposons often take on a regulatory role in the genome [16], however, our data do not yet allow us to explore downstream effects on gene expression. We for now predict that variation at this genomic region will be associated with different gene regulatory networks.

We uncovered no genetic markers associated with colony social structure in this sample. Given our reduced representation sequencing approach, this does not rule out a genetic basis for social group formation in this species, since we sampled only a small subset of the genome.

Another outcome of this project is the discovery that Anelosimus genomes have challenging features, including large size and low GC content. The genome size and composition are in line with other published spider genomes [46, 48]. The repetitive nature makes genome assembly challenging, although recent innovations in sequencing and long-range genome scaffolding approaches will facilitate substantial improvements. Despite the highly fragmented nature of the genome assembly presented here, which was based entirely on short insert Illumina sequencing, being able to align RADtags to a genome for the linkage mapping and association studies improved our ability to compare maps (Additional file 1: Figure S1) and to determine where SNPs of interest were found in other species’ genomes (Table 1). In general, the outcome of population genetic analyses using reduced representation genomic sequencing methods is enhanced by alignment to genome assemblies, even low quality ones [47]. The linkage maps will, in turn, be useful in scaffolding a more contiguous genome sequence in the future. The linkage mapping project also revealed the unexpected pattern of multiple paternity in two families (Fig. 1), which contradicts a previous study in this species that detected first-male sperm precedence [28]. Their study investigated the outcome of captive matings when different males were introduced over the course of several days. Although our sample is small, we sequenced field-collected brood, and differences in the timing and outcome of mating in natural conditions may account for the conflicting results.


Here, we take initial steps in developing genomic resources for a broadly useful behavioral model for the study of animal personality and social evolution. We propose several novel studies that could follow up on our preliminary results. The discovery of a genomic region that is strongly associated with boldness opens a variety of new research directions at the intersection of animal behavior and population genomics. There remains much to be discovered.


Genome and linkage map

In May 2014, we sampled an A. studiosus female subadult from a nest containing multiple females at Melton Hill Dam, Tennessee (35°89′N, 84°30′W) for whole genome sequencing. In June 2014, we sampled two families including a mother and her offspring (37 and 43 spiderlings, respectively) from solitary nests, also from Melton Hill Dam, in order to construct high-density genetic linkage maps.

Association study

In May 2014, we sampled 226 A. studiosus individuals from single- and multi-female colonies at four populations in Tennessee: Chilhowee Dam, Melton Hill, Melton Lake and Boone Lake (Fig. 3). Anelosimus studiosus spiders generally nest in vegetation along the water’s edge, and they tend to have a patchy distribution [29, 45]. Thus, we sampled nests opportunistically by boat and/or on foot at each site. These populations harbor spiders using two alternative social strategies. Some individuals form multi-female social groups, while others exhibit the more classical subsocial strategy of dispersing and founding independent nests before maturity [45]. These spiders usually live for approximately 1 year, with mature individuals producing offspring in the late spring or early summer. We targeted a mix of solitary and social nests, where we sampled the solitary female or a group of 2–5 individuals from a multi-female colony, respectively. We recorded the locality of each nest using a portable GPS. At each population, we took care to sample from multiple clusters of colonies separated by several meters or more.


We assessed the social structure of each colony in the field by counting the number of adult and subadult females present in the nest. Females were deemed to exhibit the subsocial strategy if only a singleton reproductive female resided in the nest, whereas females were deemed to exhibit the social strategy if two or more reproductive females cohabitated within the same next. To evaluate where females fell along the bold vs. shy axis, we transported nests back to laboratory at the University of Tennessee, Knoxville for startle tests, which are a common assay for evaluating boldness in spiders [25, 31, 44, 52]. In brief, spiders were taken from their nests and placed within a square open field and given 30 s to acclimate. Following this acclimation period two rapid jets of air were administered to the anterior prosoma of the spider using an infant nose-cleaning bulb. This stimulus elicits a death feign response in A. studiosus where females draw their legs tightly against their body and remain in a huddled position for a variable amount of time. We recorded the latency for each female to emerge from the huddle and move one body length. Individuals that hold a huddle for a longer period are deemed to be shy, whereas individuals that emerge from huddles swiftly are deemed to be bold. Following the phenotypic assessment, we placed spiders in 2 mL tubes containing 95% ethanol and transported them to the molecular laboratory. Open fields were cleaned with 95% ethanol between behavioral assays.

Sequencing and bioinformatics

For whole genome sequencing, we extracted DNA from the focal individual using a modified Insect Tissue protocol (based on Qiagen’s DNeasy protocol). In particular, we used alternative spin columns ( and eluted in 30 uL of AE buffer. We then sent the DNA to the University of Lausanne Genomics Core facility for library preparation with the TruSeq DNA kit and insert size of 500 bp and carried out 100 bp paired end sequencing on two Illumina 2500 lanes. We obtained 164 M unique reads or 33 Gbp, which amounted to approximately 14x coverage (the latter estimated through the function in BBMap, [14]).

We assembled the genome using the following pipeline: we used a custom shell command to remove reads that failed the Illumina chastity filter, we removed PCR duplicates with FastUniq [55], we combined overlapping paired end reads and removed adapters with PEAR [56], and we trimmed low quality bases with Sickle [30]. With these processed sequences, we assembled reads using SOAPdenovo2 [36] with a range of k values from 29 to 61, as well as with the multi k function. Small k values tended to yield larger but more fragmented assemblies, and we selected k43 as the best compromise between contiguity and completeness. We then scaffolded this assembly with SSPACE [10] and filled gaps with GapFiller [11].

For linkage map and association study processing, we extracted DNA following the protocol described above. We then used a customized ddRADseq protocol to amplify restriction fragments distributed across the genome [12, 42]. Briefly, we digested the genomic DNA with SbfI and MseI enzymes, ligated unique barcoded adapters to the SbfI cut site, removed small DNA fragments from each individual using AMPure beads (Agencourt), carried out PCR in three parallel reactions, size-selected fragments in the 400–500 bp range through a gel extraction, and purified the resulting DNA using a gel extraction kit (Qiagen). Finally, we sequenced each library on two Illumina 2500 lanes at the University of Lausanne genomics core facility. For the linkage map library, we used 100 bp single end reads, while we used 100 bp paired end sequencing for the association study.

For each library, we demultiplexed the sequence data using the process_radtags function in Stacks [15], aligned RADtags to our draft genome using Bowtie2 [33], and called variants using Samtools [35]. We filtered the resulting markers in VCFtools [18]. For the linkage maps, we initially separated the two families and carried out filtering independently; we retained genotype calls with a quality score greater than 20 (‘mingq = 20’), we retained markers present in at least 80% of individuals (geno = 80), and we removed markers with a minor allele frequency less than 15% (maf = 0.15). We examined missing data by individual, and we ultimately retained 32 and 36 offspring that had high quality genotype calls for at least 50% of loci; with these smaller families, we redid the geno and maf filtering.

For the association study, we retained markers that had a minimum genotype quality score of 20 or more (mingq = 20), were missing in less than 80% (geno = 0.8) of individuals, and that exhibited a minor allele frequency over 5% (maf = 0.05). We then examined the level of missing data in individuals and we removed 16 samples for which 25% or more loci were missing. Moving forward with 210 individuals, we reran the filters described above and generated a list of 4179 high-quality markers that were present in the majority of samples (mean read depth per site = 114x). In total, the sociality association test included 210 individuals at 4179 loci; four additional individuals were excluded from the boldness assay due to early mortality, resulting in a comparison of 206 individuals at 4179 loci.

Linkage map assembly

After filtering the data, we determined which markers were heterozygous in each maternal sample. Using VCFtools, we retained the list of maternally informative markers for each group of siblings for subsequent processing, in order to develop a maternal map. We then surveyed the sibling groups for markers that were variable among siblings, but homozygous in the mother, and retained these markers in order to develop a paternal map for each family. In order to construct the linkage maps, we exported the four datasets as 012 text files for manual conversion to MSTmap input format [54]. To account for the unknown phase of maternal and paternal genotypes, we duplicated each marker in the dataset, calling each genotype ‘A’ in one copy and ‘B’ in the other, following the procedure of Gadau [21]. We tested different parameter values in MSTmap. In particular, we tried cut_off_p_values ranging from 0.0005 to 0.000005. In general, larger p-value cut-offs resulted in too many markers being lumped together in single linkage groups, while small cut-offs resulted in greater fragmentation. For each map, we selected an intermediate p-value cut-off that provided a compromise between these extremes. We show the MSTmap parameter values in Additional file 1: Table S1. Before inferring male genotypes to construct paternal linkage maps, we inspected relationships among sib-groups. Based on principal component analyses (Fig. 1—1436 loci for P1 family, 7930 loci for P2 family) and assessment of relatedness (Additional file 1: Figure S2), we suspected that each female had mated with multiple males, though this contradicts a prior study that detected exclusive first-male sperm precedence in this species [28]. As a result, we could not construct paternal linkage maps for either family.

Association study analysis

We imputed individual genotypes across all our initial ddRADseq loci using BEAGLE [13]. We then used the list of high-quality markers determined through the filtering process described in the previous section to retain only the best markers from our dataset, now with missing data added through imputation. We used PLINK [43] to generate files in the bim and fam formats needed to assess associations and correct for underlying population genetic structure in GEMMA [57]. For each trait (boldness and sociality), we added the phenotype scores manually to the fam file. We ran GEMMA to calculate the association between each response variable and all the polymorphic markers, including underlying population genetic structure as a random effect in the model. To do this, we first ran GEMMA on the imputed VCF files with the -gk flag to calculate the relatedness matrix (this was done separately for each trait due to missing data for boldness) and then with the -lmm flag to perform the mixed effects model. We identified statistical outliers using the Bonferroni correction to obtain a conservative alpha value (α = 3.38 × 10− 5); when statistical outliers were present, we compared the observed p-value distribution to the expected distribution by generating a quantile-quantile plot in R. We assessed the proportion of variance in boldness that was explained by the most significantly associated locus (Fig. 2) by comparing the sum of squares for genotype to the sum of squares total.

Availability of data and materials

The genetic data generated during this study are available in the NCBI GenBank repository BioProject PRJNA557053, which includes genome raw reads, association study RADseq reads, and linkage map RADseq reads. This Whole Genome Shotgun project has been deposited at DDBJ/ENA/GenBank under the accession VSFD00000000. The version described in this paper is version VSFD01000000. The phenotypic data are provided in the supplementary information file associated with this article.


  1. 1.

    Agnarsson I. A revision of the New World eximius lineage of Anelosimus (Araneae, Theridiidae) and a phylogenetic analysis using worldwide exemplars. Zool J Linnean Soc. 2006;141:453–593.

    Article  Google Scholar 

  2. 2.

    Agnarsson I, Avilés L, Coddington JA, Maddison WP. Social theridiid spiders: repeated origins of an evolutionary dead-end. Evolution. 2006;60:2342–51.

    PubMed  Article  PubMed Central  Google Scholar 

  3. 3.

    Avilés L. Causes and consequences of cooperation and permanent sociality in spiders. In: Choe JC, Crespi BJ, editors. editors The Evolution of Social Behavior in Insects and Arachnids. Cambridge: Cambridge University Press; 1997. p. 476–98.

  4. 4.

    Avilés L, Agnarsson I, Salazar P, Purcell J, Iturralde G, Yip E, et al. Altitudinal patterns of spider sociality and the biology of a new mid-elevation social Anelosimus species in Ecuador. Am Nat. 2007;170:783–92.

    PubMed  Article  PubMed Central  Google Scholar 

  5. 5.

    Avilés L, Madison W. When is the sex ratio biased in social spiders?: chromosome studies of embryos and male meiosis in Anelosimus species (Araneae: Theridiidae). J Arachnol. 1991;19:126–35.

    Google Scholar 

  6. 6.

    Avilés L, McCormack J, Cutter A, Bukowski T. Precise, highly female-biased sex ratios in a social spider. Proc R Soc B. 2000;267:1445–9.

    PubMed  Article  PubMed Central  Google Scholar 

  7. 7.

    Avilés L, Purcell J. Anelosimus oritoyacu, a cloud forest social spider with only slightly female-biased primary sex ratios. J Arachnol. 2011;39:178–82.

    Article  Google Scholar 

  8. 8.

    Avilés L, Purcell J. The evolution of inbred social systems in spiders and other organisms: from short-term gains to long-term evolutionary dead-ends? Adv Study Behav. 2012;44:99–133.

    Article  Google Scholar 

  9. 9.

    Bengston SE, Dahan RA, Donaldson Z, Phelps SM, van Oers K, Sih A, et al. Genomic tools for behavioural ecologists to understand repeatable individual differences in behaviour. Nature Ecol Evol. 2018;2:944–55.

    Article  Google Scholar 

  10. 10.

    Boetzer M, Henkel CV, Jansen HJ, Butler D, Pirovano W. Scaffolding pre-assembled contigs using SSPACE. Bioinformatics. 2011;27:578–9.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  11. 11.

    Boetzer M, Pirovano W. Toward almost closed genomes with GapFiller. Genome Biol. 2012;13:R56.

    PubMed  PubMed Central  Article  Google Scholar 

  12. 12.

    Brelsford A, Dufresnes C, Perrin N. High-density sex-specific linkage maps of a European tree frog (Hyla arborea) identify the sex chromosome without information on offspring sex. Heredity. 2016;116:177–81.

    CAS  PubMed  Article  Google Scholar 

  13. 13.

    Browning SR, Browning BL. Rapid and accurate haplotype phasing and missing data inference for whole genome association studies by use of localized haplotype clustering. Am J Hum Genet. 2007;81:1084–97.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  14. 14.

    Bushnell B. BBMap: a fast, Accurate, Splice-Aware Aligner. 2014.

    Google Scholar 

  15. 15.

    Catchen J, Amores A, Hohenlohe P, Cresko W, Postlethwait J. Stacks: building and genotyping loci de novo from short-read sequences. G3: genes, genomes. Genetics. 2011;1:171–82.

    CAS  Google Scholar 

  16. 16.

    Chuong EB, Elde NC, Feschotte C. Regulatory activities of transposable elements: from conflicts to benefits. Nature Rev Genet. 2017;18:71–86.

    CAS  PubMed  Article  Google Scholar 

  17. 17.

    Dall SRX, Houston AI, McNamara JM. The behavioural ecology of personality: consistent individual differences from an adaptive perspective. Ecol Lett. 2004;7:734–9.

    Article  Google Scholar 

  18. 18.

    Danecek P, Auton A, Abecasis G, Albers CA, Banks E, DePristo MA, et al. The variant call format and VCFtools. Bioinformatics. 2011;27:2156–8.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  19. 19.

    Dochtermann NA, Schwab T, Sih A. The contribution of additive genetic variation to personality variation: heritability of personality. Proc R Soc B. 2015;282:20142201.

    PubMed  Article  Google Scholar 

  20. 20.

    Furey RE. Two cooperatively social populations of the theridiid spider Anelosimus studiosus in a temperate region. Anim Behav. 1998;55:727–35.

    CAS  PubMed  Article  Google Scholar 

  21. 21.

    Gadau J. Phase-unknown linkage mapping in ants. Cold Spring Harb Protoc. 2009;2009:5251.

    Google Scholar 

  22. 22.

    Grinsted L, Pruitt JN, Settepani V, Bilde T. Individual personalities shape task differentiation in a social spider. Proc R Soc B. 2013;280:1407.

    Article  Google Scholar 

  23. 23.

    Henschel JR, Lubin YD, Schneider J. Sexual competition in an inbreeding social spider, Stegodyphus dumicola (Araneae: Eresidae). Insect Soc. 1995;42:419–26.

    Article  Google Scholar 

  24. 24.

    Johannesen J, Lubin Y, Smith DR, Bilde T, Schneider JM. The age and evolution of sociality in Stegodyphys spiders: a molecular phylogenetic perspective. Proc R Soc B. 2006;274:231–7.

    Article  CAS  Google Scholar 

  25. 25.

    Jones TC, Akoury TS, Hauser CK, Neblett MF, Linville BJ, Edge AA, et al. Octopamine and serotonin have opposite effects on antipredator behavior in the orb-weaving spider, Larinioides cornutus. J Comp Phys A. 2011;197:819–25.

    CAS  Article  Google Scholar 

  26. 26.

    Jones TC, Parker PG. Costs and benefits of foraging associated with delayed juvenile dispersal in the spider Anelosimus studiosus (Araneae: Theridiidae). J Arachnol. 2000;28:61–9.

    CAS  Article  Google Scholar 

  27. 27.

    Jones TC, Parker PG. Delayed juvenile dispersal benefits both mother and offspring in the cooperative spider Anelosimus studiosus (Araneae: Theridiidae). Behav Ecol. 2002;13:142–8.

    Article  Google Scholar 

  28. 28.

    Jones TC, Parker PG. First male sperm precedence in multiply-mated females of the cooperative spider Anelosimus studiosus (Araneae, Theridiidae). J Arachnol. 2008;36:527–32.

    Article  Google Scholar 

  29. 29.

    Jones TC, Riechert SE. Patterns of reproductive success associated with social structure and microclimate in a spider system. Anim Behav. 2008;76:2011–9.

    Article  Google Scholar 

  30. 30.

    Joshi NA, Fass JN. Sickle: A sliding-window, adaptive, quality-based trimming tool for FastQ files. 2011.

    Google Scholar 

  31. 31.

    Kralj-Fiser S, Schneider JM. Individual behavioural consistency and plasticity in an urban spider. Anim Behav. 2012;84:197–204.

    Article  Google Scholar 

  32. 32.

    Laine VN, van Oers K. The quantitative and molecular genetics of individual differences in animal personality. In: Vonk J, Weiss A, Kuczaj S, editors. editors Personality in nonhuman animals. Springer: Cham; 2017. p. 55–72.

    Google Scholar 

  33. 33.

    Langmead B, Salzberg SL. Fast gapped-read alignment with bowtie 2. Nat Methods. 2012;9:357–9.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  34. 34.

    Ledón-Rettig CC, Richards CL, Martin LB. Epigenetics for behavioral ecologists. Behav Ecol. 2013;24:311–24.

    Article  Google Scholar 

  35. 35.

    Li H, Handsaker B, Wysoker A, Fennell T, Ruan J, Homer N, et al. The sequence alignment/map format and SAMtools. Bioinformatics. 2009;25:2078–9.

    PubMed  PubMed Central  Article  CAS  Google Scholar 

  36. 36.

    Luo R, Liu B, Xie Y, Li Z, Huang W, Yuan J, et al. SOAPdenovo2: an empirically improved memory-efficient short-read de novo assembler. GigaScience. 2012;1:18.

    PubMed  PubMed Central  Article  Google Scholar 

  37. 37.

    Mattila TM, Bechsgaard JS, Hansen TT, Schierup MH, Bilde T. Orthologous genes identified by transcriptome sequencing in the spider genus Stegodyphus. BMC Genomics. 2012;13:70.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  38. 38.

    Pruitt JN, Goodnight CJ. Site-specific group selection drives locally adapted group compositions. Nature. 2014;514:359.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  39. 39.

    Pruitt JN, Iturralde G, Avilés L, Riechert SE. Amazonian social spiders share similar within-colony behavioural variation and behavioural syndromes. Anim Behav. 2011;82:1449–55.

    Article  Google Scholar 

  40. 40.

    Pruitt JN, Oufiero CE, Avilés L, Riechert SE. Iterative evolution of increased behavioral variation characterizes the transition to sociality in spiders and proves advantageous. Am Nat. 2012;180:496–510.

    PubMed  Article  Google Scholar 

  41. 41.

    Pruitt JN, Riechert SE, Jones TC. Behavioural syndromes and their fitness consequences in a socially polymorphic spider, Anelosimus studiosus. Anim Behav. 2008;76:871–9.

    Article  Google Scholar 

  42. 42.

    Purcell J, Brelsford A, Wurm Y, Perrin N, Chapuisat M. Convergent genetic architecture underlies social organization in ants. Curr Biol. 2014;24:2728–32.

    CAS  PubMed  Article  Google Scholar 

  43. 43.

    Purcell S, Neale B, Todd-Brown K, Thomas L, Ferreira MAR, Bender D, et al. PLINK: a toolset for whole-genome association and population-based linkage analysis. Am J Hum Genet. 2007;81:559–75.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  44. 44.

    Riechert SE, Hedrick AV. A test for correlations among fitness-linked behavioural traits in the spider Agelenopsis aperta (Araneae, Agelenidae). Anim Behav. 1993;46:669–75.

    Article  Google Scholar 

  45. 45.

    Riechert SE, Jones TC. Phenotypic variation in the social behaviour of the spider Anelosimus studiosus along a latitudinal gradient. Anim Behav. 2008;75:1893–902.

    Article  Google Scholar 

  46. 46.

    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.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  47. 47.

    Schafer ABA, Peart CR, Tusso S, Maayan I, Brelsford A, Wheat CW, Wolf JBW. Bioinformatic processing of RAD-seq data dramatically impacts downstream population genetic inference. Methods Ecol Evol. 2017;8:907–17.

    Article  Google Scholar 

  48. 48.

    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.

    PubMed  PubMed Central  Article  CAS  Google Scholar 

  49. 49.

    Settepani V, Schou MF, Greve M, Grinsted L, Bechsgaard J, Bilde T. Evolution of sociality in spiders leads to depleted genomic diversity at both population and species levels. Mol Ecol. 2017;26:4197–210.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  50. 50.

    Sih A, Mathot KJ, Moirón M, Montiglio P-O, Wolf M, Dingemanse NJ. Animal personality and state-behaviour feedbacks: a review and guide for empiricists. TREE. 2015;30:50–60.

    PubMed  Google Scholar 

  51. 51.

    van Oers K, Sinn DL. The quantitative and molecular genetics of animal personality. In: Carere C, Maestripieri D, editors. editors Animal Personalities: Behavior, Physiology, and Evolution. Chicago: University of Chicago Press; 2013. p. 148–200.

  52. 52.

    Watts JC, Herrig A, Allen WD, Jones TC. Diel patterns of foraging aggression and antipredator behaviour in the trashline orb-weaving spider Cyclosa turbinate. Anim Behav. 2014;94:79–86.

    Article  Google Scholar 

  53. 53.

    Wolf M, Weissing FJ. Animal personalities: consequences for ecology and evolution. TREE. 2012;27:452–61.

    PubMed  Google Scholar 

  54. 54.

    Wu Y, Bhat PR, Close TJ, Lonardi S. Efficient and accurate construction of genetic linkage maps from the minimum spanning tree of a graph. PLoS Genet. 2008;4:e1000212.

    PubMed  PubMed Central  Article  CAS  Google Scholar 

  55. 55.

    Xu H, Luo X, Qian J, Pang X, Song J, Qian G, et al. FastUniq: a fast de novo duplicates removal tool for paired short reads. PLoS One. 2012;7:e52249.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  56. 56.

    Zhang J, Kobert K, Flouri T, Stamatakis A. PEAR: a fast and accurate Illumina paired-end reAd mergeR. Bioinformatics. 2014;30:614–20.

    CAS  Article  Google Scholar 

  57. 57.

    Zhou X, Stephens M. Genome-wide efficient mixed-model analysis for association studies. Nature Genet. 2012;44:821–4.

    CAS  PubMed  Article  Google Scholar 

Download references


The authors thank S.E. Riechert for her hospitality and F. Armagost, A.E. DeMarco, and T.C. Jones for their assistance during the field portion of this study. Thanks also to A. Brelsford, M. Sankovitz, A. de la Mora, and two anonymous reviewers for their feedback on the manuscript.


The authors thank the Fondation pour l’Université de Lausanne (450eAnniversaire, grant to JP) and the NSF (IOS # 1352705, to JNP) for funding this research. These agencies did not contribute to the collection, analysis, interpretation, or writing of this manuscript.

Author information




Both authors designed the study and carried out field work. JNP oversaw the phenotypic portion of the study, while JP performed the genetic data collection and analysis. JP wrote the first draft of the manuscript; JNP provided editorial feedback on an early draft. Both authors read and approved the final manuscript.

Corresponding author

Correspondence to Jessica Purcell.

Ethics declarations

Ethics approval and consent to participate

Not applicable.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary information

Additional file 1.

Detailed information about linkage maps is shown, including a comparison between linkage maps for P1 family and P2 family (Figure S1), methods and results of linkage map family relatedness calculations (Figure S2), MSTmap parameters (Table S1), and full linkage map alignments and SNP positions for P1 family (Table S2) and P2 family (Table S3). In addition, the full list of individuals used in the association study along with their social and boldness phenotypes is provided (Table S4).

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (, which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Purcell, J., Pruitt, J.N. Are personalities genetically determined? Inferences from subsocial spiders. BMC Genomics 20, 867 (2019).

Download citation


  • Social spiders
  • Boldness
  • Animal personality
  • RADseq
  • Genome assembly
  • Linkage map