Deciphering a potentially hyperdiverse diet of wandering spider,
(Phoneutria boliviensis; Araneae: Ctenidae) by DNA metabarcoding
of gut contents
Diego Sierra Ramírez1, Giovany
Guevara2, Lida Marcela Franco
Pérez3, Arie van der Meijden1,4,
Julio César González-Gómez1,3, Juan Carlos
Valenzuela-Rojas1, Carlos Fernando Prada
Quiroga1,*
1 Grupo de Investigación Biología y Ecología de
Artrópodos (BEA), Corporación Huiltur, Neiva, Huila y Facultad de
Ciencias, Universidad del Tolima, Altos de Santa Helena, Ibagué 730001,
Colombia
2 Grupo de Investigación en Zoología (GIZ), Facultad
de Ciencias, Universidad del Tolima, Altos de Santa Helena, Ibagué
730001, Colombia
3 Facultad de Ciencias Naturales y Matemáticas,
Universidad de Ibagué, Carrera 22 calle 67, Ibagué 730001, Colombia
4CIBIO Research Centre in Biodiversity and Genetic
Resources, InBIO, Universidade do Porto, Campus Agrário de Vairão, Rua
Padre Armando Quintas 7, 4485-661 Vairão, Vila do Conde, Portugal
Correspondence: Carlos F Prada Universidad del Tolima, Altos de Santa
Helena, Ibagué 730001, Colombia. Email:cfpradaq@ut.edu.co
Abstract
Arachnids are the most abundant land predators. Despite the importance
of their functional roles as predators and the of necessity to
understand their diet for conservation and nutrient fluxes, the trophic
ecology of many arachnid species is not fully understood. In the case of
the wandering spider, Phoneutria boliviensis F. O.
Pickard-Cambridge, 1897, only selected field and laboratory
observational studies about their diet exist. By using a DNA
metabarcoding approach, we compared the prey found in the gut content of
males and females from three distant Colombian populations of P.
boliviensis . By DNA metabarcoding of the cytochrome c oxidase
subunit I (COI), we detected and identified 234 prey records belonging to
96 operational taxonomic units (OTUs), as prey for this wandering
predator. Our results broaden the known diet of P. boliviensiswith at least 75 prey taxa not previously registered in fieldwork or
laboratory experimental trials. These results suggest that P.
boliviensis feeds predominantly on invertebrates (Diptera, Lepidoptera,
Coleoptera and Orthoptera) and opportunistically on small squamates.
Intersex and interpopulation differences are observed. Assuming that
prey preference does not vary between populations, these differences are
likely associated with a higher local prey availability. Finally, we
suggest that DNA metabarcoding can be used for evaluating subtle
differences in the diet of distinct populations of P.
boliviensis, particularly when predation records in the field cannot be
established or quantified using direct observation.
Keywords: Predator–prey interaction, gut content metabarcoding,
molecular diet analysis, diet, Phoneutria, prey detection.
INTRODUCTION
Understanding the contribution of predators in shaping the structure of
ecological communities is a central issue in ecology
(Lima 1998;
Schmitz 2007;
Seibold, et al. 2018). Therefore,
predator-prey interactions are a main driver of natural selection,
population dynamics, food web structure, community assemblages, and
ecosystem function (Portalier, et al.
2019; Severtsov and Shubkina 2015;
Start, et al. 2020). Spiders are among
the most abundant predators in terrestrial ecosystems, playing an
important role in the controlling prey species populations
(Betz and Tscharntke 2017;
Michalko, et al. 2019;
Pekár, et al. 2011;
Pusceddu, et al. 2018). Almost all
spiders are carnivores, feeding predominantly on arthropods including,
to a lesser extent, other spiders
(Birkhofer and Wolters 2012;
Nyffeler 1999;
Pekár and Toft 2015). Very rarely
non-arthropod prey are consumed as
a supplement to the arthropod diet (Foelix
2011; Nyffeler, et al. 2016;
Symondson, et al. 2002). However, this
appears to be relatively frequent in some spider families which include
large-sized species, such as Theraphosidae, Ctenidae, Lycosidae, and
Pisauridae, among others (Hazzi 2014;
Valdez 2020).
Predators in their search for food sources can develop narrow
(stenophagous) or broad (euryphagous) eating habits
(Pekár and Toft 2015). In the case of
spiders, the stenophagous specialists possess adaptations for the
capture of large focal prey, minimizing handling time, and live in
proximity to their prey so they can minimize foraging time. As a result,
the capture time is much shorter than in euryphagous species
(García, et al. 2018;
Michálek, et al. 2017;
Pekár, et al. 2011;
Pompozzi, et al. 2019). On the other
hand, euryphagous spiders, which select prey smaller than their body and
must thus capture more items of prey, could minimize foraging time by
shortening the duration of consumption of each item of prey
(García, et al. 2016;
Pompozzi, et al. 2019). Indeed, recent
studies showed that cursorial obligatory stenophagous species selected
larger prey and fed for a significantly longer time, extracting
significantly more mass than euryphagous spiders
(García, et al. 2018;
Michálek, et al. 2017).
Spiders of the genus Phoneutria , popularly known as “Banana
spiders” or “wandering spiders”, are restricted to South America and
they are essentially wandering, nocturnal spiders. This genus represents
one of the main groups of medically important spiders in South America
because of their defensive behavior, anthropogenic habits, and potent
venom (Hazzi 2014;
Valenzuela-Rojas, et al. 2019;
Vetter and Isbister 2008). Recent records
on the diet of P. boliviensis suggest these spiders prey on
several arthropod species but also consume vertebrates, mainly reptiles
and anurans (Valenzuela-Rojas, et al.
2019). In addition, mammals and birds have occasionally been reported
as prey in other Phoneutria species, suggesting they are likely a
generally euryphagous genus (Bücherl, et
al. 2013). Intersexual differences have been reported in the production
of venom of Phoneutria. For example, experiments with P.
nigriventer and P. boliviensis showed that there is greater
venom production in females than in males
(Estrada-Gomez, et al. 2015;
Herzig, et al. 2002;
Valenzuela-Rojas, et al. 2019), which
could indicate a possible difference in prey between the sexes. However,
recent laboratory observations on P. boliviensis showed no
difference in prey acceptance between males and females
(Valenzuela-Rojas, et al. 2019). In this
same study venom of males and females was found to be more effective
against vertebrate (geckos) and spider prey than other prey types.
Spiders feed on the pre-digested fluids of their prey through external
digestion and ingest nutrients only in liquid form. Following ingestion
of liquefied material through the esophagus and sucking stomach, the
food enters the midgut which branches into highly complex diverticula
extending throughout the prosoma and into the legs
(Foelix 2011;
Macías-Hernández, et al. 2018).
Consequently, digestion takes place in many different parts of the body
and dissecting the whole gut is nearly impossible due to its complexity
and extent (Macías-Hernández, et al.
2018). This condition restricts field studies of spider diet to direct
observation, stable isotope analysis, and molecular analysis
(Birkhofer, et al. 2017;
Jackson, et al. 2001;
Lafage, et al. 2020;
Pompanon, et al. 2012).
Molecular gut content analysis is a valuable tool for characterizing
trophic interactions. DNA-based methods possess a range of advantages
over other approaches in diet analysis
(Lafage, et al. 2020;
Sheppard and Harwood 2005). The inherent
bias and arduous process of taxonomic determination by field work
observation has led to an increased use of molecular techniques to
investigate the diet of spiders, of which DNA metabarcoding is currently
one of the most accurate and efficient for the analysis of generalist
polyphagous diets, achieving very good results in ecological studies on
multiple species of arachnids with euryphagous diets
(Piñol, et al. 2019). Likewise, DNA
metabarcoding workflows are intrinsically different from conventional
invertebrate identification approaches in that they use short DNA
sequences as a proxy for species detection
(Liu, et al. 2020).
We recently investigated the feeding preferences and related role of the
venom of the wandering spider P. boliviensis(Valenzuela-Rojas, et al. 2020;
Valenzuela-Rojas, et al. 2019), a widely
distributed spider species on the American continent inhabiting several
ecosystems in South America, spread across Bolivia, Paraguay, Perú,
Ecuador, Brazil, and Colombia
(Estrada-Gomez, et al. 2015;
Hazzi 2014). To further understand the
ecological and evolutionary interactions among trophic ecology, venom
composition and predatory behavior in this species, we study the diet.
The objective of this study, beyond knowing what P. boliviensiseats, was to test whether there are interpopulation and intersexual
differences in diet, through the technique of DNA metabarcoding of gut
contents. To our knowledge, this is the first study using this molecular
technique for detecting the diet of P. boliviensis .
MATERIALS AND METHODS
Collection and locations
Sixty adult specimens of P. boliviensis were used for DNA
metabarcoding of the entire gut contents. Twenty individuals (ten
females and ten males) of the spiders came from three Colombian
localities each, distant by an average of 300 km: Barbosa (Antioquia;
6°40’ 54.7”N, 75°41’ 10.4”W), Oporapa (Huila; 2°01’40.5”N,
75°59’43”W) and Ibagué (Tolima; 4°32’22.3”N, 75°05’37.1”W). Spiders
were collected from July 2019 to September 2019 (Figure 1). Locations
were selected based on the previous distribution reports for the species
in Colombia, and accessibility
(Estrada-Gomez, et al. 2015;
Valenzuela-Rojas, et al. 2020;
Valenzuela-Rojas, et al. 2019).
Selected adult individuals were opportunistically collected during one
night by randomly following cardinal transects of 500m between 20:00 and
04:00 h for three nights. However, most collection was done in what we
found to be the “peak hours” of activity of the species (between 00.00
– 03.00 h; JCV-R and JCG-G, Pers. Obs.). For each individual,
elevation, temperature, relative humidity and mass (g) was recorded
(Supplementary Table S1). At each locality, in addition to the 20
specimens, three more individuals were taken for standardization of the
DNA metabarcoding technique. After euthanized by freezing, the collected
individuals were stored in separate Falcon© tubes in 96% alcohol and
transported to the Biology Laboratory of the University of Ibagué
(Ibagué, Colombia). Subsequently, 96% alcohol washes were performed to
remove impurities. Later, the distal parts of the legs (tarsus and
metatarsus) were removed, since they do not contain gut diverticula
(Macías-Hernández, et al. 2018). In all
cases, dissection procedures were conducted using forceps and scissors,
flame sterilized after each dissection to prevent cross contamination.
We performed all laboratory activities on a clean and sterilized
laboratory bench. The dissected parts were then subjected to DNA
extraction separately (Macías-Hernández,
et al. 2018). Tissue samples were submitted to the standard procedures
of the Canadian Centre for DNA Barcoding (CCDB) for sequencing the
cytochrome c oxidase subunit I gene (COI), the standard marker
for the identification of most animals
(Pentinsaari, et al. 2020). We used the
Qiagen DNEasy Tissue kit for DNA extraction under manufacturer’s
conditions and used Bioline Plant TAQ to amplify PCRs.
Preliminary assays and samples processing
Since dissection of the highly diverticulated gut is difficult, we
performed three preliminary assays to determine which sampling of the
spider’s body (prosoma, prosoma+opisthosoma or the entire individual
except tarsus and metatarsus) yielded the most prey sequences
(Macías-Hernández, et al. 2018). We used
the primers dgHCO and mlCOIntF for PCR and sequencing
(Leray, et al. 2013). Once these
conditions were determined, DNA extraction was performed on the 60P. boliviensis samples from the three selected locations, using
the blocking primer pair noSPI blocking Primer
5‘TACACGACGCTCTTCCGATCTTCATTTYCCHCGWATAAAYAAYATAAG3‘ and dgHCO1
5‘CAGACGTGTGCTCTTCCGATCAGGAGTAAACTTCAGGGTGACCAAARAAYCA3‘.
Sequencing and DNA Metabarcoding
The metabarcoding of the P. boliviensis samples was performed
with the Illumina platform, where the pair-end-merging was done using
USEARCH v11.0.667 (Edgar 2010) and the
trimming processing was done in a custom program using Python 2.7.15 to
obtain the reads. The process of sequencing and quality control of the
reads was carried out by the company AIM (www.aimethods-lab.com)
following standard protocols (Kress and
Erickson 2007; Sang, et al. 1997).
Sequences with a length of less than 300pb were eliminated (with a
maximum expected error of 1%) using FastQC version 0.11.8
(de Sena Brandine and Smith 2019). In
addition, singletons, unique sequences and chimera sequences were
filtered using VSEARCH 2.9.1 (Rognes, et
al. 2016) to generate the final FASTQ files by sample, following the
protocols proposed by Leidenfrost et al and Liu et al
(Leidenfrost, et al. 2020;
Liu, et al. 2020). After quality
filtering, a total of 2,410,269 pair-end reads were obtained.
Sequences were blasted against the complete sequence database of the
Barcode of Life Data systems (BOLD) in order to find the closest matches
using the BOLD Identification Engine (http://www.boldsystems.org)
(Ratnasingham and Hebert 2007). Taxon
nomenclature follows the catalogue used in the BOLD and NCBI databases
(Federhen 2012). Based on the FASTQ files
for each individual, 256 operative taxonomic units (OTUs) were
identified.
Based on these OTUs, different filters according to standard exclusion
criteria for this technique were applied
(Deagle, et al. 2019;
Lafage, et al. 2020), in which sequences
with the following characteristics were eliminated: a) all reads
representing fewer than 0.01% of the total number of reads per sample
were removed from the analysis, b) sequences that corresponded to
environmental DNA or intestinal microbiota, and c) OTUs that did not
correspond to the geographical distribution of P. boliviensis .
Once these filters were applied, the OTUs were identified at the species
taxonomic level to those sequences with identity at ≥ 97%, genus with ≥
95%, family with 90% and order ≥ 75%. Additionally, BINs (Barcode
Index Numbers) were used to identify sequence clusters within the
database, correlating with species in 98% of all cases
(Lafage, et al. 2020). After OTU
filtering, a total of 105,583 sequences were obtained, corresponding to
96 OTUs. In order to not underestimate the total reads, a
standardization was made for each sample (percentage of reads per
sample).
Data analyses
To test differences in body mass (g) between the sexes and populations,
the mass of males and females were initially compared using the
Wilcoxon-Mann-Whitney Test with the package “coin”, and the mass
between populations using the Kruskal-Wallis Test with the package
“dplyr” in R version 4.0.2. In addition, a Wilcoxon post-hoc analysis
with Bonferroni adjustment was used (R
Development Core 2018), and the mass of the different spiders was
correlated with the estimated size of the different potential prey
obtained using the bugGuide database
(https://bugguide.net/node/view/15740 )
(Bartlett 2003). Table Supplementary S1).
The data on the diet of Phoneutria was analyzed separately by sex
and population. A Chi-square goodness-of-fit test was applied to compare
the ratio of number of reads between sexes and populations with Minitab
17 (Arend 2010). Additionally, the
confidence intervals between sexes under the binomial distribution were
corroborated with Statpages (https://statpages.info). Bipartite
prey–spider reads interaction network diagrams by sex and populations
were created with Sankeymatic
(http://sankeymatic.com/);
accessed on March 2020.
We also compared the richness of prey (i.e., OTUs that result after of
applying all filtering procedures for the reads) detected for each of
the three Colombian populations of P. boliviensis using
sample-based rarefaction curves (Hsieh, et
al. 2016) by calculating the Hill’s numbers by means the q parameter
(q0 = species richness) (Chao, et al.
2014). The confidence intervals for the obtained curves were determined
using the bootstrap method with 1000 resampling replicates. All analyses
and visualizations were performed in R version 3.5.2
(R Development Core 2018), using the
iNEXT package (Hsieh, et al. 2016).
- RESULTS
- Mass of spiders
Our results show that the females have a mass of 3.79 ± 0.21g (mean ±
standard error), with females weighing from 1.86 to 6.08g, while males
weigh 2.25 ± 0.17g, with individuals weighing from 1.01 to 4.64g. The
difference in mass between the sexes was highly significant (P
< 0.01) (Supplementary Table S1 and Figure S1, section a).
When analyzing the weight by locality, also significant differences were
found (P= 0.027). We found that spiders from the Ibagué sampling zone
tend to be larger (mean of 3.61g) compared to those from the Oporapa and
Barbosa populations (mean of 2.62 and 2.82g, respectively). These
results are summarized in the Supplementary Figure S1, section b. No
correlation (R2=0,003) was observed between fresh
spider mass (g) and prey size (mm).
DNA metabarcoding
analysis in Phoneutria boliviensis
Preliminary analyses with the blocking primers identified that the
optimal annealing temperature was 48°C. We also found that the
Prosoma+Opisthosoma region contained the higher prey stomach content. We
used these conditions for the processing of the 60 P. boliviensissamples from the three sampled locations for metabarcoding.
Of the 60 samples analyzed, and after applying all the filters, 57 had
at least one prey (for two Barbosa spiders and one Ibagué spider, no
prey was observed). We identified 96 different OTUs belonging to 10
different taxonomic orders (9 invertebrate and 1 vertebrate). Arthropods
contributed the largest portion of the diet (97.9% of prey DNA
recovered). The taxonomic orders with the greatest number of species
represented in the diet (77% of prey DNA recovered) included Diptera
(23), Coleoptera (20), Lepidoptera (17) and Orthoptera (14). The less
diversely represented orders were Phasmatodea (2), Hemiptera (2),
Araneae (2) and Dermaptera (1). By frequency of occurrence, Diptera made
up 23.9%, while Coleoptera made up (20.8%), Lepidoptera (17.7%) and
Orthoptera (14.6%). Vertebrates (Squamata) represented a small portion
of the diet, accounting for only 2% of the sequences obtained from gut
content analysis (Table 1).
The relative abundance analysis of reads indicates that between 2 and a
maximum of 64,287 reads were observed in each OTU. Of a total of 105,583
reads, 78.1% (82,537 reads) correspond to the order Orthoptera,
followed by Diptera with 7.3% (7,776) and Blattodea with 6.4% (6,841).
The remaining seven taxonomic orders represent only the remaining 8% of
the reads in the diet of P. boliviensis . However, our results
show a high dominance of reads for certain species. For example, the
cricket of the genus Neoconocephalus (sp5) had the highest number
of reads (64,287), which represented 60.8% of all reads (Supplementary
Table S2 and Figure S2).
A total of 234 prey were identified, belonging to 96 different species
in the 57 spiders sampled. The most abundant orders found as prey were
Diptera (62) and Orthoptera (59). Six orthopteran species belonging to
the genus Neoconocephalus were found, which were detected in 45
spiders, representing 19.2% of all registered prey items. For example,
the OTU Neoconocephalus sp.5 (Orthoptera) is present in the diet
of 32 of 57 spiders sampled. The preys with the second highest frequency
were the two species of flies belonging to the genus Delia , which
had a frequency of 9.8% of all prey (Table 1).
Intersexual differences in prey composition ofPhoneutria boliviensis
In this work, differences were identified in the consumption of prey
between sexes of P. boliviensis . Of the 105,583 reads analyzed,
28.9% (30,516) correspond to prey detected in females and 71.1%
(75,067) of the prey detected in males. For most orders, the number of
reads of prey in males was more than double that of prey detected in
females. Exclusive orders of a single genus were observed, as in the
case of Squamata in females or Hemiptera in males (Supplementary Table
S2). The prey–spider reads interaction network by sex is summarized in
the Supplementary Figure S2, section a.
The greatest difference in the number of reads by prey between the sexes
was found in the order Orthoptera. This order presents of 68.7% of
reads in females compared with 82% in males. On the other hand, the
results show us that of the total reads for this order (82,537), the
females represent 25.4% (20,961), while the males represent 74.6%
(61,576) and Ibagué represents 18.8% (15,517). When comparing the
proportion of reads in the other taxonomic orders between the sexes,
although they present significant differences between them (P= 0,000),
these are not as marked as those observed in Orthoptera. The
prey–spider reads interaction network by sex is summarized in the
Supplementary Figure S2, section a.
A total of 234 prey belonging to 96 different OTUs were observed, of
which 108 prey were detected in females and 126 in males. The
differences are mainly in the orders Orthoptera (females: 26 and males:
33), Lepidoptera (females: 7 and males: 17), Hemiptera (females: 0 and
males: 2) and Squamata (females: 2 and males: 0). We identified 60
different prey species in females, while 65 different prey species were
identified in males. Half (30) of the identified prey in the females are
unique, while the other half (30) are shared by both sexes. Thirty-five
species were observed as unique in males (Table 1). These results of the
prey–spider species interaction network by sex is summarized in Figure
2.
Interpopulation differences in prey composition ofPhoneutria boliviensis.
Differences in the composition of diet in P. boliviensis were
observed among the different populations. Of the 105,583 reads analyzed,
63.4% (66,970) correspond to prey detected in Barbosa, 13.1% (13,829)
in Oporapa and 23.5% (24,784) of the prey detected in Ibagué
(Supplementary Table S2). The relative abundance of reads per sample and
locality are shown in Figure 3.
The greatest difference between populations in reads per prey is
observed in the order Orthoptera. For example, with respect to the total
reads for each population, this order presented 90.3% in Barbosa,
46.9% in Oporapa and 62.5% in Ibagué. On the other hand, the results
showed that of the total reads for this order (82,537), the population
of Barbosa represents 73.4% (60,537), while Oporapa presented 7.8%
(6,487) and Ibagué 18.8%. (15,513).
Chi-square goodness-of-fit test
show significant differences between reads of populations (P= 0,000).
The prey–spider reads interaction network by population is summarized
in the Supplementary Figure S2, section b.
A total of 234 prey belonging to 96 different OTUs were observed, of
which 70 were detected in Barbosa, 83 in Oporapa and 81 in Ibagué. In
terms of prey richness, Oporapa had the highest number of prey OTUs
across individuals, followed by Ibagué, and Barbosa, which was supported
by the rarefaction curves (Figure 4). The differences were mainly
related to the orders Diptera (n=14, 20 and 28, respectively),
Orthoptera (n=17, 17 and 25) and Coleoptera (n=11, 24 and 13). In
Barbosa, 35.7% (n=25) of the species identified as prey are unique,
while in Oporapa it has 30.1% (n=25) and Ibagué with 23.5% (n=19). Our
results show that 21 to 26 spider prey species are shared between two or
three spider populations (Table 1). The results of the prey–spider
species interaction network by population is summarized in Figure 5.
- DISCUSSION
- Diet of Phoneutria boliviensis
In trophic ecology studies, it is not always possible to track trophic
links between predators and prey by direct observation. This is
especially critical when observing small, wandering, or elusive animals
with nocturnal or cryptic food‐web ecology
(Sheppard and Harwood 2005). In this
study, we developed a blocking primer pair (noSPI/dgHCO1) specifically
designed for metabarcoding of gut contents of the spider P.
boliviensis and phylogenetically closely related members of the
Ctenidae family. Most studies based on molecular gut-content analyses
have focused on small-sized spider families such as the Linyphiidae
(Agustí, et al. 2003;
Macías-Hernández, et al. 2018), Lycosidae
(Eitzinger, et al. 2019;
Lafage, et al. 2020;
Zhong, et al. 2019), Theridiidae and
Salticidae (Furlong, et al. 2014;
Whitaker, et al. 2019), Tetragnathidae
(Chapman, et al. 2010;
Toju and Baba 2018) and Oxyopidae
(Greenstone, et al. 2014). In these
studies, the procedure of extraction was either homogenizing the whole
spider or crushing the abdomen, taking into account the amount of sample
and the size of the spider. The performance of the primer pair
noSPI/dgHCO1 was compared in different body section samples of the
spider, with the Prosoma+Opisthosoma sample showing the highest content
of prey in P. boliviensis . Similar results were found in previous
studies of other spider species, showing that this sampling of the body
generally results in the largest prey content
(Lafage, et al. 2020;
Macías-Hernández, et al. 2018). However,
due to the use of blocking primes in our study, ctenid prey that were
previously reported as prey for Phoneutria , such asSpinoctenus , Ctenus and Phoneutria are blocked
(Valenzuela-Rojas, et al. 2020). This
might explain the low number of spiders we found as prey in our analyses
(see results).
Studies focusing on spider diet in natural habitats using metabarcoding
are scarce, particularly for spider species that do not build an aerial
web. Wirta et al. (Wirta, et al. 2015)
found Pardosa glacialis to consume mainly Diptera and Lepidoptera
in the High Arctic. Hambäck, et al.
(Hambäck, et al. 2016) also found Diptera
as the main prey of Pardosa prativaga followed by Lepidoptera,
Coleoptera, and Heteroptera in the Baltic shoreline. Lafage et al.
(Lafage, et al. 2020) found similar
results to those observed here (up to 12 orders and 117 species) in
Lycosidae species (genera tested: Pardosa, Trochosa, Hygrolycosa,
Xerolycosa , and Pirata ), using the same molecular technique.
They also found Diptera, Hemiptera, Hymenoptera, and Lepidoptera as the
main prey of different lycosid spiders. Similar results were also found
in the wolf spiders of the genus Hygrolycosa and in multiple
species of Pardosa , where Diptera are also the main component of
the diet (Lafage, et al. 2020; Zhong, et al. 2019).
In a recent study of gut content analysis of four spider species of the
genus Tetragnatha using DNA metabarcoding, Diptera, Hemiptera and
Lepidoptera were observed to be the most frequent in the diet of these
species in order of importance (Kennedy,
et al. 2019). These results are strongly related with their capture
strategy since the genus Tetragnatha builds webs to obtain its
food. However, this kind of spiders is very small and web-building and
therefore difficult to compare to the large, wandering P.
boliviensis .
In this work, we determined a wide diversity of potential prey (up to 10
orders and 96 species) consumed by P. boliviensis through DNA
metabarcoding analysis of gut contents; identifying prey belonging to
the orders Diptera, Coleoptera, Lepidoptera and Orthoptera, principally.
Our results show that P. boliviensis can consume, in addition to
invertebrates, vertebrate species such as Anolis sp. (lizard) andStenorrhina sp. (snake). Some vertebrates have previously
been reported as part of the P. boliviensis diet
(Valenzuela-Rojas, et al. 2020;
Valenzuela-Rojas, et al. 2019;
Viera and Gonzaga 2017). For example, of
the 57 spiders with gut contents identified, 52.6% (30) of them has at
least one Diptera species as prey. According to our results, frequent
consumption of some Diptera species could be considered as the basis of
the diet of this spider. However, the study by Valenzuela-Rojas et al.
(Valenzuela-Rojas, et al. 2020) reported
no captures of Diptera, hypothesizing that P. boliviensis has a
preference for larger prey. They identified prey up to three times
larger than the size of the spider itself, which is contrary to our
observations of a high preference for Diptera prey, which are small
species between 4 to 13 mm long (Bartlett
2003).
We propose that at least three hypotheses that could explain the
presence of Diptera as a prey for P. boliviensis .
a) Diptera as parasites of prey : Our analyses detected the
presence of Diptera parasitoids, as in the case of Anisia(Stireman, et al. 2019) or Ogcodessp. (Kerr and Winterton 2008) andTachina sp. (Sullivan and
Ozman-Sullivan 2012) including spider parasites as in the case ofTachina sp. (Gillung and Borkent
2017; Sullivan and Ozman-Sullivan 2012).
However, in our analysis, out of the 62 Diptera detected, only 3
(Anisia sp. (1), Ogcodes sp. (1) and Tachina sp.
(1); see Table 1) would be parasitoids of other insects. Therefore, it
does not lend support to a possible hypothesis that most of the Diptera
species we detected were erroneously interpreted as prey but were
actually present as parasitoids of arthropods preyed on by the spider.
Similarly, in our analyses, several species of Hymenoptera were detected
as prey. However, some of these species are reported as parasitoids,
such as Mesochorus sp. (Ashfaq, et
al. 2005) or predators of different insect species such asVespula sp. (Pusceddu, et al.
2018). Among the parasitic Hymenoptera families detected from the
spider samples, only Ichneumonidae includes well-characterized species
known to parasitize on spiders (Eberhard
2000; Gauld and Dubois 2006).
b) Diptera and or Hymenoptera as parasites of the Phoneutriathemselves : Despite the frequent presence of parasitic insects in
spiders, since most parasitoid hymenopterans and other parasitic Diptera
(Gillung and Borkent 2017) detected in
this study may be prey of the examined spiders; our spider samples did
not seem to carry eggs/larvae/pupae of spider-specific Ichneumonidae
wasps, which are easily recognizable ectoparasites
(Eberhard 2000;
Gauld and Dubois 2006).
c) Diptera (another insects) in the diet of the prey : We found a
small number of prey species that are predators, as in the case of the
spider Gonatium sp., which preys mainly on Coleoptera
(Tamaddoni-Nezhad, et al. 2013); or the
spider Neoscona sp. which preys mainly on Hemiptera (Xaaceph &
Butt, 2014). We also found predatory Coleoptera, as in the case ofAgonum sp., Aleochara lata , Bembidion sp.,Saprinus sp., which prey on different insects
(Battán Horenstein and Linhares 2011;
Caron, et al. 2008;
Riddick and Mills 1994;
von Berg, et al. 2009). Similarly, we
detected a species of Hemiptera of the genus Arma , which is
reported as a predator of several insects
(Zou, et al. 2012). Another explanation
is the possible cannibalistic behavior of P. boliviensis towards
juveniles or his own offspring, behavior reported in some scorpion
species. This enables adults to collect food from prey that are usually
too small for them (Polis 1980). If
Phoneutria do the same, this could account for the higher number of
small Diptera. Finally, within our samples of possible P.
boliviensis prey, two reptile prey (Anolis sp. andStenorrhina sp.) were detected, which are reported mainly as
insectivores (Barragán-Contreras and
Calderón-Espinosa 2013; Shewchuk and
Austin 2001). Of the 234 prey identified in P. boliviensis , only
24 (see Table 1) are insect predators or parasitoids; therefore, we
consider that a possible false prey event would be low, but cannot be
ruled out. Unfortunately, DNA metabarcoding analysis alone cannot
indicate whether such parasitoids are directly preyed on by spiders or
were just parasitizing on insect hosts caught by spiders
(Toju and Baba 2018). According to Toju
and Baba (Toju and Baba 2018), existing
DNA barcoding methods, including those presented in this work, are not
designed to distinguish between signs of cannibalism, intra-guild
predation, scavenging, and secondary predation (i.e., DNA of a prey’s
prey). Despite these limitations, DNA metabarcoding is a powerful tool
for systematically investigating prey compositions of a whole spider
community because many spider species are nocturnal, as P.
boliviensis (Bucaretchi, et al. 2016),
making their predation behavior difficult to observe directly in the
field (Fenk, et al. 2010;
Toju and Baba 2018).
d) Bias in previous studies : In a recent observational study on
the diet of P. boliviensis in the field, diet was composed
predominantly of arthropods (86%) and to a lesser extent, small
vertebrates (14%), consuming 21 different species (Valenzuela-Rojas, et
al. 2020). According to Valenzuela-Rojas, et al. (2020) the natural diet
of P. boliviensis included individuals of the orders Araneae,
Blattodea, Coleoptera, Hymenoptera, Lepidoptera, Mantodea, Orthoptera,
and Phasmatodea. Furthermore, the diet of this spider is shown to
include lizards (Gekkonidae and Sphaerodactylidae) and frogs (Hylidae).
According to the authors, the prey of P. boliviensis was larger
than the spiders themself, and predation was not observed of prey items
larger than three times the body length of P. boliviensis . Our
results, as well as the previous reports, confirm that P.
boliviensis has a generalist polyphagous diet. However, in the study
conducted by Valenzuela-Rojas et al.
(Valenzuela-Rojas, et al. 2020), contrary
to our results, no captures of Diptera were observed as a prey of this
spider. The difference between the observational studies and the studies
of analysis of gut content by DNA metabarcoding appears to be the
latter’s greater sensitivity in the detection of prey
(Liu, et al. 2020;
Pompanon, et al. 2012). This was
exemplified again here in the case of Diptera. DNA metabarcoding
presents the same limitations as observational fieldwork studies in
species identification, especially in tropical regions with high
biodiversity (Barsoum, et al. 2019;
Liu, et al. 2020). However, in DNA
metabarcoding, species identification can be automated in a
high-throughput way, given a sufficiently diverse reference database.
Likewise, Diptera (and other small-sized insects such as some
coleopteran species), may be missing from previous studies because they
were all carried out at night, while diptera are mostly diurnal.
Additionally, these insects are small, and consumed fast, reducing the
likelihood that they are encountered while being consumed.
Despite all the advantages of prey detection using DNA metabarcoding,
and the large number of arthropod species with DNA barcodes (246,069) in
the Barcode of Life Data Systems (BOLD)
(Ratnasingham and Hebert 2007), the low
representativity of sequences from South America (Colombia with only
39,741 of 9,265,546 Specimen Records; accessed June 2020) reduces the
likelihood that sequences are correctly assigned, which may lead to not
identifying sequences of 77 species (identified only at the order level)
of prey that were in the gut contents of the spider.
According to recent studies, DNA metabarcoding is mostly used to detect
the presence of species in samples rather than their relative abundance,
although sequence frequency is sometimes used as a proxy of species
abundance (Aizpurua, et al. 2018;
Deagle, et al. 2019). Several studies
have demonstrated positive relationships between input species biomass
and output sequencing reads. For example, Bista et al.
(Bista, et al. 2018) detected a strong
positive correlation (R 2 = 0.83) between
biomass input and read abundance by DNA metabarcoding for the beetleGyrinus marinus . Therefore, based on our analyses of relative
abundance of reads of prey, the most abundant taxonomic groups were
Orthoptera, Diptera and Blattodea, respectively (see Supplementary Table
S2 and Figure S2), which could mean that the highest biomass associated
with stomach content in P. boliviensis would be associated mainly
with the cricket of the genus Neoconocephalus and the species of
Diptera detected as prey. On the other hand, according to Sint et al.
(Sint, et al. 2011;
Sint, et al. 2015), the molecular
detection of prey in spiders has a high sensitivity up to 84-h
post-feeding; which would explain both the high number of prey per
spider detected in this work (234), and the variation in the reads of
the same prey species (i.e., the cricket of the genusNeoconocephalus , detected in 37/57 spider’ guts).
Differences between males and females
In most animals, males and females show marked differences in primary
and secondary sexual traits. The sexual dimorphism literature pertaining
to invertebrates is fragmented, particularly for arachnids
(McLean, et al. 2018). Spiders are
sexually dimorphic in various morphological, behavioral, and life
history traits (Cordellier, et al. 2020).
Our study provides evidence for both sexual dimorphism (significant
differences in mass between females and males), and a different trophic
ecology in males and females of P. boliviensis . We identified
that the males on average are smaller than the females; however, we
found that some males are larger than some females. Adult male and
female spiders generally differ in body weight, but they can be of
similar size and shape or differ markedly, depending on species
(Cordellier, et al. 2020). These analyses would suggest that a
differential trophic environment could affect this character (see above
for our analysis by location).
We found sexual differences in diet in our analysis of gut contents ofP. boliviensis . Our results show that of the total reads, 28.9%
are present in females and 71.1% detected in males. We also found a
lower number of preys in females than in males (108 and 126,
respectively). In a study on wolf spiders, differences in predatory
behavior were observed between males and females
(Persons and Uetz 1999). In this study,
it was observed that male wolf spiders travel greater distances, with
greater hunting efforts compared to females. In addition, other factors
such as the amount of venom produced and toxicity of the venom differ
between the sexes in this species
(Valenzuela-Rojas, et al. 2019). This in
turn could influence, or be an adaptation to, differences in the diet,
for example the consumption of vertebrates by females (Table 1).
Our results show that wild female P. boliviensis prey on
vertebrates. Different studies under laboratory conditions have shown
the tendency of P. boliviensis to be a euryphagus given its great
voracity, as well as its capacity to consume small vertebrates
(Valenzuela-Rojas, et al. 2019). Through
the analysis of gut contents by metabarcoding, these findings are
corroborated in the wild.
The difference in the number of intakes of prey or variety of food items
is usually attributed to differences in size, either total size or
specific structures such as pedipalps or carapace
(Foellmer and Moya-Larano 2007). Body
size dimorphism may be the result of selection for many factors, such as
reproductive success, hyper-predation, dispersal capacity among others
(Crawley 2009). In our results, males are smaller than females and also
have a greater number of prey items than females, which supports the
hypothesis that in spiders with sexual dimorphism a larger size does not
necessarily imply that a greater range of prey can be captured
(McLean, et al. 2018). Similarly, the
capacity to consume certain prey is influenced by multiple factors, not
only morphological differences between sexes, but also differential
predatory behavior and the active search for females by males to
reproduce (Kotiaho, et al. 1998). The
toxicity of P. boliviensis venom is higher in males than in
females, which could ensure greater success in hunting, without ruling
out the possibility that it is used defensively (Valenzuela-Rojas, et
al. 2019).
Differences between populations
Although the number of prey items among the Barbosa, Oporapa and Ibagué
populations are very similar (70, 83 and 81, respectively), differences
exist in the type of prey. When comparing the diet in P.
boliviensis among populations, the diversity (number of species) of
Coleoptera was higher in Oporapa, representing twice as many items as
the other two populations studied. Furthermore, it was the only
population where the order Dermaptera was consumed, while in Ibagué the
most consumed order was Diptera followed by Orthoptera. Finally, Barbosa
was the locality with the lowest number of prey species. The three
localities present differences in the prey species that compose the
diet, which may indicate that the selection of prey species may be in
good part be determined by their availability
(Birkhofer and Wolters 2012;
Eitzinger, et al. 2019;
Hambäck, et al. 2016). Nevertheless, the
analysis of relative read counts by population shows that spiders from
Barbosa have two or three times more reads than the populations of
Oporapa and Ibagué. One hypothesis that could explain this behavior
would be that larger spiders have a greater prey capacity. However, our
results do not identify a clear correlation (R2=0,003)
between read count per population and the body size of the spider (in
term of fresh mass).
The proportions of each order of prey show marked differences between
populations, especially in the order Orthoptera. These differences
between populations may be influenced by the availability of prey in
each population, and these in turn are influenced by biogeographical
characteristics such as the Andes mountain range that crosses the
populations and could affect the ranges of certain prey species, as well
as factors such as temperature, humidity and anthropogenic pressures
(Ramirez-Villegas, et al. 2014).
According to the environmental data of relative humidity, temperature
and altitude, some differences are observed between the populations (see
Supplementary Table S1), which would justify the differences between
them at the trophic level.
CONCLUSIONS
This project contributed to the knowledge of the trophic ecology ofP. boliviensis using the DNA metabarcoding approach, confirming
its Euryphagous feeding behavior. We show a wider range of prey species
than reported previously for this species in Colombia. Some prey species
are reported here for the first time. Our results also could be the
first evidence that males of P. boliviensis apparently have a
greater hunting effort, as indicated by an increased number of prey and
reads compared to females. This suggests different predatory strategies
between the sexes, perhaps based on different energy requirements.
Similarly, there is a small difference in trophic ecology depending on
the locality (comparing our records of diet composition and taxa
richness), again confirming the generalist nature and flexibility of the
diet. However, our results differ from previous findings using different
diet assessment techniques in the field and the laboratory
(Valenzuela-Rojas et al., 2019, 2020). Particularly the high fraction of
Diptera in the diet we found could not have been anticipated based on
these previous findings and requires further corroboration. In addition,
it is necessary to extend the studies with DNA barcoding in South
America, in order to identify, at the species level, the possible prey
and the ecological impact of this predatory spider.