| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
|
|
|
|||||||
|
|||||||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
1 Institut Pasteur, UPRE Lyssavirus Dynamics and Host Adaptation, National Reference Centre for Rabies, WHO Collaborating Centre for Reference and Research on Rabies, 25 rue du Docteur Roux, 75724 Paris Cedex 15, France
2 Center for Infectious Disease Dynamics, Department of Biology, The Pennsylvania State University, Mueller Laboratory, University Park, PA 16802, USA
3 Istituto Zooprofilattico Sperimentale delle Venezie, Research and Development Department, Rabies Unit, Viale dell'Università 10, Legnaro, Padova, Italy
4 Institut Pasteur de Dakar, Laboratoire Arbovirologie, 36 Avenue Pasteur, BP 220, Dakar, Senegal
5 Institut Pasteur de Bangui, Laboratoire des Arbovirus et Fièvres Hémorragiques Virales, BP 923, Bangui, Central African Republic
6 Direction des Laboratoires Vétérinaires, Niamey, Niger
7 Laboratoire Central Vétérinaire, km 8 route de Koulikoro, BP 2295, Bamako, Mali
8 Centre National d'Elevage et de Recherche Vétérinaires, Nouakchott, Mauritania
9 Virology Department, Laboratoire National d'Elevage, BP 7026, Ouagadougou 03, Burkina Faso
10 Institut Pasteur de Côte d'Ivoire, Unité des Virus du Système Nerveux, Département des Virus Epidémiques, 01 BP 490 Abidjan 01, Ivory Coast
Correspondence
Hervé Bourhy
herve.bourhy{at}pasteur.fr
 |
ABSTRACT |
|---|
 TOP ABSTRACT INTRODUCTIONMETHODSRESULTSDISCUSSIONREFERENCES |
|---|
Published online ahead of print on 9 January 2009 as DOI 10.1099/vir.0.007765-0.
The GenBank/EMBL/DDBJ accession numbers for newly acquired sequences are EU853566–EU853657 and FJ545659–FJ545685.
A supplementary table showing epidemiological information for all RABV isolates analysed in this study is available with the online version of this paper.
|
INTRODUCTION |
|---|
|
TOPABSTRACTINTRODUCTIONMETHODSRESULTSDISCUSSIONREFERENCES |
|---|
; Gould et al., 1998; Knobel et al., 2005). Dogs (Canis familiaris L.) have always been the principal host species of rabies throughout this area (Nel & Rupprecht, 2007). In most of Africa, and specifically western and central African countries, notification of rabies disease is not mandatory, so epidemiological data are scarce (Dodet et al., 2008).
The aetiological agent, rabies virus (RABV; genus Lyssavirus, family Rhabdoviridae), is an RNA virus with a single-stranded genome of negative polarity (around 12 kb in length) that contains five genes. Phylogenetic analyses have determined the existence of seven genotypes of lyssavirus, although this number is likely to increase with more intensive sampling (Bourhy et al., 1993; Delmas et al., 2008; Gould et al., 1998; Kuzmin et al., 2005). Among this biodiversity, three genotypes (1, 2 and 3) circulate in Africa, with dog-adapted viruses belonging to genotype 1 (i.e. RABV). Previous molecular epidemiological studies of rabies in Africa (Bourhy et al., 2008; David et al., 2007; Kissi et al., 1995) have revealed that four clades of genotype 1 circulate in this continent. The Africa 1 clade, adapted to dogs, is the most similar to current Eurasian RABV lineages and was therefore grouped into a larger ‘Cosmopolitan’ clade (Kissi et al., 1995; Swanepoel et al., 1993), together with a novel clade (Africa 4) identified recently in Egypt (Bourhy et al., 2008; David et al., 2007). In contrast, the Africa 3 clade is restricted to the Republic of South Africa and is adapted to mongoose, so that it constitutes an epidemiological cycle distinct from that of dog RABV (Bourhy et al., 2008; Davis et al., 2007; Kissi et al., 1995; Nel & Rupprecht, 2007). Finally, the Africa 2 clade includes RABV strains that circulate in dogs in several central and western African countries.
Although rabies has been studied extensively in wildlife populations in Europe and North America, the dynamics of RABV in domestic dog populations have been largely ignored. In particular, data concerning the diversity, distribution and origin of dog RABV circulating in western and central Africa are scarce. Prior studies were either restricted geographically or limited to a very small region of the viral genome, and compromised by a small sample size and a lack of exact spatial co-ordinates (Durr et al., 2008; Kissi et al., 1995; Sacramento et al., 1992; Smith et al., 1993). A recent study based on high-resolution temporal and spatial data demonstrated that rabies epidemics in southern and eastern Africa cycle with a periodicity of 3–6 years and also show significant synchrony across this region (Hampson et al., 2007). This suggests that considerable dispersal occurs from endemic rabies foci in Africa, perhaps due to the presence of ‘superspreader’ dogs that transmit the disease over large distances (Hampson et al., 2007).
The aim of the present study was to use a combination of phylogenetic and coalescent approaches to infer key aspects of the phylogeography, especially spatio-temporal dynamics, of dog RABV over a large geographical area covering western and central Africa. Most importantly, we addressed the question of the timescale of RABV evolution and its pattern of migration to understand how RABV spread within this geographical region. We therefore analysed a total of 182 isolates sampled from 27 African countries over 29 years; 92 of these isolates were newly collected and sequenced in the present study. To enhance the power of our phylogenetic analysis, we investigated evolutionary patterns and dynamics by using sequences of the complete nucleoprotein (N) and glycoprotein (G) genes from some samples.
|
METHODS |
|---|
|
TOPABSTRACTINTRODUCTIONMETHODSRESULTSDISCUSSIONREFERENCES |
|---|
Total RNA from the original brain samples or after the first passage in mice was extracted by using Trizol reagent (Invitrogen) according to the manufacturer's instructions. RT-PCRs and sequencing reactions were performed as described previously (Bourhy et al., 2008; Kissi et al., 1995, 1999). Additional primers used in this study were N127 (5'-ATGTAACACCTCTACAATGG-3'), N8m (5'-CAGTCTCYTCNGCCATCTC-3'), N42 (5'-CACATYTTGRTGAGTTGTCA-3'), N577 (5'-AAGATGTGYGCYAAYTGGAG-3'), GAS (5'-TCCTCAATTGAACCCCTTCT-3'), G310AS (5'-GTGTAGGGCGAAAGTGCTTT-3') and G1311AS (5'-TGCACATCAGGAAGGTGCA-3'), at positions 55–74, 1568–1586, 633–651, 644–663, 5850–5869, 3624–3643 and 4625–4643, respectively, of the lyssavirus genome (Delmas et al., 2008).
Evolutionary analysis.
The gene sequences newly described here were combined with relevant sequences from GenBank, resulting in three datasets of (i) 97 N gene sequences from the Africa 1 and Africa 2 lineages combined, (ii) 134 N gene sequences from the Africa 2 lineage, and (iii) 34 G gene sequences from the Africa 2 lineage. All alignments are available from the authors on request. For each dataset, we inferred the maximum clade credibility (MCC) phylogenetic tree by using the Bayesian Markov chain Monte Carlo (MCMC) method available in the BEAST package (Drummond & Rambaut, 2007), thereby incorporating information on sampling time. Posterior probability values provide an assessment of the degree of support for each node on the tree. This analysis utilized a relaxed (uncorrelated log-normal) molecular clock and the HKY85+
4 model of nucleotide substitution, although highly similar results were obtained under more complex substitution models. As demographic history can be considered a nuisance parameter in our study, we utilized the Bayesian skyline model as a coalescent prior. All chains were run for a sufficient length to ensure convergence, with 10 % removed as burn-in. As well as generating the MCMC tree, this analysis also allowed us to estimate both the rate of nucleotide substitution per site (substitutions per site year–1) and the time to most recent common ancestor (TMRCA) in years. The degree of uncertainty in each parameter estimate is provided by 95 % highest posterior density (HPD) values.
We used a parsimony-based approach to determine the geographical structure of dog RABV of Africa 2, based on the MCC tree of the N gene and with the country from where each sequence was collected as the unit of analysis (although this broad-scale categorization clearly hides a great deal of geographical structure). Further, because the MCC tree is automatically rooted through the assumption of a (relaxed) molecular clock, we were also able to infer the direction of migration events. Overall, we collected sequences from at least 30 different cities, encompassing 14 countries, in western and central Africa. Each RABV sequence was first assigned a character state reflecting its country of origin. The number of unambiguous changes in character state observed among each country, as determined from the MCMC tree, was then recorded and compared with the number expected under the null hypothesis of entirely random mixing by repeating this analysis on 1000 randomized trees. A matrix of the observed-minus-expected character-state changes was then constructed to determine the strength of migration (positive values) or population substitution (negative values). All of these analyses were conducted by using the PAUP* package (Swofford, 2003).
Finally, to determine the nature of the selection pressures acting on the Africa 2 RABV sequences, we estimated the mean ratio of non-synonymous (dN) to synonymous (dS) substitution per site (ratio dN/dS) by using the maximum-likelihood single likelihood ancestor counting (SLAC) method available through the DATAMONKEY web interface of the HyPhy package (http://www.datamonkey.org).
|
RESULTS |
|---|
|
TOPABSTRACTINTRODUCTIONMETHODSRESULTSDISCUSSIONREFERENCES |
|---|
. The topology is similar to those of earlier phylogenetic analyses of the N gene (Durr et al., 2008; Kissi et al., 1995). In particular, the sequences were clearly divided into two major clusters, identified as the ‘Cosmopolitan’ and ‘Africa 2’ clades. The Cosmopolitan clade includes isolates that are distributed in northern, central, eastern and southern Africa, as shown previously (Bourhy et al., 2008; Kissi et al., 1995; Nel & Rupprecht, 2007). More notably, we provide the first evidence for a widely distributed Africa 2 clade, comprising dog isolates that have a wide geographical range in western and central Africa, including Guinea, Sierra Leone, Senegal, Niger, Nigeria, Mauritania, Ivory Coast, Burkina Faso, Cameroon, Benin, Chad, Mali, Gambia and the Central African Republic (CAR) (Fig. 2).
|
|
The Africa 2 sequences could also be grouped arbitrarily on the basis of phylogenetic placement into eight main groups (denoted A–H), all of which were supported by strong Bayesian posterior probability values (Figs 2
and 3). The only exceptions were groups A and D, which include only one isolate from Chad (Chad 73/2006) and one from Benin (8697BEN), respectively. Group B consists of isolates originating in Guinea, Sierra Leone, Gambia and the Casamance region of southern Senegal (isolate Nb 07023SEN). Group C consists of one isolate originating in Chad (isolate 9218TCH) and another from the north-western part of CAR (Bossoangoa) (isolate Nb 07149RCA). Group E contains viruses sampled from Cameroon, Chad, Nigeria and Niger (9014NIG). Group F includes isolates originating in Niger and Burkina Faso, whilst group G consists of viruses from Mauritania, Burkina Faso and Senegal. Finally, group H contains viruses from Mauritania, Mali, Burkina Faso and the Ivory Coast.
|
), with equivalent phylogenetic grouping (B–H; no group A sequences were available for analysis) and with matching geographical distributions.
|
; Davis et al., 2007), as well as a more global analysis of dog RABV sampled worldwide (Bourhy et al., 2008). By using the same approach, we were able to estimate the TMRCA of the Africa 2 clade to be 163 years (95 % HPD, 72–288 years), which corresponds to the year 1845 (1720–1936). These data also depict a gradual diversification of the Africa 2 clade into the eight phylogenetic groups that may have occupied this entire geographical range by the year 1955 (1935–1972), perhaps 100 years after the TMRCA of the clade as a whole (Fig. 3). A similar timescale, with overlapping 95 % HPD values, was observed for the smaller sample of G gene sequences (mean, 162 years; 95 % HPD, 112–225 years), suggesting that these estimates of the age of genetic diversity are robust.
Finally, the overall (mean) dN/dS values of the N and G genes of Africa 2 RABV were very low (0.070 and 0.080, respectively); this is indicative of strong selective constraints, as documented previously in RABV (Bourhy et al., 2008; Holmes et al., 2002).
Population and spatial dynamics of dog RABV
We also conducted a more detailed investigation of the patterns and dynamics of the spatial diffusion of the Africa 2 clade, using the larger N gene dataset. This revealed a very strong population subdivision (P
0.001), similar to that observed for RABV in other geographical areas (Bourhy et al., 2008). Indeed, of 182 pairwise comparisons, only 11 (6 %) exhibited a positive correlation between countries (>0); this is indicative of migration between them (Table 1). The strongest evidence for inter-country migration is for from Chad to Nigeria (1.06), Chad to Benin (0.903), Chad to CAR (0.929), Chad to Cameroon (0.777), Chad to Niger (0.333), Niger to Burkina Faso (0.813) and Burkina Faso to Mali (0.636); this is indicative of a general westward movement across west/central Africa. However, in all cases, the strength of migration was weak, with population subdivision being by far the strongest signal in these data.
|
|
DISCUSSION |
|---|
|
TOPABSTRACTINTRODUCTIONMETHODSRESULTSDISCUSSIONREFERENCES |
|---|
; Fevre et al., 2006; Windiyaningsih et al., 2004). Ancient texts describe the existence of this disease in Mesopotamia. From this location, the virus is proposed to have spread to Europe and then, following patterns of human colonization, to Africa (Nel & Rupprecht, 2007). More recently, we were able to determine a common origin of all dog viruses circulating globally and to propose that the ancestor of these viruses existed <1500 years ago, perhaps in the Indian subcontinent (Bourhy et al., 2008).
The results described here provide the first evidence for the spread of dog RABV in western and central Africa, associated with the emergence of an Africa 2 clade, which is the dominant lineage of RABV in this area. Notably, this clade is surrounded to the north, the east and the south by the Cosmopolitan clade (Fig. 1). Such a spatial and phylogenetic distinction strongly supports the idea that these clades represent independent introductions into Africa.
More notably, the recent TMRCA for this Africa 2 clade, with a mean date (from the N gene) of 1845 (95 % HPD, 1720–1936), is consistent with the documented timescale of the expanding European colonial influence in western and central Africa. At the beginning of the 19th century, the French had contact in different areas of coastal west Africa, but their efforts were most clearly focused on the Senegal River area and its hinterland (Crowder, 1990). From 1880, French control was established over much of northern, western and central Africa, such that, by the early years of the 20th century, the French held most of what would become their colonial territory in western and central Africa (including present-day Mauritania, Senegal, Mali, Burkina Faso, Benin, Guinea, Ivory Coast, Chad, CAR, Niger and the Democratic Republic of Congo).
Our analysis of phylogeographical structure also provides information on the pattern of RABV spread in this part of Africa. In particular, the facts that (i) sequences sampled from Chad fall into multiple clades, notably those that tend to be basal in the tree, and (ii) those sequences sampled from countries bordering the Atlantic Ocean tend to fall at distant locations, are suggestive of an initial introduction of RABV into the eastern part of west/central Africa, with Chad as a potential source population, followed by a general westward (and southward) diffusion. However, the pattern is not strong and will need to be confirmed with a larger sample of sequences and more refined spatial data. Indeed, the main spatial pattern present in our data is that of population subdivision, with only limited evidence of viral movement among countries. As such, the Africa 2 clade satisfies the model proposed for dog RABV in general; that is, of a series of spatially distinct clusters that experience relatively little contact among them (Bourhy et al., 2008).
Overall, our study illustrates how the establishment and intensification of travel and trade routes between African countries following colonization and during the first half of the 20th century have been accompanied by the spread of rabies in dogs among a large part of west/central Africa. It is also possible that parallel urbanization facilitated the spread and maintenance of dog RABV in the same region. In addition, our detailed phylogeographical analysis reveals that the Africa 2 clade spread gradually – perhaps taking longer than a century – to western and central Africa, with no evidence for positive selection acting on either the N or the G gene.
Our data sit in marked contrast to a recent study based on high-resolution temporal and spatial data, which suggested that RABV spreads rapidly and continually from endemic rabies regions in Africa (Hampson et al., 2007). Specifically, the authors explained the occurrence of rabies epidemics in southern and eastern Africa by suggesting that some ‘superspreader’ dogs may transmit the disease over large distances (Hampson et al., 2007). However, the very strong geographical clustering of RABV sequences depicted in our study, as well as a timescale of spatial diffusion measured in decades, argue against any model in which superspreaders transmit their virus rapidly over a large territory. Similarly, the possible transportation of latent, or infectious, dogs by people travelling over this vast region does not seem to have a large epidemiological impact, as demonstrated by the limited movement of virus among localities.
More detailed phylogenetic and epidemiological analysis of higher-resolution data from across temporal and spatial scales is clearly necessary to determine the mechanisms underlying circulation patterns of canine rabies and to develop a more predictive understanding of the spatio-temporal dynamics of rabies in Africa and other localities. As such, this study also highlights the need for improved inter-country collaboration to better describe the spread of dog RABV among countries in western and central Africa and to obtain a reliable picture of rabies epidemiology in this region (Dodet et al., 2008). Finally, our study is also of importance when trying to design an effective strategy for the control and elimination of dog rabies in western and central Africa. In particular, we demonstrate, for the first time, that there is no spread of rabies between the countries of northern Africa and those of the sub-Saharan region and, further, that the exchange of viruses between countries of western and central Africa is limited. As a consequence, a progressive strategy of rabies elimination from western and central Africa is conceivable.
|
ACKNOWLEDGEMENTS |
|---|
. |
REFERENCES |
|---|
|
TOPABSTRACTINTRODUCTIONMETHODSRESULTSDISCUSSIONREFERENCES |
|---|
Bourhy, H., Kissi, B., Audry, L., Smreczak, M., Sadkowska-Todys, M., Kulonen, K., Tordo, N., Zmudzinski, J. F. & Holmes, E. C. (1999). Ecology and evolution of rabies virus in Europe. J Gen Virol 80, 2545–2557.
Bourhy, H., Cowley, J. A., Larrous, F., Holmes, E. C. & Walker, P. J. (2005). Phylogenetic relationships among rhabdoviruses inferred using the L polymerase gene. J Gen Virol 86, 2849–2858.
Bourhy, H., Reynes, J. M., Dunham, E. J., Dacheux, L., Larrous, F., Huong, V. T., Xu, G., Yan, J., Miranda, M. E. & Holmes, E. C. (2008). The origin and phylogeography of dog rabies virus. J Gen Virol 89, 2673–2681.
Crowder, M. (1990). History of French West Africa until independence. In Africa South of the Sahara 1990. Edited by O. Crowder, R. Crowder & M. Crowder. London: Europa Publications.
David, D., Hughes, G. J., Yakobson, B. A., Davidson, I., Un, H., Aylan, O., Kuzmin, I. V. & Rupprecht, C. E. (2007). Identification of novel canine rabies virus clades in the Middle East and North Africa. J Gen Virol 88, 967–980.
Davis, P. L., Rambaut, A., Bourhy, H. & Holmes, E. C. (2007). The evolutionary dynamics of canid and mongoose rabies virus in southern Africa. Arch Virol 152, 1251–1258.[CrossRef][Medline]
Delmas, O., Holmes, E. C., Talbi, C., Larrous, F., Dacheux, L., Bouchier, C. & Bourhy, H. (2008). Genomic diversity and evolution of the lyssaviruses. PLoS One 3, e2057[CrossRef][Medline]
Dodet, B., Adjogoua, E. V., Aguemon, A. R., Amadou, O. H., Atipo, A. L., Baba, B. A., Bara Ada, S., Boumandouki, P., Bourhy, H. & other authors (2008). Fighting rabies in Africa: the Africa Rabies Expert Bureau (AfroREB). Vaccine 26, 6295–6298.[CrossRef][Medline]
Drummond, A. J. & Rambaut, A. (2007). BEAST: Bayesian evolutionary analysis by sampling trees. BMC Evol Biol 7, 214[CrossRef][Medline]
Durr, S., Naissengar, S., Mindekem, R., Diguimbye, C., Niezgoda, M., Kuzmin, I., Rupprecht, C. E. & Zinsstag, J. (2008). Rabies diagnosis for developing countries. PLoS Negl Trop Dis 2, e206[CrossRef][Medline]
Fevre, E. M., Bronsvoort, B. M., Hamilton, K. A. & Cleaveland, S. (2006). Animal movements and the spread of infectious diseases. Trends Microbiol 14, 125–131.[CrossRef][Medline]
Gould, A. R., Hyatt, A. D., Lunt, R., Kattenbelt, J. A., Hengstberger, S. & Blacksell, S. D. (1998). Characterisation of a novel lyssavirus isolated from Pteropid bats in Australia. Virus Res 54, 165–187.[CrossRef][Medline]
Hampson, K., Dushoff, J., Bingham, J., Bruckner, G., Ali, Y. H. & Dobson, A. (2007). Synchronous cycles of domestic dog rabies in sub-Saharan Africa and the impact of control efforts. Proc Natl Acad Sci U S A 104, 7717–7722.
Holmes, E. C., Woelk, C. H., Kassis, R. & Bourhy, H. (2002). Genetic constraints and the adaptive evolution of rabies virus in nature. Virology 292, 247–257.[CrossRef][Medline]
Kissi, B., Tordo, N. & Bourhy, H. (1995). Genetic polymorphism in the rabies virus nucleoprotein gene. Virology 209, 526–537.[CrossRef][Medline]
Kissi, B., Badrane, H., Audry, L., Lavenu, A., Tordo, N., Brahimi, M. & Bourhy, H. (1999). Dynamics of rabies virus quasispecies during serial passages in heterologous hosts. J Gen Virol 80, 2041–2050.
Knobel, D. L., Cleaveland, S., Coleman, P. G., Fevre, E. M., Meltzer, M. I., Miranda, M. E., Shaw, A., Zinsstag, J. & Meslin, F. X. (2005). Re-evaluating the burden of rabies in Africa and Asia. Bull World Health Organ 83, 360–368.[Medline]
Kuzmin, I. V., Hughes, G. J., Botvinkin, A. D., Orciari, L. A. & Rupprecht, C. E. (2005). Phylogenetic relationships of Irkut and West Caucasian bat viruses within the Lyssavirus genus and suggested quantitative criteria based on the N gene sequence for lyssavirus genotype definition. Virus Res 111, 28–43.[CrossRef][Medline]
Nel, L. H. & Rupprecht, C. E. (2007). Emergence of lyssaviruses in the Old World: the case of Africa. Curr Top Microbiol Immunol 315, 161–193.[Medline]
Sacramento, D., Badrane, H., Bourhy, H. & Tordo, N. (1992). Molecular epidemiology of rabies virus in France: comparison with vaccine strains. J Gen Virol 73, 1149–1158.
Smith, J. S., Yager, P. A. & Orciari, L. A. (1993). Rabies in wild and domestic carnivores of Africa: epidemiological and historical associations determined by limited sequence analysis. Onderstepoort J Vet Res 60, 307–314.[Medline]
Swanepoel, R., Barnard, B. J., Meredith, C. D., Bishop, G. C., Bruckner, G. K., Foggin, C. M. & Hubschle, O. J. (1993). Rabies in southern Africa. Onderstepoort J Vet Res 60, 325–346.[Medline]
Swofford, D. L. (2003). PAUP*: Phylogenetic analysis using parsimony (and other methods), version 4.0b10. Sunderland, MA: Sinauer Associates.
Windiyaningsih, C., Wilde, H., Meslin, F. X., Suroso, T. & Widarso, H. S. (2004). The rabies epidemic on Flores Island, Indonesia (1998–2003). J Med Assoc Thai 87, 1389–1393.[Medline]
Received 7 October 2008;
accepted 21 December 2008.
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
| INT J SYST EVOL MICROBIOL | MICROBIOLOGY | J GEN VIROL |
| J MED MICROBIOL | ALL SGM JOURNALS | |
