Abstract
Bubonic plague is widely regarded as a disease of mainly historical importance; however, with increasing reports of incidence1,2,3 and the discovery of antibiotic-resistant strains of the plague bacterium Yersinia pestis4, it is re-emerging as a significant health concern5,6. Here we bypass the conventional human-disease models, and propose that bubonic plague is driven by the dynamics of the disease in the rat population. Using a stochastic, spatial metapopulation model, we show that bubonic plague can persist in relatively small rodent populations from which occasional human epidemics arise, without the need for external imports. This explains why historically the plague persisted despite long disease-free periods, and how the disease re-occurred in cities with tight quarantine control. In a contemporary setting, we show that human vaccination cannot eradicate the plague, and that culling of rats may prevent or exacerbate human epidemics, depending on the timing of the cull. The existence of plague reservoirs in wild rodent populations has important public-health implications for the transmission to urban rats and the subsequent risk of human outbreaks.
Similar content being viewed by others
Main
Large-scale human epidemics of bubonic plague have been recorded throughout history, from Roman times to the pandemic in the early 1900s. This disease has had a major social and demographic effect7,8,9,10; its arrival in Europe in 1348 led to the death of around one-third of the human population, and even today bubonic plague kills people in many areas of the world1,2,3. Historical data, from a variety of locations, show occasional large outbreaks of plague separated by long disease-free periods, and yet the disease clearly persists7,8,9. Understanding persistence is a common problem in general epidemic modelling11,12,13, and for bubonic plague it is a central historical question. Previous models of bubonic plague have been highly anthropocentric, modelling the disease as if it were transmitted solely within human populations10,14,15. But consideration of the biology shows that bubonic plague is primarily a disease of rodents that is spread by fleas and only occasionally infects humans; such a disease is termed a zoonosis. From this perspective, we formulate an epizootic (animal-based disease) model for the rat and flea populations, and by coupling this with a standard epidemic (human disease) model, we identify epidemic patterns and the circumstances in which the disease causes a large number of human cases.
The life cycle of the plague can be partitioned into four stages. (1) Fleas feeding on an infected rat ingest the bacteria causing bubonic plague, and soon become infectious. (2) When an infected rat dies, its fleas leave to search for a new host. (3) The fleas usually find other rats, infect them, and so spread the disease through the rodent community. (4) Only when the density of rats is low are the fleas forced to feed on alternative hosts such as humans, and a human epidemic occurs. (Although the pneumonic form of the disease can be passed directly from human to human, historical evidence suggests that this is rare8,16, and is only triggered as a consequence of many rat-related cases.) This informal set of rules can be made precise via a set of differential equations, given in the Methods. We identify the potential force of infection for humans, λH, which allows us to assess the expected number of human cases. Owing to the slow localized spread of plague in rodents8,9, the dynamics and persistence of the disease in both humans and rats is appropriately considered using a spatially structured version of this model with stochastic dynamics.
Using a metapopulation framework17,18, the community is divided into a set of subpopulations, each of which has its own stochastic dynamics19 (see Methods). Two phenomena at the local subpopulation scale control the global dynamics of the model (Fig. 1 ). When the proportion of susceptible rats is low (between 25% and 50%), an import of bubonic plague can often lead to highly persistent endemic behaviour in rats and a very low force of infection to humans. In contrast, when the proportion of susceptible rats is high (greater than 80%) the arrival of infection is followed by a short-lived rat epizootic and a large force of infection to humans (Fig. 1). The global persistence of the disease within the metapopulation depends on a small number of rat subpopulations which are in the endemic state and do not cause human cases; these endemic populations occasionally trigger large but short-lived outbreaks in neighbouring rat subpopulations which in turn lead to major human epidemics.
The potential force of infection for humans, λH, from the stochastic metapopulation model can be compared with historical data for mortality in London20 (Fig. 2a and b). Although the predicted number of human cases is erratic with long inter-epidemic periods, plague persists in the rat population (Fig. 2c). In fact, relatively small rodent metapopulations (50,000 individuals) allow the disease to persist globally for many years, although localized extinctions are common. This is attributable to the spatial heterogeneity that develops21 and the nonlinear density-dependent behaviour of the disease in rats. The persistence of the disease within the rodent population presents a new perspective on the historical data, as a sudden outbreak of bubonic plague in humans does not have to coincide with fresh imports of the disease from an external source. This resolves a controversy in the historical records for some European cities during the fifteenth and sixteenth centuries, when despite very tight quarantine controls frequent human epidemics still occurred8,22.
Contemporary control or eradication of the plague can be considered using two main approaches, vaccination or culling. However, both methods have associated complications. Because bubonic plague is a zoonosis, primarily passed from rats (via fleas) to humans, there is no ‘herd immunity’ effect23,24 associated with vaccinating humans. For most contagious diseases, once the vaccination coverage has exceeded the vaccination threshold (70–80% for smallpox, 90–95% for measles and 99% for malaria) the disease cannot persist in the population23. For bubonic plague, only those individuals vaccinated are protected, and vaccination of humans cannot provide a means of eradication because the disease is driven by the rat population. An alternative approach is to reduce the rodent population; but our model predicts that there may be difficulties with the timing of any cull. If the rats are kept at a permanently low level, then the risk of a large outbreak in rodents—and therefore the risk of human cases—is reduced. But if the cull is only brought into effect after the first human cases have been reported, then this action can create a far larger force of infection for humans (Fig. 3a). By the time human cases arise, there is already a large reservoir of infection in the rodent population; a cull therefore releases many infected fleas, which may bite humans in the absence of a suitable rodent host. This highlights the importance of considering bubonic plague as an epizootic infection and explicitly modelling the rat dynamics.
Given the difficulties in control, it is important to calculate the conditions under which a large human outbreak is likely. The recent discovery of antibiotic-resistant strains4 increases the importance of the prevention of human cases in the absence of a reliable cure. For many cities in both developing and developed countries, rats are still a serious problem25 and in many of these areas, including the United States, southern Africa, southern Asia and South America, there are large reservoirs of plague in wild rural rodents26. Hence there exists the potential for the disease to enter the highly susceptible urban rat population, with clear public-health consequences.
From the metapopulation model, the potential number of human cases associated with an outbreak is found to depend primarily on the basic reproductive ratio of the disease in rats, R0 (Fig. 3b), with human cases much more likely above some threshold value. Consideration of the disease in other wild rodent species allows us to place bounds on the flea searching efficiency a, which in turn means that the R0 threshold corresponds to a rat density of about 3,000 km-2. Many urban populations exceed this critical density, and identifying those communities at high risk is of obvious public-health importance.
Methods
Deterministic model
The natural history of the epizootic infection leads to five differential equations describing the dynamics of rats and fleas,
Here SR is the number of susceptible rats, IR is the number of infectious rats, RR is the number of resistant rats and TR = SR+ I R+ RR is the total rat population. In the absence of plague, the rat population experiences logistic growth, with slowly waning resistance to the disease at the population level governed by the parameter f. N gives the average number of fleas on a rat (the flea index), and F is the total number of infected fleas currently searching for a host. N is again modelled by logistic growth, although it is increased whenever a searching flea finds a new host. Given that fleas search randomly for a new rat host in a limited area, the probability of success is given by 1 - exp(-aTR), where a is a measure of searching efficiency.
λH is the potential force of infection to humans, and as such gives a measure of the number of expected human cases. In fact λH is calculated as the number of infected fleas that fail to find a rat host, and therefore is proportional to (but greater than) the rate at which human cases occur. This set of differential equations always shows oscillatory convergence to an equilibrium point with a very low λH.
Parameter values
See Table 1.
Stochastic metapopulation
Stochastic spatial dynamics are becoming an increasingly common form of epidemiological modelling. We move from a deterministic framework to a stochastic (or Monte Carlo) simulation, by allowing each possible event (for example, birth, death, infection) to occur randomly but at a rate given by the differential equations19. To introduce spatial structure, the entire rodent and flea population is discretized into 25 smaller subpopulations in a square lattice arrangement. Each subpopulation has its own stochastic dynamics derived from the differential equations, and adjacent subpopulations are coupled by the random movement of rats and fleas at rates µR = 0.03 and µF = 0.008, respectively. These movement rates are consistent with the observed slow spread of the disease through communities8,9.
References
Barreto, A., Aragon, M. & Epstein, R. Bubonic plague outbreak in Mozambique 1994. Lancet 345, 983–984 ( 1994).
Kumar, S. Bubonic plague in Surat. Lancet 345, 714 (1995).
Boisier, P. et al. Urban epidemic of bubonic plague in Majunga, Madagascar: epidemiological aspects. Trop. Med. Int. Health 2, 422– 427 (1997).
Galimand, M. et al. Multidrug resistance in Yersinia pestis mediated by a transferable plasmid. New Engl. J. Med. 337, 677–680 (1997).
Pinner, R. W. Addressing the challenges of emerging infectious disease. Am. J. Med. Sci. 311, 3–8 ( 1996).
McCormick, J. B. Epidemiology of emerging/re-emerging antimicrobial-resistant bacterial pathogens. Curr. Opin. Microbiol. 1, 125– 129 (1998).
Shrewsbury, J. F. D. A History of Bubonic Plague in the British Isles (Cambridge Univ. Press, Cambridge, 1970).
Appleby, A. B. The disappearance of plague: a continuing puzzle. Econ. Hist. Rev. 33, 161–173 ( 1980).
Curson, P. & McCracken, K. Plague in Sydney. The Anatomy of an Epidemic (New South Wales Univ. Press, Kensington, 1986).
Scott, S., Duncan, C. J. & Duncan, S. R. The plague in Penrith, Cumbria, 1597/8: its causes, biology and consequences. Ann. Hum. Biol. 23, 1–21 (1996).
Bartlett, M. S. Measles periodicity and community size. J. R. Statist. Soc. A 120, 48–70 (1957).
Keeling, M. J. Modelling the persistence of measles. Trends Microbiol. 5, 513–518 (1997).
Grenfell, B. T. & Dobson, A. (eds) Ecology of Infectious Diseases in Natural Populations (Cambridge Univ. Press, Cambridge, 1995).
Noble, J. V. Geographic and temporal development of plagues. Nature 250, 726–728 (1974).
Raggett, G. F. Modelling the Eyam plague. Inst. Math. Appl. 18, 221–226 (1982).
Busvine, J. R. Disease Transmission by Insects (Springer, Berlin, 1993).
Grenfell, B. & Harwood, J. (Meta)population dynamics of infectious diseases. Trends Ecol. Evol. 12, 395– 399 (1997).
Hanski, I. A. & Gilpin, M. E. (eds) Metapopulation Biology. Ecology, Genetics and Evolution (Academic, San Diego, 1997).
Renshaw, E. Modelling Biological Populations in Space and Time (Cambridge Univ. Press, Cambridge, 1991).
Graunt, J. Natural and Political Observations Mentioned In A Following Index, and Made Upon the Bills of Mortality (Roycroft, London, 1662; available at 〈http://www.ac.wwu.edu/∼stephan/Graunt/bills. html〉 (1996).
Bolker, B. M. & Grenfell, B. T. Space, persistence and dynamics of measles epidemics. Philos. Trans. R. Soc. Lond. B 348, 309–320 (1995).
Slack, P. The disappearance of plague: an alternative view. Econ. Hist. Rev. 33, 469–476 ( 1980).
Anderson, R. M. & May, R. M. Infectious Diseases of Humans (Oxford Univ. Press, London, 1992).
Fox, J. P., Elveback, L., Scott, W., Gatewood, L. & Ackerman, E. Herd-immunity—Basic concept and relevance to public-health immunization practices. Am. J. Epidemiol. 94, 179–189 (1971).
Buckle, A. P. & Smith, R. H. (eds) Rodent Pests and Their Control (Oxford Univ. Press, London, 1994).
Poland, J. D. & Barnes, A. M. in CRC Handbook Series in Zoonoses Sect. A, Vol. 1 (ed Steele, J. F.) 515– 556 (CRC, Boca Raton) (1979).
Hirst, L. F. Plague. Br. Encyclop. Med. Practice 9, 675 –698 (1938).
Macchiavello, A. Reservoirs and vectors of plague. J. Trop. Med. Hyg. 57, 1–68 (1954).
Bacot, A. W. Observations on the length of time that fleas (Ceratophyllus fasciatus ) carrying Bacillus pestis in the alimentary canals are able to survive in the absence of a host and retain the power to re-infect with the plague. J. Hyg. 6, 770– 773 (1915).
Keeling, M. J. & Gilligan, C. A. Bubonic plague: A metapopulation model of a zoonosis. Proc. R. Soc. Lond. B 267, 2219–2230 (2000).
Acknowledgements
This work was funded by the Royal Society (M.J.K. and C.A.G.), BBSRC and the Leverhulme Trust (C.A.G.), and arose from the King's College Research Centre project on spatially extended dynamics. We thank S. Betts, B. Grenfell, P. Jones, K. MacCracken, D. MacDonald, F. Nichols, J. Oeppen, R. Smith and J. Swinton for their comments on the manuscript and model parametrization.
Author information
Authors and Affiliations
Corresponding author
Rights and permissions
About this article
Cite this article
Keeling, M., Gilligan, C. Metapopulation dynamics of bubonic plague. Nature 407, 903–906 (2000). https://doi.org/10.1038/35038073
Received:
Accepted:
Issue Date:
DOI: https://doi.org/10.1038/35038073
This article is cited by
-
Population dynamics of the Multimammate rat (Mastomys natalensis) and Makundi’s brush fur rat (Lophuromys makundii) and their implications in disease persistence in Mbulu District, Tanzania
European Journal of Wildlife Research (2024)
-
Endemic infectious states below the epidemic threshold and beyond herd immunity
Communications Physics (2023)
-
The Role of Mobility in the Dynamics of the COVID-19 Epidemic in Andalusia
Bulletin of Mathematical Biology (2023)
-
Spatially structured eco-evolutionary dynamics in a host-pathogen interaction render isolated populations vulnerable to disease
Nature Communications (2022)
-
A digital reconstruction of the 1630–1631 large plague outbreak in Venice
Scientific Reports (2020)
Comments
By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.