PMCCPMCCPMCC

Search tips
Search criteria 

Advanced

 
Logo of peerjLatest ArticlesFor AuthorsEditorial BoardPeerJPeerJ
 
PeerJ. 2017; 5: e3176.
Published online 2017 April 12. doi:  10.7717/peerj.3176
PMCID: PMC5391789

Elucidating the diet of the island flying fox (Pteropus hypomelanus) in Peninsular Malaysia through Illumina Next-Generation Sequencing

Academic Editor: David Coltman

Abstract

There is an urgent need to identify and understand the ecosystem services of pollination and seed dispersal provided by threatened mammals such as flying foxes. The first step towards this is to obtain comprehensive data on their diet. However, the volant and nocturnal nature of bats presents a particularly challenging situation, and conventional microhistological approaches to studying their diet can be laborious and time-consuming, and provide incomplete information. We used Illumina Next-Generation Sequencing (NGS) as a novel, non-invasive method for analysing the diet of the island flying fox (Pteropus hypomelanus) on Tioman Island, Peninsular Malaysia. Through DNA metabarcoding of plants in flying fox droppings, using primers targeting the rbcL gene, we identified at least 29 Operationally Taxonomic Units (OTUs) comprising the diet of this giant pteropodid. OTU sequences matched at least four genera and 14 plant families from online reference databases based on a conservative Least Common Ancestor approach, and eight species from our site-specific plant reference collection. NGS was just as successful as conventional microhistological analysis in detecting plant taxa from droppings, but also uncovered six additional plant taxa. The island flying fox’s diet appeared to be dominated by figs (Ficus sp.), which was the most abundant plant taxon detected in the droppings every single month. Our study has shown that NGS can add value to the conventional microhistological approach in identifying food plant species from flying fox droppings. At this point in time, more accurate genus- and species-level identification of OTUs not only requires support from databases with more representative sequences of relevant plant DNA, but probably necessitates in situ collection of plant specimens to create a reference collection. Although this method cannot be used to quantify true abundance or proportion of plant species, nor plant parts consumed, it ultimately provides a very important first step towards identifying plant taxa and spatio-temporal patterns in flying fox diets.

Keywords: Amplicon, rbcL, Pteropodid, Fruit bat, Metabarcoding, Phytophagous, Frugivory, Nectarivory

Introduction

Understanding the contribution of animals to the functioning of rainforests has become an important issue in conservation biology. Conservation studies are now recognizing the need to collect qualitative and quantitative information on trophic relationships between animals and plants, not only to identify potential ecosystem service providers (Pompanon et al., 2012; Hibert et al., 2013), but also to inform management interventions for threatened species (Valentini et al., 2009; Ando et al., 2013).

Bats (Order: Chiroptera) provide important ecosystem services such as insect pest suppression, pollination, and seed dispersal (Fujita & Tuttle, 1991; Kunz et al., 2011). Characterising their diet is a fundamental step towards understanding their ecological roles. In the Old World, fruit bats such as flying foxes (Pteropodidae: Pteropus spp., Acerodon spp.; Kingston, 2010) have become increasingly threatened by hunting for bushmeat and medicine (Mildenstein, Tanshi & Racey, 2016). Identifying their diet and roles as ecosystem service providers can help strengthen arguments for their protection. It will also help us understand the wider implications of large-scale flying fox extinctions, as these giant bats are known to interact with plants on a large landscape scale, performing ecological roles over vast transboundary areas (Epstein et al., 2009). Flying foxes are likely to be particularly important players in island ecosystems where they often serve as keystone pollinators and seed dispersers both within and between islands (Cox et al., 1991; Banack, 1998; McConkey & Drake, 2007; McConkey & Drake, 2015), and where maintaining their numbers at high densities is necessary for the survival of plant communities (McConkey & Drake, 2006). Such data are also important to understand the drivers and potential mitigation strategies for conflicts between fruit bats and humans (Aziz et al., 2016; Aziz et al., in press).

Whilst in-depth, comprehensive dietary/foraging studies have been conducted for certain flying fox species, particularly in Australia (e.g., Boulter et al., 2005; Williams et al., 2006), Oceania (e.g., McConkey & Drake, 2006; Luskin, 2010), Japan (e.g., Nakamoto, Kinjo & Izawa, 2007; Nakamoto, Kinjo & Izawa, 2009; Lee et al., 2009), South Asia (e.g., Mahmood-ul-Hassan et al., 2010; Sudhakaran & Doss, 2012), and Indian Ocean islands (e.g., Nyhagen et al., 2005; Oleksy, Racey & Jones, 2015), the diets of Southeast Asian species, which are some of the most threatened due to the additional threat of commercial hunting (Mildenstein, Tanshi & Racey, 2016), remain largely unknown. Indeed, apart from a few studies in the Philippines (Reiter & Curio, 2001; Mildenstein et al., 2005; Stier & Mildenstein, 2005), Thailand (Weber et al., 2015), and Myanmar (Win & Mya, 2015), all other dietary and foraging studies on Southeast Asian Pteropodidae have focused on the smaller pteropodids (e.g., Hodgkison et al., 2004; Fletcher, Zubaid & Kunz, 2012; Bumrungsri et al., 2013; Stewart, Makowsky & Dudash, 2014). This is of particular concern given that out of the 67 flying fox species listed on the IUCN Red List, almost half (30 species i.e., 45%) are found in Southeast Asia (IUCN, 2016).

Due to the nocturnal and volant nature of bats, invasive analyses (by capturing individuals) or indirect methods (by collecting droppings) have traditionally been used to study their diets. Microscope analyses of pteropodid faeces have provided insights into their diet (Bumrungsri, Leelapaibul & Racey, 2007), as well as their roles in pollination (Bumrungsri et al., 2013) and seed dispersal (Sritongchuay et al., 2014). However, all these studies have relied on physical identification of food plant species—either through direct observations of foraging bats, or microhistological identification of seeds, pollen, fruit fibres and leaf fragments in faeces and ejecta (chewed-up pellets of plant parts spat out by bats after swallowing the juice or soft pulp; Nyhagen et al., 2005; Long & Racey, 2007). The successful use of these approaches relies on several important factors such as accessibility and visibility of foraging bats for the former method, and also the availability of expert botanical knowledge or resources such as reference collections. Another limitation of these conventional approaches is that they require physically identifiable remains to be expelled by the bats; any plant parts that were consumed or expelled solely in liquid form will be missed out in the analysis (Pompanon et al., 2012). Studies on the foraging ecology of wide-ranging species such as flying foxes also require the use of expensive, hi-tech equipment such as GPS collars, which is often not feasible.

Although molecular analysis of pteropodid diets can potentially be used to overcome the obstacles outlined above, this approach has yet to be applied. DNA analyses of faeces collected non-invasively have already been conducted to determine the herbivorous diets of animals such as primates (Bradley et al., 2007; Quéméré et al., 2013; Srivathsan et al., 2016), marmots, bears, capercaillies, grasshoppers, molluscs, slugs (Valentini et al., 2009a), pigeons (Ando et al., 2013) and tapirs (Hibert et al., 2013), but this has never before been attempted for pteropodids in the Palaeotropics. To date, molecular analyses of bat diets in the Old World have only been conducted for insectivorous species (e.g., Clare et al., 2009; Razgour et al., 2011; Zeale et al., 2011). In fact, to our knowledge, the only successful attempt to identify the diet of plant-visiting bats through molecular analysis has been done by one study in the Neotropics (Hayward, 2013).

We evaluated the utility of Illumina Next Generation Sequencing (NGS) to identify plant species present in the droppings of the island flying fox (Pteropus hypomelanus) from Tioman Island in Peninsular Malaysia. In addition, we evaluated the potential of NGS analysis in complementing or even replacing conventional microhistological analysis to elucidate flying fox diets. First, we created a site-specific reference collection of potential flying fox food plants—DNA sequencing of these plants facilitated the construction of a phylogenetic tree to support morphology-based identification of these food plant species and primer design for NGS. Next, we assessed the feasibility of extracting plant DNA from flying fox droppings, and evaluated whether DNA sequences obtained from NGS could be matched with those from online and site-specific DNA reference databases. Finally, we compared the performance of NGS with conventional microhistological analyses to identify food plant species from flying fox droppings.

Materials and Methods

Study site

We conducted this study on Tioman Island (2°48′38″N, 104°10′38″E; 136 km2; Fig. 1A), located 32 km off the east coast of Peninsular Malaysia in the State of Pahang. This research was approved by the Economic Planning Unit of Malaysia (Permit number: 3242). Much of the island inland is still covered by primary tropical rainforest, which has been designated as Pulau Tioman Wildlife Reserve (83 km2). It has a hilly topography, with flat areas only along the coast (Abdul, 1999). The area designated as a wildlife reserve is composed of lowland mixed dipterocarp forest and hill dipterocarp forest. Most forested areas are still inaccessible due to the rugged topography, with many steep slopes and rocky outcrops (Latiff et al., 1999). The climate is tropical, uniformly warm and humid throughout the year (Hasan Basyri et al., 2001), but the island experiences the northeast monsoon from November to March (Bullock & Medway, 1966).

Figure 1
Map of study area and images of sampling site and method.

There are currently seven villages on the island, situated along the coastline (Fig. 1A). The majority of the local people are Muslim, and therefore due to religious dietary restrictions do not hunt bats for food or medicine (Aziz, 2016). As the island’s marine area is also a designated Marine Park and a popular tourist destination, many of the local people are heavily involved in the tourism industry (Abdul, 1999).

Study species

The island flying fox (Pteropus hypomelanus), also known as the variable flying fox and the small flying fox, roosts gregariously, forming colonies of up to 5,000 individuals. It is a widespread insular species, considered to be abundant throughout a distribution range that extends from the Maldives and Indian islands in the west to Melanesia in the east. Because of this, it is considered to be Least Concern on a global scale by the IUCN Red List; however, populations are decreasing (Francis et al., 2008; Olival, 2008), and the species is now listed as endangered on the Malaysian Red list (DWNP , 2010).

On Tioman, the island flying fox can be found roosting permanently in two villages: Tekek, on the west coast, and Juara, on the east coast (Fig. 1A), and forages throughout the island (Medway, 1966; Ong, 2000). Monthly roost counts conducted during March–October 2015 using a thermalscope (Pulsar Quantum HD38S) yielded estimated ranges of 2,178–5,385 individuals for the entire island (see Fig. S1 for further details). Local people have reported that the flying foxes do forage in other villages on the island.

Flying fox dropping collection

Collection of flying fox droppings took place once a month during March-October 2016 (i.e., eight months). Samples of droppings consisting of faeces and ejecta were collected for three mornings in the last week of each month from three separate day roosts in Juara (east coast) and two separate day roosts in Tekek (west coast). The number of roosts and sampling days were determined based on species accumulation curves of pollen morphospecies that were detected through preliminary microhistological analysis in June 2014. Program EstimateS (version 9.1.0; http://viceroy.eeb.uconn.edu/estimates/) indicated that sampling completeness (i.e., observed/estimated number of species; Soberón, Llorente & Onate, 2000) was around 97% using this sampling regime.

In Juara, three suitable roost trees (Fig. 1B) for sampling were selected based on accessibility and also on the highest/largest amount of faecal/ejecta splatter produced under the roost, in order to maximise sample yield. As flying foxes often shifted roosts or temporarily abandoned degraded roosts, this meant that sometimes different roosts were sampled in each location every month or even every morning, although most roosts were consistently sampled each month due to their constant high occupancy and best accessibility.

Plastic sheets measuring 0.8 × 1.0 m were placed under each roost after dark, once the bats had exited the roost to forage. The roosts were then visited the next morning for collection starting at 0700 h and ending at 1200 h (bats typically returned to the roosts around 0500–0600h); the plastic sheets were carefully moved away from the roost to a clear area for processing (Fig. 1C). As it was often difficult to differentiate faeces from ejecta, both were collected and analysed equally as ‘droppings’ (Fig. 1D). Droppings collected for processing were selected based on unique colour and texture, as this was assumed to be representative of plant diversity in the bats’ diet. Following the approach used by Stier & Mildenstein (2005) based on short gut-passage time for flying foxes (12–34 min; Tedman & Hall, 1985), we assumed that each bat voided its last meal once, and therefore each dropping represented a different individual’s food choice. We devised our own novel collection protocol where droppings were collected by swabbing them with a cotton bud, then placing each individual dropping into a 5 ml Eppendorf tube containing ~1,000 μl of 95% ethanol. These tubes were then kept cool in the field, either by storing in a conventional freezer or by using a portable cooler box with ice packs, for 1–3 days before being transported off the island and then stored in a −80 °C freezer.

In order to simultaneously test the utility of NGS and compare it with conventional approaches, we collected two duplicate sets of 10 individual droppings from one single roost in Juara village during a single morning on 6 May 2015. One sample set was then kept in a conventional fridge for microscope analysis, whilst the other set was stored in the −80 °C freezer for molecular analysis.

Site-specific plant reference collection

We first checked a published list of genera of known food plants for Pteropus across its range (Marshall, 1985), cross-checked this against a preliminary checklist of seed plants for Tioman (Latiff et al., 1999), and also obtained information on possible flying fox food plants through talking to local people in Juara. We then searched for genera of similar plants in and around the two villages with the aid of a local plant expert. The botanical identification of plants (at least to genus) were subsequently verified by a trained botanist familiar with plants from the region. When we found an individual plant from one of these genera, we recorded its GPS location and collected pollen, fruit, and/or seeds if it was flowering or fruiting. We also collected 3–5 mature leaves for DNA extraction. The leaves were stored in Ziploc bags with silica gel under cool conditions to retard decomposition rates. Leaf samples from 19 different plant species were obtained for this purpose, constituting a preliminary library (Table S1).

Genomic DNA was extracted from approximately 25 mg of one leaf from each plant species using DNAeasy Plant Mini Kit (Qiagen, Halden, Germany) according to the manufacturer’s protocols. Primers targeted the rbcL gene, a protein-coding gene associated with the chloroplast genome of all living plants. DNA amplifications were performed in a mastermix containing 1 µL of DNA, 25 µL of OneTaq Quick-Load 2X Master Mix with Standard Buffer, (New England Biolab, Ipswich, MA, USA), 1 µL of 10mM forward primer rbcLaf-M13, 1 µL of 10 mM reverse primer rbcLa-revM13 (Table S2), and 22 µL of nuclease-free water. The PCR protocol was started with an initial denaturation step for 30 s at 94 °C, followed by 30 cycles of 30 s at 94 °C, 30 s at 48 °C, 40 s at 68 °C, and final elongation for 2 min at 68 °C. The PCR products were purified using 0.8× volume ratio of Agencourt Ampure XP beads (Beckman Coulter, Inc). The purified samples were sent to 1st BASE laboratories (http://www.base-asia.com/) for Sanger sequencing. The sequencing results were quality trimmed using CodonCode TraceViewer (http://www.codoncode.com/TraceViewer/) and aligned using MAFFT version 7.0 (Katoh & Standley, 2013).

A phylogenetic tree (Fig. S4) was constructed to support morphology-based identification of site-specific reference plants. Additional rbcL gene sequences that displayed high similarity to the sequenced plant specimens were mined from the BOLD/NCBI databases, aligned with MAFFT version 7.0 and trimmed using TrimAl version 1.9 optimized for maximum likelihood tree construction (Capella-Gutiérrez, Silla-Martínez & Gabaldón, 2009). The trimmed alignment was subsequently used to construct a maximum likelihood tree with FastTree using the –nt (nucleotide) and –gtr (generalised time-reversible model) setting.

Laboratory procedure

Based on the alignment of rbcL gene fragments from the site-specific reference plant DNA (Fig. S2), we designed a primer pair targeting 198 bp of the rbcL gene (minus primer sequence) using Primer3 (http://bioinfo.ut.ee/primer3-0.4.0/) on default settings (Table S2). Prior to primer synthesis, partial Illumina adapter sequences were added to the 5′ end of the designed primers, rbcL-357F and rbcL-556R, to allow barcoding and sequencing on the Illumina platform.

Individual droppings were pooled according to roost (n = 5, 2 in Tekek and 3 in Juara) and month (n = 8), creating 40 separate mixtures for analysis. The tubes containing the daily samples were first vortexed for 2 min to homogenise the content and subsequently, 1,000 µL of the sample was pipetted into another tube to form the mixture. Next, 100 µL of the mixture was used for DNA extraction similarly using DNAeasy Plant Mini Kit (Qiagen, Halden, Germany) instead of a stool-specific DNA extraction kit to improve the recovery of plant-derived DNA from faecal samples.

PCR reaction was performed using IlluM_rbcLF and IlluM_rbcLR. The 20 µL PCR cocktail consists of 10 µL Q5 Hot Start High-Fidelity 2× Master Mix (New England Biolab, Ipswich, MA, USA), 1 µL each of 10 µM forward and reverse primer, 1 µL gDNA and 7 µL nuclease-free water. All reactions were performed in a Veriti® 96-Well Fast Thermal Cycler with the following protocol: initial denaturation for 30 s at 98 °C, 25 cycles of 10 s at 98 °C, 30 s at 55 °C and 10 s at 65 °C, with a final 1 min extension at 65 °C. The PCR product was purified using 0.8× vol. ratio Agencourt Ampure XP beads (Beckman Coulter, Inc). Then, 1 uL of Index 1 and Index 2 primers from Nextera XT kit were added to 3 uL of purified PCR product and combined with 5 uL of Q5 Hot Start High-Fidelity 2X M aster Mix (New England Biolabs, Ipswich, MA). The PCR protocol was as followed: initial denaturation for 30 s at 98 °C, 8 cycles of 10 s at 98 °C and 1 min at 65 °C, with a final 1 min extension at 65 °C.

Sequence analysis and taxon assignation

The purified amplicons containing the full length Illumina adapter and appropriate unique barcode were then quantified using KAPA Library Quantification kit (Kapa Biosystems, CapeTown, South Africa) on the EcoRealTime PCR system (Illumina, San Diego, CA, USA). Based on the qPCR data, the amplicons were normalised, pooled and subsequently sequenced on the MiSeq (2 × 250 bp paired-end run) located at the Monash University Malaysia Genomics Facility.

Illumina Nextera adapters and primer sequences of the reads were trimmed off using Trimmomatic v 0.33 and FastX trimmer, respectively (Bolger, Lohse & Usadel, 2014; http://hannonlab.cshl.edu/fastx_toolkit/). The trimmed paired-end reads were then merged using PEAR (Zhang et al., 2014) using default settings. Dereplication, singleton removal and Operationally Taxonomic Unit (OTU) clustering (default setting of 97% identity clustering; -id 0.97) were performed using the pipeline implemented in UPARSE (Edgar, 2013). The filtered OTUs were translated into protein sequence and manually inspected for translated sequence containing stop codon(s), which were removed from subsequent analyses. Taxonomic assignment of OTUs was conducted by matching sequences against those from the BOLD/NCBI non-redundant nucleotide database as of December 2016. The BLASTn output for each OTU was downloaded in “XML” format and imported into MEGAN6 (Huson et al., 2016) to calculate their Lowest Common Ancestor (LCA). LCA predicts the taxonomic rank of OTUs based on matches to a set of reference lineages such that the predicted taxonomic rank is shared by all matches within user pre-defined parameters. The implemented parameters for MEGAN6 were: Minimum Blast bit score = 200, Max expected E-value = 0.01, Min Percent Identity = 90.0, Top percent (the percentage of the best score a match needs to lie within) = 10 and Weighted LCA% = 75. To estimate the OTU relative abundance in each sample, reads were mapped to the filtered OTUs via USEARCH (97% similarity cut-off) and normalised to 10,000 reads (Caporaso et al., 2010).

Microhistological analysis

For the 10 dropping samples collected in May, we sent one set for NGS analysis (following the protocol above) and used another set for microscope analysis. For the latter, we first manually broke up the dropping contents in the tube to produce a relatively more representative liquid sample. We then droppped 1–3 drops of this liquid onto a microscope slide using a pipette. Fuchsin jelly was added to this in order to stain pollen grains within the dropping, a slip cover was placed on top, and the jelly was then melted over an open flame, sealing the slip cover to the slide. The slide was then cooled down in a conventional fridge in order to allow the jelly to solidify again before examination.

Once cooled, we placed the slide under a conventional light microscope (Leica DM E) and first examined it using 10/0.25 magnification in order to detect pollen grains and other plant parts. We used a self-made reference collection as well as photos from Start (1974), S. Bumrungsri (http://www.seabcru.org/seabcru-resources) and Mohamed (2014) to identify pollen and other plant parts. When necessary, we used higher magnification (40/0.65) to view the pollen grains. We used just ‘presence/absence’ to assess pollen to avoid quantification biases towards species that naturally produce greater amounts of pollen. Following Thomas (2009), we considered a species present if we found three or more pollen grains on a single slide.

Results

Site-specific plant reference database

For the 19 specimens in the site-specific plant reference database, the phylogenetic tree (Fig. S4) based on a 689 bp aligned rbcL gene region recovered the monophyly of several plant genera with high SH-like local support (>0.90). However, the sister grouping among genera in general was not well-supported, possibly due to the lack of comprehensive plant taxon sampling, and the use of a single genetic marker. The rbcL-based phylogeny generally supported the morphology-based identification of reference specimens of potential flying fox food plants, with 15 specimens correctly identified at least to genus level as shown by their monophyletic clustering with plant species of the same genus (Table 1). Sample PTMN02 (identified as Streblus asper based on morphology) was an interesting exception as it formed a sister group with the Ficus clade (SH-like local support = 0.938), but not with other congeners. rbcL sequences representing the plant species collected on Tioman were scarce in the BOLD/NCBI databases, suggesting these databases have insufficient sequence representation that can affect taxonomic assignment of plant DNA in the droppings.

Table 1
Summary information of 29 OTUs detected in flying fox droppings over eight months (Mar–Oct 2015) in Tioman Island, Malaysia: OTU identities (genus in bold) based on matches with plant sequences from online reference database using the Least Common ...

Feasibility of using NGS to study flying fox diet

With our newly designed rbcL primer, we were able to successfully extract, amplify, and subsequently identify plant DNA from all of the collected flying fox droppings. Initially, a total of 160 Operationally Taxonomic Units (OTUs) were recovered from the sequencing reads, of which 29 OTUs (Table 1) were retained after filtering based on cumulative relative abundance (>0.5%) and presence of stop codon(s) in reading frame. Using a conservative LCA approach, we identified at least three different plant genera and at least 14 plant families from the droppings (Table 1). In addition, 8 OTUs matched with specimens from the site-specific plant reference collection.

Based on sampling completeness (calculated using EstimateS) for OTU relative abundance data from five roosts (data pooled over three days) per month using Chao 1 species richness estimator (good for datasets skewed towards low abundance classes; Chao, 1984), sampling completeness was relatively high for the months March (99%), April (100%), June (100%), August (100%), September (88%), and October (96%). However, sampling completeness could be improved for May (55%) and July (79%). The month of May, which had the lowest overall sampling completeness, also had the highest number of droppings collected (Table S3).

The results from our NGS analysis of island flying fox droppings over eight months suggest that the diet at both Juara and Tekek during this time was dominated by four different plant taxa (Table 1; Fig. 2): OTU 1 (Ficus sp.), OTU 3 (likely to be Mangifera indica based on sequence match with site-specific plant reference collection) and OTUs 4 & 5 (Rubiaceae). Spatio-temporal patterns in the relative abundance of these four taxa in the diet were observed during the sampling period (Fig. 2). For example, OTU 5 appeared to be consumed in similar proportions at both Juara and Tekek across all months whereas OTU 4 was consistently consumed in low proportions in Tekek yet consumed irregularly in Juara over the same period (Fig. 3). Even between different roosts in the same site, spatio-temporal differences were observed, such as for OTU 7 (Fig. S5).

Figure 2
Relative abundance of 8 OTUs detected in flying fox droppings.
Figure 3
Spatio-temporal trends in predicted consumption of the top four most dominant plant taxa.

Performance of NGS vs. microhistological analysis

Microscope analysis identified two plant taxa in flying fox droppings (Table 2). Out of 10 individual droppings, three contained durian (Durio sp.) pollen. All the other droppings contained fig parts exclusively; no other plant parts were detected. Durian pollen occurred at extremely low abundance; in all cases, only 3–4 grains were detected per slide. No other pollen or plant parts were detected. NGS identified the same two plant taxa detected by microhistological analysis, and further identified an additional six plant taxa. However, Durio was not detected in the same samples as those identified via microscope.

Table 2
Comparative effectiveness of microhistological vs. NGS analyses in identifying plants from two duplicate sets of 10 samples of flying fox droppings collected on 6 May 2015.

Discussion

Our study is the first to describe the diet of the island flying fox, which was previously unknown. To our knowledge, this is also the first use of NGS to identify plant taxa in the diet of a pteropodid, which has been difficult to characterise due to these animals’ volant nature, large home ranges and nocturnal foraging behaviour. Figs consistently formed the highest amount of plant taxa detected in the droppings each month, at both sampling sites. This strongly suggests that figs compose the core diet of flying foxes on the island. It is thus highly likely that the island flying fox plays a key role in dispersing fig seeds throughout Tioman, making these bats important keystone species for the island (Cox et al., 1991; McConkey & Drake, 2015); future studies on seed dispersal and germination are required to confirm this.

NGS is a reliable tool to study flying fox diet

We have demonstrated that identification of plant taxa to family level is generally possible based on the partial sequence of rbcL using the LCA approach. In addition, some OTUs in our study were successfully assigned to the genus level. In order to be conservative, however, we avoided assigning most OTUs to species level, unless there were matches with BOLD/NCBI database sequences and site-specific reference plant sequences. Asspecies-level plant identification based solely on NGS is not straightforward, the use of a site-specific food plant DNA reference database is vital to help identify plant species in flying fox diets.

Other genes have been successfully used to identify plant species in animal diets. For example, Valentini et al. (2009a) found the trnL intronic region to be effective for Asian mammals, birds, and invertebrates, identifying 50% of the plant taxa found in the diets of these animals to species level. The same approach has been used for European bison (Kowalczyk et al., 2011), alpine chamois (Raye’ et al., 2011), and red-headed wood pigeons (Ando et al., 2013). The trnL intron has been previously reported to evolve more than three times faster than the protein-coding rbcL gene, thus potentially harbouring more variation and phylogenetic signal per base pair (Gielly & Taberlet, 1994). However, we chose the rbcL gene instead of the P6 loop of the trnL intron to study flying fox diet, as the rbcL gene is currently one of the two genes (the other being matK) that is increasingly being used for plant species identification. This is useful for studying plants originating from a less-studied region of high biodiversity such as Southeast Asia. Unlike the trnL database, the rbcL database is consistently growing in the BOLD database, and as a result is also likely to be represented by more sequences from properly identified and vouchered specimens. In addition, given that rbcL is a protein-coding gene, it enables the screening of erroneous OTU resulting from sequencing and/or amplification errors based on the presence of stop codon(s) in the translated reading frame. Instead of using recently developed and robustly tested primers (e.g., Little, 2014), we designed a new set of rbcL minibarcode primers due to the current lack of rbcL sequence representation for plant species from Malaysia or specifically, Tioman. One caveat of this approach would be the possibility of preferential primers binding to known diet items instead of unexpected ones. However, our results seem to indicate that our newly designed primers were, to a certain extent, capable of recovering OTUs belonging to a wide variety of plant families, in part due to the diverse representation of plant reference specimens that contributed to the expanded taxon coverage of our primers. Future studies aiming to achieve greater power of identification of plants in flying fox diets could consider using more than one target region (e.g., Hibert et al., 2013; Clare, 2014) coupled with an additional in silico PCR optimisation step using ecoprimers to improve the reliability and universality of the newly designed primers (Riaz et al., 2011). Also, in order to completely eliminate bias associated with PCR, metagenomic shotgun sequencing could be performed, albeit at a relatively higher cost depending on the required sequencing coverage. (e.g., Srivathsan et al., 2015; Srivathsan et al., 2016).

NGS can complement microhistological analysis

Our results showed that NGS can provide greater insights into the diet of flying foxes than conventional microhistological approaches by detecting a wider range of plant taxa, thus highlighting the utility and discriminatory potential of the newly designed rbcL primers to study flying fox diets. More importantly, the use of NGS allowed us to identify plant species even when no physical plant parts were found in the flying fox droppings. The plant genera and families detected from NGS have also been recorded by botanists as being present on Tioman, including the top four genera/families detected most abundantly in the droppings (Latiff et al., 1999; Mohd. Norfaizal et al., 2014).

In our study, attempts to use microscope analysis to identify plant parts in droppings proved to be challenging, as no pre-existing reference collection was available. Building our own microhistological reference collection for Tioman was time-consuming and labour-intensive—and the resulting collection often did not match up with the plant parts found in the flying fox droppings. However, obtaining DNA from plant specimens is still necessary to narrow down the identity of OTUs to species level. Indeed, 8 out of 29 OTUs had 100% matches to the sequences of plant specimens collected from the study site, highlighting the importance of building a comprehensive local sequence library beforehand, preferably specific to one’s particular study site.

It is important to note that NGS did not detect Durio in the same individual droppings as those identified via microscope. This is likely due to the low abundance of this plant taxon in the droppings affecting detection probability, especially since the NGS analysis used a more general primer that was not specific to Durio. This pollen detection probability is another caveat to be aware of; Scanlon et al. (2014) have cautioned that faecal subsampling methods can potentially lead to inaccurate detection of pollen in dietary studies, regardless of which method is used.

Caveats

Our sample collection method in the field, selecting only for droppings with unique colour and texture, may have introduced a bias that could result in underestimating the relative abundance of OTUs in the droppings. In particular, sampling completeness for the months May, June and July were relatively low, showing that more roosts and/or days needed to be sampled in order to obtain a complete representation of diet for these months. However, roost count data (Fig. S1) show that this high diversity in diet was not influenced by population abundance. Instead, perhaps diet choice and/or food resource diversity were relatively higher during this period compared to the rest of the year. Future studies should aim to collect every single dropping found underneath a roost to improve representativeness.

Given the potentially short flying fox gut passage times (Tedman & Hall, 1985), droppings collected from day roosts in the morning may also bias the analysis results towards food items that were consumed only at the end of the foraging period (Schmelitschek, French & Parry-Jones, 2009). Food plants that were consumed during the start or middle of the evening may not have been detected by our methods. Although Banack & Grant (2002) observed flying foxes returning to food resources that were foraged upon earlier, before then returning to day roosts, this is still a potential caveat to bear in mind. For example, primates are known to exhibit temporal patterning in diet choice, structuring their diet throughout their foraging period with different food items; it is believed that this is due to how different foods are processed, and give energy, at different rates, and therefore helps to ensure that the animals maintain high energy levels (Robinson, 1984; Ganzhorn & Wright, 1994; Chapman & Chapman, 1991). Given the size of Tioman, and the logistical challenges of observing flying foxes foraging, the best way to overcome this possible information gap is to conduct GPS tracking studies.

We also acknowledge that NGS approaches to diet identification are semi-quantitative because chloroplast abundance is variable in different plant species and different parts of the leaf. Ultimately, the ability of NGS to accurately identify food plants will always depend on sequence specificity of the primers. While the NGS approach has proven to be useful in elucidating the island flying fox’s varied diet on Tioman, for animals with such a diverse phytophagous diet, primer specificity will always be a limiting factor and there is a chance that unknown plant species will not be detected due to primer mispriming. Also, identical chloroplast DNA sequences can be present in different but related species, making it impossible to distinguish closely related plant species from each other in the diet. This could be one possible factor why several OTUs could only be identified to order/family/subfamily levels, suggesting that they require further phylogenetic investigation and/or may benefit from identification based on more rapidly-evolving plastid-coding genes such as matK. It is worth noting that the sequenced rbcL gene of some plant specimens collected in this study did not exhibit 100% identity matches to species in the BOLD/NCBI databases, which may be attributed to genetic diversity at the intra-species level or gaps in the database i.e., certain plant species consumed by the flying foxes may not yet have their corresponding sequences deposited in the database. Nevertheless, there is an urgent need for the BOLD/NCBI databases to have more representation of plant sequences from Peninsular Malaysia and Southeast Asia in general.

Another limitation of the NGS approach for generalist diets is that it does not identify which part of the plant was consumed. For specialised frugivores, nectarivores, or herbivores (e.g., folivorous), this may not be an issue. Flying foxes, however, are generalists which consume fruits, flowers, nectar, and even leaves (Marshall, 1985). It is this dietary plasticity which allows them to perform more than one ecological role in tropical landscapes. Therefore, identifying which plant parts are actually consumed is a crucial step towards identifying the ecosystem services that these bats provide. NGS can provide a first step towards identifying flying fox diet but should not be viewed as a replacement for microhistological analysis. Nevertheless, this approach has shed new light on flying fox diet by discovering plant taxa that were entirely missed out by the conventional approach. Ideally, studies using NGS should be combined with micro-histological analysis in order to fill in the gaps and broaden our understanding of pteropodid diet and foraging ecology. NGS can also be used in combination with comprehensive and long-term data on plant phenology, to observe which food resources are available at which time. Following on from this preliminary study, the identification of specific food plants via NGS can now help guide more in-depth plant sample collection and phenological observations.

Conclusion

Our study is the first to use NGS to identify potential plant species in flying fox diet, paving the way for a new approach to studying flying fox diets. Since our NGS analysis of flying fox diet was semi-quantitative, it is not yet possible to make any definite conclusions regarding food preference vs. food availability; ultimately it is unclear to what extent sampling bias and detection probability may have influenced the type and relative abundance of plant taxa detected in our study. Yet some of the interesting patterns we observed are worth investigating in greater detail, particularly in combination with microhistological analysis. The results will also help to guide us in conducting more accurate and expanded phenology monitoring, and further collection of botanical samples. Further and more rigorous sampling, especially at the level of the individual animal, is required to understand the dietary patterns of this particular flying fox population, expand on the information provided here and build on our understanding of how flying foxes provide ecosystem services on Tioman Island and elsewhere.

Supplemental Information

10.7717/peerj.3176/supp-1

Figure S1

Monthly roost counts of flying foxes in Tekek and Juara villages:

Monthly roost counts of flying foxes in Tekek and Juara villages between Mar and Oct 2015, Tioman Island. Due to the accessibility of the roosts, these counts were found to be more accurate than exit/emergence/fly-out counts. Exit counts were found to be impractical, as the bats did not fly out in steady streams towards specific directions, and no appropriate vantage points existed to be used as counting stations.

10.7717/peerj.3176/supp-2

Figure S2

Sequence alignment of partial rbcL genes for plants collected from Tioman Island:

Residues that are divergent across all sequences are highlighted in green (>90% similarity) and yellow (<90% similarity).

10.7717/peerj.3176/supp-3

Figure S3

Overview of our newly designed primers for this study and expected construct consisting of the complete Illumina adapter, dual index barcode and partial rbcL gene.
10.7717/peerj.3176/supp-4

Figure S4

Maximum likelihood phylogenetic tree:

Maximum likelihood phylogenetic tree depicting the evolutionary relationship among rbcL sequences obtained from individually collected leaf samples (red-coloured tips) of potential flying fox food plants and rbcL sequences from online reference databases. PTMN codes refer to barcodes generated in this study and accession numbers are shown next to each taxon. Values in nodes indicate SH-like local support values. Scale bar indicates number of substitution per site.

10.7717/peerj.3176/supp-5

Figure S5

Spatio-temporal trends in consumption of OTU 7 (Moraceae) by flying foxes, suggesting possible inter-roost variation at Juara and Tekek villages between March and October 2016 on Tioman Island:
10.7717/peerj.3176/supp-6

Table S1

Summary information of 19 potential flying fox food plant reference specimens:

Summary information of 19 potential flying fox food plant reference specimens obtained from Tioman Island, Peninsular Malaysia: identification of plant specimen based on morphology, genus match(tick)/non-match(cross) based on rbcL sequences (Fig. 3), approximate GPS coordinates where the plant specimens were collected, their specimen code and GenBank accession codes.

10.7717/peerj.3176/supp-7

Table S2

Primers used in this study for the amplification of rbcL gene fragment from flying fox droppings:
10.7717/peerj.3176/supp-8

Table S3

Number of droppings per roost sampled from Juara (1406 droppings) and Tekek (1,336 droppings) over 8 months on Tioman Island:
10.7717/peerj.3176/supp-9

Data S1

Fasta files:

Acknowledgments

We thank the people of Juara for their hospitality and support. We are grateful to the Economic Planning Unit of Malaysia (Permit No. 3242) and the Department of Wildlife and National Parks (DWNP) for allowing this research to be conducted. In particular, we thank Dr. Jeffrine Rovie Ryan Japning and Frankie Thomas Sitam for the use of DWNP lab facilities. Special thanks to Lam Wai Yee, Esteban Brenes-Mora, Anna Deasey, Noraisah Majri, Mahfuzatul Izyan, Khatijah Haji Hussin, Joanne Tong, Yek Sze Huei, Liz Moleski, Sri Rao Venkateswara, Foon Junn Kitt, Jackie Tan May Li, David Bickford, Mary Rose Posa, Lim Lee Sim, Jasdev Sohanpal and Manpreet Kaur for assisting with the collection of droppings. We are indebted to Sara Bumrungsri for training in dropping collection and microscope analysis, and Kartini Mohamed and Muhamad Hamirul for facilitating the use of light microscopes. We are also grateful to two anonymous reviewers and Johan Ponsu for providing constructive feedback to improve the manuscript.

Funding Statement

This work was supported by The Rufford Foundation (grant number 17325-1), Bat Conservation International (Scholarships and Grassroots Grants for Bat-Centric Projects Focused on Critical Conservation Needs), Muséum National d’Histoire Naturelle, University of Nottingham Malaysia Campus, Dr. Yeoh Suat Hui from Universiti Malaya, and Ms. Jade Ong. SA Aziz was the recipient of a PhD scholarship from the French Government issued by the French Embassy in Malaysia. GR Clements was supported by a LABEX BCDiv grant while in France. There was no additional external funding received for this study. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.

The following grant information was disclosed by the authors:

The Rufford Foundation 17325-1.
Bat Conservation International .
Muséum National d’Histoire Naturelle, University of Nottingham Malaysia Campus .
Universiti Malaya .
French Embassy in Malaysia .
LABEX BCDiv .

Additional Information and Declarations

Competing Interests

Sheema Abdul Aziz and Gopalasamy Reuben Clements are co-founders of Rimba, Selangor, Malaysia.

Author Contributions

Sheema Abdul Aziz and Han Ming Gan conceived and designed the experiments, performed the experiments, analyzed the data, contributed reagents/materials/analysis tools, wrote the paper, prepared figures and/or tables, reviewed drafts of the paper.

Gopalasamy Reuben Clements conceived and designed the experiments, performed the experiments, wrote the paper, prepared figures and/or tables, reviewed drafts of the paper.

Lee Yin Peng performed the experiments, analyzed the data, wrote the paper, prepared figures and/or tables.

Ahimsa Campos-Arceiz, Kim R. McConkey and Pierre-Michel Forget reviewed drafts of the paper.

Field Study Permissions

The following information was supplied relating to field study approvals (i.e., approving body and any reference numbers):

This research was approved by the Economic Planning Unit of Malaysia (Permit number: 3242).

DNA Deposition

The following information was supplied regarding the deposition of DNA sequences:

Sanger sequencing results for the rbcL fragment of collected plant specimens were assigned accession numbers in GenBank/NCBI (see Table S1). NGS data were also deposited in the SRA database as mentioned in the Data Availability Section.

Accession numbers:

SRX1988609, SRX1988608, SRX1988607, SRX1988606, SRX1988605, SRX1988604, SRX1988603, SRX1988602, SRX1988601, SRX1988600, SRX1988599, SRX1988598, SRX1988597, SRX1988596, SRX1988595, SRX1988594, SRX1988593, SRX1988592, SRX1988591, SRX1988590, SRX1988589, SRX1988588, SRX1988587, SRX1988586, SRX1988585, SRX1988584, SRX1988583, SRX1988582, SRX1988581, SRX1988580, SRX1988579, SRX1988578, SRX1988577, SRX1988576, SRX1988575, SRX1988574, SRX1988573, SRX1988572, https://www.ncbi.nlm.nih.gov/sra?term=SRX1988571, SRX1988570, SRX1988569, SRX1988568, SRX1988567, SRX1988566, SRX1988565, SRX1988564, SRX1988563, SRX1988562, SRX1988561, SRX1988560

https://www.ncbi.nlm.nih.gov/sra/SRP080299.

Data Availability

The following information was supplied regarding data availability:

The research in this article only generated DNA sequences, which have already been listed in detail in the relevant section above. No other raw data or code were generated.

References

Abdul (1999) Abdul J. An introduction to Pulau Tioman. Raffles Bulletin of Zoology. 1999;6:3–4.
Ando et al. (2013) Ando H, Setsuko S, Horikoshi K, Suzuki H, Umehara S, Inoue-Murayama M, Isagi Y. Diet analysis by next-generation sequencing indicates the frequent consumption of introduced plants by the critically endangered red-headed wood pigeon Columba janthina nitens in oceanic island habitats. Ecology and Evolution. 2013;3:4057–4069. doi: 10.1002/ece3.773. [PMC free article] [PubMed] [Cross Ref]
Aziz (2016) Aziz SA. PhD thesis. 2016. Interactions between flying foxes, plants and people: implications for conservation.
Aziz et al. (in press) Aziz SA, Clements GR, Giam X, Forget P-M, Campos-Arceiz A. Coexistence and conflict between the Island Flying Fox (Pteropus hypomelanus) and local people on a tropical island in Peninsular Malaysia. Human Ecology. 2017 In Press.
Aziz et al. (2016) Aziz SA, Olival KJ, Bumrungsri S, Richards GC, Racey PA. Bats in the anthropocene: conservation of bats in a changing world. Springer International Publishing; 2016. The conflict between pteropodid bats and fruit growers: species, legislation and mitigation; pp. 377–426.
Banack (1998) Banack SA. Diet selection and resource use by flying foxes (genus Pteropus) Ecology. 1998;79:1949–1967.
Banack & Grant (2002) Banack SA, Grant GS. Spatial and temporal movement patterns of the flying fox, Pteropus tonganus, in American Samoa. The Journal of Wildlife Management. 2002;66:1154–1163. doi: 10.2307/3802948. [Cross Ref]
Bolger, Lohse & Usadel (2014) Bolger AM, Lohse M, Usadel B. Trimmomatic: a flexible trimmer for Illumina sequence data. Bioinformatics. 2014;30(15):2114–2120. doi: 10.1093/bioinformatics/btu170. [PMC free article] [PubMed] [Cross Ref]
Boulter et al. (2005) Boulter SL, Kitching RL, Howlett BG, Goodall K. Any which way will do—the pollination biology of a northern Australian rainforest canopy tree (Syzygium sayeri; Myrtaceae) Botanical Journal of the Linnean Society. 2005;149:69–84. doi: 10.1111/j.1095-8339.2005.00430.x. [Cross Ref]
Bradley et al. (2007) Bradley BJ, Stiller M, Doran-Sheehy DM, Harris T, Chapman CA, Vigilant L, Poinar H. Plant DNA sequences from feces: potential means for assessing diets of wild primates. American Journal of Primatology. 2007;69:699–705. doi: 10.1002/ajp.20384. [PubMed] [Cross Ref]
Bullock & Medway (1966) Bullock JA, Medway L. The fauna of Pulau Tioman. I. General information. Bulletin of the Natural Museum. 1966;34:1–8.
Bumrungsri et al. (2013) Bumrungsri S, Lang D, Harrower C, Sripaoraya E, Kitpipit K, Racey PA. The dawn bat, Eonycteris spelaea Dobson (Chiroptera: Pteropodidae) feeds mainly on pollen of economically important food plants in Thailand. Acta Chiropterologica. 2013;15:95–104. doi: 10.3161/150811013X667894. [Cross Ref]
Bumrungsri, Leelapaibul & Racey (2007) Bumrungsri S, Leelapaibul W, Racey PA. Resource partitioning in sympatric Cynopterus bats in lowland tropical rainforest, Thailand. Biotropica. 2007;39:241–248. doi: 10.1111/j.1744-7429.2006.00245.x. [Cross Ref]
Capella-Gutiérrez, Silla-Martínez & Gabaldón (2009) Capella-Gutiérrez S, Silla-Martínez JM, Gabaldón T. trimAl: a tool for automated alignment trimming in large-scale phylogenetic analyses. Bioinformatics. 2009;25:1972–1973. doi: 10.1093/bioinformatics/btp348. [PMC free article] [PubMed] [Cross Ref]
Caporaso et al. (2010) Caporaso JG, Kuczynski J, Stombaugh J, Bittinger K, Bushman FD, Costello EK, Fierer N, Pena AG, Goodrich JK, Gordon JI, Huttley GA. QIIME allows analysis of high-throughput community sequencing data. Nature Methods. 2010;7:335–336. doi: 10.1038/nmeth.f.303. [PMC free article] [PubMed] [Cross Ref]
Chao (1984) Chao A. Non-parametric estimation of the number of classes in a population. Scandinavian Journal of Statistics. 1984;11:265–270.
Chapman & Chapman (1991) Chapman CA, Chapman LJ. The foraging itinerary of spider monkeys: when to eat leaves? Primatologica. 1991;56:162–166.
Clare (2014) Clare EL. Molecular detection of trophic interactions: emerging trends, distinct advantages, significant considerations and conservation applications. Evolutionary Applications. 2014;7:1144–1157. doi: 10.1111/eva.12225. [PMC free article] [PubMed] [Cross Ref]
Clare et al. (2009) Clare EL, Fraser EE, Braid HE, Fenton MB, Hebert PD. Species on the menu of a generalist predator, the eastern red bat (Lasiurus borealis): using a molecular approach to detect arthropod prey. Molecular Ecology. 2009;18:2532–2542. doi: 10.1111/j.1365-294X.2009.04184.x. [PubMed] [Cross Ref]
Cox et al. (1991) Cox PA, Elmqvist T, Pierson ED, Rainey WE. Flying foxes as strong interactors in South Pacific island ecosystems: a conservation hypothesis. Conservation Biology. 1991;5:448–454. doi: 10.1111/j.1523-1739.1991.tb00351.x. [Cross Ref]
DWNP (2010) DWNP (Department of Wildlife and National Parks Peninsular Malaysia) Red List of mammals for Peninsular Malaysia. DWNP; Cheras, Kuala Lumpur: 2010.
Edgar (2013) Edgar RC. UPARSE: highly accurate OUT sequences from microbial amplicon reads. Nature Methods. 2013;10:996–998. doi: 10.1038/nmeth.2604. [PubMed] [Cross Ref]
Epstein et al. (2009) Epstein JH, Olival KJ, Pulliam JRC, Smith C, Westrum J, Hughes T, Dobson AP, Zubaid A, Sohayati AR, Misliah MB, Field HE, Daszak P. Pteropus vampyrus, a hunted migratory species with a multinational home-range and a need for regional management. Journal of Applied Ecology. 2009;46:991–1002. doi: 10.1111/j.1365-2664.2009.01699.x. [Cross Ref]
Fletcher, Zubaid & Kunz (2012) Fletcher C, Zubaid A, Kunz TH. Fruit diet of frugivorous bats (Cynopterus brachyotis and Cynopterus horsfieldii) in tropical hill forests of Peninsular Malaysia. Mammalia. 2012;76:389–397.
Francis et al. (2008) Francis C, Rosell-Ambal G, Bonaccorso FA, Heaney L, Molur S, Srinivasulu C. Pteropus hypomelanus. 2008. [19 April 2016]. The IUCN Red List of Threatened Species 2008: e.T18729A8511932. [Cross Ref]
Fujita & Tuttle (1991) Fujita MS, Tuttle MD. Flying foxes (Chiroptera: Pteropodidae): threatened animals of key ecological and economic importance. Conservation Biology. 1991;5:455–463. doi: 10.1111/j.1523-1739.1991.tb00352.x. [Cross Ref]
Ganzhorn & Wright (1994) Ganzhorn JU, Wright PC. Temporal patterns in primate leaf eating: the possible role of leaf chemistry. Folia Primatologica. 1994;63:203–208. [PubMed]
Gielly & Taberlet (1994) Gielly L, Taberlet P. The use of chloroplast DNA to resolve plant phylogenies: noncoding versus rbcL sequences. Molecular Biology and Evolution. 1994;11:769–777. [PubMed]
Hasan Basyri et al. (2001) Hasan Basyri D, Norhan AR, Amir Hashim MK, Kamarul Azlan MN. Study on water resources in Tioman Island. Jurnal Teknologi. 2001;34(F):51–64.
Hayward (2013) Hayward CE. MSc thesis. 2013. DNA barcoding expands dietary identification and reveals dietary similarity in Jamaican frugivorous bats.
Hibert et al. (2013) Hibert F, Taberlet P, Chave J, Scotti-Saintagne C, Sabatier D, Richard-Hansen C. Unveiling the diet of elusive rainforest herbivores in next generation sequencing era? The tapir as a case study. PLOS ONE. 2013;8(4):e60799. doi: 10.1371/journal.pone.0060799. [PMC free article] [PubMed] [Cross Ref]
Hodgkison et al. (2004) Hodgkison R, Balding ST, Zubaid A, Kunz TH. Temporal variation in the relative abundance of fruit bats (Megachiroptera: Pteropodidae) in relation to the availability of food in a lowland Malaysian rain forest. Biotropica. 2004;36:522–533.
Huson et al. (2016) Huson DH, Beier S, Flade I, Górska A, El-Hadidi M, Mitra S, Ruscheweyh HJ, Tappu R. MEGAN community edition-interactive exploration and analysis of large-scale microbiome sequencing data. PLOS Computational Biology. 2016;12(6):e1004957. doi: 10.1371/journal.pcbi.1004957. [PMC free article] [PubMed] [Cross Ref]
IUCN (2016) IUCN IUCN Red List of Threatened Species. Version 2016-1. http://www.iucnredlist.org. [4 July 2016];2016
Katoh & Standley (2013) Katoh K, Standley DM. MAFFT multiple sequence alignment software version 7: improvements in performance and usability. Molecular Biology and Evolution. 2013;30:772–780. doi: 10.1093/molbev/mst010. [PMC free article] [PubMed] [Cross Ref]
Kingston (2010) Kingston T. Research priorities for bat conservation in Southeast Asia: a consensus approach. Biodiversity Conservation. 2010;19:471–484. doi: 10.1007/s10531-008-9458-5. [Cross Ref]
Kowalczyk et al. (2011) Kowalczyk R, Taberlet P, Coissac E, Valentini A, Miquel C, Kamiński T, Wójcik JM. Influence of management practices on large herbivore diet—Case of European bison in Białowieza Primeval Forest (Poland) Forest Ecology and Management. 2011;261:821–828. doi: 10.1016/j.foreco.2010.11.026. [Cross Ref]
Kunz et al. (2011) Kunz TH, De Torrez EB, Bauer D, Lobova T, Fleming TH. Ecosystem services provided by bats. Annals of the New York Academy of Sciences. 2011;1223:1–38. doi: 10.1111/j.1749-6632.2011.06004.x. [PubMed] [Cross Ref]
Latiff et al. (1999) Latiff A, Faridah Hanum I, Zainudin Ibrahim A, Goh MWK, Loo AHB, Tan HTW. On the vegetation and flora of Pulau Tioman, Peninsular Malaysia. Raffles Bulletin of Zoology. 1999;6:11–72.
Lee et al. (2009) Lee Y-F, Takaso T, Chiang T-Y, Kuo YM, Nakanishi N, Tzeng HY, Yasuda K. Variation in the nocturnal foraging distribution of and resource use by endangered Ryukyu flying foxes (Pteropus dasymallus) on Iriomotejima Island, Japan. Contributions to Zoology. 2009;78:51–64.
Little (2014) Little DP. A DNA mini-barcode for land plants. Molecular Ecology Resources. 2014;14:437–446. doi: 10.1111/1755-0998.12194. [PubMed] [Cross Ref]
Long & Racey (2007) Long E, Racey PA. An exotic plantation crop as a keystone resource for an endemic megachiropteran, Pteropus rufus in Madagascar. Journal of Tropical Ecology. 2007;23:1–11. doi: 10.1017/S0266467406003774. [Cross Ref]
Luskin (2010) Luskin MS. Flying foxes prefer to forage in farmland in a tropical dry forest landscape mosaic in Fiji. Biotropica. 2010;42:246–250. doi: 10.1111/j.1744-7429.2009.00577.x. [Cross Ref]
Mahmood-ul-Hassan et al. (2010) Mahmood-ul-Hassan M, Gulraiz TL, Rana SA, Javid A. The diet of Indian flying-foxes (Pteropus giganteus) in urban habitats of Pakistan. Acta Chiropterologica. 2010;12:341–347. doi: 10.3161/150811010X537927. [Cross Ref]
Marshall (1985) Marshall AG. Old World phytophagous bats and their food plants: a survey. Zoological Journal of the Linnean Society. 1985;83:351–369.
McConkey & Drake (2006) McConkey KR, Drake DR. Flying foxes cease to function as seed dispersers long before they become rare. Ecology. 2006;87:271–276. doi: 10.1890/05-0386. [PubMed] [Cross Ref]
McConkey & Drake (2007) McConkey KR, Drake DR. Indirect evidence that flying foxes track food resources among islands in a Pacific archipelago. Biotropica. 2007;39:436–440. doi: 10.1111/j.1744-7429.2007.00269.x. [Cross Ref]
McConkey & Drake (2015) McConkey KR, Drake DR. Low redundancy in seed dispersal within an island frugivore community. AoB Plants. 2015;7 doi: 10.1093/aobpla/plv088. plv088. [PMC free article] [PubMed] [Cross Ref]
Medway (1966) Medway L. The mammals. Observations on the fauna of Pulau Tioman and Pulau TulaiBulletin of the National Museum, Singapore. 1966;34:9–32.
Mildenstein et al. (2005) Mildenstein TL, Stier SC, Nuevo-Diego CE, Mills LS. Habitat selection of endangered and endemic large flying-foxes in Subic Bay, Philippines. Biological Conservation. 2005;126:93–102.
Mildenstein, Tanshi & Racey (2016) Mildenstein T, Tanshi I, Racey PA. Bats in the Anthropocene: conservation of bats in a changing world. Berlin Heidelberg: Springer International Publishing; 2016. Exploitation of bats for bushmeat and medicine; pp. 325–375.
Mohamed (2014) Mohamed NZ. PhD thesis. 2014. The role of nectar-feeding bats (Pteropodidae) in the pollination ecology of the genus Sonneratia at Setiu mangrove areas, Malaysia.
Mohd. Norfaizal et al. (2014) Mohd. Norfaizal G, Masrom H, Omar Y, Aisha-Farhana S. A preliminary flora survey in Gunung Kajang, Pulau Tioman, Pahang Darul Makmur, Peninsular Malaysia. Malaysian Applied Biology. 2014;43:17–23.
Nakamoto, Kinjo & Izawa (2007) Nakamoto A, Kinjo K, Izawa M. Food habits of Orii’s flying-fox, Pteropus dasymallus inopinatus, in relation to food availability in an urban area of Okinawa-jima Island, the Ryukyu Archipelago, Japan. Acta Chiropterologica. 2007;9:237–249. doi: 10.3161/1733-5329(2007)9[237:FHOOFP]2.0.CO;2. [Cross Ref]
Nakamoto, Kinjo & Izawa (2009) Nakamoto A, Kinjo K, Izawa M. The role of Orii’s flying-fox (Pteropus dasymallus inopinatus) as a pollinator and seed disperser on Okinawa-jima Island, the Ryukyu Archipelago, Japan. Ecological Research. 2009;24:405–414. doi: 10.1007/s11284-008-0516-y. [Cross Ref]
Nyhagen et al. (2005) Nyhagen DF, Turnbull SD, Olesen JM, Jones CG. An investigation into the role of the Mauritian flying fox, Pteropus niger, in forest regeneration. Biological Conservation. 2005;122:491–497. doi: 10.1016/j.biocon.2004.08.012. [Cross Ref]
Oleksy, Racey & Jones (2015) Oleksy R, Racey PA, Jones G. High-resolution GPS tracking reveals habitat selection and the potential for long-distance seed dispersal by Madagascan flying foxes Pteropus rufus. Global Ecology and Conservation. 2015;3:678–692. doi: 10.1016/j.gecco.2015.02.012. [Cross Ref]
Olival (2008) Olival KJ. PhD dissertation. 2008. Population genetic structure and phylogeography of southeast asian flying foxes: implications for conservation and disease ecology.
Ong (2000) Ong SY. BSc thesis. 2000. Kajian Ekologi Pteropus hypomelanus di Kampung Ayer Batang, Pulau Tioman, Pahang.
Pompanon et al. (2012) Pompanon F, Deagle BE, Symondson WOC, Brown DS, Jarman SN, Taberlet P. Who is eating what: diet assessment using next generation sequencing. Molecular Ecology. 2012;21:1931–1950. doi: 10.1111/j.1365-294X.2011.05403.x. [PubMed] [Cross Ref]
Quéméré et al. (2013) Quéméré E, Hibert F, Miquel C, Lhuillier E, Rasolondraibe E, Champeau J, Rabarivola C, Nusbaumer L, Chatelain C, Gautier L, Ranirison P. A DNA metabarcoding study of a primate dietary diversity and plasticity across its entire fragmented range. PLOS ONE. 2013;8(3):e58971 doi: 10.1371/journal.pone.0058971. [PMC free article] [PubMed] [Cross Ref]
Raye’ et al. (2011) Raye’ G, Miquel C, Coissac E, Redjadj C, Loison A, Taberlet P. New insights on diet variability revealed by DNA barcoding and high-throughput pyrosequencing: chamois diet in autumn as a case study. Ecological Research. 2011;26:265–276. doi: 10.1007/s11284-010-0780-5. [Cross Ref]
Razgour et al. (2011) Razgour O, Clare EL, Zeale MR, Hanmer J, Schnell IB, Rasmussen M, Gilbert TP, Jones G. High-throughput sequencing offers insight into mechanisms of resource partitioning in cryptic bat species. Ecology and Evolution. 2011;1:556–570. doi: 10.1002/ece3.49. [PMC free article] [PubMed] [Cross Ref]
Reiter & Curio (2001) Reiter J, Curio E. Home range, roost switching, and foraging area in a Philippine fruit bat (Ptenochirus jagori) Ecotropica. 2001;7:109–113.
Riaz et al. (2011) Riaz T, Shehzad W, Viari A, Pompanon F, Taberlet P, Coissac E. ecoPrimers: inference of new DNA barcode markers from whole genome sequence analysis. Nucleic Acids Research. 2011;39:e145. doi: 10.1093/nar/gkr732. [PMC free article] [PubMed] [Cross Ref]
Robinson (1984) Robinson JG. Diurnal variation in foraging and diet in the wedge-capped capuchin Cebus olivaceus. Folia Primatologica. 1984;43:216–228.
Scanlon et al. (2014) Scanlon AT, Petit S, Tuiwawa M, Naikatini A. High similarity between a bat-serviced plant assemblage and that used by humans. Biological Conservation. 2014;174:111–119.
Schmelitschek, French & Parry-Jones (2009) Schmelitschek E, French KO, Parry-Jones K. Fruit availability and utilisation by grey-headed flying-foxes (Pteropodidae: Pteropus poliocephalus) in a human-modified environment on the south coast of New South Wales, Australia. Wildlife Research. 2009;36:592–600.
Soberón, Llorente & Onate (2000) Soberón JM, Llorente JB, Onate L. The use of specimen-label databases for conservation purposes: an example using Mexican Papilionid and Pierid butterflies. Biodiversity and Conservation. 2000;9:1441–1466. doi: 10.1023/A:1008987010383. [Cross Ref]
Sritongchuay et al. (2014) Sritongchuay T, Gale GA, Stewart A, Kerdkaew T, Bumrungsri S. Seed rain in abandoned clearings in a lowland evergreen rain forest in southern Thailand. Tropical Conservation Science. 2014;7:572–585. doi: 10.1177/194008291400700314. [Cross Ref]
Srivathsan et al. (2016) Srivathsan A, Ang A, Vogler AP, Meier R. Fecal metagenomics for the simultaneous assessment of diet, parasites, and population genetics of an understudied primate. Frontiers in Zoology. 2016;13:17. doi: 10.1186/s12983-016-0150-4. [PMC free article] [PubMed] [Cross Ref]
Srivathsan et al. (2015) Srivathsan A, Sha J, Vogler AP, Meier R. Comparing the effectiveness of metagenomics and metabarcoding for diet analysis of a leaf-feeding monkey (Pygathrix nemaeus) Molecular Ecology Resources. 2015;15:250–261. doi: 10.1111/1755-0998.12302. [PubMed] [Cross Ref]
Start (1974) Start AN. PhD thesis. 1974. The feeding biology in relation to food sources of nectarivorous bats (Chiroptera: Macroglossinae) in Malaysia.
Stewart, Makowsky & Dudash (2014) Stewart AB, Makowsky R, Dudash MR. Differences in foraging times between two feeding guilds within Old World fruit bats (Pteropodidae) in southern Thailand. Journal of Tropical Ecology. 2014;30:249–257. doi: 10.1017/S0266467414000042. [Cross Ref]
Stier & Mildenstein (2005) Stier SC, Mildenstein TL. Dietary habits of the world’s largest bats: the Philippine flying foxes, Acerodon jubatus and Pteropus vampyrus lanensis. Journal of Mammalogy. 2005;86:719–728. doi: 10.1644/1545-1542(2005)086[0719:DHOTWL]2.0.CO;2. [Cross Ref]
Sudhakaran & Doss (2012) Sudhakaran MR, Doss PS. Food and foraging preferences of three pteropodid bats in southern India. Journal of Threatened Taxa. 2012;4:2295–2303. doi: 10.11609/JoTT.o2227.2295-303. [Cross Ref]
Tedman & Hall (1985) Tedman RA, Hall LS. The morphology of the gastrointestinal tract and food transit time in the fruit bats Pteropus alecto and P. poliocephalus (Megachiroptera) Australian Journal of Zoology. 1985;33:625–640. doi: 10.1071/ZO9850625. [Cross Ref]
Thomas (2009) Thomas DW. Analysis of diets of plant-visiting bats. In: Kunz TH, Parsons S, editors. Ecological and behavioural methods for the study of bats. Second edition The Johns Hopkins University Press; Baltimore: 2009.
Valentini et al. (2009) Valentini A, Miquel C, Nawaz MA, Bellemain EVA, Coissac E, Pompanon F, Gielly L, Cruaud C, Nascetti G, Wincker P, Swenson JE. New perspectives in diet analysis based on DNA barcoding and parallel pyrosequencing: the trnL approach. Molecular Ecology Resources. 2009;9:51–60. doi: 10.1111/j.1755-0998.2008.02352.x. [PubMed] [Cross Ref]
Weber et al. (2015) Weber N, Duengkae P, Fahr J, Dechmann D. KN, Phengsakul P, Khumbucha W, Siriaroonrat B, Wacharapluesadee S, Maneeorn P, Wikelski M, Newman S. High-resolution GPS tracking of Lyle’s flying fox between temples and orchards in central Thailand. Journal of Wildlife Management. 2015;79:957–968. doi: 10.1002/jwmg.904. [Cross Ref]
Williams et al. (2006) Williams N. SG, McDonnell MJ, Phelan GK, Keim LD, Van der Ree R. Range expansion due to urbanization: increased food resources attract Grey-headed Flying-foxes (Pteropus poliocephalus) to Melbourne. Austral Ecology. 2006;31:190–198. doi: 10.1111/j.1442-9993.2006.01590.x. [Cross Ref]
Win & Mya (2015) Win SS, Mya KM. The diet of the Indian Flying Fox Pteropus giganteus (Brünnich, 1782) (Chiroptera: Pteropodidae) in Myanmar—conflicts with local people? Journal of Threatened Taxa. 2015;7:7568–7572. doi: 10.11609/JoTT.o4178.7568-72. [Cross Ref]
Zeale et al. (2011) Zeale MR, Butlin RK, Barker GL, Lees DC, Jones G. Taxon-specific PCR for DNA barcoding arthropod prey in bat faeces. Molecular Ecology Resources. 2011;11:236–244. doi: 10.1111/j.1755-0998.2010.02920.x. [PubMed] [Cross Ref]
Zhang et al. (2014) Zhang J, Kobert K, Flouri T, Stamatakis A. PEAR: a fast and accurate Illumina Paired-End reAd MergeR. Bioinformatics. 2014;30:614–620. doi: 10.1093/bioinformatics/btt593. [PMC free article] [PubMed] [Cross Ref]

Articles from PeerJ are provided here courtesy of PeerJ, Inc