Antigenic characterisation of lyssaviruses in South Africa

Copyright: © 2014. The Authors. Licensee: AOSIS OpenJournals. This work is licensed under the Creative Commons Attribution License. There are at least six Lyssavirus species that have been isolated in Africa, which include classical rabies virus, Lagos bat virus, Mokola virus, Duvenhage virus, Shimoni bat virus and Ikoma lyssavirus. In this retrospective study, an analysis of the antigenic reactivity patterns of lyssaviruses in South Africa against a panel of 15 anti-nucleoprotein monoclonal antibodies was undertaken. A total of 624 brain specimens, collected between 2005 and 2009, confirmed as containing lyssavirus antigen by direct fluorescent antibody test, were subjected to antigenic differentiation. The lyssaviruses were differentiated into two species, namely rabies virus (99.5%) and Mokola virus (0.5%). Furthermore, rabies virus was further delineated into two common rabies biotypes in South Africa: canid and mongoose. Initially, it was found that the canid rabies biotype had two reactivity patterns; differential staining was observed with just one monoclonal antibody. This difference was likely to have been an artefact related to sample quality, as passage in cell culture restored staining. Mongoose rabies viruses were more heterogeneous, with seven antigenic reactivity patterns detected. Although Mokola viruses were identified in this study, prevalence and reservoir host species are yet to be established. These data demonstrate the usefulness of monoclonal antibody typing panels in lyssavirus surveillance with reference to emergence of new species or spread of rabies biotypes to new geographic zones.


Introduction
Rabies is a viral zoonotic disease caused by lyssaviruses that are capable of infecting all mammalian species.These viruses belong to the Lyssavirus genus of the Rhabdoviridae family.The Lyssavirus genus consists of 14 recognised viral species based on serological assays, genetic distances of the nucleoprotein gene, topology and consistency of the phylogenetic trees, antigenic patterns in reaction with anti-nucleocapsid monoclonal antibodies and/or other additional characters such as ecological properties, host and geographic range and pathological features (Dietzgen et al. 2011;ICTV 2014).These includes classical rabies virus (RABV) and the other viral species are referred to as rabies-related viruses.Of the 14 species, only six have been identified and isolated in Africa to date, namely: RABV, Lagos bat virus (LBV), Mokola virus (MOKV), Duvenhage virus (DUVV), Shimoni bat virus (SHIBV) and Ikoma lyssavirus (IKOV).In addition, Lleida bat lyssavirus was recently proposed as a putative species of the Lyssavirus genus (Freuling et al. 2011;Marston et al. 2012).Within the classical rabies virus two common variants have been identified in southern Africa: the canid rabies and mongoose rabies biotypes (King, Meredith & Thomson 1993;Nel et al. 2005;Von Teichman et al. 1995).Rabies in mongooses was first identified in southern Africa in the 1800s (Snyman 1940;Swanepoel 2004).In contrast, rabies in dogs appears to have spread into southern Africa via the north-west from Angola in the 1940s, and by the 1950s dog rabies had entered Zimbabwe and South Africa, where it became established in various canid host species in subsequent years (Nadin-Davis & Bingham 2004).In South Africa, rabies still remains a significant public and veterinary health threat, as the disease is well maintained by various host species.
Genetic typing by nucleotide sequence analysis and antigenic typing using discriminatory monoclonal antibody (mAb) panels in an indirect immunofluorescence assay (IFA) can be used to differentiate Lyssavirus species, as well as variants within the species (Favoretto et al. 2002;Favoretto et al. 2006;Smith 2002).The discriminatory capability of IFA is generally determined by the number of mAbs used in a panel (Nadin-Davis et al. 2010;Smith 2002).Antigenic typing panels generally depend on mAbs specific for the nucleoprotein (N) and phosphoprotein (P); glycoprotein (G) mAbs can also be used.In addition, specific mAb panels have also been generated and used to map antigenic sites on these proteins to characterise lyssaviruses, and also to distinguish RABV variants persisting in geographically restricted host reservoirs (Nadin-Davis et al. 2010;Rupprecht et al. 1991;Smith & King 1996).Lyssavirus typing with mAbs provides a rapid and inexpensive screening method for large-scale surveillance of both new lyssaviruses (Gould et al. 1998) and common or emerging variants of the rabies virus (Nadin-

Antigenic characterisation of lyssaviruses in South Africa
Davis et al. 2000).This technology was used to differentiate between street virus and the rabies vaccine strain during the oral vaccination campaigns against fox rabies in Western Europe and in Ontario, Canada (Fehlner-Gardiner et al. 2008;Schneider et al. 1988).In previous studies, a panel of 80 anti-nucleoprotein mAbs was used to differentiate rabies viruses from Namibia and South Africa (King et al. 1993).It was shown that canid rabies viruses conformed to a single reactivity pattern, whereas mongoose rabies viruses had several variations in their reactivity patterns (King et al. 1993).This retrospective study describes the mAb typing of lyssavirus isolates collected in South Africa from different animal species between 2005 and 2009 is reported on.

Brain specimens
The specimens were routinely submitted to the Agricultural Research Council-Ondestepoort Veterinary Institute (ARC-OVI) for rabies diagnosis from various provinces of South Africa except the Western Cape and Kwazulu-Natal (KZN), which submit samples to Allerton Veterinary Laboratory, KZN and Western Cape provincial veterinary laboratory (Figure 1).All the samples were confirmed to be positive for lyssavirus antigen by direct fluorescent antibody test (dFAT) using fluorescently labelled goat antiserum raised against purified RABV and MOKV ribonucleoproteins (Dean, Abelseth & Atanasiu 1996).A total of 624 brain specimens, collected between 2005 and 2009, from different animal species were analysed by mAb-IFA.

Viruses included in the analysis
Further phylogenetic analyses and comparison with phylogroups determined from previous studies (Nel et al. 2005;Van Zyl, Markotter & Nel 2010) were carried out on a panel of 21 representative isolates selected from the seven mongoose biotype groups detected with antigenic typing and from six isolates selected from the canid biotype groups.The year of isolation, host species and geographical origin of Page 2 of 9

Lagos bat virus Mokola virus
Duvenhage virus Note: -, no reactivity; +, reactivity observed; variable, some isolates within the species or biotype react with the mAb and others do not react.
the mongoose biotype isolates are shown in Table 2.

Viral RNA extraction and reverse transcription polymerase chain reaction
Total viral RNA was extracted from original brain material using Trizol ® reagent (Sigma-Aldrich, USA) according to the manufacturer's instructions.For the cDNA synthesis of the G-L intergenic region and nucleoprotein genes, 5 µL of total RNA (approximately 2 µg) and 20 pmol of the positive sense primers (2 µL), 001lys (N-gene) and G (+) (G-L intergenic region) (Table 3) were denatured at 65 °C for 5 min.The reaction mixtures were immediately cooled on ice and RNA was reverse transcribed at 42 °C for 60 min in a 20 µL reaction mixture containing 200 U Superscript III reverse transcriptase (Invitrogen), 4 µL of reaction buffer, 10 mM deoxynucleotide triphosphates (dNTPs) and 40 U of RNase inhibitor (Promega, USA).The polymerase chain reactions (PCRs) were performed in a total volume of 50 µL containing 1X PCR reaction buffer, 1.5 mM MgCl 2 , 250 µM of each dNTP, 5 µL of cDNA, 1.25 U of Taq DNA polymerase (Applied Biosystem, USA) and 40 pmol of each of the primers (G [+] and L [-], 001lys and 304) (Table 3).The amplification reactions were performed in a Geneamp 2400 thermocycler with an initial denaturation of 94 °C for 2 min, followed by 40 cycles (25 cycles) of denaturation at 94 °C, primer annealing at 37 °C for 30 s (42 °C for 90 s) and primer extension at 72 °C for 90 s followed by a final extension at 72 °C for 7 min (values in parenthesis indicate conditions for G-L intergenic region).The amplicons were subsequently analysed by ethidium bromide-stained agarose gel (1%) electrophoresis.

DNA purification, nucleotide sequencing and phylogenetic analysis
The PCR products were purified using the DNA purification kit according to the manufacturer's instructions (Promega, USA).The purified products were sequenced in both directions using both primer sets used in the reverse transcription polymerase chain reaction (RT-PCR) step together with the sequencing kit according to the manufacturer's instructions (Applied Biosystems, USA).The sequencing products were electrophoresed on an ABI 377 automated sequencer (Applied Biosystems) and the consensus sequences of 1353 nucleotides (nt) of the N gene and 592 nt (G-L intergenic region) were obtained by using MEGA 3.1 sequence editor (Kumar, Tamura & Nei 2004).The phylogenetic trees were constructed using the neighbour joining (NJ) method in Clustal X (Saitou & Nei 1987).Two sets of 18 full-length nucleoprotein and partial G-L intergenic region nucleotide sequences representing previous mongoose RABV phylogroups from southern Africa were retrieved from Genbank and included for phylogenetic analysis.The branching order of the trees was statistically evaluated by 1000 bootstrap replicates; values of more than 700 (70%) were considered to be statistically significant (Hills & Bull 1993).
The canid RABV biotype isolates fell into two categories, with differences observed only with mAb 32GD12 (Table 4).The majority of the canid viruses (403/438, 92%) were stained by this mAb (reactivity pattern I).Although viruses that conformed to canid biotype reactivity pattern II were less common, they exhibited similar geographic   distribution to viruses that conformed to canid biotype reactivity pattern I (Figure 1).Among the mongoose RABV biotypes, seven antigenic reactivity patterns (A-G) were detected (Table 5) and these viruses were recovered from a wide range of different animal species (Table 6).Mongoose rabies biotype viruses that conformed to group B and D staining patterns were most common (Table 5).However, the antigenic reactivity patterns do not appear to have any species-specific association with the different hosts (Table 6).

Galerella sanguinea
The PCR amplicons generated from the selected isolates were of the expected size of approximately 850 bp (G-L region) and 1680 bp (N-gene).Sequencing yielded about 740 bases (G-L region) and 1400 bases (N-gene) of sequence on average, which were trimmed to a consensus of 592 and 1353 nucleotides (nt) for phylogenetic analysis (data not shown).Representative canid RABV isolates conforming to both reactivity patterns (n = 3 for each group) had 98% nucleotide sequence identity in the nucleoprotein-encoding region (data not shown).
For the mongoose biotype, phylogenetic analysis of representative isolates from each group (n = 21) revealed trees with similar topologies when analysing either the  partial G-L intergenic region or the full-length nucleoprotein nucleotide sequences (Figure 2a and Figure 2b).Isolates corresponding to staining pattern A (n = 3) all clustered within the previously described phylogroup IV, whereas those with staining pattern E (n = 3) and F (n = 3) clustered within phylogroup III.These isolates also clustered geographically with other isolates that fell within phylogroups III and IV (Figure 2).In contrast, the isolates corresponding to the other staining patterns (B, C, D, F and G) did not group together under a single phylogroup, and did not exhibit geographical clustering, as observed for isolates with staining patterns A, E and F. Thus, two of the three isolates with staining patterns B, C and G clustered with phylogroup III, whilst the third isolate examined fell within phylogroup II and IV.The three viruses with staining pattern D fell within phylogroups II, III and IV.Whilst the isolates from staining pattern groups B, C, D, and G did not cluster geographically, they did map to similar geographic areas with other isolates belonging to the same phylogroups (Figure 2).

Discussion
The geographical distribution of the rabies cases detected between 2005 and 2009 reflects the distribution of the susceptible host species in South Africa.For example, bat-eared foxes maintain and transmit rabies in the western regions of the country (Sabeta et al. 2007a), whereas blackbacked jackals maintain the disease cycles in the northern region of South Africa (Swanepoel 2004;Zulu, Sabeta & Nel 2009).Canid rabies biotype cases were detected in all provinces included in the present study, and although the provinces of KZN and Western Cape were not represented, the canid RABV biotype is known to be endemic in these regions, with principal reservoirs in domestic dogs and bateared foxes, respectively (Swanepoel 2004).The majority of mongoose RABV biotype cases were found in North-West, Free State and Mpumalanga provinces, which is consistent with the distribution of Herpestidae species (Nel et al. 2005).Two mongoose RABV biotype variants were detected in a domestic dog and black-backed jackal in Limpopo province, whereas all other cases from this region were of the canid RABV biotype.These detections may have been the result of sporadic introduction from neighbouring Zimbabwe and Botswana, where mongoose biotype RABV is endemic (Bingham et al. 2001;Foggin 1988;Sabeta et al. 2008;Tremlett , Wibberley & King 1994).However, as the sampling method was only passive, relying on the state authorities to submit animal rabies suspected cases to the laboratory, it is possible that the mongoose RABV biotype is also established in this region of South Africa, but its prevalence  underestimated.Spill-over of the mongoose rabies biotype from the host species is common; for example, mongoose RABVs were recovered from domestic dogs, sheep and cattle in the present study.In contrast, spillover of the canid RABV biotype into non-canid wildlife species is rarely observed (Von Teichman et al. 1995).Understanding the movement of RABV biotype and particularly spillover of wildlife variants into domestic animals and vice versa is a key aspect in the rabies control in the region.
The mAb panel used in this study was originally selected to differentiate amongst canid and mongoose biotypes of the rabies virus, as well as the non-rabies Lyssavirus species of MOKV, LBV and DUVV circulating in the region (Bingham et al. 2001).Several different staining patterns were observed for the rabies viruses in the present study, notably two for the canid biotype and seven for the mongoose biotype.A similar variation in mAb reactivity patterns was observed amongst the mongoose viruses in previous studies (King et al. 1993;Tremlett et al. 1994).In contrast, a single antigenic reactivity pattern of the canid RABV biotype was previously observed in isolates from South Africa and Namibia (King et al. 1993).
In previous phylogenetic studies, whilst seven canid RABV biotype variants or clusters were observed in the Northern provinces of South Africa, they conformed to a single antigenic typing pattern (Agricultural Research Council-Onderstepoort Veterinary Institute 2007;Zulu 2007).Similar observations were made in a study from Botswana, where antigenic typing could not differentiate smaller clades or subgroupings established by phylogenetic analyses (Johnson et al. 2004).It is unlikely that the detection of two reactivity patterns for the canid RABV biotype in the present study is due to a superior performance of the current mAb panel distinguishing between two variants that have long co-existed in the region.However, if this scenario were the case, observation of some geographic segregation of the variants would be expected, and it was not.Given that the observed canid RABV reactivity patterns differed by only one mAb (32GD12), it is more likely that the epitope recognised by mAb 32GD12 is unstable and the observed staining variability was an artefact of the quality of the tissue that was sampled.When the five canid rabies viruses that had antigenic reactivity pattern II were grown in cell culture, reactivity with mAb 32GD12 was restored, suggesting that the epitope recognised by this mAb is indeed unstable (data not shown).This is further supported by the observation that comparison of the N-gene nucleotide sequence of multiple isolates conforming to each reactivity pattern indicated that they were 98% identical, and predicted amino acid sequences of N protein revealed no sequence motif that could explain the specific binding of mAb 32GD12 to reactivity pattern indicated that they I isolates alone (data not shown).A phylogenetic study of canid biotype viruses from Free State province and surrounding areas, collected between 1995 and 2007, showed similar high nucleotide sequence identity (Ngoepe, Sabeta & Nel 2009).
Molecular phylogeny studies using both G-L intergenic region and full-length nucleoprotein nucleotide sequences showed that the mongoose RABV biotype from southern Africa (Zimbabwe and South Africa) could be differentiated into five main clades, which are consistent with their geographic localities (Nel et. al., 2005;Van Zyl et. al. 2010).
In contrast, in the present study, mAb typing revealed seven antigenic reactivity patterns within the mongoose RABV biotype.However, the selected isolates from each antigenic group clustered within just three of the five clades previously determined (Nel et al. 2005;Van Zyl et al. 2010).Concordant results were obtained whether the G-L intergenic region or the nucleoprotein nucleotide sequence was used for phylogenetic analysis.A high degree of diversity in mongoose rabies isolates examined with mAbs using a different mAb panel was also observed in studies in Botswana (Tremlett et al. 1994) The heterogeneity in mAb reactivity patterns observed in the present study is somewhat puzzling; however, not all of the different reaction patterns are associated with a given host species or geographic origin.In addition, only the antigenic types A, E and F are found on distinct branches of the sequence-based phylogenetic tree.Antigenic types B, C, D, and G are scattered over three of the five branches of the tree, with the three group D isolates falling into three different clades.It might well be that these are artefacts related to sample quality, as observed with the canid biotypes.These observations highlight the limitation of using the current mAb panel for higher resolution studies of RABV sub-variants.
The retrospective analysis of the present study confirms not only the presence of the two RABV biotypes but also the rare MOKV in South Africa.The true prevalence and the host species of MOKV are yet to be established.In southern Africa, MOKV has only been isolated from domestic cats, with one exception from a dog (Sabeta et al. 2007b).Since domestic cats are unknown to be a reservoir species for any lyssavirus, it is likely that they exist in close association with the true MOKV reservoir and are incidental hosts for this Lyssavirus species.The ability to identify MOKV quickly using mAbs is of great importance since the current vaccines do not protect against MOKV infections (Badrane et al. 2001;Von Teichman et al. 1998).

Conclusion
Rabies case surveillance using only pan-lyssavirus reactive antisera in dFAT does not provide sufficient information to fully understand the epidemiology of rabies in a region where multiple lyssaviruses are circulating.Whilst virus characterisation by genetic sequencing provides definitive information with regard to virus type, mAb typing is a useful screening tool that can be done quickly and at comparatively lower cost.This study has illustrated the value of mAb typing to identify species of the Lyssavirus genus and has contributed to an understanding of the epidemiology of rabies and rabies-related viruses in South Africa; however, evaluation of detailed relationships amongst isolates requires the use of more sophisticated molecular techniques that are not impacted to the same degree by sample quality.Routine lyssavirus differentiation using mAbs should be continued in synergy with nucleotide sequencing in this region.Such studies can inform the design of rabies control strategies, may assist in the identification of reservoir species (as in the case of MOKV), and may also lead to the identification of new lyssavirus species (Kuzmin et al. 2010).
Source: Authors' own creation Note: Bootstrap values are indicated at the branch nodes and branch lengths are drawn to scale.Both RABV9196FX and PV were used as out groups.

FIGURE 2 :
FIGURE 2: The neighbour-joining tree based on 1353 nucleotide sequences of nucleoprotein (a) and 592 nucleotide sequences of the G-L intergenic region (b) of the 21 mongoose rabies viruses selected from antigenic groups.

TABLE 1 :
Typical reactivity patterns of African lyssaviruses against a panel of 16 monoclonal antibodies.

TABLE 2 :
Representative mongoose rabies variants from each antigenic group that were included in phylogenetic analyses.
Note: The accession numbers of the targeted genes are also shown, whereby, G-L denotes G-L intergenic region and N denotes full nucleoprotein nucleotide sequences.

TABLE 3 :
Oligonucleotide primers used in the study showing annealing positions and nucleotide sequences.

TABLE 4 :
Indirect immunofluorescence assay reactivity patterns of canid rabies viruses obtained from South Africa.
mAb, monoclonal antibody; -, no reactivity; +, reactivity observed.†,This is the only mAB (32GD12) that differentiates between reactivity patterns I and II of canid rabies viruses.Source: Authors' own creationFIGURE 1: The geographical distribution of lyssavirus variants identified and received in South Africa between 2005 and 2009 at the Agricultural Research Council-Onderstepoort Veterinary Institute.Mokola virus Mongoose rabies biotype Canid rabies biotype

TABLE 6 :
The distribution of mongoose rabies virus variant antigenic reactivity patterns by host species.

TABLE 5 :
Indirect immunofluorescence assay reactivity patterns obtained for mongoose rabies viruses.
(Van Zyl et al. 2010)eactivity patterns observed within the mongoose RABV variant have been interpreted as supporting the hypothesis that the mongoose RABV variant has existed in this region for a longer period than the canid RABV variant and has therefore adapted well into the host species of this specific geographic locality (Nadin-Davis & Bingham 2004).Molecular clock analyses have supported this hypothesis and estimate the age of the mongoose RABV to be in the range of 200 years, consistent with anecdotal evidence from the historical record(Van Zyl et al. 2010).