Biodiversity Data Journal 9: 64842 CO) doi: 10.3897/BDJ.9.e64842 open access Research Article DNA barcodes of birds from northern Colombia Paulo Cesar Pulgarin-R*, Martha Olivera-Angel$, Luisa Ortiz8, Duvan Nanclares§, Sara Velasquez-Restrepo!, Juan Fernando Diaz-Nieto! + Facultad de Ciencias y Biotecnologia, Universidad CES, Medellin, Colombia § Biogénesis, Facultad de Ciencias Agrarias, Universidad de Antioquia, Cl. 73 #73A-79, Medellin, Colombia | Grupo Biodiversidad, Evolucién y Conservacién (BEC), Departamento de Ciencias Biolégicas, Escuela de Ciencias, Universidad EAFIT, Carrera 49 No. 7 sur-50, Medellin, Colombia Corresponding author: Juan Fernando Diaz-Nieto (jdiazni@eafit.edu.co) Academic editor: Caio J. Carlos Received: 22 Feb 2021 | Accepted: 02 May 2021 | Published: 21 May 2021 Citation: Pulgarin-R PC, Olivera-Angel M, Ortiz L, Nanclares D, Velasquez-Restrepo S, Diaz-Nieto JF (2021) DNA barcodes of birds from northern Colombia. Biodiversity Data Journal 9: e64842. https://doi.org/10.3897/BDJ.9.e64842 Abstract DNA barcode datasets are a useful tool for conservation and aid in taxonomic identification, particularly in megadiverse tropical countries seeking to document and describe its biota, which is dropping at an alarming rate during recent decades. Here we report the barcodes for several low elevation bird species from northern Colombia with the goal to provide tools for species identification in this region of South America. We blood- sampled birds in a lowland tropical forest with various degrees of intervention using standard 3 x 12 m mist-nets. We extracted DNA and sequenced the COI barcode gene using standard primers and laboratory methods. We obtained 26 COI sequences from 18 species, 10 families and three orders and found that barcodes largely matched (but not always) phenotypic identification (> 90%) and they also facilitated the identification of several challenging passerine species. Despite our reduced sampling, our study represents the first attempt to document COI barcodes for birds (from blood samples) in this part of Colombia, which fills a considerable gap of sampling in this part of South America. Keywords aves, lowland tropical forest, mtDNA, northern Colombia © Pulgarin-R P et al. This is an open access article distributed under the terms of the Creative Commons Attribution License (CC BY 4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited. 2 Pulgarin-R P et al Introduction DNA barcode reference libraries are a useful tool for conservation and aid in taxonomic identification (Gonzalez et al. 2009, Waugh 2007) for many biological groups (Hebert et al. 2003). Megadiverse countries, such as Colombia, are desperately in need of documenting and describing its biota, which is declining at an alarming rate during recent decades (Shaw et al. 2013), with emphasis on the putative cryptic diversity present in tropical areas (Crawford et al. 2012, Lohman et al. 2010, Stefan et al. 2018). Despite efforts to encourage sequence data collection and sharing through local and global initiatives (e.g. Barcode Life Data System, BOLD), most taxa are under-represented for barcodes (de Kerdrel et al. 2020, Gaytan et al. 2020, Ko et al. 2013). Naturally, although barcodes are particularly useful for advancing on the recognition of unknown diversity (in groups where alpha taxonomy is still developing or for highly diverse groups where many species remain to be described), it is also extremely useful for species identification in groups with better resolution in their taxonomy (Collins and Cruickshank 2012, Hebert and Gregory 2005, Hebert et al. 2003). Birds are one of the most well-known groups in terms of their taxonomy and systematics (Jarvis et al. 2014); nonetheless, new species are being described almost every year, particularly in the Neotropics (Avendano et al. 2015, Avendano et al. 2017) and some challenges still remain in the identification of species groups with very little phenotypic differentiation (Lara et al. 2012, Tavares et al. 2011, Cadena et al. 2016). Consequently, birds are an excellent group for implementing DNA barcoding for both species-identification and species-recognition purposes. Barcode studies in Neotropical birds are on the increase, particularly in Brazil and Argentina, where studies have focused on testing species limits and biogeographic patterns (Chaves et al. 2015, Kerr et al. 2009, Tavares et al. 2011, Vilaca et al. 2006). Despite that progress, a huge gap in information remains to be filled in northern South America, where very few studies have been completed (but see Mendoza et al. 2016). Here, we report the barcodes for 19 low elevation bird species from northern Colombia with the goal to provide tools for species identification and add to the existing gap in animal COI data in this part of South America. Materials and Methods Sample collection and processing We sampled birds at “Hacienda Universidad de Antioquia”, in the Municipality of Caucasia, Department of Antioquia, Colombia (8.003143 N, -75.400716 W; 70 m a.s.l., Fig. 1), from the 26th to 29th of October 2017. The landscape at the study site is composed of remnants of low land riparian forests, immersed in a matrix of pastures, secondary vegetation, shrubs and small streams (IDEAM 2012). Birds were caught at forest edges and in open areas between forest fragments using standard 3 x 12 m mist-nets and were blood- sampled from the brachial vein, using small gauge needles and non-heparinised capillary tubes (Pulgarin-R et al. 2018). All captured birds were processed, identified using field DNA barcodes of birds from northern Colombia 3 guides (Hilty et al. 1986, Donegan and McMullan 2014) and finally released in place. Blood samples were stored in 90-95% ethanol and kept at room temperature (Pulgarin-R et al. 2018). Legend i eS La Candelaria Corine Land Cover for Colombia GE 1.1.1. Continuous urban fabric MB 1.1.2. Discontinuous urban fabric IY) 2.3.1. Pastures 2.4.4, Heterogeneous agricultural areas GE 3.1.1. Dense forest 3.1.3. Fragmented forest GN 3.1.4. Gallery forest QO 3.1.5. Forest plantation ME 3.2.3. Secondary vegetation (9) 4.1.1. inland marshes MO 5.1.1. Rivers (50m) | | 75.250°W 75.035°W Figure 1. EESl Study area in the lowlands of northern Colombia. Laboratory Procedures We extracted total DNA from blood using the PureLink Genomic DNA Mini Kit (Invitrogen) according to the manufacturer’s specifications. For blood samples, 20 ul of Proteinase K, 20 ul of RNase and 200 ul of PureLink® Genomic Lysis/Binding buffer were added during the digestion phase. Later, each sample was transferred to a spin column and two washes were performed with Wash Buffer 1 and Wash Buffer 2 to perform a final elusion, dividing the total volume into two consecutive sets of 50 ul with Elution Buffer. For molecular typing, we targeted the Cytochrome c oxidase subunit 1 (COI) barcode region, using the primer combination from Ivanova et al. (2007) with the unique difference that all primers were M13-tailed to facilitate the sequencing process (Table 1). PCR amplifications were performed in 35 ul reactions that contained: 2 mM of MgClo, 1 * of buffer PCR 10 x with KCI, 0.2 mM of each dNTP, 0.14 ul of each primer cocktail, 1U of Tag DNA Polymerase (Fermentas) and 100 ng of DNA template. Thermal cycling conditions involved an initial denaturation at 95°C for 2 min followed by a single stage of 28 cycles that included denaturation at 95°C for 30 s, annealing at 52°C for 40 s, extension at 72°C for 1 min and a final 10 min extension at 72°C. PCR products were visualised on a 1.5% agarose gel, using a MiniBIS Pro-DNR Bio Imaging Systems. All amplification products 4 Pulgarin-R P et al were purified using Shrimp Alkaline Phosphatase and sent to Macrogen (Seoul, Korea) to be sequenced on an ABI PRISM 3100 Genetic Analyzer (Applied Biosystems, Carlsbad, CA, USA). Table 1. Primers used for the amplification of COI sequences obtained in this study. Name Sequence + M13 Ratio Source LepF1_t1- GTAAAACGACGGCCAGTATTCAACCAATCATAAAGATAT IGG 4 Ivanova et al. 2007 M13FWD VF1_t1- GTAAAACGACGGCCAGTTTCTCAACCAACCACAAAGACATTGG 1 Ivanova et al. 2007 M13FWD VF 1d_t1- GTAAAACGACGGCCAGTTTCTCAACCAACCACAARGAYATYGG 1 Ivanova et al. 2007 M13FWD VF 1i_t1- GTAAAACGACGGCCAGTTTCTCAACCAACCAIAAIGAIATIGG 3 Ivanova et al. 2007 M13FWD LepRI_t1- CAGGAAACAGCTATGACCTAAACTTCTGGATGTCCAAAAAATCA 1 Ivanova et al. 2007 M13REV VR1id_t1- CAGGAAACAGCTATGACCTAGACTTCTGGGTGGCCRAARAAYCA 1 Ivanova et al. 2007 M13REV VR1_t1- CAGGAAACAGCTATGACCTAGACTTCTGGGTGGCCAAAGAATCA | 1 Ivanova et al. 2007 M13REV VR1i_t1- CAGGAAACAGCTATGACCTAGACTTCTGGGTGICCIAAIAAICA 3 Ivanova et al. 2007 M13REV M13REV CAGGAAACAGCTATGACC NA Beckman Coulter, Inc 2020 M13FWD GTAAAACGACGGCCAGT NA Beckman Coulter, Inc 2020 Data analysis Sequences were edited, assembled and examined with reference to translated amino-acid sequences, using Geneious PRO 6.1.6. Nucleotide-sequences and complementary information were deposited in BOLD (www.barcodinglife.org) with the accession number dataset CANDE030-20 to CANDE055-20. For an initial sequence quality check and provisionary identification, all assembled sequences were searched in the National Centre for Biotechnology Information (NCBI) database through BLAST (http://BLAST.ncbi. nim.nih.gov/BLAST.cgi), using the Geneious Pro 6.1.6 match tool. We used the top- matching hit having the highest (> 98%) maximal percent identity score as criteria for successful conspecific/congeneric identification. After the initial BLAST-based identification on the NCBI database, we used the Animal Identification (COI) tool from the BOLD Identification System (IDS), using the Species Level Barcode Records database. For all Our Sequences, we recovered the species identification, closest matching BIN (Table 2) and a Neighbour-Joining topology, using Kimura-2-Parameter (K2P) substitution model as implemented in the BOLD portal (Suppl. material 1). Table 2. Individuals sampled and barcoded in this study. Individuals with * represent boreal migrants. Bolded taxa represent inconsistencies between our identification methods (see text). Code LCA35 LCA9 LCA12 LCA30 LCA3 LCA26 LCA27 LCA28 LCA24 LCA33 LCA20 LCA6 LCA18 LCA31 LCA36 LCA38 LCA4 LCA7 LCA21 LCA22 LCA13 LCA15 Field ID Automolus ochrolaemus Cantorchilus leucotis Catharus minimus Catharus minimus Catharus ustulatus Chaetura sp Chaetura sp Chaetura sp Coereba flaveola Dendrocincla fuliginosa Elaenia flavogaster Elaenia flavogaster Galbula ruficauda Manacus manacus Manacus manacus Manacus manacus Momotus subrufescens Momotus subrufescens Myiodinastes maculatus Myiozetetes cayanensis Parkesia noveboracensis Ramphocelus dimidiatus DNA barcodes of birds from northern Colombia BOLD ID A. ochrolaemus C. leucotis C. minimus C. minimus C. ustulatus* C. brachyura C. brachyura C. brachyura C. flaveola D. fuliginosa E. flavogaster E. flavogaster G. ruficauda M. aurantiacus M. aurantiacus M. aurantiacus M. momota M. momota M. luteiventris M. cayanensis P. noveboracensis R. dimidiatus BOLD NCBIID hit’ (%) 100 100 100 100 100 100 100 100 100 99.85 99.85 100 100 100 100 100 100 100 100 99.85 100 100 A. ochrolaemus C. leucotis C. minimus C. minimus* C. ustulatus* C. brachyura C. brachyura C. brachyura C. flaveola D. fuliginosa E. flavogaster E. flavogaster G. ruficauda M. manacus M. manacus M. manacus M. momota M. momota M. luteiventris * M. cayanensis P noveboracensis R. carbo NCBI hit"(%) 99.10 95.55 100 100 100 100 100 100 100 99.15 98.93 99.39 97.55 100 100 100 96.92 97.41 100 98.77 99.85 99.39 Consensus sp BOLD ID A. ochrolaemus C. leucotis C. minimus* C. minimus* C. ustulatus* C. brachyura C. brachyura C. brachyura C. flaveola D. fuliginosa E. flavogaster E. flavogaster G. ruficauda M. aurantiacus M. aurantiacus M. aurantiacus M. momota M. momota M. luteiventris* M. cayanensis P. noveboracensis* R. dimidiatus Seq length (bp) 671 690 657 660 702 644 642 652 651 673 681 696 675 667 663 667 681 657 651 660 658 681 BOLD:ADM4531 BOLD:ABX4224 BOLD:AAA9441 BOLD:AAA9441 BOLD:AAA9440 BOLD:AAK0488 BOLD:AAK0488 BOLD:AAK0488 BOLD:AAA4006 BOLD:ABZ6107 BOLD:AAB3859 BOLD:AAB3859 BOLD:ABX4491 BOLD:AAB9291 BOLD:AAB9291 BOLD:AAB9291 BOLD:ABX4186 BOLD:ABX4186 BOLD:AAF5348 BOLD:AAE6211 BOLD:AAB0401 BOLD:AAD5047 6 Pulgarin-R P et al Code __ Field ID BOLD ID BOLD NCBIID NCBI Consensus sp Seq BIN hit! hit"(%) BOLDID length (%) (bp) LCA40 Ramphocelus R. dimidiatus 100 R. carbo 99.23 R. dimidiatus 654 BOLD:AAD5047 dimidiatus LCA1 Sporophila S. funerea 100 S. angloennsis 98.92 S. funerea 687 BOLD:AAE5360 funerea LCA19 Tolmomyias E 99.85 T. 97.89 TT. sulphurescens 666 BOLD:ACI3658 sulphurescens sulphurescens sulphurescens LCA8 Xiphorhynchus’ X.susurrans 99.54 X. guttatus 98.15 xX. susurrans 670 BOLD:ACF 1637 susurrans Results We obtained 26 COI sequences from 18 species, 10 families and three orders and, when analysed by BOLD, the species were grouped into 18 existing BINS (access numbers in Table 2). Most bird species were residents, but four species (Catharus minimus, Catharus ustulatus, Myiodynastes luteiventris and Parkesia noveboracensis), were boreal migrants (Ayerbe 2018, Avendano et al. 2015, Avendano et al. 2017). For all species, sequence lengths varied from 642 to 702 bp (Table 2). Since most bird species are under-sampled for DNA barcodes (Mendoza et al. 2016) in this part of the tropics, our report represents an important contribution to expand the geographic sampling (for COI] sequences) of several species in South America and it also includes the first sequences for Colombia for the following species: Cantorchilus leucotis, Chaetura brachyura, Galbula_ruficauda, Myiodynastes_ luteiventris, Myiozetetes cayanensis, Tolmomyias sulphurescens and Xiphorhynchus susurrans. Most COI barcodes matched our initial phenotypic identification; however, for six (6) species, (10 individuals), we found differences between our field identification, the query hits from BOLD’s IDS and the NCBI BLAST search (Table 2). One bird species originally identified in the field as Myiodynastes maculatus, field ID LCA21) was positively identified as Myiodynastes luteiventris, (Fig. 2) by BOLD and NCBI analyses (but see Discussion). Another passerine species correctly identified in the field (field ID's LCA31, LCA36, LCA38) and by the NCBI BLAST as Manacus manacus (Fig. 3), was recovered as the Central American restricted species, Manacus aurantiacus by BOLD’s IDS. A third species was identified in the field (and NCBI BLAST search) as Momotus subrufescens (Fig. 4), but BOLD IDS recovered its former nominal assignation, Momotus momota, the name of a widely distributed form of motmot before it was split into five species-level taxa (Stiles 2009). Additionally, other three species (Ramphocelus dimidiatus, Sporophila funerea and Xiphorhynchus susurrans) were positively identified in the field and by BOLD, but exhibited erroneous identifications by the NCBI BLAST apparently because of the absence of COl sequences for either species in the latter portal. Finally, in six instances, DNA sequences helped to confirm the identification of Automolus ochrolaemus, Chaetura brachyura and Tolmomyias sulphurescens, which are all species difficult to identify in the field, even in hand, particularly the swifts. DNA barcodes of birds from northern Colombia Figure 2. EES] Myiodynastes luteiventris (but see Discussion), a boreal migrant, was initially identified in the field as Myiodynastes maculatus and was subsequently re-identified with the help of its COI barcode (BOLD ID). Figure 3. EES Manacus manacus was identified as a different manakin species according to BOLD. 8 Pulgarin-R P et al Momotus momota (43): Bolivia, La Paz Momotus momota (44): Bolivia, La Paz Momotus momota (45): Bolivia, La Paz Momotus momota [46]: Indet Momotus momota [47]: Bolivia, La Paz Momotus momota [48]: Bolivia, Santa Cruz Momotus momota [49}: Bolivia, La Paz Momotus momota (50): Peru, Ucayali Momotus momota (51): Peru, Ucayali Momotus momota (52): Peru, Ucayali Momotus momota [53]: Peru, Madre de Dios Momotus momota (54): Peru, Ucayali Momotus momota [55]: Brazil, Roraima Momotus momota (56): French Guiana Momotus momota (57): French Guiana Momotus momota (58): French Guiana Momotus momota (59): French Guiana Momotus momota (60): French Guiana Momotus momota [61]: French Guiana Momotus momota [62): Colombia, Vichada Momotus momota [63]: Colombia, Vichada Momotus momota (64): Ecuador, Sucumbios Momotus momota (65): Panama, Colon Momotus momota (66): Panama, Darien Momotus momota [67]: Panama, Colon Momotus subrufescens Momotus momota (68): Panama, Darien LCA7/LCA4: Colombia, Antioquia Momotus momota (69): Mexico, Chiapas Momotus momota [70]: Mexico, Tabasco Momotus momota {71}: Guatemala Momotus momota [72]: Guatemala, |zabal Momotus lessonii Momotus momota [73]: Panama, Veraguas Momotus momota (74): Panama, Chiriqui Momotus momota [75]: Panama, Chiriqui Momotus aequatorialis (76): Peru, Ayacucho | Momotus mexicanus [77]: Mexico Momotus mexicanus [78]: Mexico Momotus mexicanus [79]: Mexico Momotus mexicanus [80]: Mexico Momotus mexicanus (81): Mexico Momotus mexicanus [82]: Mexico Momotus mexicanus (83): Mexico Momotus mexicanus (84): Mexico Momotus mexicanus (85): Mexico Momotus mexicanus (86): Mexico Momotus mexicanus (87): Mexico Momotus mexicanus Momotus mexicanus (88): Mexico Momotus mexicanus [89]: Mexico Momotus mexicanus [90]: Mexico Momotus mexicanus {91}: Mexico Momotus mexicanus [92]: Mexico Momotus mexicanus [93]: Mexico Momotus mexicanus (94): Mexico Momotus mexicanus (95): Mexico Momotus mexicanus [96]: Mexico Momotus mexicanus [97]: Mexico Momotus mexicanus (98): Mexico Momotus mexicanus [99]: Mexico Momotus momota To Baryphthengus Momotus aequatorialis 0.02 substitutions/site Figure 4. EES] Kimura-2-parameter tree (obtained from the BOLD portal) of the "Womotus momota complex" showing the updated taxonomic arrangement of this clade. Numbers in brackets in the terminals correspond to BOLD numeric descriptors for each sample. Discussion Our assessment of species identification, using the COI barcodes, shows a strong correspondence (90%) with field identification, based on research expertise and photo ID (Table 2, Suppl. material 1). However, DNA barcodes were able to help with the identification of challenging species that can be problematic even for trained neotropical ornithologists. This was the case of field ID LCA21, identified initially as Myiodynastes maculatus (Hilty et al. 1986, Donegan and McMullan 2014), but for which both NCBI DNA barcodes of birds from northern Colombia 9 BLAST and the BOLD identification tool recovered it as M. /uteiventris (Fig. 2). According to traditional and recent literature, some relevant diagnostic characters to identify MV. maculatus include: insinuation or presence of rufous colouration on the margins of the primary feathers and coverts, a broadly pink lower-mandibular base (dark only the distal half), rufous tail with dark central stripe, narrow dusky malar stripe that does not meet under the bill (usually paler than that of M. /uteiventris) and rufous or buffy supercilium (Shah 2020, Lowther and Stotz 2020, Donegan and McMullan 2014, Hilty et al. 1986, Ayerbe 2018). Although phenotypic characteristics of our specimen match those of MV. maculatus (see Fig. 2), its COl sequence was grouped within the unique BIN containing /. luteiventris sequences (see Table 2). Moreover, its nearest neighbour is a M. maculatus BIN, with a strikingly large COI distance of 7.22%, reducing the possibility of a misidentification problem by BOLD's database. Considering this contradictory evidence (i.e. the phenotypic resemblance of our specimen with M. maculatus and the strong mitochondrial association with M. /uteiventris), we cannot rule out the possibility of an introgression event, a phenomenon that has been documented amongst closely-related species with sympatric distributions in the family Tyrannidae (Ottenburghs et al. 2017, Rheindt et al. 2009, Rheindt et al. 2013, Winger 2017). Although evaluating a possible introgression scenario is outside of the scope of this study, it is important to highlight that barcoding studies can give us clues to understand these events. Similarly, barcodes might help to identify the breeding areas or population origin for species exhibiting migratory divide or genetic structure, as happened with passing through northern South America species, Catharus minimus and Catharus ustulatus (Topp et al. 2013, Pulgarin-R et al. 2018). Additionally, barcodes can be of great help in resident species with little phenotypic variation, such as the swifts in the genus Chaetura, which are hard to capture in mist-nets and hard to identify in the field. We also found some discrepancies between IDs recovered by the NCBI BLAST tool, those recovered by BOLD and our initial identifications made in the field. For example, three specimens identified in the field (Fig. 3) and by the NCBI BLAST tool as Manacus manacus were recovered by BOLD as M. aurantiacus. The BIN containing our sequences (Table 2) groups several phenotypes that, in the past, have been treated as the same species (e.g. Snow 1975) and also as a Ssuperspecies with up to four species (IV. aurantiacus, M. candel, M. manacus and WM. vitellinus) (Snow et al. 2004). Taxonomy within this group is not fully resolved so far that M. aurantiacus has been considered a subspecies of MW. vitellinus (Snow 1975), an independent allopatrically-distributed species of the genus (Brumfield and Braun 2001, Brumfield et al. 2001,Brumfield et al. 2008) and even as a paraphyletic clade, based on mtDNA (Brumfield and Braun 2001). Moreover, it has been found that species of Manacus can hybridise in areas of sympatric distribution with other species of the genus and even the family (Brumfield and Braun 2001, Brumfield et al. 2001, H6glund and Shorey 2004). All the above-mentioned scenarios indicate that, although the phenotype of all our sequences corresponds to what is known as M. manacus (Fig. 3), in the absence of a clear phylogenetic arrangement and poor knowledge on the species limits within the genus Manacus, the DNA barcode, by itself, is not able to reconcile the morphological and molecular information and is only the reflection of a poorly understood taxonomy. 10 Pulgarin-R P et al Another result that showed some inconsistencies was the identification of LCA4 and LCA7 sequences, which were recovered by BOLD as Momotus momota (Fig. 4). This used to be a widely-distributed species in Central and South America, until it was split into five species-level taxa (M. aequatorialis, M. bahamensis, M. lessonii, M. momota and WM. subrufescens), using a combined analysis of plumage, biometrics and voice (Stiles 2009). Currently, MM. momota is considered a cis-Andean distributed species from eastern Colombia to southern Venezuela, Guianas, north-western Argentina and most of Brazil (Stiles 2009). Particularly, the specific epithet, associated with the populations and phenotype obtained in this study, corresponds to M. subrufescens; however, despite the presence of eight different BINs that span much of the distribution of all the mentioned species within the genus, the taxonomy within the BOLD portal has yet to be updated and, consequently, our sequence is part of a BIN based on a haplotype with geographical proximity that bears the outdated MW. momota taxon name. A final group of inconsistencies between identification methods corresponds to three (3) species for which no COI sequence data are available at the NCBI portal and, consequently, their closest matching sequences are inconsistent with their correct field- and BOLD-based identifications. In the case of the genera Ramphocelus and Xiphorhynchus, the BLAST search tool identified our samples as the cis-Andean distributed congeneric species (R. carbo and X. guttatus) and not as the correct trans- Andean species, R. dimidiatus and xX. sussurrans (Ayerbe 2018). For the genus Sporophila, although S. angolensis and S. funerea can show sympatric distributions (Ayerbe 2018), the morphology exhibited by their males is strikingly different and leaves no room for discussion on their morphological identification. Even though we found some discrepancies between our identification methods compared to BOLD’s IDS, a close inspection to the K2P trees from BOLD (Fig. 4) showed that individuals across all sampled species are closely related to other individuals from nearby populations/areas. This is an important fact because, even in the presence of outdated or incorrect assignment of names to a barcode sequence (and, consequently, to its corresponding BIN), the K2P topologies are able to group individuals that, based even on geography itself, can putatively represent the current taxonomical treatment of the species (as is the case with the Vomotus subrufescens in Fig. 4). Conclusions Despite our reduced sampling, this study represents the first attempt to document COl barcodes for birds (from blood samples) in this part of Colombia, which fills a considerable gap in sampling in north-western South America. Particularly, a call for broader sampling for barcodes might provide hints on cryptic species across barriers (Barreira et al. 2016) or might facilitate the identification of highly-traded species in Colombia, such as parrots (Mendoza et al. 2016, Restrepo-R and Pulgarin-R 2017). DNA barcodes of birds from northern Colombia 11 Acknowledgements This research was partially funded by Universidad de Antioquia through the grant "Sostenibilidad 2016-2017 Codi". Fieldwork was developed under the permit "Resolucion 0524 del 27 de mayo de 2014" issued by the Colombian National Authority of Environmental Licences (ANLA, by its initials in Spanish) to Universidad de Antioquia. We are grateful to Diego Calderén-Franco and two anonymous reviewers for their comments that helped to improve the present manuscript. References ° Avendano E, Bohorquez |, Rosselli L, Arzuza-Buelvas D, Estela FA, Cuervo AM, Renjifo M (2017) Lista de chequeo de las aves de Colombia: una sintesis del estado del conocimiento desde Hilty & Brown (1986). Ornitologia Colombiana 16. ° Avendano JE, Cuervo A, Lpez-O JP, Gutirrez-Pinto N, Corts-Diago A, Cadena CD (2015) Anew species of tapaculo (Rhinocryptidae: Scytalopus) from the Serrania de Perija of Colombia and Venezuela. The Auk 132 (2): 450-466. https://doi.org/10.1642/ AUK-14-166.1 ° Ayerbe F (2018) Guia Ilustrada de la Avifauna colombiana. [An ilustrated field guide to the Birds of Colombia]. 1. WCS, Bogota. ° Barreira A, Lijtmaer D, Tubaro P (2016) The multiple applications of DNA barcodes in avian evolutionary studies. Genome 59 (11): 899-911. https://doi.org/10.1139/ gen-2016-0086 ° Beckman Coulter, Inc (2020) Beckman Coulter, Inc. Definitions https://doi.org/10.32388/ kfgmki ° Brumfield R, Braun M (2001) Phylogenetic relationships in bearded manakins (Pipridae: Manacus) indicate that male plumage color is a misleading taxonomic marker. The Condor 103 (2): 248-258. https://doi.org/10.1093/condor/103.2.248 ° Brumfield R, Jernigan R, McDonald D, Braun M (2001) Evolutionary implications of divergent clines in avian (Manacus: Aves) hybrid zone. Evolution 55 (10). httos://doi.org/ 10.1554/0014-3820(2001)055[2070:eiodci]2.0.co;2 ° Brumfield RT, Liu L, Lum DE, Edwards SV (2008) Comparison of species tree methods for reconstructing the phylogeny of bearded manakins (Aves: Pipridae, Manacus) from multilocus sequence data. Systematic Biology 57 (5): 719-31. https://doi.org/ 10.1080/10635150802422290 ° Cadena CD, Caro ML, Caycedo-Rosales PC, Cuervo AM, Bowie RCK, Slabbekoorn H (2016) Henicorhina anachoreta (Troglodytidae), another endemic bird species for the Sierra Nevada de Santa Marta, Colombia. Ornitologia Colombiana 15 (May): 82-89. ° Chaves BN, Chaves A, Nascimento AA, Chevitarese J, Vasconcelos M, Santos F (2015) Barcoding neotropical birds: Assessing the impact of nonmonophyly in a highly diverse group. Molecular Ecology Resources 15 (4): 921-931. https://doi.org/ 10.1111/1755-0998.12344. ° Collins RA, Cruickshank RH (2012) The seven deadly sins of DNA barcoding. Molecular Ecology Resources 13 (6): 969-975. https://doi.org/10.1111/1755-0998.12046 12 Pulgarin-R P et al Crawford A, Cruz C, Griffith E, Ross H, lbafez R, Lips K, Driskell A, Bermingham E, Crump P (2012) DNA barcoding applied toex situtropical amphibian conservation programme reveals cryptic diversity in captive populations. Molecular Ecology Resources 8 (6). https://doi.org/10.1111/1755-0998.12054 de Kerdrel GA, Andersen JC, Kennedy SR, Gillespie R, Krehenwinkel H (2020) Rapid and cost-effective generation of single specimen multilocus barcoding data from whole arthropod communities by multiple levels of multiplexing. Scientific Reports 10 (1): 78. https://doi.org/10.1038/s41598-019-54927-z Donegan T, McMullan M (2014) Guia de las aves de Colombia. 2. Proaves Gaytan A, Bergsten J, Canelo T, Pérez-Izquierdo C, Santoro M, Bonal R (2020) DNA Barcoding and geographical scale effect: The problems of undersampling genetic diversity hotspots. Ecology and Evolution 10 (19): 10754-10772. https://doi.org/10.1002/ ece3.6733 Gonzalez MA, Baraloto C, Engel J, Mori S, Pétronelli P, Riera B, Roger A, Théebaud C, Chave J (2009) Identification of Amazonian trees with DNA barcodes. PLOS One 4 (10). https://doi.org/10.1371/journal.pone.0007483 Hebert PN, Cywinska A, Ball S, deWaard J (2003) Biological identifications through DNA barcodes. Proceedings of the Royal Society of London. Series B: Biological Sciences 270 (1512): 313-321. https://doi.org/10.1098/rspb.2002.2218 Hebert PN, Gregory TR (2005) The promise of DNA barcoding for taxonomy. Systematic Biology 54 (5): 852-859. https://doi.org/10.1080/10635150500354886 Hilty SL, Brown WL, Brown B (1986) A guide to the birds of Colombia. 1. Princeton University Press Hdglund J, Shorey L (2004) Genetic divergence in the superspecies Manacus. Biological Journal of the Linnean Society 81 (3): 439-447. https://doi.org/10.1111/j. 1095-8312.2003.00297.x IDEAM (2012) Leyenda Nacional de Coberturas de la Tierra. Metodologia CORINE Land Cover adaptada para Colombia. Escala 1:100.000. Ivanova N, Zemlak T, Hanner R, Hebert PN (2007) Universal primer cocktails for fish DNA barcoding. Molecular Ecology Notes 7 (4): 544-548. https://doi.org/10.1111/j. 1471-8286.2007.01748.x Jarvis ED, Mirarab S, Aberer AJ, Li B, Houde P, Li C, et al. (2014) Whole-genome analyses resolve early branches in the tree of life of modern birds. Science 346 (6215): 1320-31. https://doi.org/10.1126/science.1253451 Kerr KR, Lijtmaer D, Barreira A, Hebert PN, Tubaro P (2009) Probing evolutionary patterns in neotropical birds through DNA barcodes. PLOS One 4 (2). https://doi.org/ 10.1371/journal.pone.0004379 Ko H, Wang Y, Chiu T, Lee M, Leu M, Chang K, Chen W, Shao K (2013) Evaluating the accuracy of morphological identification of larval fishes by applying DNA barcoding. PLoS ONE 8 (1). https://doi.org/10.1371/journal.pone.0053451 Lara CE, Cuervo AM, Valderrama SV, Calderon D, Cadena D (2012) A new species of wren (Troglodytidae: Thryophilus) from the dry Cauca river canyon, northwestern Colombia. The Auk 129 (3): 537-550. https://doi.org/10.1525/auk.2012.12028 Lohman D, Ingram K, Prawiradilaga D, Winker K, Sheldon F, Moyle R, Ng PL, Ong P, Wang LK, Braile T, Astuti D, Meier R, et al. (2010) Cryptic genetic diversity in “widespread” Southeast Asian bird species suggests that Philippine avian endemism is DNA barcodes of birds from northern Colombia 13 gravely underestimated. Biological Conservation 143 (8): 1885-1890. htips://doi.org/ 10.1016/j.biocon.2010.04.042 Lowther PE, Stotz DF (2020) Sulphur-bellied Flycatcher (Myiodynastes luteiventris), version 1.0. In: Poole AF, Gill FB (Eds) Birds of the World. Cornell Lab of Ornithology, Itaca, New York. https://doi.org/10.2173/bow.subfly.01 Mendoza AM, Torres MF, Paz A, Trujillo-Arias N, Lopez-Alvarez D, Sierra S, Forero F, Gonzalez M (2016) Cryptic diversity revealed by DNA barcoding in Colombian illegally traded bird species. Molecular Ecology Resources 16 (4): 862-873. htips://doi.org/ 10.1111/1755-0998.12515 Ottenburghs J, Kraus RS, van Hooft P, van Wieren S, Ydenberg R, Prins HT (2017) Avian introgression in the genomic era. Avian Research 8 (1). https://doi.org/10.1186/ $40657-017-0088-z Pulgarin-R P, Gomez C, Bayly N, Bensch S, FitzGerald A, Starkloff N, Kirchman J, Gonzalez-Prieto A, Hobson K, Ungvari-Martin J, Skeen H, Castaho MI, Cadena CD (2018) Migratory birds as vehicles for parasite dispersal? Infection by avian haemosporidians over the year and throughout the range of a long-distance migrant. Journal of Biogeography 46 (1): 83-96. https://doi.org/10.1111/jbi.13453 Restrepo-R DC, Pulgarin-R PC (2017) Dinamicas de los loros en cautiverio en Colombia: trafico, mortalidad y liberacion. Ornitologia Colombiana 16 (December): eA1-eA23. Rheindt FE, Christidis L, Norman JA (2009) Genetic introgression, incomplete lineage sorting and faulty taxonomy create multiple cases of polyphyly in a montane clade of tyrant flycatchers (Elaenia, Tyrannidae). Zoologica Scripta 38: 143-145. https://doi.org/ 10.1111/j.1463-6409.2008.00369.x Rheindt FE, Fujita MK, Wilton PR, Edwards SV (2013) Introgression and phenotypic assimilation in Zimmerius flycatchers (Tyrannidae): population genetic and phylogenetic inferences from genome-wide SNPs. Systematic Biology 63 (2): 134-52. https://doi.org/ 10.1093/sysbio/syt070 Shah SS (2020) Streaked Flycatcher (Myiodynastes maculatus), version 1.0. In: Schulenberg T (Ed.) Birds of the World. Cornell Lab of Ornithology, Ithaca, New York.. https://doi.org/10.2173/bow.strfly1.01 Shaw D, Escalante P, Rappole J, Ramos M, Oehlenschlager R, Warner D, Winker K (2013) Decadal changes and delayed avian species losses due to deforestation in the northern Neotropics. PeerJ 1 https://doi.org/10.7717/peerj.179 Snow D (1975) The classification of the manakins. Bulletin of the British Ornithological Club 95: 20-27. Snow D, Del Hoyo J, Elliot A, Christie D (2004) Familly Pipridae (Manakins). Bulletin of British Ornithological Club 95. Stefan L, GOmez-Diaz E, Mironov S, Gonzalez-Solis J, McCoy K (2018) “More than meets the eye”: Cryptic diversity and contrasting patterns of host-specificity in feather mites inhabiting seabirds. Frontiers in Ecology and Evolution 6: 1-97 hitps://doi.org/ 10.3389/fevo.2018.00097 Stiles FG (2009) A review of the genus Momotus (Coraciiformes: Momotidae) in northern South America and adjacent areas. Ornitologia Colombiana 8: 29-75. URL: http://zoobank.org/c135a921-ca07-4c6a-be68-07ff3a710d17 14 Pulgarin-R P et al ° Tavares ES, Gongalves P, Miyaki CY, Baker A (2011) DNA Barcode detects high genetic structure within neotropical bird species. PLOS One 6 (12). httos://doi.org/10.1371/ journal.pone.0028543 ° Topp C, Pruett C, McCracken K, Winker K (2013) How migratory thrushes conquered northern North America: A comparative phylogeography approach. PeerJ 1 https://doi.org/10.7717/peerj.206 ° Vilaca TS, Lacerda DR, Sari HR, Santos FR (2006) DNA-based identification applied to Thamnophilidae (Passeriformes) species: The first barcodes of neotropical birds. Revista Brasileira de Ornitologia 14 (1): 7-13. ° Waugh J (2007) DNA barcoding in animal species: Progress, potential and pitfalls. BioEssays 29 (2): 188-197. https://doi.org/10.1002/bies.20529 ° Winger BM (2017) Consequences of divergence and introgression for speciation in Andean cloud forest birds. Evolution; International Journal of Organic Evolution 71 (7): 1815-1831. https://doi.org/10.1111/evo.13251 Supplementary material Suppl. material 1: Suplementary Information EE) Authors: Pulgarin et al. Data type: DNA sequences, data Tables. Download file (31.54 kb)