Streamlining DNA Barcoding Protocols: Automated DNA Extraction and a New cox1 Primer in Arachnid Systematics

Background DNA barcoding is a popular tool in taxonomic and phylogenetic studies, but for most animal lineages protocols for obtaining the barcoding sequences—mitochondrial cytochrome C oxidase subunit I (cox1 AKA CO1)—are not standardized. Our aim was to explore an optimal strategy for arachnids, focusing on the species-richest lineage, spiders by (1) improving an automated DNA extraction protocol, (2) testing the performance of commonly used primer combinations, and (3) developing a new cox1 primer suitable for more efficient alignment and phylogenetic analyses. Methodology We used exemplars of 15 species from all major spider clades, processed a range of spider tissues of varying size and quality, optimized genomic DNA extraction using the MagMAX Express magnetic particle processor—an automated high throughput DNA extraction system—and tested cox1 amplification protocols emphasizing the standard barcoding region using ten routinely employed primer pairs. Results The best results were obtained with the commonly used Folmer primers (LCO1490/HCO2198) that capture the standard barcode region, and with the C1-J-2183/C1-N-2776 primer pair that amplifies its extension. However, C1-J-2183 is designed too close to HCO2198 for well-interpreted, continuous sequence data, and in practice the resulting sequences from the two primer pairs rarely overlap. We therefore designed a new forward primer C1-J-2123 60 base pairs upstream of the C1-J-2183 binding site. The success rate of this new primer (93%) matched that of C1-J-2183. Conclusions The use of C1-J-2123 allows full, indel-free overlap of sequences obtained with the standard Folmer primers and with C1-J-2123 primer pair. Our preliminary tests suggest that in addition to spiders, C1-J-2123 will also perform in other arachnids and several other invertebrates. We provide optimal PCR protocols for these primer sets, and recommend using them for systematic efforts beyond DNA barcoding.

preliminary tests suggest that in addition to spiders, C1-J-2123 will also perform in other arachnids and several other invertebrates. We provide optimal PCR protocols for these primer sets, and recommend using them for systematic efforts beyond DNA barcoding.
The objective of our study was to explore an optimal strategy for extracting and analyzing arachnid DNA focusing on the barcoding and adjacent cox1 regions. Our work focused on the species richest arachnid lineage, spiders.
Our first goal was to improve an automated DNA extraction protocol. Compared with manual extraction procedures using kits, robotic DNA extraction methods often yield lower quantity of extracted DNA [34]. To maximize its efficiency, we experimentally adjusted an internal robotic DNA extraction program and improved it for acquisition of high concentration of genomic DNA from different quality tissues. Our second goal was to test the performance of commonly used cox1 primer combinations and to identify the optimal primer set over the major phylogenetic lineages of spiders. We screened and tested the high throughput utility with a single PCR program of ten cox1 primer pairs. Our third goal was to develop a new cox1 primer that would produce an indel-free alignment resulting in more accurate phylogenetic analyses.

Specimens and taxonomic coverage
Fifteen spider species were selected to represent all major spider clades [14,35] ( Table 2; Fig. 2). Specimens were obtained from the EZ Lab tissue bank (http:// Table 1. Common cox1 primers used in arachnid systematics, and tested in this study.
doi:10.1371/journal.pone.0113030.g002 pH 7.4); into row B: 120 mL of Wash Solution I, into row C, E, F: 120 mL of Wash Solution II, into row D: 38 mL of nuclease-free water, into row G: 40 mL of Elution Buffer I. During the run in the MagMAX Express Magnetic Particle Processor the magnetic beads solution (6.4 mL DNA binding beads with 9.6 mL nuclease-free water) into row A and 2 mL of RNase A into row D were added. In second pause 40 mL of Multisample DNA Lysis Buffer and 48 mL of Isopropanol were added into row D. During the third pause a step of incubation in thermoblock at 70˚C for 5 min was made. After the incubation, 40 mL of Elution Buffer II was added into row G (from 70 mL down to 30 mL minimum is allowed) and the run continued in the instrument. Samples of purified DNA were transferred to cryovials for storage from row G (see Appendix S1). The additional step of overnight incubation of starting material with Proteinase K was added to the protocol improving extraction efficiency. Differently sized tissue was cut and thoroughly homogenized with a pestle in a tube with 73.6 mL PK Buffer and 6.4 mL Proteinase K (100 mg/mL), shortly centrifuged and incubated over night at 55˚C on a shaker. All reagents and buffers (with exception of PBS) used for DNA extraction are components of MagMAX DNA Multi-Sample Kit (Life Technologies).
During the optimization step with MagMAX DNA Multi-Sample Kit also the comparison with MagMAX Total Nucleic Acid Isolation Kit (Life Technologies) was done (data not shown). After the quantification of extracted nucleic acids with NanoDrop Lite (Thermo Scientific) the amount of extracted DNA was up to 5-fold higher in comparison to sample concentration prepared with MagMAX DNA Multi-Sample Kit. However, better extraction of nucleic acids could also be related to the MagMAX Total Nucleic Acid Isolation Kit specifications where the isolation of all the nucleic acids (ssDNA, dsDNA and RNA) is performed at once. The comparison of material costs showed substantial differences between the kits used in this study, with the MagMAX DNA Multi-Sample Kit being the most economic.

PCR optimization strategy
For PCR amplification reaction we achieved the best results using GoTaq Flexi Polymerase Kit (Promega) and the cycling program for cox1 gene. All of the PCR reaction mixtures had a total volume of 25 mL and 1 mL of DNA template (in the range between 3.0 to 50 ng) per reaction was generally used. Each reaction included 5.1 mL of Promega's GoTaq Flexi Buffer, 0.14 mL of GoTaq Flexi Polymerase, 2.5 mL dNTP's (2 mM each), 2.3 mL MgCl 2 (25 nM), 0.5 mL of each primer (forward and reverse, 20 mM), 0.14 mL BSA (10 mg/mL) and 12.8 mL of sterile distilled water. The thermocycle program for cox1 amplification consisted of 95˚C for 1 min, 5 cycles of 94˚C for 40 sec, 45˚C for 40 sec and 72˚C for 1 min. Additional 35 cycles of 94˚C for 40 sec, 51˚C for 40 sec and 72˚C for 1 min with a final extension at 72˚C for 5 min were used. Two samples were amplified with a second round of PCR using the excised band as a template. All reagents and buffers used for PCR were supplied from Promega.

PCR amplification verification and sequencing
PCR products were stained with Sybr Safe (Invitrogen), separated by standard 1.5% agarose gel using Owl B2 EasyCast Mini Gel Electrophoresis System and visualized by UV gel imager E-box VX2/20LM (Vilber Lourmat). All of the PCRs were repeated three times to verify strong signals obtained in every reaction where specific PCR fragment was generated. Primers (Table 1)  . Sequence data were edited and assembled using Geneious Pro version 5.4 [36] and further handled in Mesquite version 2.74 [37].

Primer design
In order to construct a new primer that would bind upstream of the C1-J-2183 primer binding site, we searched for the most conserved region in that area. We used the NCBI nucleotide search facility within Geneious Pro to gather all cox1 sequences of the order Araneae that contained the keyword ''BARCODE''. This search tagged sequences meeting all CBOL criteria [38]. All sequences longer or shorter than 658 bp were deleted and the remaining 1672 sequences, already aligned, were used for primer design (see Appendix S2). Potential primers were evaluated using the program Primer3 [39].

Results
We optimized and improved the manufacturer's protocol for extraction of genomic DNA using the MagMAX Express Magnetic Particle Processor, an automated high throughput DNA extraction system. We processed a wide range of spider tissue of different taxonomic affiliation, size and quality, to improve the protocol and increase the efficiency of the procedure. The manufacturer's protocol only specifies the use of the kit with MagMAXExpress-96 Standard Magnetic Particle Processor or manually. Our procedure describes the use of the MagMAX Express Magnetic Particle Processor in combination with the MagMAX DNA Multi-Sample Kit. Additionally, we optimized the workflow for smaller quantities of starting material and accordingly adjusted and modified the internal program. Changes were made to the volume of reagents used and timing of specific steps (see Appendix S1).
To assess the amplification success of ten primer pairs routinely employed for barcoding cox1 gene, we tested fifteen target species throughout the spider phylogeny [14,35] (Fig. 2; Table 2) and performed PCR reactions using ten primer pairs (Table 2). Our goal was to compare the performance of each primer pair against every representative across the tree and evaluate their success rate. The success rate of different primer pairs varied from 30 to 100% ( Fig. 3; Table 2). To target solely the short barcoding cox1 region, we recommend using the combination of Folmer primers (LCO1490/HCO2198) [33].
To maximize sequence data for genetic, genomic, and phylogenetic analyses, however, a longer stretch of cox1 is desired. Existing primers can fail to provide a continuous stretch if used in a single pair or in combinations of pairs [12,[31][32]. The forward primer C1-J-2183, for example, is designed too close to the reverse  HCO2198 for accurate chromatogram reads, and in practice the resulting interpretations of base pairs result in two partial cox1 sequences (Fig. 1). We therefore designed, via analysis of the consensus alignment of 1672 arachnid cox1 sequences a new forward DNA primer situated 60 base pairs upstream of the C1-J-2183 binding site. The sequence GATCGAAATTTTAATACTTCTTTTTTTGA was chosen as the most conserved and appropriate for the binding of a new primer, named C1-J-2123 (Fig. 4). Our preliminary tests (data not shown) suggest that this primer will work not only in spiders, but also other arachnids (scorpions, mites and ticks) and other invertebrates (bivalves, gastropods, tunicates and others).
C1-J-2123 performed with the same success rate (93%) as the alternative primer C1-J-2183, amplifying in 14 out of 15 spider species (Fig. 3; Table 2). We recommend using the C1-J-2123/C1-N-2776 primer pair extended in upstream direction, which will allow for full overlap of this extended sequence with that obtained with the standard Folmer primers LCO1490 and HCO2198. Cox1 sequences obtained with the C1-J-2123/C1-N-2776 primer pair can fully sequence both regions (Fig. 5).

Conclusions
This study assessed the usefulness, measured as amplification success, of ten primer pairs routinely employed for targeting the barcoding and other cox1 gene regions in arachnids. Aiming to optimize the efforts in pursuing a longer stretch of cox1 that would maximize the data versus effort ratio for phylogenetic use of the barcode data, we sought an ideal protocol for automatic, reliable and fast extraction of genomic DNA, and developed a new cox1 primer for routine spider systematic work. Our newly designed cox1 primer C1-J-2123 replaces C1-J-2183 to avoid creating an indel region after the Folmer region. This may be especially useful to obtain more complete cox1 data for phylogenetic analyses. We also improved the robotic DNA extraction protocol from the manufacturer's version, adapting it for spider tissue. We are the first to convey usage and set protocol for DNA extraction with MagMAX Express Magnetic Particle Processor in combination with MagMAX DNA Multi-Sample Kit. Our protocol allows for higher DNA concentration output compared with manual DNA extraction using commercial kits. It is thus suitable for semi-high throughput preparation of arachnid DNA. With the 1.5 hour DNA extraction run, the system can be loaded about 8 times per day providing DNA isolated from 192 samples. Following our protocol, PCR amplification of 96 samples is possible in only two hours using a single PCR program. Our protocol is fast and effective and able to provide up to 1000 amplifications per week. Using an even higher throughput system such as the MagMAX Express-96 Magnetic Particle Processor that processes 96 samples at a time, this time could be further cut in half.

Supporting Information
Appendix S1.