- Methodology article
- Open Access
Enhanced primers for amplification of DNA barcodes from a broad range of marine metazoans
BMC Ecology volume 13, Article number: 34 (2013)
Building reference libraries of DNA barcodes is relatively straightforward when specifically designed primers are available to amplify the COI-5P region from a relatively narrow taxonomic group (e.g. single class or single order). DNA barcoding marine communities have been comparatively harder to accomplish due to the broad taxonomic diversity and lack of consistently efficient primers. Although some of the so-called “universal” primers have been relatively successful, they still fail to amplify COI-5P of many marine animal groups, while displaying random success even among species within each group. Here we propose a new pair of primers designed to enhance amplification of the COI-5P region in a wide range of marine organisms.
Amplification tests conducted on a wide range of marine animal taxa, rendered possible the first–time sequencing of DNA barcodes from eight separated phyla (Annelida, Arthropoda, Chordata, Cnidaria, Echinodermata, Mollusca, Nemertea and Platyhelminthes), comprising a total of 14 classes, 28 orders, 57 families, 68 genus and 76 species.
These primers demonstrated to be highly cost-effective, which is of key importance for DNA barcoding procedures, such as for building comprehensive DNA barcode libraries of marine communities, where the processing of a large numbers of specimens from a wide variety of marine taxa is compulsory.
The so-called “universal” PCR primers are research tools of great utility for molecular ecological studies where the same locus is analysed across a broad range of taxa from different phyla. The universal primers designed by Folmer and colleagues (HCO2198-LCO1490, henceforth named “Folmer primers”) for amplification of a 658 base pair (bp) fragment of the 5′ end of the mitochondrial gene cytochrome c oxidase subunit I (COI-5P), have shown to be very successful in the amplification of this gene fragment in a broad range of marine metazoan phyla. The Folmer primers have been probably the most widely used primer pair for amplification of COI in many animal groups in addition to marine organisms (3166 citations recorded in ISI, 25/10/2012) (e.g.[2, 3]). Indeed, this was the primer pair selected by Hebert and colleagues for their proof-of concept study on Canadian moths, where they propose the DNA barcoding approach for species identification. The mitochondrial DNA region delimited by Folmer’s primers (COI-5P) became the established DNA barcode region for animal life. With a growing number of studies attempting to examine DNA barcodes from different animal taxa, it quickly became apparent that the primer pair HCO2198 and LCO1490 was not so “universal” as thought before, as it would still fail to amplify some taxa.
Simison, for instance, verified that the Folmer primers and several combinations with degenerate primers were successful for only a small number of gastropod taxa. Also, Lohman et al. tried to amplify COI-5P of perching birds (Aves: Passeriformes) with the universal LCO1490/HCO2198 primer pair, but the efforts were similarly unsuccessful, as they were with different combinations of primers. In another example, Blankenship et al., failed to amplify COI-5P using the Folmer primers from the remnants of big-eye tuna (identified only through the mitochondrial 16S gene) inside the guts of deep-sea amphipods (Scopelocheirus schellenbergi). Also, the usefulness of the Folmer primers may be limited for decapod Crustacea because they are not optimized. In echinoderms, amplification of COI is often challenging, either due to low amplification successor to the amplification of pseudogenes[9, 10]. Given the vast geographic distribution and the important commercial value of holothurians, further work should focus on developing alternative primers for these species, because PCR (polymerase chain reaction) amplification could not be achieved with the primers available.
To overcome the limitations of the presumed universal primers, new primer pairs have been developed targeting specific large assemblages, such as birds, lepidopterans, or fish. Still, the limited amplification success in some groups led to the development of alternative approaches, namely the design of degenerate primer pairs or even primer cocktails. The later approach has been very successful in COI-5P amplification in fish (e.g.), although one of the primer cocktails was originally designed and tested in mammals, and the same occurred for the alternative primer pairs designed for birds or Lepidoptera. Despite the success of these group-specific PCR primers, to date no alternatives to Folmer primers have been proposed that are effective in a broad range of marine animals and particularly for marine invertebrates. Within the latter group, most PCR primers developed were phylum or class specific (see), like primers designed for Echinoderms, Crustacea (, D. Steinke unpublished in), Gastropoda and Annelids. Here we propose a new pair of enhanced primers specifically designed to amplify the COI-5P barcode region from a broad taxonomic range of marine organisms. We compared its amplification potential with other broadly used primer pairs, and tested amplification success in 76 species from both vertebrates and invertebrates and a total of 8 animal phyla. The success of amplification in a broad-range of taxa indicates that these primers can be particularly valuable for building up reference DNA barcode libraries of complex marine communities.
Marine specimens were collected from various locations along the west coast of Portugal and preserved in 96% ethanol. A total of 130 specimens and 76 species belonging to 8 among the most common phyla of marine metazoa was identified by traditional methods based on morphological characters. The list of specimens investigated and respective collection data is displayed in Table 1.
Muscle tissue was used to extract DNA from the specimens of all surveyed groups. DNA extracts were obtained using the E.Z.N.A. Mollusc DNA Kit (Omega Bio-tek), following the manufacturer’s instructions.
Degenerate primers LoboF1 (forward) and LoboR1 (reverse) for the COI-5P were designed based on publically available COI sequences obtained from GenBank that were compatible with the COI–5P sequence that included complete 5′ and 3′ ends. They were analyzed and aligned using MEGA software version 5. To assist the design, three pairs of available primers were added to the alignment: the forward primers LCO1490, CrustDF1 (D. Steinke unpublished in), CrusF1 and CrusF2 and the reverse primers HCO2198 and CrustDR1 (D. Steinke unpublished in) (Figure 1).
The amplifications were carried out in an iCycler™ (Bio-Rad) thermal cycler, using the new primers LoboF1 and LoboR1, and a pre-made PCR mix from Invitrogen™. The mix contained 1× PCR buffer, 1.5 mM of MgCl2, 0.2 mM of the dNTP mixture, 1 U of DNA Taq polymerase plus 0.5 μM of each primer and 4 μL of DNA template and completed with sterile milliQ-grade water to make up a total volume of 25 μL.
DNA barcode amplification
The optimal annealing temperatures of primers were tested before setting up the PCR thermal cycling conditions based on the commonly used protocol in DNA barcoding studies: 1) denaturation at 94°C for 1 min; 2) denaturation at 94°C (30 s), annealing at 45°C (1 min 30 s), extension as 72°C for 1 min (5 cycles); 3) denaturation at 94°C (30 s), annealing at 54°C (1 min 30 s), extension as 72°C for 1 min (45 cycles) and ended with a final extension of 5 min at 72°C. The PCR products together with a 100 bp DNA ladder (Invitrogen™) were separated by electrophoresis in 1.5% agarose gel in TAE buffer, and subsequently stained with ethidium bromide for visualization in GelDoc 2000 equipment (Bio-Rad™). PCR products were cleaned up using a three-time precipitation with isopropanol, sequenced bidirectionally using the BigDye Terminator 3 kit, and run on an ABI 3730XL DNA analyser (all from Applied Biosystems™).
An initial test was conducted to compare the potential of newly-designed primers to amplify the COI barcode fragment, with two of the most common primer pairs used for marine organisms. The primers compared were LCO1490 and HCO2198, CrustDF1 and CrustDR1 and the new LoboF1 and LoboR1. DNA templates from six specimens of three crustacean species were tested in parallel for the 3 primer pairs employing the same PCR conditions.
Sequence alignment and tree reconstruction
All sequence data were carefully checked to detect possible nuclear mitochondrial pseudogenes (numts) (see): chromatogram quality examination and sequence editing to detect ambiguities, double peaks and noise, GenBank’s BLAST search for homologies with public sequences (Table 1), and translation into amino acid for inspection for any indels, stop codons, or unusual amino acid sequence patterns. The GenBank accession numbers for sequences obtained in this study are included between KF369103 and KF369196, and specimen and sequence data is compiled in the Barcode of Life Data Systems (BOLD) project titled “Enhanced primers for amplification of DNA barcodes from marine metazoans”.
All sequence alignments and tree reconstruction were performed using MEGA5 software. The neighbour-joining (NJ) method was used, applying the Kimura-2-parameter (K2P) model for the nucleotide–base tree and the Jones-Taylor-Thornton (JTT) model for the amino acid sequence tree. For assessment of node support 1000 bootstrap iterations were run in both cases.
The primers LoboF1 and LoboR1 amplified the COI-5P (Figure 2) from a wide diversity of species from different animal phyla (Table 1), both vertebrates and invertebrates, comprising 130 specimens belonging to 76 species from 8 phyla. In addition, one brown alga species (Fucus spiralis) was successfully amplified (Table 1). All PCR products were successfully sequenced (forward and reverse sequences), with 87% of the cases returning the expected ≈ 658 bp length after sequence adition and primer trimming. No indication of the presence of numts was detected in any of the sequences, although the bivalve Anomia ephippium displayed two stop codons (TAA) on positions 3 and 158 of the amino acid sequence. This particular exception is further analysed in the Discussion section.
Contrasting the resulting sequences with those deposited at GenBank and BoldSystems databases confirmed the accordance between the morphological and genetic classification of most species. Although some specimens referred to species not yet accessioned in the databases, the specimen similarity obtained was always within the same taxonomic group, albeit at a higher taxonomic level.
Comparing the primers used in the initial test performed on six specimens, we verified that LoboF1 and LoboR1 were the only primers that provide DNA amplifications for all specimens (Figure 3). Folmer primers failed to amplify COI-5P for 4 out of 6 decapod samples and the crustacean-specific primers CrustDF1 and CrustDR1 amplified five of the six samples.
The newly-designed primers permitted amplification of COI-5P from all marine species tested, which comprised a wide phylogenetic diversity of taxa belonging to eight different phyla (Annelida, Arthropoda, Chordata, Cnidaria, Echinodermata, Mollusca, Nemertea and Platyhelminthes), and including some of the organisms considered more recalcitrant to amplification such as echinoderms and flatworms. To our best knowledge, this is the first DNA barcoding study reporting both successful amplification and sequencing from such a diverse array of taxa using a single pair of primers. Although the primers here proposed are degenerate, we did not find any indication of the presence of pseudogenes after detailed inspection of the sequences. The sequence chromatograms were generally of high quality, and neither indels or stop codons were observed, nor unusual substitution patterns in more conserved regions of the nucleotide or amino acid sequences. BLAST searches returned the expected matches or nearest neighbours, and we also did not find any evidence of accidental amplification of microbial contaminants (see).
The only exception to this pattern was the presence of two stop codons (TAA) in the bivalve Anomia ephippium. Plazzi et al. previously reported the occurrence of exactly the same stop codon (on position 158 of the amino acid sequence) on a specimen of the same genus (Anomia sp.), but still considered their sequence to have genuine mtDNA origin. Indeed, the chromatograms from our three specimens of A. ephippium had very high quality scores and did not show any other abnormal substitution except for these exact codons. Such a precise codon-specific substitution, occurring in three organisms of the same genus (diverging by 30%), should not be compatible with the expected random substitution pattern typical of a pseudogene (see for instance). Similarly to Plazzi et al., we suggest that this is a bona fide mtDNA COI sequence, and that the putative stop codon could be an exception to the mitochondrial code of this taxon that deserves further investigation.
Among the phyla that yielded successful COI sequences, echinoderms and flatworms (Platyhelminthes) have long been receiving particular attention in DNA barcoding and related studies due to acknowledged difficulties in obtaining quality sequences. Three species of three different classes of echinoderms had COI-5P successfully amplified and sequenced, albeit other authors acknowledge that amplifying COI-5P from echinoderms remains a challenge, probably due to the unsuitability of the primers used or to DNA degradation. On the other hand, according to León-Règagnon et al., the generation of primers that amplify the standard barcode region for a wide variety of flatworm groups is still at the experimental stage and the COI-5P traditionally used in platyhelminths is shorter than the standard barcodes. In the present study, only one species of platyhelminth (Leptoplana tremellaris) was attempted for amplification and sequencing, however, the standard COI-5P was successfully obtained, which may reveal a promising trend for the primers’ application to other flatworms. Additionally, even though the primer design was primarily aimed at invertebrate sequences, five distinct species of telosts (Dicentrarchus sp., Perciformes: Moronidae, Lepadogaster lepadogaster, Gobiesociformes: Gobiesocidae, Pomatoschistus lozanoi, Perciformes: Gobiidae, Solea senegalensis, Pleuronectiformes: Soleidae and Thunnus sp., Perciformes: Scombridae) were amplified and successfully sequenced. It should be noted, at this point, that several primer pairs or even multiple primer cocktails are usually used to amplify COI-5P. In spite of these promising preliminary results, research has still to be performed to confirm the widespread suitability of the LoboF1/R1 for teleosts.
Currently, COI-5P is a very useful marker to identify species and is generally regarded as a powerful tool for molecular taxonomy. Furthermore, thanks to the BOLD multiple functionalities, sequences of COI-5P can already be contrasted with those introduced by other researchers. While COI-5P gives high support values for species identification, it does not resolve well deep molecular diversity, and shows a low level of phylogenetic informative characters, but the combination of nuclear and mitochondrial genes is useful for phylogenetic relationship (see[33, 34]). Additionally, it must be stressed that DNA barcoding has not been technically conceived to recover phylogenetic relationships but rather to identify known species and to aid the discovery of new ones, which contributes to explain the misallocation of higher taxa hereby observed, in accordance to the observations by Hajibabaei et al.. Still, regardless of technical limitations, all specimens of the same species were grouped in the same clade when either amino acid– or nucleotide–based trees were reconstructed, as expected.
The newly designed primers LoboF1 and LoboR1 proved to be a rapid, practical and cost-effective tool for DNA barcoding based on COI sequencing. Given their high success rate to amplifying COI-5P from a wide and diverse range of tested taxa, it may well be that these new primers have a broader taxonomic range than all those currently used for barcoding the kingdom Animalia. These primers thus have high potential for the accelerated build–up of a global DNA–based biodiversity library, particularly with regard to the marine component.
Folmer O, Black M, Hoeh W, Lutz R, Vrijenhoek R: DNA primers for amplification of mitochondrial cytochrome c oxidase subunit I from diverse metazoan invertebrates. Mol Mar Biol Biotechnol. 1994, 3: 294-299.
Dawnay N, Ogden R, McEwing R, Carvalho GR, Thorpe RS: Validation of the barcoding gene CO1 for use in forensic genetic species identification. Forensic Sci Int. 2007, 173: 1-6. 10.1016/j.forsciint.2006.09.013.
Ekrem T, Willassen E, Stur E: A comprehensive DNA sequence library is essential for identification with DNA barcodes. Mol Phylogenet Evol. 2007, 43: 530-542. 10.1016/j.ympev.2006.11.021.
Hebert PDN, Cywinska A, Ball SL, deWaard JR: Biological identifications through DNA barcodes. Proc R Soc Lond B Biol Sci. 2003, 270: 313-321. 10.1098/rspb.2002.2218.
Simison WB: Evolution and phylogeography of new world gastropod faunas. 2000, Berkeley: University of California
Lohman DJ, Prawiradilaga DM, Meier R: Improved COI barcoding primers for Southeast Asian perching birds (Aves: Passeriformes). Mol Ecol Resour. 2009, 9: 37-40. 10.1111/j.1755-0998.2008.02221.x.
Blankenship LE, Yayanos AA: Universal primers and PCR of gut contents to study marine invertebrate diets. Mol Ecol. 2005, 14: 891-899. 10.1111/j.1365-294X.2005.02448.x.
Schubart CD: Mitochondrial DNA and decapod phylogenies: the importance of pseudogenes and primer optimization. Decapod Crustacean Phylogenetics. Edited by: Martin JW, Crandall KA, Felder DL. 2009, Boca Raton, FL: CRC Press, Taylor & Francis Group, 47-65.
Hoareau TB, Boissin E: Design of phylum-specific hybrid primers for DNA barcoding: addressing the need for efficient COI amplification in the Echinodermata. Mol Ecol Resour. 2010, 10: 960-967. 10.1111/j.1755-0998.2010.02848.x.
Ward RD, Holmes BH, O’Hara TD: DNA barcoding discriminates echinoderm species. Mol Ecol Resour. 2008, 8: 1202-1211. 10.1111/j.1755-0998.2008.02332.x.
Uthicke S, Byrne M, Conand C: Genetic barcoding of commercial Bêche-de-mer species (Echinodermata: Holothuroidea). Mol Ecol Resour. 2010, 10: 634-646.
Hebert PDN, Stoeckle MY, Zemlak TS, Francis CM: Identification of birds through DNA barcodes. PLoS Biol. 2004, 2: e312-10.1371/journal.pbio.0020312.
Hajibabaei M, Janzen DH, Burns JM, Hallwachs W, Hebert PDN: DNA barcodes distinguish species of tropical Lepidoptera. Proc Natl Acad Sci USA. 2006, 103: 968-971. 10.1073/pnas.0510466103.
Ward RD, Zemlak TS, Innes BH, Last PR, Hebert PDN: Barcoding Australia’s fish species. Philos Trans R Soc Lond B Biol Sci. 2005, 360: 1847-1857. 10.1098/rstb.2005.1716.
Ivanova NV, Zemlak TS, Hanner RH, Hebert PDN: Universal primer cocktails for fish DNA barcoding. Mol Ecol Notes. 2007, 7: 544-548. 10.1111/j.1471-8286.2007.01748.x.
Costa FO, Landi M, Martins R, Costa MH, Costa ME, Carneiro M, Alves MJ, Steinke D, Carvalho GR: A ranking system for reference libraries of DNA barcodes: application to marine fish species from Portugal. PLoS ONE. 2012, 7: e35858-10.1371/journal.pone.0035858.
Bucklin A, Steinke D, Blanco-Bercial L: DNA barcoding of marine metazoa. Ann Rev Mar Sci. 2011, 3: 471-508. 10.1146/annurev-marine-120308-080950.
Costa FO, de Waard JR, Boutillier J, Ratnasingham S, Dooh RT, Hajibabaei M, Hebert PDN: Biological identifications through DNA barcodes: the case of the Crustacea. Can J Fish Aquat Sci. 2007, 64: 272-295. 10.1139/f07-008.
Matzen da Silva J, Creer S, Dos Santos A, Costa AC, Cunha MR, Costa FO, Carvalho GR: Systematic and evolutionary insights derived from mtDNA COI barcode diversity in the decapoda (crustacea: malacostraca). PLoS ONE. 2011, 6: e19449-10.1371/journal.pone.0019449.
Zou S, Li Q, Kong L: Multigene Barcoding and phylogeny of geographically widespread muricids (gastropoda: neogastropoda) along the coast of china. Mar Biotechnol (NY). 2010, 14: 21-34.
Carr CM, Hardy SM, Brown TM, Macdonald TA, Hebert PDN: A Tri-Oceanic perspective: DNA barcoding reveals geographic structure and cryptic diversity in Canadian polychaetes. PLoS ONE. 2011, 6: e22232-10.1371/journal.pone.0022232.
Tamura K, Peterson D, Peterson N, Stecher G, Nei M, Kumar S: MEGA5: molecular evolutionary genetics analysis using maximum likelihood. Evolutionary distance and maximum parsimony methods. Mol Biol Evol. 2011, 28: 2731-2739. 10.1093/molbev/msr121.
Song H, Buhay JE, Whiting MF, Crandall KA: Many species in one: DNA barcoding overestimates the number of species when nuclear mitochondrial pseudogenes are co amplified. Proc Natl Acad Sci USA. 2008, 105: 13486-13491. 10.1073/pnas.0803076105.
McGinnis S, Madden TL: BLAST: at the core of a powerful and diverse set of sequence analysis tools. Nucleic Acids Res. 2004, 32: W20-W25. 10.1093/nar/gkh435.
Kimura M: A simple model for estimating evolutionary rates of base substitutions through comparative studies of nucleotide sequences. J Mol Evol. 1980, 16: 111-120. 10.1007/BF01731581.
Jones DT, Taylor WR, Thornton JM: The rapid generation of mutation data matrices from protein sequences. Comput Appl Biosci. 1992, 8: 275-282.
Saitou N, Nei M: The neighbor-joining method: a new method for reconstructing phylogenetic trees. Mol Biol Evol. 1987, 4: 406-425.
Siddall ME, Fontanella FM, Watson SC, Kvist S, Eresus C: Barcoding bamboozled by bacteria: convergence to Metazoan Mitochon drial primer targets by marine microbes. Syst Biol. 2009, 58: 445-451. 10.1093/sysbio/syp033.
Plazzi F, Passamonti M: Towards a molecular phylogeny of Mollusks: Bivalves’ early evolution as revealed by mitochondrial genes. Mol Phylogenet Evol. 2010, 57: 641-657. 10.1016/j.ympev.2010.08.032.
León-Règagnon V: Evidence of new species of Haematoloechus (Platyhelminthes: Digenea) using partial cox1 sequences. Mitochondrial DNA. 2010, 21: 12-17. 10.3109/19401736.2010.523700.
Ratnasingham S, Hebert PDN: BOLD: The Barcode of Life Data System. Mol Ecol notes. 2007, 7: 355-364. 10.1111/j.1471-8286.2007.01678.x.http://www.barcodinglife.org,
Zarowiecki MZ, Huyse T, Littlewood DTJ: Making the most of mitochondrial genomes – Markers for phylogeny, molecular ecology and barcodes in Schistosoma (Platyhelminthes: Digenea). Int J Parasitol. 2007, 37: 1401-1418. 10.1016/j.ijpara.2007.04.014.
Simon C, Frati F, Beckenbach A, Crespi B, Liu H, Flook P: Evolution, weighting and phylogenetic utility of mitochondrial gene sequences and a compilation of conserved polymerase chain reaction primers. Ann Entomol Soc Am. 1994, 87: 651-701.
Simon C, Buckley TR, Frati F, Stewart JB, Beckenbach A: Incorporating molecular evolution into phylogenetic analysis, and a new compilation of conserved polymerase chain reaction primers for animal mitochondrial DNA. Ann Rev Ecol Evol Syst. 2006, 37: 545-579. 10.1146/annurev.ecolsys.37.091305.110018.
Hebert PDN, Gregory TR: The Promise of DNA Barcoding for Taxonomy. Syst Biol. 2005, 54: 852-859. 10.1080/10635150500354886.
Hajibabaei M, Singer G, Hickey D: Benchmarking DNA barcodes. Genome. 2006, 49: 851-854. 10.1139/G06-025.
This work was supported by FEDER through POFC-COMPETE and by national funds from "Fundação para a Ciência e a Tecnologia (FCT)" in the scope of the grants FCOMP-01-0124-FEDER-015429 and PEst-C/BIA/UI4050/2011. Jorge Lobo is supported by a FCT PhD grant (SFRH/BD/69750/2010). Pedro M Costa is supported by a FCT post-doc grant (SFRH/BPD/72564/2010). Filipe O Costa benefitted from the Marie Curie European Reintegration Grant PERG02-GA–2007-224890, provided by the European Commission. The authors would also like to thank Carolina Madeira and Monica Landi for the initial molecular experiments, Luisa Borges and Ronaldo Sousa for the support on collections, Cátia Bartilotti (IPMA), Nuno Gomes (CBMA) for the offer of crustacean specimens and Artur Oliveira (CBMA) for the alga specimen.
The authors declare that they have no competing interests.
JL participated in sample collections and morphological identifications, performed the DNA extraction, amplification and purification, designed the primers and wrote the manuscript. PMC provided support on molecular components, participated in the design of the primers and edited the manuscript. MALT and MSGF contributed on annelid and amphipod collections, morphological identifications, molecular and bioinformatic analysis. MHC participated in sample collections, morphological identifications and edited the manuscript. FOC provided supervision and training on molecular and bioinformatics components and edited the manuscript. All authors read and approved the final manuscript.
About this article
Cite this article
Lobo, J., Costa, P.M., Teixeira, M.A. et al. Enhanced primers for amplification of DNA barcodes from a broad range of marine metazoans. BMC Ecol 13, 34 (2013) doi:10.1186/1472-6785-13-34
- Universal primers
- DNA barcoding
- Marine metazoans