The first complete mitochondrial genome sequence of Cynopterus brachyotis (Chiroptera, Pteropodidae) from the Philippines

Paul Lorenzo Gaite; Wilson Aala; Jr.; Michael Bacus; Christian Labrador; April Mae Numeron; Lief Erikson Gamalo; Lyre Anni Murao

doi:10.3897/BDJ.10.e72768

Biodiversity Data Journal : Research Article

PDF

Research Article

The first complete mitochondrial genome sequence of Cynopterus brachyotis (Chiroptera, Pteropodidae) from the Philippines

Paul Lorenzo A. Gaite^‡, Wilson F. Aala, Jr.^‡, Michael G. Bacus^‡, Christian C. Labrador^‡, April Mae M. Numeron^‡, Lief Erikson D. Gamalo^§,‡,|, Lyre Anni E. Murao^§,‡,|

‡ Philippine Genome Center Mindanao, University of the Philippines Mindanao, Davao City, Philippines

§ Department of Biological Sciences and Environmental Studies, University of the Philippines Mindanao, Davao City, Philippines

| Wildlife-Human Interaction Studies, Ecological Research and Biodiversity Conservation Laboratory, University of the Philippines Mindanao, Davao City, Philippines

Corresponding author: Lyre Anni E. Murao (lemurao@up.edu.ph)

Academic editor: Ricardo Moratelli

Received: 09 Aug 2021 | Accepted: 28 Sep 2021 | Published: 14 Feb 2022

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.

Citation: Gaite PLA, Aala, Jr. WF, Bacus MG, Labrador CC, Numeron AMM, Gamalo LED, Murao LAE (2022) The first complete mitochondrial genome sequence of Cynopterus brachyotis (Chiroptera, Pteropodidae) from the Philippines. Biodiversity Data Journal 10: e72768. https://doi.org/10.3897/BDJ.10.e72768

Abstract

The technical limitations of capillary sequencing in providing insights on phylogeny have been greatly aided in recent years by the implementation of next generation sequencing platforms which can generate whole mitochondrial genome (mitogenome) sequences. In this study, enriched mitochondrial DNA of Cynopterus brachyotis from Mindanao, Philippines was sequenced using the Illumina MiSeq platform. A total of 653,967 clean paired-end reads was assembled using a MIRA-MITObim pipeline, resulting in a consensus mitogenome sequence length of 17,382 bases and a GC content of 41.48%, which is consistent with other published mitogenomes in fruit bats. The assembled C. brachyotis mitogenome was annotated using the MITOS online server and was able to resolve all mitochondrial genes, except for one transfer RNA gene (trnT) which may be further resolved by additional capillary sequencing of the region. Sequence analysis showed that the Philippine C. brachyotis is only 90%-91% homologous with other Cynopterus spp., based on its full mitogenome sequence. Phylogenetic analysis of fruit bat mitogenomes, deposited in online repositories, revealed that the Philippine C. brachyotis in this study has diverged from Asian Cynopterus, namely Cynopterus brachyotis and Cynopterus sphinx from other parts of Asia (100% bootstrap support) with the latter two forming a separate clade. This divergence at the species level was consistent with phylogentic inference using cytochrome oxidase 1 (CO1) and cytochrome B (cytb) gene markers. Our results strengthen the previously reported hypothesis that the Cynopterus cf. brachyotis in the Philippines is distinct from its Asian counterparts and should be, therefore, treated as a new species.

Keywords

NGS, mitogenome, bats, taxonomy, phylogenetic analysis

Introduction

Cynopterus brachyotis (Müller, 1838) is a fruit bat widespread in Southeast Asia, ranging from Sri Lanka, Taiwan, Vietnam, Indonesia, Malaysia and the Philippines (Jayaraj et al. 2012). The species is believed to have reached the Philippine islands within the past 1 million years (Rowe 2018). Currently, it is abundant in different types of habitats in the country with differing degrees of disturbances, ranging from urban residential, agricultural to forested areas (Rowe 2018). Due to high population abundance, widespread geographical extent and the relatively stable population, the species is considered of least concern under the updated IUCN Red List (Csorba et al. 2019, Rosell-Ambal et al. 2019).

The status of C. brachyotis taxonomy in the Philippines is uncertain and challenging. Genetic studies suggest that the Philippine population has high levels of divergence in comparison to other C. brachyotis populations in other Southeast Asian countries (Bacus et al. 2021, Luczon et al. 2019). Moreover, other authorities suggested that the Philippine Cynopterus should be assigned as Cynopterus luzoniensis, a taxon believed to exist throughout the Philippine archipelago and in the Sulawesi Island of Indonesia (Campbell et al. 2004, Rosell-Ambal et al. 2019, Rowe 2018, Simmons 2005). However, these hypotheses should be further evaluated as the exact identification of C. luzoniensis in Sulawesi may be different from "C. luzoniensis" in the Philippines or in some islands in the country (Campbell et al. 2004, Jayaraj et al. 2012). Due to this confusion on the distribution and extent of the taxa mentioned above, we will be consistently using "Cynopterus cf. brachyotis" in this paper to represent the Philippine Cynopterus to avoid confusion until taxonomy is resolved.

The analyses of mitochondrial genes are crucial to understand the challenging and confusing taxonomy and evolutionary history of bats in the Philippines. However, information on the mitochondrial genome (mitogenome) of Philippine bats is clearly fragmented despite the wide application of this data for nomenclature, evolution and phylogeny of many species. Fortunately, the technical limitations of capillary sequencing to provide phylogenetic evidence have been greatly aided by the use of Next Generation Sequencing (NGS) platforms that can generate whole mitogenome sequences.

As part of the larger initiative to sequence the mitochondrial DNA of native species in the country, we provide the complete mitochondrial genome (mitogenome) of Cynopterus cf. brachyotis from Mindanao Island. Moreover, phylogenetic trees, based on the sequences generated from this study and from those deposited in GenBank, were provided to show the evolutionary relationship of the Philippine population to other Cynopterus spp. from other countries in Asia. Results from this study will contribute to the taxonomic status and identity of the Philippine Cynopterus which, in turn, can be used as additional supporting information for other phylogenetic, evolution and ecological studies of bat species in the world.

Material and methods

Samples

The sample used in this study was a cardiac tissue from a voucher specimen, previously identified as Cynopterus brachyotis using morphological characteristics and measurements (e.g. forearm length, head length, body length, tail length and ear length) with the aid of the taxonomic key developed by Ingle and Heaney (1992). Moreover, a ~ 202 bp barcode for the mitochondrial cytochrome oxidase 1 (CO1) gene (Walker et al. 2016) was determined in a separate study (Bacus et al. 2021) that further confirmed the identity of the specimen. An IACUC approval from the University of the Philippines Manila (Protocol No.: 2018-109) and a Gratuitous permit from the Philippine Department of Environment and Natural Resources Region XI (WGP No.: 2018-07) were obtained for the collection of voucher specimens.

DNA Extraction and Bat Mitogenome Enrichment

Around 20 mg of cardiac tissue was obtained from a voucher specimen of C. brachyotis. Genomic DNA was extracted using the DNEasy Blood and Tissue kit (Qiagen, Hilden, Germany) according to the manufacturer’s instructions. Mitogenome enrichment was performed using two primer pairs, namely LR6F (5’-GCC CAT ACC CCG AAA ATG TTG-3’), LR6R (5’-CGG CGG GAG AAG TAG ATT GAA-3’), LR15F (5’-CCA CAG GAA AAT CAG CCC AAT T-3’) and LR15R (5'-GCT GTT GCT GTG TCT GAT GTG-3’) as previously described (Jebb et al. 2017). PCR amplicons were visualised with 0.7% agarose gel and stored at -20°C prior to the library preparation step.

Library Preparation and Illumina Sequencing

The two PCR amplicons were quantified using Qubit dsDNA HS assay and normalised in a resuspension buffer to a final volume of 2.5 µl and a final concentration of 10 ng/µl (25 ng total yield). Both amplicons were pooled in a new PCR tube. DNA tagmentation, amplification and library clean-up were performed using the Nextera XT Library Prep Kit according to the manufacturer’s instructions. Cleaned libraries were then sent to the Philippine Genome Center Diliman-DNA Sequencing Core Facility for normalisation and next-generation sequencing using the Illumina MiSeq platform and a paired-end read format at 2 x 150 bp for 500 cycles.

Mitogenome Assembly and Analysis

Raw sequencing reads were subjected to initial quality assessment with the FastQC software (Andrews 2010). Raw reads were trimmed of adapters and low-quality bases at both ends using Trimmomatic (Bolger et al. 2014). The filtered reads were used to generate an initial mapping assembly using MIRA (Chevreux et al. 2013) and were iteratively mapped back to the initial mapping assembly and assembled using MITObim (Hahn et al. 2013). The reference sequence that was used for the assembly was the NCBI Reference Sequence entry for Cynopterus brachyotis mitochondrion, complete genome (GenBank Accession ID: NC_026465). The different genes/regions in the resulting assembly were annotated using the MITOS web server (Bernt et al. 2013). Both the entire assembled genome and the identified genes/regions were annotated individually using the nucleotide BLAST tool at the NCBI website (NCBI Resource Coordinators et al. 2018). Default settings were used in the parameters for running all of the aforementioned procedures. The assembled mitogenome sequence was submitted to NCBI GenBank and subsequently assigned the accession number OK163841.

Phylogenetic Analysis

The assembled mitogenome sequence of Cynopterus brachyotis from Mindanao, Philippines was aligned with available fruit bat mitogenome sequences in the NCBI GenBank. The best phylogenetic model was determined using MEGA 7, based on the lowest BIC and AIC values (Kumar et al. 2016). A phylogenetic tree was subsequently inferred using the GTR+G+I DNA substitution and site heterogeneity model using a Maximum-Likelihood estimation with 1000 bootstrap replicates in MEGA 7 (Kumar et al. 2016). A similar analysis was performed for the full cytochrome oxidase 1 (CO1) and cytochrome B (cytb) genes. For the CO1 gene, a multiple sequence alignment of 1,533 bp was used for the phylogenetic analysis using the Maximum-Likelihood method with the HKY+G+I model in MEGA 7 with 1000 bootstrap replicates. For the cytb gene, a multiple sequence alignment of 1,134 bp was used for the phylogenetic analysis using the Maximum-Likelihood method with the GTR+G+I model in MEGA 7 with 1000 bootstrap replicates. A distant vertebrate group, fish, was selected as the outgroup.

Results

Mitogenome

A total of 2,180,947 reads was initially generated by next generation sequencing of the mitochondrial DNA-enriched tissue sample (Table 1). The quality control steps reduced this number of reads to 1,784,834. From this readpool, a total of 653,967 reads was used for the assembly of the enriched mitochondrial DNA by the MIRA-MITObim workflow. The generated mitochondrial DNA assembly had a total length of 17,382 bases with a GC content of 41.48%. The average depth of coverage of this assembly was calculated to be 10,158.30.

Table 1.

Download as

CSV

XLSX

Summary of mitochondrial genome assembly statistics.

Total number of generated raw sequencing reads	2,180,947
Total number of reads after quality control processing	1,784,834
Total number of reads used for mitogenome assembly	653,967
Total length of assembled mitogenome sequence	17,382
% GC of assembled mitogenome sequence	41.48%
Average depth of coverage of assembled mitogenome	10,158.30

The assembled sequence was subjected to online nucleotide BLAST search, with the top two hits being Cynopterus sphinx voucher MZF1967 mitochondrion, complete genome (90.86% identity, 96% query coverage, 82 sequence gaps) and Cynopterus brachyotis mitochondrion, complete genome (90.56% identity, 95% query coverage, 83 sequence gaps) (Table 2). Further annotation of the assembled sequence by MITOS identified a total of 36 genes: 21 tRNA genes, 2 rRNA genes and 13 protein-coding genes (Fig. 1). All known vertebrate mitochondrial genes have been found in this mitogenome, except trnT. Performing online nucleotide BLAST on individual gene regions identified by MITOS resulted in mostly either C. brachyotis or C. sphinx top hits with percent identities within the range of 90%-97%. It is also noted that the top hit for trnM(atg) has no significant similarity found with any NCBI-BLAST database entry, a few top hits are non-Cynopterus related entries (e.g. Megaerops niphanae and Myonycteris torquata) and the top hits for trnY(tac) and trnD(gac) are non-bat related entries, such as Maxomys surifer and Equus asinus, respectively (Table 3).

Table 2.

Download as

CSV

XLSX

Top five NCBI-BLAST hits of the assembled mitogenome sequence.

BLAST hit description	query cover	e-value	% identity	accession ID
Cynopterus sphinx voucher MZF1967 mitochondrion, complete genome	96%	0	90.86	MT259587
Cynopterus brachyotis mitochondrion, complete genome	95%	0	90.56	KM659865
Cynopterus brachyotis isolate 140810x2 mitochondrion, complete genome	95%	0	91.29	MN816304
Cynopterus sphinx isolate CKM35 mitochondrion, complete genome	96%	0	90.94	NC_046902
Megaerops niphanae isolate NLN25 mitochondrion, complete genome	96%	0	88.27	NC_046915

Table 3.

Download as

CSV

XLSX

Top NCBI-BLAST hits for the MITOS-annotated mitochondrial genes.

mitochondrial gene (MITOS)	top hit annotation	query cover	e-value	% identity	accession ID
trnF(ttc)	Cynopterus sphinx voucher MZF1967 mitochondrion, complete genome	100%	2.00E-22	97.1	MT259587.1
rrnS	Cynopterus sphinx voucher MZF1967 mitochondrion, complete genome	99%	0	94.21	MT259587.1
trnV	Cynopterus sphinx voucher MZF1967 mitochondrion, complete genome	100%	8.00E-26	100	MT259587.1
rrnL	Cynopterus brachyotis isolate 140810x2 mitochondrion, complete genome	100%	0	96.31	MN816304.1
trnL2(tta)	Cynopterus sphinx voucher MZF1967 mitochondrion, complete genome	96%	6.00E-28	100	MT259587.1
nad1	Cynopterus sphinx mitochondrial ND1 gene for NADH dehydrogenase subunit 1, complete cds	100%	0	91.27	AB079802.1
trnI(atc)	Megaerops niphanae isolate NLN25 mitochondrion, complete genome	100%	1.00E-24	98.55	NC_046915.1
trnQ(caa)	Cynopterus sphinx voucher MZF1967 mitochondrion, complete genome	100%	2.00E-23	95.89	MT259587.1
trnM(atg)	(no significant similarity found)	-	-	-	-
nad2	Cynopterus brachyotis isolate 140810x2 mitochondrion, complete genome	97%	0	89.33	MN816304.1
trnW(tga)	Cynopterus brachyotis isolate 140810x2 mitochondrion, complete genome	100%	1.00E-23	97.14	MN816304.1
trnA(gca)	Myonycteris torquata isolate 1732 mitochondrion, complete genome	100%	2.00E-26	100	NC_046922.1
trnN(aac)	Cynopterus sphinx isolate CKM35 mitochondrion, complete genome	100%	2.00E-18	91.78	NC_046902.1
trnC(tgc)	Cynopterus sphinx isolate CKM35 mitochondrion, complete genome	98%	9.00E-21	95.59	NC_046902.1
trnY(tac)	Maxomys surifer voucher CBGP_R4223 mitochondrion, complete genome	100%	5.00E-23	98.51	KY464181.1
cox1	Cynopterus brachyotis mitochondrion, complete genome	100%	0	91	KM659865.1
trnS2(tca)	Cynopterus sphinx voucher MZF1967 mitochondrion, complete genome	100%	2.00E-26	100	MT259587.1
trnD(gac)	Equus asinus isolate 9348 mitochondrion, complete genome	98%	1.00E-24	100	MK896308.1
cox2	Cynopterus brachyotis isolate 140810x2 mitochondrion, complete genome	100%	0	91.48	MN816304.1
trnK(aaa)	Rousettus amplexicaudatus voucher MRMD 2101 mitochondrion, complete genome	100%	4.00E-19	92.96	NC_045044.1
atp8	Cynopterus brachyotis mitochondrion, complete genome	100%	2.00E-77	94.36	KM659865.1
atp6	Cynopterus brachyotis mitochondrion, complete genome	100%	0	90.96	KM659865.1
cox3	Cynopterus brachyotis isolate 140810x2 mitochondrion, complete genome	100%	0	93.1	MN816304.1
trnG(gga)	Megaerops niphanae isolate NLN25 mitochondrion, complete genome	100%	2.00E-22	97.06	NC_046915.1
nad3	Cynopterus sphinx voucher MZF1967 mitochondrion, complete genome	100%	8.00E-129	91.3	MT259587.1
trnR(cga)	Megaerops niphanae isolate NLN25 mitochondrion, complete genome	100%	4.00E-24	98.55	NC_046915.1
nad4l	Cynopterus brachyotis mitochondrion, complete genome	100%	1.00E-105	91.58	KM659865.1
nad4	Cynopterus sphinx isolate CKM35 mitochondrion, complete genome	100%	0	92.25	NC_046902.1
trnH(cac)	Cynopterus sphinx voucher MZF1967 mitochondrion, complete genome	100%	4.00E-24	98.53	MT259587.1
trnS1(agc)	Cynopterus sphinx voucher MZF1967 mitochondrion, complete genome	100%	6.00E-21	100	MT259587.1
trnL1(cta)	Rousettus amplexicaudatus voucher MZF1072 mitochondrion, complete genome	100%	7.00E-27	100	MT259592.1
nad5	Cynopterus brachyotis isolate 140810x2 mitochondrion, complete genome	99%	0	89.88	MN816304.1
nad6	Cynopterus brachyotis isolate 140810x2 mitochondrion, complete genome	100%	0	89.83	MN816304.1
trnE(gaa)	Cynopterus brachyotis mitochondrion, complete genome	100%	2.00E-21	95.65	KM659865.1
cob	Cynopterus brachyotis mitochondrial cytb gene for cytochrome b, complete cds, country: Philippines: Davao del Norte, Panabo, Tibungol	100%	0	87.59	AB046321.1
trnP(cca)	Megaerops niphanae isolate NLN25 mitochondrion, complete genome	100%	1.00E-19	95.52	NC_046915.1
trnF(ttc)	Cynopterus sphinx voucher MZF1967 mitochondrion, complete genome	100%	2.00E-22	97.1	MT259587.1

Figure 1.

Graphical representation of the assembled Cynopterus cf. brachyotis mitochondrial genome with MITOS-annotated genes (figure generated with Circos software by Krzywinski et al. 2009).

Phylogenetics

A phylogenetic tree, based on the complete mitogenome sequences of selected fruit bats, was generated (Fig. 2). The cynopterine fruit bats formed a monophyletic clade with 100% bootstrap support. Interestingly, the Cynopterus brachyotis from Mindanao, Philippines appeared to have diverged from the rest of the Asian Cynopterus (C. brachyotis & C. sphinx) with a robust bootstrap support (100%). This branching pattern also appeared to be consistent with phylogenetic trees generated using the two commonly-used markers, CO1 and cytb genes (Fig. 3).

Figure 2.

Maximum Likelihood phylogeny of selected fruit bat mitogenomes obtained from GenBank. The tree was generated using MEGA 7 with 1000 bootstrap replicates. The DNA substitution model used was GTR+G+I. The scale bar denotes the number of substitutions per site which is reflected in the branch lengths (i.e. 0.5 value means 5 nucleotide changes for every 10 nucleotides). Bootstrap values for clade support are written on the nodes of the branches. An outgroup consisting of fish mitogenomes was included in the phylogenetic tree inference.

Figure 3.

Maximum Likelihood phylogeny of fruit bat A) CO1 and B) cytb gene obtained from GenBank. The tree was generated using MEGA 7 with 1000 bootstrap replicates. For the CO1 gene, a multiple sequence alignment of 1,533 bp was used for the phylogenetic analysis using the Maximum-Likelihood method with the HKY+G+I model in MEGA 7 with 1000 bootstrap replicates. For the cytb gene, a multiple sequence alignment of 1,134 bp was used for the phylogenetic analysis using the Maximum-Likelihood method with the GTR+G+I model in MEGA 7 with 1000 bootstrap replicates. The scale bar denotes the number of substitutions per site which is reflected in the branch lengths (i.e. 0.5 value means 5 nucleotide changes for every 10 nucleotides). Bootstrap values for clade support are written on the nodes of the branches. An outgroup consisting of fish mitogenomes was included in the phylogenetic tree inference.

Discussion

An almost complete mitogenome (99.86%, contains a single 25-base "N" gap in the final assembled sequence) was generated using an Illumina NGS and MIRA-MITObim assembly workflow. Nucleotide BLAST of the assembled mitogenome and individual genes confirmed its identity to be Cynopterus sp., the first reported mitogenome of this genus from the Philippines. Further annotation of individual mitochondrial genes by MITOS determined that almost all vertebrate mitochondrial genes were present on the assembled mitogenome, except for trnT. The length (17,382 bases) and GC content (41.48%) of the assembled mitogenome were similar and hence concordant with other published fruit bat mitogenomes (Jahari et al. 2020, Jahari et al. 2021, Yue 2019). These steps and resulting statistics confirm that a true, high-quality bat mitogenome has been assembled from the sequence reads and the identity of this assembled mitogenome can be attributed to Cynopterus sp.. However, the relatively low similarity of the whole mitogenome (90.86% and 90.56% to C. sphinx and C. brachyotis, respectively) and individual gene regions (90%-97%) to both C. brachyotis and C. sphinx entries suggest that the assembled mitogenome originated from a species that is genetically unique from the aforementioned. The presence of non-Cynopterus and non-bat BLAST hit results and even a hit result with no significant similarity, for a number of the mitochondrial genes due to low similarity to Cynopterus-related database entries, provide further evidence to this proposition.

Phylogenetic analysis of the mitogenomes and individual CO1 and cytb genes confirmed that the sequenced bat in this study formed a genetically-distinct clade from the Asian C. brachyotis and C. sphinx. Similar findings have been reported in recent years for the Philippine Cynopterus cf. brachyotis through sequence analysis of the standard (500 bp) and even shorter (200 bp) CO1, cytb and control region (Bacus et al. 2021, Campbell et al. 2004, Luczon et al. 2019). In fact, the standard CO1 barcode depicted a high genetic distance (11.5%) between the Philippine and Southeast Asian C. brachyotis (Luczon et al. 2019), supporting its classification as a distinct species (Francis et al. 2010). Confirmation of these findings through our mitogenome analysis further strengthens this proposition.

Cynopterus brachyotis is hypothesised to be a species complex, based on sequence analysis of the mitochondrial cytb and control region sequences (Campbell et al. 2004). In the Philippines, two distinct lineages are present (Campbell et al. 2004), suggesting multiple colonisation events from multiple origins. The Sunda C. brachyotis lineage is present in Palawan Island and the rest of the Philippines form a distinct lineage which is proposed for taxonomic reassignment to "C. luzoniensis". However, this reassignment is not consistently accepted by some authors. Thus, we recommend further investigation using more genetic samples from different islands in the Philippines that includes representatives from the Sulawesi population to fully resolve the taxonomic issues revolving around the Philippine Cynopterus. Moreover, inconsistencies on the taxonomic identity involving this species complex suggests the importance of including morphological and ecological data in resolving taxonomic problems (Campbell et al. 2004).

Phylogenetic analysis, using both whole mitogenome and individual gene markers, such as CO1 and cytb, produced similar phylogenetic trees with high branch support values, suggesting that these gene markers can be used to reliably resolve taxonomic groups in fruit bats. While this holds true at the genus level, careful consideration should be given when using single gene markers for species delineation as exemplified by the discordance between the CO1 and cytb phylogenetic trees for the placement of C. sphinx within the cynopterine clade. Such contrasting branching patterns might arise especially when only a few representative samples are used in the analysis.

This study focused on the mitochondrial genome due to its utility in taxonomic resolution. Mitochondrial DNA evolves faster compared to nuclear genetic markers, generating a higher degree of sequence variation which makes it useful to resolve lower taxonomic levels for organisms (Chan et al. 2021, Hassanin et al. 2013). Furthermore, the relatively small size of the mitogenome (~ 16,000 base pairs long) makes it a practical alternative as a molecular marker. Despite these characteristics, the sole use of the mitochondrial genome in taxonomy remains to have certain limitations. Alternatively, nuclear markers represent neutral and paternal inheritance. In particular, single nucleotide polymorphisms (SNPs) and microsatellites exhibit high polymorphism amongst populations or individuals (Galbusera et al. 2000, Jiang et al. 2016). In addition, the highly-conserved nature of nuclear rRNA gene sequences makes it suitable for resolving higher taxonomic levels for organisms (Chan et al. 2021). Other studies have shown that combining mitochondrial and nuclear DNA analysis provided better taxonomic resolution (Abreu et al. 2020, Lambret-Frotté et al. 2012). A similar approach can be used in future studies to provide a more comprehensive and accurate perspective to ascertain the taxonomic resolution of the Mindanao or Philippine Cynopterus cf. brachyotis.

Conclusions

To date, the discussion on Cynopterus taxonomy in the Philippines is still fraught with conflicting evidence and, as a result, perspectives as well. Hence, this study was conducted so as to address part of these conflicts by attempting to provide further genetic evidence on Philippine Cynopterus through the sequencing of the mitogenome of Cynopterus brachyotis. This study was, indeed, able to generate an assembled Cynopterus brachyotis mitogenome from PCR-enriched bat tissue DNA, which was subsequently proven to be a high-quality and complete mitogenome sequence as determined through further quality assessment and annotation steps. The annotation steps revealed the relatively low similarity of the assembled mitogenome and individual genes with C. brachyotis and C. sphinx mitochondrial genomes and genes, respectively, suggesting that the mitogenome originated from a species that is genetically different to C. brachyotis (as originally identified); hence, phylogenetic analysis was performed for further confirmation. Phylogenetic analysis of whole mitogenome and CO1/cytb genes revealed and confirmed divergence of the Philippine Cynopterus cf. brachyotis from Asian Cynopterus (C. brachyotis and C. sphinx) with high branch support values. The results of this study altogether provide further evidence to confirm a previous proposition that Philippine Cynopterus cf. brachyotis is a unique species. Although part of the conflicting discussion on Philippine Cynopterus has been resolved by the results of this study, future studies to obtain nuclear marker/genome data, complementing the mitogenome data as well as morphological, behaviour and ecological data, are still needed so as to further establish the taxonomic resolution and the geographical extent of the Philippine Cynopterus cf. brachyotis.

Acknowledgements

This study is an output of the Philippine Genome Center Mindanao, which is funded by the Philippine Department of Science and Technology (DOST) with the support of the Philippine Council for Health Research and Development of DOST (DOST-PCHRD) as a monitoring agency. The study was also co-funded by the University of the Philippines Mindanao In-House Research Grant. The team would also like to extend its gratitude to the Philippine Department of Environment and Natural Resources (DENR) - Region XI for granting the researchers a Gratuitous permit to collect biological samples for research and to the Philippine Local Government Unit of Malagos for its support.

Hosting institution

Philippine Genome Center Mindanao, University of the Philippines Mindanao, Davao City, Philippines

Ethics and security

IACUC approval from the University of the Philippines Manila (Protocol No.: 2018-109)

Gratuitous permit from the Philippine Department of Environment and Natural Resources (DENR) - Region XI (WGP No.: 2018-07)

Conflicts of interest

The authors declare no conflicts of interest.

References

Abreu J, Salvi D, Perera A, Harris DJ (2020)

Mitochondrial lineages or discrete species? Assessing diversity within Timon tangitanus (Lacertidae) using mitochondrial and nuclear DNA sequences

Amphibia-Reptilia

(

123

‑

132

. https://doi.org/10.1163/15685381-20191186

Andrews S (2010)

FastQC: A quality control tool for high throughput sequence data

. http://www.bioinformatics.babraham.ac.uk/projects/fastqc/

Bacus MG, Burgos SC, Elizagaque HG, Malbog KMA, Responte MA, Gamalo LED, Achondo MJMM, Murao LAE (2021)

Pilot fecal DNA barcoding on selected fruit bats in Davao City, Philippines

Philippine Journal of Science

150

(

545

‑

555

Bernt M, Donath A, Jühling F, Externbrink F, Florentz C, Fritzsch G, Pütz J, Middendorf M, Stadler P (2013)

MITOS: Improved de novo metazoan mitochondrial genome annotation

Molecular Phylogenetics and Evolution

(

313

‑

319

. https://doi.org/10.1016/j.ympev.2012.08.023

Bolger A, Lohse M, Usadel B (2014)

Trimmomatic: a flexible trimmer for Illumina sequence data

Bioinformatics

(

2114

‑

2120

. https://doi.org/10.1093/bioinformatics/btu170

Campbell P, Schneider C, Adnan A, Zubaid A, Kunz T (2004)

Phylogeny and phylogeography of Old World fruit bats in the Cynopterus brachyotis complex

Molecular Phylogenetics and Evolution

(

764

‑

781

. https://doi.org/10.1016/j.ympev.2004.06.019

Chan AHE, Chaisiri K, Saralamba S, Morand S, Thaenkham U (2021)

Assessing the suitability of mitochondrial and nuclear DNA genetic markers for molecular systematics and species identification of helminths

Parasites & Vectors

(

233

‑

. https://doi.org/10.1186/s13071-021-04737-y

Chevreux B, Weber J, Hörster A, Dlugosch K (2013)

Sequence assembly with MIRA 4

. http://mira-assembler.sourceforge.net/docs/

Csorba G, Bumrungsri S, Bates P, Gumal M, Kingston T, Molur S, Srinivasulu C (2019)

IUCN Red List of threatened species: Cynopterus brachyotis

The IUCN Red List of Threatened Species

https://doi.org/10.2305/IUCN.UK.2019-3.RLTS.T6103A22113381.en

Francis C, Borisenko A, Ivanova N, Eger J, Lim B, Guillén-Servent A, Kruskop S, Mackie I, Hebert PN (2010)

The role of DNA barcodes in understanding and conservation of mammal diversity in Southeast Asia

PLOS One

(

e12575

. https://doi.org/10.1371/journal.pone.0012575

Galbusera P, van Dongen S, Matthysen E (2000)

Cross-species amplification of microsatellite primers in passerine birds

Conservation Genetics

(

163

‑

168

. https://doi.org/10.1023/A:1026587024065

Hahn C, Bachmann L, Chevreux B (2013)

Reconstructing mitochondrial genomes directly from genomic next-generation sequencing reads— a baiting and iterative mapping approach

Nucleic Acids Research

(

e129

. https://doi.org/10.1093/nar/gkt371

Hassanin A, An J, Ropiquet A, Nguyen TT, Couloux A (2013)

Combining multiple autosomal introns for studying shallow phylogeny and taxonomy of Laurasiatherian mammals: Application to the tribe Bovini (Cetartiodactyla, Bovidae)

Molecular Phylogenetics and Evolution

(

766

‑

775

. https://doi.org/10.1016/j.ympev.2012.11.003

Ingle NR, Heaney LR (1992)

A key to the bats of the Philippine Islands

Publication (USA)

1440

‑

. URL: https://agris.fao.org/agris-search/search.do?recordID=US9326218

Jahari PNS, Mohd Azman S, Munian K, M. Fauzi NF, Shamsir MS, Richter S, Mohd Salleh F (2020)

The first complete mitochondrial genome data of Geoffroy's rousette, Rousettus amplexicaudatus originating from Malaysia

Mitochondrial DNA Part B

(

3262

‑

3264

. https://doi.org/10.1080/23802359.2020.1812449

Jahari PNS, Mohd Azman S, Munian K, Zakaria NA, Omar MSS, Richter S, Mohd Salleh F (2021)

The first mitochondrial genome data of an Old World fruit bat, Cynopterus sphinx from Malaysia

Mitochondrial DNA Part B

(

‑

. https://doi.org/10.1080/23802359.2020.1846472

Jayaraj VK, Laman CJ, Abdullah MT (2012)

A predictive model to differentiate the fruit bats Cynopterus brachyotis and C. cf. brachyotis Forest (Chiroptera: Pteropodidae) from Malaysia using multivariate analysis

Zoological Studies

(

259

‑

271

Jebb D, Foley N, Puechmaille S, Teeling E (2017)

The complete mitochondrial genome of the greater mouse-eared bat, Myotis myotis (Chiroptera: Vespertilionidae)

Mitochondrial DNA Part A

(

347

‑

349

. https://doi.org/10.3109/19401736.2015.1122775

Jiang J, Yu J, Li J, Li P, Fan Z, Niu L, Deng J, Yue B, Li J (2016)

Mitochondrial genome and nuclear markers provide new insight into the evolutionary history of macaques

PLOS One

(

e0154665

. https://doi.org/10.1371/journal.pone.0154665

Krzywinski M, Schein J, Birol I, Connors J, Gascoyne R, Horsman D, Jones SJ, Marra MA (2009)

Circos: An information aesthetic for comparative genomics

Genome Research

(

1639

‑

1645

. https://doi.org/10.1101/gr.092759.109

Kumar S, Stecher G, Tamura K (2016)

MEGA7: Molecular evolutionary genetics analysis version 7.0 for bigger datasets

Molecular Biology and Evolution

(

1870

‑

1874

. https://doi.org/10.1093/molbev/msw054

Lambret-Frotté J, Perini FA, de Moraes Russo CA (2012)

Efficiency of nuclear and mitochondrial markers recovering and supporting known amniote groups

Evolutionary Bioinformatics

463

‑

473

. https://doi.org/10.4137/EBO.S9656

Luczon AU, Ampo SAMM, Roño JGA, Duya MRM, Ong PS, Fontanilla IKC (2019)

DNA barcodes reveal high genetic diversity in Philippine fruit bats

Philippine Journal of Science

148

(

133

‑

140

NCBI Resource Coordinators, Agarwala R, Barrett T, Beck J, Benson DA, Bollin C, Bolton E, Bourexis D, Brister JR, Bryant SH, Canese K, Cavanaugh M, Charowhas C, Clark K, Dondoshansky I, Feolo M, Fitzpatrick L, Funk K, Geer LY, Gorelenkov V, Graeff A, Hlavina W, Holmes B, Johnson M, Kattman B, Khotomlianski V, Kimchi A, Kimelman M, Kimura M, Kitts P, Klimke W, Kotliarov A, Krasnov S, Kuznetsov A, Landrum MJ, Landsman D, Lathrop S, Lee JM, Leubsdorf C, Lu Z, Madden TL, Marchler-Bauer A, Malheiro A, Meric P, Karsch-Mizrachi I, Mnev A, Murphy T, Orris R, Ostell J, O'Sullivan C, Palanigobu V, Panchenko AR, Phan L, Pierov B, Pruitt KD, Rodarmer K, Sayers EW, Schneider V, Schoch CL, Schuler GD, Sherry ST, Siyan K, Soboleva A, Soussov V, Starchenko G, Tatusova TA, Thibaud-Nissen F, Todorov K, Trawick BW, Vakatov D, Ward M, Yaschenko E, Zasypkin A, Zbicz K (2018)

Database resources of the National Center for Biotechnology Information

Nucleic Acids Research

(

‑

D13

. https://doi.org/10.1093/nar/gkx1095

Rosell-Ambal G, Kingston T, Maryanto I (2019)

IUCN Red List of threatened species: Cynopterus luzoniensis

The IUCN Red List of Threatened Species

https://doi.org/10.2305/IUCN.UK.2019-3.RLTS.T136798A22035092.en

Rowe K (2018)

The mammals of Luzon Island: Biogeography and natural history of a Philippine fauna. By Lawrence R. Heaney, Danilo S. Balete, and Eric A. Rickart; illustrated by Velizar Simeonovski and Andria Niedzielski. Baltimore (Maryland): Johns Hopkins University Press. $75.00. xi + 287 p.; ill.; index. ISBN: 978-1-4214-1837-7 (hc); 978-1-4214-1838-4 (eb). 2016.

The Quarterly Review of Biology

(

‑

. https://doi.org/10.1086/696789

Simmons NB (2005)

Order Chiroptera

Mammal Species of the World: A Taxonomic and Geographic Reference

312

‑

529

. URL: https://ci.nii.ac.jp/naid/20001149042

Walker F, Williamson CD, Sanchez D, Sobek C, Chambers C (2016)

Species from feces: Order-wide identification of Chiroptera from guano and other non-invasive genetic samples

PLOS One

(

e0162342

. https://doi.org/10.1371/journal.pone.0162342

Yue X (2019)

The complete mitochondrial genome of an old world fruit bat Sphaerias blanfordi (Pteropodidae)

Mitochondrial DNA Part B

(

1761

‑

1762

. https://doi.org/10.1080/23802359.2019.1611392

Supplementary material

Endnotes