Bat coronaviruses (CoV) are putative precursors of the severe acute respiratory syndrome (SARS) CoV and other CoV that crossed the species barrier from zoonotic reservoirs into the human population. To determine the presence and distribution of CoV in Iberian bats, 576 individuals of 26 different bat species were captured in 13 locations in Spain. We report for the first time the presence of 14 coronaviruses in 9 Iberian bat species. Phylogenetic analysis of a conserved CoV genome region (RdRp gene) shows a wide diversity and distribution of alpha and betacoronavirus in Spain. Interestingly, although some of these viruses are related to other European BatCoV, or to Asian CoV, some of the viruses found in Spain cluster in new groups of α and β CoV.

Papillomaviruses (PVs) are widespread pathogens. However, the extent of PV infections in bats remains largely unknown. This work represents the first comprehensive study of PVs in Iberian bats. We identified four novel PVs in the mucosa of free-ranging Eptesicus serotinus (EserPV1, EserPV2, and EserPV3) and Rhinolophus ferrumequinum (RferPV1) individuals and analyzed their phylogenetic relationships within the viral family. We further assessed their prevalence in different populations of E. serotinus and its close relative E. isabellinus. Although it is frequent to read that PVs co-evolve with their host, that PVs are highly species-specific, and that PVs do not usually recombine, our results suggest otherwise. First, strict virus–host co-evolution is rejected by the existence of five, distantly related bat PV lineages and by the lack of congruence between bats and bat PVs phylogenies. Second, the ability of EserPV2 and EserPV3 to infect two different bat species (E. serotinus and E. isabellinus) argues against strict host specificity. Finally, the description of a second noncoding region in the RferPV1 genome reinforces the view of an increased susceptibility to recombination in the E2-L2 genomic region. These findings prompt the question of whether the prevailing paradigms regarding PVs evolution should be reconsidered.

Molecular analysis of diet overcomes the considerable limitations of traditional techniques for identifying prey remains in bat faeces. We collected faeces from individual Mountain Long-eared Bats Plecotus macrobullaris trapped using mist nets during the summers of 2009 and 2010 in the Pyrenees. We analysed their diet using DNA mini-barcodes to identify prey species. In addition, we inferred some basic features of the bat's foraging ecology that had not yet been addressed. P. macrobullaris fed almost exclusively on moths (97.8%). As prey we detected one dipteran genus (Tipulidae) and 29 moth taxa: 28 were identified at species level (23 Noctuidae, 1 Crambidae, 1 Geometridae, 1 Pyralidae, 1 Sphingidae, 1 Tortricidae), and one at genus level (Rhyacia sp., Noctuidae). Known ecological information about the prey species allowed us to determine that bats had foraged at elevations between 1,500 and 2,500 m amsl (above mean sea level), mostly in subalpine meadows, followed by other open habitats such as orophilous grasslands and alpine meadows. No forest prey species were identified in the diet. As 96.4% of identified prey species were tympanate moths and no evidence of gleaning behaviour was revealed, we suggest P. macrobullaris probably forages by aerial hawking using faint echolocation pulses to avoid detection by hearing moths. As we could identify 87.8% of the analysed sequences (64.1% of the MOTUs, Molecular Operational Taxonomic Units) at species level, we conclude that DNA mini-barcodes are a very useful tool to analyse the diet of moth-specialist bats.

Niche partitioning through foraging is a mechanism likely involved in facilitating the coexistence of ecologically similar and co-occurring animal species by separating their use of resources. Yet, this mechanism is not well understood in flying insectivorous animals. This is particularly true of bats, where many ecologically similar or cryptic species coexist. The detailed analysis of the foraging niche in sympatric, cryptic sibling species provides an excellent framework to disentangle the role of specific niche factors likely involved in facilitating coexistence. We used DNA metabarcoding to determine the prey species consumed by a population of sympatric sibling Rhinolophus euryale and Rhinolophus mehelyi whose use of habitat in both sympatric and allopatric ranges has been well established through radio tracking. Although some subtle dietary differences exist in prey species composition, the diet of both bats greatly overlapped (O = 0.83) due to the consumption of the same common and widespread moths. Those dietary differences we did detect might be related to divergences in prey availabilities among foraging habitats, which prior radio tracking on the same population showed are differentially used and selected when both species co-occur. This minor dietary segregation in sympatry may be the result of foraging on the same prey-types and could contribute to reduce potential competitive interactions (e.g., for prey, acoustic space). Our results highlight the need to evaluate the spatial niche dimension in mediating the co-occurrence of similar insectivorous bat species, a niche factor likely involved in processes of bat species coexistence.

Inferences of the interactions between species' ecological niches and spatial distribution have been historically based on simple metrics such as low-resolution dietary breadth and range size, which might have impeded the identification of meaningful links between niche features and spatial patterns. We analysed the relationship between dietary niche breadth and spatial distribution features of European bats, by combining continent-wide DNA metabarcoding of faecal samples with species distribution modelling. Our results show that while range size is not correlated with dietary features of bats, the homogeneity of the spatial distribution of species exhibits a strong correlation with dietary breadth. We also found that dietary breadth is correlated with bats' hunting flexibility. However, these two patterns only stand when the phylogenetic relations between prey are accounted for when measuring dietary breadth. Our results suggest that the capacity to exploit different prey types enables species to thrive in more distinct environments and therefore exhibit more homogeneous distributions within their ranges.