Skip to main content

Fungal communities in feces of the frugivorous bat Ectophylla alba and its highly specialized Ficus colubrinae diet



Bats are important long-distance dispersers of many tropical plants, yet, by consuming fruits, they may disperse not only the plant’s seeds, but also the mycobiota within those fruits. We characterized the culture-dependent and independent fungal communities in fruits of Ficus colubrinae and feces of Ectophylla alba to determine if passage through the digestive tract of bats affected the total mycobiota.


Using presence/absence and normalized abundance data from fruits and feces, we demonstrate that the fungal communities were significantly different, even though there was an overlap of ca. 38% of Amplicon Sequence Variants (ASVs). We show that some of the fungi from fruits were also present and grew from fecal samples. Fecal fungal communities were dominated by Agaricomycetes, followed by Dothideomycetes, Sordariomycetes, Eurotiomycetes, and Malasseziomycetes, while fruit samples were dominated by Dothideomycetes, followed by Sordariomycetes, Agaricomycetes, Eurotiomycetes, and Laboulbeniomycetes. Linear discriminant analyses (LDA) show that, for bat feces, the indicator taxa include Basidiomycota (i.e., Agaricomycetes: Polyporales and Agaricales), and the ascomycetous class Eurotiomycetes (i.e., Eurotiales, Aspergillaceae). For fruits, indicator taxa are in the Ascomycota (i.e., Dothideomycetes: Botryosphaeriales; Laboulbeniomycetes: Pyxidiophorales; and Sordariomycetes: Glomerellales). In our study, the differences in fungal species composition between the two communities (fruits vs. feces) reflected on the changes in the functional diversity. For example, the core community in bat feces is constituted by saprobes and animal commensals, while that of fruits is composed mostly of phytopathogens and arthropod-associated fungi.


Our study provides the groundwork to continue disentangling the direct and indirect symbiotic relationships in an ecological network that has not received enough attention: fungi-plants-bats. Findings also suggest that the role of frugivores in plant-animal mutualistic networks may extend beyond seed dispersal: they may also promote the dispersal of potentially beneficial microbial symbionts while, for example, hindering those that can cause plant disease.


Ecological networks have been the topic of extensive research, yet it seems that we still have important missing pieces of these complex puzzles [1]. For example, few studies have addressed the role of parasites in networks [2], while indirect interactions—those that occur when the interaction between two species is modified by a third one [3, 4]—are seldom addressed and poorly understood [5, 6]. Notwithstanding, indirect interactions are ubiquitous in ecological networks and play a critical role in shaping the bonds among species within communities [7,8,9,10,11]. While we typically believe that the survival of a species depends solely on the protection of its direct interactions, realistically it may also depend on how those interactions are shaped by species indirectly linked to each node in the network; the loss of any of these indirect links may have unforeseen effects on network functioning. Therefore, in considering the vulnerability of many species due to human activities, such as introduction of exotic species, habitat loss, introduction of parasites, and climate change, we need to acknowledge both the direct and indirect interactions of a taxon within ecological networks to fully understand the potential effects of its loss.

A type of ecological network that has received much attention is that which encompasses interactions between plants and their animal pollinators and seed dispersers; the interaction is mutually beneficial because animals help transport pollen and seeds, and in exchange obtain food (reviewed in [12]). Studies of plant-animal mutualistic networks have focused on understanding the relationship among interacting species [13, 14], the importance of particular taxa and network topology for maintaining network resilience [15,16,17,18], and how the topology of mutualistic networks may influence resilience [15,16,17] and even species diversity [15, 19]. However, these studies have neglected an important component of plant communities: their microbial symbionts. One such group of important microbial symbionts are endophytic fungi, which live within aerial tissues of plants without causing any visible negative impact. Even though some endophytes may be latent pathogens or saprotrophs [20], in many cases these endosymbionts provide benefits to the plant including protection against diseases and pests, plant growth, and reduction of drought stress [21]. As a result, endophytic fungi are regarded as critical components of any healthy plant community. Less is known about epiphytic or phyllosphere fungi, but there is some evidence of benefits to the host plant [22, 23].

The ability of a fungal species to disperse its spores (or other propagules such as hyphae, chlamydospores, and fruiting structures) is one of the factors that influence fungal diversity in a natural ecosystem [24, 25]. Fungi, including those that can become endophytic, can only disperse their spores short distances (a few centimeters at the most) using their own means (e.g., forcible ejection or discharge) [26, 27]. Consequently, they rely on other mechanisms for long distance dispersal, e.g., water, wind, and animals. Wind has been reported as the most efficient long-distance spore dispersal mechanism. However, in natural forests, especially old-growth, wind may not have a large influence in spore dispersal because trees and understorey vegetation provide a barrier to wind movement [28, 29]. Therefore, it is expected that other factors besides wind are influencing long-distance dispersal of fungi in natural tropical forests. Considering that animals are capable of long-distance dispersal of plant seeds, it is possible that their role extends to the dispersal of fungal spores and other propagules.

Many animals may disperse fungi directly by eating mushrooms and then defecating the spores; or indirectly, by eating other plant parts that contain these fungi. The direct consumption of fungal fruiting bodies, or mycophagy, and spore dispersal has been described several times in insects, rodents, marmosets, and other mammals [30,31,32]. In some cases it was reported that fungal spores survive and their germination is improved after passing the digestive tract of truffle-eating rodents or other ground-dwelling animals [32]. However, indirect fungal propagule dispersal is woefully unknown. We use bats as a model to understand this interaction because these mammals are important long-distance dispersers of many tropical plants [33]. Yet by consuming fruits, bats may disperse not only the plant’s seeds, but also the fungi that are contained in those fruits. Bats may be particularly good dispersers of fungi because they fly long distances each night, defecate during flight, and may retain viable propagules for long periods of time. Unlike birds, fruit-eating bats also venture frequently into deforested areas that may otherwise lack input of beneficial fungal spores [34,35,36,37]. This study represents a first step towards identifying an interaction that may have consequences for the preservation of healthy tropical ecosystems.

In this study we aimed to explore the hypothesis of an indirect mutualistic relationship between fruit-eating animals, specifically bats, and symbiotic fungi (with emphasis on endophytes) that grow within the tissues of fruits that bats eat. Since fungi can develop in any plant tissue as endophytes, including fruits [38, 39], it is presumed that bats may disperse fungal propagules that are consumed from these structures. To begin to investigate the poorly examined premise that bats are also long-distance dispersers of fungi, the main objectives of this study were to (i) characterize the fungal communities from the fruits of Ficus colubrinae and determine whether the same species are present in the feces of Ectophylla alba; and (ii) determine if at least some fungi survive the digestive tract of bats. In this project we aim to address these basic objectives, yet many questions will remain regarding the relationship between frugivores, plants, and fungi. We hope our answers will begin to shed light on this potential interaction and hopefully foster further scrutiny.


Ectophylla alba’s main diet

Out of the total plant ITS nrDNA sequences amplified from the fecal samples, 82% (σ = 17) matched to several Ficus spp. with percent similarities of 90–99%, including F. colubrinae with 96% (GenBank accession number EU081760). However, the only two matching species that are present in La Selva Biological Station are F. colubrinae and F. costaricana (La Selva Florula Digital, The remaining 18% of plant ITS sequences corresponded to various species of microscopic green algae (e.g., Chlamydomonas spp., Parachlorella spp., and Trentepohlia spp., among others). Most (i.e., 99.8%; σ = 0.17) of the ITS sequences from fruits matched to the same species as those found in feces.

Comparison of fruit and fecal culture-independent fungal communities

The total number of fungal ASVs identified from metabarcoding of 9 pooled fruit samples (18 total fruits; sequencing failed for one of the fruit samples) and 13 bats was 460 and 1025, respectively (Additional file 1: Tables S1–S3). The phylum with the highest number of ASVs in both feces and fruits was Ascomycota, followed by Basidiomycota (Additional file 1: Fig. S1); a large percentage of ASVs (ca. 30%) did not match to any known fungal phylum. Figure 1 shows the most abundant (percent relative abundance of ASVs) classes and orders in the fecal and fruit samples. Excluding unidentified ASVs, fecal and fruit samples were dominated by the classes Agaricomycetes and Dothideomycetes, respectively (Fig. 1a, c, and Additional file 1: Fig. S1). The most frequent orders in feces (excluding unclassified ASVs) were Pleosporales, followed by Polyporales, Chaetothyriales, Hypocreales, and Agaricales (Fig. 1b, d, and Additional file 1: Fig. S1). Fruits were dominated by Pleosporales, followed by Hypocreales, Chaetothyriales, Capnodiales, and Glomerellales. Lastly, the genera with the greatest number of ASVs in the fecal samples were Malassezia, followed by Wallemia, Aspergillus, Fusarium, and Pyxidiophora (Additional file 1: Fig. S1). In contrast, Malassezia, Pyxidiophora, Colletotrichum, Ochroconis, and Diaporthe dominated the fruit mycobiota.

Fig. 1
figure 1

Barplots of fungal ASVs by class and order where panels (a) and (b) represent individual sample composition and panels (c) and (d) reflect overall group composition. Averaged taxa abundance per group is shown in c and d. Unidentified fungal taxa which could not be assigned to any other taxonomical category were aggregated as “unid” (light blue in all graphs); fungal taxa that could be assigned to Ascomycota but not further are aggregated as “unid_Ascomycota.” Only the 10 most abundant taxa are shown; uncommon taxa are combined in the category “others.”

While a larger number of fungal taxa was observed in bat feces than in fruit, alpha diversity did not differ significantly between the two communities (Fig. 2a) even when considering samples collected by tree (Additional file 1: Fig. S2a). However, fungal communities of fruits and feces were significantly different (presence/absence data: R = 0.07, F1,20 = 1.51, P < 0.001; normalized abundance data: R = 0.08, F1,20 = 1.84, P = 0.05) and there was no significant association between the fungal communities of bats and fruits collected from the same trees (Additional file 1: Figs. S2b–d). There were 866 unique fungal ASVs present in bat feces and 301 in fruits (Fig. 2b), resulting in overall distinct communities as visualized in the NMDS ordination (Fig. 2c, d); 159 ASVs overlapped. Both communities include Cystofilobasidium, Fusarium, Geranomyces, Malassezia, Pyxidiophora, and Wallemia (Additional file 1: Tables S1–S3).

Fig. 2
figure 2

Alpha diversity (a), Venn diagram (b), and NMDS based on Jaccard distance for presence/absence data (c) and NMDS based on Bray–Curtis distance for normalized abundance data (d). a: ns non-significant difference. c, d the values within the plot represent the stress values, and ellipses the 95% confidence interval. Bat and fruit illustrations provided by Silvia Chaves Ramírez

According to the LDA scores for presence/absence data (Fig. 3a), indicator taxa for bat feces belong to the phylum Basidiomycota (i.e., Agaricomycetes: Polyporales and Agaricales) and the ascomycetous class Eurotiomycetes (i.e., Eurotiales, Aspergillaceae). For fruits, there are several unidentified taxa which constitute the larger contributors to differentiating their fungal communities from those found in feces. Other indicator taxa in fruits belong in the phylum Ascomycota (i.e., Dothideomycetes: Botryosphaeriales, Phyllostictaceae, Phyllosticta; Laboulbeniomycetes: Pyxidiophorales, Pyxidiophora; and Sordariomycetes: Glomerellales, Colletotrichum). In contrast, when using normalized abundance (Fig. 3b), the analyses resulted in more and somewhat different indicator taxa in bat feces. For example, many unidentified Basidiomycota and Ascomycota groups, in addition to Malasseziomycetes (i.e., Malasseziales, Malassezia), Wallemiomycetes (i.e., Wallemiales, Wallemia), Sordariomycetes (i.e., Hypocreales, Fusarium), Laboulbeniomycetes (i.e., Pyxidiophorales, Pyxidiophora), and Dothideomycetes (i.e., Pleosporales), among others. In fruits, indicator taxa belong in Tremellomycetes (i.e., Cystofilobasidiales, Cystofilobasidium), Saccharomycetes (i.e., Saccharomycetales, Wickerhamomyces), and Laboulbeniomycetes (i.e., Pyxidiophorales, Pyxidiophora).

Fig. 3
figure 3

LDA score is the linear discriminant analysis score in LEfSe (a: results based on presence/absence data, b: results based on normalized abundance data). The letters represent taxonomic classifications: phylum (p), class (c), order (o), family (f), genus (g) and species (s). Unidentified fungal taxa that could be assigned to a specific taxonomic level are aggregated as “x_unid”. Bat and fruit illustrations provided by Silvia Chaves Ramírez

Comparison of fruit and fecal culture-dependent fungal communities

Fungal colonies grew successfully from 4 fruits samples (11 isolates) and 10 fecal samples (30 isolates). The taxa present in fruit isolates included Colletotrichum fruticola, C. siamense, Diaporthe sp., D. cf. hongkongensis, Fusarium concentricum, and Neopestalotiopsis saprophytica; all matching ASVs that were also found in the culture-independent (metabarcoding) diversity analyses (Additional file 1: Table S3 and Additional file 1: Fig. S3). In fecal samples, the culture-dependent diversity was represented by Fusarium concentricum, F. waltergamsii, Mucor irregularis, Neopestalotiopsis saprophytica, and Pseudopestalotiopsis simitheae, also corresponding to ASVs obtained from the metabarcoding analyses. Colletotrichum spp. and Diaporthe spp. were only present in fruits in both culture-dependent and independent analyses. Neopestalotiopsis saprophytica was found in both fruits and feces. Fusarium concentricum dominated bat feces; however, some ASVs and cultures that matched this taxon were also observed in fruits. Fusarium waltergamsii grew only from fecal samples. However, metabarcoding data suggest that this taxon was also present in fruits. Mucor irregularis was only found in fecal samples in both culture-dependent and independent analyses.

Analyses of putative ecological roles of the fungal taxa

Based on results from FUNGuild annotation tool with modified assignments (see Methods section), in Fig. 4 we highlight the relative contribution of ASVs per sample according to putative ecological guild and based on presence/absence (Fig. 4a) and normalized abundance (Fig. 4b). Overall, excluding the unclassified/unidentified taxa, fecal samples were dominated by saprobes (e.g., Agaricales, Polyporales, and Eurotiales), phytopathogens (i.e., Fusarium concentricum), and animal (not insect) commensals (e.g., Malasseziales and Wallemiales), while fruit mycobiota was dominated by phytopathogens (e.g., Botryosphaeriales and Glomerellales), saprobes, and arthropod-associated taxa (i.e., Pyxidiophora and Wickerhamomyces). The LDA scores using presence/absence and normalized abundance (Fig. 3) showed that the indicator taxa in feces were mostly groups classified as saprobes and animal commensals. In contrast, the indicator taxa in fruits include phytopathogens and arthropod-associated fungi.

Fig. 4
figure 4

Heatmaps showing the relative contribution of taxa with a known function within each sample (columns) in bat feces and fruit communities based on presence/absence data (a) and normalized abundance (b). Sample abundance for (a) is based on the number of taxa found within each sample per function, whereas sample abundance for (b) is based on the sum of proportions of taxa for all taxa of a given function within samples. We also include a total column for each community, which shows the sum of all identified taxa within samples (a) or the sum of proportions (b) for a given function; grey shades provide an estimate of relative abundance. Bat and fruit illustrations provided by Silvia Chaves Ramírez


One of the hypotheses we posed in this study was that bats may disperse, through fecal deposition, the fungi that inhabit the fruits they consume. The combination of qualitative (presence/absence) and quantitative (normalized abundance) data suggest that approximately 38% of the original fruit mycobiota (i.e., species composition) remains in the fecal samples, even though the overall comparison of fruit and fecal fungal assemblages yielded marginally significant differences (Fig. 2). Through culture-dependent techniques we show that some of these fungal species originating from the fruit are still viable in feces, suggesting bats may be capable of dispersing fungi over long distances. Our characterization of the mycobiota through metabarcoding indicates the presence of fungal DNA (i.e., ITS nrDNA) in fruits that is similarly observed in bat feces, but we cannot conclude that all these species are still viable after passing through the digestive system. Our study found that out of the five culturable fungal species in F. colubrinae fruits, only two of those were recovered in cultures of bat feces. However, some of the ASVs identified by metabarcoding that were present in the fruit mycobiota but were not represented in cultured isolates did successfully grow from fecal samples. A noteworthy example is Fusarium waltergamsii, which was present in both fruits and feces, but was only cultured from feces, suggesting that passage through the digestive tract may help some fungi germinate and grow [30,31,32]. However, even though ITS nrDNA continues to be the most utilised option in fungal metabarcoding studies, caution should be placed in taxonomy assignments using this marker, as studies point to the limitations in species circumscription [40].

Many studies which compare both culturable and unculturable mycobiota from environmental samples have found similar trends where only a very small portion (< 5%) of the total fungal diversity is recovered in cultures [41,42,43]. This large difference is mainly due to competition (e.g., presence of low-abundant and slow-growing microorganisms that may be outcompeted by high-abundant and fast-growing species), obligate biotrophy (e.g., the fungus can only grow on a living host or partner), and substrate specificity (e.g., failure to grow on conventional media because of inappropriate conditions of pH, temperature, redox state, or availability of essential nutrients) [44,45,46]. Notwithstanding, our results and those of others indicate that several fungal taxa can grow after passing through the digestive tract of vertebrates [32, 47]. Additional RNA, transcriptomic, or proteomic analyses may help us better understand the viable mycobiota and the functions and interactions between microbial species [48,49,50,51]. Understanding which species can successfully grow after passage through animal guts will provide a clearer picture of the role of frugivores in the dispersal of fungal endophytes.

Comparisons of the mycobiota of fruits and bat feces provide preliminary evidence on how frugivores may be affecting plants exposed to microbial communities within feces, how surviving fungal taxa may affect frugivores, and how fungi may affect the interaction between fruits and frugivores. In our study, the differences in fungal species composition between the two communities (fruits vs. feces, Fig. 2) also reflected changes in functional diversity as estimated by ecological guild analyses (Fig. 4). Indicator taxa analyses (Fig. 3) revealed that the core fecal mycobiota is comprised of saprobe fungi, while that of fruits is constituted mostly of phytopathogens and arthropod-associated fungi. Therefore, Ficus colubrinae seeds dispersed by frugivores may benefit from a reduction in the number of potentially plant pathogenic taxa inherent to the fruit body. These interactions could increase seed survival which is often very low, primarily due to pathogenesis [52]. We can then infer that if F. colubrinae fruits are not consumed by frugivores and they just fall directly from the tree, potential plant pathogens contained in fruits will remain near the parent tree and may cause disease in seedlings. These results support the Janzen–Connell hypothesis [53, 54] or the Theory of pest pressure [55], which suggest that specialized natural enemies (i.e., plant pathogens) decrease survival of seedlings that are in high densities beneath the parent tree, thus giving locally rare species an advantage. An alternative scenario is that if fruit-eating bats defecate under or near F. colubrinae trees, or under or near their roosts (leaves of Heliconia spp.) [56], passage through the digestive tract may still result in decreased abundance of potential plant pathogens, which would reduce disease. Therefore, consumption of fruits by animals (i.e., bats) could not only benefit the plant by dispersing its seeds, but also by reducing the amount of pathogen inoculum and thus escape disease pressure [57].

A noteworthy and novel finding was Pyxidiophora, an obligate mycoparasitic fungal genus [58] found in both fruit and fecal mycobiota. The biology of this genus is poorly studied, but a few studies show that the sexual spores (ascospores) attach to phoretic mites of bark beetles [59]. The fungal spores are vectored by the insects, which then land on other fungi to become mycoparasitic [58]. As we found Pyxidiophora to be a core constituent of Ficus colubrinae fruits (Fig. 3), we hypothesize that this fungus is associated with the fig pollinating wasps (Pegoscapus spp., Agaonidae, Chalcidoidea) which were abundant inside the fruits we collected [60]. We also demonstrate that Pyxidiophora DNA is present in bat feces, suggesting that Ectophylla alba may also disperse these fungi. Previous studies report that Pyxidiophora is a common viable fungus parasitizing dung fungi [58]. Our study expands the list of potential arthropod hosts, suggesting that Pyxidiophora may not only be restricted to mites and beetles. Another fungus that may be associated with fig pollinating wasps is Wickerhamomyces, which has been reported from the midgut and gonads of several insects [61, 62].

While we still lack sufficient information to predict how certain fungal taxa found in the tissues of fruits may be affecting fruit-eating animals, some of the fungi we recorded in feces that were also present in fruits are considered common animal endosymbionts associated with healthy guts [63,64,65], such as Malasseziales and Wallemiales. In addition, some of the fungi that were found in feces are known to aid in digestion by improving the metabolizable energy of plant-based diets due to their ability to produce CAZymes and endo-beta(1,4)-xylanases (e.g., Aspergillus and related taxa, [66]) and plant cell-wall degrading enzymes (e.g., Basidiomycota: Polyporales, or Ascomycota: Hypocreales, Fusarium spp. [67]), in addition to the production and release of carotenoids, lipids, and coenzyme Q10 into the intestine (e.g., Cystofilobasidium and Sporidiobolales, [68,69,70,71]).

Finally, fungi may also affect the interaction between plants and frugivores. Fungi can produce secondary metabolites and volatile organic compounds (VOCs) which can affect fruit palatability and attraction [72,73,74,75,76]. For example, Fusarium verticillioides (related to F. concentricum, a species found in our fruit samples) and other Fusarium spp. have been found to produce VOCs that attract insects [77, 78]. By extension, we hypothesize that endophytic fungal communities may change the chemical composition of fruits and hence preferences in frugivores. While unquestionably relevant for understanding mutualistic networks, this topic remains completely unexplored.

This poorly studied interaction among fungal endophytes, fruits, and frugivores suggests there is a critical component of plant-animal ecological networks that urgently requires further scrutiny. This additional link complicates our understanding of mutualistic network dynamics and has implications for models developed thus far. For example, we often regard frugivores as somewhat equally responsible for promoting seed dispersal, yet differing foraging styles and physiological conditions among fruit-eating species can potentially affect dispersal distance and viability of propagules, fungal or otherwise [79, 80]. By affecting fruit palatability and overall plant fitness, the mycobiota may also be highly responsible for the success, or failure, of certain plant species, which ultimately modifies the composition of many ecological networks and interactions therein.

Study of the mycobiota has experienced a major increase in recent years, predominantly with the advent of genetic tools, and evidence is accumulating on the many roles that fungi play in natural ecosystems. Many of the studies conducted so far have shown a diverse fungal community in plants, yet surprisingly, very little is known about mycosymbionts in fruits (but see [38, 81,82,83]) despite the obvious role of fruits for plant fitness, and no studies to date have assessed how vertebrate consumption can affect fungal endophytes. Consequently, the role of endophytic fungi in mutualistic networks has been, until now, largely ignored. Our study shows that fungal endophytes are ubiquitous within fruits, and as such may be important components of plant-animal networks. Their ubiquity in plant tissues and the potential role that plant-eating organisms can play in dispersing the fungal propagules, suggest that further studies into mutualistic interactions should consider greater focus on endosymbionts.


We conclude that fungal communities in Ectophylla alba feces and Ficus colubrinae fruits are fundamentally different, with about 38% of the ASVs shared between the two sample groups. As suggested by previous studies, we confirmed that a combination of qualitative (presence/absence) and quantitative (normalized abundance) data provides a more complete depiction of the fungal communities studied. We also show that there are several viable fungi and the presence of fruit ASVs in fecal samples, suggesting bats may be important dispersers for those fungi. The fruit mycobiota is dominated mostly by fungi with potential plant pathogenic activity, whereas fecal samples are dominated by saprobes. We also established, through metabarcoding, that Ectophylla alba’s main diet is based on Ficus colubrinae fruits. Our findings indicate that the role of frugivores in plant-animal mutualistic networks may extend beyond seed dispersal: they may also promote the dispersal of potentially beneficial microbial symbionts while simultaneously hindering those that can cause plant disease.


Study system

Fresh ripe fig fruits (Ficus colubrinae, Moraceae) and Honduran White Bat (Ectophylla alba, Chiroptera: Phyllostomidae) fecal samples were collected for fungal community analyses. It has been reported that E. alba feeds almost exclusively from F. colubrinae fruits [84]. Ectophylla alba is known only from Honduras, Nicaragua, Costa Rica, and western Panama [85]. In 2008, IUCN elevated this bat species to a near threatened Red List category [86]. Populations of this bat species have been declining due to urbanization and strong habitat (i.e., Heliconia leaves from intermediate secondary succession forests that are used for roosting) and diet (i.e., F. colubrinae fruits) preferences [86, 87]. Ficus colubrinae is an understory tree distributed from Mexico to Colombia, fruits throughout the year (La Selva Florula Digital,, and is a food source for many bats and birds [88].

Fruit and fecal sample collection

Samples were collected in La Selva Biological Station (Sarapiquí, Heredia, Costa Rica). Two trees that were fruiting at the time of the fieldwork (February 2015) were chosen for fruit and fecal sample collection. We captured bats and collected fruits from the same trees as to increase the chances that the fecal samples came from the fruits consumed in that tree. Both sample types were subjected to gene-amplicon targeted sequencing (metabarcoding) and culture analyses. Each collection consisted of four adjacent ripe fruits (i.e., from the same cluster) placed into individual Ziploc bags. The clusters were picked from five randomly selected but reachable areas of each tree. Out of the four fruits per cluster, two were placed in 2-mL Eppendorf microtubes with silica gel and frozen at -20 °C for posterior DNA extraction. The remaining two were used for isolation into pure culture (see “Culturable fungi isolation and identification” section). Hereafter, trees are labeled as “Fc1” and “Fc2”, respectively, and each fruit sample from each tree as e.g., Fc1_1, Fc1_2, and so on. In total, we obtained 10 pooled fruit samples (representing 20 individual fruits) for metabarcoding and 20 fruits for culture analyses; 40 fruits for the entire study (2 trees × 5 clusters × 4 fruits = 40 fruits).

From trees Fc1 and Fc2, five and eight bats, respectively, were captured with mist nets (Ecotone, Poland) and immediately placed in sterilized cloth bags. To avoid cross-contamination, we cleaned our hands with an alcohol-based hand gel before releasing every bat from the net. When bats defecated, sterilized cotton swabs were used to obtain the fecal sample from the cloth bag. The approximate volume collected was a 4–5-mm-diam. pellet. Half of each sample was placed in 2-mL Eppendorf tubes, in Ziploc bags with silica gel, and then placed in a -20 °C freezer for subsequent metabarcoding analyses. The other half was placed in sterile Eppendorf tubes for later same-day culturing.

Culturable fungi isolation and identification

Fruits were divided into five equal pieces and placed onto CMD + (BBL™ corn-meal-agar + 2% dextrose + antibiotic) 9-mm Petri dishes [89]. An antibiotic solution was added to the media to eliminate bacteria. From each fecal sample swab, five points of inoculation were made onto each CMD + Petri plate. The plates were incubated for several days (up to 2 weeks) at room temperature and the emerging colonies were subcultured to obtain pure isolates.

Genomic DNA from pure fungal cultures was extracted with PrepMan Ultra (Life Technologies, Waltham, MA, U.S.A.). The Internal Transcribed Spacer (ITS) and a region of the Large Subunit (28S) of the nuclear ribosomal DNA were amplified in one reaction, using the ITS5-forward (GGAAGTAAAAGTCGTAACAAGG) and LR5-reverse (TCCTGAGGGAAACTTCG) primers [90]. ITS is the official fungal barcode [40] and gives an approximate species identification. Polymerase chain reaction (PCR) conditions and protocols are described in previous publications [89]. PCR products were purified and sequenced at Macrogen U.S.A. Assembly of forward and reverse strands and sequence alignment were done in Geneious v 10.2.3 ( BLASTn algorithm was performed in Geneious with retrieve from the UNITE v 2020 database. Best hits were compared and sequence identity of > 99% was used for taxonomy assignment.

Metabarcoding of fruit and fecal samples

Genomic DNA from whole fruits and feces was extracted with the following protocol: fruits or feces were placed in the -20 °C freezer for at least one day and then, when ready to extract, placed into new tubes prefilled with 500 µm garnet beads and a 6 mm zirconium grinding satellite bead (OPS Diagnostics LLC, NJ, U.S.A.). For sample homogenization, a FastPrep® instrument (Zymo Research, Irvine, CA, U.S.A.) was used at maximum speed (6.5 m/s) for 1 min. 750 µl of Qiagen® Lysis Buffer AP1 and 6 ul of QiaGen® RNase-A were added to each tube and incubated overnight at 65 °C. Total DNA was extracted using the Qiagen® DNeasy Plant Mini Kit according to manufacturer’s instructions.

PCR amplicons of the ITS2 nrDNA region using fungal-specific primers fITS7-forward (GTGARTCATCGAATCTTTG) and ITS4-reverse (TCCTCCGCTTATTGATATGC) [91] were tagged and multiplexed for paired-end sequencing on the Illumina MiSeq 2 × 300 platform at MRDNA (, Shallowater, TX, U.S.A.). The same ITS2 primers can also amplify some plant DNA [91] and thus can be used to confirm the bat’s main diet. Three PCR stochastic replicates were pooled and purified using calibrated AMPure XP beads (Beckman Coulter Life Sciences, Indianapolis, IN, U.S.A.). PCR from a pure fungal culture of a Trichoderma koningiopsis and UltraPure™ DNase/RNase-Free Distilled Water (Thermo Fisher Scientific, Waltham, MA, U.S.A.) were used as positive and negative/blank controls, respectively, in quality control for the downstream bioinformatics. All raw sequencing data have been submitted to the NCBI Sequence Read Archive (SRA) database under the BioProject ID PRJNA759639.

Metabarcoding bioinformatics and fungal species identification

Cutadapt v 2.3 [92] in Python v 3.7.10 was used to remove primers. Dada2 v 1.21.0 [93] in RStudio v 1.4.1717 was used for quality inspection and filtering, trimming, merging paired-ends, sample or Amplicon Sequence Variant (ASV) inference, and chimera removal. Forward and reverse sequences were filtered and trimmed where the quality score dropped to < 20 (i.e., forward trimmed at nucleotide position 270 and reverse at 210), and with a maximum number of expected errors (maxEE) set to 2. Sequences were clustered into ASVs [94] and then filtered for chimeras. To assign taxonomy, ASVs were subjected to similarity searches in the UNITE v 2020 curated and quality-checked database [95] using DECIPHER v 2.0 [96]. Each name was verified manually for nomenclatural accuracy either in Index Fungorum or Mycobank.

Fungal diversity and community analyses

We analyzed the microbial communities found in fruits and bat feces using the package microeco v 0.6.5 [97] in R v 4.1.1. We transformed the ASV table to a matrix of presence/absence data (qualitative method) and estimated alpha (observed) and beta (Jaccard distance) diversity. Studies have suggested that using both quantitative and qualitative diversity measures will often be critical for understanding the factors that affect microbial diversity [98]. Therefore, to complement the results generated from presence/absence data, we also performed beta diversity analyses using abundance data (quantitative method) transformed to proportions [99]. This was accomplished by dividing the number of reads for each ASV in a sample (bat feces or fruit) by the total number of reads in that sample [99]. The distance matrix was constructed using the Bray–Curtis index. Alpha diversity (observed) was compared between the two communities using a t-test. Differences between the communities were determined based on ordination, using a non-metric multidimensional scaling (NMDS; [100]), and group distance, using a permutational multivariate analysis of variance (perMANOVA; [101]). Additionally, we identified which taxa might help us explain the differences between bat and fruit communities using the linear discriminant analysis (LDA) effect size (LEfSe) method [102]. We retained taxa with LDA scores > 4. In general, we placed more weight on the results from presence/absence data [100, 103, 104] over that from normalized abundance because we sought to determine whether the same fungal taxa in fruits were present in feces, and because of the intrinsic issues with misestimation of abundance in microbiome studies [99, 104,105,106,107,108].

Estimation of abundance of fungal taxa according to their ecological roles

We explored if there was a trend in the relative abundance and presence/absence of ASVs with specific ecological roles in fruits and feces. A putative ecological or functional role was assigned by first parsing fungal community datasets by ecological guild using FUNGuild annotation tool [109] implemented in Python v 3.6 through the supercomputer Kabré (CNCA-CONARE, Costa Rica), and then manually checking, refining, and modifying the assignments (following the approach in [110]). Since many of the assignments provided by FUNGuild were ambiguous or incorrect, and other taxa had no assignments at all, we used the following modified putative roles: saprobe, plant pathogen, entomopathogen, animal (not arthropod) commensal, insect symbiont (not pathogen), mycotroph or fungicolous, lichen-forming, epiphyte, animal (not arthropod) pathogen. At the time we collected the fruits, there were no disease symptoms or necrotrophy. Therefore, all the inferred ecological guilds refer to a hypothesized cryptic role [111, 112].

With the data on ecological guilds, presence/absence, and normalized abundance for taxa, we then constructed two separate heatmaps. In these maps, we only included taxa for which more than 10 hits were recorded overall. For data on presence/absence we estimated abundance of a given function based on the number of ASVs with that function that were found in a specific sample (i.e., bat feces or fruits). The heatmap was created by plotting abundance of each function for each sample. For normalized abundance, the heatmap was created by plotting the sum of proportions (frequency) for all taxa of a given function within samples.

Availability of data and materials

All code and raw data have been stored in the GitHub repository ( and the NCBI Sequence Read Archive (SRA) database under the BioProject ID PRJNA759639.



Amplicon sequence variant


Corn-meal agar + dextrose


Internal transcribed spacer


Linear discriminant analysis


LDA effect size method


Maximum number of expected errors


Non-metric multidimensional scaling


Polymerase chain reaction


Permutational multivariate analysis of variance


Volatile organic compound


  1. Olesen JM, Bascompte J, Dupont YL, Elberling H, Rasmussen C, Jordano P. Missing and forbidden links in mutualistic networks. Proc Biol Sci. 2011;278:725–32.

    PubMed  Google Scholar 

  2. Lafferty KD, Allesina S, Arim M, Briggs CJ, De Leo G, Dobson AP, et al. Parasites in food webs: the ultimate missing links. Ecol Lett. 2008;11:533–46.

    PubMed  PubMed Central  Article  Google Scholar 

  3. Levine SH. Competitive interactions in ecosystems. Am Nat. 1976;110:903–10.

    Article  Google Scholar 

  4. Holt RD. Predation, apparent competition, and the structure of prey communities. Theor Popul Biol. 1977;12:197–229.

    CAS  PubMed  Article  Google Scholar 

  5. Miller TE, TerHorst CP. Indirect effects in communities and ecosystems. In: Gibson D, (ed). Oxford Bibliogr Ecol. New York: Oxford University Press; 2012.

  6. Michalet R, Chen SY, An LZ, Wang XT, Wang YX, Guo P, et al. Communities: Are they groups of hidden interactions? J Veg Sci. 2015;26:207–18.

    Article  Google Scholar 

  7. Blanc LA, Walters JR. Cavity excavation and enlargement as mechanisms for indirect interactions in an avian community. Ecology. 2008;89:506–14.

    PubMed  Article  Google Scholar 

  8. Irwin RE. The consequences of direct versus indirect species interactions to selection on traits: Pollination and nectar robbing in Ipomopsis aggregata. Am Nat. 2006;167:315–28.

    PubMed  Article  Google Scholar 

  9. Poelman EH, Gols R, Snoeren TAL, Muru D, Smid HM, Dicke M. Indirect plant-mediated interactions among parasitoid larvae. Ecol Lett. 2011;14:670–6.

    PubMed  Article  Google Scholar 

  10. Wootton JT. The nature and consequences of indirect effects in ecological communities. Annu Rev Ecol Syst. 1994;25:443–66.

    Article  Google Scholar 

  11. Strauss SY. Indirect effects in community ecology: their definition, study and importance. Trends Ecol Evol Evol. 1991;6:206–10.

    CAS  Article  Google Scholar 

  12. Bascompte J. Mutualistic networks. Front Ecol Environ. 2009;7:429–36.

    Article  Google Scholar 

  13. Petanidou T, Kallimanis AS, Tzanopoulos J, Sgardelis SP, Pantis JD. Long-term observation of a pollination network: fluctuation in species and interactions, relative invariance of network structure and implications for estimates of specialization. Ecol Lett. 2008;11:564–75.

    PubMed  Article  Google Scholar 

  14. Carnicer J, Jordano P, Melian CJ. The temporal dynamics of resource use by frugivorous birds: a network approach. Ecology. 2009;90:1958–70.

    PubMed  Article  Google Scholar 

  15. Allesina S, Tang S. Stability criteria for complex ecosystems. Nature. 2012;483:205–8.

    CAS  PubMed  Article  Google Scholar 

  16. Bascompte J, Stouffer DB. The assembly and disassembly of ecological networks. Philos Trans R Soc B-Biol Sci. 2009;364:1781–7.

    Article  Google Scholar 

  17. Okuyama T, Holland JN. Network structural properties mediate the stability of mutualistic communities. Ecol Lett. 2008;11:208–16.

    PubMed  Article  Google Scholar 

  18. Mello MAR, Marquitti FMD, Guimaraes PR, Kalko EKV, Jordano P, de Aguiar MAM. The modularity of seed dispersal: differences in structure and robustness between bat- and bird-fruit networks. Oecologia. 2011;167:131–40.

    PubMed  Article  Google Scholar 

  19. Bascompte J, Jordano P. Plant–animal mutualistic networks: the architecture of biodiversity. Annu Rev Ecol Evol Syst. 2007;38:567–93.

    Article  Google Scholar 

  20. Porras-Alfaro A, Bayman P. Hidden fungi, emergent properties: endophytes and microbiomes. Annu Rev Phytopathol. 2011;49:291–315.

    CAS  PubMed  Article  Google Scholar 

  21. Rodriguez RJ, White JF, Arnold AE, Redman RS. Fungal endophytes: diversity and functional roles. New Phytol. 2009;182:314–30.

    CAS  Article  PubMed  Google Scholar 

  22. Andrews JH. Biological control in the phyllosphere. Annu Rev Phytopathol. 1992;30:603–35.

    CAS  PubMed  Article  Google Scholar 

  23. Lindow SE. Phyllosphere microbiology: a perspective. In: Bailey M, Lilley A, Timms-Wilson T, Spencer-Phillips P, editors. Microb Ecol Aer Plant Surfaces. Oxfordshire: CAB International; 2006. p. 1–20.

    Google Scholar 

  24. Kohn LM. Mechanisms of fungal speciation. Annu Rev Phytopathol. 2005;43:279–308.

    CAS  PubMed  Article  Google Scholar 

  25. Persoh D. Plant-associated fungal communities in the light of meta’omics. Fungal Divers. 2015;75:1–25.

    Article  Google Scholar 

  26. Roper M, Pepper RE, Brenner MP, Pringle A. Explosively launched spores of ascomycete fungi have drag-minimizing shapes. Proc Natl Acad Sci USA. 2008;105:20583–8.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  27. Galante TE, Horton TR, Swaney DP. 95% of basidiospores fall within 1 m of the cap: a field- and modeling-based study. Mycologia. 2011;103:1175–83.

    PubMed  Article  Google Scholar 

  28. Chen JQ, Franklin JF, Spies TA. Contrasting microclimates among clear-cut, edge, and interior of old-growth douglas-fir forest. Agric For Meteorol. 1993;63:219–37.

    Article  Google Scholar 

  29. Milleron M, de Heredia U, Lorenzo Z, Perea R, Dounavi A, Alonso J, et al. Effect of canopy closure on pollen dispersal in a wind-pollinated species (Fagus sylvatica L.). Plant Ecol. 2012;213:1715–28.

    Article  Google Scholar 

  30. Epps MJ, Arnold AE. Diversity, abundance and community network structure in sporocarp-associated beetle communities of the central Appalachian Mountains. Mycologia. 2010;102:785–802.

    PubMed  Article  Google Scholar 

  31. Hilario RR, Ferrari SF. Feeding ecology of a group of buffy-headed marmosets (Callithrix flaviceps): fungi as a preferred resource. Am J Primatol. 2010;72:515–21.

    PubMed  Google Scholar 

  32. Johnson CN. Interactions between mammals and ectomycorrhizal fungi. Trends Ecol Evol. 1996;11:503–7.

    CAS  PubMed  Article  Google Scholar 

  33. Kunz TH, Braun de Torrez E, Bauer D, Lobova T, Fleming TH. Ecosystem services provided by bats. Ann N Y Acad Sci. 2011;1223:1–38.

  34. Shilton LA, Altringham JD, Compton SG, Whittaker RJ. Old World fruit bats can be long-distance seed dispersers through extended retention of viable seeds in the gut. Proc R Soc B-Biol Sci. 1999;266:219–23.

    Article  Google Scholar 

  35. Dumont ER. Bats and fruit: an ecomorphological approach. In: Kunz TH, Fenton MB, editors. Bat Ecol. Chicago: The University of Chicago Press; 2003. p. 398–429.

    Google Scholar 

  36. Medellin RA, Gaona O. Seed dispersal by bats and birds in forest and disturbed habitats of Chiapas. Mexico Biotropica. 1999;31:478–85.

    Article  Google Scholar 

  37. Cardoso DaSilva JM, Uhl C, Murray G. Plant succession, landscape management, and the ecology of frugivorous birds in abandoned Amazonian pastures. Conserv Biol. 1996;10:491–503.

    Article  Google Scholar 

  38. Tiscornia S, Ruiz R, Bettucci L. Fungal endophytes from vegetative and reproductive tissues of Eugenia uruguayensis in Uruguay. Sydowia. 2012;64:313–28.

    Google Scholar 

  39. Martinson EO, Herre EA, Machado CA, Arnold AE. Culture-free survey reveals diverse and distinctive fungal communities associated with seveloping figs (Ficus spp.) in Panama. Microb Ecol. 2012;64:1073–84.

    PubMed  Article  Google Scholar 

  40. Schoch CL, Seifert KA, Huhndorf S, Robert V, Spouge JL, Levesque C, et al. Nuclear ribosomal internal transcribed spacer (ITS) region as a universal DNA barcode marker for Fungi. Proc Natl Acad Sci USA. 2012;109:1–6.

    Article  Google Scholar 

  41. Arnold AE, Henk DA, Eells RL, Lutzoni F, Vilgalys R. Diversity and phylogenetic affinities of foliar fungal endophytes in lobolly pine inferred by culturing and environmental PCR. Mycologia. 2007;99:185–206.

    CAS  PubMed  Article  Google Scholar 

  42. Martinson EO, Herre EA, Machado CA, Arnold AE. Culture-free survey reveals diverse and distinctive fungal communities associated with developing figs (Ficus spp.) in Panama. Microb Ecol. 2012;64:1073–84.

    PubMed  Article  Google Scholar 

  43. Unterseher M, Persoh D, Schnittler M. Leaf-inhabiting endophytic fungi of European Beech (Fagus sylvatica L.) co-occur in leaf litter but are rare on decaying wood of the same host. Fungal Divers. 2013;60:43–54.

    Article  Google Scholar 

  44. Hiergeist A, Gläsner J, Reischl U, Gessner A. Analyses of intestinal microbiota: culture versus sequencing. ILAR J. 2015;56:228–40.

    CAS  PubMed  Article  Google Scholar 

  45. Ward DM, Weller R, Bateson MM. 16S rRNA sequences reveal numerous uncultured microorganisms in a natural community. Lett Nat. 1990;345:183–7.

    Article  Google Scholar 

  46. Stefani FOP, Bell TH, Marchand C, de la Providencia IE, El Yassimi A, St-Arnaud M, et al. Culture-dependent and -independent methods capture different microbial community fractions in hydrocarbon-contaminated soils. Hu S, editor. PLoS ONE. 2015;10:e0128272.

  47. Bertolino S, Vizzini A, Wauters LA, Tosi G. Consumption of hypogeous and epigeous fungi by the red squirrel (Sciurus vulgaris) in subalpine conifer forests. For Ecol Manag. 2004;202:227–33.

    Article  Google Scholar 

  48. Spain AM, Elshahed MS, Najar FZ, Krumholz LR. Metatranscriptomic analysis of a high-sulfide aquatic spring reveals insights into sulfur cycling and unexpected aerobic metabolism. PeerJ. 2015;3:e1259.

    PubMed  PubMed Central  Article  CAS  Google Scholar 

  49. Zhao M, Zhang D, Su X, Duan S, Wan J, Yuan W, et al. An Integrated Metagenomics/ Metaproteomics investigation of the microbial communities and enzymes in solid-state fermentation of Pu-erh tea. Sci Rep. 2015;5.

  50. Kirschner R, Hsu T, Tuan NN, Chen C-L, Huang S-L. Characterization of fungal and bacterial components in gut/fecal microbiome. Curr Drug Metab. 2015;16:272–83.

    CAS  PubMed  Article  Google Scholar 

  51. Sørensen J, Nicolaisen MH, Ron E, Simonet P. Molecular tools in rhizosphere microbiology-from single-cell to whole-community analysis. Plant Soil. 2009;321:483–512.

    Article  CAS  Google Scholar 

  52. Dostál P. Post-dispersal seed mortality of exotic and native species: Effects of fungal pathogens and seed predators. Basic Appl Ecol. 2010;11:676–84.

    Article  Google Scholar 

  53. Connell JH. On the role of natural enemies in preventing competitive exclusion in some marine animals and in rain forest trees. In: Den BPJ, Gradwell GR, editors. Dyn Popul. Wageningen: Pudoc; 1971. p. 298–312.

    Google Scholar 

  54. Janzen DH. Seed predation by animals. Annu Rev Ecol Syst. 1971;2:465–92.

    Article  Google Scholar 

  55. Gillett JB. Pest pressure, an underestimated factor in evolution. Syst Assoc Publ Number. 1962;4:37–46.

    Google Scholar 

  56. Morrison DW. Foraging and day-roosting dynamics of canopy fruit bats in Panama. J Mammal. 1980;61:20–9.

    Article  Google Scholar 

  57. Gilbert GS. Dimensions of plant disease in tropical forests. In: Burslem D, Pinard M, Hartley S, editors. Biot Interact Trop their Role Maint Species Divers. Cambridge: Cambridge; 2005. p. 141–64.

    Chapter  Google Scholar 

  58. Haelewaters D, Gorczak M, Kaishian P, De Kesel A, Blackwell M. Laboulbeniomycetes, enigmatic fungi with a turbulent taxonomic history. In: Zaragoza Ó, Casadevall ABT-E of M (eds). Encycl Mycol Vol 1. Oxford: Elsevier; 2021. p. 263–83.

  59. Blackwell M. Minute mycological mysteries: the influence of arthropods on the lives of fungi. Mycologia. 1994;86:1–17.

    Article  Google Scholar 

  60. Machado CA, Robbins N, Gilbert MTP, Herre EA. Critical review of host specificity and its coevolutionary implications in the fig/fig-wasp mutualism. Proc Natl Acad Sci USA. 2005;102:6558–65.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  61. Cappelli A, Ulissi U, Valzano M, Damiani C, Epis S, Gabrielli MG, et al. A Wickerhamomyces anomalus killer strain in the malaria vector Anopheles stephensi. PLoS ONE. 2014;9:95988.

    Article  CAS  Google Scholar 

  62. Malassigné S, Minard G, Vallon L, Martin E, Valiente Moro C, Luis P. Diversity and functions of yeast communities associated with insects. Microorg. 2021;9:1552.

    Article  CAS  Google Scholar 

  63. Kunčič MK, Kogej T, Drobne D, Gunde-Cimerman N. Morphological response of the halophilic fungal genus Wallemia to high salinity. Appl Environ Microbiol. 2010;76:329–37.

    Article  CAS  Google Scholar 

  64. Amend A. From dandruff to deep-sea vents: Malassezia-like fungi are ecologically hyper-diverse. PLoS Pathog. 2014;10:e1004277.

  65. Hallen-Adams HE, Suhr MJ. Fungi in the healthy human gastrointestinal tract. Virulence. 2017;8:352–8.

    CAS  PubMed  Article  Google Scholar 

  66. Lafond M, Bouza B, Eyrichine S, Rouffineau F, Saulnier L, Giardina T, et al. In vitro gastrointestinal digestion study of two wheat cultivars and evaluation of xylanase supplementation. J Anim Sci Biotechnol. 2015;6:5.

    PubMed  PubMed Central  Article  CAS  Google Scholar 

  67. Geib SM, Filley TR, Hatcher PG, Hoover K, Carlson JE, Jimenez-Gasco M del M, et al. Lignin degradation in wood-feeding insects. Proc Natl Acad Sci U S A. 2008;105:12932–7.

  68. Yurkov A, Krüger D, Begerow D, Arnold N, Tarkka MT. Basidiomycetous yeasts from boletales fruiting bodies and their interactions with the mycoparasite Sepedonium chrysospermum and the host fungus Paxillus. Microb Ecol. 2012;63:295–303.

    PubMed  Article  Google Scholar 

  69. Petrik S, Marova I, Haronikova A, Kostovova I, Breierova E. Production of biomass, carotenoid and other lipid metabolites by several red yeast strains cultivated on waste glycerol from biofuel production—a comparative screening study. Ann Microbiol. 2013;63:1537–51.

    CAS  Article  Google Scholar 

  70. Yurkov AM, Vustin MM, Tyaglov BV, Maksimova IA, Sineokiy SP. Pigmented basidiomycetous yeasts are a promising source of carotenoids and ubiquinone Q10. Microbiology. 2008;77:1–6.

    CAS  Article  Google Scholar 

  71. de Melo PGV, Beux M, Pagnoncelli MGB, Soccol VT, Rodrigues C, Soccol CR. Isolation, selection and evaluation of antagonistic yeasts and lactic acid bacteria against ochratoxigenic fungus Aspergillus westerdijkiae on coffee beans. Lett Appl Microbiol. 2016;62:96–101.

    Article  CAS  Google Scholar 

  72. Venugopalan A, Srivastava S. Endophytes as in vitro production platforms of high value plant secondary metabolites. Biotechnol Adv. 2015;33:873–87.

    PubMed  Article  Google Scholar 

  73. Kandasamy D, Gershenzon J, Hammerbacher A. Volatile organic compounds emitted by fungal associates of conifer bark beetles and their potential in bark beetle control. J Chem Ecol. 2016;42:952–69.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  74. Whitehead SR, Poveda K. Herbivore-induced changes in fruit-frugivore interactions. J Ecol. 2011;99:964–9.

    Article  Google Scholar 

  75. Bennett JW, Hung R, Lee S, Padhi S. 18 Fungal and bacterial volatile organic compounds: an overview and their role as ecological signaling agents Bt - fungal Associations. In: Hock B, editor. Mycota IX Fungal Assoc. Berlin, Heidelberg: Springer Berlin Heidelberg; 2012. p. 373–93.

  76. Bloss J, Acree TE, Bloss JM, Hood WR, Kunz TH. Potential use of chemical cues for colony-mate recognition in the big brown bat. Eptesicus fuscus J Chem Ecol. 2002;28:819–34.

    CAS  PubMed  Article  Google Scholar 

  77. Bartelt RJ, Wicklow DT. Volatiles from Fusarium verticillioides (Sacc.) Nirenb. and their attractiveness to nitidulid beetles. J Agric Food Chem. 1999;47:2447–54.

  78. Rangel LI, Hamilton O, de Jonge R, Bolton MD. Fungal social influencers: secondary metabolites as a platform for shaping the plant-associated community. Plant J. 2021;108:632–45.

    CAS  PubMed  Article  Google Scholar 

  79. Colgan W, Claridge AW. Mycorrhizal effectiveness of Rhizopogon spores recovered from faecal pellets of small forest-dwelling mammals. Mycol Res. 2002;106:314–20.

    Article  Google Scholar 

  80. Charalambidou I, Santamaria L, Langevoord O. Effect of ingestion by five avian dispersers on the retention time, retrieval and germination of Ruppia maritima seeds. Funct Ecol. 2003;17:747–53.

    Article  Google Scholar 

  81. Vega FE, Simpkins A, Aime MC, Posada F, Peterson SW, Rehner SA, et al. Fungal endophyte diversity in coffee plants from Colombia, Hawai’i Mexico and Puerto Rico. Fungal Ecol. 2010;3:122–38.

    Article  Google Scholar 

  82. Paul NC, Lee HB, Lee JH, Shin KS, Ryu TH, Kwon HR, et al. Endophytic fungi from Lycium chinense Mill and characterization of two new korean records of Colletotrichum. Int J Mol Sci. 2014;15:15272–86.

    PubMed  PubMed Central  Article  CAS  Google Scholar 

  83. Taylor MW, Tsai P, Anfang N, Ross HA, Goddard MR. Pyrosequencing reveals regional differences in fruit-associated fungal communities. Environ Microbiol. 2014;16:2848–58.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  84. Brooke AP. Tent selection, roosting ecology and social organization of the tent-making bat, Ectophylla alba, in Costa Rica. J Zool. 1990;221:11–9.

    Article  Google Scholar 

  85. Reid FA. A field guide to the mammals of Central America and southeast Mexico. New York: Oxford University Press; 1997.

    Google Scholar 

  86. Rodríguez-Herrera B, Pineda W. The IUCN Red List of Threatened Species 2015. 2015.

  87. Rodríguez-Herrera B, Medellín RA, Gamba-Rios M. Roosting requirements of white tent-making bat Ectophylla alba (Chiroptera: Phyllostomidae). Acta Chiropterologica. 2008;10:89–95.

    Article  Google Scholar 

  88. De la Llata Quiroga E, Ruedas LA, Mora JM. A comparison of fruit removal in Ficus colubrinae between birds and Ectophylla alba (Chiroptera: Phyllostomidae) in a Costa Rican rain forest. Stud Neotrop Fauna Environ. 2021;1–8.

  89. Gazis R, Chaverri P. Diversity of fungal endophytes in leaves and stems of wild rubber trees (Hevea brasiliensis) in Peru. Fungal Ecol. 2010;3:240–54.

    Article  Google Scholar 

  90. White TJ, Bruns T, Lee S, Taylor J. Amplification and direct sequencing of fungal ribosomal RNA genes for phylogenetics. In: White TJ, editor. PCR Protoc A Guid to Methods Appl. San Diego: Academic Press; 1990. p. 315–22.

    Google Scholar 

  91. Ihrmark K, Bodeker IT, Cruz-Martinez K, Friberg H, Kubartova A, Schenck J, et al. New primers to amplify the fungal ITS2 region–evaluation by 454-sequencing of artificial and natural communities. FEMS Microbiol Ecol. 2012;82:666–77.

    CAS  PubMed  Article  Google Scholar 

  92. Martin M. Cutadapt removes adapter sequences from high-throughput sequencing reads. EMBnet.journal. 2011;17:10–2.

  93. Callahan BJ, McMurdie P, Rosen M, Han AW, Johnson AJA, Holmes SP. DADA2: High-resolution sample inference from Illumina amplicon data. Nat Methods. 2016;13:581–3.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  94. Callahan BJ, McMurdie PJ, Holmes SP. Exact sequence variants should replace operational taxonomic units in marker-gene data analysis. ISME J. 2017;11:2639–43.

    PubMed  PubMed Central  Article  Google Scholar 

  95. Kõljalg U, Nilsson RH, Abarenkov K, Tedersoo L, Taylor AFS, Bahram M, et al. Towards a unified paradigm for sequence-based identification of fungi. Mol. Ecol. 2013. p. 5271–7.

  96. Wright ES. Using DECIPHER v2.0 to analyze big biological sequence data in R. R J. 2016;8:352–359.

  97. Liu C, Cui Y, Li X, Yao M. microeco : an R package for data mining in microbial community ecology. FEMS Microbiol Ecol. 2021;97.

  98. Lozupone CA, Hamady M, Kelley ST, Knight R. Quantitative and qualitative beta diversity measures lead to different insights into factors that structure microbial communities. Appl Environ Microbiol. 2007;73:1576–85.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  99. McKnight DT, Huerlimann R, Bower DS, Schwarzkopf L, Alford RA, Zenger KR. Methods for normalizing microbiome data: an ecological perspective. Methods Ecol Evol. 2019;10:389–400.

    Article  Google Scholar 

  100. Ramette A. Multivariate analyses in microbial ecology. FEMS Microbiol Ecol. 2007;62:142–60.

    CAS  PubMed  Article  Google Scholar 

  101. Anderson MJ. A new method for non-parametric multivariate analysis of variance. Austral Ecol. 2001;26:32–46.

    Google Scholar 

  102. Segata N, Izard J, Waldron L, Gevers D, Miropolsky L, Garrett WS, et al. Metagenomic biomarker discovery and explanation. Genome Biol. 2011;12:R60.

    PubMed  PubMed Central  Article  Google Scholar 

  103. Royle JA, Nichols JD. Estimating abundance from repeated presence–absence data or point counts. Ecology. 2003;84:777–90.

    Article  Google Scholar 

  104. Porter TM, Hajibabaei M. Scaling up: A guide to high-throughput genomic approaches for biodiversity analysis. Mol Ecol. 2018;27:313–38.

    PubMed  Article  Google Scholar 

  105. Nguyen NH, Smith D, Peay K, Kennedy P. Parsing ecological signal from noise in next generation amplicon sequencing. New Phytol. 2015;205:1389–93.

    CAS  PubMed  Article  Google Scholar 

  106. McLaren MR, Willis AD, Callahan BJ. Consistent and correctable bias in metagenomic sequencing experiments. Turnbaugh P, Garrett WS, Turnbaugh P, Quince C, Gibbons S, editors. Elife. 2019;8:e46923.

  107. Hu Y-J, Lane A, Satten GA. A rarefaction-based extension of the LDM for testing presence–absence associations in the microbiome. Bioinformatics. 2021;37:1652–7.

    CAS  PubMed Central  Article  Google Scholar 

  108. McMurdie PJ, Holmes S. Waste not, want not: why rarefying microbiome data is inadmissible. PLoS Comput Biol. Public Library of Science; 2014;10:e1003531.

  109. Nguyen NH, Song Z, Bates ST, Branco S, Tedersoo L, Menke J, et al. FUNGuild: An open annotation tool for parsing fungal community datasets by ecological guild. Fungal Ecol. 2016;20:241–8.

    Article  Google Scholar 

  110. Gazis R, Chaverri P. Wild trees in the Amazon basin harbor a great diversity of beneficial endosymbiotic fungi: Is this evidence of protective mutualism? Fungal Ecol. 2015;17:18–29.

    Article  Google Scholar 

  111. Rodriguez RJ, Jr. JFW, Arnold AE, Redman RS, White Jr JF, Arnold AE, et al. Fungal endophytes: Diversity and functional roles: Tansley review. New Phytol. 2009;182:314–30.

  112. Parfitt D, Hunt J, Dockrell D, Rogers HJ, Boddy L. Do all trees carry the seeds of their own destruction? PCR reveals numerous wood decay fungi latently present in sapwood of a wide range of angiosperm trees. Fungal Ecol. 2010;3:338–46.

    Article  Google Scholar 

Download references


We deeply acknowledge J.P. Barrantes, C. Castillo-Salazar, and E. Hellman for their help in laboratory and fieldwork. J. Aihartza, I. Garin, L. Jiménes, E. Paniagua, and E. Rojas also provided valuable field assistance. We thank the staff at La Selva Biological Station for their help with logistics and accommodation, and Lourdes Vargas from SINAC for her help with research permits. We would also like to thank Silvia Chaves Ramírez for providing the bat and fruit illustrations used in Figs. 2, 3, 4.


This research was supported by a Grant from Conservation, Food and Health Foundation.

Author information




PC and GC contributed equally to conceiving, collecting and analyzing data, and writing the paper. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Priscila Chaverri.

Ethics declarations

Ethics approval and consent to participate

Not applicable.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher's Note

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

Supplementary Information

Additional file 1

. Supplementary Information (Tables and Figures).

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Chaverri, P., Chaverri, G. Fungal communities in feces of the frugivorous bat Ectophylla alba and its highly specialized Ficus colubrinae diet. anim microbiome 4, 24 (2022).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI:


  • Dispersal
  • Ecological network
  • Endophytes
  • Epiphytes
  • Gut microbiome
  • Janzen–Connell hypothesis
  • Metabarcoding
  • Mycobiota
  • Pyxidiophora
  • Theory of pest pressure