Introduction
Bat flies (Diptera: Hippoboscoidea) are obligate, blood-feeding ectoparasites that exclusively parasitize bats and are among the most abundant and frequent hematophagous parasites in this mammalian group (Hrycyna et al., Reference Hrycyna, Martins and Graciolli2019). These flies typically exhibit a high degree of host specificity, often associated with long-term evolutionary relationships between parasite and host lineages (Dick et al., Reference Dick, Miller, Brown, Borkent, Cumming, Wood and Woodley2010). However, host specificity can also be influenced by intrinsic factors such as host behaviour, health and body size, as well as extrinsic environmental factors (Palheta et al., Reference Palheta, Urbieta, Brasil, Dias-Silva, Da Silva, Graciolli, Aguiar and Vieira2020). Studying these tightly linked parasite–host interactions offers valuable insight into co-evolutionary processes and the mechanisms shaping host specialization (Dick and Patterson, Reference Dick, Patterson, Morand, Krasnov and Poulin2006; Hiller et al., Reference Hiller, Vollstädt, Brändel, Page and Tschapka2021). Globally, bat flies are divided into two main families with distinct biogeographic distributions: Nycteribiidae, more diverse in the Eastern Hemisphere, and Streblidae, predominantly found in the Western Hemisphere, especially in the Neotropics (Soares et al., Reference Soares, Graciolli, Alcântara, Ribeiro, Valença and Ferrari2013; Graciolli and Dick, Reference Graciolli and Dick2018; Barbier et al., Reference Barbier, Nunes, Rocha, Rocha and Cordeiro-Estrela2019; Graciolli et al., Reference Graciolli, Guerrero and Catzeflis2019). In South America, representatives of both families coexist, parasitizing a wide range of bat species (Biz et al., Reference Biz, Bastazini, Carvalho and Ramos-Pereira2023; Zapata-Mesa et al., Reference Zapata-Mesa, Montoya-Bustamante, Hoyos, Peña, Galindo-González, Chacón-Pacheco, Ballesteros-Correa, Graciolli, Nogueira and Mello2024). Despite their ecological relevance, detailed studies on bat–fly associations remain scarce in many parts of the Neotropics, including Peru.
Peru harbours remarkable bat diversity, with 196 species currently recorded (Pacheco et al., Reference Pacheco, Diaz, Graham-Angeles, Flores-Quispe, Calizaya-Mamani, Ruelas and Sánchez-Vendizú2021; Velazco, Reference Velazco2023; Diaz et al., Reference Díaz, Medina, Arias, López and Carrión2025), of which at least 75 are known to host ectoparasites, including 66 species of bat flies (Minaya et al., Reference Minaya, Silva and Iannacone2021). However, most records of bat–fly diversity in Peru are concentrated in the lowland Amazonian rainforests of Loreto and Madre de Dios (Theodor, Reference Theodor1967; Guerrero, Reference Guerrero, Wilson and Sandoval1996a; Graciolli et al., Reference Graciolli, Autino and Claps2007; Autino et al., Reference Autino, Claps, Barquez and Díaz2011; Gettinger, Reference Gettinger2018; Gettinger et al., Reference Gettinger, Epperson, Hermasillo and Gardner2020; Morales-Malacara and Guerrero, Reference Morales-Malacara and Guerrero2020). In contrast, montane forests – despite being among the most bat–diverse ecosystems in the Neotropics (Chaverri et al., Reference Chaverri, Garin, Alberdi, Jimenez, Castillo-Salazar and Aihartza2016; Bogoni et al., Reference Bogoni, Peres and Kmpmb2021) – remain poorly studied in terms of bat ectoparasite associations. Available information from these ecosystems is sparse and typically limited to isolated records of ectoparasite presence, rather than comprehensive analyses of interactions (Biz et al., Reference Biz, Bastazini, Carvalho and Ramos-Pereira2023; Zapata-Mesa et al., Reference Zapata-Mesa, Montoya-Bustamante, Hoyos, Peña, Galindo-González, Chacón-Pacheco, Ballesteros-Correa, Graciolli, Nogueira and Mello2024). For example, in the montane forests of Amazonas, northern Peru, a key region within the Andean forest belt, only a single study on bat flies has been published (Ibáñez and Jara, Reference Ibáñez and Jara2008), highlighting a substantial gap in our understanding of host–parasite relationships in these high-elevation systems.
The use of ecological network analysis has become increasingly important in advancing our understanding of host–parasite systems, particularly in bat–ectoparasite relationships (Runghen et al., Reference Runghen, Poulin, Monlleó-Borrull and Llopis-Belenguer2021; Biz et al., Reference Biz, Bastazini, Carvalho and Ramos-Pereira2023; Zapata-Mesa et al., Reference Zapata-Mesa, Montoya-Bustamante, Hoyos, Peña, Galindo-González, Chacón-Pacheco, Ballesteros-Correa, Graciolli, Nogueira and Mello2024). Unlike traditional species inventories, network approaches allow researchers to explore structural properties such as modularity, nestedness and interaction specialization within ecological communities (Bezerra and Bocchiglieri, Reference Bezerra and Bocchiglieri2022). Bats and their ectoparasitic flies represent an ideal model for such studies due to their high species richness and long coevolutionary history. Network-based analyses have consistently revealed high levels of specialization and modularity in bat–fly associations, suggesting that both ecological and evolutionary factors shape these interaction patterns (Falcão et al., Reference Falcão, Araújo, Leite, Fagundes, Espírito-Santo, Zazá-Borges, Vasconcelos, Fernandes and Paglia2022). Moreover, understanding the structure of host–parasite networks is increasingly relevant for public health, as ectoparasites may act as vectors of zoonotic pathogens, potentially facilitating transmission between wildlife and humans (Szentiványi et al., Reference Szentiványi, Christe and Glaizot2019). Integrating network analysis into parasite–host studies therefore offers valuable insights for both ecological theory and applied conservation and health strategies.
Given this context and the need to generate local-scale data that contribute to a broader understanding of bat–ectoparasite interactions, the aim of our study was to assess the dipteran ectoparasites associated with bats in the montane forests of Amazonas, northern Peru. Specifically, we sought to examine patterns of species distribution to explore the structure of parasite–host association using interaction networks.
Materials and methods
Study area
The study was conducted in the hamlet of Nueva Esperanza on the Numparket Waterfall Tourist Route located in Aramango, Bagua (Amazonas, Peru). The area is characterized by low montane forest vegetation and falls within the ecoregions of Very Humid Montane Forest (Bosque Montano Muy Húmedo, BMHM) and Very Humid Premontane Forest (Bosque Premontano Muy Húmedo, BMHP) (Britto, Reference Britto2017). The landscape includes both well-preserved primary forest and zones subjected to selective logging.
Fieldwork was conducted between July 2023, February 2024 and August 2024 across three sites: Numparket (1800 m a.s.l.), Chontas (1560 m a.s.l.) and Higuerón (1480 m a.s.l.). Numparket is located around the Numparket Waterfall within the conservation concession Cerro El Adobe. This area also forms part of the buffer zone of both the national sanctuary Cordillera de Colán and the communal reserve Chayu Nain (Figure 1). Due to its proximity to the waterfall and its tributary rivers, Numparket maintains high humidity throughout most of the year. The site is predominantly covered by well-preserved primary forest, with the exception of some disturbed zones near the road. In contrast, Chontas and Higuerón are located within the area of influence of Cordillera de Colán. These sites are characterized by a mosaic of preserved forest, patches of secondary growth and areas affected by selective logging. The vegetation includes species such as Ficus paraensis, Cecropia spp., various species of Araceae and Rubiaceae, and abundant pteridophytes (Authors’ observation).
Geographic location of the sampling sites in the Amazonas region, northern Peru, where bats were captured and their ectoparasitic flies collected during the 2023–2024 field campaign.
Bats and flies sampling
At every sampling station, 10 understory mist nets (12 × 2.5 m) spaced ∼20 m apart were used for 12 nights, distributed as three blocks of four consecutive nights (July 2023, February 2024 and August 2024) (MINAM, 2015). The nets were opened from 18:00 to 00:00 h to target species with peak foraging activity during that period (Jones et al., Reference Jones, McShea, Conroy, Kunz, Wilson, Cole, Nichols, Rudran and Foster1996). Individuals that could not be identified in the field were collected, preserved in alcohol and deposited at the Museo Vera Alleman de la Universidad Ricardo Palma (MURP) in Lima, Peru.
All individuals were checked while alive. Bat flies were removed using entomological forceps, then fixed and preserved in polypropylene cryovials containing 70% ethanol. Specimens were cleared in 10% KOH and examined under a Nikon SMZ745 stereomicroscope for taxonomic identification using the keys of Wenzel et al. (Reference Wenzel, Tipton, Kiewlicz, Wenzel and Tipton1966); Wenzel (Reference Wenzel1976) and Guerrero (Reference Guerrero1993; Reference Guerrero1994a, Reference Guerrero1994b, Reference Guerrero1995a, Reference Guerrero1995b, Reference Guerrero1996b, Reference Guerrero1996c). Macrophotographs of external anatomy and taxonomically important structures were taken using a TOUPCAM camera mounted on a Nikon Eclipse Si microscope with a Nikon Nii LED illumination system. Image stacking was performed with ToupView software. All ectoparasite specimens were deposited in the entomological collection of the Natural History Museum, Faculty of Natural Sciences and Mathematics, Universidad Nacional Federico Villarreal (MUFV).
Specialization of bat-fly interactions
Host specificity of bat flies was classified as follows: monoxenous (ectoparasitic flies utilizing only a single host species), oligoxenous (utilizing two or more congeneric species), pleioxenous (utilizing two or more host genera within the same family) and polyxenous (utilizing multiple hosts from different families) (Marshall, Reference Marshall and Kunz1982; Seneviratne et al., Reference Seneviratne, Fernando and Udagama-Randeniya2009). The parasite population component was analysed using standard ecological parasitological indices: prevalence (P%) and mean intensity of infection (MI), following Bush et al. (Reference Bush, Lafferty, Lotz and Shostak1997) and Bautista-Hernández et al. (Reference Bautista-Hernández, Monks, Pulido-Flores and Rodríguez-Ibarra2015).
Using the bat–fly interaction encounters, we constructed weighted bipartite networks for each of the three sites as well as an aggregated regional network. To assess the specificity of bat–fly interactions in these networks, we applied modularity and specialization metrics. We first tested for modularity (Qw) using the weighted DIRTLPAwb + algorithm (Beckett, Reference Beckett2016) and then evaluated low-level nestedness (within modules) using the WNODAsm metric (Pinheiro et al., Reference Pinheiro, Felix, Dormann and Mello2019). Following Pinheiro et al. (Reference Pinheiro, Felix and Lewinsohn2022), we did not test for nestedness in the overall network, as all networks were significantly modular (see Results). Modularity measures the extent to which species and their interactions can be divided into subgroups (modules) that are more interconnected within themselves than with others (Newman, Reference Newman2006). Nestedness reflects the pattern in which interactions of species with fewer connections (specialists) form subsets of the interactions of species with more connections (generalists) (Mariani et al., Reference Mariani, Ren, Bascompte and Tessone2019). Additionally, specialization was quantified using the H2’ metric (complementary specialization), which captures how selective the network is beyond what would be expected based on species relative abundances, approximated by the matrix’s marginal totals (Blüthgen et al., Reference Blüthgen, Menzel and Blüthgen2006). To test the significance of these metrics, we used the equiprobable (preserving species richness and total number of interactions) and proportional (same as equiprobable but also preserving marginal sums) null models described in Pinheiro et al. (Reference Pinheiro, Felix and Lewinsohn2022). Specifically, the restricted version of the equiprobable null model, which also maintains the modular structure during randomizations, was used to test WNODAsm, while the proportional null model was employed for Q w and H2’.
Species-level metrics were used to explore variation in the specialization of bat and fly species across sites, focusing only on species present at all sites. Species-level specialization (d’) was employed to describe how selectively a bat or fly interacts with available species from the opposite group within the network, based on the frequency of their interactions (Blüthgen et al., Reference Blüthgen, Menzel and Blüthgen2006). For flies, higher d′ values indicate higher host specificity. For bats, which do not choose their parasites, higher d′ values indicate that their assemblage of flies is composed of more host–specific parasites, whereas lower values indicate association with more generalist parasites. Additionally, species strength was calculated as the total sum of interaction proportions across all partners for a given species, reflecting how dependent bats or flies are on that species (Bascompte et al., Reference Bascompte, Jordano and Olesen2006). These species-level metrics, along with the network-level metrics mentioned above, were calculated using the package ‘bipartite’ (Dormann et al., Reference Dormann, Gruber and Fruend2008) in the software R 4.4.1.
Sampling coverage of networks was also analysed following the suggestion of Chiu et al. (Reference Chiu, Chao, Vogel, Kriegel and Thorn2023). This metric indicates the proportion of the total number of interaction events represented by the detected interactions. For this assessment, the ‘iNext.link’ package (Hsieh et al., Reference Hsieh, Ma and Chao2016) was used in the software R 4.4.1.
Finally, to quantify differences between interaction networks among sites, we followed the approach of Fründ (Reference Fründ2021), which decomposes total link dissimilarity into additive components. For each pair of sites, we calculated: βWN, the overall dissimilarity between the two interaction networks; βOS, the dissimilarity attributable to changes in interactions among species shared between sites (rewiring); βST, the dissimilarity attributable to species turnover, i.e. interactions that differ because one or both interacting species are present at only one site. Following Novotny (Reference Novotny2009), βST was further partitioned into turnover caused exclusively by flies (βST.f), exclusively by bats (βST.b), or jointly by both (βST.fb). This analysis was performed using the betalinkr_multi function of the package ‘bipartite’ (Dormann et al., Reference Dormann, Gruber and Fruend2008) in the software R 4.4.1., specifically with the ‘commondenom’ partition method (Fründ, Reference Fründ2021).
Results
Bats and flies
A total of 160 bats were captured, including 152 individuals from the family Phyllostomidae and 8 from Vespertilionidae. The bats belonged to 23 species, of which 71 individuals from 14 species were parasitized by bat flies (Table 1). The most abundant bat species were Carollia brevicauda (n = 48), C. perspicillata (n = 36) and Sturnira oporaphilum (n = 19). The species with the highest number of parasitized individuals were C. brevicauda (n = 30) and C. perspicillata (n = 18). Among the three sampling areas, Numparket had the highest number of captured bats (n = 73) and parasitized individuals (n = 33) (Table 1). C. brevicauda presented the highest abundance of ectoparasites across all three sites. In Numparket and Higuerón, C. perspicillata ranked second in parasite abundance, while in Chontas, the second most parasitized species was Myotis nigricans.
Bats captured along the Nueva Esperanza Trail to Numparket Falls, Amazonas, Peru, and specific characterization based on their ectoparasitic flies
Bf, bat fly abundance; E(P), no. of examined bats (parasitized); MI, mean intensity; P%, prevalence; S, bat fly species richness.
A total of 155 ectoparasitic flies were collected, representing 19 species from the families Streblidae (17 species) and Nycteribiidae (2 species) (Figures 2 –3). The most abundant fly species were Paraeuctenodes similis Wenzel, 1976 (n = 44), Trichobius joblingi Wenzel, 1966 (n = 43) and Megistopoda proxima (Séguy, 1926) (n = 13). Numparket and Higuerón exhibited the highest bat fly species richness (s = 11) and abundance (n = 56), followed by Chontas (s = 8, n = 43). Only P. similis and T. joblingi were recorded in all three areas. In terms of host specificity, most bat fly species were classified as monoxenous (n = 10), followed by oligoxenous (n = 8), and pleioxenous (n = 1) (Table 2).
Species of ectoparasitic diptera from bats captured along the Nueva Esperanza Trail to Numparket Falls, Amazonas, Peru (first part). (A) Anastrebla caudiferae, (B) Aspidoptera falcata, (C) Exastinion oculatum, (D) Megistopoda proxima, (E) Metelasmus pseudopterus, (F) Neotrichobius bisetosus, (G) Paraeuctenodes similis, (H) Paratrichobius longicrus complex, (I) Strebla guajiro.
Species of ectoparasitic diptera from bats captured along the Nueva Esperanza Trail to Numparket Falls, Amazonas, Peru (second part). (A) Anatrichobius scorzai, (B) Basilia anceps, (C) Anatrichobius sp., (D) Speiseria ambigua, (E) Trichobius joblingi, (F) Paratrichobius salvini complex.
Associations and characterization of ectoparasitic flies from bat captured along the Nueva Esperanza Trail to Numparket Falls, Amazonas, Peru
C, Chontas (1560 m); H, Higuerón (1480 m); N, Numparket (1800 m); P%, Prevalence; Spe, Specificity (Mon, Monoxenous; Oli, Oligoxenous; Ple, Pleioxenous).
** New species record for Peru; *New record only for Amazonas.
Specialization of bat–fly interactions
Sampling coverage of bat–fly networks was always above 0.85 (Table 3), indicating that they are a good representation of bat–fly interactions at each site as well as at the regional scale. All networks exhibited a modular topology, but nestedness within modules was observed only in the regional network (Table 3). Specificity of interactions was intermediate to high across all networks, as indicated by Qw (≥0.48) and H2’ (≥0.75). The highest values for these metrics were observed in Chontas (middle elevation), while Higuerón and Numparket showed similar values, still reflecting high specificity. The modular structure showed a clear partition based on the phylogenetic relationships of bats, with modules never including unrelated taxa of bats (Figure 4). Anoura, Carollia, Myotis and Sturnira species were always grouped within the same module, except for Carollia at higher altitudes (Numparket). In this latter case, C. brevicauda and C. perspicillata formed their own module, although they still shared the same flies.
Modular structure of bat–fly interaction networks along the Nueva Esperanza Trail to Numparket Falls, Amazonas, Peru. The regional network is the result of the aggregation of interactions of the three other local networks. Interactions and species of the same module share specific colours and interactions between species of different modules are in grey.
Sampling coverage and structural properties of bat–fly interaction networks along the Nueva Esperanza Trail to Numparket Falls, Amazonas, Peru. The regional network is the result of the aggregation of interactions of the three local networks
SC, sampling coverage; Qw, weighted modularity; WNODAsm, within-module nestedness, H2’: complementary specialization. Statistically significant values (P < 0.05) are in bold.
Species showed different patterns across sites in terms of species-level specialization and species strength (Figure 5, Supplementary table S1-S2). Carollia species exhibited higher specialization at middle elevations (Chontas), with C. brevicauda consistently showing higher values than C. perspicillata. Myotis riparius displayed the highest specialization values at lower (Higuerón) and higher (Numparket) elevations. However, species strength for C. brevicauda decreased continuously from lower to higher elevations, while values for C. perspicillata increased at higher elevations, eventually surpassing C. brevicauda. Myotis riparius showed the same pattern observed in specialization, with lower species strength values at middle elevations. Among flies, T. joblingi was more specialized than P. similis at lower and middle elevations, but roles reversed at higher elevations. However, T. joblingi had higher species strength values than P. similis at the lower-elevation site, while the opposite occurred at middle- and higher-elevation sites.
Species-level metrics of bat–fly interaction networks along the Nueva Esperanza Trail to Numparket Falls, Amazonas, Peru. Species-level specialization d’ (A and B) and species strength (C and D). Only species present in the three evaluated sites are assessed.
Network dissimilarity between elevations was moderate to high, with βWN values ranging from 0.60 (Higuerón–Chontas) to 0.70 (Chontas–Numparket) (Table 4). Across all pairwise comparisons, the contribution of rewiring among shared species (βOS) was consistently low (0.03–0.04), representing only 5–7% of total link dissimilarity. In contrast, species turnover (βST) accounted for the vast majority of network differences (93–95%). Within the turnover component, co-turnover of both bats and flies (βST.fb) was consistently the dominant factor, contributing approximately 50–60% of βST in all comparisons. Turnover restricted to flies (βST.f) also made a substantial contribution (33–40%), whereas turnover restricted to bats (βST.b) was smaller (0–21%) and, in one comparison, absent.
Pairwise network dissimilarity between sites. ΒOS represents the dissimilarity attributable to rewiring among species shared between sites. ΒWN is the overall dissimilarity between interaction networks. ΒST corresponds to the component of dissimilarity explained by species turnover. ΒST.F, βST.B and βST.Fb indicate the portions of βST attributable to turnover restricted to fly species, restricted to bat species, and jointly to both trophic levels, respectively
Discussion
Bats and flies
This study expands the current knowledge of bat–ectoparasite interactions in the montane forests of northern Peru by documenting 19 species of parasitic flies associated with 14 bat species. All bat fly species reported represent new records for Amazonas (Minaya et al., Reference Minaya, Silva and Iannacone2021), and Anastrebla caudiferae, Exastinion deceptivum, Exastinion oculatum, P. similis and the Paratrichobius salvini complex are documented for the first time in Peru. An important observation concerns the record of T. joblingi. Although this species was previously reported from Condorcanqui, Amazonas, by Ibáñez and Jara (Reference Ibáñez and Jara2008), we noted inconsistencies between their figure and the diagnostic characters of T. joblingi. Based on our specimens, we provide the first confirmed record of T. joblingi for the department of Amazonas.
Two bat fly species, P. similis and T. joblingi, represented the dominant and most abundant core ectoparasites in the bats sampled in Amazonas. Both species parasitized more than 50% of the C. brevicauda and C. perspicillata individuals examined. A similar pattern was observed in the Magdalena River basin (López-Rivera et al., Reference López-Rivera, Robayo-Sánchez, Ramírez-Hernández, Cuéllar-Saénz, Villar, Cortés-Vecino, Rivera-Páez, Ossa-López, Ospina-Pérez, Henao-Osorio, Cardona-Giraldo, Racero-Casarrubia, Rodríguez-Posada, Morales-Martinez, Hidalgo and Ramírez-Chaves2024) and in Caldas (Raigosa et al., Reference Raigosa, García, Autino and Gomes2020), both in Colombia, where approximately 50% of individuals were primarily parasitized by T. joblingi and P. similis, consistent with our observations in Amazonas. These findings suggest that both dipteran species exhibit strong host specificity toward Carollia bats, maintaining a stable host–parasite association across distinct Neotropical ecosystems. This stability may be further reinforced by the high sociality and frequent sharing of roosts and foraging resources among Carollia bats, which facilitate parasite transmission (Altizer et al., Reference Altizer, Nunn, Thrall, Gittleman, Antonovics, Cunningham, Dobson, Ezenwa, Jones, Pedersen and Poss2003; McLellan and Koopman, Reference McLellan, Koopman and Gardner2008; Rifkin et al., Reference Rifkin, Nunn and Garamszegi2012; Webber et al., Reference Webber, McGuire, Smith and Willis2015; Medina and Torres, Reference Medina and Torres2018).
Among other streblid flies recorded, M. proxima and N. bisetosus stand out for exhibiting the broadest host associations, though for distinct biological reasons. M. proxima was the only species parasitizing more than one congeneric host within Sturnira, whereas N. bisetosus exploited two hosts from different genera – A. glaucus and V. caraccioli. Under classical host-specificity categories (Seneviratne et al., Reference Seneviratne, Fernando and Udagama-Randeniya2009), M. proxima qualifies as oligoxenous and N. bisetosus as pleioxenous, the latter being an uncommon pattern in Streblidae, a group known for strong phylogenetic fidelity (Dick and Patterson, Reference Dick, Patterson, Morand, Krasnov and Poulin2006; Autino et al., Reference Autino, Claps, Barquez and Díaz2011).
These broader host associations likely reflect ecological opportunities for cross–host transmission. Although direct evidence for multispecies roost sharing among Sturnira in Andean forests is limited, phyllostomid bats commonly use diverse natural shelters, where mixed-species roosts can occur (Kunz and Lumsden, Reference Kunz, Lumsden, Kunz and Fenton2003; Patterson et al., Reference Patterson, Dick and Dittmar2007). Such conditions plausibly increase contact opportunities among sympatric hosts and may facilitate the movement of M. proxima among closely related Sturnira species. Similarly, rare pleioxenous patterns like that of N. bisetosus have been reported in other Neotropical systems (e.g. Neotrichobius delicatus in Loreto, Peru; Autino et al., Reference Autino, Claps, Barquez and Díaz2011), typically involving flies associated with ecologically overlapping phyllostomid bats (Fagundes et al., Reference Fagundes, Antonini and Aguiar2017).
Host–parasite associations between bat flies and their chiropteran hosts are generally characterized by strong specificity, as seen in Basilia and Anatrichobius, which primarily parasitize Myotis species (Guerrero, Reference Guerrero1995b; Ospina-Pérez et al., Reference Ospina-Pérez, Rivera-Páez and Ramírez-Chaves2023), or Exastinion, apparently restricted to Anoura (Guerrero, Reference Guerrero1995b). Similar patterns have been documented in Peru (Minaya et al., Reference Minaya, Silva and Iannacone2021), and our findings corroborate these associations while extending their known geographic distributions into the montane forests of northern Peru.
Specialization of bat–fly interactions
This study indicates a great specialization at the community level among parasitic bat–flies and their hosts in the montane forests near Nueva Esperanza, Amazonas. These specific associations are well documented in different parts of the world (Lim et al., Reference Lim, Hitch, Lee, Low, Neves, Borthwick, Smith and Mendenhall2020; Poon et al., Reference Poon, Chen, Tsang, Shek, Tsui, Zhao, Guénard and Sin2023); the Neotropical region (Guerrero, Reference Guerrero2019; Ospina-Pérez et al., Reference Ospina-Pérez, Rivera-Páez and Ramírez-Chaves2023; Ramírez-Martínez and Tlapaya-Romero, Reference Ramírez-Martínez and Tlapaya-Romero2023; França et al., Reference França, Alexandre, Correia, Souza, Graciolli, Moura de Souza Aguiar and Vieira2024) and specifically Peru are no exception (Autino et al., Reference Autino, Claps, Barquez and Díaz2011).
Bat–fly interaction networks showed high specificity (high Qw and H2’), a pattern frequently observed in bat–fly interactions at other locations (Fagundes et al., Reference Fagundes, Antonini and Aguiar2017; Urbieta et al., Reference Urbieta, Graciolli and da Cunha Tavares2022; Hiller et al., Reference Hiller, Vollstädt, Brändel, Page and Tschapka2021; Ramalho et al., Reference Ramalho, Diniz and Aguiar2021; Ospina-Pérez et al., Reference Ospina-Pérez, Rivera-Páez and Ramírez-Chaves2023; Ramírez-Martínez and Tlapaya-Romero, Reference Ramírez-Martínez and Tlapaya-Romero2023). This high specificity is mainly driven by the parasitic nature of these interactions, where parasites typically depend strongly on specific hosts to maximize their fitness (Runghen et al., Reference Runghen, Poulin, Monlleó-Borrull and Llopis-Belenguer2021). Such a high degree of dependency often results in parasite–host networks forming modules composed of phylogenetically related species (Felix et al., Reference Felix, Pinheiro, Poulin, Krasnov and Mello2022), as was also observed in all our bat–fly networks. The regional network showed internally nested modules due to the aggregation of interactions that were uniquely observed at specific sites. For example, this pattern is evident in the module of Myotis species: at Higuerón, only M. riparius is present, interacting with A. scorzai and Basilia sp.; at Chontas, M. nigricans appears along with B. anceps; and at Numparket, Anatrichobius sp. is present. When aggregating all these interactions, nestedness increases within the Myotis module, and a similar pattern occurs in other modules, resulting in a compound structure (internally nested modules). This is consistent with the integrative hypothesis of specialization proposed for parasitic networks, which suggests that at larger scales, internally nested modules are more likely to emerge due to the aggregation of different allopatric species and interactions (Felix et al., Reference Felix, Pinheiro, Poulin, Krasnov and Mello2022).
In addition to the evolutionary component behind interactions between bats and flies, these associations have also been particularly discussed in the context of roost-sharing among bat species (Reckardt and Kerth, Reference Reckardt and Kerth2006; Patterson et al., Reference Patterson, Dick and Dittmar2007; Fagundes et al., Reference Fagundes, Antonini and Aguiar2017; Urbieta et al., Reference Urbieta, Graciolli and da Cunha Tavares2022). Logically, species that share roosts are more susceptible to sharing flies, as mentioned in the previous section. This could also be a factor driving the interactions observed in this study, although specific information on roost-sharing is not available for most species. Various bat species share roost with congeneric species, especially in caves (Tanalgo et al., Reference Tanalgo, Tabora and de Oliveira2022), which could have contributed to the independent module aggregation observed in our study for Anoura, Carollia, Myotis and Sturnira, which have been reported to roost in caves frequently (Tanalgo et al., Reference Tanalgo, Tabora and de Oliveira2022). The Stenodermatini species recorded in our study apparently prefer different kinds of roosts (Garbino and Tavares, Reference Garbino and Tavares2018), and as far as we know, there are no records of them roosting together. However, A. glaucus and V. caraccioli were both hosts of P. salvini complex. Artibeus glaucus is a strictly tent-making bat (Ortega et al., Reference Ortega, Arroyo-Cabrales, Martínez-Mendez, Del Real-Monroy, Moreno-Santillán and Velazco2015), and V. caraccioli is also suggested to be a tent-making bat (Page and Dechmann, Reference Page and Dechmann2022). This may suggest they could potentially share roosts; however, tents are usually inhabited by only one species (Rodríguez-Herrera et al., Reference Rodríguez-Herrera, Medellín and Timm2007). Nevertheless, considering that bats can colonize a tent previously used (Rodríguez-Herrera et al., Reference Rodríguez-Herrera, Medellín and Timm2007) or may possibly try to exclude bats from an existent tent (Kunz and McCracken, Reference Kunz and McCracken1996), there is a possibility that flies can be transmitted through tents. This could represent a strategy by flies to spread among populations and species, taking advantage of the complex roosting dynamics of tent-making bats (Chaverri and Kunz, Reference Chaverri and Kunz2006, Reference Chaverri and Kunz2010; Fernandez et al., Reference Fernandez, Schmidt, Schmidt, Rodríguez-Herrera and Knörnschild2021). In summary, roosting behaviour of bats may be closely related to their interactions with flies and should be explored in detail to better understand these relationships.
Specialization and species strength of species varied across sites but followed different patterns, which may indicate that bat–fly relationships change according to specific properties of each site, even though the overall network structure can remain similar (Nielsen and Totland, Reference Nielsen and Totland2014). The observed changes in specialization and species strength for both Carollia species, T. joblingi and P. similis, are consistent with the module separation observed for Carollia at Numparket. At this site, fly species associated with the genus Carollia were much more dependent on C. perspicillata (as indicated by a disproportionately high species strength), while the specialization of P. similis at this site surpassed that of T. joblingi. Although these metrics do not causally drive modularity, they highlight complementary information about how these species structure the network. This suggests that not only species turnover or richness differences can modify networks, but also that changes in how bat and fly species interact can be a driver of subtle network structural variations (Jordán et al., Reference Jordán, Okey, Bauer and Libralato2008; Fründ, Reference Fründ2021).
Despite of the evident species-level variations, rewiring contributed only a small fraction of total dissimilarity in all pairwise comparisons, indicating that species occurring at multiple sites tended to retain similar partners. In contrast, species turnover accounted for more than 90% of link dissimilarity, with the largest contribution coming from the joint turnover of bats and flies, followed by turnover restricted to flies. This pattern reflects both the inherent specificity of flies (Runghen et al., Reference Runghen, Poulin, Monlleó-Borrull and Llopis-Belenguer2021) and the considerable variation in bat and fly assemblages across elevations. These results show that differences in community composition were the main driver of the to the observed variation in interactions, although network structure remained broadly similar among sites as has been observed in other studies (e.g. Kemp et al., Reference Kemp, Evans, Augustyn and Ellis2017; White et al., Reference White, Collier and Stout2022).
In conclusion, we provide novel insights into the diversity and structure of bat–fly interactions in the montane forests of northern Peru and represents the first in the country to apply a network-based approach to these associations. Our records reveal new distributional records at both local and national levels. Bat–fly relationships were highly specialized at both local and regional scales, with slight structural variation across sites. Network structure appears to be shaped by phylogenetic constraints and the roosting behaviour of bat hosts. Species turnover was the major factor behind interaction differences along the elevational gradient. However, species-level roles of bats and flies varied across sites, suggesting that specific interaction dynamics, rather than species turnover alone, contributed to the observed differences in interactions. This also points to a possible interplay between environmental factors and bat–fly relationships. Overall, our findings in this important but previously unexplored region of the Peruvian Andes contribute substantially to the broader ecological understanding of bat–fly interactions in Neotropical ecosystems.
Supplementary material
The supplementary material for this paper can be found at https://doi.org/10.1017/S0031182025101479.
Acknowledgements
The authors would like to thank Christian Olivera Tarifeño, Jhonny Ramos Sandoval and Marlon Hoyos Cerna from the Cordillera Colán National Sanctuary for their support in logistics and advice. We would also like to thank the authorities from the Nueva Esperanza Population Centre for their trust and for allowing us to carry out the work within their jurisdiction. Ricardo Ruiz Chumpitazi, Yanira Zárate Pantoja and Karla Ramos Inga for their support in the field. To our friends at the Nueva Esperanza Population Center, who were very kind and hospitable to us. Especially to the Medina-Leyva family, with special mention to Dalila who at her young age always showed curiosity and interest in the work and supported us in the camps. GG to CNPq (#308119/2022-3). Finally, to IDEA WILD for donating equipment and materials that were of great help.
Author contributions
DM, JJP, CY, KP, BC, JP, GG and JI conceived and designed the study. DM, JJP, CY, KP, BC and JP conducted data gathering. DM and JJP performed statistical analyses. DM, JJP, CY, KP, BC, JP, GG and JI wrote the article.
Financial support
This research received no specific grant from any funding agency, commercial or not-for-profit sectors.
Competing interests
The authors declare there are no conflicts of interest.
Ethical standards
The animals were captured with permission granted by National Forest and Wildlife Service – SERFOR (No. AUT-IFS-2023-048).
Open access
