Next Article in Journal
Telomere Homeostasis: Interplay with Magnesium
Previous Article in Journal
Individualizing Treatment Approaches for Epileptic Patients with Glucose Transporter Type1 (GLUT-1) Deficiency
Previous Article in Special Issue
Hepatitis C Virus: Viral Quasispecies and Genotypes
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
IJMS
Volume 19
Issue 1
10.3390/ijms19010156
ijms-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Review

Molecular Epidemiology and Evolution of European Bat Lyssavirus 2

by
Lorraine M. McElhinney
1,2,*,†,
Denise A. Marston
1,†,
Emma L. Wise
1,
Conrad M. Freuling
3,
Hervé Bourhy
4,
Reto Zanoni
5,
Torfinn Moldal
6,
Engbert A. Kooi
7,
Antonie Neubauer-Juric
8,
Tiina Nokireki
9,
Thomas Müller
3 and
Anthony R. Fooks
1,2
1
Wildlife Zoonoses and Vector Borne Disease Research Group, Animal and Plant Health Agency (APHA), Surrey KT15 3NB, UK
2
Institute of Global Health, University of Liverpool, Liverpool L69 3BX, UK
3
Institute of Molecular Virology and Cell Biology, Friedrich Loeffler Institute, (FLI), 17493 Greifswald, Germany
4
Institut Pasteur, Unit Lyssavirus Dynamics and Host Adaptation, 75015 Paris, France
5
Institute of Virology and Immunology, University of Berne, 3012 Berne, Switzerland
6
Norwegian Veterinary Institute, 0454 Oslo, Norway
7
Wageningen Bioveterinary Research, 8221 RA Lelystad, The Netherlands
8
Bavarian Health and Food Safety Authority, 85764 Oberschleissheim, Germany
9
Finnish Food Safety Authority Evira, 00790 Helsinki, Finland
*
Author to whom correspondence should be addressed.
†
These authors contributed equally to this work.
Int. J. Mol. Sci. 2018, 19(1), 156; https://doi.org/10.3390/ijms19010156
Submission received: 29 November 2017 / Revised: 26 December 2017 / Accepted: 29 December 2017 / Published: 5 January 2018
(This article belongs to the Special Issue Virus Comparative Genomics)
Browse Figures
Versions Notes

Abstract

:
Bat rabies cases in Europe are mainly attributed to two lyssaviruses, namely European Bat Lyssavirus 1 (EBLV-1) and European Bat Lyssavirus 2 (EBLV-2). Prior to the death of a bat worker in Finland in 1985, very few bat rabies cases were reported. Enhanced surveillance in the two subsequent years (1986–1987) identified 263 cases (more than a fifth of all reported cases to date). Between 1977 and 2016, 1183 cases of bat rabies were reported, with the vast majority (>97%) being attributed to EBLV-1. In contrast, there have been only 39 suspected cases of EBLV-2, of which 34 have been confirmed by virus typing and presently restricted to just two bat species; Myotis daubentonii and Myotis dasycneme. The limited number of EBLV-2 cases in Europe prompted the establishment of a network of European reference laboratories to collate all available viruses and data. Despite the relatively low number of EBLV-2 cases, a large amount of anomalous data has been published in the scientific literature, which we have here reviewed and clarified. In this review, 29 EBLV-2 full genome sequences have been analysed to further our understanding of the diversity and molecular evolution of EBLV-2 in Europe. Analysis of the 29 complete EBLV-2 genome sequences clearly corroborated geographical relationships with all EBLV-2 sequences clustering at the country level irrespective of the gene studied. Further geographical clustering was also observed at a local level. There are high levels of homogeneity within the EBLV-2 species with nucleotide identities ranging from 95.5–100% and amino acid identities between 98.7% and 100%, despite the widespread distribution of the isolates both geographically and chronologically. The mean substitution rate for EBLV-2 across the five concatenated genes was 1.65 × 10−5, and evolutionary clock analysis confirms the slow evolution of EBLV-2 both between and within countries in Europe. This is further supported by the first detailed EBLV-2 intra-roost genomic analysis whereby a relatively high sequence homogeneity was found across the genomes of three EBLV-2 isolates obtained several years apart (2007, 2008, and 2014) from M. daubentonii at the same site (Stokesay Castle, Shropshire, UK).

Graphical Abstract

1. Introduction

There are currently 14 recognised species in the lyssavirus genus of the family Rhabdoviridae [1]. These species are comprised of the following viruses: rabies virus (RABV), Lagos bat virus (LBV), Mokola virus (MOKV), Duvenhage virus (DUVV), European bat lyssavirus 1 (EBLV-1), European bat lyssavirus 2 (EBLV-2), Australian bat lyssavirus (ABLV), Aravan virus (ARAV), Khujand virus (KHUV), Irkut virus (IRKV), West Caucasian bat virus (WCBV) Shimoni bat virus (SHIBV), Ikoma virus (IKOV), and Bokeloh bat lyssavirus (BBLV). A further two putative bat lyssaviruses await classification: Lleida lyssavirus (LLEBV) [2] and Gannoruwa bat lyssavirus (GBLV) [3]. With the exception of MOKV and IKOV, all lyssaviruses have been isolated from bats.
Cases of bat rabies occur globally, although different bat lyssavirus species are present in different regions [4,5]. Bat rabies cases in Europe are principally attributed to EBLV-1 and EBLV-2. However, three other lyssaviruses have been detected in European insectivorous bats, albeit in very few cases. WCBV was isolated in 2002 from a bent-winged bat (Miniopterus schreibersii) in the Caucasus Mountains and LLEBV from the same bat species in Spain in 2011 [2,6]. Antigenic and phylogenetic profiling of BBLV, isolated from Myotis nattererii in Germany and France, confirms the presence of a fifth lyssavirus in European bats [7,8,9].
The first record of a lyssavirus infected insectivorous bat in Europe was in Hamburg, Germany, in 1954 [10]. The reports described a biting incident from a bat that was caught and tested but not formally identified. Negri bodies were detected following mouse inoculation tests, but no tissue was retained for virus typing, so the bat and virus species involved are unknown. Prior to the death of a bat worker in Finland in 1985 [11], very few bat rabies cases were reported. Interestingly, 23% (263 cases) of all subsequent reported cases were identified through enhanced surveillance in the subsequent two years (1986 and 1987). With the advent of molecular PCR assays (polymerase chain reaction) and subsequent sequencing, the viruses were genetically characterised as two separate lyssavirus species, EBLV-1 and -2 [12].
Between 1977 and 2016, 1183 infected bats were reported in Europe with an average of 34 infected bats reported per year to Rabies Bulletin Europe in the last decade. More than 80% of the reported cases were identified in three countries the Netherlands (33%), Germany (27%), and Denmark (21%). Differences in surveillance systems, bat ecology, public awareness, and available resources may account for the heterogeneous spatial pattern of confirmed bat rabies cases [13]. However, the apparent host specificity may also limit the prevalence of bat lyssaviruses when considering the wide distribution, high frequency, and sympatric habitation of some European insectivorous bat species. Whilst lyssaviruses have the ability to cause clinical disease in bats, they can also circulate within healthy bat populations detectable only via active surveillance screening [14,15].
Cross species transmission (CST) of bat lyssaviruses to other mammals has been documented but is extremely rare [16]. Two CST events involving EBLV-2 from bats to humans have been reported, in Finland in 1985 and in the UK in 2002. Unlike for EBLV-1, EBLV-2 spillover into other wildlife or domestic species has not been documented and no other lyssavirus has yet been isolated from M. Daubentonii or M. dasycneme.
Evolutionary analyses, sub-clinical infections and restricted virus–host relationships suggest that all lyssaviruses, including RABV, most likely originated in bats [16,17,18]. Previous attempts to study the phylogeography, molecular epidemiology, and evolutionary rates for EBLV-2 have been restricted to only a small number of virus isolates or partial genomic sequences [19,20,21]. The historic sharing of viruses between European rabies laboratories for confirmation and research was not standardised. This is particularly evident for the EBLV-2 viruses due to the small number involved. However, differences in virus handling (multiple passages in vivo and in vitro), various sequencing approaches (direct sequencing from amplicons, sequencing from cloned products, and deep sequencing), and improvements in the proof reading fidelity of molecular enzymes have resulted in anomalous published sequences being submitted for the same isolate [20]. To address this, a network of rabies reference laboratories involved in EBLV-2 surveillance and diagnosis was established to review, collate, and analyse all available EBLV-2 submissions. The aim of this review was to collate and analyse a comprehensive and accurate genomic sequence dataset for EBLV-2 providing researchers with both epidemiological and genetic data for every confirmed case. Furthermore, whole genome sequences from the same site were analysed to understand the diversity of viruses within the same maternity roost.

2. Results and Discussion

2.1. Review of Reported Cases

Data from all potential cases of EBLV-2 were collated. The distribution and submission details are summarised in Figure 1 and Table 1 and include confirmed cases of EBLV-2.
EBLV-2 was first isolated in 1985 from a Swiss bat biologist who had been working with bats in Finland, Switzerland, and Malaysia. In 1986, EBLV-2 was isolated in Switzerland from M. daubentonii (Daubenton’s bats). Between 1985 and 2017, a total of 39 cases are suspected to have been associated with EBLV-2, of which 34 have been confirmed (Table 1). Full genome sequences have been previously published or obtained in this study for 29 EBLV-2. In addition, EBLV-2 RNA has been detected by RT-PCR in oro-pharyngeal swabs of two apparently healthy M. daubentonii bats during active surveillance (Scotland and Switzerland), yielding partial N gene sequences (Table 1). Unfortunately, sequence data obtained for three outstanding EBLV-2 cases (UK RV3385, Germany F860/LGL and 2Finland2016) were obtained too late to be included in this analysis. Assuming that rabies in M. daubentonii and M. dasycneme (Pond bats) is restricted to EBLV-2 infection, there were three further probable EBLV-2 cases in Denmark from single M. dasycneme (1986, 1987) and M. daubentonii (1986) bats [22]. In Germany, bat rabies was diagnosed in two M. daubentonii in 1986 [23]. However, no material is available for confirmation of an EBLV-2 infection from these probable cases so they remain unconfirmed. The total number of confirmed (typed) EBLV-2 cases is relatively small compared to EBLV-1 isolated from Eptesicus bat species. This discrepancy could be due to less effective surveillance in non-synanthropic bat species like Daubenton’s bats compared to Eptesicus bats or the low transmission rates of EBLV-2 among conspecifics as indicated by experimental studies in Daubenton’s bats [24].
Fourteen cases of EBLV-2 in M. daubentonii were reported through passive surveillance in the UK between 1996 and 2017 and a single human case in a male in Scotland in 2002 [29]. Unfortunately, sequence data from the latest UK case in September 2017 (Table 1, Derbyshire, England) was obtained too late to be included in the phylogenetic analysis in this review. However, preliminary sequence analysis (N1356) suggests a closer identity (99.6%) to the 1996 EBLV-2 RV628 and 2006 EBLV-2 RV2159 cases (Table 1, Sussex and Oxfordshire, respectively). In addition, EBLV-2 viral RNA was detected in a single oro-pharyngeal swab specimen collected from a healthy Daubenton’s bat sampled as part of an active surveillance programme in Scotland in 2008 [20].
The discovery of EBLV-2 in M. daubentonii in Finland (2009) provided evidence for the possible origin of human infection in 1985 [40]. Prior to this latter case, surveillance had failed to detect EBLV-2 in Finland despite intensive efforts in 1986 by active surveillance [45] and continued passive surveillance since (albeit for low submissions averaging 10/year). Subsequent active surveillance has failed to detect viral RNA in the oro-pharyngeal swabs of healthy Daubenton’s bats [46]. A second case of EBLV-2 in M. daubentonii has more recently been reported in Inkoo, Southern Finland (October 2016) but genome sequence case was not available for this review [44]. However, partial N gene data (N253, accession number MF326269) confirms a closer sequence identity (97.6–98%) to the previous 2 Finland EBLV-2 cases than to any other EBLV-2 (91.7–95.7%).
In Germany, six rabies cases in Daubenton’s bats have been typed as EBLV-2 [33,34,35,41], and five genome sequences were available to be included in this study. Notably, only one case was detected in the frame of routine rabies surveillance, whilst the other five isolates were only found during enhanced surveillance programmes targeting all dead found bats [34]. Genomic sequence data for the outstanding German case (Table 1, Bavaria 2013 860/LGL) were obtained too late to be included in this study, but preliminary sequence analysis (N1336) confirms the closest sequence identity to the 2012 Bavarian German EBLV-2 F500/LGL (99.8%).
The Netherlands is the only country where EBLV-2 was detected in M. dasycneme bats [47], of which sequences are available for four isolates (Table 1). Interestingly, passive surveillance has not yielded any further EBLV-2 cases since 1993.
In Switzerland, passive surveillance between 1976 and 2009 revealed three cases of EBLV-2, all in M. daubentonii and a single PCR positive oro-pharyngeal swab sample via active surveillance [26]. The three Swiss EBLV-2 isolates identified through passive surveillance have previously been analysed [19,20,26,48], but discrepancies in the sequence obtained for the three isolates were noted between different studies. In addition, it is unclear if the partial N gene sequence published as U89479, which shows a number of mismatches compared to all other Swiss EBLV-2 sequences, represents an additional Swiss virus or is a result of erroneous sequencing. In order to solve these issues, original virus isolates for the Swiss EBLV-2 held at APHA and University of Berne were subject to NGS and the correct genome sequences were obtained in this study (Table 1).
The first case of EBLV-2 in Norway was reported in a Daubenton’s bat in 2015. The bat was found by a member of the public in Oppland, in the southern part of Norway [42].

2.2. Phylogenetic Analysis

The phylogenetic analysis confirmed the close relationship of EBLV-2 sequences to the Phylogroup I lyssaviruses, in particular KHUV and BBLV, irrespective of the gene studied (Figure 2a,b, maximum likelihood trees for other genes not shown). Where available, a number of representative lineages for the lyssaviruses were included. The genetic diversity observed for the other lyssavirus species is in contrast to the highly homogeneous sequence identities observed for EBLV-2. However, despite this homogeneity, the EBLV-2 sequences are resolved geographically for all analyses, with clear clustering at the country level. Even within country, resolution can be observed when sufficient numbers allow. The Scotland and Northern England EBLV-2 isolates can be seen to cluster separately to the Southern England and Wales EBLV-2 isolates, although bootstrap support is not consistent (Figure 2 and Figure 3). The human EBLV case (RV1333) in Scotland in 2002, resolved closest to the Northern England and Scotland cases. In the partial N400 analysis (not shown), human UK EBLV-2 was most similar to the partial sequence obtained from a Daubenton’s saliva sample during active surveillance in Scotland [20]. Although originating from disparate regions in the south and north of Germany, the German EBLV-2 viruses cluster together and are most closely associated with the 2002 Swiss isolate (TW0118/02) (Figure 2). All of the Dutch EBLV-2 sequences cluster separately from the other sequences possibly reflecting either geographical and/or host species relationships. The latter may be likely as the EBLV-2 isolates from the Netherlands are the only ones sampled from M. dasycneme. The Norwegian EBLV-2 is distinct but is most closely related to the Swiss EBLV-2.
The two available EBLV-2 genome sequences from Finland are clearly distinct from the other EBLV-2 studied and have previously been proposed as a separate lineage, although depending upon the gene studied, the bootstrap support is not always consistent. However, the relatively low number of available EBLV-2 sequences makes it difficult to confidently further divide into monophyletic sub lineages. Although, the subtypes 2a and 2b of EBLV-2 were proposed in earlier studies [48], it has more recently been reported that the highest nucleotide sequence divergence was between two isolates (AY863408—Switzerland isolate and AY863406—Finland human case), previously assigned as EBLV-2b [49]. The two EBLV-2 sequences from Finland [21] and the erroneous sequence U89479 consistently formed a well-supported sub lineage (EBLV-2b?) irrespective of the model used (bootstrap support 100% for Neighbour Joining (p-values) and 99% for Maximum Likelihood (GTR) models). In contrast, the 1992 & 1993 Swiss EBLV-2 sequences were not consistently associated with either sub-lineage when individual sequences were removed or additional lyssavirus sequences added.
Intra-genotypic comparisons were reported previously for the five viral proteins of EBLV-2 [49] and demonstrated the limited genetic diversity of these lyssaviruses compared to RABV at both the nucleotides and amino acid levels. In these analyses, using concatenated gene sequences, the similarity of the EBLV-2 ranged from 95.5–100% for nucleotides and from 98.7–100% for amino acids.
With respect to the anomalous sequence U89479, no evidence has been found to support it as an additional EBLV-2 sequence, hence we conclude it is erroneous and recommend it is not used in future EBLV-2 analyses.

2.3. Intra-Roost Variation

The three EBLV2 isolates from Daubenton’s bats found at Stokesay castle (Shropshire, UK) were studied in greater detail as they represented the only cluster of EBLV-2 with three separate isolations at the same site in 2007, 2008, and 2014. The percentage of viral reads assembled from the total reads varied between isolates, ranging from 0.98% (RV2473) to 2.08% (RV2336). Depth of coverage (average and maximum) was variable and correlated with the number of viral reads assembled (Table S1).
The three genomes were aligned and found to be 99.9% identical to one another. Positions at which single nucleotide polymorphisms (SNPs) occurred were identified and investigated to determine whether they denote a change in amino acid within a region of the genome coding for a viral protein (Figure 4). Surprisingly, only one such SNP was identified within the P gene at Position 2015, despite up to a 7-year interval between isolations. RV2336 and RV2974 code for lysine at this position, whereas the change in base from A to G in RV2473 results in the translation of arginine instead (Table S2). This SNP is located within the LC8 domain (cytoplasmic dynein light chain) at residue 164. The LC8 domain has been previously identified as a region at which the P protein interacts with LC8, a cytoplasmic protein associated with axonal transport of RABV within neurons [49,50], but the SNP is not within the key residues of 145–149, the motif to which LC8 binds. In addition, the substitution is considered conservative due to the similar biochemical properties of the amino acids.
To investigate heterogeneity of the 3 Stokesay Castle isolates across the length of each genome, each position in the genome was analysed to determine the number of reads showing each Base A, T, C, or G. The proportion of reads showing each base was calculated and plotted to give a graphical representation of heterogeneity across the genome for each of the three EBLV-2 isolates. Data was filtered to remove positions that showed less than 100 reads, and less than or equal to 1% heterogeneity. Across the length of the genome, RV2336 has 23 positions at which heterogeneity equals or exceeds 1%, RV2473 has 257 positions and RV2974 has 212. The similar level of heterogeneity between RV2473 and RV2974 gives us confidence that this data is not merely a reflection of the varying depth of coverage between the genomes; RV2473 has the lowest average depth of coverage (245 reads) whilst RV2974 has the highest (970 reads). There are three positions at which heterogeneity equals or exceeds 2% in RV2336, 50 positions in RV2473, and 86 positions in RV2974. Additional analysis and figures can be found in the Supplementary data, Supplementary Tables S1 and S2 and Supplementary Figure S1.

2.4. Evolutionary Analysis of EBLV-2

Evolutionary analysis was undertaken using an alignment of concatenated genes for the 29 EBLV-2 sequences. The analysis supports the previously observed relatively low evolutionary rates for EBLV-2 compared to other RNA viruses [19,21]. The ancestral root of the available EBLV-2 concatenated dataset was likely to have diverged as far back as 436AD, but the 95% HPD intervals are wide (−9553 to 1655). The recent ancestors of the geographical clusters of EBLV-2 are estimated as follows (Figure 3); Finland in 1880 (range 1133 to1969), Norway in 706 (single isolate −7535 to 1664); Switzerland in 890 (−6280 to 1734), Germany in 1567 (−1326 to 1910); Netherlands in 1862 (1042 to 1961); UK in 1639 (−727 to 1923). The Stokesay EBLV-2 ancestors are estimated to have existed much more recently in 1970 (1725 to 2003). The long TMRCAs observed (Figure 3) is suggestive of minimal and discrete evolution at the local level (Figure 3). Migration of infected bats, whilst possible, does not appear to be supported in our dataset. Mean rate of evolution for the concatenated genes for the EBLV-2 dataset (n = 29) was 1.7 × 10−5 (4.65 × 10−8, 3.63 × 10−5).

3. Materials and Methods

Firstly, a literature review was undertaken of all reported cases of EBLV-2 between 1985 and 2017. As well as peer-reviewed literature and published sequences on GenBank, other sources included Rabies Bulletin Europe Reports, ProMed reports, and reference laboratory archives. When available and where gaps existed, original brain material was sought to clarify anomalies in published data and obtain complete genome sequences (Table 1).

3.1. Next Generation Sequencing (NGS)

Total RNA (from original brain material) was depleted of host genomic DNA (gDNA) and ribosomal RNA (rRNA) following methods described previously [51,52]. Briefly, gDNA was depleted using the on-column DNase digestion protocol in RNeasy plus mini kit (Qiagen, Manchester, UK) following manufacturer’s instructions, eluting in 30 µL in molecular grade water. Subsequently, rRNA was depleted, using Terminator 5′-phosphate-dependent exonuclease (CamBio, Cambridge, UK). Briefly, 30 µL of gDNA depleted RNA was mixed with 3 µL of Buffer A, 0.5 µL of RNAsin Ribonuclease inhibitor (20–40 U/µL) (Promega, Southampton, UK), and incubated at 30 °C for 60 min. The depleted RNA was purified to remove the enzyme using the RNeasy plus mini kit as above, without the DNase digestion, eluting in 30 µL of molecular grade water. ds cDNA (double-stranded copy DNA) was synthesised using random hexamers and a cDNA synthesis kit (Roche, Burgess Hill, UK) following manufacturer’s instructions. The resulting ds-cDNA was purified using AMPure XP magnetic beads (Beckman Coulter, High Wycombe, UK), quantified using Quantifluor (Promega) and approximately 1 ng of each sample was used in a “tagmentation” reaction mix using a Nextera XT DNA sample preparation kit (Illumina, San Diago, CA, USA) following the manufacturer’s instructions—without the bead normalization step. DNA libraries were quantified using Quantifluor (Promega). Individual libraries were pooled (approximately 8–12 samples per MiSeq run) and normalized to equimolar concentrations. Libraries were sequenced as 2 × 150 bp paired-end reads on an Illumina MiSeq platform (Illumina).

3.2. Mapping to Reference Sequence to Obtain Complete Genome Sequences

Reads were mapped using the Burrow–Wheeler Aligner (BWA version 0.7.5a–r405) and were visualized in Tablet [53]. A modified SAMtools/vcfutils [54] script was used to generate an intermediate consensus sequence in which any indels and SNPs relative to the original reference sequence were appropriately called. The intermediate consensus sequence was used as the reference for four subsequent iterations of mapping and consensus calling. Sequencing resulted in 100% coverage of the five genes for all samples. Consensus sequences were submitted to Genbank (accession numbers provided in Table 1).

3.3. Evolutionary Analysis

Sequences generated in this study were combined with published EBLV-2 sequences and aligned using ClustalX2 (version 1.2, www.clustal.org/clustal2). Nucleotide and amino acid identities were obtained using MegAlign, v11.1.0 (Lasergene DNASTAR. Madison, WI, USA).
A phylogenetic analysis of all 29 EBLV-2 concatenated coding gene sequences (Accession numbers provided in Table 1) was undertaken using Bayesian Markov Chain Monte Carlo simulation in the BEAST package v1.8.0 [55]. A GTR nucleotide substitution model with rate variation and a proportion of invariant sites (Gamma) were determined to best fit the data using Akaike Information Criterion (AIC) in MEGA 6.0 [56]. A relaxed molecular clock model with constant population prior were used in Markov Chain Monte Carlo (MCMC) simulation for 100,000,000 iterations, sampling every 10,000 states to give effective sample sizes of over 200. Molecular clock and population coalescent models were compared using a modified AIC (AICM) in Tracer v1.5, as described previously [57], and maximum clade credibility trees were annotated using TreeAnnotator (v1.8.4) after 10% of trees (burn-in) were discarded [55]. The maximum clade credibility (MCC) trees were then visualised using Fig Tree (v1.4.3, http://tree.bio.ed.ac.uk/software/figtree/). Maximum likelihood phylogenetic reconstructions were also undertaken on the EBLV-2 dataset (n = 29) in MEGA 6.0 [56] for each individual gene and for partial N gene sequences (400 nucleotides).

4. Conclusions

This review represents the most comprehensive and accurate genomic analysis of EBLV-2 to date and includes the first detailed intra-roost genetic analysis of EBLV-2 isolated from the same site between 2007 and 2014.
The chiropteran host range of confirmed EBLV-2 appears to be restricted to two bat species M. Daubentonii (n = 28) and M. dasycneme (n = 5) (Table 1). The distribution of M. dasycneme is confined to Central and Eastern Europe, whereas M. daubentonii is widely distributed throughout Europe [58]. M. dasycneme is considered to be absent in the UK. However, given the relationship between bat species and their lyssaviruses [4], it remains unclear as to whether M. dasycneme is a true reservoir for EBLV-2 or whether these cases represent spill over infections. The genetic structure of M. daubentonii in Western Europe is relatively homogenous and data suggests migratory movements between the UK and the mainland Europe by Daubenton’s bats [59]. Many factors affecting the behaviour of Daubenton’s bats, may allow EBLV-2 to become established and rapidly spread throughout the population. This may account for the greater detection of EBLV-2, albeit at low levels, throughout the UK [60]. However, despite the long distance migratory behaviour of Daubenton’s bats [61] and the extended distribution of EBLV-2 across Europe, there is a striking spatial clustering of EBLV-2 sequences, irrespective of the genomic region analysed.
Within EBLV-2 there is a high level of sequence homogeneity and low evolutionary rate compared to other lyssaviruses. The intra-roost analysis of the three Stokesay EBLV-2 genomes sampled within a 7-year period demonstrates extremely few SNPs and insignificant levels of heterogeneity in the viral reads. However, partial sequence analysis would have been limited in differentiating the isolates, thus confirming the importance of generating full genome/concatenated gene data to improve resolution. The low mean substitution rate of 1.7 × 10−5 of the EBLV-2 concatenated gene dataset is similar to the 4.19 × 10−5 rate recently determined for EBLV-1 [62] but contrasts starkly with rates published for similar concatenated gene or complete genome datasets for RABV, e.g., 3.10 × 10−4 for Turkey canine RABV dataset [63], a comprehensive canine RABV dataset 2.44 × 10−4 [64] and a fox RABV dataset from Greenland 2.5 × 10−4 [65]. The long TMRCAs and low substitution rate suggests limited and localised evolution has occurred over very long periods of time. The distribution and evolutionary analysis suggest an east–west axis possibly in an east (Finland) to west (UK) direction (Figure 1 and Figure 3).
In recent years, multiple novel lyssaviruses have been isolated from bats. Passive bat surveillance, complemented by active surveillance where appropriate, will facilitate the further isolation of known and novel bat lyssaviruses. To fully assess the diversity and evolution of bat lyssaviruses in Europe, all isolates from bat rabies cases should be archived and further characterized. Molecular typing of wildlife rabies cases will facilitate the early detection of CST events which could give rise to the adaption and maintenance of bat lyssaviruses in non-flying mammals. As confirmed in this study, NGS offers unique insights into viral evolution. For a virus with relatively few cases, it is prudent to obtain full genome sequence and study the possible causes for such limited viral evolution and restricted host range [16].

Supplementary Materials

Supplementary materials can be found at www.mdpi.com/1422-0067/19/1/156/s1.

Acknowledgments

This work was supported by grants from Defra, the Scottish Government, and the Welsh Government through grants SE0431 and SV3500, the EU Horizon 2020 research and innovation programme under grant agreement No. 653316, and the German Ministry for Education and Research (grant 01KI1016A). The authors acknowledge support from the Research and Policy for Infectious Disease Dynamics program of the Science and Technology Directorate, Department of Homeland Security, USA. The authors would also like to acknowledge diagnostic staff within the Wildlife Zoonoses and Vector Borne Diseases Research Group at APHA and Richard Ellis (Central Sequencing Unit, APHA). Reproduced with permission of the Controller of Her Majesty’s Stationary Office/Queen’s Printer for Scotland and Animal and Plant Health Agency.

Author Contributions

Lorraine M. McElhinney, Denise A. Marston, Emma L. Wise, Conrad M. Freuling, Antonie Neubauer-Juric, and Tiina Nokireki performed the genome sequencing and analysed the data. All authors (Lorraine M. McElhinney, Denise A. Marston, Emma L. Wise, Conrad M. Freuling, Antonie Neubauer-Juric, Tiina Nokireki, Herve Bourhy, Reto Zanoni, Torfinn Moldal, Engbert A. Kooi, Thomas Müller and Anthony R. Fooks) contributed to the review of submission data and material from previous EBLV-2 cases. All authors contributed to writing the manuscript and gave final approval for publication.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. International Committee on Taxonomy of Viruses (ICTV). Virus Taxonomy: The Classification and Nomenclature of Viruses. The Online (10th) Report of the ICTV. Available online: https://talk.ictvonline.org/ictv-reports/ictv_online_report/ (accessed on 28 November 2017).
  2. Arechiga Ceballos, N.; Vazquez Moron, S.; Berciano, J.M.; Nicolas, O.; Aznar Lopez, C.; Juste, J.; Rodriguez Nevado, C.; Aguilar Setien, A.; Echevarria, J.E. Novel lyssavirus in bat, Spain. Emerg. Infect. Dis. 2013, 19, 793–795. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  3. Gunawardena, P.S.; Marston, D.A.; Ellis, R.J.; Wise, E.L.; Karawita, A.C.; Breed, A.C.; McElhinney, L.M.; Johnson, N.; Banyard, A.C.; Fooks, A.R. Lyssavirus in Indian flying foxes, Sri Lanka. Emerg. Infect. Dis. 2016, 22, 1456–1459. [Google Scholar] [CrossRef] [PubMed]
  4. Streicker, D.G.; Turmelle, A.S.; Vonhof, M.J.; Kuzmin, I.V.; McCracken, G.F.; Rupprecht, C.E. Host phylogeny constrains cross-species emergence and establishment of rabies virus in bats. Science 2010, 329, 676–679. [Google Scholar] [CrossRef] [PubMed]
  5. Banyard, A.C.; Hayman, D.; Johnson, N.; McElhinney, L.; Fooks, A.R. Bats and lyssaviruses. Adv. Virus Res. 2011, 79, 239–289. [Google Scholar] [PubMed]
  6. Botvinkin, A.D.; Poleschuk, E.M.; Kuzmin, I.V.; Borisova, T.I.; Gazaryan, S.V.; Yager, P.; Rupprecht, C.E. Novel lyssaviruses isolated from bats in Russia. Emerg. Infect. Dis. 2003, 9, 1623–1625. [Google Scholar] [CrossRef] [PubMed]
  7. Freuling, C.M.; Beer, M.; Conraths, F.J.; Finke, S.; Hoffmann, B.; Keller, B.; Kliemt, J.; Mettenleiter, T.C.; Muhlbach, E.; Teifke, J.P.; et al. Novel lyssavirus in Natterer’s bat, Germany. Emerg. Infect. Dis. 2011, 17, 1519–1522. [Google Scholar] [CrossRef] [PubMed]
  8. Picard-Meyer, E.; Servat, A.; Robardet, E.; Moinet, M.; Borel, C.; Cliquet, F. Isolation of Bokeloh bat lyssavirus in Myotis nattereri in France. Arch. Virol. 2013, 158, 2333–2340. [Google Scholar] [CrossRef] [PubMed]
  9. Eggerbauer, E.; Troupin, C.; Passior, K.; Pfaff, F.; Hoper, D.; Neubauer-Juric, A.; Haberl, S.; Bouchier, C.; Mettenleiter, T.C.; Bourhy, H.; et al. The recently discovered bokeloh bat lyssavirus: Insights into its genetic heterogeneity and spatial distribution in Europe and the population genetics of its primary host. Adv. Virus Res. 2017, 99, 199–232. [Google Scholar] [PubMed]
  10. Mohr, W. Die Tollwut. Med. Klin. 1957, 52, 1057–1060. [Google Scholar] [PubMed]
  11. Lumio, J.; Hillbom, M.; Roine, R.; Ketonen, L.; Haltia, M.; Valle, M.; Neuvonen, E.; Lahdevirta, J. Human rabies of bat origin in Europe. Lancet 1986, 1. [Google Scholar] [CrossRef]
  12. Bourhy, H.; Kissi, B.; Lafon, M.; Sacramento, D.; Tordo, N. Antigenic and molecular characterization of bat rabies virus in Europe. J. Clin. Microbiol. 1992, 30, 2419–2426. [Google Scholar] [PubMed]
  13. Schatz, J.; Fooks, A.R.; McElhinney, L.; Horton, D.; Echevarria, J.; Vazquez-Moron, S.; Kooi, E.A.; Rasmussen, T.B.; Muller, T.; Freuling, C.M. Bat rabies surveillance in Europe. Zoonoses Public Health 2013, 60, 22–34. [Google Scholar] [CrossRef] [PubMed]
  14. Brookes, S.M.; Aegerter, J.N.; Smith, G.C.; Healy, D.M.; Jolliffe, T.A.; Swift, S.M.; Mackie, I.J.; Pritchard, J.S.; Racey, P.A.; Moore, N.P.; et al. European bat lyssavirus in Scottish bats. Emerg. Infect. Dis. 2005, 11, 572–578. [Google Scholar] [CrossRef] [PubMed]
  15. Amengual, B.; Bourhy, H.; Lopez-Roig, M.; Serra-Cobo, J. Temporal dynamics of European bat lyssavirus type 1 and survival of Myotis myotis bats in natural colonies. PLoS ONE 2007, 2, e566. [Google Scholar] [CrossRef] [PubMed]
  16. Marston, D.A.; Banyard, A.C.; McElhinney, L.M.; Freuling, C.M.; Finke, S.; de Lamballerie, X.; Müller, T.; Fooks, A.R. The lyssavirus host-specificity conundrum—Rabies virus—The exception not the rule. Curr. Opin. Virol. 2018, 28, 68–73. [Google Scholar] [CrossRef] [PubMed]
  17. Badrane, H.; Tordo, N. Host switching in lyssavirus history from the chiroptera to the carnivora orders. J. Virol. 2001, 75, 8096–8104. [Google Scholar] [CrossRef] [PubMed]
  18. Vos, A.; Kaipf, I.; Denzinger, A.; Fooks, A.R.; Johnson, N.; Müller, T. European bat lyssaviruses—An ecological enigma. Acta Chiropterol. 2007, 9, 283–296. [Google Scholar] [CrossRef]
  19. Davis, P.L.; Holmes, E.C.; Larrous, F.; Van der Poel, W.H.; Tjornehoj, K.; Alonso, W.J.; Bourhy, H. Phylogeography, population dynamics, and molecular evolution of European bat lyssaviruses. J. Virol. 2005, 79, 10487–10497. [Google Scholar] [CrossRef] [PubMed]
  20. McElhinney, L.M.; Marston, D.A.; Leech, S.; Freuling, C.M.; van der Poel, W.H.; Echevarria, J.; Vazquez-Moron, S.; Horton, D.L.; Muller, T.; Fooks, A.R. Molecular epidemiology of bat lyssaviruses in Europe. Zoonoses Public Health 2013, 60, 35–45. [Google Scholar] [CrossRef] [PubMed]
  21. Jakava-Viljanen, M.; Nokireki, T.; Sironen, T.; Vapalahti, O.; Sihvonen, L.; HuoVilainen, A. Evolutionary trends of European bat lyssavirus type 2 including genetic characterization of Finnish strains of human and bat origin 24 years apart. Arch. Virol. 2015, 160, 1489–1498. [Google Scholar] [CrossRef] [PubMed]
  22. King, A.; Davies, P.; Lawrie, A. The rabies viruses of bats. Vet. Microbiol. 1990, 23, 165–174. [Google Scholar] [CrossRef]
  23. Kappeler, A. Bat rabies surveillance in Europe. Rabies Bull. Eur. 1989, 13, 12–13. [Google Scholar]
  24. Johnson, N.; Vos, A.; Neubert, L.; Freuling, C.; Mansfield, K.L.; Kaipf, I.; Denzinger, A.; Hicks, D.; Nunez, A.; Franka, R.; et al. Experimental study of European bat lyssavirus type-2 infection in Daubenton’s bats (Myotis daubentonii). J. Gen. Virol. 2008, 89, 2662–2672. [Google Scholar] [CrossRef] [PubMed]
  25. Rabies Bulletin Europe (RBE). Quarterly rabies cases reported in Europe. Rabies Bull. Eur. 1987, 11, 1–42. [Google Scholar]
  26. Megali, A.; Yannic, G.; Zahno, M.L.; Brugger, D.; Bertoni, G.; Christe, P.; Zanoni, R. Surveillance for European bat lyssavirus in Swiss bats. Arch. Virol. 2010, 155, 1655–1662. [Google Scholar] [CrossRef] [PubMed]
  27. Whitby, J.E.; Heaton, P.R.; Black, E.M.; Wooldridge, M.; McElhinney, L.M.; Johnstone, P. First isolation of a rabies-related virus from a Daubenton’s bat in the United Kingdom. Vet. Rec. 2000, 147, 385–388. [Google Scholar] [CrossRef] [PubMed]
  28. Johnson, N.; Selden, D.; Parsons, G.; Healy, D.; Brookes, S.M.; McElhinney, L.M.; Hutson, A.M.; Fooks, A.R. Isolation of a European bat lyssavirus type 2 from a Daubenton’s bat in the United Kingdom. Vet. Rec. 2003, 152, 383–387. [Google Scholar] [CrossRef] [PubMed]
  29. Fooks, A.R.; McElhinney, L.M.; Pounder, D.J.; Finnegan, C.J.; Mansfield, K.; Johnson, N.; Brookes, S.M.; Parsons, G.; White, K.; McIntyre, P.G.; et al. Case report: Isolation of a European bat lyssavirus type 2a from a fatal human case of rabies encephalitis. J. Med. Virol. 2003, 71, 281–289. [Google Scholar] [CrossRef] [PubMed]
  30. Fooks, A.R.; McElhinney, L.M.; Marston, D.A.; Selden, D.; Jolliffe, T.A.; Wakeley, P.R.; Johnson, N.; Brookes, S.M. Identification of a European bat lyssavirus type 2 in a Daubenton’s bat found in Staines, Surrey, UK. Vet. Rec. 2004, 155, 434–435. [Google Scholar] [PubMed]
  31. Fooks, A.R.; Selden, D.; Brookes, S.M.; Johnson, N.; Marston, D.A.; Jolliffe, T.A.; Wakeley, P.R.; McElhinney, L.M. Identification of a European bat lyssavirus type 2 in a Daubenton’s bat found in Lancashire. Vet. Rec. 2004, 155, 606–607. [Google Scholar] [PubMed]
  32. Fooks, A.R.; Marston, D.; Parsons, G.; Earl, D.; Dicker, A.; Brookes, S.M. Isolation of EBLV-2 in a Daubenton’s bat (Myotis daubentonii) found in Oxfordshire. Vet. Rec. 2006, 159, 534–535. [Google Scholar] [CrossRef] [PubMed]
  33. Freuling, C.M.; Kliemt, J.; Schares, S.; Heidecke, D.; Driechciarz, R.; Schatz, J.; Muller, T. Detection of European bat lyssavirus 2 (EBLV-2) in a Daubenton’s bat (Myotis daubentonii) from Magdeburg, Germany. Berl. Munchener Tierarztliche Wochenschr. 2012, 125, 255–258. [Google Scholar]
  34. Schatz, J.; Ohlendorf, B.; Busse, P.; Pelz, G.; Dolch, D.; Teubner, J.; Encarnacao, J.A.; Muhle, R.U.; Fischer, M.; Hoffmann, B.; et al. Twenty years of active bat rabies surveillance in Germany: A detailed analysis and future perspectives. Epidemiol. Infect. 2014, 142, 1155–1166. [Google Scholar] [CrossRef] [PubMed]
  35. Freuling, C.; Grossmann, E.; Conraths, F.J.; Schameitat, A.; Kliemt, J.; Auer, E.; Greiser-Wilke, I.; Muller, T. First isolation of EBLV-2 in Germany. Vet. Microbiol. 2008, 131, 26–34. [Google Scholar] [CrossRef] [PubMed]
  36. Harris, S.L.; Mansfield, K.; Marston, D.A.; Johnson, N.; Pajamo, K.; O’Brien, N.; Black, C.; McElhinney, L.M.; Fooks, A.R. Isolation of european bat lyssavirus type 2 from a Daubenton’s bat (Myotis daubentonii) in Shropshire. Vet. Rec. 2007, 161, 384–386. [Google Scholar] [CrossRef] [PubMed]
  37. Pajermo, K.; Harkess, G.; Goddard, T.; Marston, D.; McElhinney, L.; Johnson, N.; Fooks, A.R. Isolation of european bat lyssavirus type 2 (EBLV-2) in a Daubenton’s bat in the UK with a minimum in Cubation period of 9 months. Rabies Bull. Eur. 2008, 32, 6–8. [Google Scholar]
  38. Banyard, A.C.; Johnson, N.; Voller, K.; Hicks, D.; Nunez, A.; Hartley, M.; Fooks, A.R. Repeated detection of european bat lyssavirus type 2 in dead bats found at a single roost site in the UK. Arch. Virol. 2009, 154, 1847–1850. [Google Scholar] [CrossRef] [PubMed]
  39. Horton, D.L.; Voller, K.; Haxton, B.; Johnson, N.; Leech, S.; Goddard, T.; Wilson, C.; McElhinney, L.M.; Fooks, A.R. European bat lyssavirus type 2 in a Daubenton’s bat in Scotland. Vet. Rec. 2009, 165, 383–384. [Google Scholar] [CrossRef] [PubMed]
  40. Jakava-Viljanen, M.; Lilley, T.; Kyheroinen, E.M.; HuoVilainen, A. First encounter of European bat lyssavirus type 2 (EBLV-2) in a bat in Finland. Epidemiol. Infect. 2010, 138, 1581–1585. [Google Scholar] [CrossRef] [PubMed]
  41. Neubauer-Juric, A.; Krebs, S.; Schille, A.-U.; Schürmann, E.-M.; Blum, H. Tollwut in bayern. In Jahrestagung der DVG-FG AVID; Kaiser Medien GmbH: Nürnberg, Germany, 2014. [Google Scholar]
  42. Moldal, T.; Vikøren, T.; Cliquet, F.; Marston, D.A.; van der Kooij, J.; Madslien, K.; Ørpetveit, I. First detection of European bat lyssavirus type 2 (EBLV-2) in Norway. BMC Vet. Res. 2017, 13, 216. [Google Scholar] [CrossRef] [PubMed]
  43. Johnson, N.; Goddard, T.M.; Goharriz, H.; Wise, E.; Jennings, D.; Selden, D.; Marston, D.A.; Banyard, A.C.; McElhinney, L.M.; Fooks, A.R. Two eblv-2 infected Daubenton’s bats detected in the north of England. Vet. Rec. 2016, 179, 311–312. [Google Scholar] [CrossRef] [PubMed]
  44. Nokireki, T.; Sironen, T.; Smura, T.; Karkamo, V.; Sihvonen, L.; Gadd, T. Second case of European bat lyssavirus type 2 detected in a Daubenton’s bat in Finland. Acta Vet. Scand. 2017, 59. [Google Scholar] [CrossRef] [PubMed]
  45. Hagner, N.; Ek-Kommonen, C.; Lokki, J.; Neuvonen, E.; Stjernberg, T.; Valle, M.; Veijalainen, P. No bat rabies found in Finland 1986. In European Bat Research; Hanák, V., Horácek, I., Gaisler, J., Eds.; Charles University Press: Praha, Czech Republic, 1989; pp. 631–635. [Google Scholar]
  46. Nokireki, T.; HuoVilainen, A.; Lilley, T.; Kyheroinen, E.M.; Ek-Kommonen, C.; Sihvonen, L.; Jakava-Viljanen, M. Bat rabies surveillance in Finland. BMC Vet. Res. 2013, 9, 174. [Google Scholar] [CrossRef] [PubMed]
  47. Van der Poel, W.H.; Van der Heide, R.; Verstraten, E.R.; Takumi, K.; Lina, P.H.; Kramps, J.A. European bat lyssaviruses, The Netherlands. Emerg. Infect. Dis. 2005, 11, 1854–1859. [Google Scholar] [CrossRef] [PubMed]
  48. Amengual, B.; Whitby, J.E.; King, A.; Cobo, J.S.; Bourhy, H. Evolution of European bat lyssaviruses. J. Gen. Virol. 1997, 78, 2319–2328. [Google Scholar] [CrossRef] [PubMed]
  49. Marston, D.A.; McElhinney, L.M.; Johnson, N.; Muller, T.; Conzelmann, K.K.; Tordo, N.; Fooks, A.R. Comparative analysis of the full genome sequence of European bat lyssavirus type 1 and type 2 with other lyssaviruses and evidence for a conserved transcription termination and polyadenylation motif in the G-L 3′ non-translated region. J. Gen. Virol. 2007, 88, 1302–1314. [Google Scholar] [CrossRef] [PubMed]
  50. Lo, K.W.; Naisbitt, S.; Fan, J.S.; Sheng, M.; Zhang, M. The 8-kda dynein light chain binds to its targets via a conserved (K/R)XTQT motif. J. Biol. Chem. 2001, 276, 14059–14066. [Google Scholar] [CrossRef] [PubMed]
  51. Marston, D.A.; McElhinney, L.M.; Ellis, R.J.; Horton, D.L.; Wise, E.L.; Leech, S.L.; David, D.; de Lamballerie, X.; Fooks, A.R. Next generation sequencing of viral RNA genomes. BMC Genom. 2013, 14. [Google Scholar] [CrossRef] [PubMed]
  52. Brunker, K.; Marston, D.A.; Horton, D.L.; Cleaveland, S.; Fooks, A.R.; Kazwala, R.; Ngeleja, C.; Lembo, T.; Sambo, M.; Mtema, Z.J.; et al. Elucidating the phylodynamics of endemic rabies virus in Eastern Africa using whole-genome sequencing. Virus Evol. 2015, 1. [Google Scholar] [CrossRef] [PubMed]
  53. Milne, I.; Stephen, G.; Bayer, M.; Cock, P.J.; Pritchard, L.; Cardle, L.; Shaw, P.D.; Marshall, D. Using tablet for visual exploration of second-generation sequencing data. Brief. Bioinform. 2013, 14, 193–202. [Google Scholar] [CrossRef] [PubMed]
  54. Li, H.; Handsaker, B.; Wysoker, A.; Fennell, T.; Ruan, J.; Homer, N.; Marth, G.; Abecasis, G.; Durbin, R.; Genome Project Data Processing, S. The sequence alignment/map format and samtools. Bioinformatics 2009, 25, 2078–2079. [Google Scholar] [CrossRef] [PubMed]
  55. Lemey, P.; Rambaut, A.; Drummond, A.J.; Suchard, M.A. Bayesian phylogeography finds its roots. PLoS Comput. Biol. 2009, 5, e1000520. [Google Scholar] [CrossRef] [PubMed]
  56. Tamura, K.; Stecher, G.; Peterson, D.; Filipski, A.; Kumar, S. Mega6: Molecular evolutionary genetics analysis version 6.0. Mol. Biol. Evol. 2013, 30, 2725–2729. [Google Scholar] [CrossRef] [PubMed]
  57. Baele, G.; Lemey, P.; Bedford, T.; Rambaut, A.; Suchard, M.A.; Alekseyenko, A.V. Improving the accuracy of demographic and molecular clock model comParison while accommodating phylogenetic uncertainty. Mol. Biol. Evol. 2012, 29, 2157–2167. [Google Scholar] [CrossRef] [PubMed]
  58. Stebbings, R.E.; Griffith, F. Distribution and Status of Bats in Europe; Institute of Terrestrial Ecology: Huntingdon, UK, 1986; p. 142. [Google Scholar]
  59. Smith, G.C.; Aegerter, J.N.; Allnutt, T.R.; MacNicoll, A.D.; Learmount, J.; Hutson, A.M.; Atterby, H. Bat population genetics and lyssavirus presence in Great Britain. Epidemiol. Infect. 2011, 139, 1463–1469. [Google Scholar] [CrossRef] [PubMed]
  60. Harris, S.L.; Aegerter, J.N.; Brookes, S.M.; McElhinney, L.M.; Jones, G.; Smith, G.C.; Fooks, A.R. Targeted surveillance for European bat lyssaviruses in English bats (2003–06). J. Wildl. Dis. 2009, 45, 1030–1041. [Google Scholar] [CrossRef] [PubMed]
  61. Hutterer, R.; Ivanova, T.; Meyer-Cords, C.; Rodrigues, L. Bat migrations in Europe: A review of banding data and literature. Ger. Agency Nat. Conserv. 2005, 28, 83–84. [Google Scholar]
  62. Troupin, C.; Picard-Meyer, E.; Dellicour, S.; Casademont, I.; Kergoat, L.; Lepelletier, A.; Dacheux, L.; Baele, G.; Monchatre-Leroy, E.; Cliquet, F.; et al. Host genetic variation does not determine spatio-temporal patterns of European bat 1 lyssavirus. Genome Boil. Evol. 2017, 9, 3202–3213. [Google Scholar] [CrossRef] [PubMed]
  63. Marston, D.A.; Horton, D.L.; Nunez, J.; Ellis, R.J.; Orton, R.J.; Johnson, N.; Banyard, A.; McElhinney, L.; Freuling, C.; Firat, M.; et al. Genetic analysis of a rabies virus host shift event reveals within-host viral dynamics in a new host. Virus Evol. 2017. [Google Scholar] [CrossRef] [PubMed]
  64. Troupin, C.; Dacheux, L.; Tanguy, M.; Sabeta, C.; Blanc, H.; Bouchier, C.; Vignuzzi, M.; Duchene, S.; Holmes, E.C.; Bourhy, H. Large-scale phylogenomic analysis reveals the complex evolutionary history of rabies virus in multiple carnivore hosts. PLoS Pathog. 2016, 12, e1006041. [Google Scholar] [CrossRef] [PubMed]
  65. Hanke, D.; Freuling, C.M.; Fischer, S.; Hueffer, K.; Hundertmark, K.; Nadin-Davis, S.; Marston, D.; Fooks, A.R.; Botner, A.; Mettenleiter, T.C.; et al. Spatio-temporal analysis of the genetic diversity of arctic rabies viruses and their reservoir hosts in Greenland. PLoS Negl. Trop. Dis. 2016, 10, e0004779. [Google Scholar] [CrossRef] [PubMed]
Ijms 19 00156 g001 550
Figure 1. Distribution of 34 confirmed cases of EBLV-2 (Table 1) in Europe between 1985 and 2017 (WHO Rabies Bulletin Europe).
Figure 1. Distribution of 34 confirmed cases of EBLV-2 (Table 1) in Europe between 1985 and 2017 (WHO Rabies Bulletin Europe).
Ijms 19 00156 g001
Ijms 19 00156 g002a 550Ijms 19 00156 g002b 550
Figure 2. Molecular phylogenetic analysis of (a) nucleoprotein sequences (1353 nucleotides) and (b) glycoprotein sequences (1569 nucleotides) from 29 EBLV-2 cases and 36 non EBLV-2 lyssavirus sequences (n = 65). Sequences were analysed by the Maximum Likelihood method. Bootstrap values > 70% are shown for each branch and analysis was performed using the GTR + G + I model in MEGA6 with 1000 bootstrap samples. The EBLV-2 sequences are coloured depending upon sampling location (see 3 letter codes within case names) (a) N1353; (b) G1569.
Figure 2. Molecular phylogenetic analysis of (a) nucleoprotein sequences (1353 nucleotides) and (b) glycoprotein sequences (1569 nucleotides) from 29 EBLV-2 cases and 36 non EBLV-2 lyssavirus sequences (n = 65). Sequences were analysed by the Maximum Likelihood method. Bootstrap values > 70% are shown for each branch and analysis was performed using the GTR + G + I model in MEGA6 with 1000 bootstrap samples. The EBLV-2 sequences are coloured depending upon sampling location (see 3 letter codes within case names) (a) N1353; (b) G1569.
Ijms 19 00156 g002aIjms 19 00156 g002b
Ijms 19 00156 g003 550
Figure 3. Bayesian maximum clade credibility (MCC) tree of 29 EBLV-2 concatenated gene sequences (10818 nucleotide sites comprising the nucleoprotein, matrix, phosphoprotein, glycoprotein, and polymerase genes) with branch lengths scaled in time by enforcing a relaxed molecular clock. The EBLV-2 strains are coloured to indicate the location of sampling. Mean substitution rate of the dataset was 1.7 × 10−5.
Figure 3. Bayesian maximum clade credibility (MCC) tree of 29 EBLV-2 concatenated gene sequences (10818 nucleotide sites comprising the nucleoprotein, matrix, phosphoprotein, glycoprotein, and polymerase genes) with branch lengths scaled in time by enforcing a relaxed molecular clock. The EBLV-2 strains are coloured to indicate the location of sampling. Mean substitution rate of the dataset was 1.7 × 10−5.
Ijms 19 00156 g003
Ijms 19 00156 g004 550
Figure 4. A circular map of the EBLV2 genome. Numbers indicate genomic positions at which heterogeneity is seen between the three genomes from Stokesay Castle. Black positions are located within coding regions, teal positions are in non-coding regions. Position 2015 (asterisk) has a base change corresponding to an amino acid substitution. Orange blocks indicate known domains within each protein where heterogeneity between the three genomes is observed. Underlined positions indicate those at which >1% heterogeneity is seen within the host. Domains key: A = T-cell epitopes; B = L binding region; C = soluble N binding region; D = LC8 protein binding region; E = Domain III; F = Domain IV; G = Domain V.
Figure 4. A circular map of the EBLV2 genome. Numbers indicate genomic positions at which heterogeneity is seen between the three genomes from Stokesay Castle. Black positions are located within coding regions, teal positions are in non-coding regions. Position 2015 (asterisk) has a base change corresponding to an amino acid substitution. Orange blocks indicate known domains within each protein where heterogeneity between the three genomes is observed. Underlined positions indicate those at which >1% heterogeneity is seen within the host. Domains key: A = T-cell epitopes; B = L binding region; C = soluble N binding region; D = LC8 protein binding region; E = Domain III; F = Domain IV; G = Domain V.
Ijms 19 00156 g004
Table
Table 1. Confirmed cases of EBLV-2 in Europe 1985–2017 and sequences used to generate phylogenetic trees.
Table 1. Confirmed cases of EBLV-2 in Europe 1985–2017 and sequences used to generate phylogenetic trees.
YearOriginal ReferenceAPHA Virus Ref.PI Virus Ref.LocationCountrySpeciesN400 Data Accession No.Genome Accession No.Case or Sequence Ref.
1985FINMANRV89007FINHelsinkiFinlandHumanAY062091JX129233, KY688151[11]
198741729RV309018HOLWommelsNetherlandsM.dasycEU293114EU293114[19]
198747072RV299482HOLTjerkwerdNetherlandsM.dasycU89480KY688145[25]
198992666RV22894112HOLAndijkNetherlandsM.dasycAY062089KY688146[19]
1992TW1814/92RV594, RV2478-PlaffeienSwitzerlandM.daub.AY212117KY688133[26]
1993--9375HOLRodenNetherlandsM.dasycAY863404KY688152[19]
1993TW1392/93RV621, RV24799337SWIVersoixSwitzerlandM.daubAY212118KY688140[19]
199618/96RV628EBL2GBSussexEnglandM.daubU89478KY688136[27]
2002105/02RV1332-LancashireEnglandM.daubAY212120KY688150[28]
2002Human caseRV1333-AngusScotlandHumanEF157977EF157977[29]
2002TW0118/02RV248002053SWIGenevaSwitzerlandM.daubAY863408KY688132[19]
2004603/04RV1787-SurreyEnglandM.daubJQ796807KY688142[30]
2003696/04RV1788-LancashireEnglandM.daubJQ796808KF155004[31]
2006672/06RV2159-OxfordshireEnglandM.daubJQ796809KY688144[32]
200618856RV2506-MagdeburgGermanyM.daubJQ796805KY688135[33]
200625538--ThuringiaGermanyM.daubKF826115KY688137[34]
200716618RV2505-Bad BuchauGermanyM.daubGU227648KY688138[35]
2007762/07RV2336-ShropshireEnglandM.daubJQ796810KY688139[36]
2008166/08RV2418-SurreyEnglandM.daubJQ796811KY688141[37]
20081218/08RV2473-ShropshireEnglandM.daubJQ796812KY688147[38]
2009882/09RV2482-West LothianScotlandM.daubJQ796806KY688143[39]
2009ra3278/09--Nr TurkuFinlandM.daubGU002399JX129232[40]
2012F500/LGL--Ingolstadt, BavariaGermanyM.daub-KY688149[34,41]
201331452--Gießen, HesseGermanyM.daubKF826149KY688134[34]
2013F860/LGL--Ingolstadt, BavariaGermanyM.daub-Not availableThis study
2014281/14RV2974-ShropshireEnglandM.daub-KY688148This study
2015421/15RV3158-Newtown, PowysWalesM.daub-KY688153This study
20152015-04-18832--OpplandNorwayM.daub-KY688154[42]
2016307/16RV3369-North YorkshireEnglandM.daub-KY688156This study
2016432/16RV3370-NorthumberlandEnglandM.daub-KY688155[43]
20162Finland2016--Inkoo, UusimaaFinlandM.daubMF326269Not available[44]
2017591/17RV3385-DerbyshireEnglandM.daub-Not availableThis study
Active Surveillance Detection of Viral RNA in Oro-Pharyngeal Swab Only:
200880063M08/09-PerthshireScotlandM.daubJQ796804Not available[20]
200970--GenthodSwitzerlandM.daubHM067110Not available[26]
Key—not known. M.daub—Myotis daubentonii; M.dasyc—Myotis dasycneme.

Share and Cite

MDPI and ACS Style

McElhinney, L.M.; Marston, D.A.; Wise, E.L.; Freuling, C.M.; Bourhy, H.; Zanoni, R.; Moldal, T.; Kooi, E.A.; Neubauer-Juric, A.; Nokireki, T.; et al. Molecular Epidemiology and Evolution of European Bat Lyssavirus 2. Int. J. Mol. Sci. 2018, 19, 156. https://doi.org/10.3390/ijms19010156

AMA Style

McElhinney LM, Marston DA, Wise EL, Freuling CM, Bourhy H, Zanoni R, Moldal T, Kooi EA, Neubauer-Juric A, Nokireki T, et al. Molecular Epidemiology and Evolution of European Bat Lyssavirus 2. International Journal of Molecular Sciences. 2018; 19(1):156. https://doi.org/10.3390/ijms19010156

Chicago/Turabian Style

McElhinney, Lorraine M., Denise A. Marston, Emma L. Wise, Conrad M. Freuling, Hervé Bourhy, Reto Zanoni, Torfinn Moldal, Engbert A. Kooi, Antonie Neubauer-Juric, Tiina Nokireki, and et al. 2018. "Molecular Epidemiology and Evolution of European Bat Lyssavirus 2" International Journal of Molecular Sciences 19, no. 1: 156. https://doi.org/10.3390/ijms19010156

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop