Skip to main content

DNA barcoding a unique avifauna: an important tool for evolution, systematics and conservation



DNA barcoding utilises a standardised region of the cytochrome c oxidase I (COI) gene to identify specimens to the species level. It has proven to be an effective tool for identification of avian samples. The unique island avifauna of New Zealand is taxonomically and evolutionarily distinct. We analysed COI sequence data in order to determine if DNA barcoding could accurately identify New Zealand birds.


We sequenced 928 specimens from 180 species. Additional Genbank sequences expanded the dataset to 1416 sequences from 211 of the estimated 236 New Zealand species. Furthermore, to improve the assessment of genetic variation in non-endemic species, and to assess the overall accuracy of our approach, sequences from 404 specimens collected outside of New Zealand were also included in our analyses. Of the 191 species represented by multiple sequences, 88.5% could be successfully identified by their DNA barcodes. This is likely a conservative estimate of the power of DNA barcoding in New Zealand, given our extensive geographic sampling. The majority of the 13 groups that could not be distinguished contain recently diverged taxa, indicating incomplete lineage sorting and in some cases hybridisation. In contrast, 16 species showed evidence of distinct intra-species lineages, some of these corresponding to recognised subspecies. For species identification purposes a character-based method was more successful than distance and phylogenetic tree-based methods.


DNA barcodes accurately identify most New Zealand bird species. However, low levels of COI sequence divergence in some recently diverged taxa limit the identification power of DNA barcoding. A small number of currently recognised species would benefit from further systematic investigations. The reference database and analysis presented will provide valuable insights into the evolution, systematics and conservation of New Zealand birds.


DNA barcoding sensu Hebert et al. [1] has been suggested as a means of species identification through comparison of a standardised segment of the mitochondrial genome. In the case of animals, the ‘barcode’ is a 648 bp region of the 5′ end of the cytochrome c oxidase I (COI) gene. Since its proposal, DNA barcoding has become a large scale and well-supported global enterprise [2]. DNA barcoding has two distinct goals; species discovery and specimen identification [1, 3, 4]. The former, which involves using DNA barcodes to delimit species boundaries or identify novel species has been criticised, for among other reasons, being a form of DNA taxonomy and for relying on a single gene to infer species relationships [3, 5, 6]. Although DNA barcoding does not provide a way of defining new species, the results of such studies can highlight taxa that require further investigation. When applied to the latter problem of identifying specimens within taxonomically well-resolved groups, DNA barcoding has proven to be a very useful tool [7, 8].

Traditional taxonomic identification requires increasingly rare expert knowledge and is often difficult or impossible for degraded specimens or incomplete remains. As only a small amount of DNA is required, samples that would usually be difficult or impossible to identify morphologically such as blood, eggs, embryos, feathers and faeces can be accurately identified by DNA barcoding. DNA barcoding has been successfully applied to a variety of issues, such as the identification of historic specimens [9, 10], wildlife forensics [11,12,13], diet analysis [14, 15], identification of species involved in birdstrike (a collision between a bird and an aircraft) [16, 17] and conservation biology (reviewed in Krishnamurthy et al., [18]). In cases where DNA is highly degraded, a shorter “mini-barcode” may still enable specimen identification [19]. Furthermore, where DNA barcodes have highlighted inconsistences with established taxonomy, more detailed studies using a range of approaches have been undertaken and in many cases have been able to inform the processes of molecular evolution, biogeography and speciation (reviewed in Barreira et al., [20]).

Debate has centred on the best way to use DNA barcodes for species identification. Early studies analysed barcodes exclusively using distance based methods that numerically quantify the degree of genetic divergence between taxa e.g. [1, 21]. However, character-based methods that rely on the presence or absence of diagnostic characters (in this case nucleotides), are considered more consistent with modern taxonomy [22]. Many early studies also reported the existence of a global ‘barcode gap’, a discontinuity between intra- and interspecific genetic divergences. However, most of these studies had limited congeneric and geographic sampling resulting in underestimation of intraspecific variation and overestimation of interspecific divergence [8]. Subsequent studies have found that within well-sampled groups, intra- and interspecific distances usually overlap significantly so that no global barcode gap exists [8]. However, when used in combination with character-based methods, distance based analyses can still provide useful insights [4].

Avian taxonomy is relatively well-resolved making it an ideal group with which to test the efficacy of DNA barcoding for specimen identification [23]. The All Birds Barcoding Initiative ( was launched in 2005 and so far the avifauna of a variety of different geographic regions has been successfully DNA barcoded including North America, the eastern Palearctic, the Neotropics, Scandinavia, the Netherlands, Japan and Turkey [21, 23,24,25,26,27,28,29,30]. While methodology differs between each study, generally they report high success rates for species identification between 93% (520 species) [26] and 96.6% (226 species) [29].

The avifauna of New Zealand is evolutionarily and taxonomically distinct. After the continent of Zealandia split from Gondwana approximately 83 million years ago [31], it became the largest landmass free from ground-dwelling mammals allowing the avifauna to flourish [32]. Today, New Zealand is an archipelago of two main islands and over 330 smaller ones, with a total land area of approximately 270,000 km2, separated from any other significant land mass by almost 1500 km [33, 34]. Despite this geographic isolation, the region has not been completely isolated biologically, as demonstrated by the heterogeneous composition of the modern avifauna which consists of representatives from globally diverse taxa [35]. Although there is strong evidence for vicariant speciation in some groups, other taxa dispersed to New Zealand following the break-up of Gondwana with the majority arriving from Australia or the Pacific [35]. There is a high degree of endemism (of 168 contemporary native bird species, 93 are endemic [36]), which is also indicative of isolation.

Many features of the New Zealand avifauna are reflective of the country being an archipelago. As with other islands, representation of groups is highly variable and the overall diversity of some groups is low [35]. The numerous offshore islands have facilitated allopatric divergence, with some island taxa being recognised as separate species from their mainland New Zealand relatives [37]. These islands provide breeding grounds for many seabird species, and as a result New Zealand is often referred to as the ‘Seabird Capital of the World’ [38]. Nearly a quarter of the world’s 359 seabird species breed in New Zealand and almost 10% breed exclusively in New Zealand [38]. Unfortunately, 80% of New Zealand’s native birds are now either ‘threatened’ or ‘at risk’, mostly as a consequence of predation by introduced mammalian predators [36]. Native birds are a large part of New Zealand’s national identity and the country’s strong conservation ethos has established it as a world-leader in avian conservation [39].

The composition and evolutionary history of the New Zealand avifauna is very different from that of other regions where DNA barcoding of birds has been successful. With the exception of Saitoh et al. [29], most studies have focused on continental regions. New Zealand however, is a continental island [40], and its avifauna has characteristics of both a continental remnant and an isolated archipelago [32]. Additionally, seabirds which make up a large portion of native species, have very different life history traits and population dynamics than land birds [41]. These features make it difficult to predict the success of DNA barcodes for species identification in New Zealand. The present study aims to: 1) develop a working DNA barcoding database for the birds of New Zealand; 2) determine the percentage of currently recognized species that can be discriminated by DNA barcoding; 3) test the potential of DNA barcodes to correctly assign specimens to their nominal species; 4) identify taxa that could benefit from further investigation.


COI sequence data was obtained from 1416 specimens representing 211 avian species found within the New Zealand region. Over 90% of these species were represented by > 2 specimens (Table 1). Where available, the sister species or a close relative of all New Zealand species were included in the analyses (an additional 404 sequences from 107 species). Data was analysed using three methods. The first two methods are based on analysis of genetic distances (pairwise distance analysis and neighbour-joining tree building). The third is a diagnostic character assignment method implemented in the program CAOS [42,43,44]. The mean intraspecific uncorrected p-distance was 0.32% (range 0.00–7.94%) and the mean nearest neighbour (i.e. minimum interspecific) p-distance was 4.24% (range 0.00–13.27%). There was substantial overlap between these values (see Additional file 1). The optimised threshold was 0.25% with a cumulative error rate of 15.8% (see Additional file 2).

Table 1 Summary of species used in this study, including sequences obtained from Genbank (New Zealand endemic species are also by definition New Zealand native species)

The local barcode gap (Fig. 1a) reflects whether or not, within each species, the genetic distance between each conspecific individual is smaller than to any allospecific individual [4]. For New Zealand species with > 1 specimen, 17.8% did not have a local barcoding gap meaning that the difference between the maximum intraspecific and the minimum interspecific distances for that species was ≤0 (Table 2). There was no correlation between the number of specimens per species and the maximum intraspecific distance (Pearsons correlation coefficient 0.13; p-value = 0.07) (Fig. 1b).

Fig. 1

Distance analysis of COI data. a Comparison of nearest neighbour (minimum interspecific distance) and maximum intraspecific distances of the COI sequences from each New Zealand species represented by > 1 specimen (n = 191). Equal intra– and inter–specific variation is marked by the black line. Points above the black line indicate species with ‘local barcode gaps’. b Comparison of maximum intraspecific distance and sampling effort (number of specimens) for each species. There is no observable sampling bias in levels of intraspecific variation. In both scatterplots, green points represent species with a local barcode gap, while red points represent those with no barcode gap

Table 2 Groups of New Zealand bird species with limited COI divergence. For each species the number of specimens analysed is indicated, as is the neighbour-joining tree profile (≥ 95 = monophyletic with greater than or equal to 95% bootstrap support, < 95 = monophyletic with less than 95% bootstrap support). Whether a species had a local barcode gap and could be reliably identified by CAOS is indicated. Maximum intraspecific as well as both minimum and mean distances between the species are given in percentages. The potential reason(s) for the observed similarity in barcodes is provided along with supporting references [47,48,49,50,51,52, 54,55,56,57, 60,61,62, 64, 65, 67, 68, 98,99,100,101,102,103,104]

Of the 191 New Zealand species represented by > 1 specimen, 134 (70.2%) formed well supported monophyletic groups (≥95% bootstrap support) and 29 (15.2%) were monophyletic but with < 95% support (Additional file 3). Nine species (4.7%) were paraphyletic and the remaining 19 species (9.9%) were polyphyletic (Table 2 and Additional file 3). For species represented by only one specimen, no bootstrap support could be calculated. However, with the exception of the fulmar prion (Pachyptila crassirostris) and the Chatham Island pigeon (Hemiphaga chathamensis), sequences from these single specimens formed distinct branches in the tree and did not interfere with other groupings (Additional file 3). Species with two or more distinct clusters in the neighbour-joining tree, supported by high bootstrap values, were identified as candidates for further investigation. Sixteen species showed evidence of two or more divergent lineages (> 1.7% divergence with > 82% bootstrap support) (Table 3). The majority of these groupings corresponded to recognised subspecies and/or populations separated by large geographic distances.

Table 3 List of New Zealand species that show geographically structured populations or divergent lineages. Indicated are New Zealand specific lineages, splits within New Zealand and the presence of subspecies recognised by Clements [87]. Individual clusters are separated by the symbol /

Of the 25 groups that were problematic to distinguish using neighbour-joining trees and/ or distance methods, species within 13 groups could be correctly identified using diagnostic characters in CAOS (Table 2). Although CAOS distinguished the pacific black duck (Anas superciliosa) from the mallard (A. platyrhynchos) based on two nucleotides at positions 315 (A/G) and 402 (C/T), two mallard sequences had ambiguous calls at these positions (R and Y respectively) indicating the occurrence of heteroplasmy in these individuals. As such, these characters were not truly diagnostic and the species were considered to be indistinguishable. In total 169 out of 191 species with > 1 specimen (88.5%) could be successfully identified from their COI barcodes. Fifteen species had COI sequences that were difficult or impossible to distinguish from their respective closest relatives which do not occur in New Zealand highlighting the importance of thorough within genera sampling (Table 2).


DNA barcoding using the COI region has proven to be an effective tool for identifying New Zealand birds to species level, correctly identifying 88.5% of species represented by multiple specimens. Our success rate is slightly lower than other avian DNA barcoding studies which have reported upwards of a 93% success rate [23, 24, 26,27,28,29]. This is likely a reflection of our comprehensive dataset in which intraspecific variation was determined through the inclusion of conspecific individuals from throughout their world-wide distribution and through the inclusion of other closely related species that do not occur in New Zealand. The average intraspecific distance of 0.32% was slightly larger than the values reported for the avifauna of Scandinavia (0.24%), North America (0.23%), Argentina (0.24%) and the Netherlands (0.29%) [23, 25, 27, 28] though smaller than for the Japanese (0.46%) [29] and Turkish (0.62%) [30] avifaunas. This result is also likely a reflection of intraspecific sampling from a wide geographic distribution. While these earlier studies have used Kimura-2-Parameter (K2P) genetic distances, this does not affect comparisons as the average intraspecific K2P distance for this study is only 0.01% higher (0.33%) than the uncorrected p-distance. Nearest neighbour divergence varied from 0 to 13.27%, similar to the range found in eastern Palearctic birds by Kerr et al. [24]. Despite our best efforts, this is likely an inflated estimate due to the under sampling of some groups.

Evolutionary and systematic applications

Low levels of genetic divergence, particularly at just one locus, do not invalidate established taxonomy [45]. In cases of recent divergence, phenotypic differentiation can occur more rapidly than the complete sorting of mtDNA [45] while hybridisation and back-crossing can result in genetic introgression from one species to another [46]. Similar genetic patterns can also result from misidentification of specimens, although we made all efforts to minimise this issue. In this study, the majority of the 13 species pairs and triads that could not be distinguished by their COI barcodes, represent well-studied, valid species. For example, the species status of the extant New Zealand snipes (Coenocorypha spp.) are supported by reciprocal monophyly in both nuclear and mitochondrial markers as well as morphometric and plumage data [47]. Divergence is estimated to have occurred only about 96,000 years ago [47] suggesting incomplete lineage sorting as the most likely explanation for COI similarity. This is likely also the case for the masked gulls (Chroicocephalus scopulinus and C. bulleri) which diverged about 240,000 years ago [48]. Occasional hybridisation between these species has also been observed [49] and slow mutation rates have also been implicated [50]. Many gull species within the closely related genus Larus have indistinguishable COI barcodes e.g. [24, 27,28,29] which is attributed to recent speciation and hybridisation. The brown and south polar skua (Stercorarius antarcticus and S. maccormicki respectively) diverged only about 200,000 years ago and speciation is considered incomplete with hybridisation common [51, 52]. Originally considered a single species, the northern and southern giant petrels (Macronectes halli and M. giganteus) were split on the basis of morphological and behavioural differences [53]. This taxonomy is supported by nuclear and mitochondrial markers though genetic divergence levels are low [54], a reflection of recent divergence (about 200,000 years ago [54]) and hybridisation [55]. The low divergence between parakeets (Cyanoramphus spp.) is also likely a reflection of recent speciation [56]. Mitochondrial control region sequences were used to measure the divergence of the Antipodes parakeet (C. unicolor) which was estimated to have occurred ~ 270,000 years ago [56]. The sympatric Reischek’s parakeet (C. hochstetteri) colonised the Antipodes Islands much more recently, diverging from the red-crowned parakeet (C. novaezelandiae) ~ 100,000 years ago [56]. This is consistent with our finding that the Antipodes parakeet could be distinguished from the other two species by two diagnostic nucleotides. The introduced mallard and native pacific black duck are known to hybridise extensively with mtDNA introgression being bidirectional [57]. In two mallards, there was evidence of heteroplasmy at two nucleotide sites. Mitochondrial heteroplasmy, the occurrence of more than one haplotype within an individual, can occur as a result of mutation, recombination or paternal leakage [58]. The general assumption that mtDNA is uniparentally inherited and homoplasmic is being questioned by the accumulating evidence of paternal leakage in a variety of taxa (reviewed in Barr et al., [59]).

For other species that could not be identified by their COI barcodes, in-depth studies are lacking and further investigation is required. For example, the royal penguin (Eudyptes schlegeli) and macaroni penguin (E. chrysolophus) are considered conspecific by some [60, 61] and the mitochondrial hypervariable control region and the COI barcoding region show very low levels of divergence [62, 63]. Here we show that the COI sequences of specimens from royal and macaroni penguins generated by Baker et al. [63] could not be distinguished by CAOS. Taxonomic uncertainty also surrounds the New Zealand oystercatchers. A preliminary genetic study using mtDNA found no differences between the mainland species the South Island pied oystercatcher (Haematopus finschi) and the variable oystercatcher (H. unicolor) [64] which occasionally hybridise [65], but vary substantially in morphology. However, the Chatham Island species (H. chathamensis), which is considered by some to be a subspecies of the variable oystercatcher [66], was found to be distinct from the mainland species [64]. While none of the three species could be distinguished by COI barcodes, the Chatham Island oystercatcher showed the highest genetic divergence. Consistent plumage differences are currently the only basis for the separation of the Australasian bittern (Botaurus poiciloptilus) and the great bittern (B. stellaris) [61]. Our results suggest limited genetic divergence at the COI locus between the two species (mean distance of 0.12%) indicating their taxonomy may require further investigation.

Since we only have one Chatham Island pigeon (Hemiphaga chathamensis) specimen we can conclude little from the similarity between it and its sister taxa the New Zealand pigeon (H. novaeseelandiae). The Chatham Island pigeon was only recently elevated to species status on the basis of morphometric differences [67]. In a larger study, Goldberg et al. [68] found no differences in the COI region, and low divergence in cytochrome b and D-loop sequences (1.2 and 2.8% respectively) which was attributed to recent widespread dispersal. In the present study, the fulmar prion (Pachyptila crassirostris) was represented by a single sequence which showed low divergence when compared to the fairy prion (P. turtur). These species are sometimes considered conspecific [69].

Divergent COI lineages were evident within 16 species (Table 3). There is no level of genetic distance that can be used as a cut-off for species status, as speciation results in genetic divergence but is not caused by it [70]. However, divergence in COI barcodes can identify taxa in which further scrutiny may be required [24]. In seven globally distributed species, divergent lineages corresponded to recognised subspecies separated by large geographic distances; gentoo penguin (Pygoscelis papua), Wilson’s storm petrel (Oceanites oceanicus), great egret (Ardea alba), common pheasant (Phasianus colchicus), whimbrel (Numenius phaeopus), little owl (Athene noctua) and Eurasian skylark (Alauda arvensis). Within five of the six species that showed divergent lineages between New Zealand and Australia there are recognised subspecies; spur-winged plover (Vanellus miles), morepork (Ninox novaeseelandiae), little penguin (Eudyptula minor), New Zealand pipit (Anthus novaeseelandiae) and white-faced storm petrel (Pelagodroma marina). However, detailed geographic sampling would be required for each species to determine if COI barcodes could distinguish subspecies. Furthermore, it has been suggested that some of these currently recognised subspecies warrant separate species status. For example, there is evidence in the form of nuclear and mitochondrial markers [70, 71], behavioural [72] and plumage differences [71] that the New Zealand and Australian populations of little penguins (E. minor) should be recognised as separate species. Additionally, a recent study of the genus Ninox recommended that the mainland Australian population be treated as a separate species from Tasmanian and New Zealand populations [73].

The taxonomy of diving petrels remains unresolved and is the subject of debate [74, 75]. Within New Zealand, there was evidence of divergent lineages corresponding to recognised subspecies within four species. Populations of both little shearwater (Puffinus assimilis) and white-faced storm petrel (Pelagodroma marina) were divergent between the Mokohinau and Kermadec Islands. The white-faced storm petrels collected from the Mokohinau Islands and a beach wrecked individual found in the North Island (presumably from the Mokohinau Islands population) were over 5% divergent from the specimens collected in Australia and the Kermadec Islands. The Kermadec Island population is regarded as a distinct species (P. albiclunis) by Birds New Zealand [76]. Though we cannot be certain of which breeding population the Australian sample originated from, given its collection location it was probably a member of the Australian subspecies P. m. dulciae. Genetic comparisons of these populations are lacking, however, Silva et al. [77] found that the Mokohinau Island population was highly differentiated from North and South Atlantic populations using mitochondrial and nuclear markers. Rifleman (Acanthisitta chloris) and New Zealand robin (Petroica australis) showed divergent lineages between the North and South Islands. New Zealand robin lineages showed a divergence of 4.12%, similar to control region sequences which showed 5.9% divergence suggesting long-term isolation [78]. Indeed, Birds New Zealand recognises the North Island robin as a separate species P. longipes [76].

There was no clear pattern in the success rates of DNA barcoding despite the unique composition of the New Zealand avifauna. High levels of endemism had no obvious effect on success rates. For New Zealand species represented by > 1 specimen, 85.5% of endemic, 88.75% of native and 87.1% of introduced species could be identified by DNA barcodes. Divergent lineages were evident in a similar proportion of native and introduced species (6.9 and 9.7% respectively). The high prevalence of seabirds did not appear to influence success rates with 90.4% of seabirds successfully identified compared to 88% of land birds and 7.2% of seabirds showing evidence of divergent lineages compared to 9.2% of land birds. Importantly, this demonstrates that DNA barcoding can be successfully applied to species discrimination in fauna with a wide range of evolutionary patterns and life history traits.

Conservation management applications

The 928 COI sequences from 180 New Zealand bird species generated from this study form a substantial reference database that will be a valuable tool for specimen identification and the conservation of New Zealand birds. DNA barcoding has many advantages over morphological identification when applied to conservation management [14]. DNA barcoding can utilise non-invasive samples such as feathers or faeces [11] which is beneficial when dealing with rare and endangered birds or elusive predators [79]. Invasive mammalian predators are the biggest threat to the survival of New Zealand birds, responsible for the majority of the 26.6 million chick and egg losses of native bird species each year [80]. Diet studies using DNA barcodes can be used to assess predator impact on prey populations and provide superior detection and identification of prey species when compared to morphological analysis [14].

Performance of different methods of analysis

CAOS was found to be the most successful method for identifying specimens to the species level. All species that were distinguishable using neighbour-joining tree or other distance-based methods were also successfully identified by CAOS and an additional 14 species could only be identified using CAOS. While previously the application of CAOS has been limited by scalability issues [24], this has now been overcome and large datasets such as ours can be successfully analysed. We found that there were no differences in output when CAOS was run using smaller datasets consisting of species from one order. When sequences from species not found in our database were queried, CAOS correctly identified these individuals to the genus level. While this issue highlights the importance of thorough sampling in the reference database, genus level identification is more useful than no identification at all.

It is well-established that phylogenetic trees may perform poorly for the purpose of specimen identification [81, 82]. It is not possible to determine if a query sequence belongs to the species which it is topologically closest to unless it is nested within a monophyletic cluster [83]. Additionally, when either speciation is recent and individual genes are still incompletely sorted, or when introgressive hybridisation is occurring, non-monophyly is to be expected [84, 85]. Despite these limitations, quantifying the level of monophyly is still a useful descriptor of the data [4]. In this study, non-monophyly was observed in 14.7% of species, that is similar to values reported in other studies of Aves, between 10.4% [86] and 16.7% [46]. While distance and phylogenetic tree-based methods do not have the same level of success as CAOS, they reveal interesting features of the data which character-based methods do not. For example, evidence of divergent lineages can be quickly observed in a phylogenetic tree while large intraspecific variation may also indicate divergence. For species that show small interspecific distances and/or non-monophyly, we should be more cautious about identifications provided by CAOS as discussed above.


This study demonstrates that DNA barcoding can identify the majority of New Zealand birds to the species level. DNA barcoding has proved effective in ‘the Seabird Capital of the world’, a region where the unique avifauna has characteristics of both a continent and an island and is of mixed evolutionary origin. COI barcodes have highlighted species groups with limited divergence and other species that show evidence of divergent lineages that require further taxonomic scrutiny. Widespread geographic sampling means that the reported success rates are more conservative than they would have been had we only included specimens from the New Zealand region. The reference database generated by this study will provide a powerful tool for the conservation management of New Zealand birds.



We generated COI sequences from 928 specimens representing 180 species from the New Zealand region (Fig. 2a and b). Samples included voucher specimens from the Auckland War Memorial Museum, Museum Victoria, Museum of Natural Sciences at Louisiana State University and the Royal Ontario Museum. Other specimens were collected in the field by a large number of people over the last 35 years. Where possible, individual birds were sampled from across the species’ geographic range in order to determine levels of geographic variation. Additional Genbank sequences were also included in the analysis (see Additional file 4). Taxonomy was based upon Clements [87], including corrections and updates up to 7 March 2017.

Fig. 2

Map of New Zealand region as defined by this study including (a) New Zealand and its outlying islands and (b) the Ross Dependency, Antarctica

DNA sequencing

For the majority of samples, the DNA extraction protocol, PCR conditions, sequencing methodology and primer details were as previously described by Patel et al. [88]. For the remaining samples the methodology is outlined in Tavares and Baker [89]. Sequences shorter than 519 bp or which contained ten or more ambiguous base calls were excluded from analysis. Specimen information, sequences and trace files can be accessed on the Barcode of Life Data Systems website (BOLD, GenBank accession numbers MK261779 - MK262706) [90].

Additional data

The data gathered in this study were supplemented by 892 sequences from GenBank that fell into two categories. Firstly, 488 sequences from species that occur in the New Zealand region. For non-endemic species, specimens from across their geographic distribution were preferentially included, to capture the most geographic variation. Secondly, in instances where a species’ closest relative did not occur in New Zealand, sequences from the most closely related species available were included to increase within genera sampling. These additional 404 sequences are referred to as related species and while they were included in all analyses, only the success rates and divergence levels of New Zealand species are reported. A number of GenBank sequences followed outdated taxonomic classifications and were renamed to follow Clements [87] (see Additional file 4).

In total 1820 sequences were included in the analysis, of which 1416 were from 211 species that occur in the New Zealand region. A complete list of GenBank accession numbers of the sequences used in this study is available in Additional file 4.


Three DNA barcoding analysis methods were used; tree building, distance analysis and diagnostic character assignment. Tree building was conducted in MEGA version 7 [91]. Sequence alignment was performed with MUSCLE [92] and a neighbour-joining tree was produced based on uncorrected p-distances. P-distances have been shown to produce higher or similar levels of correct identification than Kimura-2-Parameter (K2P) distances which are commonly employed in barcoding studies [93, 94]. Support for monophyletic clades was measured using bootstrap values with 1000 replicates. Patterns of divergence were classified as either monophyletic with either greater than or less than 95% bootstrap support, paraphyletic or polyphyletic.

Local barcode gap analysis was conducted by calculating maximum intraspecific and minimum interspecific (nearest neighbour) genetic distances for each species using the Spider package [95] for RStudio [96]. For each species, these values were plotted against each other to visualise ‘local’ barcoding gaps; discontinuity between levels of intraspecific and interspecific distances [4]. Nearest neighbour distances were used in preference to average interspecific distances because species identification is ultimately dependent upon how different a sequence is from its closest allospecific sequence, as opposed to the distance to the “average” sequence [97]. An optimised global distance threshold was also calculated from the data, minimising the cumulative error rate [8].

Character-based identification was implemented in CAOS [42,43,44]. CAOS identifies diagnostic characters, termed ‘character attributes’ (CA’s) from a tree of pre-defined species. Single CA’s may be either pure (sPu’s) if they are shared by all members of a clade and are absent from the other clades or private (sPr’s) if they are shared only by some members of a clade [7]. Detailed methodology can be found in the Additional file 5. In brief, CAOS barcoding is comprised of three steps, each performed by a separate program. Firstly, the CAOS-Analyzer extracts CA’s from the input nexus file that consists of the sequence alignment and tree file fused together. Next, the outputs of the CAOS-Analyzer are converted into an easily interpretable character-based barcode matrix using the CAOS-Barcoder. Finally, the CAOS-Classifier tests the efficacy of this matrix by attempting to assign a new query specimen to the correct species in the reference dataset. For species with multiple representatives, the shortest sequence was excluded from the reference database and used as a query sequence.



Characteristic Attribute Organization System


Cytochrome c oxidase subunit I


mitochondrial DNA


  1. 1.

    Hebert PDN, Cywinska A, Ball SL, DeWaard JR. Biological identifications through DNA barcodes. Proc R Soc B Biol Sci. 2003;270:313–21.

    CAS  Article  Google Scholar 

  2. 2.

    Taylor HR, Harris WE. An emergent science on the brink of irrelevance: a review of the past 8 years of DNA barcoding. Mol Ecol Resour. 2012;12:377–88.

    CAS  PubMed  Article  Google Scholar 

  3. 3.

    DeSalle R. Species discovery versus species identification in DNA barcoding efforts: response to Rubinoff. Conserv Biol. 2006;20:1545–7.

    PubMed  Article  Google Scholar 

  4. 4.

    Collins RA, Cruickshank RH. The seven deadly sins of DNA barcoding. Mol Ecol Resour. 2012;13:969–75.

    PubMed  Google Scholar 

  5. 5.

    Dunn CP. Keeping taxonomy based in morphology. Trends Ecol Evol. 2003;18:269–70.

    Article  Google Scholar 

  6. 6.

    Seberg O, Humphries CJ, Knapp S, Stevenson DW, Petersen G, Scharff N, et al. Shortcuts in systematics? A commentary on DNA-based taxonomy. Trends Ecol Evol. 2003;18:63–5.

    Article  Google Scholar 

  7. 7.

    Rach J, DeSalle R, Sarkar I, Schierwater B, Hadrys H. Character-based DNA barcoding allows discrimination of genera, species and populations in Odonata. Proc R Soc B Biol Sci. 2008;275:237–47.

    CAS  Article  Google Scholar 

  8. 8.

    Meyer CP, Paulay G. DNA barcoding: error rates based on comprehensive sampling. PLoS Biol. 2005;3:e422.

    PubMed  PubMed Central  Article  CAS  Google Scholar 

  9. 9.

    Price BW, Henry CS, Hall AC, Mochizuki A, Duelli P, Brooks SJ. Singing from the grave: DNA from a 180 year old type specimen confirms the identity of Chrysoperla carnea (Stephens). PLoS One. 2015;10:e0121127.

    PubMed  PubMed Central  Article  CAS  Google Scholar 

  10. 10.

    Areta JI, Piacentini VDQ, Haring E, Gamauf A, Silveira LF, Machado E, et al. Tiny bird, huge mystery—the possibly extinct hooded seedeater (Sporophila melanops) is a Capuchino with a melanistic cap. PLoS One. 2016;11:e0154231.

    PubMed  PubMed Central  Article  CAS  Google Scholar 

  11. 11.

    Li J, Cui Y, Jiang J, Yu J, Niu L, Deng J, et al. Applying DNA barcoding to conservation practice: a case study of endangered birds and large mammals in China. Biodivers Conserv. 2017;26:653–68.

    Article  Google Scholar 

  12. 12.

    Dalton DL, Kotze A. DNA barcoding as a tool for species identification in three forensic wildlife cases in South Africa. Forensic Sci Int. 2011;207:e51–4.

    CAS  PubMed  Article  Google Scholar 

  13. 13.

    Ardura A, Linde AR, Moreira JC, Garcia-Vazquez E. DNA barcoding for conservation and management of Amazonian commercial fish. Biol Conserv. 2010;143:1438–43.

    Article  Google Scholar 

  14. 14.

    Zarzoso-Lacoste D, Bonnaud E, Corse E, Gilles A, Meglecz E, Costedoat C, et al. Improving morphological diet studies with molecular ecology: an application for invasive mammal predation on island birds. Biol Conserv. 2016;193:134–42.

    Article  Google Scholar 

  15. 15.

    Méheust E, Alfonsi E, Le Ménec P, Hassani S, Jung JL. DNA barcoding for the identification of soft remains of prey in the stomach contents of grey seals (Halichoerus grypus) and harbour porpoises (Phocoena phocoena). Mar Biol Res. 2015;11:385–95.

  16. 16.

    Waugh J, Evans MW, Millar CD, Lambert DM. Birdstrikes and barcoding: can DNA methods help make the airways safer? Mol Ecol Resour. 2011;11:38–45.

    PubMed  Article  Google Scholar 

  17. 17.

    Dove CJ, Rotzel NC, Heacker M, Weigt LA. Using DNA barcodes to identify bird species involved in birdstrikes. J Wildl Manag. 2008;72:1231–6.

    Article  Google Scholar 

  18. 18.

    Krishnamurthy PK, Francis RA. A critical review on the utility of DNA barcoding in biodiversity conservation. Biodivers Conserv. 2012;21:1901–19.

  19. 19.

    Hajibabaei M, Smith MA, Janzen DH, Rodriguez JJ, Whitfield JB, Hebert PDN. A minimalist barcode can identify a specimen whose DNA is degraded. Mol Ecol Notes. 2006;6:959–64.

    CAS  Article  Google Scholar 

  20. 20.

    Barreira AS, Lijtmaer DA, Tubaro PL. The multiple applications of DNA barcodes in avian evolutionary studies. Genome. 2016;13 October:1–13.

  21. 21.

    Hebert PDN, Stoeckle MY, Zemlak TS, Francis CM. Identification of birds through DNA barcodes. PLoS Biol. 2004;2:e312.

    PubMed  PubMed Central  Article  CAS  Google Scholar 

  22. 22.

    DeSalle R, Egan MG, Siddall M. The unholy trinity: taxonomy, species delimitation and DNA barcoding. Philos Trans R Soc B Biol Sci. 2005;360:1905–16.

    CAS  Article  Google Scholar 

  23. 23.

    Kerr KCR, Stoeckle MY, Dove CJ, Weigt LA, Francis CM, Hebert PDN. Comprehensive DNA barcode coverage of north American birds. Mol Ecol Notes. 2007;7:535–43.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  24. 24.

    Kerr KC, Birks SM, Kalyakin MV, Red’kin YA, Koblik EA, Hebert PD. Filling the gap - COI barcode resolution in eastern Palearctic birds. Front Zool. 2009;6:29.

    PubMed  PubMed Central  Article  CAS  Google Scholar 

  25. 25.

    Kerr KCR, Lijtmaer DA, Barreira AS, Hebert PDN, Tubaro PL. Probing evolutionary patterns in Neotropical birds through DNA barcodes. PLoS One. 2009;4:e4379.

    PubMed  PubMed Central  Article  CAS  Google Scholar 

  26. 26.

    Tavares ES, Gonçalves P, Miyaki CY, Baker AJ. DNA barcode detects high genetic structure within Neotropical bird species. PLoS One. 2011;6:e28543.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  27. 27.

    Johnsen A, Rindal E, Ericson PGP, Zuccon D, Kerr KCR, Stoeckle MY, et al. DNA barcoding of Scandinavian birds reveals divergent lineages in trans-Atlantic species. J Ornithol. 2010;151:565–78.

    Article  Google Scholar 

  28. 28.

    Aliabadian M, Beentjes K, Roselaar K, van Brandwijk H, Nijman V, Vonk R. DNA barcoding of Dutch birds. ZooKeys. 2013;365:25–48.

    Article  Google Scholar 

  29. 29.

    Saitoh T, Sugita N, Someya S, Iwami Y, Kobayashi S, Kamigaichi H, et al. DNA barcoding reveals 24 distinct lineages as cryptic bird species candidates in and around the Japanese archipelago. Mol Ecol Resour. 2015;15:177–86.

    CAS  PubMed  Article  Google Scholar 

  30. 30.

    Bilgin R, Ebeoğlu N, İnak S, Kırpık MA, Horns JJ, Şekercioğlu ÇH. DNA barcoding of birds at a migratory hotspot in eastern Turkey highlights continental phylogeographic relationships. PLoS One. 2016;11:e0154454.

    PubMed  PubMed Central  Article  CAS  Google Scholar 

  31. 31.

    McLoughlin S. The breakup history of Gondwana and its impact on pre-Cenozoic floristic provincialism. Aust J Bot. 2001;49.

  32. 32.

    Daugherty CH, Gibbs GW, Hitchmough RA. Mega-Island or micro-continent? New Zealand and its Fauna. Trends Ecol Evol. 1993;8:437–42.

    CAS  PubMed  Article  Google Scholar 

  33. 33.

    Goldberg J, Trewick SA, Paterson AM. Evolution of New Zealand’s terrestrial fauna: a review of molecular evidence. Philos Trans R Soc B Biol Sci. 2008;363:3319–34.

    Article  Google Scholar 

  34. 34.

    Atkinson IAE, Taylor RH. Distribution of alien mammals on New Zealand islands. Wellington: DSIR Land Resources Report No. 91/50; 1991.

  35. 35.

    Trewick S, Gibb G. Vicars, tramps and assembly of the New Zealand avifauna: a review of molecular phylogenetic evidence. Ibis. 2010;152:226–53.

    Article  Google Scholar 

  36. 36.

    Robertson HA, Baird K, Dowding JE, Elliott GP, Hitchmough RA, Miskelly CM, et al. Conservation status of New Zealand birds, 2016. Wellington: New Zealand Threat Classification Series 19; 2017

  37. 37.

    Fleming CA. New Zealand as a minor source of terrestrial plants and animals in the Pacific. Tuatara: Journal of the Biological Society 1976;22.

  38. 38.

    Taylor GA. Action plan for seabird conservation in New Zealand. Part A, Threatened seabirds. Threatened species occasional publication. Department of Conservation; 2000.

  39. 39.

    Jones HP, Holmes ND, Butchart SHM, Tershy BR, Kappes PJ, Corkery I, et al. Invasive mammal eradication on islands results in substantial conservation gains. Proc Natl Acad Sci. 2016;113:4033–8.

    CAS  PubMed  Article  Google Scholar 

  40. 40.

    Cowie RH, Holland BS. Dispersal is fundamental to biogeography and the evolution of biodiversity on oceanic islands. J Biogeogr. 2006;33:193–8.

    Article  Google Scholar 

  41. 41.

    Schippers P, Snep RPH, Schotman AGM, Jochem R, Stienen EWM, Slim PA. Seabird metapopulations: searching for alternative breeding habitats. Popul Ecol. 2009;51:459–70.

    Article  Google Scholar 

  42. 42.

    Sarkar NI, Thornton JW, Planet PJ, Figurski DH, Schierwater B, DeSalle R. An automated phylogenetic key for classifying homeoboxes. Mol Phylogenet Evol. 2002;24:388–99.

  43. 43.

    Sarkar IN, Planet PJ, Desalle R. CAOS software for use in character-based DNA barcoding. Mol Ecol Resour. 2008;8:1256–9.

    CAS  PubMed  Article  Google Scholar 

  44. 44.

    Bergmann T, Hadrys H, Breves G. Character-based DNA barcoding: a superior tool for species classification. Berl Münch Tierärztl. 2009;12:446–50.

    Google Scholar 

  45. 45.

    Joseph L, Omland KE. Phylogeography: its development and impact in Australo-Papuan ornithology with special reference to paraphyly in Australian birds. Emu Austral Ornithol. 2009;109:1–23.

    Article  Google Scholar 

  46. 46.

    Funk DJ, Omland KE. Species-leve paraphyly and polyphyly: frequency, causes, and consequences, with insights from animal mitochondrial DNA. Annu Rev Ecol Evol Syst. 2003;34:397–423.

    Article  Google Scholar 

  47. 47.

    Baker AJ, Miskelly CM, Haddrath O. Species limits and population differentiation in New Zealand snipes (Scolopacidae: Coenocorypha). Conserv Genet. 2010;11:1363–74.

    Article  Google Scholar 

  48. 48.

    Given AD, Mills JA, Baker AJ. Molecular evidence for recent radiation in southern hemisphere masked gulls. Auk. 2005;122:268–79.

    Article  Google Scholar 

  49. 49.

    Gurr L. Inbreeding of Larus novaehollandiae scopulinus and Larus bulleri in the wild in New Zealand. Ibis. 1966;109:552–5.

    Article  Google Scholar 

  50. 50.

    Crochet P-A, Desmarais E. Slow rate of evolution in the mitochondrial control region of gulls (Aves: Laridae). Mol Biol Evol. 2000;17:1797–806.

    CAS  PubMed  Article  Google Scholar 

  51. 51.

    Ritz MS, Millar C, Miller GD, Phillips RA, Ryan P, Sternkopf V, et al. Phylogeography of the southern skua complex; rapid colonization of the southern hemisphere during a glacial period and reticulate evolution. Mol Phylogenet Evol. 2008;49:292–303.

    CAS  PubMed  Article  Google Scholar 

  52. 52.

    Ritz MS, Hahn S, Janicke T, Peter H-U. Hybridisation between south polar skua (Catharacta maccormicki) and Brown skua (C. antarctica lonnbergi) in the Antarctic peninsula region. Polar Biol. 2006;29:153–9.

    Article  Google Scholar 

  53. 53.

    Bourne W, Warham J. Geographical variation in the giant petrels of the genus Macronectes. Ardea. 1966;54:45–67.

    Google Scholar 

  54. 54.

    Techow N, O’Ryan C, Phillips R, Gales R, Marin M, Patterson-Fraser D, et al. Speciation and phylogeography of giant petrels Macronectes. Mol Phylogenet Evol. 2010;54:472–87.

  55. 55.

    Brown RM, Techow NMSM, Wood AG, Phillips RA. Hybridization and Back-crossing in Giant petrels (Macronectes giganteus and M. halli) at Bird Island, South Georgia, and a summary of hybridization in seabirds. PLoS One. 2015;10:e0121688.

    PubMed  PubMed Central  Article  CAS  Google Scholar 

  56. 56.

    Boon WM, Kearvell JC, Daugherty CH, Chambers GK. Molecular systematics and conservation of kakariki (Cyanorampbus spp.). Wellington: Sciences for Conservation 176; 2001.

  57. 57.

    Rhymer JM, Williams MJ, Braun MJ. Mitchondrial analysis of gene flow between New Zealand mallards (Anas platyrhynchos) and Grey ducks (A. superciliosa). Auk. 1994;111:970–8.

    Article  Google Scholar 

  58. 58.

    Kvist L, Martens J, Nazarenko AA, Orell M. Paternal leakage of mitochondrial DNA in the great tit (Parus major). Mol Biol Evol. 2003;20:243–7.

    CAS  PubMed  Article  Google Scholar 

  59. 59.

    Barr CM, Neiman M, Taylor DR. Inheritance and recombination of mitochondrial genomes in plants, fungi and animals. New Phytol. 2005;168:39–50.

    CAS  PubMed  Article  Google Scholar 

  60. 60.

    Falla RA, Mougin RL. Sphenisciformes. In: Mayr E, Cottrell GW, editors. Checklist of the Birds of the World. 2nd ed. Cambridge: Mass., Mass.; 1979. p. 121–34.

    Google Scholar 

  61. 61.

    Christidis L, Boles WE. Systematics and taxonomy of Australian birds. Collingwood: CSIRO Publishing; 2008.

    Book  Google Scholar 

  62. 62.

    Ritchie P. The evolution of the mitochondrial DNA control region in the Adelie penguins of Antarctica. PhD thesis. New Zealand: Massey University; 2001.

  63. 63.

    Baker AJ, Tavares ES, Elbourne RF. Countering criticisms of single mitochondrial DNA gene barcoding in birds. Mol Ecol Resour. 2009;9(SUPPL. 1):257–68.

    CAS  PubMed  Article  Google Scholar 

  64. 64.

    Banks JC, Paterson AM. A preliminary study of the genetic differences in New Zealand oystercatcher species. N Z J Zool. 2007;34:141–4.

    Article  Google Scholar 

  65. 65.

    Crocker T, Petch S, Sagar P. Hybridisation by South Island pied oystercatcher (Haematopus finschi) and variable oystercatcher (H. unicolor) in Canterbury. Notornis. 2010;57:27–32.

    Google Scholar 

  66. 66.

    Dickinson EC. The Howard and Moore complete checklist of birds of the world. London: Christopher Helm; 2003.

    Google Scholar 

  67. 67.

    Millener PR, Powlesland RG. The Chatham Islands pigeon (Parea) deserves full species status; Hemiphaga chathamensis (Rothschild 1891); Aves: Columbidae. J R Soc N Z. 2001;31:365–83.

    Article  Google Scholar 

  68. 68.

    Goldberg J, Trewick SA, Powlesland RG. Population structure and biogeography of Hemiphaga pigeons (Aves: Columbidae) on islands in the New Zealand region. J Biogeogr. 2011;38:285–98.

    Article  Google Scholar 

  69. 69.

    Cox JB. Some remarks on the breeding distribution and taxonomy of the prions (Procellariidae: Pachyptila). Rec South Aust Mus. 1980;18:91–121.

  70. 70.

    Waugh J. DNA barcodes highlight two clusters within the little penguin (Eudyptula minor): time to reassess species delineation? Notornis. 2016;63:66–72.

    Google Scholar 

  71. 71.

    Grosser S, Burridge CP, Peucker AJ, Waters JM. Coalescent modelling suggests recent secondary-contact of cryptic penguin species. PLoS One. 2015;10:e0144966.

    PubMed  PubMed Central  Article  CAS  Google Scholar 

  72. 72.

    Daniel TA, Chiaradia A, Logan M, Quinn GP, Reina RD. Synchronized group association in little penguins, Eudyptula minor. Anim Behav. 2007;74:1241–8.

    Article  Google Scholar 

  73. 73.

    Gwee CY, Christidis L, Eaton JA, Norman JA, Trainor CR, Verbelen P, et al. Bioacoustic and multi-locus DNA data of Ninox owls support high incidence of extinction and recolonisation on small, low-lying islands across Wallacea. Mol Phylogenet Evol. 2017;109:246–58.

    CAS  PubMed  Article  Google Scholar 

  74. 74.

    Scofield P, Stephenson B. Birds of New Zealand: a photographic guide. Auckland: Auckland University Press; 2013.

    Google Scholar 

  75. 75.

    Fischer JH, Debski I, Miskelly CM, Bost CA, Fromant A, Tennyson AJD, et al. Analyses of phenotypic differentiations among South Georgian Diving Petrel (Pelecanoides georgicus) populations reveal an undescribed and highly endangered species from New Zealand. PLoS ONE. 2018;13:e0197766.

  76. 76.

    Gill B. Ornithological Society of new Zealand, checklist committee. Checklist of the birds of New Zealand, Norfolk and Macquarie Islands, and the Ross dependency, Antarctica. 4th ed. Wellington: Te Papa Press in association with the Ornithological Society of New Zealand; 2010.

    Google Scholar 

  77. 77.

    Silva MC, Matias R, Wanless RM, Ryan PG, Stephenson BM, Bolton M, et al. Understanding the mechanisms of antitropical divergence in the seabird white-faced storm-petrel (Procellariiformes: Pelagodroma marina) using a multilocus approach. Mol Ecol. 2015;24:3122–37.

    PubMed  Article  Google Scholar 

  78. 78.

    Miller HC, Lambert DM. A molecular phylogeny of New Zealand’s Petroica (Aves: Petroicidae) species based on mitochondrial DNA sequences. Mol Phylogenet Evol. 2006;40:844–55.

    CAS  PubMed  Article  Google Scholar 

  79. 79.

    Waits LP, Paetkau D. Noninvasive genetic sampling tools for wildlife biologists: a review of applications and recommendations for accurate data collection. J Wildl Manag. 2005;69:1419–33.

    Article  Google Scholar 

  80. 80.

    Russell JC, Innes JG, Brown PH, Byrom AE. Predator-free New Zealand: conservation country. Bioscience. 2015;65:520–5.

    PubMed  PubMed Central  Article  Google Scholar 

  81. 81.

    Meier R, Shiyang K, Vaidya G, Ng P. DNA barcoding and taxonomy in Diptera: a tale of high intraspecific variability and low identification success. Syst Biol. 2006;55:715–28.

    PubMed  Article  Google Scholar 

  82. 82.

    Zhang A-B, Muster C, Liang H-B, Zhu C-D, Crozier R, Wan P, et al. A fuzzy-set-theory-based approach to analyse species membership in DNA barcoding. Mol Ecol. 2012;21:1848–63.

    CAS  PubMed  Article  Google Scholar 

  83. 83.

    Will KW, Mishler BD, Wheeler QD. The perils of DNA barcoding and the need for integrative taxonomy. Syst Biol. 2005;54:844–51.

    PubMed  Article  Google Scholar 

  84. 84.

    Avise JC. Phylogeography: the history and formation of species. Cambridge (Massachusetts): Harvard University Press; 2000.

    Google Scholar 

  85. 85.

    McKay BD, Zink RM. The causes of mitochondrial DNA gene tree paraphyly in birds. Mol Phylogenet Evol. 2010;54:647–50.

    PubMed  Article  Google Scholar 

  86. 86.

    Chaves BRN, Chaves AV, Nascimento ACA, Chevitarese J, Vasconcelos MF, Santos FR. Barcoding Neotropical birds: assessing the impact of nonmonophyly in a highly diverse group. Mol Ecol Resour. 2015;15:921–31.

    PubMed  Article  Google Scholar 

  87. 87.

    Clements J. The Clements checklist of the birds of the world. 6th ed. New York: Cornell University Press; 2007.

    Google Scholar 

  88. 88.

    Patel S, Waugh J, Millar CD, Lambert DM. Conserved primers for DNA barcoding historical and modern samples from New Zealand and Antarctic birds. Mol Ecol Resour. 2010;10:431–8.

    CAS  PubMed  Article  Google Scholar 

  89. 89.

    Tavares ES, Baker AJ. Single mitochondrial gene barcodes reliably identify sister-species in diverse clades of birds. BMC Evol Biol. 2008;8:81.

    PubMed  PubMed Central  Article  CAS  Google Scholar 

  90. 90.

    Ratnasingham S, Hebert PDN. BOLD: the barcode of life data system ( Mol Ecol Notes 2007;7:355–64.

  91. 91.

    Kumar S, Stecher G, Tamura K. MEGA7: molecular evolutionary genetics analysis version 7.0 for bigger datasets. Mol Biol Evol. 2016;33:1870–4.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  92. 92.

    Edgar RC. MUSCLE: a multiple sequence alignment method with reduced time and space complexity. BMC Bioinformatics. 2004;5:113.

    PubMed  PubMed Central  Article  CAS  Google Scholar 

  93. 93.

    Srivathsan A, Meier R. On the inappropriate use of Kimura-2-parameter (K2P) divergences in the DNA-barcoding literature. Cladistics. 2012;28:190–4.

    Article  Google Scholar 

  94. 94.

    Collins RA, Boykin LM, Cruickshank RH, Armstrong KF. Barcoding’s next top model: an evaluation of nucleotide substitution models for specimen identification. Methods Ecol Evol. 2012;3:457–65.

    Article  Google Scholar 

  95. 95.

    Brown S, Collins R, Boyer S, Lefort M-C, Malumbres-Olarte J, Vink C, et al. Spider: an R package for the analysis of species identity and evolution, with particular reference to DNA barcoding. Mol Ecol Resour. 2012;12:562–5.

    PubMed  Article  Google Scholar 

  96. 96.

    RStudio Team. RStudio: integrated development environment for R. 2016.

  97. 97.

    Meier R, Zhang G, Ali F. The use of mean instead of smallest interspecific distances exaggerates the size of the “barcoding gap” and leads to misidentification. Syst Biol. 2008;57:809–13.

    PubMed  Article  Google Scholar 

  98. 98.

    Baker AJ, Pereira SL, Haddrath OP, Edge K-A. Multiple gene evidence for expansion of extant penguins out of Antarctica due to global cooling. Proc R Soc B Biol Sci. 2006;273:11–7.

    Article  Google Scholar 

  99. 99.

    Borboroglu PG, Boersma PD. Penguins: natural history and conservation. Seattle: University of Washington Press; 2013.

    Google Scholar 

  100. 100.

    Gonzalez J, Düttmann H, Wink M. Phylogenetic relationships based on two mitochondrial genes and hybridization patterns in Anatidae. J Zool. 2009;279:310–8.

    Article  Google Scholar 

  101. 101.

    Sun Z, Pan T, Hu C, Sun L, Ding H, Wang H, et al. Rapid and recent diversification patterns in Anseriformes birds: inferred from molecular phylogeny and diversification analyses. PLoS One. 2017;12:e0184529.

    PubMed  PubMed Central  Article  CAS  Google Scholar 

  102. 102.

    Alström P, Barnes KN, Olsson U, Barker FK, Bloomer P, Khan AA, et al. Multilocus phylogeny of the avian family Alaudidae (larks) reveals complex morphological evolution, non-monophyletic genera and hidden species diversity. Mol Phylogenet Evol. 2013;69:1043–56.

    PubMed  Article  Google Scholar 

  103. 103.

    Irwin DE, Rubtsov AS, Panov EN. Mitochondrial introgression and replacement between yellowhammers (Emberiza citrinella) and pine buntings (Emberiza leucocephalos) (Aves: Passeriformes). Biol J Linn Soc. 2009;98:422–38.

    Article  Google Scholar 

  104. 104.

    Marthinsen G, Wennerberg L, Lifjeld JT. Low support for separate species within the redpoll complex (Carduelis flammea–hornemanni–cabaret) from analyses of mtDNA and microsatellite markers. Mol Phylogenet Evol. 2008;47:1005–17.

    CAS  PubMed  Article  Google Scholar 

  105. 105.

    Peucker AJ, Dann P, Burridge CP. Range-wide Phylogeography of the little penguin (Eudyptula minor): evidence of long-distance dispersal. Auk. 2009;126:397–408.

    Article  Google Scholar 

  106. 106.

    Paterson AM, Wallis LJ, Wallis GP. Preliminary molecular analysis of Pelecanoides georgicus (Procellariiformes: Pelecanoididae) on Whenua Hou (Codfish Island): implications for its taxonomic status. N Z J Zool. 2000;27:415–23.

  107. 107.

    Austin JJ, Bretagnolle V, Pasquet E. A global molecular phylogeny of the small Puffinus shearwaters and implications for systematics of the little-Audubon’s shearwater complex. Auk. 2004;121:647–864.

    Article  Google Scholar 

  108. 108.

    Pellegrino I, Negri A, Cucco M, Mucci N, Pavia M, Šálek M, et al. Phylogeography and Pleistocene refugia of the little owl Athene noctua inferred from mtDNA sequence data. Ibis. 2014;156:639–57.

    Article  Google Scholar 

  109. 109.

    Marchant S, Higgins PJ. Handbook of Australian, New Zealand and Antarctic birds. Volume 5: tyrant-flycatchers to chats. Melbourne: Oxford University Press; 2001.

    Google Scholar 

  110. 110.

    Schodde R, Mason I. The directory of Australian birds: a taxonomic and zoogeographic atlas of the biodiversity of birds in Australia and its territories. Victoria: CSIRO Publishing; 1999.

    Google Scholar 

Download references


This project would not have been possible without the support of a large number of people who very generously provided samples. We wish to thank the following people and institutions for providing samples: Murray Potter (Massey University), Bruce Robertson (University of Otago), Todd Landers (Auckland Council), James Russell (University of Auckland), Matt Rayner (Auckland War Memorial Museum), Gerry Kooyman (Scripps Institution of Oceanography), Emma Marks (University of Auckland), Graeme Taylor (New Zealand Department of Conservation), Kate McInnes (New Zealand Department of Conservation), Mark Hauber (University of Illinois), Steffanie Ismar (Helmholtz Centre for Ocean Research), Murray Williams (Victoria University of Wellington), Steven Lawrence, Colin Miskelly (Museum of New Zealand Te Papa Tongarewa), Sara Treadgold (New Zealand Department of Conservation), Mick Clout (University of Auckland), Bill Peacock (Northcote College), Shaun O’Connor (New Zealand Department of Conservation) Lawson Davey (Nelson Marlborough Fish and Game), Sandra Anderson (University of Auckland), Daryl Eason (New Zealand Department of Conservation), Andrew Fidler (Cawthron Institute), Rosemary Barraclough (University of Auckland), Shinichi Nakagawa (University of New South Wales), Michael Anderson (Massey University), Erica Sommer (New Zealand Department of Conservation), Jo Hiscock (New Zealand Department of Conservation), Stuart Cockburn (New Zealand Department of Conservation), Peter Moore (New Zealand Department of Conservation), David Thompson (National Institute of Water and Atmospheric Research), the late Ian Jamieson (University of Otago), Sylvia Durant (Bird Rescue) and Rosemary Tulley (Bird Rescue).


This research was funded by the University of Auckland’s Vice-Chancellor’s University Development Fund; the University of Auckland Faculty Research Development Funds; the Allan Wilson Centre for Molecular Ecology and Evolution; the Canadian Barcode of Life Network from Genome Canada through the Ontario Genomics Institute, the Natural Sciences and Engineering Research Council of Canada; and the Royal Ontario Museum Governors’ Fund. Tjard Bergmann was financed by the H. Wilhelm Schaumann Stiftung.

Availability of data and materials

Specimen information, sequences and trace files can be accessed on the Barcode of Life Data Systems website (BOLD, GenBank accession numbers MK261779 - MK262706).

Author information




CM, DL, AB, LC, JN and BG conceived the study and participated in its design and coordination. SP, JW, ET, JT and OH generated the sequence data and curated the database. JT and TB conducted the analysis. JT, SP, CD, JW, PS and JN interpreted the data. JT, SP and CM drafted the manuscript. All authors contributed to revisions and approved the final manuscript.

Corresponding author

Correspondence to Craig Millar.

Ethics declarations

Ethics approval and consent to participate

Many of the samples in this study were obtained from dead specimens held in museums and other institutions for which no ethics approval is required. The remainder of the samples were collected by the New Zealand Department of Conservation under their animal ethics committee approval.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Additional files

Additional file 1:

Frequency distribution of maximum intraspecific and minimum interspecific genetic distances measured using a standardised 648 bp region of the cytochrome c oxidase gene for all New Zealand bird species with > 1 specimen obtained during the study. The dashed line indicates the calculated optimised distance threshold (0.025%). (DOCX 96 kb)

Additional file 2:

Cumulative error plot of type I (false positive) and type II (false negative) errors for different divergence thresholds of maximum intraspecific and minimum interspecific genetic distances measured using a standardised 648 bp region of the cytochrome c oxidase gene for all New Zealand bird species with > 1 specimen obtained during the study. The optimal threshold occurs at 0.25%. (DOCX 171 kb)

Additional file 3:

Neighbour Joining tree of sequences of a standardised 648 bp region of the cytochrome c oxidase gene obtained from New Zealand and closely related bird species in this study. Bootstrap support values ≥0.5 are indicated. Monophyletic clades have been collapsed. Branches are coloured by Order. (PDF 6097 kb)

Additional file 4:

List of all sequences of a standardised 648 bp region of the cytochrome c oxidase gene obtained from New Zealand and closely related bird species used in this study including Genbank accession numbers. (DOCX 63 kb)

Additional file 5:

Detailed methodology of CAOS analysis. (DOCX 29 kb)

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (, which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Tizard, J., Patel, S., Waugh, J. et al. DNA barcoding a unique avifauna: an important tool for evolution, systematics and conservation. BMC Evol Biol 19, 52 (2019).

Download citation


  • New Zealand birds
  • Cytochrome c oxidase subunit I
  • COI
  • Specimen identification
  • Conservation
  • DNA barcodes