Cryptic hybridization between the ancient lineages of Natterer's bat (Myotis nattereri).

Studying hybrid zones that form between morphologically cryptic taxa offers valuable insights into the mechanisms of cryptic speciation and the evolution of reproductive barriers. Although hybrid zones have long been the focus of evolutionary studies, the awareness of cryptic hybrid zones increased recently due to rapidly growing evidence of biological diversity lacking obvious phenotypic differentiation. The characterization of cryptic hybrid zones with genome-wide analysis is in its early stages and offers new perspectives for studying population admixture and thus the impact of gene flow. In this study, we investigate the population genomics of the Myotis nattereri complex in one of its secondary contact zones, where a putative hybrid zone is formed between two of its cryptic lineages. By utilizing a whole-genome shotgun sequencing approach, we aim to characterize this cryptic hybrid zone in detail. Demographic analysis suggests that the cryptic lineages diverged during the Pliocene, c. 3.6 million years ago. Despite this ancient separation, the populations in the contact zone exhibit mitochondrial introgression and a considerable amount of mixing in nuclear genomes. The genomic structure of the populations corresponds to geographic locations and the genomic admixture changes along a geographic gradient. These findings suggest that there is no effective hybridization barrier between both lineages, nevertheless, their population structure is shaped by dispersal barriers. Our findings highlight how such deeply diverged cryptic lineages can still readily hybridize in secondary contact.

Agriculture shapes the trophic niche of a bat preying on multiple pest arthropods across Europe: Evidence from DNA metabarcoding.

The interaction between agricultural production and wildlife can shape, and even condition, the functioning of both systems. In this study, we i) explored the degree to which a widespread European bat, namely the common bent-wing bat Miniopterus schreibersii, consumes crop-damaging insects at a continental scale, and ii) tested whether its dietary niche is shaped by the extension and type of agricultural fields. We employed a dual-primer DNA metabarcoding approach to characterize arthropod 16S and COI DNA sequences within bat faecal pellets collected across 16 Southern European localities, to first characterize the bat species' dietary niche, second measure the incidence of agricultural pests across their ranges and third assess whether geographical dietary variation responds to climatic, landscape diversity, agriculture type and vegetation productivity factors. We detected 12 arthropod orders, among which lepidopterans were predominant. We identified >200 species, 44 of which are known to cause agricultural damage. Pest species were detected at all but one sampling site and in 94% of the analysed samples. Furthermore, the dietary diversity of M. schreibersii exhibited a negative linear relation with the area of intensive agricultural fields, thus suggesting crops restrict the dietary niche of bats to prey taxa associated with agricultural production within their foraging range. Overall, our results imply that M. schreibersii might be a valuable asset for biological pest suppression in a variety of agricultural productions and highlight the dynamic interplay between wildlife and agricultural systems.

Molecular and serological survey of lyssaviruses in Croatian bat populations.

BACKGROUND: Rabies is the only known zoonotic disease of bat origin in Europe. The disease is caused by species belonging to the genus Lyssavirus. Five Lyssavirus species, i.e., European bat lyssavirus (EBLV)-1, EBLV-2, Bokeloh bat lyssavirus, Lleida bat lyssavirus, and West Caucasian bat virus, have been identified in European bats. More recently, a proposed sixth species, Kotalahti bat lyssavirus, was detected. Thus, in this study, active surveillance was initiated in order to obtain insights into the prevalence of lyssaviruses in Croatian bat populations and to improve our understanding of the public health threat of infected bats. RESULTS: In total, 455 bats were caught throughout Continental and Mediterranean Croatia. Antibodies were found in 20 of 350 bats (5.71%, 95% confidence interval 3.73-8.66). The majority of seropositive bats were found in Trbusnjak cave (Continental Croatia, Eastern part), and most seropositive bats belonged to Myotis myotis (13/20). All oropharyngeal swabs were negative for the presence of Lyssavirus. CONCLUSIONS: The presence of lyssaviruses in bat populations was confirmed for the first time in Croatia and Southeastern Europe. The results of this study suggest the need for further comprehensive analyses of lyssaviruses in bats in this part of Europe.

Circulation of SARS-CoV-Related Coronaviruses and Alphacoronaviruses in Bats from Croatia.

Bats are natural hosts of various coronaviruses (CoVs), including human CoVs, via an assumed direct zoonotic spillover or intermediate animal host. The present study aimed to investigate the circulation of CoVs in a bat colony in the Mediterranean region of Croatia. Guano and individual droppings from four bat species were sampled and tested with the E-gene sarbecovirus RT-qPCR, the pan-CoV semi-nested RT-PCR targeting the RdRp gene and NGS. Furthermore, bat blood samples were investigated for the presence of sarbecovirus-specific antibodies with the surrogate virus neutralization test (sVNT). The initial testing showed E-gene Sarebeco RT-qPCR reactivity in 26% of guano samples while the bat droppings tested negative. The application of RdRp semi-nested RT-PCR and NGS revealed the circulation of bat alpha- and betaCoVs. Phylogenetic analysis confirmed the clustering of betaCoV sequence with SARS-CoV-related bat sarbecoviruses and alpha-CoV sequences with representatives of the Minunacovirus subgenus. The results of sVNT show that 29% of bat sera originated from all four species that tested positive. Our results are the first evidence of the circulation of SARS-CoV-related coronaviruses in bats from Croatia.

Viral Metagenomic Profiling of Croatian Bat Population Reveals Sample and Habitat Dependent Diversity.

To date, the microbiome, as well as the virome of the Croatian populations of bats, was unknown. Here, we present the results of the first viral metagenomic analysis of guano, feces and saliva (oral swabs) of seven bat species (Myotis myotis, Miniopterus schreibersii, Rhinolophus ferrumequinum, Eptesicus serotinus, Myotis blythii, Myotis nattereri and Myotis emarginatus) conducted in Mediterranean and continental Croatia. Viral nucleic acids were extracted from sample pools, and analyzed using Illumina sequencing. The presence of 63 different viral families representing all seven Baltimore groups were confirmed, most commonly insect viruses likely reflecting the diet of insectivorous bats. Virome compositions of our samples were largely impacted by the sample type: invertebrate-infecting viruses were most frequently found in feces, bacterial viruses in guano, whereas vertebrate-infecting viruses were most common in swabs. Most vertebrate-infecting virus sequences were assigned to retroviruses, parvoviruses, iridoviruses, and poxviruses. We further report the complete genome sequence of a novel adeno-associated virus, densovirus and a near complete length genome sequence of a novel iflavirus. Additionally, one of the most interesting findings in this study was the difference in viromes between two contrasting habitats, the continental and Mediterranean Croatia.

DNA metabarcoding and spatial modelling link diet diversification with distribution homogeneity in European bats.

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.