- Research article
- Open Access
Evidence for genetic variation in Natterer’s bats (Myotis nattereri) across three regions in Germany but no evidence for co-variation with their associated astroviruses
BMC Evolutionary Biology volume 17, Article number: 5 (2017)
As bats have recently been described to harbor many different viruses, several studies have investigated the genetic co-variation between viruses and different bat species. However, little is known about the genetic co-variation of viruses and different populations of the same bat species, although such information is needed for an understanding of virus transmission dynamics within a given host species. We hypothesized that if virus transmission between host populations depends on events linked to gene flow in the bats, genetic co-variation should exist between host populations and astroviruses.
We used 19 nuclear and one mitochondrial microsatellite loci to analyze the genetic population structure of the Natterer’s bat (Myotis nattereri) within and among populations at different geographical scales in Germany. Further, we correlated the observed bat population structure to that of partial astrovirus sequences (323–394 nt fragments of the RNA-dependent RNA polymerase gene) obtained from the same bat populations. Our analyses revealed that the studied bat colonies can be grouped into three distinct genetic clusters, corresponding to the three geographic regions sampled. Furthermore, we observed an overall isolation-by-distance pattern, while no significant pattern was observed within a geographic region. Moreover, we found no correlation between the genetic distances among the bat populations and the astrovirus sequences they harbored. Even though high genetic similarity of some of the astrovirus haplotypes found in several different regions was detected, identical astrovirus haplotypes were not shared between different sampled regions.
The genetic population structure of the bat host suggests that mating sites where several local breeding colonies meet act as stepping-stones for gene flow. Identical astrovirus haplotypes were not shared between different sampled regions suggesting that astroviruses are mostly transmitted among host colonies at the local scale. Nevertheless, high genetic similarity of some of the astrovirus haplotypes found in several different regions implies that occasional transmission across regions with subsequent mutations of the virus haplotypes does occur.
With their particular social, ecological, physiological and immunological traits, bats provide unique hosts for many viruses to co-evolve with (e.g. [1, 2]). Indeed, bats are increasingly recognized as reservoirs for a wide range of viruses, some of which carry a zoonotic potential, for example rabies and other viruses of the genus Lyssavirus, SARS-like, MERS-like and other coronaviruses [2–5]. Consequently, several studies have investigated the genetic co-variation between different bat species and their associated viruses [6, 7]. However, much less is known about genetic co-variation of viruses and different populations of the same bat species, although such information is required to gain a better understanding of the transmission dynamics within a given host species (e.g. [8, 9]).
The Natterer’s bat (Myotis nattereri Kuhl 1817 sensu lato ) is a non-migratory vespertilionid bat that is widespread throughout Europe with the exception of the Iberian and Italian peninsula and the South of France . While this species uses underground sites for hibernation during the winter months, it mostly roosts in trees and buildings during summer forming maternity colonies that consist of female bats and their juveniles as well as occasionally some males [12, 13]. Males typically roost either individually or in small groups in the vicinity of the maternity colonies . Male and female Natterer’s bats have been found to exhibit philopatry even though males leave their natal colony but stay in its vicinity . Mating takes place during autumn at swarming sites that are typically up to 50 km away from the summer colony .
Recently, various viruses have been reported to be harbored by M. nattereri , including herpes- , lyssa-  and astroviruses (e.g. ). The Astroviridae form a large family of non-enveloped, positive-sense, single-stranded RNA viruses . Astroviruses are mostly transmitted via the fecal-oral route  and may cause diarrhea in many animal species, including humans . However, the route of transmission in bats has not yet been elucidated. Even though astroviruses have been detected in a variety of species , bats have been hypothesized to be a potential reservoir host in Europe and Asia [7, 23]. The high prevalence and diversity of astroviruses harbored by bats is remarkable [18, 21, 23] and their capability to cross species barriers and become adapted to new hosts, including spill-over to other taxa, has been suggested . Due to the occurrence of astroviruses in animals in close contact to humans, e.g. livestock and also bats using human habitation as roosting, it has been argued that astroviruses should be considered as potential candidates for zoonotic infections (e.g. ). However, almost nothing is known about the transmission of astrovirus among different populations of their bat hosts.
We studied patterns of population genetic structure and dispersal of Myotis nattereri within and among three geographic regions of Germany using both nuclear and mitochondrial microsatellite markers. The population genetic structure of M. nattereri has previously only been investigated in the United Kingdom (UK; ). However, the population genetic structure of bats occurring in the UK may be affected by their insular status, and some important differences between insular and continental populations have been described (e.g. ). Thus, our study adds new important insights into the dispersal behavior of M. nattereri in mainland Europe. Moreover, we investigated for the first time if genetic co-variation occurs between populations of a bat host and its harboured astroviruses. Fischer et al.  reported distinctly higher similarities in astrovirus sequences of samples collected from the same bat species in different geographic localities than between samples from different species sampled at the same locality, whereas different results were obtained in the Czech Republic for some other European bat species . Because our analyses are based on the sequences found by Fischer et al.  in M. nattereri, we assume that, similarly to coronaviruses in Chinese bats , astroviruses are mostly transmitted within Natterer’s bats rather than within the local bat community as a whole. We hypothesized that if virus transmission between host populations were associated with events linked to gene flow in the bats, e.g. mating , genetic co-variation should be detectable between host populations and astroviruses on a larger scale (e.g. between geographic regions), even though not necessarily within a certain region.
Many bat astroviruses form distinct phylogenetic clusters , but little is known whether astrovirus population structure matches that of their host species. As astroviruses are hypothesized to be transmitted via the fecal-oral route , both direct and indirect virus transmission within the breeding colonies of bats and at swarming sites during mating seem possible . In comparison, for bat ectoparasites such as bat flies (Nycteribiidae) that are also transmitted both by direct body contact and indirectly through the bats’ roosts, a more efficient transmission among different bat populations has been detected at swarming sites as compared to breeding colonies . As swarming sites represent the main mating sites for M. nattereri, where gene flow occurs, it is expected that if viruses are also mainly transmitted between conspecifics during that period, the transmission route of the viruses should resemble the pattern of host gene flow. Since transmission of astroviruses at swarming sites would cause viruses to be transmitted between members of different colonies visiting a given swarming site, local differentiation should not occur at the colony level but rather on a larger scale. Viral genetic patterns should therefore follow an isolation-by-distance pattern using swarming sites as stepping-stones for gene flow, as suggested for the bat host .
Study area and sample collection
Sampling occurred between 2010 and 2014 from maternity colonies in three regions in Germany: Bavaria (BY, N = 92 adult females, one colony), North Rhine Westphalia (NRW, N = 100 adult females, five colonies) and Mecklenburg Western Pomerania (MV, N = 172 adult females, ten colonies; Fig. 1). Bats were either taken directly out of the bat boxes provided for the colonies or captured using mist nets or harp traps when leaving the colonies’ roosts. Three millimeter wing tissue samples were collected and stored in 90% ethanol until DNA extraction. In addition, between 2012 and 2014, samples of bat saliva, feces and urine were opportunistically taken within the three regions of interest (BY, N = 177; NRW, N = 74; MV, N = 19) and screened for the presence of astrovirus-related RNA (Additional file 1: Table S3).
Bat DNA extraction and amplification
Genomic DNA was extracted using an ammonium acetate precipitation method . Individuals were genotyped using 19 nuclear microsatellite markers and one mitochondrial microsatellite marker. The following nuclear markers were used: EF15 ; b22 ; A2-Mluc, A13-Mluc, E5-Mluc, G6-Mluc, G30-Mluc, G31-Mluc , D15, H19, H23, H29 ; Mnatt-1, Mnatt-2 ; Mschreib3 ; Mnatt-8, Mnatt-11 ; FV5AP  and GZBYR (5′-TCCTTGTCACTATAAGCTCAGTGG-3′ (forward); 5′-CCAGGCAATAGTCTCCTAGCAC-3′ (reverse)). The 5′ end of the reverse primers FV5AP and G30-Mluc were PIG-tailed  with the sequence GTTT and the 5′ end of the reverse primer G31-Mluc with the sequence GTTTT to facilitate adenylation. These 19 autosomal microsatellite markers were amplified in two multiplex polymerase chain reactions (PCRs, Table 1). PCRs were carried out in 9 μl reaction volumes using the Qiagen multiplex Kit (Qiagen, Hilden, Germany). Each multiplex reaction contained 1 x Qiagen Multiplex Master Mix and between 0.11 μM and 1.06 μM of each primer. After drying 1 μl of DNA (approximately 10 ng) for 15 min at 52 °C in a 96-well PCR plate (VWR), multiplex reactions were performed. The PCR amplification was carried out in a 2720 Thermal Cycler (Applied Biosystems), with an initial 5 min denaturation at 95 °C, followed by 30 cycles with denaturation at 95 °C for 30 s, annealing at 60 °C for 90 s and extension at 72 °C for 1 min. Final incubation occurred at 60 °C for 30 min.
The mitochondrial DNA marker AT-2  was amplified in a separate PCR. After drying 1 μl of DNA (approximately 10 ng) for 15 min at 52 °C in a 96-well PCR plate (VWR), this PCR was carried out in 10 μl reaction volume which contained 0.2 μM of Primer AT-2, 1.0 mM dNTPs, 0.8 mM of MgCl2, 1.0 μl 10 x Taq Buffer B1 (Solis, BioDyne, Tartu, Estonia) and 1 unit of Taq Hot FirePol® DNA Polymerase (Solis, BioDyne, Tartu, Estonia). This PCR amplification was carried out in a 2720 Thermal Cycler (Applied Biosystems), with an initial 15 min denaturation at 95 °C, followed by 40 cycles with denaturation at 95 °C for 30 s, annealing at 53 °C for 60 s and extension at 72 °C for 1 min. Final incubation occurred at 72 °C for 7 min.
PCR products were separated using an ABI 3130 Genetic Analyzer (Applied Biosystems) together with the internal size standard Genescan 500 liz (Applied Biosystems) and analyzed using Genemapper v 5.0 (Applied Biosystems).
Bat population genetic structure
Since it has been reported that the presence of closely related individuals within populations can bias Bayesian multi-locus clustering methods , we removed closely related individuals from the dataset before conducting population genetic structure analyses using the program Structure . For this purpose, we first determined the degree of relatedness between all pairs of individuals within a population using TrioML , as implemented in the Coancestry 220.127.116.11 software package . Relatedness densities were further plotted using R (R Core ) in order to determine the relatedness threshold for excluding individuals from the analyses. This threshold was selected manually by best separating the first peak of the plotted distribution (i.e. closely related individuals) from the rest (unrelated individuals). From every pair of individuals with a relatedness value exceeding the determined threshold (0.3), one individual was randomly removed, respectively. All other analyses, except for Structure, were performed on the whole data set. As uneven sampling can bias inferences on the number of clusters in the program Structure , efforts were made to have comparable number of individuals from the three regions investigated after the removal of closely related individuals (71 from BY, 73 from NRW and 104 from MV; cf. Results).
As preliminary runs using the original Structure model showed limited population structure (Additional file 1: Figure S2), Structure  was run on the nuclear DNA dataset assuming admixture and correlated allele frequencies using the LOCPRIOR model that allows for the use of sample group information (here the colony) in the clustering process . Thus, twenty independent runs of K = 1–10 were conducted for the whole dataset after removing closely related individuals as mentioned above. Additionally, twenty independent runs of K = 1–10 were run for each dataset within a sampling region (NRW, MV and BY), respectively. All runs used 106 iterations after a burn-in period of 105.
For each of the genetically distinct populations inferred by Structure the significance of deviations from Hardy-Weinberg equilibrium (HWG, ) and linkage disequilibrium between loci was tested in Genepop 4.0.7 . The False Discovery Rate (FDR) correction method was used to deal with multiple testing .
To assess the level of genetic diversity, the observed (H o ) and expected heterozygosity (H ex ) for each locus as well as for each population inferred by Structure for the complete data set were calculated using Genetix 4.05.2 . The mean number of alleles (A) and the allelic richness (A R ) were calculated for each locus and each subpopulation using Fstat v.2.9.3  (Table 1). Differences in the number of alleles per locus, allelic richness and expected and observed heterozygosities between the populations inferred by Structure were tested for significance using the Wilcoxon signed-rank test in R. MICRO-CHECKER 2.2.3 , set for 10 000 iterations and a 95% confidence interval, was used to test for null alleles.
Population pairwise F ST values  on the whole dataset including closely related individuals, were used to measure the level of genetic differentiation between the populations inferred by Structure. For both the nuclear and the mitochondrial DNA data set, a hierarchical population structure was assumed where colonies were clustered within regions. F ST values were thus calculated using hierarchical analyses within the R-package HierFstat [56, 57]. Genetic structure was tested among colonies within sampling regions with more than 1 colony (MV and NRW) and among sampling regions. The significance of the F-statistics was tested by 10,000 permutations.
Isolation-by-distance for the entire set of 17 colonies was tested via Mantel’s test  from the comparison of all pairwise F ST/(1-F ST) values with pairwise geographic distances using the R package ecodist  with 10,000 permutations. The test was performed for the whole dataset as well as within regions consisting of several colonies (NRW and MV).
Correlation of bat and virus genetic distances
A total of 270 samples obtained from saliva, feces and urine of Myotis nattereri were screened for the presence of astroviruses by using a published hemi-nested PCR protocol  for the amplification of a highly conserved region of the RNA dependent RNA polymerase gene (RdRp). This PCR assay enables the detection of a variety of bat-associated astroviruses by using degenerated primers. Further details of the virus-related sampling and laboratory are presented in Fischer et al. . A total 57 sequences representing 19 different astrovirus haplotypes (=sequences having 100% identity) (N1-19, Fig. 4) were identified from the screening of 270 samples. Three haplotypes were excluded from the analyses (N17-19, Fig. 4) as they were phylogenetically too distant from the remaining 16 and too rarely encountered (only once per haplotype) to make biologically meaningful inferences about their correlation to the bat host’s population genetic structure. Out of the 270 samples, Fischer et al.  were able to assign 73 individual sequences to a specific haplotype (N1-19, Fig. 4), whereas the remaining astrovirus positive samples were shorter than 279 nt and could therefore not be assigned to a haplotype. As in this study we only analyzed individually marked adult female bats that were clearly identifiable via their RFID-tag, in total 57 individual sequences, representing these 16 different haplotypes, were used (N1-16, Fig. 4).
The genetic distances between the different astrovirus haplotypes were calculated using the Tamura-Nei model implemented in Mega 6.0 . Further, the sequences were translated into amino acids and amino acid genetic distances were computed using the p-dist method implemented in Mega 6.0 . This latter measure was tentatively used to differentiate viral species following the recommendation of the Astroviridae working group of the International Committee for Taxonomy of Viruses (ICTV Astroviridae Working ), proposing that mean amino acid genetic distances (p-dist) of the full length ORF2 larger than 0.284 distinguish different species.
As we found no shared identical sequences between regions (BY, MV and NRW), a permutation test was conducted using R (R Core ) to test how likely such a pattern would be expected by chance (Additional file 1: Database S1).
Moreover, the program Poptree2  was used to generate genetic distance matrices for the nuclear and mitochondrial DNA datasets based on the DA distances  of population microsatellite allele frequencies within a priori populations. For the astrovirus dataset, we first used the program jModeltest 2.1.7  to calculate likelihood scores for substitution model selection. Genetic and amino acid distances were then calculated using Mega 6.0  based on the Maximum Composite Likelihood substitution model with gamma correction for among-site rate heterogeneity and an estimated proportion of invariable sites.
Associations between astrovirus and bat host genetic distances (both mitochondrial and nuclear) were first analyzed via a Mantel Test  using the software GenAlex 6.501 . To control for the possible effect of geographic distance, we also carried out partial Mantel tests in PASSaGE v 18.104.22.168 (9999 permutations, ). The geographic distance matrix used was calculated from the GPS coordinates of the different sampling sites using the Create option in PASSaGE.
For the genetic correlations between host and astrovirus sequences the NRW dataset was modified. Due to significant population genetic structuring on the basis of mtDNA and the fact that the sampled colonies in NRW are up to 46 km apart from each other compared to maximally 6.5 km in MV (see Fig. 1) we split the NRW data set in four separate sampling units (NRW2, NRW3, NRW5 and NRW6; Fig. 1). Here, we only used genetic data from the colonies within NRW where virus sequences were detected. Together with MV and BY, the total dataset for comparing host and virus population structures now consisted of six populations between which pairwise genetic distances were computed as mentioned above.
Bat population genetic structure
After removing closely related individuals at the threshold of 0.3, the nuclear DNA dataset consisted of samples from 248 bats in total, including 71 from BY, 73 from NRW and 104 from MV. Using the 19 autosomal microsatellite loci, Structure inferred the presence of three distinct genetic clusters (Fig. 2; Additional file 1: Figure S1), splitting our data set into the three sampled regions NRW, MV and BY. No additional sub-structuring was detected by Structure within any of these three sampling regions (data not shown).
The 19 autosomal microsatellite loci had between 4 and 24 alleles and an average of 10.2–11.7 alleles per each of the three genetically distinct populations inferred by the program Structure (Table 1). No significant deviations from Hardy–Weinberg equilibrium were detected in these three populations. Deviations from linkage disequilibrium at the α = 0.05 level (after FDR correction) did occur consistently in all three regions between two loci (Mnatt-2 and Mnatt-11, Table 1) but were not detected in all respective colonies (linkage occurred in BY, in three out of ten colonies in MV and in three out of five colonies in NRW). Because of the inconsistencies at the colony level, we decided to nevertheless keep both loci for further analyses. No marker with consistently appearing null alleles was found within NRW and MV, whereas in BY the two loci EF15 and FV5AP showed the presence of null alleles. However, since the estimated frequency of null alleles per locus was low (<0.1), we kept those loci.
Mean expected and observed heterozygosity were globally similar across colonies and regions (Tables 1 and 2). The number of alleles found in MV was significantly higher than in BY (P = 0.0020) and NRW (P = 0.0284), but BY and NRW did not differ significantly from each other. Moreover, allelic richness was significantly higher in MV as compared to BY (P = 0.0124) but no significant differences were found between BY and NRW and between NRW and MV. The overall level of differentiation among the three regions was weak but significant (F ST = 0.0088, P = 0.0001) based on nuclear DNA estimated using hierarchical F-statistics. Significant genetic differentiation was further identified among colonies within regions (F ST = 0.0194, P = 0.0001). When analyzing the data obtained by the mitochondrial DNA marker AT-2, pairwise genetic differentiation was found to be much higher at both geographic scales (among regions: F ST = 0.4979, P = 0.0001; among colonies within regions: F ST = 0.3657, P = 0.0001).
When correlating genetic differentiation F ST/(1-F ST) with ln of geographic distances between all colonies within our study area, a significant pattern of isolation-by-distance was detected (Mantel r = 0.2989, P = 0.0160; Fig. 3). In contrast, the Isolation-by-distance patterns within the NRW (Mantel r = -0.2762, P = 0.7869) and the MV population (Mantel r = 0.2275, P = 0.0535) were not significant, even though in MV there was a similar trend visible as in the entire data set.
Correlation of bat and virus genetic distances
We found no overlap in the detected haplotypes between the three regions (Fig. 4). The permutation test revealed that the probability of having no such overlap by chance was very low (P = 0.007; 100,000 permutations). However, a closer look at the genetic relationship of the different astrovirus haplotypes (Fig. 4) revealed that sequences of high similarity occur in different geographic regions. Genetic distances in astrovirus haplotypes that cluster together (N1-3, N4-6, N7-9, N10-11, N12-15; Fig. 4) ranged between 0.0033 and 0.0313, whereas genetic distances among clusters were distinctly higher (0.2093–0.6895). As for the amino acid genetic distances (p-dist), distances ranged between 0.0221–0.0588 and 0.5732–0.6786 within and between clusters, respectively, the latter being typical of species differences. Rough estimations of divergence times using an astrovirus mutation rate of 3 × 10−3  indicate that differences within clusters have occurred within 1–10 years, versus 70–230 years for between-cluster divergence time.
Furthermore, the astrovirus genetic distances neither correlated with those of the host’s nuclear DNA nor with those of the host’s mitochondrial DNA. This was true both for the Mantel Tests and for the partial Mantel Test correcting for geographic distance (Table 3).
In this study, we analyzed the population genetic structure of M. nattereri within and among summer colonies at different geographical scales in Germany and correlated it with that of the astrovirus sequences found in the respective host colonies. Our aim was to assess whether the population structure and dispersal patterns of host populations can explain the genetic structure of astrovirus sequences and thus ultimately shed light on the virus transmission dynamics within a given bat species [8, 9].
Our findings show significant population structure in M. nattereri with the detection of three genetic clusters which correspond to the three regions sampled (NRW, MV and BY). Both the results obtained using the Bayesian clustering approach and the F ST values between the three genetic clusters show the existence of significant population genetic structure. The detection of an isolation-by-distance pattern over the whole study area combined with the continuous presence of the species across Germany suggests that the levels of gene flow are not high enough to prevent some population differentiation .
The observed strong mitochondrial substructure and weak but significant structure at the nuclear level, likely results from differences in effective population size and female philopatry combined with male-mediated gene flow. Differences in patterns of genetic structure in mitochondrial and nuclear DNA have been reported in many European bat species and were associated with male-biased dispersal [41, 69]. According to previous studies in M. nattereri in the UK [13, 30], both sexes appear to be highly philopatric to their natal area but visit central swarming sites during autumn for mating. According to Rivers et al. , this pattern results in the same genetic pattern as permanent male dispersal .
In connection with the existing overall population isolation-by-distance pattern detected and the absence of a significant pattern on a local scale, i.e. within a region, we suggest that individuals from different summer colonies meet and mate at swarming sites within each region (e.g. [13, 70]). This would result in gene flow following a stepping-stone model and would both lead to a significant isolation-by-distance over larger scales  and to the absence or a weak signal on a local scale [30, 33]. Within the UK, isolation-by-distance was not detected between summer colonies of M. nattereri unless distances exceeded 100 km . Our results agree with those obtained by Rivers et al.  as we also did not detect significant genetic isolation-by-distance between summer colonies within a given region in Germany (even though there is a trend in MV), but over the study area as a whole. Therefore, we suggest that the isolation-by-distance pattern observed here is generated by swarming sites acting as stepping-stones for gene flow.
In our study area, identical astrovirus haplotypes harbored by M. nattereri do not overlap between geographic regions (BY, MV and NRW, Fig. 4). Based on the partial sequence of the conserved RdRp gene analysed in this study , where the mean amino acid distances ranged from 0.022 to 0.059 within and 0.573–0.679 between clusters, the analyzed haplotypes form six different groups which might represent six different viral species  (C1-6, Fig. 4). The detection of the same putative viral species in the different regions combined with the estimated divergence times (<10 years for within species) suggests occasional transmission between host populations. Both the observed genetic population structure of M. nattereri and existing data from ringing studies in Germany  show that Natterer’s bats rarely if ever move over long distances (more than 100 km). Thus, it is unlikely that individual bats directly transmit a certain virus haplotype between the three geographic regions analyzed (NRW, BY and MV). Since viruses have a considerably higher mutation rate  compared to the bat host, mutations in the virus sequences do occur within much shorter time scales than mutation in the bats’ genome [74, 75]. As bats need to transmit the astroviruses directly, the movement of viruses across the landscape should mirror the movement of the bats and hence occur successively over large distances following the stepping stone model of the host. As a consequence, virus transmission over large distances is likely to take multiple years. In the course of these successive transmissions events, mutations will occur in the virus which will lead to viral population differentiation as we observe within the putative viral species.
No association was found when correlating the genetic distances of the different astrovirus sequences with their bat host genetic distances. We had originally expected that if the transmission route of astrovirus sequences resembles the pattern of host gene flow, genetic co-variation between astroviruses and host populations should be detected, especially across regions. The reason why no such correlation was found for M. nattereri and its associated astroviruses could be due to strong differences in population size (hence genetic drift) and mutation rate between bats and viruses. In contrast to higher Eukaryotes such as the bat host, RNA viruses are subject to higher selective pressures and combined with a high mutation rate allow continuous and rapid adaptation to changing environmental conditions [73, 76, 77]. Coupled with large population sizes, virus evolution can thus already be observed within very short time scales of weeks to months [74, 75]. The frequent fluctuations in the prevalence of viral populations (e.g. bottlenecks) and hence the higher genetic drift they face might prevent these populations from showing patterns of isolation-by-distance . Moreover, our virus sampling could only be performed during the summer period since only some of the autumn swarming sites are known so far. Additionally, due to logistic reasons exact sampling dates differed between the sampling localities, which could blur the signal if some haplotypes or putative viral species are more abundant in different periods (see Additional file 1: Table S1). Finally, although large from a virology perspective, the number of samples with viral material was relatively limited to perform population genetics analyses. This combined with the variations in sampling times could confound our analysis if viruses show quick temporal variation in prevalence and/or turnover. A further possible explanation for the lack of a genetic correlation between hosts and viruses is that at swarming sites bats may not only interact in ways that lead to gene flow. Multiple mating and contact between individuals that does not result in successful mating might also represent transmission opportunities for viruses that are not reflected within the host genetic pattern.
In summary, our findings suggest that for M. nattereri within Germany, swarming sites act as stepping-stones for gene flow, as indicated by an overall isolation-by-distance pattern and the absence of such a significant pattern on a local scale. The observed population genetic structure indicates that no apparent strong barriers to gene flow exist within our study area for the bat host. While putative viral species were mostly shared between geographic regions, no haplotypes were shared for any putative viral species. Despite the observed genetic differentiation between the three geographic regions in the bat host and to a certain extent also in the harbored astroviruses, we did not detect a correlation between host and virus genetic distances. This could potentially be due to differences in genetic drift, selective pressure, population size and mutation rate between bats and viruses. Further studies with a higher astrovirus sample size and with specific simultaneous sampling during autumn mating at swarming sites are required to shed further light on the host-virus relationship between bats and their astroviruses.
Wilkinson GS, South JM. Life history, ecology and longevity in bats. Aging Cell. 2002;1(2):124–31.
Han HJ, Wen HL, Zhou CM, Chen FF, Luo LM, Liu JW, Yu XJ. Bats as reservoirs of severe emerging infectious diseases. Virus Res. 2015;205:1–6.
Freuling CM, Binger T, Beer M, Adu-Sarkodie Y, Schatz J, Fischer M, Hanke D, Hoffmann B, Hoper D, Mettenleiter TC, et al. Lagos bat virus transmission in an Eidolon helvum bat colony, Ghana. Virus Res. 2015;210:42–5.
Calisher CH, Childs JE, Field HE, Holmes KV, Schountz T. Bats: Important reservoir hosts of emerging viruses. Clin Microbiol Rev. 2006;19(3):531–45.
Hu B, Ge X, Wang L-F, Shi Z. Bat origin of human coronaviruses. Virol J. 2015;12(1):221.
Streicker DG, Turmelle AS, Vonhof MJ, Kuzmin IV, McCracken GF, Rupprecht CE. Host phylogeny constrains cross-species emergence and establishment of rabies virus in bats. Science. 2010;329(5992):676–9.
Zhu HC, Chu DKW, Liu W, Dong BO, Zhang SY, Zhang JX, Li LF, Vijaykrishna D, Smith GJD, Chen HL, et al. Detection of diverse astroviruses from bats in China. J Gen Virol. 2009;90:883–7.
Atterby H, Aegerter JN, Smith GC, Conyers CM, Allnutt TR, Ruedi M, MacNicoll AD. Population genetic structure of the Daubenton’s bat (Myotis daubentonii) in western Europe and the associated occurrence of rabies. Eur J Wildl Res. 2010;56(1):67–81.
Ostfeld RS, Glass GE, Keesing F. Spatial epidemiology: an emerging (or re-emerging) discipline. Trends Ecol Evol. 2005;20(6):328–36.
Puechmaille SJ, Allegrini B, Boston ESM, Dubourg-Savage MJ, Evin A, Knochel A, Le Bris Y, Lecoq V, Lemaire M, Rist D, et al. Genetic analyses reveal further cryptic lineages within the Myotis nattereri species complex. Mamm Biol. 2012;77(3):224–8.
Dietz C, Kiefer A: Die Fledermäuse Europas: kennen, bestimmen, schützen: Kosmos; 2014
Altringham JD. British Bats (Collins New Naturalist Library, Book 93), vol. 93. UK: HarperCollins; 2014.
Rivers NM, Butlin RK, Altringham JD. Autumn swarming behaviour of Natterer’s bats in the UK: Population size, catchment area and dispersal. Biol Conserv. 2006;127(2):215–26.
Smith PG: Habitat preference, range use and roosting ecology of Natterer’s bats (Myotis nattereri) in a grassland-woodland landscape. University of Aberdeen; 2001
Kohl C, Kurth A. European bats as carriers of viruses with zoonotic potential. Viruses. 2014;6(8):3110–28.
Wibbelt G, Kurth A, Yasmum N, Bannert M, Nagel S, Nitsche A, Ehlers B. Discovery of herpesviruses in bats. J Gen Virol. 2007;88:2651–5.
Garcia-Perez R, Ibanez C, Godinez JM, Arechiga N, Garin I, Perez-Suarez G, de Paz O, Juste J, Echevarria JE, Bravo IG. Novel Papillomaviruses in Free-Ranging Iberian Bats: No Virus-Host Co-evolution, No Strict Host Specificity, and Hints for Recombination. Genome Biol Evol. 2014;6(1):94–104.
Fischer K, Zeus V, Kwasnitschka L, Kerth G, Haase M, Groschup MH, Balkema-Buschmann A. Insectivorous bats carry host specific astroviruses and coronaviruses across different regions in Germany. Infect Genet Evol. 2016;37:108–16.
Schultz-Cherry S: Astrovirus Research: Essential Ideas, Everyday Impacts, Future Directions: Springer Science & Business Media; 2012
Mendenhall IH, Smith GJD, Vijaykrishna D. Ecological Drivers of Virus Evolution: Astrovirus as a Case Study. J Virol. 2015;89(14):6978–81.
Chu DKW, Poon LLM, Guan Y, Peiris JSM. Novel astroviruses in insectivorous bats. J Virol. 2008;82(18):9107–14.
Moser LA, Schultz-Cherry S. Pathogenesis of astrovirus infection. Viral Immunol. 2005;18(1):4–10.
Drexler JF, Corman VM, Wegner T, Tateno AF, Zerbinati RM, Gloza-Rausch F, Seebens A, Muller MA, Drosten C. Amplification of Emerging Viruses in a Bat Colony. Emerg Infect Dis. 2011;17(3):449–56.
De Benedictis P, Schultz-Cherry S, Burnham A, Cattoli G. Astrovirus infections in humans and animals - Molecular biology, genetic diversity, and interspecies transmissions. Infect Genet Evol. 2011;11(7):1529–44.
Durrant CJ, Beebee TJC, Greenaway F, Hill DA. Evidence of recent population bottlenecks and inbreeding in British populations of Bechstein’s bat, Myotis bechsteinii. Conserv Genet. 2009;10(2):489–96.
Dufkova L, Strakova P, Sirmarova J, Salat J, Moutelikova R, Chrudimsky T, Bartonicka T, Nowotny N, Ruzek D. Detection of Diverse Novel Bat Astrovirus Sequences in the Czech Republic. Vector Borne Zoonotic Dis. 2015;15(8):518–21.
Tang X, Zhang J, Zhang S, Wang P, Fan X, Li L, Li G, Dong B, Liu W, Cheung C. Prevalence and genetic diversity of coronaviruses in bats from China. J Virol. 2006;80(15):7481–90.
Criscione CD. Parasite co-structure: Broad and local scale approaches. Parasite. 2008;15(3):439–43.
van Schaik J, Dekeukeleire D, Kerth G. Host and parasite life history interplay to yield divergent population genetic structures in two ectoparasites living on the same bat species. Mol Ecol. 2015;24(10):2324–35.
Rivers NM, Butlin RK, Altringham JD. Genetic population structure of Natterer’s bats explained by mating at swarming sites and philopatry. Mol Ecol. 2005;14(14):4299–312.
Nicholls J, Double M, Rowell D, Magrath R. The evolution of cooperative and pair breeding in thornbills Acanthiza (Pardalotidae). J Avian Biol. 2000;31(2):165–76.
Vonhof MJ, Davis CS, Fenton MB, Strobeck C. Characterization of dinucleotide microsatellite loci in big brown bats (Eptesicus fuscus), and their use in other North American vespertilionid bats. Mol Ecol Notes. 2002;2(2):167–9.
Kerth G, Safi K, Konig B. Mean colony relatedness is a poor predictor of colony structure and female philopatry in the communally breeding Bechstein’s bat (Myotis bechsteinii). Behav Ecol Sociobiol. 2002;52(3):203–10.
Jan C, Dawson DA, Altringham JD, Burke T, Butlin RK. Development of conserved microsatellite markers of high cross-species utility in bat species (Vespertilionidae, Chiroptera, Mammalia). Mol Ecol Resour. 2012;12(3):532–48.
Castella V, Ruedi M. Characterization of highly variable microsatellite loci in the bat Myotis myotis (Chiroptera : Vespertilionidae). Mol Ecol. 2000;9(7):1000–2.
Boston ESM, Puechmaille SJ, Scott DD, Buckley DJ, Lundy MG, Montgomery IW, Prodoehl PA, Teeling EC. Empirical assessment of non-invasive population genetics in bats: comparison of DNA quality from faecal and tissue samples. Acta Chiropterologica. 2012;14(1):45–52.
Miller-Butterworth CM, Jacobs DS, Harley EH. Isolation and characterization of highly polymorphic microsatellite loci in Schreibers’ long-fingered bat, Miniopterus schreibersii (Chiroptera : Vespertilionidae). Mol Ecol Notes. 2002;2(2):139–41.
Scott DD, Fitzpatrick S, Bailie DA, Boston ESM, Lundy MG, Buckley DJ, Teeling EC, Montgomery WI, Prodoehl PA. Isolation and characterization of eight polymorphic microsatellite loci for Natterer’s bat, Myotis nattereri (Vespertilionidae, Chiroptera). Conserv Genet Resour. 2013;5(3):643–5.
O’Donnell CF, Richter S, Dool S, Monks JM, Kerth G: Genetic diversity is maintained in the endangered New Zealand long-tailed bat (Chalinolobus tuberculatus) despite a closed social structure and regular population crashes. Conserv Genet. 2015;17(1):1-12.
Brownstein MJ, Carpten JD, Smith JR. Modulation of non-templated nucleotide addition by taq DNA polymerase: Primer modifications that facilitate genotyping. Biotechniques. 1996;20(6):1004–6. 1008-1010.
Kerth G, Mayer F, Petit E. Extreme sex-biased dispersal in the communally breeding, nonmigratory Bechstein’s bat (Myotis bechsteinii). Mol Ecol. 2002;11(8):1491–8.
Rodríguez-Ramilo ST, Wang J. The effect of close relatives on unsupervised Bayesian clustering algorithms in population genetic structure analysis. Mol Ecol Resour. 2012;12(5):873–84.
Pritchard J, Stephens M, Donnelly P. Inference of population structure using multilocus genotype data. Genetics. 2000;155(2):945–59.
Wang J. Triadic IBD coefficients and applications to estimating pairwise relatedness. Genet Res. 2007;89(03):135–53.
Wang J. COANCESTRY: a program for simulating, estimating and analysing relatedness and inbreeding coefficients. Mol Ecol Resour. 2011;11(1):141–5.
Team RC: R: A language and environment for statistical computing. R Foundation for Statistical Computing, Vienna, Austria, 2012. In.: ISBN 3-900051-07-0; 2013
Puechmaille SJ. The program structure does not reliably recover the correct population structure when sampling is uneven: subsampling and new estimators alleviate the problem. Mol Ecol Resour. 2016;16(3):608–27.
Hubisz MJ, Falush D, Stephens M, Pritchard JK. Inferring weak population structure with the assistance of sample group information. Mol Ecol Resour. 2009;9(5):1322–32.
Nei M. Estimation of average heterozygosity and genetic distance from a small number of individuals. Genetics. 1978;89(3):583–90.
Raymond M, Rousset F: GENEPOP (version-1.2) - Population-genetics software for exact tests and ecumenicism. J Hered. 1995;86(3):248-49.
Benjamini Y, Hochberg Y: Controlling the false discovery rate - A practical and powerful approach to multiple testing. J R Stat Soc Ser B Methodol. 1995;57(1):289-300.
Belkhir K, Borsa P, Chikhi L, Raufaste N, Bonhomme F. GENETIX 4.05, logiciel sous Windows pour la génétique des populations. France: Laboratoire génome populations interactions, CNRS, Université de Montpellier; 2004.
Goudet J: FSTAT (Version 1.2): A computer program to calculate F-statistics. J Hered. 1995;86(6):485-6.
Van Oosterhout C, Hutchinson WF, Wills DPM, Shipley P. MICRO-CHECKER: software for identifying and correcting genotyping errors in microsatellite data. Mol Ecol Notes. 2004;4(3):535–8.
Wright S. The genetical structure of populations. Ann Eugenics. 1951;15(4):323–54.
Goudet J. HIERFSTAT, a package for R to compute and test hierarchical F-statistics. Mol Ecol Notes. 2005;5(1):184–6.
De Meeûs T, Goudet J. A step-by-step tutorial to use HierFstat to analyse populations hierarchically structured at multiple levels. Infect Genet Evol. 2007;7(6):731–5.
Mantel N. The detection of disease clustering and a generalized regression approach. Cancer Res. 1967;27(2):209–20.
Goslee SC, Urban DL. The ecodist package for dissimilarity-based analysis of ecological data. J Stat Softw. 2007;22(7):1–19.
Tamura K, Stecher G, Peterson D, Filipski A, Kumar S. MEGA6: Molecular Evolutionary Genetics Analysis Version 6.0. Mol Biol Evol. 2013;30(12):2725–9.
Group IAW. ICTV Official Taxonomy: Updates since the 8th Report. 2012.
Takezaki N, Nei M, Tamura K. POPTREE2: Software for Constructing Population Trees from Allele Frequency Data and Computing Other Population Statistics with Windows Interface. Mol Biol Evol. 2010;27(4):747–52.
Nei M, Tajima F, Tateno Y. Accuracy of extimated phylogenetic trees from molecular data, II. Gene-Frequency data. J Mol Evol. 1983;19(2):153–70.
Posada D. jModelTest: Phylogenetic model averaging. Mol Biol Evol. 2008;25(7):1253–6.
Peakall R, Smouse PE. GenAlEx 6.5: genetic analysis in Excel. Population genetic software for teaching and research-an update. Bioinformatics. 2012;28(19):2537–9.
Rosenberg MS, Anderson CD. PASSaGE: Pattern Analysis, Spatial Statistics and Geographic Exegesis. Version 2. Methods Ecol Evol. 2011;2(3):229–32.
Babkin IV, Tikunov AY, Zhirakovskaia EV, Netesov SV, Tikunova NV. High evolutionary rate of human astrovirus. Infect Genet Evol. 2012;12(2):435–42.
Frantz AC, Cellina S, Krier A, Schley L, Burke T. Using spatial Bayesian methods to determine the genetic structure of a continuously distributed population: clusters or isolation by distance? J Appl Ecol. 2009;46(2):493-505.
Castella V, Ruedi M, Excoffier L. Contrasted patterns of mitochondrial and nuclear structure among nursery colonies of the bat Myotis myotis. J Evol Biol. 2001;14(5):708–20.
Kerth G, Kiefer A, Trappmann C, Weishaar M. High gene diversity at swarming sites suggest hot spots for gene flow in the endangered Bechstein’s bat. Conserv Genet. 2003;4(4):491–9.
Kimura M, Weiss GH. The stepping stone model of population structure and the decrease of genetic correlation with distance. Genetics. 1964;49(4):561–76.
Steffens R, Zöphel U, Brockmann D, Dresden F. 40 Jahre Fledermausmarkierungszentrale Dresden: methodische Hinweise und Ergebnisübersicht: Saxoprint; 2004
van Hemert FJ, Lukashov VV, Berkhout B. Different rates of (non-)synonymous mutations in astrovirus genes; correlation with gene function. Virol J. 2007;4:12.
Domingo E, Holland JJ. RNA virus mutations and fitness for survival. Annu Rev Microbiol. 1997;51:151–78.
Ojosnegros S, Beerenwinkel N. Models of RNA virus evolution and their roles in vaccine design. Immunome Res. 2010;6 Suppl 2:S5.
Holmes EC. Future issues in RNA virus evolution. Futur Virol. 2006;1(2):243–9.
Drake JW, Charlesworth B, Charlesworth D, Crow JF. Rates of spontaneous mutation. Genetics. 1998;148(4):1667–86.
Hutchison DW, Templeton AR. Correlation of pairwise genetic and geographic distance measures: Inferring the relative influences of gene flow and drift on the distribution of genetic variability. Evolution. 1999;53(6):1898–914.
We thank the local forestry and conservation departments as well as the “Fledermausforschungsprojekt Wooster Teerofen e.V.” for support and the following people for help in the field: Daniela Fleischmann, Lena Grosche, Ralf Koch, Markus Melber, and Jaap van Schaik. For general support during the lab work we thank Ina Römer, and Katja Böhm and Franziska Neitzel who helped to design the Primers FV5AP and GZBYR. We are grateful to Stefan Braendel, Kathleen Drescher, Nicola Fischer, Alain Frantz, Darren Obbard, Jaap van Schaik, Caroline Schöner, Eric Witt and one anonymous referee for providing comments on the manuscript.
This work was funded by the Deutsche Forschungsgemeinschaft (KE 746/6-1 and KE 746/7-1) within the priority programm “Ecology and species barriers in emerging viral diseases (SPP 1596)”.
Availability of data and material
The datasets supporting the conclusions of this article are included within the article (and its Additional files 1 and 2). The GenBank accession numbers for the astrovirus sequences detected to be harboured by M. nattereri are KT894898 - KT894916.
TH and GK designed the study with input from AB-B, SJP and KF. GK, VZ, FM and CT provided samples. TH and RG performed the laboratory work. TH performed data analyses on bats and viral sequences and drafted the manuscript with input from AB-B, KF, SJP and GK. SJP contributed to analysis. All authors have commented on the manuscript, read and approved the final version of the manuscript.
The authors declare that they have no competing interests.
Consent for publication
Ethics approval and consent to participate
Tissue sampling was carried out under the following permits: BY: 55.2-2531.01-79/10 and 55.2-2531.01-47/11; MV: 7221.3-1.1-023/12; NRW: 70.2-0197/08 and 70.2-2012/0254.
Supplementary material. Figure S1. Summary of the log-likelihood values from the 20 independent runs conducted with Structure for the number of genetic clusters (K) set to a minimum of 1 and a maximum of 10. The left graph shows the log-likelihood results of the runs for each K, whereas the right graph shows Delta K plotted against K. The most likely number of genetic clusters is three using both methods. Figure S2. Comparisons of the STRUCTURE runs for K = 3 with (top) or without (bottom) the LOCPRIOR option. Table S1. Table of sampling times and associated screening for astroviruses. AstV = astrovirus, BY = Bavaria, MV = Mecklenburg Western Pomerania, NRW = North Rhine Westphalia. AstV positive samples with a length of less than 279 nt could not be assigned to a specific haplotype, but are nevertheless included in the number of AstV positive samples presented here. Database S1. R-script used for the permutation test. With this script we tested how likely it is to have no overlap in virus haplotypes across regions and whether virus haplotypes are more different between regions than expected by chance. (DOCX 318 kb)
Genotype table Myotis nattereri. Table of alleles per locus for each individual bat used for the population genetic analyses on the Natterer’s bat (Myotis nattereri) determined using Genemapper v 5.0 (Applied Biosystems). (XLSX 75 kb)
About this article
Cite this article
Halczok, T.K., Fischer, K., Gierke, R. et al. Evidence for genetic variation in Natterer’s bats (Myotis nattereri) across three regions in Germany but no evidence for co-variation with their associated astroviruses. BMC Evol Biol 17, 5 (2017) doi:10.1186/s12862-016-0856-0
- Population structure
- Myotis nattereri
- Gene flow