Environmental DNA (eDNA) COI metabarcoding Illumina MiSeq NGS PCR Protocol V2
Collin Closek1, Anni Djurhuus2, Katie Pitz3, Ryan Kelly4, Reiko Michisaki3, Kristine Walz3, Hilary Starks1, Francisco Chavez3, Alexandria Boehm1, Mya Breitbart2
1Center for Ocean Solutions, Stanford University, CA; 2University of South Florida, College of Marine Science, St Petersburg, FL; 3Monterey Bay Aquarium Research Institute, Moss Landing, CA; 4University of Washington, Seattle, WA
Mar 20, 2018
ABSTRACT
This protocol is aimed at amplifying the cytochrome c oxidase subunit I (COI) mitochondrial gene in eukaryotes. The primers (forward: mlCOIintF, reverse: HCO2198) utilized in this protocol are based on the primers utilized in Leray et al. 2013 (forward) and Folmer et al. 1994 (reverse).
Amplicons generated using this protocol can then be sequenced using the Illumina platform.
Primers used: Fluidigm CS1+mlCOIinfF Fluidigm CS2+HCO2198
Secondary COI PCR primers PE1-BC-CS1 PE2-BC-CS2
Primer Citations:
PCR
- PCR reactions for COI were run with Fluidigm two-step amplification protocol for each sample
Primary PCR
- Primary PCR amplifications were carried out in single 75-μl reactions using
- 3 μl DNA extract
- 37.5 μl AmpliTaq Gold Fast PCR master mix (Applied Biosystems)
- 3 μl each of forward and reverse primers (5 μM)
- 28.5 μl molecular-biology grade water.
- PCR reactions were run in 96-well plates with a NTC run in singleton for each plate.
- Primary COI cycling parameters - Thse parameters use a normal ramp speed:
95 °C for 10 minutes 16 cycles of the following three steps: - 94 °C for 10 seconds - 62 °C for 30 seconds (this changes -1°C for each subsequent cycle) - 68 °C for 60 seconds Then 25 cycles of the following three steps: - 94 °C for 10 seconds - 46 °C for 30 seconds - 68 °C for 60 seconds A final elongation step of 72 °C for 10 minutes Hold at 4 °C
- COI Primary PCR primers (primers listed in 5’ to 3’ direction)
- Fluidigm CS1+mlCOIinfF (forward):
ACACTGACGACATGGTTCTACAGGWACWGGWTGAACWGTWTAYCCYCC
- Fluidigm CS2+HCO2198 (reverse):
TACGGTAGCAGAGACTTGGTCTTAAACTTCAGGGTGACCAAAAAATCA
- Fluidigm CS1+mlCOIinfF (forward):
Primary PCR clean-up (Bead Clean-up Protocol)
- After primary PCR amplification of the marker region, the PCR products were run through an agarose gel to confirm the presence of target bands and absense of non-specific amplification across environmental samples as well as the absence of amplification in no-template controls (NTCs).
- Primary PCR products were purified and size selected using the Agencourt AMPure XP bead system (Beckman Coulter, USA).
- A second agarose gel was run to confirm primer removal and retention of target amplicons after purication.
The following steps are performed by MSU’s RTSF Genomics Core
Secondary PCR
- An aliquot of 20 μl from each purified primary PCR product was sent to RTSF Genomics Core at MSU for secondary PCR amplification with primers which targeted the CS1/CS2 ends of the primary PCR products and added dual indexed, Illumina compatible adapters with barcodes.
- Secondary PCR amplifications were carried out as single 15-μl reactions using:
- 1 μl template of primary PCR product (no dilution)
- 6 μl OneTaq Hot Start 2X master mix with standard buffer (NEB)
- 1 μl of forward and reverse primer mix (6 μM)
- 7 μl molecular-biology grade water
- Secondary COI cycling parameters:
Secondary COI cycling parameters: 95 °C for 3 minutes 15 cycles of the following three steps: - 95 °C for 15 seconds - 60 °C for 30 seconds - 72 °C for 60 seconds Then a final elongation step of 72 °C for 3 minutes Hold at 25 °C
- Secondary Fluidigm PCR primers (primers listed in 5’ to 3’ direction)
- PE1-BC-CS1 (forward):
AATGATACGGCGACCACCGAGATCT-[i5-BC(index 2)]-ACACTGACGACATGGTTCTACA
- PE2-BC-CS2 (reverse):
CAAGCAGAAGACGGCATACGAGAT-[i7-BC(index 1)]-TACGGTAGCAGAGACTTGGTCT
- PE1-BC-CS1 (forward):
Quality control, PCR clean-up and sequencing parameters
- An agarose gel was run after secondary PCR to confirm the presence of target bands and absense of non-specific amplification across environmental samples as well as the absence of amplification in no-template controls (NTCs).
- After secondary PCR, products were run through Invitrogen SequalPrep Normalization Plate (ThermoFisher Scientific) using manufacturer’s protocol to create pooled library.
Sequencing
- The pooled product for the genetic locus was loaded on a standard MiSeq v2 flow cell and sequenced in a 2x250bp paired end format using a v2 500-cycle MiSeq reagent cartridge.
- The MiSeq run was performed with a 10% PhiX spike added.
- Primers complementary to the Fluidigm CS1 & CS2 oligomers were added to appropriate wells of the reagent cartridge to server as sequencing and index read primers.
- Base calling was done by Illumina Real Time Analysis (RTA) v1.18.54 and output of RTA was demultiplexed and converted to FastQ format with Illumina Bcl2fastq v2.18.0
- COI Sequencing primers (5’ to 3’ direction):
- FL1-CS1(read1)
A+CA+CTG+ACGACATGGTTCTACA
- FL1-CS2(read2)
T+AC+GGT+AGCAGAGACTTGGTCT
- FL2-CS1rc
T+GT+AG+AACCATGTCGTCAGTGT
- FL2-CS2rc(index)
A+GAC+CA+AGTCTCTGCTACCGTA
- FL1-CS1(read1)
- Sequencing is performed at the Research Technology Support Facility (RTSF) Genomics Core at Michigan State University (MSU).