Research Article
Print
Research Article
Invasive terrestrial invertebrate detection in water and soil using a targeted eDNA approach
expand article infoCecilia Villacorta-Rath, Lori Lach, Natalia Andrade-Rodriguez, Damien Burrows, Dianne Gleeson§, Alejandro Trujillo-González§
‡ James Cook University, Townsville, Australia
§ University of Canberra, Canberra, Australia
Open Access

Abstract

Terrestrial invasive invertebrates can rapidly colonise new areas, causing detrimental effects on biodiversity, economy and lifestyle. Targeted environmental DNA (eDNA) methods could constitute an early detection tool given their sensitivity to small numbers of individuals. We hypothesised that terrestrial runoff would transport eDNA from the land into adjacent waterbodies and used the invasive yellow crazy ant (Anoplolepis gracilipes) as a model species to test this hypothesis. We collected water samples from four waterbodies adjacent (< 10 m from the creek edge) to infestations following rainfall events for eDNA analysis. We also collected soil samples from areas of known infestations and tested five eDNA extraction methods to determine their efficiency to extract eDNA from soil. Water samples resulted in positive yellow crazy ant eDNA amplification (20–100% field replicates across all sites), even at one site located 300 m away from where ants had been detected visually. Soil samples resulted in a higher percentage of false negatives when sampled from ant transit areas than from nest entrances. Unpurified DNA extracts from soil also resulted in false negative detections and only after applying a purification step of DNA extracts, did we detect yellow crazy ant eDNA in 40–100% of field replicates across all methods and sites. This is the first study to empirically show that eDNA from a terrestrial invertebrate can be successfully isolated and amplified from adjacent or downstream waterbodies. Our results indicate that eDNA has the potential to be a useful method for detecting terrestrial invertebrates from soil and water.

Keywords

Biosecurity, eDNA runoff, false negative, invasive species detection, qPCR inhibition, terrestrial eDNA, yellow crazy ant

Introduction

Over the past 15 years, environmental DNA (eDNA) analysis has gained momentum for biomonitoring of both marine and freshwater aquatic systems (Zaiko et al. 2018; Sepulveda et al. 2020; Trujillo-González et al. 2021). Targeted-eDNA methods (i.e. quantitative real-time polymerase chain reaction [qPCR]-based) are now considered a sensitive, specific and robust tool for detection of aquatic or semi-aquatic invasive species (Ficetola et al. 2008; Piaggio et al. 2014; Smart et al. 2015; Villacorta-Rath et al. 2020). This has resulted in eDNA methods being increasingly adopted into monitoring programmes by natural resource management agencies, consultancy companies and citizen-science groups (Darling and Mahon 2011; Larson et al. 2020). However, advances in the use of eDNA to detect aquatic species have not been matched by developments for terrestrial species.

A key challenge in using targeted eDNA as a method to detect terrestrial species is determining an effective sampling strategy (Taberlet et al. 2012; Van Der Heyde et al. 2020). This needs to consider where eDNA is likely to be deposited by the target organism, the effect of the substrate on eDNA detectability and ease of sampling. To date, a few studies have used soil samples for targeted eDNA detection (Kucherenko et al. 2018; Katz et al. 2021; Yasashimoto et al. 2021) and most of them have reported false negatives arising from sampling soil areas that the species did not occupy (Kucherenko et al. 2018; Yasashimoto et al. 2021). Similarly, DNA binds to soil to a varying level depending on its physicochemical composition and inappropriate eDNA extraction methods can lead to qPCR inhibition, affecting eDNA detectability and data interpretation (Andersen et al. 2012; Katz et al. 2021; Yasashimoto et al. 2021).

Invasive ants are amongst the most harmful invasive species globally (Kenis et al. 2009; Siddiqui et al. 2021); they threaten the environment, human health and livelihoods and social amenities (Holway et al. 2002; Lach and Hooper-Bùi 2010; Gruber et al. 2022). Millions of dollars are spent on prevention, treatment and control of ant invasions globally to avoid these impacts (Zenni et al. 2021). Current methods for invasive ant detection (i.e. baited traps or cards, pitfall traps and detection dogs) rely on luring, trapping, smelling or sighting active individuals or colonies, which can be labour-intensive, costly and reliant on species behaviour and weather conditions (Hoffmann et al. 2010; Lach and Barker 2013; L. Lach, unpublished data). These methods have low detection sensitivity to low numbers of individuals, which can lead to false negatives (Stringer et al. 2011; L. Lach, pers obs.). Environmental DNA analysis could improve detectability of invasive ants, since it does not require sighting the target species (Jerde et al. 2011). However, studies investigating terrestrial insect eDNA capture for species detection are seldomly available and are designed in agricultural contexts, wherein their sampling strategy consisted of spraying water over crops to aggregate the available eDNA deposited there (Valentin et al. 2018, 2020; Allen et al. 2021).

In this study, we test different methods to capture and detect terrestrial invertebrate eDNA. We use the yellow crazy ant, Anoplolepis gracilipes (Smith, 1857), one of the most environmentally and socioeconomically damaging invasive insect species in the world (Clarke et al. 2021; Gruber et al. 2022), as a model species. We collected water and soil samples in field conditions to: (1) investigate whether yellow crazy ant eDNA can be detected from creeks or rivers adjacent to existing infestations; and (2) to compare the efficiency of laboratory extraction methods for eDNA from soil samples. We hypothesised that terrestrial runoff would transport eDNA from the land into adjacent waterbodies.

Materials and methods

Study system

The yellow crazy ant is a widespread invader in tropical regions, particularly the Indo-Pacific (Janicki et al. 2016). Colonies consist of multiple interacting nests, usually termed ‘supercolonies’ (Abbott 2005; Hoffmann 2014). In forested areas, nests are typically at the base of trees or under rocks, leaf litter or logs, but the ants will nest virtually anywhere with the right temperature, humidity and protection from sunlight (e.g. discarded car engines and soft drink cans) and can readily relocate nests when disturbed (Hoffmann 2015; Lach, pers. obs.). Workers may also shelter temporarily in nest-like sites without brood or queens, for unknown periods of time (Lach, pers. obs.). Their nesting and foraging behavioural plasticity make choosing appropriate sites from which to sample soil for eDNA more challenging than it may be for other species, such as Argentine ants, Linepithema humile (Yasashimoto et al. 2021). The first recorded infestation of the yellow crazy ant in mainland Australia was in the Northern Territory in the early 1990s (Majer 1984). Yellow crazy ant incursions have been reported in Queensland, Australia, since 2001, including in the Cairns and Townsville regions (Lach and Hoskin 2015).

Townsville is in the “dry tropics” region of Australia and is characterised by a wet (November to April) and a dry (May to October) season. During the wet season, ambient temperature ranges between 23.2 °C (± 1.4) and 30.9 °C (± 0.7), there is a mean humidity of 65.6% (± 5.1) and a mean rainfall of 169.3 mm (± 102.7) (http://www.bom.gov.au/climate/averages/tables/cw_032040.shtml). During the dry season, ambient temperature ranges between 16.6 °C (± 2.6) and 27 °C (± 1.6), the mean humidity is 58.2% (± 5.7) and mean precipitation is 19.9 mm (± 8.6) (http://www.bom.gov.au/climate/averages/tables/cw_032040.shtml).

Environmental DNA sample collection

Water samples for eDNA analysis were collected and preserved from waterbodies adjacent (< 10 m from the creek edge) to yellow crazy ant infestations in Townsville, Queensland, Australia, with “high activity” (M. Green, Townsville City Council, pers. obs.) in February and March 2021 (Ross River, unnamed creek adjacent to Chauncy Crescent, Stuart Creek) (Table 1, Fig. 1), during or immediately after rainfall events. We also sampled from an unnamed creek (Palmetum site) in which yellow crazy ants had been detected 300 m from the creek edge (M. Green, Townsville City Council, pers. obs.). At each site, five replicate 30 ml surface water samples were collected using a sterile 50 ml Falcon tube and decanted into another 50 ml Falcon tube containing 10 ml of Longmire’s preservative solution (Longmire and Baker 1997) as per Villacorta-Rath et al. (2021). Replicate water samples were taken approximately 10 m apart. At every site, a field blank was also carried out to ensure that the process of sample collection did not introduce contamination. The field blank consisted of decanting 30 ml of laboratory-grade water into a Falcon tube containing 10 ml of preservative solution.

Table 1.

Field sites where water and soil samples were collected for eDNA analyses in the Townsville area, north Queensland. All samples were collected by the authors, except for water samples from Stuart Creek, which were collected by the Townsville City Council.

Site Latitude, Longitude # eDNA replicates Collection date Collection time Sample volume per replicate
Water eDNA samples
Ross River 19.3127°S, 146.7563°E 5 26/02/2021 morning 30 mL
Creek at Chauncy Crescent 19.3126°S, 146.7575°E 5 26/02/2021 morning 30 mL
Creek at Palmetum 19.3114°S, 146.7631°E 5 26/02/2021 morning 30 mL
Stuart Creek 19.3405°S, 146.8533°E 5 8/03/2021 morning 30 mL
Soil eDNA samples – Ant transit sampling event
Gieseman Road 19.2680°S, 146.5784°E 7 10/12/2020 morning 1 g
Copper Refinery 19.3411°S, 146.8530°E 7 10/12/2020 morning 1 g
Soil eDNA samples – Nest entrance sampling event
Gieseman Road 19.2680°S, 146.5784°E 10 27/04/2021 morning 1 g
Douglas 19.3126°S, 146.7575°E 10 26/04/2021 dusk 1 g
Figure 1.

Sampling sites where eDNA samples were collected in the Townsville region, Queensland, Australia. Water sampling sites included: (1) Ross River; (2) creek at Chauncy Crescent; (3) creek at Palmetum; and (4) Stuart Creek. Soil sampling sites included: (5) Gieseman Road; (6) Copper Refinery; and (7) Douglas. Map Data: Google 2023 CNES / Airbus, Maxar Technologies.

Soil samples were collected during two other sampling events in Townsville at known infestation sites (Table 1, Fig. 1). Soil samples were not collected at the same time as water samples because ants tend to retreat into their nests during rainfall events, when water samples were collected. Soil sampling took place in the morning and dusk and between the months of December and April, when, due to the environmental conditions of the sampling region, yellow crazy ant activity is likely to be high (Hoffmann 2015). During the first sampling event (December 2020; Table 1), samples were collected from areas where ants could be observed transiting (hereafter referred to as ant transit samples). Seven replicate samples were collected from two sites along 50–70 m transects, starting where ant activity was observed. One of the sites (Gieseman Road site; Table 1, Fig. 1) was adjacent (< 10 m away) to a water body (Stuart Creek, Table 1, Fig. 1). Sampling consisted of collecting 2 ml of soil into a sterile 50 ml Falcon collecting tube. Subsequently, 1 ml of this soil was transferred into a tube containing 700 µl cetyltrimethylammonium bromide (CTAB) buffer from the CTAB protocol (Adamkewicz and Harasewych 1996) and the other 1 ml into a tube containing 2 ml of Biomeme lysis buffer from the M1 Bulk Sample Prep Kit for DNA – High Concentration (Philadelphia, Pennsylvania). Two field controls were collected to assess potential field contamination: one was a tube containing 700 µl CTAB buffer and the other was a Biomeme homogenisation tube containing 2 ml of Biomeme lysis buffer.

In the second sampling event (April 2021; Table 1), we collected soil from suitable microhabitats (bases of trees and under fallen logs) that had visible evidence of high yellow crazy ant traffic and high incidence of dead ants, which were considered putative nest entrances (hereafter referred to as ant nest entrance samples). In this case, ten 50 ml replicate soil samples were collected from two sites. One of the sites (Gieseman Road, Table 1) was resampled; however, the other previously sampled site (Copper Refinery) had recently been treated with insecticide, so samples were collected from a different site (Douglas: Table 1). Similar to the Gieseman Road site, the Douglas site was adjacent (< 10 m away) to a water body (Ross River, Table 1). Sampling consisted of collecting the equivalent of 50 ml of soil in sterile plastic containers, carefully avoiding visible dead ants. Containers were stored in ice and immediately transported to the laboratory at James Cook University. Each sample was then shaken to homogenise and approximately 1 g of soil was partitioned into four sub-samples for extraction method comparisons and stored in DNA LoBind 2 ml tubes at -80 °C until eDNA was extracted. Field controls consisted of five tubes containing the lysis buffers of each eDNA extraction method (described below).

DNA extraction, quantitative real-time polymerase chain reaction and sequencing

Environmental DNA extraction from water

We extracted eDNA from water samples at the dedicated TropWATER eDNA laboratory, James Cook University, Townsville, Australia. We followed the preserve, precipitate, lyse, precipitate and purify (PPLPP) method in Edmunds and Burrows (2020). Briefly, the PPLPP workflow uses a glycogen-aided isopropanol-based precipitation, followed by a guanidinium hydrochloride and TritonX-based lysis and a subsequent glycogen-aided polyethylene glycol (PEG)-based precipitation. For each eDNA extraction batch, an extraction control was added to ensure that no contamination was introduced during laboratory procedures. This consisted of decanting 30 ml of laboratory-grade water into a Falcon tube containing 10 ml of preservative solution in the laboratory. DNA extracts were then purified to remove environmental inhibitors using the Qiagen DNeasy PowerClean Pro Cleanup Kit (Germantown, Maryland), as per the usual procedure of eDNA water sample analysis (Villacorta-Rath et al. 2020; Villacorta‐Rath et al. 2021).

Environmental DNA extraction from soil

Ant transit soil samples were extracted using two methods: the field-based extraction and the laboratory-based method. The field-based extraction method involved using the Biomeme M1 Bulk Sample Prep Kit for DNA – High Concentration (Philadelphia, Pennsylvania) and the laboratory-based method involved using a chloroform-based extraction protocol (CTAB; Adamkewicz and Harasewych (1996)). For the field-based extraction, the equivalent of 1 ml of soil sample was added to a Biomeme homogenising tube with 2 ml of Biomeme lysis buffer and mixed vigorously using vortex mixing for one minute. After this time, the supernatant was drawn from the tube using the Biomeme extraction column and syringe and pumped out 20 times. Extraction methods were then followed as per the manufacturer’s instructions and eluted in 100 µl of Biomeme Elution buffer. For the laboratory-based extraction, the equivalent of 1 ml of soil sample was directly transferred into a microcentrifuge tube containing 700 µl CTAB buffer and 10 µl proteinase K was added to samples upon arrival to the laboratory following Adamkewicz and Harasewych (1996). Samples were vortexed, crushed with plastic pestles and left to lyse at 65 °C in a hibernation oven on a rocking platform for 20 hours. After lysis was performed, 700 µl chloroform-isoamyl was added and samples were centrifuged at 16,000 g for 10 min. The supernatant was then transferred into a new tube containing 600 µl chloroform-isoamyl and centrifuged again at 16,000 g for 10 min. The resulting supernatant was transferred into a new tube containing 600 µl of cold isopropanol, inverted to mix and stored at -20 °C overnight. After freezing, samples were centrifuged at 16,000 g at 4 °C for a total of one hour and the supernatant was pipetted off taking care to not lose the formed pellet. The pellet was then washed with 1,000 µl of 70% ethanol and centrifuged at 16,000 g for 10 min. Finally, the ethanol was pipetted off and pellets were allowed to air dry for 15 min, before re-suspending them in 100 µl TE buffer and stored at 4 °C. DNA extracts from both extraction methods were tested for the presence of contaminants using a NanoDrop spectrophotometer.

Ant nest entrance soil samples were extracted using five methods: (1) the Biomeme M1 Bulk Sample Prep Kit for DNA – High Concentration (Philadelphia, Pennsylvania); (2) the CTAB method; (3) the PPLPP method (Edmunds and Burrows 2020); (4) the Qiagen Dneasy PowerSoil kit (Germantown, Maryland); and (5) the modular-universal DNA extraction (Mu-DNA) method (Sellers et al. 2018).

For the Biomeme M1 Bulk Sample Prep Kit for DNA – High Concentration (Philadelphia, Pennsylvania), we modified the previously-described procedure in the first and last steps: the equivalent of 1 ml of soil samples were initially transferred into 3 ml of Biomeme lysis buffer and DNA extracts were eluted in 400 µl Biomeme Elution buffer. Similarly, the first step of the CTAB method was modified, wherein each 1 ml-field replicate was split into two 2 ml DNA LoBind tubes (the equivalent of approximately 0.5 ml of soil sample/tube) containing 1000 µl CTAB buffer.

For the PPLPP method, each replicate of 1 ml soil sample was transferred into a 50 ml DNA LoBind Falcon tube containing 10 ml Longmire’s buffer (Longmire and Baker 1997) and 10 ml of MilliQ water. Environmental DNA was extracted following Edmunds and Burrows (2020) with eDNA eluted in 100 µl elution buffer.

For the Qiagen Dneasy PowerSoil kit (Germantown, Maryland) (from hereon referred to as ‘Qiagen method’), each field replicate consisting of the equivalent of 1 ml soil was partitioned into four tubes with the equivalent of 250 ml of soil and mixed with 60 µl Solution C1. The bead beating step was not performed given the target was not bacterial DNA from the soil samples. We then followed the manufacturer’s protocol handling each field replicate in four separate tubes, sequentially passed through a single spin column and eluted in 100 µl elution buffer.

Finally, we followed the soil sample workflow of the Mu-DNA protocol without the bead beating step. Each equivalent of 1 ml soil replicate was split into two 2 ml DNA LoBind tubes and mixed with 550 µl lysis solution, 200 µl soil lysis additive and 20 µl proteinase K. Samples were then vortexed and incubated for 3 h at 55 °C. Subsequently, samples were centrifuged at 4,000 g for 1 min, the supernatant was transferred into a new tube, centrifuged at 10,000 g for 1 min and the supernatant was again transferred into a new tube containing 300 µl flocculant solution. Extraction was then carried out as published in Sellers et al. (2018). A negative extraction control was added to each batch of eDNA extractions to ensure that no contamination was introduced during laboratory procedures.

Quantitative real-time polymerase chain reaction (qPCR) and cycle sequencing

We screened samples for yellow crazy ant eDNA presence using two probe-based, species-specific eDNA assays developed and optimised by EcoDNA, targeting two separate sections of the yellow crazy ant Cytochrome Oxidase 1 (COI) gene region: Agra1 assay (112 base pair [bp] long) and Agra2 assay (131 bp long) (Suppl. material 1).

qPCR plates were set-up using the Arise Biotech EzMate 401 Automated Pipetting System (Taipei, Taiwan) and run in a Thermo Fisher Scientific QuantStudio 5 Real-Time PCR System (Singapore). We tested four technical replicates of each sample and each site, including field and extraction blanks, three no-template controls and genomic DNA positive controls. Each qPCR reaction and cycling conditions were as explained in the assay development section of this study (Suppl. material 1). Inhibition was tested using a Thermo Fisher Scientific TaqMan Exogenous Internal Positive Constrol (IPC) qPCR assay (Burlington, Ontario). A total of 3 µl IPC was applied to duplicate samples and three reactions containing only IPC were included as controls. A departure of three or more cycles (ΔCt) would indicate sample inhibition (Hartman et al. 2005). Samples were subsequently purified using the Qiagen Dneasy PowerClean Pro Cleanup Kit (Germantown, Maryland) and another qPCR was carried out. Ct values between unpurified and purified samples were compared to assess the level of sample inhibition. A subset of amplicons with positive detections were Sanger-sequenced for confirmation of results at the Australian Genome Research Facility (AGRF). Resulting sequences were deposited in GenBank (accession numbers: MZ330820MZ330832).

Data analysis

Differences in yellow crazy ant eDNA capture sensitivity (number of DNA copies per assay) across different methods were assessed with a generalised linear mixed model using a template model builder (TMB) computed in the R package glmmTMB version 1.7.19 (Brooks et al. 2021). The response variable was the DNA copy number and the explanatory variables were the eDNA extraction method (fixed effect) and field replicate/technical replicate (nested random effects, with technical replicates nested within field replicates). Two models were run: the first one testing replicates as an additive fixed factor and the second one testing replicates as an interaction. The best performing model was chosen, based on the corrected Akaike Information Criterion (AICc). We tested for overdispersion with the DHARMa R package version 0.4.4 (Hartig 2021). Post-hoc paired comparisons of means were performed using Tukey’s HSD. Given that the differences in mean eDNA concentration across laboratory extraction methods were large, data were log-transformed solely for comparison (Fig. 3). Statistical analyses were completed in R statistical environment (R Development Core Team 2021).

Data accessibility

All data supporting the findings of this study can be found under the Suppl. materials and archived in the James Cook University Research Data Hub.

Results

Yellow crazy ant eDNA detection in water and soil

Water samples collected adjacent or in the vicinity of yellow crazy ant infestations showed positive eDNA amplification with both assays. The highest percentage of eDNA detections were observed at Stuart Creek (100% of field and technical replicates using the Agra2 assay), followed by Ross River (80% and 60% of field and technical replicates using the Agra2 assay, respectively), Palmetum Creek (100% and 50% of field and technical replicates using the Agra2 assay, respectively) and Chauncy Crescent Creek (20% and 10% of field and technical replicates using the Agra2 assay, respectively) (Fig. 2). The Agra2 eDNA assay amplified DNA extracts from water samples at a higher percentage than the Agra1 assay, which failed to detect yellow crazy ant eDNA at Chauncy Crescent Creek (Fig. 2).

Figure 2.

Percentage of positive yellow crazy ant eDNA detections from purified water samples using the Agra1 and Agra2 eDNA assays, targeting two different fragments of the COI gene.

In soil samples, the Agra1 assay was more efficient at amplifying DNA extracts. Samples collected from ant transit sites showed high concentration of contaminants (Suppl. material 2: table S1) and complete inhibition was observed in eDNA samples tested using both eDNA assays (Suppl. material 2: table S2). After purification of DNA extracts, both assays successfully amplified 43–64% of the laboratory-extracted samples (Table 2, Suppl. material 2: table S3). However, in-situ extracted samples had a low percentage of positive detections (14–18%) (Table 2, Suppl. material 2: table S3).

Table 2.

Yellow crazy ant eDNA detection in purified soil samples using Agra1 and Agra2 eDNA assays, targeting two different fragments of the COI gene. Extraction methods used were field-based (Biomeme M1 Bulk Sample Prep Kit for DNA – High Concentration [Philadelphia, Pennsylvania]), CTAB (cetyltrimethylammonium bromide, Adamkewicz and Harasewych 1996), PPLPP (preserve, precipitate, lyse, precipitate, purify method, Edmunds & Burrows, 2020), Qiagen (Qiagen DNeasy PowerSoil kit [Germantown, Maryland]), and Mu-DNA (modular-universal DNA extraction method, Sellers et al. 2018).

eDNA assay Site Extraction method # Field replicates % Positive field replicates # Technical replicates % Positive technical replicates
Soil samples – Transit sites sampling event
Agra1 Gieseman Road Field-based 7 29 28 18
CTAB 7 86 28 64
Copper Refinery Field-based 7 29 28 14
CTAB 7 71 28 50
Agra2 Gieseman Road Field-based 7 4 28 14
CTAB 7 57 28 43
Copper Refinery Field-based 7 29 28 14
CTAB 7 57 28 43
Soil samples – Nest entrance sampling event
Agra1 Gieseman Road Field-based 10 70 40 45
CTAB 10 100 40 90
PPLPP 10 100 40 83
Qiagen 10 80 40 78
Mu-DNA 10 50 40 38
Douglas Field-based 10 40 40 28
CTAB 10 80 40 68
PPLPP 10 90 40 78
Qiagen 10 70 40 48
Mu-DNA 10 40 40 23
Agra2 Gieseman Road Field-based 10 50 40 35
CTAB 10 100 40 93
PPLPP 10 100 40 100
Qiagen 10 100 40 95
Mu-DNA 10 70 40 58
Douglas Field-based 10 40 40 28
CTAB 10 80 40 70
PPLPP 10 70 40 70
Qiagen 10 70 40 68
Mu-DNA 10 50 40 35

There was a greater percentage of positive technical replications from nest entrance samples compared to transit samples for both assays and for field based and CTAB purification (Table 2, Suppl. material 2: table S4). In general, the Agra1 assay was more efficient at amplifying yellow crazy ant eDNA, reflected in higher eDNA yields across all extraction methods than the Agra 2 assay (Suppl. material 2: table S4). In addition, the CTAB, PPLPP and Qiagen extraction methods had the highest percentage of positive detections when using both assays (Suppl. material 2: table S4), whereas eDNA yield (number DNA copies/assay) was variable depending on the sampling site and eDNA assay used (Fig. 3, Suppl. material 2: table S4). Field and extraction controls from both the ant transit sites and nest entrance sampling events did not show positive amplification.

Figure 3.

Mean DNA concentration (mean DNA copy number per assay ± Standard Error) yielded by each of the five eDNA extraction methods from purified soil samples collected from likely yellow crazy ant nest entrances (ant nest entrance sampling event) and run using the: (a) Agra1 assay on Gieseman Road samples, (b) Agra2 assay on Gieseman Road samples, (c) Agra1 assay on Douglas samples and (d) Agra2 assay on Douglas samples. DNA yield was log10 transformed. Note the differences in y-axis scales across panels. Methods with different letters above the bars within each panel differ significantly in post-hoc tests (at P < 0.05) using Tukey HSD. Extraction methods used were field-based (Biomeme M1 Bulk Sample Prep Kit for DNA – High Concentration [Philadelphia, Pennsylvania]), CTAB (cetyltrimethylammonium bromide, Adamkewicz and Harasewych (1996)), PPLPP (preserve, precipitate, lyse, precipitate, purify method, Edmunds and Burrows (2020)), Qiagen (Qiagen DNeasy PowerSoil kit [Germantown, Maryland]) and Mu-DNA (modular-universal DNA extraction method, Sellers et al. (2018)).

At Ross River, where we collected both water and soil samples, yellow crazy ant eDNA detections of both substrata were similar: 40–90% of positive soil field replicates (Agra1 assay) compared to 80% of positive water field replicates (Agra2 assay) and 28–78% of positive soil technical replicates (Agra1 assay) compared to 60% of positive water field replicates (Agra2 assay).

Comparison of eDNA yield across extraction methods

There were significant differences between mean number of DNA copies across all eDNA extraction methods at both sites and using both assays (Fig. 3). For Gieseman Road samples tested using Agra1 assay, the field-based method yielded significantly fewer DNA copies than the CTAB (F = -3.006, P = 0.0008), PPLPP (F = -3.888, P < 0.001) and Qiagen (F = -2.345, P = 0.0136) methods (Fig. 3). Additionally, the Mu-DNA method yielded significantly fewer DNA copies than the CTAB (F = 2.198, P = 0.0345) and PPLPP (F = -3.080, P = 0.0008) methods (Fig. 3). When using the Agra2 assay, we found that the CTAB method yielded significantly higher number of DNA copies than the field-based (F = -2.590, P = 0.0018), Mu-DNA (F = 4.295, P < 0.0001), PPLPP (F = 1.711, P = 0.0003) and Qiagen (F = 1.886, P = 0.0001) methods (Fig. 3). In addition, the Mu-DNA method had significantly higher number of DNA copies than the PPLPP (F = -2.583, P = 0.0001) and Qiagen (F = -2.409, P = 0.0004) methods (Fig. 3).

For Douglas samples, the Qiagen method exhibited significantly higher number of DNA copies than the field-based (F = -7.477, P = 0.0002), CTAB (F = -4.399, P = 0.0115) and PPLPP (F = -6.051, P < 0.0001) methods of samples tested using the Agra1 assay (Fig. 3). When using the Agra2 assay, the field-based method yielded significantly fewer DNA copies than the PPLPP (F = -2.8819, P = 0.0036) and Qiagen (F = -3.5810, P = 0.0118) methods (Fig. 3); and the Mu-DNA method had significantly fewer DNA copies than the CTAB (F = 2.4420, P = 0.0077), PPLPP (F = -2.9663, P = 0.0039) and Qiagen (F = -3.6653, P = 0.0002) methods (Fig. 3).

Discussion

Detection methods that are sensitive to small number of individuals, such as eDNA analysis, have the capacity to complement and improve the detection of invasive species (Jerde et al. 2011; Smart et al. 2015; Villacorta-Rath et al. 2020). We used yellow crazy ants as a model species to test eDNA detection of a terrestrial invasive species in water and soil, as well as to explore the effect of soil eDNA extraction methods on eDNA detectability. Both substrata yielded positive yellow crazy ant eDNA detection. Importantly, we detected yellow crazy ant eDNA in water samples from creeks and rivers directly adjacent and in the vicinity of known infestations. To the best of our knowledge, our findings are the first demonstration of the feasibility of detecting terrestrial invertebrate eDNA in natural waterways. Additionally, we found that eDNA detectability in soil is dependent on the extraction method and the area from which the samples are collected (i.e. ant transit areas vs. nest entrances) and that purification of DNA extracts is necessary to avoid false negative detections.

Many aquatic eDNA studies show that population size (Yates et al. 2019; Spear et al. 2021; Yates et al. 2021), time-elapsed since a target organism has occupied a system (Kucherenko et al. 2018; Schmidt et al. 2021) and target species behaviour (Buxton et al. 2017; Dunn et al. 2017), influence eDNA detectability in water. In the present study, qualitative data from the infestations adjacent to our water sampling sites suggested that ant activity (behaviour) was high in all sampled areas, except next to Palmetum Creek, where 80% of the field replicates exhibited yellow crazy ant eDNA. On the other hand, at Chauncy Crescent Creek, adjacent to an area of high ant activity, only 20% of the field replicates resulted in positive eDNA detections. Therefore, future research to disentangle the factors related to ant density and activity that could affect eDNA detectability in water will be useful. We also expect that the amount of rainfall prior to water sampling and the distance between the infestations and receiving waterbodies will play an important role in transporting terrestrial eDNA into the aquatic system. Once yellow crazy ant eDNA is in the aquatic system, we hypothesise that the amount of rainfall will influence eDNA detectability by increasing water flow and dilution. As with aquatic eDNA studies, we would expect that the time elapsed since eDNA transport into aquatic systems would determine eDNA detectability due to factors affecting aquatic eDNA detection (i.e. eDNA production, decay, transport, retention and resuspension; see Barnes et al. (2016)). Future studies should focus on investigating four main factors that could influence yellow crazy ant eDNA detectability in waterbodies, namely: (1) total area of the infestation; (2) ant activity; (3) time since the establishment at a site; and (4) amount of rainfall prior to water sampling.

In soil samples, eDNA detectability from areas of yellow crazy ant transit was lower than that of ant nesting areas. A recent study on Argentine ants eDNA also found the highest eDNA concentration in soil from nest entrances, as opposed to surface soil samples from an infestation area and found no relationship between eDNA concentration and distance from nests or trails (Yasashimoto et al. 2021). The authors argued that Argentine ants may move nests frequently and, therefore, strong relationships between eDNA concentration and distance from a nest were not expected (Yasashimoto et al. 2021). Yellow crazy ants also move nests and transfer brood to different locations frequently (Lach pers. obs.) and, as with most ants, move dead ants to immediately outside of the nests (Hölldobler and Wilson 1990). Therefore, it is at the nest entrances where we would expect a significant amount of eDNA to be deposited. There is also the possibility that higher detections at nest entrance sites are due to small ant parts present in soil samples, even though we avoided sampling dead ants. If eDNA methods are used to check the progress of eradication efforts, this could constitute a source of false positive detections. Therefore, it would be important to investigate how to avoid the potential of false positive detections arising from dead yellow crazy ants. If the aim is to detect presence of the species in a new area, we propose the highest likelihood of collecting yellow crazy ant eDNA is in soil from samples at the base of trees or other moist areas where they are more likely to establish long-term nests. Yasashimoto et al. (2021) also concluded that the type of ant activity and their behaviour at different areas will determine eDNA detectability, indicating the importance of understanding the ecology of the species to avoid false negative detections.

Soil type may have also affected detectability. Samples collected from Gieseman Road, which has coarse sandy soils (Murtha 1975), showed a higher percentage of positive detections than those from Douglas, which has clay soils (Murtha 1975). Regardless of the eDNA assay used, soils with higher percentage of organic matter or clay and higher pH (i.e. more negatively charged) are more likely to bind to eDNA (Allemand et al. 1997; Andersen et al. 2012) and, therefore, inhibit the qPCR reaction. Therefore, we would expect to have more effective eDNA extraction from the coarser soil from Gieseman Road compared to the more organic-rich soil from Douglas (Murtha 1975), which is shown by the higher percentage of positive eDNA detections found at the former.

Our results showed that column-based eDNA extraction methods (Qiagen and Mu-DNA) perform better at removing sample inhibition than the other three methods, which only showed eDNA amplification after a purification step. This means that the purification step could be avoided, cutting laboratory costs and shortening the sample processing time. In terms of eDNA yield, the Qiagen method was more or equally as effective in recovering eDNA from soil than the CTAB and PPLPP. Although the Mu-DNA method was less efficient than Qiagen, it can be scaled up to any starting volume of soil and it is almost ten times more cost-effective than the latter (Sellers et al. 2018).

Conclusions

In the present study, we used yellow crazy ants as a model species to explore eDNA detectability in two different substrata: water and soil. We demonstrated that terrestrial eDNA can be detected in waterbodies near yellow crazy ant infestations. Our findings suggest that there are opportunities for detecting terrestrial invertebrate eDNA across large areas given that mechanisms, such as rainfall runoff, could aggregate eDNA into nearby or downstream waterbodies. However, factors influencing terrestrial invertebrate eDNA detectability in water should be explored further. We showed that detectability of eDNA in soil is dependent on sampling location and the eDNA extraction method and that purification of DNA extracts is important to avoid false negative detections, making soil sampling less attractive than water sampling.

Acknowledgements

We acknowledge the Traditional Owners of the land where the sampling and laboratory analyses were conducted. The contribution of C.V.R. and N.A.R. were supported through funding from the Australian Government’s National Environmental Science Program—Northern Australia Environmental Resources Hub, Project 4.3 (to D.B.) and the Department of Agriculture, Water and Environment (to C.V.R. and L.L.). In addition, the contribution of A.T.G. was funded by the Biosecurity Innovation Program, from the Australian Government Department of Agriculture, Water and the Environment (to D.G. and A.T.G.). Thanks to Angela Strain (JCU) and the Wet Tropics Management Authority (WTMA) for providing yellow crazy ant specimens for assay development. We thank Melissa Green (Townsville City Council), as well as Bev Job and Janet Cross (Invasive Species Council) for assistance with sample collection and discussions on locations of yellow crazy ant infestations in Townsville.

References

  • Abbott KL (2005) Supercolonies of the invasive yellow crazy ant, Anoplolepis gracilipes, on an oceanic island: Forager activity patterns, density and biomass. Insectes Sociaux 52(3): 266–273. https://doi.org/10.1007/s00040-005-0800-6
  • Adamkewicz SL, Harasewych MG (1996) Systematics and biogeography of the genus Donax (Bivalvia: Donacidae) in Eastern North America. American Malacological Bulletin 13: 97–103.
  • Allen MC, Nielsen AL, Peterson DL, Lockwood JL (2021) Terrestrial eDNA survey outperforms conventional approach for detecting an invasive pest insect within an agricultural ecosystem. Environmental DNA 3(6): 1102–1112. https://doi.org/10.1002/edn3.231
  • Andersen K, Bird KL, Rasmussen M, Haile J, Breuning-Madsen H, Kjær KH, Orlando L, Gilbert MTP, Willerslev E (2012) Meta-barcoding of “dirt” DNA from soil reflects vertebrate biodiversity. Molecular Ecology 21(8): 1966–1979. https://doi.org/10.1111/j.1365-294X.2011.05261.x
  • Brooks ME, Kristensen K, van Benthem KJ, Magnusson A, Berg CW, Nielsen A, Skaug HJ, Mächler M, Benjamin M (2021) glmmTMB balances speed and flexibility among packages for zero-inflated generalized linear mixed modelling. The R Journal 9(2): 378–399. https://doi.org/10.32614/RJ-2017-066
  • Buxton AS, Groombridge JJ, Zakaria NB, Griffiths RA (2017) Seasonal variation in environmental DNA in relation to population size and environmental factors. Scientific Reports 7(1): 1–9. https://doi.org/10.1038/srep46294
  • Clarke DA, Palmer DJ, McGrannachan C, Burgess TI, Chown SL, Clarke RH, Kumschick S, Lach L, Liebhold AM, Roy HE, Saunders ME, Yeates DK, Zalucki M, McGeoch MA (2021) Options for reducing uncertainty in impact classification for alien species. Ecosphere 12: e03461. https://doi.org/10.1002/ecs2.3461
  • Darling JA, Mahon AR (2011) From molecules to management: Adopting DNA-based methods for monitoring biological invasions in aquatic environments. Environmental Research 111(7): 978–988. https://doi.org/10.1016/j.envres.2011.02.001
  • Dunn N, Priestley V, Herraiz A, Arnold R, Savolainen V (2017) Behavior and season affect crayfish detection and density inference using environmental DNA. Ecology and Evolution 7(19): 7777–7785. https://doi.org/10.1002/ece3.3316
  • Edmunds RC, Burrows D (2020) Got glycogen?: Development and multispecies validation of the novel preserve, precipitate, lyse, precipitate, purify (PPLPP) workflow for environmental DNA extraction from Longmire’s preserved water samples. Journal of Biomolecular Techniques 31: 125–150. https://doi.org/10.7171/jbt.20-3104-003
  • Gruber MAM, Santoro D, Cooling M, Lester PJ, Hoffmann BD, Boser C, Lach L (2022) A global review of socio-economic and environmental impacts of ants reveals new insights for risk assessment. Ecological Applications 32(4): e2577. https://doi.org/10.1002/eap.2577
  • Hoffmann BD (2014) Quantification of supercolonial traits in the yellow crazy ant, Anoplolepis gracilipes. Journal of Insect Science 14(1): 25. https://doi.org/10.1093/jis/14.1.25
  • Hoffmann BD (2015) Integrating biology into invasive species management is a key principle for eradication success: The case of yellow crazy ant Anoplolepis gracilipes in northern Australia. Bulletin of Entomological Research 105(2): 141–151. https://doi.org/10.1017/S0007485314000662
  • Janicki J, Narula N, Ziegler M, Guénard B, Economo EP (2016) Visualizing and interacting with large-volume biodiversity data using client-server web-mapping applications: The design and implementation of antmaps.org. Ecological Informatics 32: 185–193. https://doi.org/10.1016/j.ecoinf.2016.02.006
  • Katz AD, Harper LR, Sternhagen EC, Pearce SE, Melder CA, Sperry JH, Davis MA (2021) Environmental DNA is effective in detecting the federally threatened Louisiana Pinesnake (Pituophis ruthveni). Environmental DNA 3(2): 409–425. https://doi.org/10.1002/edn3.126
  • Kenis M, Auger-Rozenberg MA, Roques A, Timms L, Péré C, Cock MJW, Settele J, Augustin S, Lopez-Vaamonde C (2009) Ecological effects of invasive alien insects. Biological Invasions 11(1): 21–45. https://doi.org/10.1007/s10530-008-9318-y
  • Kucherenko A, Herman JE III, Iii EME, Urakawa H (2018) Terrestrial snake environmental DNA accumulation and degradation dynamics and its environmental application. Herpetologica 74(1): 38–49. https://doi.org/10.1655/Herpetologica-D-16-00088
  • Lach L, Hoskin CJ (2015) Too much to loose: Yellow crazy ants in the wet tropics. Wildlife Australia 52: 37–41.
  • Larson ER, Graham BM, Achury R, Coon JJ, Daniels MK, Gambrell DK, Jonasen KL, King GD, LaRacuente N, Perrin-Stowe TIN, Reed EM, Rice CJ, Ruzi SA, Thairu MW, Wilson JC, Suarez AV (2020) From eDNA to citizen science: Emerging tools for the early detection of invasive species. Frontiers in Ecology and the Environment 18(4): 194–202. https://doi.org/10.1002/fee.2162
  • Longmire JL, Baker RJ (1997) Use of “lysis buffer” in DNA isolation and its implications for museum collections. Museum of Texas Tech University, Lubbock, Texas, 1–3. https://doi.org/10.5962/bhl.title.143318
  • Majer JD (1984) Recolonisation by ants in rehabilitated opencut mines in northern Australia. Reclamation and Revegetation Review 2: 279–298.
  • Murtha GG (1975) Soils and land use on the northern section of the Townsville coastal plain, north Queensland. Soils and Land Use Series 55: 5–70.
  • Piaggio AJ, Engeman RM, Hopken MW, Humphrey JS, Keacher KL, Bruce WE, Avery ML (2014) Detecting an elusive invasive species: A diagnostic PCR to detect Burmese python in Florida waters and an assessment of persistence of environmental DNA. Molecular Ecology Resources 14(2): 374–380. https://doi.org/10.1111/1755-0998.12180
  • Schmidt KJ, Soluk DA, Maestas SEM, Britten HB (2021) Persistence and accumulation of environmental DNA from an endangered dragonfly. Scientific Reports 11(1): 1–8. https://doi.org/10.1038/s41598-021-98099-1
  • Sellers GS, Di Muri C, Gómez A, Hänfling B (2018) Mu-DNA: A modular universal DNA extraction method adap for a wide range of sample types. Metabarcoding and Metagenomics 2: 1–11. https://doi.org/10.3897/mbmg.2.24556
  • Sepulveda AJ, Nelson NM, Jerde CL, Luikart G (2020) Are environmental DNA methods ready for aquatic invasive species management? Trends in Ecology & Evolution 35(8): 668–678. https://doi.org/10.1016/j.tree.2020.03.011
  • Siddiqui A, Bamisile BS, Khan MM, Islam E, Hafeez M, Bodlah I, Xu Y (2021) Impact of invasive ant species on native fauna across similar habitats under global environmental changes. Environmental Science and Pollution Research International 28(39): 54362–54382. https://doi.org/10.1007/s11356-021-15961-5
  • Smart AS, Tingley R, Weeks AR, Van Rooyen AR, McCarthy MA (2015) Environmental DNA sampling is more sensitive than a traditional survey technique for detecting an aquatic invader. Ecological Applications 25(7): 1944–1952. https://doi.org/10.1890/14-1751.1
  • Spear MJ, Embke HS, Krysan PJ, Vander Zanden MJ (2021) Application of eDNA as a tool for assessing fish population abundance. Environmental DNA 3(1): 83–91. https://doi.org/10.1002/edn3.94
  • Stringer LD, Suckling M, Baird D, Vander Meer RK, Christian SJ, Lester PJ (2011) Sampling efficacy for the red imported fire ant Solenopsis invicta (Hymenoptera: Formicidae). Environmental Entomology 40(5): 1276–1284. https://doi.org/10.1603/EN11002
  • Taberlet P, Prud’Homme SM, Campione E, Roy J, Miquel C, Shehzad W, Gielly L, Rioux D, Choler P, Clément JC, Melodelima C, Pompanon F, Coissac E (2012) Soil sampling and isolation of extracellular DNA from large amount of starting material suitable for metabarcoding studies. Molecular Ecology 21(8): 1816–1820. https://doi.org/10.1111/j.1365-294X.2011.05317.x
  • Trujillo‐González A, Villacorta‐Rath C, White NE, Furlan EM, Sykes M, Grossel G, Divi UK, Gleeson D (2021) Considerations for future environmental DNA accreditation and proficiency testing schemes. Environmental DNA 3(6): 1049–1058. https://doi.org/10.1002/edn3.243
  • Valentin RE, Fonseca DM, Nielsen AL, Leskey TC, Lockwood JL (2018) Early detection of invasive exotic insect infestations using eDNA from crop surfaces. Frontiers in Ecology and the Environment 16(5): 265–270. https://doi.org/10.1002/fee.1811
  • Valentin RE, Fonseca DM, Gable S, Kyle KE, Hamilton GC, Nielsen AL, Lockwood JL (2020) Moving eDNA surveys onto land: Strategies for active eDNA aggregation to detect invasive forest insects. Molecular Ecology Resources 20(3): 746–755. https://doi.org/10.1111/1755-0998.13151
  • Van Der Heyde M, Bunce M, Wardell‐Johnson G, Fernandes K, White NE, Nevill P (2020) Testing multiple substrates for terrestrial biodiversity monitoring using environmental DNA metabarcoding. Molecular Ecology Resources 20(3): 732–745. https://doi.org/10.1111/1755-0998.13148
  • Villacorta-Rath C, Adekunle AI, Edmunds RC, Strugnell JM, Schwarzkopf L, Burrows D (2020) Can environmental DNA be used to detect first arrivals of the cane toad, Rhinella marina, into novel locations? Environmental DNA 2(4): 635–646. https://doi.org/10.1002/edn3.114
  • Villacorta‐Rath C, Hoskin CJ, Strugnell JM, Burrows DW (2021) Long distance (> 20 km) downstream detection of endangered stream frogs suggests an important role for eDNA in surveying for remnant amphibian populations. PeerJ 9: e12013. https://doi.org/10.7717/peerj.12013
  • Yasashimoto T, Sakata MK, Sakita T, Nakajima S, Ozaki M, Minamoto T (2021) Environmental DNA detection of an invasive ant species (Linepithema humile) from soil samples. Scientific Reports 11(1): 1–12. https://doi.org/10.1038/s41598-021-89993-9
  • Yates MC, Fraser DJ, Derry AM (2019) Meta‐analysis supports further refinement of eDNA for monitoring aquatic species‐specific abundance in nature. Environmental DNA 1(1): 5–13. https://doi.org/10.1002/edn3.7
  • Yates MC, Glaser DM, Post JR, Cristescu ME, Fraser DJ, Derry AM (2021) The relationship between eDNA particle concentration and organism abundance in nature is strengthened by allometric scaling. Molecular Ecology 30(13): 3068–3082. https://doi.org/10.1111/mec.15543
  • Zaiko A, Pochon X, Garcia-Vazquez E, Olenin S, Wood SA (2018) Advantages and limitations of environmental DNA/RNA tools for marine biosecurity: Management and surveillance of non-indigenous species. Frontiers in Marine Science 5: 322. https://doi.org/10.3389/fmars.2018.00322
  • Zenni RD, Essl F, García-Berthou E, McDermott SM (2021) The economic costs of biological invasions around the world. In: Zenni RD, McDermott S, García-Berthou E, Essl F (Eds) The economic costs of biological invasions around the world. NeoBiota 67: 1–9. https://doi.org/https://doi.org/10.3897/neobiota.67.69971

Supplementary materials

Supplementary material 1 

Invasive terrestrial invertebrate detection in water and soil using a targeted eDNA approach

Cecilia Villacorta-Rath, Lori Lach, Natalia Andrade-Rodriguez, Damien Burrows, Dianne Gleeson, Alejandro Trujillo-González

Data type: methods on assay development

Explanation note: This supplementary file contains the methods involved in yellow crazy ant eDNA assays development.

This dataset is made available under the Open Database License (http://opendatacommons.org/licenses/odbl/1.0/). The Open Database License (ODbL) is a license agreement intended to allow users to freely share, modify, and use this Dataset while maintaining this same freedom for others, provided that the original source and author(s) are credited.
Download file (34.77 kb)
Supplementary material 2 

Additional qPCR results

Cecilia Villacorta-Rath, Lori Lach, Natalia Andrade-Rodriguez, Damien Burrows, Dianne Gleeson, Alejandro Trujillo-González

Data type: Quantitative PCR results

Explanation note: This file contains additional data on qPCR results from eDNA samples.

This dataset is made available under the Open Database License (http://opendatacommons.org/licenses/odbl/1.0/). The Open Database License (ODbL) is a license agreement intended to allow users to freely share, modify, and use this Dataset while maintaining this same freedom for others, provided that the original source and author(s) are credited.
Download file (92.42 kb)
login to comment