Skip to main content

A universal DNA mini-barcode for biodiversity analysis

Abstract

Background

The goal of DNA barcoding is to develop a species-specific sequence library for all eukaryotes. A 650 bp fragment of the cytochrome c oxidase 1 (CO1) gene has been used successfully for species-level identification in several animal groups. It may be difficult in practice, however, to retrieve a 650 bp fragment from archival specimens, (because of DNA degradation) or from environmental samples (where universal primers are needed).

Results

We used a bioinformatics analysis using all CO1 barcode sequences from GenBank and calculated the probability of having species-specific barcodes for varied size fragments. This analysis established the potential of much smaller fragments, mini-barcodes, for identifying unknown specimens. We then developed a universal primer set for the amplification of mini-barcodes. We further successfully tested the utility of this primer set on a comprehensive set of taxa from all major eukaryotic groups as well as archival specimens.

Conclusion

In this study we address the important issue of minimum amount of sequence information required for identifying species in DNA barcoding. We establish a novel approach based on a much shorter barcode sequence and demonstrate its effectiveness in archival specimens. This approach will significantly broaden the application of DNA barcoding in biodiversity studies.

Background

DNA barcoding seeks to develop a comprehensive species-specific sequence library for all eukaryotes [1]. The 650 bp mitochondrial cytochrome c oxidase 1 (CO1, cox1) DNA barcode [2] is easily sequenced and provides greater than 97% species-level specificity for birds [3], mammals [4], fishes [5], and various arthropods [6]. However, conventional DNA barcoding encounters two problems. First, DNA degradation in archival specimens and processed biological material (i.e. food products) often prevents the recovery of PCR fragments longer than 200 bp, impeding barcode recovery [79]. Second, current approaches cannot be used for comprehensive analysis of environmental samples because high sequence variability necessitates the use of distinct primer sets for each major taxonomic group. In this study, we propose the use of a "mini-barcode" sequence to overcome these problems. We begin by identifying the minimum amount of sequence information required for accurate species identification. We then test the gain in amplification success for smaller fragments in specimens with degraded DNA. Finally, by targeting conserved priming sites within the barcode region we develop primers with the universality required for the analysis of all major eukaryotes.

Results and Discussion

To determine how much sequence information is required for identifications, we retrieved all CO1 barcode sequences from GenBank and calculated the probability of having species-specific barcodes for varied size fragments (Figure 1). Our analysis shows that while full-length DNA barcodes perform best (97% species resolution), 90% identification success is obtained with 100 bp regions and 95% success with 250 bp barcodes. In another words, in 90% of the species tested a DNA barcode of only 100 bp contains nucleotide substitution(s) specific to members of a particular species. Having established the potential of mini-barcodes to identify unknown specimens, we proceeded to design primers and test their performance.

Figure 1
figure 1

Short DNA barcode sequences (less than 150 bp) provide efficient taxonomic sequence tags. Increase in species resolution with increasing sequence length of cytochrome c oxidase 1.

We designed universal primers for amplifying mini-barcodes by aligning CO1 sequences from all major eukaryote groups (animals, fungi, plants, and protists) and identifying conserved amino acid strings for primers in the size range 120–150 bp. We selected a single primer set for a 130 bp amplicon (see below). We tested the ability of this universal primer set to amplify DNA extracts from 1,566 specimens derived from 691 species of mammals, fishes, birds, and insects (mainly Lepidoptera, Hymenoptera, Ephemeroptera, Plecoptera, and Trichoptera). Moreover, we examined another 330 DNA extracts from plants, fungi, and macroalgae (Additional file 1). We obtained PCR amplicons from 92% of these species (Additional file 2) and compared this success to that obtained via routine barcoding. We were able to amplify mini-barcodes with greater success than full-length sequences in all groups but one (9% lower success in Plecoptera) (Figure 2). We verified that the 130 bp amplicons were CO1 by sequencing approximately half of them (results not shown).

Figure 2
figure 2

Amplification success for a 130 bp amplicon of CO1 versus the standard 650 bp. The mini-barcodes were amplified with a single primer set whereas the full-length DNA barcode was amplified using different, taxon-specific primers.

One important application of mini-barcodes lies in obtaining sequence information from old type specimens. We tested this approach on a collection of Coleophora (Order Lepidoptera) which are difficult to identify because of their small size and cryptic morphology. We successfully sequenced all 15 dried museum specimens collected from 1871 to 1944 (including 7 type specimens) in one PCR pass with the universal mini-barcode primers and compared these specimens to 84 full-length barcodes from recently collected specimens of this genus (Figure 3). Comparison of sequence information obtained from these old specimens with recently collected samples provided excellent corroboration and species-level resolution. Mini-barcode sequences from old museum specimens formed monophyletic groups, containing zero or very low sequence divergence, with freshly collected specimens of the corresponding species (Figure 3).

Figure 3
figure 3

The 130 bp amplicon is reliably amplified from old museum specimens. Example here is from the species of the Lepidoptera genus Coleophora. The tree is reconstructed using neighbor-joining method [12] with Kimura two-parameter evolutionary distances [13]. Museum samples are shown by blue and red (for type specimens) dots. Collection dates are shown in parentheses.

Conclusion

The mini-barcode system dramatically broadens the applications of DNA barcoding. We have now demonstrated that sequence information can be reliably obtained from archival specimens or those with degraded DNA. Further, the universality of the primers enables the recovery of comprehensive barcode information from environmental mixtures. Finally, the short universally-primed amplicon is ideal for sequence characterization through new parallelized high-throughput sequencing technologies, allowing inexpensive but comprehensive studies of biodiversity to be a realistic goal.

Methods

Bioinformatics analysis

Metazoan COI sequences bearing the "BARCODE" keyword were downloaded from GenBank using the NCBI eFetch tool. Barcodes that were less than 650 bases in length were eliminated, leaving a dataset of 6,695 barcode sequences from 1,587 species. For various sizes of 5'-end minibarcodes, ranging in size from 10 bases up to the full-length of the barcode sequence, we analyzed the number of species that could be uniquely identified (to the exclusion of other species) using that sequence.

Specimens and their taxonomic coverage

All DNA extracts were obtained from different barcoding projects in the Canadian Centre for DNA Barcoding and external collaborators. We selected these samples considering maximum taxonomic range.

Primer design strategy

We selected the 5' end of the barcode region targeting a 100–150 base amplicon. By comparing a wide range of taxa in this region, we found well-conserved strings of amino acids across all taxa in priming sites. Interestingly, this high level of conservation is also evident at nucleotide level. We designed multiple oligos by using the Primer3 program [10] and considering physical and structural properties of oligos (such as annealing temperature, G+C percentage, and self-complementarity). We selected the primer Uni-MinibarR1: 5'-GAAAATCATAATGAAGGCATGAGC-3' for further testing as it represented the highest similarity – especially at the 3' end – to other taxa. A similar strategy was used for designing a forward primer: Uni-MinibarF1: 5'-TCCACTAATCACAARGATATTGGTAC-3'. This primer is positioned in the same region as other common barcoding primers are located. We attached M13 forward and reverse tails to our forward and reverse primers, respectively, to facilitate the high throughput sequencing process. These tails did not decrease the PCR success.

PCR Optimization Strategy

PCR reactions were performed using a standard PCR pre-mix [11]. We used the above mentioned universal primer set in all the reactions in a touch up PCR program: 95°C for 2 min, followed by 5 cycles of 95°C-1 min, 46°C-1 min, and 72°C-30 sec, followed by 35 cycles of 95°C-1 min, 53°C-1 min, and 72°C-30 sec, and finally a final extension at 72°C for 5 min. We used a Mastercycler ep gradient S (Eppendorf, Mississauga, ON, Canada) thermalcycler. We included two negative control reactions (no DNA template) in all our PCR 96-well plates. To compare the universal mini-barcode primer set with the specific full-length primers we amplified DNA extracts using taxonomically specific primer sets (i.e. 2 primer sets for fish species) [11].

PCR amplification verification and sequencing

PCR products were visualized on a 2% E-gel® 96 Agarose (Invitrogen, Burlington, ON, Canada). The bands on E-gel were used as a measure of PCR success. To verify the amplification of the target region, we sequenced 747 PCR products from at least 363 species. Standard BigDye kits (Applied Biosystems, Foster City, CA) were used for sequencing. Sequencing reactions were cleaned up by using Agencourt® CleanSEQ® kit (Agencourt Bioscience Corporation, Beverly, MA). The sequences were run bidirectionally on a 3730xl DNA analyzer (Applied Biosystems, Foster City, CA, USA), edited with Sequencher™ (Gene Codes Corporation, Ann Arbor, MI), and aligned using BioEdit version 7.0.5.3.

References

  1. Marshall E: Taxonomy. Will DNA bar codes breathe life into classification?. Science. 2005, 307 (5712): 1037-10.1126/science.307.5712.1037.

    Article  PubMed  CAS  Google Scholar 

  2. Hebert PD, Cywinska A, Ball SL, deWaard JR: Biological identifications through DNA barcodes. Proc Biol Sci. 2003, 270 (1512): 313-321. 10.1098/rspb.2002.2218.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  3. Hebert PD, Stoeckle MY, Zemlak TS, Francis CM: Identification of birds through DNA barcodes. PLoS Biology. 2004, 2 (10): E312-10.1371/journal.pbio.0020312.

    Article  PubMed  PubMed Central  Google Scholar 

  4. Hajibabaei M, Singer GA, Clare EL, Hebert PDN: Design and applicability of DNA arrays and DNA barcodes in biodiversity monitoring. BMC Biol. 2007, 5: 24-10.1186/1741-7007-5-24.

    Article  PubMed  PubMed Central  Google Scholar 

  5. Ward RD, Zemlak TS, Innes BH, Last PR, Hebert PD: DNA barcoding Australia's fish species. Philos Trans R Soc Lond B Biol Sci. 2005, 360 (1462): 1847-1857. 10.1098/rstb.2005.1716.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  6. Hajibabaei M, Janzen DH, Burns JM, Hallwachs W, Hebert PDN: DNA barcodes distinguish species of tropical Lepidoptera. Proceedings of the National Academy of Sciences of the United States of America. 2006, 103 (4): 968-971. 10.1073/pnas.0510466103.

    Article  PubMed  PubMed Central  Google Scholar 

  7. Goldstein PZ, Desalle R: Calibrating phylogenetic species formation in a threatened insect using DNA from historical specimens. Molecular Ecology. 2003, 12 (7): 1993-1998. 10.1046/j.1365-294X.2003.01860.x.

    Article  PubMed  CAS  Google Scholar 

  8. Hajibabaei M, Smith MA, Janzen DH, Rodriguez JJ, Whitfield JB, Hebert PDN: A minimalist barcode can identify a specimen whose DNA is degraded. Molecular Ecology Notes. 2006, 6: 959-964. 10.1111/j.1471-8286.2006.01470.x.

    Article  CAS  Google Scholar 

  9. Wandeler P, Hoeck PE, Keller LF: Back to the future: museum specimens in population genetics. Trends Ecol Evol. 2007, 22 (12): 634-642. 10.1016/j.tree.2007.08.017.

    Article  PubMed  Google Scholar 

  10. Rozen S, Skaletsky H: Primer3 on the WWW for general users and for biologist programmers. Methods Mol Biol. 2000, 132: 365-386.

    PubMed  CAS  Google Scholar 

  11. The Canadian Centre for DNA Barcoding. [http://www.dnabarcoding.ca/pa/ge/research/protocols/amplification]

  12. Saitou N, Nei M: The neighbor-joining method: a new method for reconstructing phylogenetic trees. Mol Biol Evol. 1987, 4 (4): 406-425.

    PubMed  CAS  Google Scholar 

  13. Kimura M: A simple method for estimating evolutionary rates of base substitutions through comparative studies of nucleotide sequences. Journal of molecular evolution. 1980, 16 (2): 111-120. 10.1007/BF01731581.

    Article  PubMed  CAS  Google Scholar 

Download references

Acknowledgements

We thank Elizabeth Clare, Robin Floyd, Kevin Kerr, Natalia Ivanova, Gary Saunders, Alex Smith, Magali Sole, Dirk Steinke, and Xin Zhou for providing DNA extracts. This work is supported by grants from Genome Canada (through the Ontario Genomics Institute) and Natural Sciences and Engineering Research Council of Canada to the Canadian Barcode of Life Network.

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Mehrdad Hajibabaei.

Additional information

Authors' contributions

IM performed molecular methods, assisted in project design and data analysis, and edited the manuscript. GACS gathered sequence information from GenBank, designed and conducted bioinformatics analysis, and edited the manuscript. J–FL carried out taxonomic analysis and edited the manuscript. DAH aided the study design, provided tools/reagents, and edited the manuscript. PDNH aided the study design, provided tools/reagents, and edited the manuscript. MH designed the project, performed DNA barcode data analysis, and wrote the manuscript. All authors read and approved the final manuscript.

Electronic supplementary material

Additional file 1: Specimens used in this study. (PDF 86 KB)

12864_2007_1407_MOESM2_ESM.pdf

Additional file 2: Mini-barcodes of 130 bp can be PCR amplified from the majority of specimens using a single universal primer set for all major eukaryotic lineages. (PDF 606 KB)

Authors’ original submitted files for images

Below are the links to the authors’ original submitted files for images.

Authors’ original file for figure 1

Authors’ original file for figure 2

Authors’ original file for figure 3

Rights and permissions

Open Access This article is published under license to BioMed Central Ltd. This is an Open Access article is distributed under the terms of the Creative Commons Attribution License ( https://creativecommons.org/licenses/by/2.0 ), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Reprints and permissions

About this article

Cite this article

Meusnier, I., Singer, G.A., Landry, JF. et al. A universal DNA mini-barcode for biodiversity analysis. BMC Genomics 9, 214 (2008). https://doi.org/10.1186/1471-2164-9-214

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/1471-2164-9-214

Keywords