CO1-Based DNA barcoding for assessing diversity of Pteropus giganteus from the state of Azad Jammu Kashmir, Pakistan

Código de barras de DNA à base de CO1 para avaliar a diversidade de Pteropus giganteus do estado de Azad Jammu e Caxemira, Paquistão

Abstract

The flying fox (Pteropus giganteus) also familiar with the name of the greater Indian fruit Bat belongs to the order Chiroptera and family Pteropodidae. Current research emphasis on the DNA barcoding of P. giganteus in Azad Jammu Kashmir. Bat sequences were amplified and PCR products were sequenced and examined by bioinformatics software. Congeneric and conspecific, nucleotide composition and K2P nucleotide deviation, haplotype diversity and the number of haplotypes were estimated. The analysis showed that all of the five studied samples of P. giganteus had low G contents (G 19.8%) than C (27.8%), A (25.1%) and T (27.3%) contents. The calculated haplotype diversity was 0.60% and the mean intraspecific K2P distance was 0.001% having a high number of transitional substitutions. The study suggested that P. giganteus (R=0.00) do not deviate from the neutral evolution. It was determined from the conclusion that this mtDNA gene is a better marker for identification of Bat species than nuclear genes due to its distinctive characteristics and may serve as a landmark for the identification of interconnected species at the molecular level and in the determination of population genetics.

Keywords:
Pteropodidae; evolution; mtDNA; population genetics; haplotype diversity; nucleotide diversity; flying fox

Resumo

A raposa-voadora (Pteropus giganteus), também conhecida como morcego indiano, pertence à ordem dos Chiroptera e à família Pteropodidae. A presente pesquisa dá ênfase ao código de barras de DNA de P. giganteus em Azad Jammu e Caxemira. Sequências genéticas dos morcegos foram amplificadas, e os produtos de PCR foram sequenciados e examinados por software de bioinformática. De espécies congenérica e coespecífica, foram estimados composição nucleotídica e desvio de nucleotídeos K2P, diversidade de haplótipos e número de haplótipos. A análise mostrou que todas as cinco amostras estudadas de P. giganteus apresentaram baixos teores de G (19,8%) em comparação com C (27,8%), A (25,1%) e T (27,3%). A diversidade de haplótipos calculada foi de 0,60%, e a distância média intraespecífica de K2P foi de 0,001%, com um elevado número de substituições transicionais. O estudo sugeriu que P. giganteus (R = 0,00) não se desviou da evolução neutra. É possível concluir que o gene mtDNA é um marcador favorável para identificação de espécies de morcegos do que genes nucleares por causa de suas características distintivas e pode servir como um marco para a identificação de espécies interconectadas em nível molecular e para a determinação genética de populações.

Palavras-chave:
Pteropodidae; evolução; mtDNA; genética de populações; diversidade de haplótipos; diversidade de nucleotídeos; raposa-voadora

1. Introduction

Bats (Order: Chiroptera) offer important ecosystem services such as pollination, seed dispersal and insect pest suppression (Kunz et al., 2011). Flying foxes are mainly significant players in island ecosystems where they often work as seed dispersers and keystone pollinator both within and between islands (McConkey and Drake, 2007; McConkey and Drake, 2015). Other than that, flying foxes greatly maintains their numbers at high densities which is very essential for the existence of plant communities (McConkey and Drake, 2006). Chiroptera comprises many closely linked species which are sharing very similar ecologies and morphologies (Clare et al., 2007). These sort of similarities in morphology results in overlooked diversity (Mayer et al., 2007). Indian Flying Fox (P. giganteus) is an extensively dispersed species of Bats that present in the different tropical areas of central Asia, mostly between China and Pakistan (Bates and Harrison, 1997). In Pakistan, these Indian flying foxes are widely distributed in the areas of Sindh, Punjab and Northwest Frontier Province (Sheikh and Molur, 2004).

There is an estimation that about 10 to 100 million species may exist on our earth from which at least 1.5 million species of organisms are described (Wilson, 2003). Due to this enormous gap, there is a need to search for new technologies to describe the immense biological diversity. For this, the genetic screening methods which involve a single or a small number of genes – usually called DNA barcoding – were used in diverse taxonomic groups. There are two major aims for this: (i) to use a standardized technique for the identification of species, and (ii) to detect new species efficiently to bring us closer to the true number of species (Blaxter et al., 2004; Moritz and Cicero, 2004; Hebert and Gregory, 2005). Therefore, DNA barcodes were primarily used for genetic identification of taxonomically poorly studied taxa and geographic regions like the tropics (Blaxter et al., 2004; Monaghan et al., 2005; Smith et al., 2006).

DNA barcoding is one of the best-developed technique (Hebert et al., 2004) in which a portion of the mitochondrial cytochrome c oxidase subunit 1 (COI) gene is used. DNA barcoding techniques contributed to identifying the cryptic species which increased the knowledge of immense diversity in the different species of Chiroptera (Clare et al., 2007; Mayer et al., 2007). The COI gene size has been preserved completely and used in evolutionary studies. There is a huge number of unidentified species or individuals below the species level in some higher vertebrate classes that could be documented by this best-accomplished tool of identification. DNA barcoding technique is also applicable for biodiversity analysis and certification of species (Ospina-Alvarez and Piferrer, 2008).

Pakistan is recognized as a country of very rich biodiversity but most of the species, including relatively familiar mammalian groups, are poorly identified and there is a lack of prior research and published data (Sheikh and Molur, 2004). This limitation serves as an opportunity for further research to explore unidentified species, therefore, the present study was conducted and the objectives were: (i) to assess the taxonomic identification of Pteropus giganteus by molecular sequencing data, (ii) to explore the diversity of chiropteran (Genus Pteropus) in State of Azad Jammu Kashmir by using mitochondrial COI gene, and (iii) to establish a national reference library in GenBank for Bat (Indian flying fox) species of State of Azad Jammu Kashmir, Pakistan.

2. Material and Methods

Blood samples of bats were collected from the districts Muzaffarabad (34° 21' N and 73° 28' E) and Bagh (33° 58' N and 73° 46' E) Azad Jammu Kashmir, Pakistan. Bats were randomly captured for the blood collection and anesthetized while blood collection and were released after sampling (McMichael et al., 2015). Identification of the bats was based on morphological principles described in key taxonomic references including Bates and Harrison (1997), Borissenko and Kruskop (2003), and Francis (2008), as well as describing individual species by primary literature review. Blood was collected in the 3 ml EDTA tubes with the help of 3 ml syringes and samples were labeled and were brought back to the Laboratory of Genetics, University of Azad Jammu and Kashmir. Blood samples were stored in the fridge until DNA extraction has been done. A total of 5 specimens were barcoded, from which 3 samples (M1, M2 and M3) were collected from district Muzaffarabad and 2 samples (B1 and B2) were collected from district Bagh, State of Azad Jammu Kashmir.

DNA extraction was done by using a phenol-chloroform extraction procedure followed by Sambrook et al., (1989). The extracted DNA was allowed to run on 1% agarose gel and visualized under UV transilluminator and results were recorded by using a gel documentation system. The quantity was measured by a spectrophotometer at 260/280 nm wavelength. Online software Primer 3 was used for the designing of primers sequences of Cytochrome C Oxidase 1 (CO1) gene of Pteropus giganteus. Primer sequences used were: PG-COI-F (5´-AAC GGC CCT CAG CCT ACT AA-3´), PG-COI-R (5´-GAA AAA GGT GGT ATT TAG GTT-3´) (1000 bps product size).

COI gene was amplified by using PCR (Polymerase Chain Reaction). PCR was carried out in 25μl (standard amount) reaction volumes in each tube (14 μl PCR water, 3 μl template DNA, 2.5 μl Taq buffer, 0.5 μl dNTPs (2.5 μM), primers 1 μl of each (10 μM), 2.5 μl magnesium chloride, and 0.5 U Taq polymerase. For thorough mixing, the reaction mixture kept for 30 seconds on 8,000 rpm for centrifugation. Thermal cycling comprised 95ºC for 3 min, followed by 39 cycles of 95ºC for 30 secs, annealing at 55-63 °C for 30 s, and an extension temperature of 72 °C for 1 min. This was followed by a final extension of 72 °C for 10 min. Before cycle sequencing, PCR products were purified with Exo Sap-IT (Affymetrix purification kit). Nucleotide (nt) sequencing was carried out in an ABI Prism 3100 Genetic Analyzer (PE Applied Biosystems, Foster City, CA, USA) using gene-specific forward and reverse primers for COI genes. One microliter cleaned PCR product was used for each 10 μL reaction. Indian flying fox gene sequences were edited by using the BioEdit program (Hall, 1999) for nucleotides and amino acid variations. Alignment of these nucleotide sequences was done with other already known international sequences by using Clustal W in the MegAlign program of LASERGENE package (DNASTAR Inc, Madison, WI, USA) and the international sequences were retrieved from the NCBI GenBank resource. After alignment and sequencing, these sequences were deposited in GenBank for accession numbers (Table 1).

Table 1
Details of COI gene sequences, voucher no, sequence length, Genbank accession numbers and collection locality of P. giganteus.

S.no	Voucher no.	Sequence length	NCBI Acc. No	Collection Locality
1	P. giganteus COI B1	627	MK415806	District Bagh, AJK
2	P. giganteus COI B2	627	MK415807	District Bagh, AJK
3	P. giganteus COI M1	627	MK415808	District Muzaffarabad, AJK
4	P. giganteus COI M2	627	MK415809	District Muzaffarabad, AJK
5	P. giganteus COI M3	627	MK415810	District Muzaffarabad, AJK

The identification of Bat samples was done through Barcode of Life Data Systems (BOLD) and Basic Local Alignment Search Tool (BLAST) of GenBank/NCBI. Sequence divergence between the Bat species of AJK (Pakistan) and other international reference species was determined by multiple sequence alignments. A phylogenetic tree was generated using the neighbor-joining method and Kimura 2-parameter supported by 1,000 bootstrap replicates (Akhtar and Ali, 2016; Akhtar et al., 2017; Ayesha et al., 2019) in MEGA 6.06 software. Replicate trees percentage in which all the related taxa are grouped in the bootstrap test (thousand replicates) were shown in the next branches of the MEGA 6.06 software. Considering the results of Gager et al. (2016) three Cynomops species (GenBank accession numbers JF447634, JN312044 and EF080319) and Molossus rufus (JF444936) were used as outgroups to root the tree inferred from CO1 sequences. Nucleotide composition and rate of transition/transversion ratio (R) were also calculated by MEGA 6.06 software. DNASP 5.0 program (Rozas et al., 2003) was used for the estimation of the number of polymorphic sites, no. of haplotypes, haplotype diversity (Hd), and nucleotide diversity (Pi).

3. Results

From (~1000 bps) PCR amplicons, the sequence used for comparison was the 627 bps consensus (Figure 1). Of the 627 sites, 626 (99.8%) were constant, 1 (0.15%) site was variable and also parsimony informative. The nucleotides’ translation into a total of 209 amino acids sequences produced the 188 (89.9%) constant sites, while 1 (0.47%) site was variable and also parsimony informative. The empirical base composition of the COI gene among the examined samples was T (27.3%), C (27.8%), A (25.1%) and G (19.8%). The frequencies of A+T contents (52.4%) were slightly higher than those of G+C contents (47.6%), resulting in anti-G bias sequencing which is a mitochondrial gene characteristic.

Figure 1
2% agarose gel profile of DNA amplification of the COI gene on 100V for 30 minutes. The 100 bps ladder was used.

Cytochrome oxidase I gene of bats was sequenced and aligned with other international sequences available in the BOLD and GenBank database using the Clustal W program. Because of the alignment species were grouped into definite clusters. These five Bats specimens were named as P. giganteus because they offered maximum sequence similarity with the BOLD (99.84 to 100%) and BLAST (100%) reference sequence of P. giganteus (KT291772.1) of Assam, India (Table 2).

Table 2
Results for identification based on BOLD and BLAST for Cytochrome Oxidase I (COI) sequencing for five Bat specimens collected from district Bagh and Muzaffarabad, Azad Kashmir.

No	Species Code	BOLD	Similarity (%)	Genbank/BLAST	Similarity (%)	Accession number
		Species identification		Species identification
1	PG-B1	Pteropus giganteus	100	Pteropus giganteus	100	KT291772.1
2	PG-B2	Pteropus giganteus	100	Pteropus giganteus	100	KT291772.1
3	PG-M1	Pteropus giganteus	99.84	Pteropus giganteus	100	KT291772.1
4	PG-M2	Pteropus giganteus	99.84	Pteropus giganteus	100	KT291772.1
5	PG-M3	Pteropus giganteus	99.84	Pteropus giganteus	100	KT291772.1

The total number of haplotypes in the sequences of P. giganteus was 2 (haplotype 1 is shared by B1 & B2 while, haplotype 2 is shared by M1, M2 & M3) and haplotype (gene) diversity was 0.60±0.17. The remarkable rate of transversional substitutions was observed as compared to transitional substitutions. The transition/transversion rate ratios are K1 = 0 (purines) and K2 = 0 (pyrimidines). The overall transition/transversion bias is R = 0, where R = [A*G*k1 + T*C*k2]/[(A+G)*(T+C)]. The estimated proportion of invariable sites was 0.0010%, and the configuration of the gamma parameter was anticipated as 200. The best-fit ML model for the CO1 gene data of P. giganteus was Jukes-Cantor (JC) as DNA evolution with the lowest Bayesian Information Criterion (BIC= 1811.843) scores.

3.1. Phylogenetic analysis of present studied sequences

Present study sequences of P. giganteus were formed two main clusters in the phylogenetic tree. One cluster was formed by the three specimens collected from Muzaffarabad district with 100% bootstrap value while the second cluster was formed by the two specimens collected from Bagh district with 100% bootstrap value (Figure 2). All the samples of P. giganteus had same nucleotide and protein sequence, the only difference is that the sequences of Bagh district (B1 & B2) had C (Proline) at position 290, whereas, all the sequences of Muzaffarabad district (M1, M2 & M3) had A (Histidine) at this position. The sequence divergence between the present studied samples of district Muzaffarabad and district Bagh is 0.2% and these sequences show the 0.001% overall nucleotide divergence (Table 3).

Figure 2
Unrooted phylogenetic tree of COI gene based on nucleotide sequences of present study.

Table 3
Table of Pairwise distance of present studied Bat samples.

	1	2	3	4	5
P. giganteus COI B1
P. giganteus COI M3	0.200
P. giganteus COI B2	0.000	0.200
P. giganteus COI M1	0.200	0.000	0.200
P. giganteus COI M2	0.200	0.000	0.200	0.000

3.2. Combined (studied + international) sequences analyses

Combined sequences of different Bat species available on NCBI GenBank resources were analyzed together with present studied samples to create the collective phylogenetic tree. These species made three major clades with more than 50% bootstrap values. All samples of P. giganteus were grouped in clade B with strong bootstrap support (59% posterior probabilities (pP), more than 59% from NJ). While the remaining sequences of Pteropus species form two sub-clusters (A1 & A2). One sub-cluster consists of Pteropus lylei with more than 50% posterior probabilities (pP), which support the 93% bootstrap values. Similarly, all the Pteropus vampyrus samples were clustered together in a vampyrus lineage with 86% and 91% bootstrap values. Clade C consists of outgroup species with 65% and 100% bootstrap values. These outgroup species were also positioned as a base of the tree (Figure 3). The mean intraspecific K2P genetic distance of studied samples was 0.001 ± 0.001. Similarly, the overall mean interspecific genetic distance was 0.016±0.003 while, the mean intergeneric genetic distance was 0.097±0.011 (Table 4).

Figure 3
Rooted phylogenetic tree of COI gene based on nucleotide sequences of the present study with international sequences (with accession number) and outgroup species.

Table 4
Pairwise genetic distance (nucleotide Kimura 2 parameter) for the COI gene sequences of the present study and international sequences.

		1	2	3	4	5	6	7	8	9	10	11	12	13	14	15	16	17	18	19
1	Pteropus giganteus COI B1
2	C. abrasus COI (JN312044.1)	0.36
3	C. paranus COI (JF447634.1)	0.6	0.03
4	C. planirostris COI(EF080319.1)	0.89	0.14	0.04
5	P. giganteus COI M3	0	0.33	0.56	0.84
6	P. giganteus CO1 (KT291772.1)	0	0.36	0.6	0.89	0.01
7	P. giganteus COI (JX500416.1)	0.04	0.15	0.33	0.57	0.03	0.04
8	P. lylei C01 (HM541305.1)	0.01	0.51	0.79	1.12	0.02	0.02	0.1
9	P. lylei CO1 (HM541302.1)	0.01	0.51	0.79	1.12	0.02	0.01	0.1	0
10	P. lylei CO1 (HM541303.1)	0.01	0.37	0.61	0.9	0.01	0.01	0.04	0.01	0.01
11	P. lylei CO1 (HM541304.1)	0.01	0.35	0.59	0.86	0.01	0.01	0.04	0.01	0.01	0.01
12	P. lylei CO1 (HM541307.1)	0.01	0.35	0.59	0.86	0.01	0.01	0.04	0.01	0.01	0.01	0
13	P. lylei CO1 (KP975224.1)	0.01	0.35	0.59	0.86	0.01	0.01	0.04	0.01	0.01	0.01	0	0
14	P. lylei CO1 (KP975225.1)	0.01	0.35	0.59	0.86	0.01	0.01	0.04	0.01	0.01	0.01	0	0	0
15	P. vampyrus CO1 (HM541308.1)	0.01	0.33	0.57	0.88	0.01	0.01	0.03	0.02	0.02	0.01	0.01	0.01	0.01	0.01
16	P. vampyrus CO1 (HM541309.1)	0.01	0.26	0.47	0.72	0.01	0.01	0.02	0.04	0.04	0.01	0	0.01	0.01	0.01	0.01
17	P. vampyrus CO1 (HM541309.1)	0.01	0.26	0.47	0.72	0.01	0.01	0.02	0.04	0.04	0.01	0.01	0.01	0.01	0.01	0.01	0.01
18	P. giganteus COI B2	0	0.36	0.6	0.89	0.01	0	0.04	0.01	0.01	0.01	0.01	0.01	0.01	0.01	0.01	0.01	0.01
19	P. giganteusCOI M1	0.01	0.33	0.56	0.84	0	0.01	0.03	0.02	0.02	0.01	0.01	0.01	0.01	0.01	0.01	0.01	0.01	0.01
20	P. giganteus COI M2	0.01	0.33	0.56	0.84	0	0.01	0.03	0.02	0.02	0.01	0.01	0.01	0.01	0.01	0.01	0.01	0.01	0.01	0

4. Discussion

In mitochondrial DNA some specific sequences are helpful in the satisfactory identification of unknown species (Dawnay et al., 2007). In contrast to other mitochondrial genes, cytochrome oxidase I gene produces very low changes in the sequence of Amino acids (Hebert et al., 2003). Mitochondrial CO1 gene consists of 600-800 nucleotides, an important part of animal DNA, and it helps in identification of global diversity of different species (Janzen et al., 2005; Savolainen et al., 2005; Dawnay et al., 2007). It is an important tool of precise documentation of several groups and shows the best way to identify unknown species of many animal classes (Elmeer et al., 2012; Ospina-Alvarez and Piferrer 2008) below the level of species and biodiversity study (Mayden et al., 2007). Surprisingly, it can also disclose secrets of unknown diversity and hidden species (Meyer and Paulay, 2005; Kerr et al., 2009). Researchers of wild biodiversity considered DNA barcoding an important means to open the new fields of research (Janzen et al., 2005; Savolainen et al., 2005; Dawnay et al., 2007).

The records of DNA barcoding also provide a system of identification where smashed and blemished materials are not identified by their morphological characters (Wong and Hanner, 2008; Victor et al., 2009). It might enable following the exotic insidious species (Ficetola et al., 2008) also give support to identification through analysis of fecal sample or content of gut (Kaartinen et al., 2010). Therefore, barcoding has a significant role in protection from poaching and illegal trade (e.g. International Trade of Endangered Species convention) but it also protects traders from advertising deceit (Wong and Hanner, 2008).

Current investigations include the phylogenetic analysis among the same species of Indian flying fox that are collected from Muzaffarabad and Bagh (AJK) and also includes the comparison of Pteropus giganteus (studied) with the international sequences of Pteropus giganteus and other closely related species present in the genus Pteropus with outgroup species used by Gager et al. (2016). The sequences for the international species are obtained from NCBI (GenBank). Our results confirmed that CO1 gene sequences of Bats are closely related with P. giganteus of BOLD (99.84-100%) and BLAST (100%) species in databanks. It is the fact that if the CO1 gene of the samples, under investigation, shows similarity below the 2% distance (Tamura et al., 2013) or less than 3% (Wong and Hanner, 2008) then the population belongs to the same species. The intra-specific sequence divergence possibility was may be due to ancestral polymorphisms and hybridization (Hajibabaei et al., 2006). The neighbor-joining (NJ) method also revealed similar relationships and authenticated that these Bat samples are P. giganteus. Similarly, Figure (3) demonstrated the phylogenetic tree of studied specimens made a single cluster with NCBI data of P. giganteus whereas, all other species ended in separate clusters. Sympatry of very similar-looking species is common in Bats (Von Helversen et al., 2001; Ibáñez et al., 2006) and it can make identification of species in the field very problematic. We clustered individuals successfully with the CO1 sequences. Besides, we incorporated GenBank sequences from different species and countries with the help of the COI tree.

The present findings offered very low (0.001%) intra-specific K2P genetic distance showing extremely low genetic variations between the species of two (Bagh and Muzaffarabad) district, as Ward et al. (2009) recommended that K2P distance between the entities of the same species should be very low. The current findings disclosed that there was minimum haplotype diversity of the COI gene along with low nucleotide diversity. Less Haplotype (gene) diversity in bats may be due to some environmental disturbances that could be a reason for the decrease of population. Furthermore, after the population becomes small, the genetic factors speed up the process of extinction for that population (Westemeier et al., 1998).

5. Conclusion

The efficacy of COI barcodes has strongly been validated by this study for identifying the Indian flying fox. The straight forward amplification and sequencing were proved by the COI barcode region, and this would enable the accurate identification of specimens and also aid the consequent generation of barcode database. It was confirmed through BOLD and BLAST that CO1 gene sequences of Bats were P. giganteus. This study concluded that P. giganteus had low G contents (19.8%) and mean intraspecific K2P genetic distance was very low, with a low rate of mutation. For the conservational benefits and to overcome the challenge for correctly describing and mapping biodiversity DNA barcodes will surely facilitate. There is an urgent need for this work because, despite high levels of genetic diversity within many species, conservationists and politicians still focus their effort around named species. There should be some healthy conservation actions within this region because habitat lost is greater already, we are confident that the DNA barcodes use and access of such information to the public with the help of Web portals will surely boost the intensified taxonomic effort needed to catalog and describe this diversity and ultimately help in its protection.

• (With 3 figures)

References

• AKHTAR, T. and ALI, G., 2016. DNA barcoding of Schizothorax species from the Neelum and Jhelum Rivers of Azad Jammu and Kashmir. Mitochondrial DNA. Part B, Resources, vol. 1, no. 1, pp. 934-936. http://dx.doi.org/10.1080/23802359.2016.1258337
  » http://dx.doi.org/10.1080/23802359.2016.1258337
• AKHTAR, T., SHAFI, N. and ALI, G., 2017. DNA bar-coding of snow trout (S. plagiostomus) from Neelum and Jhelum rivers of Azad Jammu and Kashmir, Pakistan. Inter. Journal of Biosciences, vol. 10, pp. 302-310.
• AYESHA, U.R., SHAFI, N., AKHTAR, T., ZAREEN, A. and AYUB, H., 2019. DNA barcoding of cyprinids (Labeo rohita, Catla catla and Cirrhinus mrigala), mitochondrial CO1based study. Mitochondrial DNA. Part B, Resources, vol. 4, no. 1, pp. 406-407. http://dx.doi.org/10.1080/23802359.2018.1547132
  » http://dx.doi.org/10.1080/23802359.2018.1547132
• BATES, P.J.J. and HARRISON, D.L., 1997. Bats of the Indian Subcontinent England: Harrison Zoological Museum Press, pp. 112-145.
• BLAXTER, M., ELSWORTH, B. and DAUB, J., 2004. DNA taxonomy of a neglected animal phylum: an unexpected diversity of tardigrades. Proc Roy Soc London. Biol Sci, vol. 4, pp. 189-192.
• BORISSENKO, A.V. and KRUSKOP, S.V., 2003. Bats of Vietnam and adjacent territories. An identification manual Korzun LP, Kalyakin MV, eds. Moscow, Geos, 203 p.
• CLARE, E.L., LIM, B.K., ENGSTROM, M.D., EGER, J.L. and HEBERT, P.D., 2007. DNA barcoding of Neotropical Bats: species identification and discovery within Guyana. Molecular Ecology Notes, vol. 7, no. 2, pp. 184-190. http://dx.doi.org/10.1111/j.1471-8286.2006.01657.x
  » http://dx.doi.org/10.1111/j.1471-8286.2006.01657.x
• DAWNAY, N., OGDEN, R., MCEWING, R., CARVALHO, G.R. and THORPE, R.S., 2007. Validation of the barcoding gene COI for use in forensic genetic species identification. Forensic Science International, vol. 173, no. 1, pp. 1-6. http://dx.doi.org/10.1016/j.forsciint.2006.09.013 PMid:17300895.
  » http://dx.doi.org/10.1016/j.forsciint.2006.09.013
• ELMEER, K., ALMALKI, A., MOHRAN, K.A., AL-QAHTANI, K.N. and ALMARRI, M., 2012. DNA barcoding of Oryx leucoryx using the mitochondrial cytochrome C oxidase gene. Genetics and Molecular Research, vol. 11, no. 1, pp. 539-547. http://dx.doi.org/10.4238/2012.March.8.2 PMid:22535389.
  » http://dx.doi.org/10.4238/2012.March.8.2
• FICETOLA, G.F., BONIN, A. and MIAUD, C., 2008. Population genetics reveals origin and number of founders in a biological invasion. Molecular Ecology, vol. 17, no. 3, pp. 773-782. http://dx.doi.org/10.1111/j.1365-294X.2007.03622.x PMid:18194168.
  » http://dx.doi.org/10.1111/j.1365-294X.2007.03622.x
• FRANCIS, C.M., 2008. A Field Guide to the Mammals of South-East Asia London, New Holland Publishers. 392.
• GAGER, Y., TARLAND, E., LIECKFELDT, D., MÉNAGE, M., BOTERO-CASTRO, F., ROSSITER, S.J., KRAUS, R.H., LUDWIG, A. and DECHMANN, D.K., 2016. The value of Molecular vs. Morphometric and acoustic information for species identification using sympatric Molossid Bats. PLoS One, vol. 11, no. 3, pp. e0150780. http://dx.doi.org/10.1371/journal.pone.0150780 PMid:26943355.
  » http://dx.doi.org/10.1371/journal.pone.0150780
• HAJIBABAEI, M.M., SMITH, D.H., JANZEN, J.J., RODRIGUEZ, J.B., WHITFIELD, J.B. and HEBERT, P.D.N., 2006. A minimalist barcode can identify a specimen whose DNA is degraded. Molecular Ecology Notes, vol. 6, no. 4, pp. 959-964. http://dx.doi.org/10.1111/j.1471-8286.2006.01470.x
  » http://dx.doi.org/10.1111/j.1471-8286.2006.01470.x
• HALL, T.A., 1999. BioEdit: a user-friendly biological sequence alignment editor and analysis program for windows 95/98/NT. Nucleic Acids Symposium Series, vol. 41, pp. 95-98.
• HEBERT, P.D. and GREGORY, T.R., 2005. The promise of DNA barcoding for taxonomy. Systematic Biology, vol. 54, no. 5, pp. 852-859. http://dx.doi.org/10.1080/10635150500354886 PMid:16243770.
  » http://dx.doi.org/10.1080/10635150500354886
• HEBERT, P.D., PENTON, E.H., BURNS, J.M., JANZEN, D.H. and HALLWACHS, W., 2004. Ten species in one: DNA barcoding reveals cryptic species in the neotropical skipper butterfly Astraptes fulgerator. Proceedings of the National Academy of Sciences of the United States of America, vol. 101, no. 41, pp. 14812-14817. http://dx.doi.org/10.1073/pnas.0406166101 PMid:15465915.
  » http://dx.doi.org/10.1073/pnas.0406166101
• HEBERT, P.D.N., CYWINSKA, A., BALL, S.L. and DEEWAARD, J.R., 2003. Biological identifications through DNA barcodes. Proceedings. Biological Sciences, vol. 270, no. 1512, pp. 313-321. http://dx.doi.org/10.1098/rspb.2002.2218 PMid:12614582.
  » http://dx.doi.org/10.1098/rspb.2002.2218
• IBÁÑEZ, C., GARCÍA-MUDARRA, J.L., RUEDI, M., STADELMANN, B. and JUSTE, J., 2006. The Iberian contribution to cryptic diversity in European Bats. Acta Chiropterologica, vol. 8, no. 2, pp. 277-297. http://dx.doi.org/10.3161/1733-5329(2006)8[277:TICTCD]2.0.CO;2
  » http://dx.doi.org/10.3161/1733-5329(2006)8[277:TICTCD]2.0.CO;2
• JANZEN, D.H., HAJIBABAEI, M., BURNS, J.M., HALLWACHS, W., REMIGIO, E. and HEBERT, P.D., 2005. Wedding biodiversity inventory of a large and complex Lepidoptera fauna with DNA barcoding. Philosophical transactions of the Royal Society of London. Series B, Biological sciences, vol. 360, no. 1462, pp. 1835-1845. http://dx.doi.org/10.1098/rstb.2005.1715 PMid:16214742.
  » http://dx.doi.org/10.1098/rstb.2005.1715
• KAARTINEN, R., STONE, G.N., HEARN, J., LOHSE, K. and ROSLIN, T., 2010. Revealing secret liaisons: DNA barcoding changes our understanding of food webs. Ecological Entomology, vol. 35, no. 5, pp. 623-638. http://dx.doi.org/10.1111/j.1365-2311.2010.01224.x
  » http://dx.doi.org/10.1111/j.1365-2311.2010.01224.x
• KERR, K.C., LIJTMAER, D.A., BARREIRA, A.S., HEBERT, P.D. and TUBARO, P.L., 2009. Probing evolutionary patterns in Neotropical birds through DNA barcodes. PLoS One, vol. 4, no. 2, pp. 4379. http://dx.doi.org/10.1371/journal.pone.0004379 PMid:19194495.
  » http://dx.doi.org/10.1371/journal.pone.0004379
• KUNZ, T.H., BRAUN DE TORREZ, E., BAUER, D., LOBOVA, T. and FLEMING, T.H., 2011. Ecosystem services provided by Bats. Annals of the New York Academy of Sciences, vol. 1223, no. 1, pp. 1-38. http://dx.doi.org/10.1111/j.1749-6632.2011.06004.x PMid:21449963.
  » http://dx.doi.org/10.1111/j.1749-6632.2011.06004.x
• MAYDEN, R.L., TANG, K.L., CONWAY, K.W., FREYHOF, J., CHAMBERLAIN, S., HASKINS, M., SCHNEIDER, L., SUDKAMP, M., WOOD, R.M., AGNEW, M., BUFALINO, A., SULAIMAN, Z., MIYA, M., SAITOH, K. and HE, S., 2007. Phylogenetic relationships of Danio within the order Cypriniformes: a framework for comparative and evolutionary studies of a model species. Journal of Experimental Zoology. Part B, Molecular and Developmental Evolution, vol. 308, no. 5, pp. 642-654. http://dx.doi.org/10.1002/jez.b.21175 PMid:17554749.
  » http://dx.doi.org/10.1002/jez.b.21175
• MAYER, F., DIETZ, C. and KIEFER, A., 2007. Molecular species identification boosts Bat diversity. Frontiers in Zoology, vol. 4, no. 1, pp. 4. http://dx.doi.org/10.1186/1742-9994-4-4 PMid:17295921.
  » http://dx.doi.org/10.1186/1742-9994-4-4
• MCCONKEY, K.R. and DRAKE, D.R., 2006. Flying foxes cease to function as seed dispersers long before they become rare. Ecology, vol. 87, no. 2, pp. 271-276. http://dx.doi.org/10.1890/05-0386 PMid:16637350.
  » http://dx.doi.org/10.1890/05-0386
• MCCONKEY, K.R. and DRAKE, D.R., 2007. Indirect evidence that flying foxes track food resources among islands in a Pacific archipelago. Biotropica, vol. 39, no. 3, pp. 436-440. http://dx.doi.org/10.1111/j.1744-7429.2007.00269.x
  » http://dx.doi.org/10.1111/j.1744-7429.2007.00269.x
• McConkey, KR, and Drake DR., 2015. Low redundancy in seed dispersal within an island frugivore community. AoB Plants, vol. 7, pp. 1-13.
• MCMICHAEL, L., EDSON, D., MCLAUGHLIN, A., MAYER, D., KOPP, S., MEERS, J. and FIELD, H., 2015. Haematology and Plasma Biochemistry of Wild Black Flying-Foxes, (Pteropus alecto) in Queensland, Australia. PLoS One, vol. 10, no. 5, pp. 0125741. http://dx.doi.org/10.1371/journal.pone.0125741 PMid:25938493.
  » http://dx.doi.org/10.1371/journal.pone.0125741
• MEYER, C.P. and PAULAY, G., 2005. DNA barcoding: error rates based on comprehensive sampling. PLoS Biology, vol. 3, no. 12, pp. 422. http://dx.doi.org/10.1371/journal.pbio.0030422 PMid:16336051.
  » http://dx.doi.org/10.1371/journal.pbio.0030422
• MONAGHAN, M.T., BALKE, M., GREGORY, T.R. and VOGLER, A.P., 2005. DNA-based species delineation in tropical beetles using mitochondrial and nuclear markers. Philosophical Transactions of the Royal Society of London. Series B, Biological Sciences, vol. 360, no. 1462, pp. 1925-1933. http://dx.doi.org/10.1098/rstb.2005.1724 PMid:16214750.
  » http://dx.doi.org/10.1098/rstb.2005.1724
• MORITZ, C. and CICERO, C., 2004. DNA barcoding: promise and pitfalls. PLoS Biology, vol. 2, no. 10, pp. e354. http://dx.doi.org/10.1371/journal.pbio.0020354 PMid:15486587.
  » http://dx.doi.org/10.1371/journal.pbio.0020354
• OSPINA-ALVAREZ, N. and PIFERRER, F., 2008. Temperature-dependent sex determination in fish revisited: prevalence, a single sex ratio response pattern, and possible effects of climate change. PLoS One, vol. 3, no. 7, pp. e2837. http://dx.doi.org/10.1371/journal.pone.0002837 PMid:18665231.
  » http://dx.doi.org/10.1371/journal.pone.0002837
• ROZAS, J., SÁNCHEZ-DELBARRIO, J.C., MESSEGUER, X. and ROZAS, R., 2003. DnaSP, DNA polymorphism analyses by the coalescent and other methods. Bioinformatics, vol. 19, no. 18, pp. 2496-2497. http://dx.doi.org/10.1093/bioinformatics/btg359 PMid:14668244.
  » http://dx.doi.org/10.1093/bioinformatics/btg359
• SAMBROOK, J., FRITSCHI, E. and MANIATIS, T. , 1989. Molecular cloning: a laboratory manual, vol. 1.New York: Cold Spring Harbor Laboratory Press, E5 p.
• SAVOLAINEN, V., COWAN, R.S., VOGLER, A.P., RODERICK, G.K. and LANE, R., 2005. Towards writing the encyclopaedia of life: an introduction to DNA barcoding. Philosophical Transactions of the Royal Society of London. Series B, Biological Sciences, vol. 360, no. 1462, pp. 1805-1811. http://dx.doi.org/10.1098/rstb.2005.1730 PMid:16214739.
  » http://dx.doi.org/10.1098/rstb.2005.1730
• SHEIKH, K.M. and MOLUR, S., 2004. Status and red list of Pakistan’s mammals, Based on the Conservation Assessment and Management Plan Workshop Pakistan: IUCN Pakistan.
• SMITH, M.A., WOODLEY, N.E., JANZEN, D.H., HALLWACHS, W. and HEBERT, P.D., 2006. DNA barcodes reveal cryptic host-specificity within the presumed polyphagous members of a genus of parasitoid flies (Diptera: tachinidae). Proceedings of the National Academy of Sciences of the United States of America, vol. 103, no. 10, pp. 3657-3662. http://dx.doi.org/10.1073/pnas.0511318103 PMid:16505365.
  » http://dx.doi.org/10.1073/pnas.0511318103
• TAMURA, K., DUDLEY, J., NEI, M. and KUMAR, S., 2013. MEGA4: Molecular Evolutionary Genetics Analysis (MEGA) software version 4.0. Molecular Biology and Evolution, vol. 24, no. 8, pp. 1596-1599. http://dx.doi.org/10.1093/molbev/msm092 PMid:17488738.
  » http://dx.doi.org/10.1093/molbev/msm092
• VICTOR, B.C., HANNER, R., SHIVJI, M., HYDE, J. and CALDOW, C., 2009. Identification of the larval and juvenile stages of the Cubera Snapper, Lutjanus cyanopterus, using DNA barcoding. Zootaxa, vol. 2215, pp. 24-36.
• VON HELVERSEN, O., HELLER, K.G., MAYER, F., NEMETH, A., VOLLETH, M. and GOMBKÖTÖ, P., 2001. Cryptic mammalian species: a new species of whiskered Bat (Myotis alcathoe n. sp.) in Europe. Naturwissenschaften, vol. 88, no. 5, pp. 217-223. http://dx.doi.org/10.1007/s001140100225 PMid:11482435.
  » http://dx.doi.org/10.1007/s001140100225
• WARD, R.D., HANNER, R. and HEBERT, P.D.N., 2009. The campaign to DNA barcode all fishes. Journal of Fish Biology, vol. 74, no. 2, pp. 329-356. http://dx.doi.org/10.1111/j.1095-8649.2008.02080.x PMid:20735564.
  » http://dx.doi.org/10.1111/j.1095-8649.2008.02080.x
• WESTEMEIER, R.L., BRAWN, J.D., SIMPSON, S.A., ESKER, T.L., JANSEN, R.W., WALK, J.W., KERSHNER, E.L., BOUZAT, J.L. and PAIGE, K.N., 1998. Tracking the long-term decline and recovery of an isolated population. Science, vol. 282, no. 5394, pp. 1695-1698. http://dx.doi.org/10.1126/science.282.5394.1695 PMid:9831558.
  » http://dx.doi.org/10.1126/science.282.5394.1695
• WILSON, E.O., 2003. The encyclopedia of life. Trends in Ecology & Evolution, vol. 18, no. 2, pp. 77-80. http://dx.doi.org/10.1016/S0169-5347(02)00040-X
  » http://dx.doi.org/10.1016/S0169-5347(02)00040-X
• WONG, E.H.K. and HANNER, R.H., 2008. DNA barcoding detects market substitution in North American seafood. Food Research International, vol. 41, no. 8, pp. 828-837. http://dx.doi.org/10.1016/j.foodres.2008.07.005
  » http://dx.doi.org/10.1016/j.foodres.2008.07.005

Publication Dates

• Publication in this collection
  10 Aug 2020
• Date of issue
  Jul-Sep 2021

History

• Received
  22 July 2019
• Accepted
  21 Feb 2020