Population genetic structure of serotine bats (Eptesicus serotinus) across Europe and implications for the potential spread of bat rabies (European bat lyssavirus EBLV-1)

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.

Genetic identification of bat species for pathogen surveillance across France

With more than 1400 chiropteran species identified to date, bats comprise one-fifth of all mammalian species worldwide. Many studies have associated viral zoonoses with 45 different species of bats in the EU, which cluster within 5 families of bats. For example, the Serotine bats are infected by European Bat 1 Lyssavirus throughout Europe while Myotis bats are shown infected by coronavirus, herpesvirus and paramyxovirus. Correct host species identification is important to increase our knowledge of the ecology and evolutionary pattern of bat viruses in the EU. Bat species identification is commonly determined using morphological keys. Morphological determination of bat species from bat carcasses can be limited in some cases, due to the state of decomposition or nearly indistinguishable morphological features in juvenile bats and can lead to misidentifications. The overall objective of our study was to identify insectivorous bat species using molecular biology tools with the amplification of the partial cytochrome b gene of mitochondrial DNA. Two types of samples were tested in this study, bat wing punches and bat faeces. A total of 163 bat wing punches representing 22 species, and 31 faecal pellets representing 7 species were included in the study. From the 163 bat wing punches tested, a total of 159 were genetically identified from amplification of the partial cyt b gene. All 31 faecal pellets were genetically identified based on the cyt b gene. A comparison between morphological and genetic determination showed 21 misidentifications from the 163 wing punches, representing ~12.5% of misidentifications of morphological determination compared with the genetic method, across 11 species. In addition, genetic determination allowed the identification of 24 out of 25 morphologically non-determined bat samples. Our findings demonstrate the importance of a genetic approach as an efficient and reliable method to identify bat species precisely.

Incursion of European Bat Lyssavirus 1 (EBLV-1) in Serotine Bats in the United Kingdom

Lyssaviruses are an important genus of zoonotic viruses which cause the disease rabies. The United Kingdom is free of classical rabies (RABV). However, bat rabies due to European bat lyssavirus 2 (EBLV-2), has been detected in Daubenton’s bats (Myotis daubentonii) in Great Britain since 1996, including a fatal human case in Scotland in 2002. Across Europe, European bat lyssavirus 1 (EBLV-1) is commonly associated with serotine bats (Eptesicus serotinus). Despite the presence of serotine bats across large parts of southern England, EBLV-1 had not previously been detected in this population. However, in 2018, EBLV-1 was detected through passive surveillance in a serotine bat from Dorset, England, using a combination of fluorescent antibody test, reverse transcription-PCR, Sanger sequencing and immunohistochemical analysis. Subsequent EBLV-1 positive serotine bats have been identified in South West England, again through passive surveillance, during 2018, 2019 and 2020. Here, we confirm details of seven cases of EBLV-1 and present similarities in genetic sequence indicating that emergence of EBLV-1 is likely to be recent, potentially associated with the natural movement of bats from the near continent

Population structure, gene flow and relatedness of Natterer’s bats in Northern England

There have been significant declines in population numbers of many bat species in the United Kingdom, including Natterer’s batsMyotis nattereri, over the last century, largely due to anthropogenic changes. The philopatry, which temperate-zone bats often exhibit to their natal landscapes, in combination with anthropogenic threats, can lead to fragmentation, isolation and sub-division of populations. This may result in bottlenecks and declines in genetic diversity. Multi-scaled research is required to disentangle how the variation in the physical traits of bat species (e.g. affecting flight), as well as their social and behavioural traits (e.g. community size, migration, breeding systems), may affect the genetic health of populations and provide a potential buffer against fragmentation. We used microsatellite markers to characterise the genetic diversity and population structure present in Natterer’s bat colonies to determine whether summer roosting bat colonies were spatially differentiated or part of a meta-population. Analyses of population structure and measures of genetic relatedness suggest spatially differentiated populations of bats exhibit long term site fidelity to summer roosting sites, whilst high genetic diversity at sites indicates gene exchange occurs via swarming sites. Natterer’s bats in northern England may travel greater distances to swarming sites than has been previously documented.

Phylogeny of Hungarian EBLV-1 strains using whole-genome sequence data

European bat lyssavirus 1 (EBLV-1) is a widespread lyssavirus across Europe, whose epizootic cycle is linked to a few bat species. Occasionally, EBLV-1 infection may occur in domestic animals and humans. EBLV-1 can be classified into two subtypes, where subtype EBLV-1a shows a wide geographic distribution between France and Russia whereas subtype EBLV-1b is distributed between Spain and Poland. In this study, we determined the genome sequence of two recent EBLV-1a strains detected in Hungary and analysed their adaptive evolution and phylodynamics. The data set that included 100 EBLV-1 genome sequences identified positive selection at selected sites in genes coding for viral proteins (N, codon 18; P, 141 and 155; G, 244 and 488; L, 168, 980, 1597 and 1754). A major genetic clade containing EBLV-1a isolates from Hungary, Slovakia, Denmark and Poland was estimated to have diverged during the 19th century whereas the divergence of the most recent ancestor of Hungarian and Slovakian isolates dates back to 1950 (time span, 1930 to 1970). Phylogeographic analysis of the EBLV-1a genomic sequences demonstrated strong evidence of viral dispersal from Poland to Hungary. This new information indicates that additional migratory flyways may help the virus spread, a finding that supplements the general theory on a west-to-east dispersal of EBLV-1a strains. Long-distance migrant bats may mediate the dispersal of EBLV-1 strains across Europe; however, structured surveillance and extended genome sequencing would be needed to better understand the epizootiology of EBLV-1 infections in Europe.

Integrating population genetics to define conservation units from the core to the edge of Rhinolophus ferrumequinum western range

The greater horseshoe bat (Rhinolophus ferrumequinum) is among the most widespread bat species in Europe but it has experienced severe declines, especially in Northern Europe. This species is listed Near Threatened in the European IUCN Red List of Threatened Animals, and it is considered to be highly sensitive to human activities and particularly to habitat fragmentation. Therefore, understanding the population boundaries and demographic history of populations of this species is of primary importance to assess relevant conservation strategies. In this study, we used 17 microsatellite markers to assess the genetic diversity, the genetic structure, and the demographic history of R. ferrumequinum colonies in the western part of its distribution. We identified one large population showing high levels of genetic diversity and large population size. Lower estimates were found in England and northern France. Analyses of clustering and isolation by distance suggested that the Channel and the Mediterranean seas could impede R. ferrumequinum gene flow. These results provide important information to improve the delineation of R. ferrumequinum management units. We suggest that a large management unit corresponding to the population ranging from Spanish Basque Country to northern France must be considered. Particular attention should be given to mating territories as they seem to play a key role in maintaining high levels of genetic mixing between colonies. Smaller management units corresponding to English and northern France colonies must also be implemented. These insular or peripheral colonies could be at higher risk of extinction in the near future.

The Gambian epauletted fruit bat shows increased genetic divergence in the Ethiopian highlands and in an area of rapid urbanization

The Gambian epauletted fruit bat (Epomophorus gambianus) is an abundant species that roosts in both urban and rural settings. The possible role of E. gambianus as a reservoir host of zoonotic diseases underlines the need to better understand the species movement patterns. So far, neither observational nor phylogenetic studies have identified the dispersal range or behavior of this species. Comparative analyses of mitochondrial and nuclear markers from 20 localities across the known distribution of E. gambianus showed population panmixia, except for the populations in Ethiopia and southern Ghana (Accra and Ve-Golokwati). The Ethiopian population may be ancestral and is highly divergent to the species across the rest of its range, possibly reflecting isolation of an ancient colonization along an east-west axis. Mitochondrial haplotypes in the Accra population display a strong signature of a past bottleneck event; evidence of either an ancient or recent bottleneck using microsatellite data, however, was not detected. Demographic analyses identified population expansion in most of the colonies, except in the female line of descent in the Accra population. The molecular analyses of the colonies from Ethiopia and southern Ghana show gender dispersal bias, with the mitochondrial DNA fixation values over ten times those of the nuclear markers. These findings indicate free mixing of the species across great distances, which should inform future epidemiological studies.

Host Genetic Variation Does Not Determine Spatio-Temporal Patterns of European Bat 1 Lyssavirus

The majority of bat rabies cases in Europe are attributed to European bat 1 lyssavirus (EBLV-1), circulating mainly in serotine bats (Eptesicus serotinus). Two subtypes have been defined (EBLV-1a and EBLV-1b), each associated with a different geographical distribution. In this study, we undertake a comprehensive sequence analysis based on 80 newly obtained EBLV-1 nearly complete genome sequences from nine European countries over a 45-year period to infer selection pressures, rates of nucleotide substitution, and evolutionary time scale of these two subtypes in Europe. Our results suggest that the current lineage of EBLV-1 arose in Europe ∼600 years ago and the virus has evolved at an estimated average substitution rate of ∼4.19×10-5 subs/site/year, which is among the lowest recorded for RNA viruses. In parallel, we investigate the genetic structure of French serotine bats at both the nuclear and mitochondrial level and find that they constitute a single genetic cluster. Furthermore, Mantel tests based on interindividual distances reveal the absence of correlation between genetic distances estimated between viruses and between host individuals. Taken together, this indicates that the genetic diversity observed in our E. serotinus samples does not account for EBLV-1a and -1b segregation and dispersal in Europe.