Next Article in Journal
Real Time Analysis of Bovine Viral Diarrhea Virus (BVDV) Infection and Its Dependence on Bovine CD46
Next Article in Special Issue
Using the LN34 Pan-Lyssavirus Real-Time RT-PCR Assay for Rabies Diagnosis and Rapid Genetic Typing from Formalin-Fixed Human Brain Tissue
Previous Article in Journal
Systematic Identification of Host Immune Key Factors Influencing Viral Infection in PBL of ALV-J Infected SPF Chicken
Previous Article in Special Issue
Rabies Surveillance Identifies Potential Risk Corridors and Enables Management Evaluation
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Viruses
Volume 12
Issue 1
10.3390/v12010115
viruses-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:
Article

Population- and Variant-Based Genome Analyses of Viruses from Vaccine-Derived Rabies Cases Demonstrate Product Specific Clusters and Unique Patterns

by
Sten Calvelage
1,
Marcin Smreczak
2,
Anna Orłowska
2,
Conrad Martin Freuling
3,*,
Thomas Müller
3,
Christine Fehlner-Gardiner
4,
Susan Nadin-Davis
4,
Dirk Höper
1 and
Paweł Trębas
2
1
Institute of Diagnostic Virology, Friedrich-Loeffler-Institut, 17493 Greifswald-Insel Riems, Germany
2
Department of Virology, National Veterinary Research Institute, 24-100 Puławy, Poland
3
Institute of Molecular Virology and Cell Biology, Friedrich-Loeffler-Institut, WHO Collaborating Centre for Rabies Surveillance and Research, OIE Reference Laboratory for Rabies, 17493 Greifswald-Insel Riems, Germany
4
National Reference Centre for Rabies, Ottawa Laboratory–Fallowfield, Canadian Food Inspection Agency, Ottawa, ON K2H 8P91, Canada
*
Author to whom correspondence should be addressed.
Viruses 2020, 12(1), 115; https://doi.org/10.3390/v12010115
Submission received: 17 December 2019 / Revised: 10 January 2020 / Accepted: 13 January 2020 / Published: 17 January 2020
(This article belongs to the Special Issue Rabies Virus: Knowledge Gaps and Challenges to Elimination)
Browse Figures
Versions Notes

Abstract

:
Rabies in wildlife has been successfully controlled in parts of Europe and North America using oral rabies vaccination, i.e., the distribution of baits containing live-attenuated virus strains. Occasionally, these vaccines caused vaccine virus-induced rabies cases. To elucidate the mechanisms of genetic selection and the effect of viral populations on these rabies cases, a next generation sequencing approach as well as comprehensive data analyses of the genetic diversity of Street Alabama Dufferin (SAD) and ERA vaccine virus strains and vaccine-induced rabies cases from Canada and several European countries were conducted. As a result, twelve newly generated sets of sequencing data from Canada and Poland were added to a pool of previously investigated samples. While the population-based analysis showed a segregation of viruses of ERA vaccine-induced rabies cases from those of SAD Bern original (SAD Bernorig)-derived rabies cases, the in-depth variant analysis revealed three distinct combinations of selected variants for the ERA vaccine-induced cases, suggesting the presence of multiple replication-competent haplotypes in the investigated ERA-BHK21 vaccine. Our findings demonstrate the potential of a deep sequencing approach in combination with comprehensive analyses on the consensus, population, and variant level.

1. Introduction

Vaccination programs are one of the most effective means of controlling infectious diseases [1] and with the development of oral vaccines and bait delivery systems, the elimination of diseases circulating in wildlife populations has become a realistic possibility. The large-scale oral rabies vaccination (ORV) campaigns that have eliminated fox-mediated rabies from Western Europe and North America [2,3] and substantially reduced disease incidence in central Europe [2] are pre-eminent examples for the success of such control programs.
ORV programs in foxes are aimed at increasing herd immunity in the target population using oral rabies vaccines distributed into the environment. Over the past four decades several oral rabies vaccines, mainly live replication-competent attenuated rabies virus vaccines, have been successfully used in ORV campaigns [4]. In Canada for example, the ERA-BHK21vaccine virus strain, a derivative of the cell culture adapted vaccine virus strain Street Alabama Dufferin (SAD) [5,6], was the only live attenuated vaccine deployed in fox ORV campaigns [7]. In Europe, with the exception of a recombinant vaccinia virus expressing the rabies virus glycoprotein, all constructs have been based on live attenuated rabies virus strains [4], derived from the SAD Bern original (SAD Bernorig) vaccine virus strain, a successor of the ERA strain [2,8]. While all these vaccines have been highly efficient in fox rabies control, the first generation of SAD-derived vaccines demonstrated residual pathogenicity in non-target species particularly in rodents [9,10,11,12]. Although several cases of vaccine virus-induced rabies were observed even in species other than rodents over the course of vaccination campaigns in a number of countries, including Germany [1], Austria [1], Slovenia, Romania [13], Poland, and Canada [7], such cases were without epidemiological relevance [1,7,8].
While previous analyses using high-throughput sequencing approaches revealed substantial genetic heterogeneity within commercial SAD-derived oral rabies virus vaccines [14], the sub-consensus genetic heterogeneity of viruses isolated from these vaccine virus-induced rabies cases on the contrary revealed nearly clonal genotypes, indicating the presence of a strong bottleneck during infection [15].
In this study, we attempted to further analyze and elucidate the mechanisms of genetic selection using a combined deep sequencing, variant, and haplotype analysis of a large set of vaccine-induced rabies cases that included additional SAD- and ERA-induced rabies cases from Poland and Canada alongside vaccine virus batches. With this comprehensive approach, we could demonstrate the utility of this approach for identity and genetic stability (revision to virulence) testing of vaccines with a heterogeneous genetic background. Additionally, we were interested to know whether viruses from vaccine-induced rabies cases differ from virulent field rabies viruses (RABV) at a population level.

2. Materials and Methods

2.1. Samples

In this study, original brain material of additional ERA-BHK21 and SAD Bernorig-derived vaccine-induced rabies cases from Canada [7] and Poland, respectively, were analyzed (Table 1). In addition, new batches of the attenuated rabies virus vaccine SAD Bern (Lysvulpen) were obtained and included in the analysis (Table 2). In order to check for validity of the PCR based and Illumina generated sequence reads for the use of the methods to be applied, we included datasets of two case samples (C/CAN/1994/MB and C/CAN/1996/B) as well as a dataset of the ERA/lot16 vaccine batch already sequenced on an IonTorrent PGM platform from a previous study [15].

2.2. Nucleic Acid Extraction, Sample Processing, and Sequencing

Samples were processed according to two different methodical approaches. Methodology 1 corresponds to the published protocol of Wylezich and colleagues [16]. Briefly, for sample disintegration prior to nucleic acid extraction, 20 mg tissue material or 200 µL of the vaccine virus preparation from the blister were cryofractured using the cryoPREP CP02 (Covaris, Brighton, United Kingdom). The pulverized material was lysed with pre-heated (56 °C) buffer AL (Qiagen, Hilden, Germany). A threefold volume of Trizol was added and RNA extracted using the RNeasy Mini Kit (Qiagen) followed by on-column DNase I digestion (Qiagen). Subsequently, the obtained RNA was converted into double stranded cDNA using a combination of SuperScript™ IV First-Strand cDNA Synthesis System (Invitrogen/Thermo Fisher Scientific, Waltham, MA, USA) and NEBNext® Ultra™ II Non-Directional RNA Second Strand Synthesis Module (New England Biolabs, Ipswich, MA, USA). The generated cDNA was fragmented using the Covaris M220 (Covaris) and subsequently converted to Ion Torrent compatible libraries using the GeneRead L Core Kit (Qiagen) and IonXpress Barcode Adaptors (Thermo Fisher Scientific, Waltham, MA, USA) followed by size selection as described [16]. After quality control (Agilent 2100 Bioanalyzer; High Sensitivity DNA Kit, Agilent Technologies, Santa Clara, CA, USA) and quantification (KAPA Library Quantification Kit—Ion Torrent PGM Uni; KAPA Biosystems/Roche, Basel, Switzerland), the libraries including the Ion S5 Calibration Standard were sequenced on an Ion 530 chip with an Ion S5 XL System (Thermo Fisher Scientific) in 400 bp-mode according to the manufacturer’s instructions.
Methodology 2 followed sample preparation as described before [17]. In short, RNA was prepared using a protocol combining Trizol and a semi-automated MagMax™ (Applied Biosystems/Thermo Fisher Scientific, Waltham, MA, USA) extraction system. Subsequently, dsDNA for input into library preparation was generated by reverse-transcription polymerase chain reaction (RT-PCR), yielding overlapping amplicons of the RABV genome (Tables S1 and S2, Supplementary File). The generated amplicons were individually quantified for each sample utilizing a Qubit fluorometer (ThermoFisher Scientific) and pooled equimolarly followed by library generation using the Nextera XT kit (Illumina, San Diego, CA, USA). Sequencing was performed with an Illumina MiSeq instrument generating 250 base paired end reads.

2.3. Population-Based Analysis, Sequence Assembly, and In-Depth Variant Analysis

Data analysis was performed essentially as described before with the respective published data included [14,15] (Tables S3–S6, Supplementary File). After trimming for quality of reads (all datasets), the dataset sizes were adjusted to match the size of the published datasets [14,15]. All sequence reads were mapped along a strict majority consensus sequence and its reverse complement (454 Sequencing System Software v3.0; Roche) that was previously generated and introduced via a MAFFT (multiple sequence alignment using fast Fourier transform) alignment of published SAD-related sequences [15]. Using R (version 3.6.0) and RStudio (version 1.2.1335) in combination with additional packages (Table 3), base frequencies of the forward and reverse mapping were calculated for each position of the reference sequence. Subsequently, the data were corrected for strand-bias and pairwise Manhattan distances were calculated for each population. After non-linear multidimensional scaling, the distances were displayed in 2-dimensional plots [15]. For reasons of simplification, previously sequenced in vitro selected attenuated rabies vaccine strains, i.e., P5/88, VA1, SAG2, SAD B19CS, and SAD B19P1 (Table S3, Supplementary File), were excluded to achieve a more convenient presentation of the SAD Bernorig-derived vaccine-induced rabies cases analyzed in this study.
Full genome sequences of all samples listed in Table 1 and Table 2 were obtained by de-novo assembly of full or partial matching reads from the forerun mapping (454 Sequencing System Software v3.0; Roche) and submitted to the European Nucleotide Archive (ENA) under the study accession PRJEB35810.
For in-depth variant analysis, strand-bias corrected data of base frequencies were used as the basis for comparison of single nucleotide variants between vaccine virus strains and their induced rabies cases. Here, all nucleotide variants with a variant frequency of at least five percent (Table S7, Supplementary File) were considered. For sample datasets generated with both IonTorrent and Illumina (C/CAN/1994/MB—C/CAN/1994N6762M, C/CAN/1996/B—C/CAN/1996N5367, and ERA/lot16), variants were taken into account if they were present in both datasets and their mean variant frequency was equal or above five percent. In parallel, fully or partially mapped sequence reads were analyzed by Geneious Prime (2019.2.3; build 2019-09-24) for an additional variation analysis using standard settings and a minimum variant frequency threshold of 0.05 to confirm positions and frequencies of single nucleotide variants found in the strand-bias corrected dataset. Noteworthy, the populations determined for samples C/CAN/1994/MB—C/CAN/1994N6762M and C/CAN/1996/B—C/CAN/1996N5367 from Ion Torrent and Illumina data were nearly identical, putting the related viruses in close proximity in the distance analysis.

3. Results

To expand the knowledge and understanding of the dynamics of ERA-BHK21 and SAD Bernorig-derived oral rabies vaccines, twelve additional datasets of vaccine-induced rabies cases and their corresponding vaccines were investigated (Table 1 and Table 2).

3.1. Viruses from Vaccines Form Product Specific Clusters in Population Analysis

Overall, the distance plot generated from the population data segregates the different oral rabies vaccine viruses into two major clusters, i.e., ERA-BHK21- and SAD Bernorig-derived vaccines (Figure 1), confirming previous observations [14,15]. When additional SAD Bern (Lysvulpen) batches were included, the SAD Bernorig-derived cluster appeared to form two sub-clusters. While the majority of SAD B19 and SAD Bern vaccine batches clearly separated into these sub-clusters, one batch of each group was dislocated. Namely, batch B/SAD/B19/958 clustered with SAD Bern (Lysvulpen) vaccine batches, whereas vaccine batch B/SAD/Bern/0213 was shifted to the vaccine-induced rabies cases (Figure 1) as previously described [15].
Population analysis of the viruses isolated from vaccine-induced rabies cases revealed two distinct patterns. Those viruses that were associated with SAD Bernorig-derived oral rabies vaccines (including the newly analyzed Polish case C/POL/2017/B) did not cluster with the related vaccines as shown before [15]. In contrast, the viruses from the ERA vaccine-induced rabies cases all clustered closely together with the ERA-BHK21 vaccine batch (B/ERA/lot16, Figure 1). To further elucidate the observed differences, we conducted an in-depth analysis of the relations between viruses from SAD Bernorig- and ERA vaccine-induced rabies cases to their progenitor vaccines.

3.2. ERA-Derived Subclusters Comprise Unique Patterns of Base Exchanges

The ERA-BHK21 vaccine virus population displayed a high genetic heterogeneity represented by diverse single nucleotide variants distributed across the N, P, G, and L genes (Figure 2), supporting previous findings [15]. Interestingly, the in-depth variant analysis revealed distinct dependencies of the ERA vaccine-induced cases from vaccine variants, with three specific combinations of single nucleotide variants (SNVs) being selected from the progenitor population at the consensus level (Figure 2). Most of these differences were derived from single nucleotide variants (SNVs) already present in the progenitor ERA-BHK21 vaccine virus strain, while a few additional mutations at the consensus level were of spontaneous origin. While all but one of the latter mutations appeared to be unique, the selected SNVs could be found across several viruses from ERA vaccine-induced rabies cases. Thus, the analyses of viruses from the ERA vaccine-induced cases confirmed distinct combinations of SNVs. According to these combinations, three groups of viruses from ERA vaccine-induced rabies cases could be defined (Table 4, Figure 2).
Each of these groups was characterized by a unique set of differences at the consensus level ranging between three and nine genome positions that were not present in any of the other groups (Table 4, Figure 2). This suggests the presence of at least three replication-competent virus haplotypes in the progenitor ERA-BHK21 vaccine virus strain that were selected in vivo by a genetic bottleneck in affected animals [15]. However, although the limited number of ERA vaccine-induced rabies cases hampers a statistical analysis, there seems to be no spatio-temporal nor species-specific correlation among the three groups of viruses from ERA vaccine-induced rabies cases.
With 8–14 differences at the consensus level, sequence diversity of the ERA vaccine-induced rabies cases and their progenitor vaccine virus strain was relatively high compared to the viruses of the SAD Bernorig rabies cases. This observation seemingly contradicts results from the population-based analysis where the ERA vaccine-induced cases clustered closely together with the ERA-BHK21 vaccine batch. However, differences at the consensus level of the viruses of ERA vaccine-induced rabies cases were already included in the progenitor vaccine population (ERA-BHK21) and likewise, in some variants, the former consensus sequence was present as a minor variant. For instance, in each of two ERA vaccine-induced rabies cases (C/CAN/1992N8944 and C/CAN/1994N6762M) one additional single nucleotide variant (A10840[G,A] and C3085[T,C]) appeared in the viral populations; in both cases, the progenitor consensus was maintained as a minor variant (Table 5).

3.3. SAD Bernorig-Derived Cases are Defined by High Sequence Identities Relative to Their Related Vaccines

Strikingly, in contrast to the ERA vaccine-induced rabies cases, for the SAD Bernorig-induced rabies cases no SNVs already present in the progenitor vaccines were selected, even though a higher number of variants were detected in the SAD Bernorig-derived vaccine virus strain (Nmean ≈ 47, see Table S8, Supplementary File) than in the ERA-BHK21 strain (N = 31, see Figure 2 and Table S7). Also, the consensus sequences of viruses from the SAD Bernorig vaccine-induced rabies cases showed a high sequence identity to their progenitor vaccine virus strain at the consensus level, ranging from zero to four nucleotide exchanges, confirming previous studies [13,15,18].

3.4. Known Antigenic Sites Are Rarely Affected by Amino Acid Exchanges in Vaccine-Induced Cases

To elucidate whether the induction of disease in these vaccine-induced rabies cases was a result of mutations at known antigenic sites [19,20,21,22,23,24,25,26,27,28,29,30,31,32], all sequences were screened for nucleotide exchanges at the respective genome positions. The analysis revealed no nucleotide exchanges between the SAD Bernorig-derived vaccine virus strains and their induced rabies cases in any of the known antigenic sites, including SNVs on these particular positions. Furthermore, none of the antigenic sites were affected by mutations in the ERA vaccine-induced rabies cases, except for one amino acid (AA) exchange Thr178Ile at an interferon-antagonist motif in the P-gene [28] found in all members of group 2 (Table 4). Whether this single AA exchange has any influence on pathogenicity remains speculative, as the entire motif contains eleven AAs and only the deletion of segments or the complete motif was studied [28]. On the other hand, if pathogenicity is increased by this presumed haplotype of group 2, more vaccine-induced rabies cases of this variant would be expected.

3.5. Viruses from Vaccine-Induced Cases Can Clearly Be Distinguished from Field RABV Viruses

When selected field RABV viruses were included in the Manhattan distance matrix analysis, the attenuated vaccine virus strains and viruses from their vaccine-induced rabies cases were clearly separated from them by high distances (Figure 3). Furthermore, dimensions of the plot increased substantially compared to Figure 1, resulting in the condensation of minor differences and the formation of a single tight cluster.

4. Discussion

Although various studies have focused on oral rabies vaccine-induced cases [1,7,13,18,33], only a few studies have comprehensively investigated those cases and their progenitor vaccine viruses at a population level [14,15]. The latter showed a high genetic heterogeneity in different vaccine virus strains as opposed to viruses from vaccine-induced rabies cases, indicating a substantial loss in viral population diversity. In this study, we could identify two different vaccine virus-dependent patterns in the genetic markup of viruses from vaccine-induced rabies cases. This was only possible by including further samples of ERA [7] and SAD Bern vaccine-induced rabies cases in the in-depth variant analysis (Table 1). Hence, in total 20 datasets from viruses of 15 cases of vaccine-induced rabies from Europe and Canada as well as 15 related batches of the ERA-BHK21, SAD Bern, and SAD B19 vaccine strains that were collected between 1991 and 2017 [7,14,15] were considered in this analysis.
While viruses from ERA vaccine-derived cases seem to descend from different replication-competent virus haplotypes that were already present within the progenitor ERA-BHK21 vaccine virus strain (Figure 2), the absence of SNV combinations in viruses from SAD Bernorig vaccine-induced rabies cases strongly suggests that only one replication-competent haplotype was present in those cases, despite a higher number of variants. One reason could be the passage history. Both ERA-BHK21 and SAD Bernorig vaccine virus strains descend from the same derivative of the cell culture-adapted field virus strain, namely SAD [5,6], but were independently further developed by multiple in vitro passaging on various cell types (see [14] for a more comprehensive overview on this topic). Genetic as well as phylogenetic analyses, however, suggest that the ERA-BHK21 vaccine virus remained closer to its progenitor SAD strain [14,15]. On the other hand, the SAD Bernorig vaccine virus strain was further developed into the commercial vaccine virus strains SAD Bern and SAD B19. Eventually, increased cell culture adaptation of SAD Bern and SAD B19 vaccine virus strains may have selected for only a single replication-competent haplotype. In parallel, this adaptation towards efficient replication in cell culture is favoring the development of defective interfering particles [34] that provide the background for the high genetic diversity as found by the variant analyses (Figure 1). Unfortunately, our analysis was restricted to viruses for which sufficient sequencing data was available, and thus, could not cover the SAD P5/88 vaccine-induced rabies cases [15].
As attenuated live vaccines consist of replication-competent virus particles, a major concern in their application is the potential reversion to virulence as demonstrated by vaccines for the severe acute respiratory syndrome coronavirus (SARS-CoV) [35], peste des petits ruminants virus (PPRV) [36], the avian metapneumovirus (aMPV) [37] and Rift Valley fever virus [38]. For the attenuated oral rabies vaccines investigated here, our analyses on alterations in antigenic sites showed no evidence of a reversion to virulence of viruses from vaccine-induced rabies cases either as result of the production process nor of the in vivo replication in affected animals. This is also supported by the population-based analysis conducted, which clearly separated field strains from vaccines and their induced cases (Figure 3). Rather, other factors including the individual immune status or genetic predisposition of the infected animals most likely have contributed to the development of disease. This is supported by the observation, that none of these sporadic cases resulted in onward transmission or was of epidemiological relevance [1,8].
The results of this study clearly demonstrate the benefits of utilizing deep sequencing followed by a comprehensive stepwise data analysis focusing on consensus, population, and variant level. While each of the analyses may lead to biased conclusions, only the combination allows for a holistic assessment, and thus should be standard in similar scenarios.

Supplementary Materials

Supplementary materials can be found at https://www.mdpi.com/1999-4915/12/1/115/s1.

Author Contributions

Conceptualization, D.H.; data curation, S.C.; formal analysis, S.C.; investigation, S.C., P.T., C.F.-G., S.N.-D., D.H., A.O., and M.S.; methodology, D.H.; resources, T.M.; supervision, T.M., D.H., and M.S.; visualization, D.H.; writing—original draft, S.C., C.M.F., T.M., A.O., and M.S.; writing—review and editing, S.C., P.T., C.M.F., T.M., S.N.-D., D.H., A.O., and M.S. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by an intramural collaborative research grant on Lyssaviruses at the Friedrich-Loeffler-Institut (Rie-0373).

Acknowledgments

We would like to thank Patrick Zitzow (FLI) for excellent technical support.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Müller, T.; Bätza, H.-J.; Beckert, A.; Bunzenthal, C.; Cox, J.H.; Freuling, C.M.; Fooks, A.R.; Frost, J.; Geue, L.; Hoeflechner, A.; et al. Analysis of vaccine-virus-associated rabies cases in red foxes (Vulpes vulpes) after oral rabies vaccination campaigns in Germany and Austria. Arch. Virol. 2009, 154, 1081–1091. [Google Scholar] [CrossRef] [PubMed]
  2. Freuling, C.M.; Hampson, K.; Selhorst, T.; Schröder, R.; Meslin, F.X.; Mettenleiter, T.C.; Müller, T. The elimination of fox rabies from Europe: Determinants of success and lessons for the future. Philos. Trans. R. Soc. B Biol. Sci. 2013, 368, 20120142. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  3. MacInnes, C.D.; Smith, S.M.; Tinline, R.R.; Ayers, N.R.; Bachmann, P.; Ball, D.G.; Calder, L.A.; Crosgrey, S.J.; Fielding, C.; Hauschildt, P.; et al. Elimination of rabies from red foxes in eastern Ontario. J. Wildl. Dis. 2001, 37, 119–132. [Google Scholar] [CrossRef] [PubMed]
  4. Rupprecht, C.E.; Plotkin, S.A. Rabies vaccines. In Vaccines, 6th ed.; Offit, P.A., Plotkin, S.A., Orenstein, W.A., Eds.; Elsevier Saunders: Saint Louis, MO, USA, 2013; pp. 646–668. ISBN 9781455700905. [Google Scholar]
  5. Abelseth, M.K. An Attenuated Rabies Vaccine for Domestic Animals Produced in Tissue Culture. Can. Vet. J. 1964, 5, 279–286. [Google Scholar]
  6. Abelseth, M.K. Propagation of Rabies Virus in Pig Kidney Cell Culture. Can. Vet. J. 1964, 5, 84–87. [Google Scholar]
  7. Fehlner-Gardiner, C.; Nadin-Davis, S.; Armstrong, J.; Muldoon, F.; Bachmann, P.; Wandeler, A. Era vaccine-derived cases of rabies in wildlife and domestic animals in Ontario, Canada, 1989–2004. J. Wildl. Dis. 2008, 44, 71–85. [Google Scholar] [CrossRef] [Green Version]
  8. Müller, T.F.; Schröder, R.; Wysocki, P.; Mettenleiter, T.C.; Freuling, C.M. Spatio-temporal Use of Oral Rabies Vaccines in Fox Rabies Elimination Programmes in Europe. PLoS Negl. Trop. Dis. 2015, 9, e0003953. [Google Scholar] [CrossRef] [Green Version]
  9. Artois, M.; Guittré, C.; Thomas, I.; Leblois, H.; Brochier, B.; Barrat, J. Potential pathogenicity for rodents of vaccines intended for oral vaccination against rabies: A comparison. Vaccine 1992, 10, 524–528. [Google Scholar] [CrossRef]
  10. Vos, A.; Neubert, A.; Aylan, O.; Schuster, P.; Pommerening, E.; Müller, T.; Chivatsi, D.C. An update on safety studies of SAD B19 rabies virus vaccine in target and non-target species. Epidemiol. Infect. 1999, 123, 165–175. [Google Scholar] [CrossRef]
  11. Schneider, L.; Cox, J.H. Ein Feldversuch zur oralen Immunisierung von Füchsen gegen die Tollwut in der Bundesrepubik Deutschland: Unschädlichkeit, Wirksamkeit und Stabilität der Vakzine SAD B 19. Tierarztliche Umschau 1983, 1983, 315–324. [Google Scholar]
  12. Güzel, T.; Aylan, O.; Vos, A. Innocuity tests with SAD B19 in Turkish non-target species. J. Etlik Vet. Microbiol. 1998, 1998, 103–112. [Google Scholar]
  13. Vuta, V.; Picard-Meyer, E.; Robardet, E.; Barboi, G.; Motiu, R.; Barbuceanu, F.; Vlagioiu, C.; Cliquet, F. Vaccine-induced rabies case in a cow (Bos taurus): Molecular characterisation of vaccine strain in brain tissue. Vaccine 2016, 34, 5021–5025. [Google Scholar] [CrossRef] [PubMed]
  14. Höper, D.; Freuling, C.M.; Müller, T.; Hanke, D.; von Messling, V.; Duchow, K.; Beer, M.; Mettenleiter, T.C. High definition viral vaccine strain identity and stability testing using full-genome population data--The next generation of vaccine quality control. Vaccine 2015, 33, 5829–5837. [Google Scholar] [CrossRef] [PubMed]
  15. Pfaff, F.; Müller, T.; Freuling, C.M.; Fehlner-Gardiner, C.; Nadin-Davis, S.; Robardet, E.; Cliquet, F.; Vuta, V.; Hostnik, P.; Mettenleiter, T.C.; et al. In-depth genome analyses of viruses from vaccine-derived rabies cases and corresponding live-attenuated oral rabies vaccines. Vaccine 2019, 37, 4758–4765. [Google Scholar] [CrossRef]
  16. Wylezich, C.; Papa, A.; Beer, M.; Höper, D. A Versatile Sample Processing Workflow for Metagenomic Pathogen Detection. Sci. Rep. 2018, 8, 13108. [Google Scholar] [CrossRef] [Green Version]
  17. Nadin-Davis, S.A.; Colville, A.; Trewby, H.; Biek, R.; Real, L. Application of high-throughput sequencing to whole rabies viral genome characterisation and its use for phylogenetic re-evaluation of a raccoon strain incursion into the province of Ontario. Virus Res. 2017, 232, 123–133. [Google Scholar] [CrossRef]
  18. Hostnik, P.; Picard-Meyer, E.; Rihtarič, D.; Toplak, I.; Cliquet, F. Vaccine-induced rabies in a red fox (Vulpes vulpes): Isolation of vaccine virus in brain tissue and salivary glands. J. Wildl. Dis. 2014, 50, 397–401. [Google Scholar] [CrossRef]
  19. Brzózka, K.; Finke, S.; Conzelmann, K.-K. Inhibition of interferon signaling by rabies virus phosphoprotein P: Activation-dependent binding of STAT1 and STAT2. J. Virol. 2006, 80, 2675–2683. [Google Scholar] [CrossRef] [Green Version]
  20. Faber, M.; Faber, M.-L.; Papaneri, A.; Bette, M.; Weihe, E.; Dietzschold, B.; Schnell, M.J. A single amino acid change in rabies virus glycoprotein increases virus spread and enhances virus pathogenicity. J. Virol. 2005, 79, 14141–14148. [Google Scholar] [CrossRef] [Green Version]
  21. Gholami, A.; Kassis, R.; Real, E.; Delmas, O.; Guadagnini, S.; Larrous, F.; Obach, D.; Prevost, M.-C.; Jacob, Y.; Bourhy, H. Mitochondrial dysfunction in lyssavirus-induced apoptosis. J. Virol. 2008, 82, 4774–4784. [Google Scholar] [CrossRef] [Green Version]
  22. Kgaladi, J.; Wright, N.; Coertse, J.; Markotter, W.; Marston, D.; Fooks, A.R.; Freuling, C.M.; Müller, T.F.; Sabeta, C.T.; Nel, L.H. Diversity and epidemiology of Mokola virus. PLoS Negl. Trop. Dis. 2013, 7, e2511. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  23. Langevin, C.; Tuffereau, C. Mutations conferring resistance to neutralization by a soluble form of the neurotrophin receptor (p75NTR) map outside of the known antigenic sites of the rabies virus glycoprotein. J. Virol. 2002, 76, 10756–10765. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  24. Lentz, T.L.; Wilson, P.T.; Hawrot, E.; Speicher, D.W. Amino acid sequence similarity between rabies virus glycoprotein and snake venom curaremimetic neurotoxins. Science 1984, 226, 847–848. [Google Scholar] [CrossRef] [PubMed]
  25. Masatani, T.; Ito, N.; Shimizu, K.; Ito, Y.; Nakagawa, K.; Abe, M.; Yamaoka, S.; Sugiyama, M. Amino acids at positions 273 and 394 in rabies virus nucleoprotein are important for both evasion of host RIG-I-mediated antiviral response and pathogenicity. Virus Res. 2011, 155, 168–174. [Google Scholar] [CrossRef]
  26. Mita, T.; Shimizu, K.; Ito, N.; Yamada, K.; Ito, Y.; Sugiyama, M.; Minamoto, N. Amino acid at position 95 of the matrix protein is a cytopathic determinant of rabies virus. Virus Res. 2008, 137, 33–39. [Google Scholar] [CrossRef]
  27. Préhaud, C.; Wolff, N.; Terrien, E.; Lafage, M.; Mégret, F.; Babault, N.; Cordier, F.; Tan, G.S.; Maitrepierre, E.; Ménager, P.; et al. Attenuation of rabies virulence: Takeover by the cytoplasmic domain of its envelope protein. Sci. Signal 2010, 3, ra5. [Google Scholar] [CrossRef]
  28. Rieder, M.; Brzózka, K.; Pfaller, C.K.; Cox, J.H.; Stitz, L.; Conzelmann, K.-K. Genetic dissection of interferon-antagonistic functions of rabies virus phosphoprotein: Inhibition of interferon regulatory factor 3 activation is important for pathogenicity. J. Virol. 2011, 85, 842–852. [Google Scholar] [CrossRef] [Green Version]
  29. Takayama-Ito, M.; Ito, N.; Yamada, K.; Sugiyama, M.; Minamoto, N. Multiple amino acids in the glycoprotein of rabies virus are responsible for pathogenicity in adult mice. Virus Res. 2006, 115, 169–175. [Google Scholar] [CrossRef]
  30. Tian, D.; Luo, Z.; Zhou, M.; Li, M.; Yu, L.; Wang, C.; Yuan, J.; Li, F.; Tian, B.; Sui, B.; et al. Critical Role of K1685 and K1829 in the Large Protein of Rabies Virus in Viral Pathogenicity and Immune Evasion. J. Virol. 2016, 90, 232–244. [Google Scholar] [CrossRef] [Green Version]
  31. Tuffereau, C.; Leblois, H.; Bénéjean, J.; Coulon, P.; Lafay, F.; Flamand, A. Arginine or lysine in position 333 of ERA and CVS glycoprotein is necessary for rabies virulence in adult mice. Virology 1989, 172, 206–212. [Google Scholar] [CrossRef]
  32. Wirblich, C.; Tan, G.S.; Papaneri, A.; Godlewski, P.J.; Orenstein, J.M.; Harty, R.N.; Schnell, M.J. PPEY motif within the rabies virus (RV) matrix protein is essential for efficient virion release and RV pathogenicity. J. Virol. 2008, 82, 9730–9738. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  33. Forró, B.; Marton, S.; Kecskeméti, S.; Hornyák, Á.; Bányai, K. Vaccine-associated rabies in red fox, Hungary. Vaccine 2019. [Google Scholar] [CrossRef] [PubMed]
  34. Rezelj, V.V.; Levi, L.I.; Vignuzzi, M. The defective component of viral populations. Curr. Opin. Virol. 2018, 33, 74–80. [Google Scholar] [CrossRef] [PubMed]
  35. Jimenez-Guardeño, J.M.; Regla-Nava, J.A.; Nieto-Torres, J.L.; DeDiego, M.L.; Castaño-Rodriguez, C.; Fernandez-Delgado, R.; Perlman, S.; Enjuanes, L. Identification of the Mechanisms Causing Reversion to Virulence in an Attenuated SARS-CoV for the Design of a Genetically Stable Vaccine. PLoS Pathog. 2015, 11, e1005215. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  36. Eloiflin, R.-J.; Boyer, M.; Kwiatek, O.; Guendouz, S.; Loire, E.; Servan de Almeida, R.; Libeau, G.; Bataille, A. Evolution of Attenuation and Risk of Reversal in Peste des Petits Ruminants Vaccine Strain Nigeria 75/1. Viruses 2019, 11, 724. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  37. Franzo, G.; Naylor, C.J.; Drigo, M.; Croville, G.; Ducatez, M.F.; Catelli, E.; Laconi, A.; Cecchinato, M. Subpopulations in aMPV vaccines are unlikely to be the only cause of reversion to virulence. Vaccine 2015, 33, 2438–2441. [Google Scholar] [CrossRef]
  38. Lokugamage, N.; Ikegami, T. Genetic stability of Rift Valley fever virus MP-12 vaccine during serial passages in culture cells. NPJ Vaccines 2017, 2. [Google Scholar] [CrossRef]
Viruses 12 00115 g001 550
Figure 1. Population-based analysis of ERA-BHK21 (left) and SAD Bernorig-derived vaccine batches (B/…) (right) and viruses of their vaccine-induced rabies cases (C/…) displayed by a fitted pairwise Manhattan distances plot. Data sets of samples that were added in this study (Table 1 and Table 2) are highlighted in bold, whereas data sets of vaccine batches and their related induced rabies cases are displayed in different colors: blue—SAD B19 vaccine batches; green—SAD Bern vaccine batches, dark blue—SAD Bernorig-derived vaccine-induced rabies cases; orange—ERA-BHK21 vaccine batch; and dark orange—ERA-derived vaccine-induced rabies cases (orange/dark orange). For the ERA vaccine-induced rabies cases, three distinct combinations of single nucleotide variants were selected that originate from the progenitor ERA-BHK21 vaccine (see Section 3.2), forming groups of samples with specific differences that are represented by the colors of the respective sample names (1–3).
Figure 1. Population-based analysis of ERA-BHK21 (left) and SAD Bernorig-derived vaccine batches (B/…) (right) and viruses of their vaccine-induced rabies cases (C/…) displayed by a fitted pairwise Manhattan distances plot. Data sets of samples that were added in this study (Table 1 and Table 2) are highlighted in bold, whereas data sets of vaccine batches and their related induced rabies cases are displayed in different colors: blue—SAD B19 vaccine batches; green—SAD Bern vaccine batches, dark blue—SAD Bernorig-derived vaccine-induced rabies cases; orange—ERA-BHK21 vaccine batch; and dark orange—ERA-derived vaccine-induced rabies cases (orange/dark orange). For the ERA vaccine-induced rabies cases, three distinct combinations of single nucleotide variants were selected that originate from the progenitor ERA-BHK21 vaccine (see Section 3.2), forming groups of samples with specific differences that are represented by the colors of the respective sample names (1–3).
Viruses 12 00115 g001
Viruses 12 00115 g002 550
Figure 2. Schematic illustration of all variants found for the ERA vaccine-induced rabies cases as well as nucleotide exchanges between viruses from vaccine-induced rabies cases and the ERA-BHK21 vaccine virus strain. Colored dotted lines display differences at the consensus level for viruses from vaccine-induced rabies cases that derived from selected combinations of single nucleotide variants (SNVs) found in the progenitor ERA-BHK21 vaccine virus (group 1–3; Table 4, for variant frequencies, see Table S7, Supplementary File). Nucleotide bases illustrated on the genome sequence represent spontaneous mutations found for viruses of vaccine-induced rabies cases. Each of these three groups were characterized by a unique set of differences at the consensus level (Table 4). The only exception was the sample C/CAN/1994N35116 which lacked one difference (position 3587, Table S7) and had one additional difference in close proximity that cannot be found in any of the other samples (position 3734, Table S7).
Figure 2. Schematic illustration of all variants found for the ERA vaccine-induced rabies cases as well as nucleotide exchanges between viruses from vaccine-induced rabies cases and the ERA-BHK21 vaccine virus strain. Colored dotted lines display differences at the consensus level for viruses from vaccine-induced rabies cases that derived from selected combinations of single nucleotide variants (SNVs) found in the progenitor ERA-BHK21 vaccine virus (group 1–3; Table 4, for variant frequencies, see Table S7, Supplementary File). Nucleotide bases illustrated on the genome sequence represent spontaneous mutations found for viruses of vaccine-induced rabies cases. Each of these three groups were characterized by a unique set of differences at the consensus level (Table 4). The only exception was the sample C/CAN/1994N35116 which lacked one difference (position 3587, Table S7) and had one additional difference in close proximity that cannot be found in any of the other samples (position 3734, Table S7).
Viruses 12 00115 g002
Viruses 12 00115 g003 550
Figure 3. Population-based analysis of ERA-BHK21- and SAD-Bernorig-derived vaccine strains and viruses of their vaccine-induced rabies cases in relation to selected field RABV strains (F/…).
Figure 3. Population-based analysis of ERA-BHK21- and SAD-Bernorig-derived vaccine strains and viruses of their vaccine-induced rabies cases in relation to selected field RABV strains (F/…).
Viruses 12 00115 g003
Table
Table 1. Overview of additional vaccine-induced rabies cases from Poland and the Ontario province of Canada analyzed in this study.
Table 1. Overview of additional vaccine-induced rabies cases from Poland and the Ontario province of Canada analyzed in this study.
Sample IDAnimalOriginDate of FindingDistributed VaccineReferenceSample MaterialSequencing MethodRABV ReadsMean Depth
C/POL/2017Red foxMałopolska, PolandDec 2017SAD BernThis StudyBrain18286170
C/CAN/1992N8944Red foxCarleton, Canada1992ERA[7]Brain249,6101138
C/CAN/1992N8991RaccoonTyendinaga, Canada1992ERA[7]Brain250,5871092
C/CAN/1994N6762-BCowGloucester, Canada1994ERA[7]Brain251,2251129
C/CAN/1994N6762-MMouseGloucester, Canada1994ERA[7]Brain250,4031165
C/CAN/1994N35116Striped skunkHullet, Canada1994ERA[7]Brain249,3191127
C/CAN/1996N5367Red foxOps, Canada1996ERA[7]Brain250,1881090
C/CAN/1998N4916RaccoonNormanby, Canada1998ERA[7]Brain251,2391092
C/CAN/1999N6567Red foxMulmur, Canada1999ERA[7]Brain251,411995
Table
Table 2. Overview of ERA and Street Alabama Dufferin (SAD) Bern vaccine batches analyzed in this study.
Table 2. Overview of ERA and Street Alabama Dufferin (SAD) Bern vaccine batches analyzed in this study.
Sample IDVaccine StrainBatch IDDescriptionSequencing MethodRABV ReadsMean Depth
B/ERA/lot16ERAlot 016ERA-BHK21 vaccine batch from 2005250,78811,301
B/SAD/Bern/4024SAD Bern17/61Lj (4024)Lysvulpen vaccine batches that were distributed in the area of Małopolska in autumn 2016, spring 2017 and autumn 2017151,0981317
B/SAD/Bern/7024SAD Bern17/76Lj (7024)152,0011333
Table
Table 3. List of additional R packages used for the calculation and graphical representation of the Manhattan distances.
Table 3. List of additional R packages used for the calculation and graphical representation of the Manhattan distances.
R PackageVersionDescription
ape5.3Analyses of phylogenetics and evolution
cluster2.0.8Extended cluster analysis
fpc2.2-3Flexible procedures for clustering
ggplot23.2.0Data visualizations
ggrepel0.8.1Non-overlapping text labels for ggplot2
grid3.6.0Grid graphics package
phyloseq1.28.0Handling and analysis of high-throughput microbiome census data
vegan2.5-5Community ecology package
Table
Table 4. Groups of viruses from ERA vaccine-induced rabies cases showing similar patterns of nucleotide exchanges compared to the progenitor ERA-BHK21 vaccine virus strain at the consensus level.
Table 4. Groups of viruses from ERA vaccine-induced rabies cases showing similar patterns of nucleotide exchanges compared to the progenitor ERA-BHK21 vaccine virus strain at the consensus level.
SamplesHost SpeciesNumber of Group-Specific Differences
Group 1C/CAN/1992N8944Red fox9
C/CAN/1992N8991Raccoon
C/CAN/1996N5367Red fox
C/CAN/1996
Group 2C/CAN/1994N35116Striped skunk7
C/CAN/1994N6762BCow
C/CAN/1994N6762MCow (M-passage) 1
C/CAN/1994
C/CAN/1998N4916Raccoon
Group 3C/CAN/1991Striped skunk3
C/CAN/1999N6567Red fox
1 Mouse-passaged virus originating from C/CAN/1994N6762-B.
Table
Table 5. Viruses of ERA vaccine-induced rabies cases that maintained the vaccine consensus in the form of an SNV.
Table 5. Viruses of ERA vaccine-induced rabies cases that maintained the vaccine consensus in the form of an SNV.
Sample IDPosition *Consensus Vaccine-Induced CaseVariant Vaccine-Induced Case (Consensus Vaccine)Frequency of SNVs in the Vaccine-Induced Case
CAN/1992N894410,840AG33.7%
CAN/1994N6762-M3085TC31.9%
* Nucleotide position relative to the sequence of the ERA-BHK21 vaccine batch B/ERA/lot16 full genome.

Share and Cite

MDPI and ACS Style

Calvelage, S.; Smreczak, M.; Orłowska, A.; Freuling, C.M.; Müller, T.; Fehlner-Gardiner, C.; Nadin-Davis, S.; Höper, D.; Trębas, P. Population- and Variant-Based Genome Analyses of Viruses from Vaccine-Derived Rabies Cases Demonstrate Product Specific Clusters and Unique Patterns. Viruses 2020, 12, 115. https://doi.org/10.3390/v12010115

AMA Style

Calvelage S, Smreczak M, Orłowska A, Freuling CM, Müller T, Fehlner-Gardiner C, Nadin-Davis S, Höper D, Trębas P. Population- and Variant-Based Genome Analyses of Viruses from Vaccine-Derived Rabies Cases Demonstrate Product Specific Clusters and Unique Patterns. Viruses. 2020; 12(1):115. https://doi.org/10.3390/v12010115

Chicago/Turabian Style

Calvelage, Sten, Marcin Smreczak, Anna Orłowska, Conrad Martin Freuling, Thomas Müller, Christine Fehlner-Gardiner, Susan Nadin-Davis, Dirk Höper, and Paweł Trębas. 2020. "Population- and Variant-Based Genome Analyses of Viruses from Vaccine-Derived Rabies Cases Demonstrate Product Specific Clusters and Unique Patterns" Viruses 12, no. 1: 115. https://doi.org/10.3390/v12010115

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