Figures
Abstract
The objective of this study was to build a DNA barcode reference library for the true butterflies of Peninsula Malaysia and assess the value of attaching subspecies names to DNA barcode records. A new DNA barcode library was constructed with butterflies from the Museum of Zoology, University of Malaya collection. The library was analysed in conjunction with publicly available DNA barcodes from other Asia-Pacific localities to test the ability of the DNA barcodes to discriminate species and subspecies. Analyses confirmed the capacity of the new DNA barcode reference library to distinguish the vast majority of species (92%) and revealed that most subspecies possessed unique DNA barcodes (84%). In some cases conspecific subspecies exhibited genetic distances between their DNA barcodes that are typically seen between species, and these were often taxa that have previously been regarded as full species. Subspecies designations as shorthand for geographically and morphologically differentiated groups provide a useful heuristic for assessing how such groups correlate with clustering patterns of DNA barcodes, especially as the number of DNA barcodes per species in reference libraries increases. Our study demonstrates the value in attaching subspecies names to DNA barcode records as they can reveal a history of taxonomic concepts and expose important units of biodiversity.
Citation: Wilson J-J, Sing K-W, Sofian-Azirun M (2013) Building a DNA Barcode Reference Library for the True Butterflies (Lepidoptera) of Peninsula Malaysia: What about the Subspecies? PLoS ONE 8(11): e79969. https://doi.org/10.1371/journal.pone.0079969
Editor: M. Alex Smith, University of Guelph, Canada
Received: April 4, 2013; Accepted: October 7, 2013; Published: November 25, 2013
Copyright: © 2013 Wilson et al. This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.
Funding: This study was supported by University of Malaya Research Grant RG158/12SUS and also through subsidised analyticals costs at the Canadian Centre for DNA Barcoding under the iBOL program. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.
Competing interests: The authors have declared that no competing interests exist.
Introduction
Surveys of butterfly species have often been considered good surrogates for surveys of total biodiversity (e.g. in Malaysia [1]). This is because of their role in food webs - caterpillars consume large quantities of plants and are themselves consumed by other animals in large numbers - and because, relative to most other animal groups, collecting and identifying adult butterflies is considered easy [1]. This is particularly so in Peninsula Malaysia where butterflies have received intensive taxonomic attention. The “Butterflies of the Malay Peninsula” have been the subject of a series of comprehensive field guides, beginning with Distant in 1882–1886 [2], and followed by four editions of Corbet and Pendlebury’s classic checklist, first published in 1934 [3] and most recently revised by Eliot in 1992 [4]. Butterflies have benefitted and suffered from intensive taxonomic attention. In many cases a preponderance of names exists for the same species and names are often used incorrectly (see list of synonyms in [4]). During a recent survey of butterflies in Southern Thailand, 150 km north of the Malaysian border, fewer than 50% of the observed butterflies were identified to species [5]. Adding to these difficulties is widespread but inconsistent use of butterfly subspecies names and concepts [6]–[7]. Butterfly surveys in Peninsula Malaysia have not been consistent in using or ignoring subspecies names [8]–[10]. This can make a big difference to biodiversity surveys - if we consider species as the biodiversity unit there are 793 units in Peninsula Malaysia, but if subspecies is considered the biodiversity unit, the number rises to 930 [11].
Butterfly trinomials have traditionally been used to recognize ‘moderate’ morphological differentiation correlated with disjunct geographical distributions [6]–[7], [12]. However, non-discrete morphological variation and the application to contiguously distributed populations, often make subspecies boundaries ambiguous [7]. Following Tobias et al. [13]’s recommendations for avian subspecies delimitation, Braby et al. [7] recently suggested standardized phenotypic criteria for subspecies delimitation in butterflies. Although considered desirable, Braby et al. [7] refrained from setting criteria based on DNA characters, citing a lack of data. However, they did acknowledge that under their concept, subspecies are genetically distinct, but not reciprocally monophyletic according to mitochondrial DNA, noting that lineages possessing a diagnostic morphological character and also showing reciprocal monophyly are probably better regarded as distinct species [7]. This criterion of concordance for species delimitation is in line with “state-of-the-art” practice in taxonomy i.e. the MTMC (Mitochondrial Tree Morphological Character congruence) of Miralles and Vences [14].
Mitochondrial DNA barcodes [15]–[16] are increasingly being used as a supplementary taxonomic identification tool in surveys of Lepidoptera (e.g. [5], [17]–[19]). However, DNA barcoders have often ignored subspecies names [18], [20]–[21], and have used personalized alphanumeric codes for biodiversity units discovered below the traditionally recognized species boundary (e.g. Hamadryas feroniaECO01 [18], [22]). These units used to account for previously overlooked (and possibly cryptic) diversity have come to be known as “dark taxa” [23] and the correspondence between subspecies, recognized by morphological differentiation, and dark taxa is often difficult to resolve (e.g. does H. feroniaECO01 = H. feronia farinulenta? [22]). Most GenBank [24] and BOLD [25] records do not include subspecies names, meaning it is impossible to tell if the authors of the DNA sequence could determine which subspecies the butterfly belonged to or not. It may be possible to narrow down subspecies identity based on locality, but locality is often missing, or imprecise, for GenBank records too.
The aim of this study was to build a DNA barcode reference library for the true butterflies (species from the families – Papilionidae, Pieridae, Nymphalidae, Lycaenidae, Riodinidae) of Peninsula Malaysia from specimens in the Museum of Zoology, University of Malaya (UMKL) collection. We tested the capacity of the library to function as an accurate identification tool for species, screening for signatures of misidentifications, of multiple species sharing identical or very similar DNA barcodes, and of currently unrecognised diversity within the collection. Given the inconsistency in using or ignoring subspecies names in surveys of butterflies, we also explored the value of attaching subspecies names to records in DNA barcode reference libraries. The new DNA barcode library for Peninsula Malaysia was analysed in conjunction with publicly available DNA barcodes from other Asia-Pacific localities to test the ability of the DNA barcodes to discriminate subspecies. Are butterfly subspecies distinctive biodiversity units that can be distinguished by their DNA barcodes and if so, what differentiates them from species? This is an important question. Twenty-eight native butterfly species are currently protected under Malaysian law [26] but in other jurisdictions subspecies can also have legal status [27].
Materials and Methods
Building a DNA Barcode Reference Library for the True Butterflies of Peninsula Malaysia
The UMKL butterfly collection comprises three thousand specimens with representatives of around 30% of the known fauna of Peninsula Malaysia. DNA barcodes were obtained by sampling dry legs from specimens in UMKL. Sampling was restricted to a few specimens per species, including morphologically and geographically diverse specimens where possible. Taxonomy and nomenclature follows our scratchpad [11], (see [28]) and reflects taxonomic decisions since Eliot [4]. The legs were sent to the Canadian Centre for DNA Barcoding for DNA barcode assembly following standard high-throughput protocols for insects [29]. Details of the specimens and DNA barcodes (including GenBank accessions) are available on BOLD [25] in the public dataset: DS-BUTMAY and in Table S1.
We performed an initial screen of the dataset by blasting each new DNA barcode against the full database of BOLD. In cases where new DNA barcodes matched DNA barcodes assigned to a different species name (with >98% similarity) we reexamined the specimens’ morphology to determine the accuracy of the original identifications (provided in the “Taxonomy Note” field of the specimen records on BOLD).
Following this initial screen we subsequently noted cases where specimens currently with different species names have identical or similar DNA barcodes (with >98% similarity) and cases where specimens currently with the same species name have dissimilar DNA barcodes (≤98% similarity). The genetic distances referred to are all K2P corrected (Kimura 2-parameter; as provided by BOLD). We used 2% as the basis for our screening following the example of previous DNA barcoding studies (e.g. [15], [17], [21]–[22], [30]) which have demonstrated that although there is no expectation for a universal threshold of genetic distances between or within species, 2% provides a useful heuristic upon which to base deeper investigation.
Testing if Subspecies can be Distinguished by their DNA Barcodes
By blasting the UMKL DNA barcodes against the full BOLD database we determined which species in the dataset have DNA barcodes on BOLD from other researchers (see Table S1). When a subspecies name was not provided we derived a subspecies name for these DNA barcodes by searching published accounts of the DNA sequences (i.e. journal articles or authors’ websites, e.g. [31]–[32]) and by making inferences based on the reported geographical distribution of the subspecies (e.g. [33], [34]). Note that many DNA barcodes come from GenBank with poorly reliable data, especially imprecise geographical origin, or are “Private” or “Early Release” on BOLD and not publicly viewable, but which nevertheless contribute to a BOLD identification. Where a species from UMKL was determined to be present on BOLD with DNA barcodes from multiple subspecies we then examined if the subspecies were distinguishable based on a “Tree Based Identification” (Neighbor-Joining) in BOLD (see Subspecies Trees S1). Specifically, we observed if each subspecies: i) shared identical DNA barcodes with another subspecies; ii) had unique DNA barcodes but which did not form an exclusive cluster on the tree provided by BOLD; iii) had unique DNA barcodes which formed an exclusive cluster (Figure 1).
Results and Discussion
A DNA Barcode Reference Library for Identification of True Butterflies in Peninsula Malaysia
A DNA barcode was obtained from 458 of 561 specimens (82%) submitted for analysis, accounting for 233 species. While similar to that reported for other Lepidoptera DNA barcoding studies (e.g. [17], [35]), considering that the oldest specimen submitted for analysis was 20 years old the success rate seemed low for a relatively recent collection. This could serve as a warning for those attempting to build a DNA barcode library from tropical museum collections (but see [18]) and has prompted a review of specimen storage conditions at UMKL. An approach that has been suggested is to freeze newly collected butterflies and store them as frozen tissue vouchers rather than the traditional pinning and drying of specimens. DNA extraction, amplification and sequencing using ‘Lep’ primers [29] was highly efficient with fresh (<3 yrs) material. Further mining of public and private collections coupled with targeted field sampling should gradually move the library to completion and increase the number of representatives per species. However, in view of the hyper-diversity of Peninsula Malaysia [4] this is a challenge compared to temperate regions (e.g. the 180 butterfly species of Romania [11]).
Screening the new DNA barcode dataset against the full BOLD database followed by reexamination of morphology revealed that about 15% of specimens in the UMKL collection were originally misidentified. Many of these were nymphalids from the subfamily Satyrinae and the tribe Heliconiini within the Heliconiinae. One noteworthy case was a pierid originally identified as Delias barcasa dives and collected at Genting Highlands, Pahang, in 2011. DNA barcoding conclusively assigned the specimen to Delias agostina (99.3% similarity with DQ082779 from Chiang Mai in northern Thailand [36]) confirming the presence of the species in Peninsula Malaysia. Delias agostina is not included in the plates of D’Abrera [4] but is featured in the Corbet and Pendlebury Delias key with “Burma” printed in bold and in the species checklist with an asterisk, indicating resident status as unconfirmed [4]. Successive screening also revealed several cases of multiple species within the same genus showing identical or similar DNA barcodes.
UMKL DNA barcodes for Danaus melanippus hegesippus shared 99.1% similarity with a “Private” D. genutia DNA barcode from Australia (subspecies not given but probably D. g. alexis [34]) which in turn was >2% distant from UMKL D. g. intermedius DNA barcodes. The Australian subspecies has previously been treated as a distinct species [34]. Interestingly, the phylogenetic sister of D. melanippus, D. affinis (according to [37]), was not the closest matching species, being >2.9% from D. melanippus and >2.6% from D. genutia.
UMKL DNA barcodes recorded under Euploea modesta modesta matched closely (<99.8%) with GenBank DNA barcodes from India recorded under E. core [38] and “Early-Release” DNA barcodes (98.8%) recorded under E. alcathoe and E. core from Australia and Papua New Guinea (Euploea Tree S1). E. m. modesta is found in India, E. m. lugens in Australia and Papua New Guinea. Similarly, the single short UMKL DNA barcode (307 bp) for E. camaralzeman malayica matched closely (99.6%) with a “Published” DNA barcode for E. core from Thailand and matched 100% to other “Early-Release” E. core DNA barcodes on BOLD. Furthermore, the UMKL DNA barcodes for E. doubledayi evalina matched 100% to E. algea (KC306717) from India and yet another “Early-Release” E. core from Australia. There was a further distinct cluster of E. core on BOLD containing DNA barcodes from Australia and Thailand which was distant from all the UMKL Euploea. One UMKL DNA barcode recorded under Euploea eunice leucogonis and collected from Genting Highlands, Pahang, in 2012 was distant (3.3%) from the two other UMKL E. eunice leucogonis DNA barcodes (Euploea Tree S1), which themselves were similar (99.2%) to E. kluji from India (KC306728) but relatively distant (98.0%) from E. kluji from Southern Thailand (HQ962260). The morphologically similar, and one time synonym [4], E. leucostictos formed a distinct sister to this cluster. As wittily noted by Corbet and Pendlebury (2nd edition) in the legend to Plate 23 [39], “it is easier to ascertain the country of origin of a (Euploea) specimen than to determine its specific identity”, the genus is notorious for being taxonomically difficult. Any taxonomic interpretation is further complicated by reports of hybrids [40] and the fact that species are commonly reared for butterfly parks (and released). There may be a tendency for collectors to assign difficult specimens to the most common species - E. core - accounting for its appearance in many places in this screening.
Identification of Eurema species, a genus found abundantly in disturbed and undisturbed habitats alike, is also notoriously difficult [4], [41]. UMKL DNA barcodes recorded under three species of Eurema (E. ada iona, E. hecabe contubernalis, E. lacteola lacteola) showed low divergence amongst themselves and also with various Eurema species from various Asia-Pacific localities. The DNA barcodes all sat within the same BIN (Barcode Index Number) (BIN S1); the system on BOLD which clusters DNA barcodes into operational taxonomic units closely corresponding to traditionally recognized species [42]. A review of Eurema in Peninsula Malaysia is currently underway by our research group. Whether Eurema as an example of ‘barcode sharing’ is actually a reflection of the difficulty assigning these small yellow butterflies to species on the basis of wing patterns remains to be seen.
Loxura atymnus fuconius and L. cassiopea cassiopea are morphologically similar [4] and the UMKL specimen of L. atymnus fuconius was originally recorded under L. cassiopea cassiopea. However, these species cannot be confused as the wing patterns, when studied carefully, and the DNA barcodes, although close (1.7% distant and in the same BIN), are characteristic for each species.
In the UMKL dataset the single representative of Polyura athamas athamas was distant from the P. a. uraeus DNA barcodes (2.1%) and closer to P. hebe (1.7%). Like Eliot [4] we are hesitant to draw conclusions about the species status of these two taxa, in our case because of the small number of specimens available in UMKL and because only a short DNA barcode (307 bp) was generated for the P. a. athamas specimen. However, these taxa are easily distinguished as the wing patterns and the DNA barcodes, although close, are characteristic for each taxon (Figure 2).
The BOLD Process ID is followed by the sequence length.
UMKL DNA barcodes recorded under four species of Tanaecia (T. aruna aruna, T. iapis puseda, T. munda waterstradti, T. pelea pelea) sat in the same BIN along with three DNA barcodes from Thailand, also representing multiple species (BIN S2). The taxonomy of this genus is difficult [43], with species specific diagnostic characters mostly from the male genitalia [4], [43] (not studied here), and needs further investigation.
Non-monophyly of Charaxes bernardus has been reported before, with C. marmax nested within C. bernardus on the molecular phylogenetic tree of Aduse-Poku et al. [44]. We found that C. durnfordi durnfordi and C. bernardus crepax shared identical DNA barcodes, despite very distinctive wing patterns (Figure 3). This interesting and rare pattern deserves further study and may reflect the complex biogeographical history of this genus [44] or mitochondrial introgression.
The BOLD Process ID or GenBank Accession (GB) is followed by locality.
UMKL DNA barcodes recorded as Mycalesis mineus macromalayana sat in a BIN with GenBank DNA barcodes for M. mineus from India, but the BIN also contained DNA barcodes from GenBank recorded under M. visala, M. intermedia and M. perseoides (BIN S3). Also present were unpublished M. mineus and M. panthaka DNA barcodes from China. Like the other genera above the Malaysian Mycalesis have a long history of taxonomic difficulty [45].
UMKL DNA barcodes recorded as Tirumala septentrionis septentrionis, the only common Tirumala in Peninsula Malaysia, sat in a BIN with “Early Release” T. hamata DNA barcodes from Australia and Papua New Guinea [35] (BIN S4). T. septentrionis septentrionis overwintering in Taiwan has previously been treated as T. hamata septentrionis [46]. T. limniace, a similar looking species, DNA barcodes from India were also in the BIN and may be misidentifications.
UMKL Troides helena cerberus DNA barcodes matched closely (>98.8%) with GenBank and BOLD T. oblongomaculatus from Indonesia. These closely related species have been treated historically as a single species [34]. T. oblongomaculatus, a “relic race of uncertain status” [47], has been reported to hybridize, including with taxonomically distant species [48].
UMKL DNA barcodes recorded under Ypthima horsfieldi humei shared close similarity (>99.5%) with Ypthima nebulosa DNA barcodes from Thailand [5]. Y. nebulosa has not been reported for Peninsula Malaysia [34] but according to Corbet and Pendlebury is likely to be found in the region [4] suggesting the specimens in UMKL require further evaluation.
Screening against BOLD highlighted 27 other species with unique DNA barcodes but which were <2% distant from other species. These represented borderline cases for the screening threshold which were nevertheless allocated to different BINs by BOLD (see Table S1; Figure 4) and cases associated with short sequence lengths or suspected misidentified DNA barcodes on GenBank/BOLD (see Table S1; Figure 4).
The DNA barcodes are arranged alphabetically by species name along the horizontal axes. Conspecific similarities below 98% and interspecific similarities above 98% that were associated with different subspecies, misidentified BOLD barcodes, potential cryptic species, short sequence length, barcode sharing and different BINs are highlighted with different coloured data points.
Within the new Peninsula Malaysia dataset, only three species showed wide (>2%) conspecific distances: Euploea eunice, Polyura athamas (see above) and Hebomoia glaucippe. DNA barcodes for H. g. anomala found on Pulau Aur, Johor, a small island off the east coast of mainland Peninsula Malaysia, were 4.2% distant from the DNA barcode for H. g. aturia from the mainland which clustered closely with BOLD DNA barcodes from Thailand, most likely H. g. aturia, and different subspecies from Taiwan and China (Figure 5). The differences in wing pattern between these two groups are readily apparent with the Pulau Aur butterflies exhibiting a deeper yellow upperside [4] (Figure 5). H. g. anomala was described as a distinct species by Pendlebury in 1939 [34].
The BOLD Process ID or GenBank Accession (GB) is followed by locality.
Compared with the levels of cryptic diversity discovered in other DNA barcoding surveys (e.g. [18], [22], [49]) three species showing wide conspecific distances is relatively few, suggesting that the long history of taxonomic study of the butterflies of Peninsula Malaysia has led to a relatively accurate account of species diversity. Furthermore, two of the three cases, Polyura and Hebomoia, were associated with subspecies, one of which had previously been treated as distinct species. There were no other cases of species being represented by more than one subspecies in the UMKL collection. Perhaps the one case of truly unrecognized diversity within the Peninsula Malaysia dataset was the distinct DNA barcodes within Euploea eunice leucogonis and this deserves further study to determine if this is truly the exception.
Following correction of morphological misidentifications in UMKL, the DNA barcodes for 78% of the 233 species were unique (with non-overlapping conspecific and interspecific distances for multiple representatives) when compared with conspecifics and closest matches on BOLD (Table S1; Figure 4). Excluding outliers - confirmed or probable misidentified DNA barcodes on BOLD and conspecific distances associated with divergent subspecies or cryptic species diversity - the number of distinct species rises to 92%, validating the capacity of the DNA barcode reference library for rapid and effective assignment of true butterflies to species. The few cases of ‘barcode sharing’ that remain provide stimulus for subsequent studies. Considering the importance of butterflies as bioindicators and conservation flagships we are particularly encouraged by the potential of DNA barcoding to enable local species inventories, without the need for lethal sampling [50], but with much higher accuracy and precision than can be achieved by observing butterflies on the wing [1], [5], or even by traditional morphological identification (considering the misidentifications in UMKL).
Can Subspecies be Distinguished by their DNA Barcodes and What Differentiates them from Species?
There were 1189 DNA barcodes on BOLD for the 233 UMKL species and we determined that 80 species were represented by multiple subspecies (Table S1). Of the 192 subspecies, 86 were represented by singletons and 87% of these singletons had unique DNA barcodes. Of the 106 subspecies represented by multiple DNA barcodes 81% had unique DNA barcodes not shared with other subspecies and 66% formed exclusive clusters on identification trees (Subspecies Trees S1; Figure 6). Because many of the subspecies were represented by singletons and “Early-Release” or “Private” DNA barcodes on BOLD, it is outside the scope of this study to examine how many of the subspecies would be reciprocally monophyletic for mtDNA in phylogenetic (maximum or statistical parsimony) analyses. However, under current levels of representation, most subspecies are genetically distinct for mtDNA (Figure 6) which is in accordance with the expectations of the butterfly subspecies concept of Brady et al. [7]. How this pattern changes or stabilizes as BOLD continues to grow will clarify the nature of the relationship between DNA barcodes and subspecies more accurately. The results suggest that as subspecies move from singletons to multiple representatives the number of subspecies with unique DNA barcodes could decrease (87% versus 81%; Figure 6). Because the butterfly DNA barcodes available on BOLD came from a range of local surveys or phylogenetic studies, the geographic coverage was patchy and no biogeographic patterns were apparent from the analysis. However, it was not uncommon for UMKL DNA barcodes to be similar to conspecific DNA barcodes from India or China, at the extremities of the Asia-Pacific region while distinct subspecies were from one of the region’s many islands.
Of the subspecies with unique DNA barcodes many were highlighted when we screened the UMKL DNA barcode dataset against the full BOLD database using the >2% conspecific distance threshold (Table 1). Conspecific genetic distances of this magnitude, i.e. distances typically seen between species, would normally warrant “dark taxon” status in DNA barcoding studies and some of these cases have in fact been highlighted by previous studies (e.g. several species from Western Ghats, India [39]). Historical studies have likewise highlighted the morphological distinctiveness of these taxa as implied through their current disparate subspecies designations. Many of these subspecies had previously been treated as distinct species (Table 1), and the DNA barcode data supports a re-evaluation of their status. Similarly, ‘unrecognized’ lepidopteran diversity revealed through DNA barcoding in the other surveys (e.g. [18], [22]) had been recognized previously, although as subspecies taxa, or as sunken or forgotten names [51]. This may reflect the challenge of meshing Linnaean taxonomy with DNA taxonomy systems [52]. In these cases above, consistent application of subspecies names in DNA barcode reference libraries would negate the need for dark taxon designation. Following a reverse MTMC [14], DNA barcoding could provide a means of testing, through concordance, if subspecies, established on the basis of moderate morphological differentiation between localities [7], are of sufficient evolutionary independence to merit species status. Note that 13 other species had specimens with unique DNA barcodes but which were >2% distant from conspecific DNA barcodes. However, these further cases were due to confirmed or suspected misidentified DNA barcodes on GenBank/BOLD (See Table S1).
DNA Barcode Reference Libraries and Subspecies
In this study we present a preliminary DNA barcode library for a major component of the true butterfly species of Peninsula Malaysia. The majority of species and subspecies sampled possessed unique DNA barcodes. Although there is no fixed threshold of genetic distances clearly differentiating conspecific from interspecific distances, BOLD identification trees generally show a discernible pattern of low conspecific distances compared to interspecific distances, so can enable effective assignment of unknown DNA barcodes to species, especially when examined in conjunction with the BIN system.
Unlike assignment to a species, assignment of an unknown DNA barcode to a subspecies using a BOLD identification tree would not be easily accomplished. The genetic distances between most conspecific subspecies are small and indistinguishable from distances between members of the same subspecies. Although the majority of subspecies with multiple representatives formed exclusive clusters on Neighbor-joining trees in our analyses, forming an exclusive cluster cannot logically guide taxonomic assignments in the absence of other discernible patterns - exclusive clusters are present at, and between, all taxonomic levels on a tree [19].
Those subspecies that show ‘large’ inter-taxa distances probably warrant full species status. Conversely, there are undoubtedly cases where subspecies names are applied to groups that probably do not warrant taxonomic recognition [7], [12]. For example, considering the similarity of Loxura atymnus and L. cassiopea the necessity for finer taxonomic divisions [53] is dubious. Subspecies designations as shorthand for geographically and morphologically differentiated groups provide a useful heuristic for assessing how such groups correlate with clustering patterns of DNA barcodes, especially as the number of DNA barcodes per species in reference libraries increases. Considering this, we feel there is significant value in attaching subspecies names to records in DNA barcode databases. A beneficial addition to BOLD would be the facility to allow data contributors to specify subspecies names while still recognising that members of different subspecies are conspecific for the purpose of progress statistics and other analytical tools.
Supporting Information
Table S1.
Details of the specimens codes, including GenBank accessions, used in this study and results of screening of DNA barcodes.
https://doi.org/10.1371/journal.pone.0079969.s001
(XLSX)
BIN S1.
Tree for Barcode Index Number - BIN6082[BOLD:AAA6082]
https://doi.org/10.1371/journal.pone.0079969.s002
(PDF)
BIN S2.
Tree for Barcode Index Number - BIN511656[BOLD:ABZ1656]
https://doi.org/10.1371/journal.pone.0079969.s003
(PDF)
BIN S3.
Tree for Barcode Index Number - BIN108780[BOLD:AAK8780]
https://doi.org/10.1371/journal.pone.0079969.s004
(PDF)
BIN S4.
Tree for Barcode Index Number - BIN21308[BOLD:AAC1308]
https://doi.org/10.1371/journal.pone.0079969.s005
(PDF)
Euploea Tree S1.
Result of BOLD tree based identification of a UMKL DNA barcode for Euploea modesta modesta
https://doi.org/10.1371/journal.pone.0079969.s006
(PDF)
Subspecies Trees S1.
Trees resulting from a BOLD tree based identification of UMKL DNA barcodes
https://doi.org/10.1371/journal.pone.0079969.s007
(PDF)
Acknowledgments
A. Sasekumar, Rosli Ramli, Thary Gazi, and Evelyn Leong provided invaluable assistance during this study. Anonymous reviewer 3 provided helpful comments which greatly improved this manuscript.
Author Contributions
Conceived and designed the experiments: JJW MSA. Performed the experiments: JJW KWS. Analyzed the data: JJW KWS. Contributed reagents/materials/analysis tools: JJW MSA. Wrote the paper: JJW.
References
- 1. Cleary DFR (2004) Assessing the use of butterflies as indicators of logging in Borneo at three taxonomic levels. Journal of Economic Entomology 97(2): 429–435.
- 2.
Distant WL (1882–1886) Rhopalocera malayana: A description of the butterflies of the Malay Peninsula. Penang: Logan. 548 p.
- 3.
Corbet AS, Pendlebury HM (1934) The butterflies of the Malay Peninsula. Kuala Lumpur: Kyle, Palmer & Co. 304 p.
- 4.
Corbet AS, Pendlebury HM, Eliot JN (1992) The butterflies of the Malay Peninsula 4th edition. Kuala Lumpur: Malayan Nature Society. 664 p.
- 5. Basset Y, Eastwood R, Sam L, Lohman DJ, Novotny V, et al. (2011) Comparison of rainforest butterfly assemblages across three biogeographical regions using standardized protocols. Journal of Research on the Lepidoptera 44: 17–28.
- 6. Vane-Wright RI (2003) Indifferent Philosophy versus Almighty Authority: On consistency, consensus and unitary taxonomy. Systematics and Biodiversity 1 (1): 3–11.
- 7. Braby MF, Eastwood R, Murray N (2012) The subspecies concept in butterflies: has its application in taxonomy and conservation biology outlived its usefulness? Biological Journal of the Linnean Society 106 (4): 699–716.
- 8. Quek KC, Sodhi NS, Liow LH (1999) New records of butterfly species for Pulau Tioman, Peninsular Malaysia. Raffles Bulletin of Zoology S6: 271–276.
- 9.
Sofian-Azirun M, Khaironizam MZ, Norma-Rashid Y, Diacus B (2005) Butterflies (Insecta: Lepidoptera) of the southwestern Endau-Rompin Park, Johor, Malaysia. In: Mohamed H, Zakaria-Ismail M, editors. The Forests and Biodiversity of Selai, Endau-Rompin. Kuala Lumpur: University of Malaya. 169–178.
- 10. Sofian-Azirun M, Normaisarah I, Khaironizam MZ (2009) A checklist of butterfly (Insecta: Lepidoptera) from Kenaboi Forest Reserve, Negeri Sembilan. Malaysian Journal of Science 28: 415–425.
- 11.
Wilson JJ (2013) Butterflies of Peninsula Malaysia. Available: http://malaysiabutterflies.myspecies.info. Accessed 16 June 2013.
- 12. Kodandaramaiah U, Weingartner E, Janz N, Leski M, Slove J, et al. (2012) Investigating concordance among genetic data, subspecies circumscriptions and hostplant use in the nymphalid butterfly Polygonia faunus. PLoS ONE 7(7): e41058.
- 13. Tobias JA, Seddon N, Spottiswoode CN, Pilgrim JD, Fishpool LDC, et al. (2010) Quantitative criteria for species delimitation. Ibis 152: 724–746.
- 14. Miralles A, Vences M (2013) New metrics for comparison of taxonomies reveal striking discrepancies among species delimitation methods in Madascincus lizards. PLoS ONE 8(7): e68242.
- 15. Hebert PDN, Cywinska A, Ball SL, deWaard JR (2003) Biological identifications through DNA barcodes. Proc Roy Soc Lond B 270: 313–321.
- 16.
Floyd R, Wilson JJ, Hebert PDN (2009) DNA barcodes and insect biodiversity. In: Foottit RG, Adler PH, editors. Insect Biodiversity: Science and Society. Oxford: Blackwell Publishing. 417–431.
- 17. Hausmann A, Haszprunar G, Hebert PDN (2011) DNA barcoding the geometrid fauna of Bavaria (Lepidoptera): successes, surprises, and questions. PLoS ONE 6(2): e17134.
- 18. Janzen D, Hallwachs W, Blandin P, Burns JM, Cadiou J-M, et al. (2009) Integration of DNA barcoding into an ongoing inventory of complex tropical biodiversity. Mol Ecol Res 9: 1–25.
- 19. Wilson JJ, Rougerie R, Schonfeld J, Janzen DH, Hallwachs W, et al. (2011) When species matches are unavailable are DNA barcodes correctly assigned to higher taxa? An assessment using sphingid moths. BMC Ecology 11: 18.
- 20. Lukhtanov VA, Sourakov A, Zakharov EV, Hebert PDN (2009) DNA barcoding Central Asian butterflies: increasing geographical dimension does not significantly reduce the success of species identification. Mol Ecol Res 9: 1302–1310.
- 21.
Dincă V, Zakharov E, Hebert PDN, Vila R (2010) Complete DNA barcode reference library for a country's butterfly fauna reveals high performance for temperate Europe. Proc R Soc Lond B doi:https://doi.org/10.1098/rspb.2010.1089 1471–2954.
- 22. Prado BR, Pozo C, Valdez-Moreno M, Hebert PDN (2011) Beyond the colours: discovering hidden diversity in the Nymphalidae of the Yucatan Peninsula in Mexico through DNA barcoding. PLoS ONE 6: e27776.
- 23. Maddison DR, Guralnick R, Hill A, Reysenbach A-L, McDade LA (2012) Ramping up biodiversity discovery via online quantum contributions. Trends Ecol Evol 27: 72–77.
- 24. Benson DA, Cavanaugh M, Clark K, Karsch-Mizrachi I, Lipman DJ, et al. (2013) GenBank. Nucl Acids Res 41 (D1): D36–D42.
- 25. Ratnasingham S, Hebert PDN (2007) BOLD: The Barcode of Life Data System (www.barcodinglife.org). Mol Ecol Notes. 7: 355–364.
- 26.
Kato K, Raman S (2005) A guide to the butterflies of Langkawi. Kuala Lumpur. 319.
- 27. Haig SM, Beever EA, Chambers SM, Draheim HM, Dugger BD, et al. (2006) Taxonomic considerations in listing subspecies under the U.S. Endangered Species Act. Conservation Biology 20: 1584–1594.
- 28. Smith VS, Rycroft SD, Brake I, Scott B, Baker E, et al. (2011) Scratchpads 2.0: a virtual research environment supporting scholarly collaboration, communication and data publication in biodiversity science. ZooKeys 150: 53–70.
- 29. Wilson JJ (2012) DNA barcodes for insects. Methods in Molecular Biology 858: 17–46.
- 30.
Hebert PDN, deWaard JR, Landry J-F (2009) DNA barcodes for 1/1000 of the animal kingdom. Biol Lett doi:https://doi.org/ 1744-957X.
- 31. Peña C, Malm T (2012) VoSeq: a voucher and DNA sequence web application. PLoS ONE 7(6): e39071.
- 32.
Nazari V (2006) E. H. Strickland Entomological Museum: DNA vouchers. Available: http://www.biology.ualberta.ca/uploads//uasm/Vouchers/Accessed 1 August 2013.
- 33.
Inayoshi Y (2012) A check list of butterflies in Indo-China. Available: http://yutaka.it-n.jp/Accessed 1 August 2013.
- 34.
Savela M (2013) Lepidoptera and some other life forms. Available: http://www.nic.funet.fi/pub/sci/bio/life/insecta/lepidoptera/Accessed 1 August 2013.
- 35. Hebert PDN, deWaard JR, Zakharov EV, Prosser SWJ, Sones JE, et al. (2013) A DNA ‘barcode blitz’: rapid digitization and sequencing of a natural history collection. PLoS ONE 8(7): e68535.
- 36. Braby MF, Pierce NE (2007) Systematics, biogeography and diversification of the Indo-Australian genus Delias Hubner (Lepidoptera: Pieridae): phylogenetic evidence supports an ‘out-of-Australia’ origin. Syst Entomol 32(1): 2–25.
- 37. Smith DAS, Lushai G, Allen JA (2005) A classification of Danaus butterflies (Lepidoptera: Nymphalidae) based upon data from morphology and DNA. Zool J Linn Soc 144: 191–212.
- 38. Gaikwad SS, Ghate HV, Ghaskadbi SS, Patole MS, Schouce YS (2012) DNA barcoding of nymphalid butterflies (Nymphalidae: Lepidoptera) from Western Ghats of India. Mol Biol Rep 39: 2375–2383.
- 39.
Corbet AS, Pendlebury HM (1956) The butterflies of the Malay Peninsula 2nd edition. London: Oliver and Boyd. 537 p.
- 40.
Scheermeyer E (1999) The crows, Euploea species, with notes on the blue tiger, Tirumala hamata (Nymphalidae: Danainae). In: Kitching RL, Scheermeyer E, Jones RE, Pierce NE, editors. Biology of Australian Butterflies. Monographs on Australian Lepidoptera. Volume 6. Melbourne: CSIRO Publishing. 191–216.
- 41. Narita S, Nomura M, Kato Y, Yata O, Kageyama D (2007) Molecular phylogeography of two species of Eurema butterflies. Genetica 131: 241–253.
- 42. Ratnasingham S, Hebert PDN (2013) A DNA-based registry for all animal species: the Barcode Index Number (BIN) system. PLoS ONE 8(7): e66213.
- 43. Corbet AS (1941) XXXIII. - A revision of the Malaysian genus Tanaecia Btlr. (Lepidoptera: Nymphalidae). Journal of Natural History Series 11 7(42): 507–520.
- 44. Aduse-Poku K, Vingerhoedt E, Wahlberg N (2009) Out-of-Africa again: a phylogenetic hypothesis of the genus Charaxes (Lepidoptera: Nymphalidae) based on five gene regions. Molecular Phylogenetics and Evolution 53: 463–478.
- 45. Corbet AS (1937) Observations on species of Satyridae and Amathusiidae from the Malay Peninsula. Proc R Entomol Soc Lond Ser B Taxonomy 6: 96–99.
- 46. Wang HY, Emmel TC (1990) Migration and overwintering aggregations of nine danaine butterfly species in Taiwan (Nymphalidae). Journal of the Lepidopterists' Society 44(4): 216–228.
- 47. Vane-Wright RI, de Jong R (2003) The butterflies of Sulawesi: annotated checklist for a critical island fauna. Zool Verh Leiden 343: 3–267.
- 48. Sands DPA, Sawyer PF (1977) An example of natural hybridization between Troides oblongomaculatus papuensis Wallace and Ornithoptera priamus poseidon Doubleday (Lepidoptera: Papilionidae). Australian Journal of Entomology 16(1): 81–82.
- 49. Hajibabaei M, Janzen DH, Burns JM, Hallwachs W, Hebert PDN (2006) DNA barcodes distinguish species of tropical Lepidoptera. PNAS 103(4): 968–971.
- 50. Koscinski D, Crawford LA, Keller HA, Keyghobadi N (2009) Effects of different methods of non-lethal tissue sampling on butterflies. Ecological Entomology 36: 301–308.
- 51. Wilson JJ (2011) Taxonomy and DNA sequence databases: a perfect match. Terrestrial Arthropod Reviews 4: 221–263.
- 52. Vogler AP, Monaghan MT (2007) Recent advances in DNA taxonomy. J Zool Syst Evol Res 45(1): 1–10.
- 53. Hayashi H (1976) Two new subspecies of Loxura cassiopia Distant from Palawan and Mindanao (Lepidoptera: Lycaenidae). Trans Lep Soc Japan 47(1): 1–3.
- 54. Eliot JN (1986) A review of the Miletini (Lepidoptera: Lycaenidae). Bull Br Mus Nat Hist 53: 1–105.
- 55. Yata O, Chainey JE, Vane-Wright RI (2010) The Golden and Mariana albatrosses, new species of pierid butterflies, with a review of subgenus Appias (Catophaga) (Lepidoptera). Systematic Entomology 35: 764–800.
- 56. Ohshima I, Tanikawa-Dodo Y, Saigusa T, Nishiyama T, Kitani M, et al. (2010) Phylogeny, biogeography, and host-plant association in the subfamily Apaturinae (Insecta: Lepidoptera: Nymphalidae) inferred from eight nuclear and seven mitochondrial genes. Mol Phylogenet Evol 57 (3): 1026–1036.
- 57. Kim MI, Wan X, Kim MJ, Jeong HC, Ahn N-H, et al. (2010) Phylogenetic relationships of true butterflies (Lepidoptera: Papilionoidea) inferred from COI, 16S rRNA and EF-1α sequences. Mol Cells 30: 409–425.
- 58. Tsao W-C, Yeh W-B (2008) DNA-based discrimination of subspecies of swallowtail butterflies (Lepidoptera: Papilioninae) from Taiwan. Zoological studies 47(5): 633–643.
- 59. Peña C, Wahlberg N, Weingartner E, Kodandaramaiah U, Nylin S, et al. (2006) Higher level phylogeny of Satyrinae butterflies (Lepidoptera: Nymphalidae) based on DNA sequence data. Mol Phylogenet Evol 40: 29–49.