Population structure of a cave-dwelling bat, Miniopterus schreibersii: does it reflect history and social organization?

Fine scale genetics reveals the subtle negative effects of roads on an endangered bat

The effective management of species with small and fragmented populations requires an in-depth understanding of how the effects of human-induced habitat disturbance shape the structure and gene flow at fine spatial scales. Identification of putative environmental barriers that affect individual exchange among subpopulations is imperative to prevent extinction risks. Here, we investigated how landscape affects the gene flow and relatedness structure of a population of the endangered lesser horseshoe bat (Rhinolophus hipposideros). We also assessed the effects of sexbiased dispersal on genetic relatedness. We genotyped 287 bat samples collected across southern Portugal and developed resistance surfaces for landscape variables hypothesized to affect gene flow. Then, we used spatially explicit models to fit relatedness distance through the resistance surfaces. We found genetic evidence of sex-biased dispersal and identified a significant fine scale structuring in the relatedness regarding females, the philopatric sex. Males displayed uniform levels of relatedness throughout the landscape. The results indicated less relatedness between the female´ from roosts located on proximity of roads than in roosts away from roads. Also, when analysing the sexes together the relatedness on roosts separated by highway were subtly less related in comparison to those occurring on the same side. Roads seem to be major shapers of the contemporary population structure of females, regardless of being relatively recent structures in the landscape. Furthermore, the relatedness patterns detected suggested that high tree density among roosts and continuity of forest patches in broader surrounding areas, promotes the relatedness among individuals. Landscape heterogeneity among roosts slightly decreases genetic relatedness. Nevertheless, those relationships are still weak, suggesting that population structuring driven by those factors is slowly ongoing. Thus, effective management measures should focus on issues for promoting safe road passages and suitable habitat corridors, allowing for the exchange of individuals and gene flow among lesser horseshoe bat roosts.

Inferring genetic structure when there is little: population genetics versus genomics of the threatened bat Miniopterus schreibersii across Europe

Despite their paramount importance in molecular ecology and conservation, genetic diversity and structure remain challenging to quantify with traditional genotyping methods. Next-generation sequencing holds great promises, but this has not been properly tested in highly mobile species. In this article, we compared microsatellite and RAD-sequencing (RAD-seq) analyses to investigate population structure in the declining bent-winged bat (Miniopterus schreibersii) across Europe. Both markers retrieved general patterns of weak range-wide differentiation, little sex-biased dispersal, and strong isolation by distance that associated with significant genetic structure between the three Mediterranean Peninsulas, which could have acted as glacial refugia. Microsatellites proved uninformative in individual-based analyses, but the resolution offered by genomic SNPs illuminated on regional substructures within several countries, with colonies sharing migrators of distinct ancestry without admixture. This finding is consistent with a marked philopatry and spatial partitioning between mating and rearing grounds in the species, which was suspected from marked-recaptured data. Our study advocates that genomic data are necessary to properly unveil the genetic footprints left by biogeographic processes and social organization in long-distant flyers, which are otherwise rapidly blurred by their high levels of gene flow.

Focus on Hyperparasites: Biotic and Abiotic Traits Affecting the Prevalence of Parasitic Microfungi on Bat Ectoparasites

The tritrophic association of bats, bat flies, and Laboulbeniales microfungi is a remarkably understudied system that may reveal patterns applicable to community ecology theory of (hyper)parasites. Laboulbeniales are biotrophic microfungi, exclusively associated with arthropods, with several species that are specialized on bat flies, which themselves are permanent ectoparasites of bats. Several hypotheses were tested on biotic and abiotic traits that may influence the presence and prevalence of hyperparasitic Laboulbeniales fungi on bat flies, based on southeastern European data. We found a wide distribution of fungal infection on bat flies, with underground-dwelling bats hosting more Laboulbeniales-infected flies compared to crevice-dwelling species. Bat host behavior, sociality, roost selection (underground versus crevice), bat fly sex, and season all have significant effects on the prevalence of fungal infection. Laboulbeniales infections are more common on bat flies that are infecting bat species with dense and long-lasting colonies (Miniopterus schreibersii, Myotis myotis, Myotis blythii), which roost primarily in underground sites. Inside these sites, elevated temperature and humidity may enhance the development and transmission of Laboulbeniales fungi. Sexual differences in bat hosts’ behavior also have an effect on fungal infection risk, with densely roosting female bat hosts harboring more Laboulbeniales-infected bat flies.

Sharing roosts but not ectoparasites: high host-specificity in bat flies and wing mites of Miniopterus schreibersii and Rhinolophus ferrumequinum (Mammalia: Chiroptera)

Schreiber’s bent-winged bat Miniopterus schreibersii and the greater horseshoe bat Rhinolophus ferrumequinum are widespread and common cavernicolous species across southern Europe that host numerous specialized ectoparasite species. The objective of this study was to characterize the species assemblage, genetic diversity, and host specificity of bat flies (Nycteribiidae, Diptera) and wing mites (Spinturnicidae, Acari) found on these bat hosts in Serbia and Bosnia and Herzegovina. Notably, while bat flies lay puparia on the cave walls and can thus be transmitted indirectly, wing mites require direct body contact for transmission. Morphological identification and sequencing of a 710-bp fragment of cytochrome oxidase I gene of 207 bat flies yielded 4 species, 3 on M. schreibersii and 1 on R. ferrumequinum. Sequencing of a 460-bp small subunit ribosomal RNA fragment, in all 190 collected wing mites revealed 2 species, 1 per host. In no case was a parasite associated with 1 host found on the other host. Species and genetic diversity of flies were higher in M. schreibersii, likely reflecting their host’s larger colony sizes and migratory potential. Mite species of both hosts showed similarly low diversity, likely due to their faster life history and lower winter survival. Our findings highlight a remarkably high host-specificity and segregation of ectoparasite species despite direct contact among their hosts in the roost, suggesting a defined host preference in the investigated ectoparasite species. Furthermore, the differences in ectoparasite genetic diversity exemplify the interplay between host and parasite life histories in shaping parasite population genetic structure.

Out of Africa: The origins of the protozoan blood parasites of the Trypanosoma cruzi clade found in bats from Africa

Understanding geographic patterns of interaction between hosts and parasites can provide useful insight into the evolutionary history of the organisms involved. However, poor taxon sampling often hinders meaningful phylogenetic descriptions of groups of parasites. Trypanosome parasites that constitute the Trypanosoma cruzi clade are worldwide distributed infecting several mammalian species, especially bats. Diversity in this clade has been recently expanded by newly discovered species, but the common ancestor and geographical origins of this group of blood parasites are still debated. We present here results based on the molecular characterization of trypanosome isolates obtained from 1493 bats representing 74 species and sampled over 16 countries across four continents. After estimating the appropriate number of hypothetical species in our data set using GMYC models in combination with Poisson Tree Processes (mPTP) and ABGD, the 18S rRNA and gGAPDH genes were used for phylogenetic analyses to infer the major evolutionary relationships in the T. cruzi clade. Then, biogeographical processes influencing the distribution of this cosmopolitan group of parasites was inferred using BioGeoBEARS. Results revealed a large lineages diversity and the presence of trypanosomes in all sampled regions which infected 344 individuals from 31 bat species. We found eight Trypanosoma species, including: five previously known; one subspecies of Trypanosoma livingstonei (Trypanosoma cf. livingstonei); and two undescribed taxa (Trypanosoma sp. 1, Trypanosoma sp. 2), which were found exclusively in bats of the genus Miniopterus from Europe and Africa. The new taxa discovered have both an unexpected position in the global phylogeny of the T. cruzi clade. Trypanosoma sp. 1 is a sister lineage of T. livingstonei which is located at the base of the tree, whereas Trypanosoma sp. 2 is a sister lineage of the Shizotrypanum subclade that contains T. c. cruzi and T. dionisii. Ancestral areas reconstruction provided evidence that trypanosomes of the T. cruzi clade have radiated from Africa through several dispersion events across the world. We discuss the impact of these findings on the biogeography and taxonomy of this important clade of parasites and question the role played by bats, especially those from the genus Miniopterus, on the dispersal of these protozoan parasites between continents.

A continent-scale study of the social structure and phylogeography of the bent-wing bat, Miniopterus schreibersii (Mammalia: Chiroptera), using new microsatellite data

Miniopterus schreibersii is a cave-dwelling bat species with a wide distribution in the western Palearctic spanning southern and central Europe, North Africa, Anatolia, the Middle East, and the Caucasus. We investigated the social structure and its effects on the genetic makeup of this species, using 10 nuclear microsatellite markers and a partial fragment of the mitochondrial cytochrome b gene. Samples were examined from the species' entire circum-Mediterranean range. Local structuring that was previously detected among populations of M. schreibersii using mitochondrial markers was not observed for microsatellite markers, indicating male-biased dispersal for the species. Some support was found for postglacial expansions in Europe, with Anatolia potentially acting as the primary refugium during the Last Glacial Maximum (LGM). However, support for this hypothesis is not as strong as that previously detected using mitochondrial DNA markers. This is likely due to the diminishing effect of male-mediated dispersal, replenishing the nuclear diversity faster than the mitochondrial diversity in regions that are relatively far from the glacial refugia.

Eco-epidemiology of Novel Bartonella Genotypes from Parasitic Flies of Insectivorous Bats

Bats are important zoonotic reservoirs for many pathogens worldwide. Although their highly specialized ectoparasites, bat flies (Diptera: Hippoboscoidea), can transmit Bartonella bacteria including human pathogens, their eco-epidemiology is unexplored. Here, we analyzed the prevalence and diversity of Bartonella strains sampled from 10 bat fly species from 14 European bat species. We found high prevalence of Bartonella spp. in most bat fly species with wide geographical distribution. Bat species explained most of the variance in Bartonella distribution with the highest prevalence of infected flies recorded in species living in dense groups exclusively in caves. Bat gender but not bat fly gender was also an important factor with the more mobile male bats giving more opportunity for the ectoparasites to access several host individuals. We detected high diversity of Bartonella strains (18 sequences, 7 genotypes, in 9 bat fly species) comparable with tropical assemblages of bat-bat fly association. Most genotypes are novel (15 out of 18 recorded strains have a similarity of 92-99%, with three sequences having 100% similarity to Bartonella spp. sequences deposited in GenBank) with currently unknown pathogenicity; however, 4 of these sequences are similar (up to 92% sequence similarity) to Bartonella spp. with known zoonotic potential. The high prevalence and diversity of Bartonella spp. suggests a long shared evolution of these bacteria with bat flies and bats providing excellent study targets for the eco-epidemiology of host-vector-pathogen cycles.

Geographical variation in the echolocation calls of bent-winged bats, Miniopterus fuliginosus

Evolutionary biologists had a long-standing interest in the evolutionary forces underlying geographical variation in the acoustic signals of animals. However, the evolutionary forces driving acoustic variation are still unclear. In this study, we quantified the geographical variation in the peak frequencies of echolocation calls in eight Miniopterus fuliginosus bat colonies, and assessed the forces that drive acoustic divergence. Our results demonstrated that seven of the colonies had very similar peak frequencies, while only one colony was significantly higher than the others. This similarity in echolocation call frequency among the seven colonies was likely due to frequent dispersal and migration, leading to male-mediated infiltration of nuclear genes. This infiltration enhances gene flow and weakens ecological selection, and also increases interactions in the presence of conspecifics. Significant correlations were not observed between acoustic distances and morphological distances, climatic differences, geographic distances or mtDNA genetic distances. However, variation in acoustic distances was significantly positive correlated with nDNA genetic distance, even after controlling for geographic distance. Interestingly, the relationship between call divergence and genetic distance was no longer significant after excluding the colony with the highest call frequency, which may be due to the minimal genetic distance among the other seven colonies. The highest frequencies of echolocation calls observed in the one colony may be shaped by selection pressure due to loud background noise in the area. Taken together, these results suggest that geographic divergence of echolocation calls may not be subject to genetic drift, but rather, that the strong selective pressure induced by background noise may lead to acoustic and genetic differentiation between JXT and the other colonies.

Combined hybridization and mitochondrial capture shape complex phylogeographic patterns in hybridogenetic Cataglyphis desert ants

Some species of Cataglyphis desert ants have evolved a hybridogenetic mode of reproduction at the social scale. In hybridogenetic populations, two distinct genetic lineages coexist. Non-reproductive offspring (workers) are hybrids of the two lineages, whereas sexual offspring (males and new queens) are produced by parthenogenesis and belong to the mother queen lineage. How this unusual reproductive system affects phylogeographic patterns and speciation processes remains completely unknown to date. Using one mitochondrial and four nuclear genes, we examined the phylogenetic relationships between three species of Cataglyphis (C. hispanica, C. humeya and C. velox) where complex DNA inheritance through social hybridogenesis may challenge phylogenetic inference. Our results bring two important insights. First, our data confirm a hybridogenetic mode of reproduction across the whole distribution range of the species C. hispanica. In contrast, they do not provide support for hybridogenesis in the populations sampled of C. humeya and C. velox. This suggests that these populations are not hybridogenetic, or that hybridogenesis is too recent to result in reciprocally monophyletic lineages on nuclear genes. Second, due to mitochondrial introgression between lineages (Darras and Aron, 2015), the faster-evolving COI marker is not lineage specific, hence, unsuitable to further investigate the segregation of lineages in the species studied. Different mitochondrial haplotypes occur in each locality sampled, resulting in strongly structured populations. This micro-allopatric structure leads to over-splitting species delimitation on mitochondrial gene, as every locality could potentially be considered a putative species; haploweb analyses of nuclear markers, however, yield species delimitations that are consistent with morphology. Overall, this study highlights how social hybridogenesis varies across species and shapes complex phylogeographic patterns.

Detection of the Severe Acute Respiratory Syndrome-Related Coronavirus and Alphacoronavirus in the Bat Population of Taiwan

Bats have been demonstrated to be natural reservoirs of severe acute respiratory syndrome coronavirus (SARS CoV) and Middle East respiratory syndrome (MERS) CoV. Faecal samples from 248 individuals of 20 bat species were tested for partial RNA‐dependent RNA polymerase gene of CoV and 57 faecal samples from eight bat species were tested positive. The highest detection rate of 44% for Scotophilus kuhlii, followed by 30% for Rhinolophus monoceros. Significantly higher detection rates of coronaviral RNA were found in female bats and Scotophilus kuhlii roosting in palm trees. Phylogenetic analysis classified the positive samples into SARS‐related (SARSr) CoV, Scotophilus bat CoV 512 close to those from China and Philippines, and Miniopterus bat CoV 1A‐related lineages. Coronaviral RNA was also detected in bat guano from Scotophilus kuhlii and Myotis formosus flavus on the ground and had potential risk for human exposure. Diverse bat CoV with zoonotic potential could be introduced by migratory bats and maintained in the endemic bat population in Taiwan.

Diversity, Host Specialization, and Geographic Structure of Filarial Nematodes Infecting Malagasy Bats

We investigated filarial infection in Malagasy bats to gain insights into the diversity of these parasites and explore the factors shaping their distribution. Samples were obtained from 947 individual bats collected from 52 sites on Madagascar and representing 31 of the 44 species currently recognized on the island. Samples were screened for the presence of micro- and macro-parasites through both molecular and morphological approaches. Phylogenetic analyses showed that filarial diversity in Malagasy bats formed three main groups, the most common represented by Litomosa spp. infecting Miniopterus spp. (Miniopteridae); a second group infecting Pipistrellus cf. hesperidus (Vespertilionidae) embedded within the Litomosoides cluster, which is recognized herein for the first time from Madagascar; and a third group composed of lineages with no clear genetic relationship to both previously described filarial nematodes and found in M. griveaudi, Myotis goudoti, Neoromicia matroka (Vespertilionidae), Otomops madagascariensis (Molossidae), and Paratriaenops furculus (Hipposideridae). We further analyzed the infection rates and distribution pattern of Litomosa spp., which was the most diverse and prevalent filarial taxon in our sample. Filarial infection was disproportionally more common in males than females in Miniopterus spp., which might be explained by some aspect of roosting behavior of these cave-dwelling bats. We also found marked geographic structure in the three Litomosa clades, mainly linked to bioclimatic conditions rather than host-parasite associations. While this study demonstrates distinct patterns of filarial nematode infection in Malagasy bats and highlights potential drivers of associated geographic distributions, future work should focus on their alpha taxonomy and characterize arthropod vectors.

Mitochondrial genetic differentiation and morphological difference of Miniopterus fuliginosus and Miniopterus magnater in China and Vietnam

Because of its complicated systematics, the bent-winged bat is one of the most frequently studied bat species groups. In China, two morphologically similar bent-winged bat species, Miniopterus fuliginosus and Miniopterus magnater were identified, but their distribution range and genetic differentiation are largely unexplored. In this study, we applied DNA bar codes and two other mitochondrial DNA genes including morphological parameters to determine the phylogeny, genetic differentiation, spatial distribution, and morphological difference of the M. fuliginosus and M. magnater sampled from China and one site in Vietnam. Mitochondrial DNA gene genealogies revealed two monophyletic lineages throughout the Tropic of Cancer. According to DNA bar code divergences, one is M. fuliginosus corresponding to the Chinese mainland and the other is M. magnater corresponding to tropical regions including Hainan and Guangdong provinces of China and Vietnam. Their most recent common ancestor was dated to the early stage of the Quaternary glacial period (ca. 2.26 million years ago [Ma] on the basis of D-loop data, and ca. 1.69-2.37 Ma according to ND2). A population expansion event was inferred for populations of M. fuliginosus at 0.14 Ma. The two species probably arose in separate Pleistocene refugia under different climate zones. They significantly differed in forearm length, maxillary third molar width, and greatest length of the skull.

How a haemosporidian parasite of bats gets around: the genetic structure of a parasite, vector and host compared

Parasite population structure is often thought to be largely shaped by that of its host. In the case of a parasite with a complex life cycle, two host species, each with their own patterns of demography and migration, spread the parasite. However, the population structure of the parasite is predicted to resemble only that of the most vagile host species. In this study, we tested this prediction in the context of a vector-transmitted parasite. We sampled the haemosporidian parasite Polychromophilus melanipherus across its European range, together with its bat fly vector Nycteribia schmidlii and its host, the bent-winged bat Miniopterus schreibersii. Based on microsatellite analyses, the wingless vector, and not the bat host, was identified as the least structured population and should therefore be considered the most vagile host. Genetic distance matrices were compared for all three species based on a mitochondrial DNA fragment. Both host and vector populations followed an isolation-by-distance pattern across the Mediterranean, but not the parasite. Mantel tests found no correlation between the parasite and either the host or vector populations. We therefore found no support for our hypothesis; the parasite population structure matched neither vector nor host. Instead, we propose a model where the parasite's gene flow is represented by the added effects of host and vector dispersal patterns.

Taxonomic evaluation of the bent-winged bat (Miniopterus) populations occurring in Iran inferred from mitochondrial cytochrome-b sequences

We studied the genetic differentiation and taxonomy of the bent-wing bat (

Detection of bat coronaviruses from Miniopterus fuliginosus in Japan

Bats have great potential as reservoirs for emerging viruses such as severe acute respiratory syndrome-coronavirus. In this study, bat coronaviruses (BtCoVs) were detected by RT-PCR from intestinal and fecal specimens of Miniopterus fuliginosus breeding colonies in Wakayama Prefecture caves, where we previously identified bat betaherpesvirus 2. Two primer sets were used for the detection of BtCoV: one was for the RNA-dependent RNA polymerase (RdRp) region and the other was for the spike (S) protein region. Eleven and 73% of intestinal and fecal specimens, respectively, were positive for RdRp region, and 2 and 40% of those were positive for S protein region. Sequencing and phylogenetic analysis showed that the detected BtCoV belonged to the group 1 (alpha) coronaviruses. These data suggest that BtCoV is endemic in M. fuliginosus in Japan.

Wings or winds: inferring bat migration in a stepping-stone archipelago

Eocene ocean currents and prevailing winds correlate with over-water dispersals of terrestrial mammals from Africa to Madagascar. Since the Early Miocene (about 23 Ma), these currents flowed in the reverse direction, from the Indian Ocean towards Africa. The Comoro Islands are equidistant between Africa and Madagascar and support an endemic land vertebrate fauna that shares recent ancestry predominantly with Madagascar. We examined whether gene flow in two Miniopterus bat species endemic to the Comoros and Madagascar correlates with the direction of current winds, using uni- and bi-parentally inherited markers with different evolutionary rates. Coalescence-based analyses of mitochondrial matrilines support a Pleistocene (approximately 180,000 years ago) colonization event from Madagascar west to the Comoros (distance: 300 km) in the predicted direction. However, nuclear microsatellites show that more recent gene flow is restricted to a few individuals flying against the wind, from Grande Comore to Anjouan (distance: 80 km).

Cryptic diversity of the bent-wing bat, Miniopterus schreibersii (Chiroptera: Vespertilionidae), in Asia Minor

Background

Two or more species are cryptic, if they are morphologically similar, biologically distinct, and misclassified as a single species. Cryptic species complexes were recently discovered within many bat species and we suspect that the bent-wing bat, Miniopterus schreibersii, found in Europe, northern Africa, and Asia Minor, could also form such a complex. Populations of M. schreibersii decline in most of the European countries and the species is currently listed as Near Threatened in the IUCN Red List. Finding that M. schreibersii is not a single species, but a species complex, would have a considerable impact on its conservation strategies, as the abundance of each component taxon would be much smaller than the one estimated for the nominal species.

Results

Miniopterus schreibersii in Asia Minor consists of two genetically diverged lineages, which are reciprocally monophyletic on three mitochondrial DNA markers, have a diagnostic set of multilocus allele frequencies, and show a marked difference in their population structures. The lineages differ slightly in their size, wing shape, and echolocation call parameters. Although these differences are sufficient to discriminate between the lineages, they are not fully diagnostic in reference to individuals. We suggest that the lineages endured the major Northern Hemisphere glaciations in different glacial refugia and colonized Asia Minor after the last glacial maximum. The lineages are allopatric, which is neither delineated by the presence of geographical barriers nor associated with the specific climatic conditions, and which we link to competitive exclusion.

Conclusions

The distinctions between the lineages comply with most of the criteria required for species delineation imposed by various species concepts. Accordingly, we conclude that M. schreibersii in Asia Minor is represented by two cryptic species. Our results imply that the distributional range of the nominal species is almost exclusively limited to Europe and the coastal zones of Asia Minor. As populations of M. schreibersii seem to be much smaller than currently assumed, conservation strategies regarding this taxon need to be revised. The exact distributional range and the vulnerability of the suggested sister species to M. schreibersii is yet to be assessed.