Human ectoparasites and the spread of plague in Europe during the Second Pandemic

Katharine R. Dean; Fabienne Krauer; Lars Walløe; Ole Christian Lingjærde; Barbara Bramanti; Nils Chr. Stenseth; Boris V. Schmid

doi:10.1073/pnas.1715640115

Contributed by Nils Chr. Stenseth, December 4, 2017 (sent for review September 4, 2017; reviewed by Xavier Didelot and Kenneth L. Gage)

Significance

Plague is infamous as the cause of the Black Death (1347–1353) and later Second Pandemic (14th to 19th centuries CE), when devastating epidemics occurred throughout Europe, the Middle East, and North Africa. Despite the historical significance of the disease, the mechanisms underlying the spread of plague in Europe are poorly understood. While it is commonly assumed that rats and their fleas spread plague during the Second Pandemic, there is little historical and archaeological support for such a claim. Here, we show that human ectoparasites, like body lice and human fleas, might be more likely than rats to have caused the rapidly developing epidemics in pre-Industrial Europe. Such an alternative transmission route explains many of the notable epidemiological differences between historical and modern plague epidemics.

Abstract

Plague, caused by the bacterium Yersinia pestis, can spread through human populations by multiple transmission pathways. Today, most human plague cases are bubonic, caused by spillover of infected fleas from rodent epizootics, or pneumonic, caused by inhalation of infectious droplets. However, little is known about the historical spread of plague in Europe during the Second Pandemic (14–19th centuries), including the Black Death, which led to high mortality and recurrent epidemics for hundreds of years. Several studies have suggested that human ectoparasite vectors, such as human fleas (Pulex irritans) or body lice (Pediculus humanus humanus), caused the rapidly spreading epidemics. Here, we describe a compartmental model for plague transmission by a human ectoparasite vector. Using Bayesian inference, we found that this model fits mortality curves from nine outbreaks in Europe better than models for pneumonic or rodent transmission. Our results support that human ectoparasites were primary vectors for plague during the Second Pandemic, including the Black Death (1346–1353), ultimately challenging the assumption that plague in Europe was predominantly spread by rats.

Plague, caused by the bacterium Yersinia pestis, has been extensively studied due to its modern and historical significance. In the past, plague has famously caused at least three pandemics in human history: the First Pandemic beginning with the Justinianic Plague (6th to 8th centuries), the Second Pandemic beginning with the “Black Death” (14th to 19th centuries), and the Third Pandemic (beginning in the 19th century) (1). Today, plague persists primarily in rodent reservoirs in Asia, Africa, and the Americas, where it poses a recurrent threat to nearby human settlements (2).

The most common forms of plague infection are bubonic and pneumonic (2). Bubonic plague occurs when bacteria enter the skin, usually from the bite of an infected flea vector. The bacteria are then transported to the lymph nodes, causing characteristic swelling, or “buboes.” Bubonic plague is typically transmitted to humans from wild or commensal rodents (3), but transmission between people is also thought to occur by human ectoparasites, such as human fleas (Pulex irritans) or body lice (Pediculus humanus humanus) (4). Primary pneumonic plague occurs when aerosolized bacteria enter and infect the lungs. Pneumonic plague can also arise as a complication of bubonic or septicemic infections (2), known as secondary pneumonic plague. Individuals with pneumonic plague can transmit the disease through the respiratory route, although outbreaks of pneumonic plague are typically small because infected persons die rapidly without treatment (5). Septicemic plague occurs when bacteria infect the bloodstream, commonly from a primary pneumonic or bubonic infection (2).

A central focus of historical plague research has been to understand the spread and persistence of plague in Europe. Little is known about the transmission of plague in Europe, the Middle East, and North Africa during the Second Pandemic, including the Black Death, when the disease killed an estimated one-third of the population. Many studies (4, 6, 7) have suggested that human ectoparasites, like human fleas and body lice, were more likely than commensal rats to have caused the rapidly spreading epidemics. Proponents of the “human ectoparasite hypothesis” argue that plague epidemics during the Second Pandemic differ from the rat-associated epidemics that occurred later, during the Third Pandemic. Specifically, the geographic spread and total mortality of the Black Death far exceeds that of modern plague epidemics (8). While contemporaneous accounts of symptoms during the Second Pandemic are consistent with those of plague (7), there are no descriptions of rat epizootics, or “rat falls,” that often precede epidemics in the Third Pandemic (7 ⇓–9). Some have noted that the climate of northern Europe could not have fostered the widespread distribution of Rattus rattus (10), a claim that is supported by a scarcity of rats in the archaeological record (6). Finally, epidemiological characteristics of plague in Europe, such as a high rate of household transmission (11), are suggestive of a more direct transmission route (12).

Despite support for human ectoparasite transmission, it has been difficult to assess their historical contribution because their role in modern plague epidemics appears to be relatively minor. Today, human ectoparasite diseases have declined in most developed countries, but they are still associated with poverty and unhygienic conditions (13). In the past, human ectoparasites have been efficient vectors for diseases such as epidemic typhus (14) and relapsing fever (15). In 1941, plague-infected body lice and human fleas were found on septicemic patients during an outbreak in Morocco (16), indicating that humans can transmit the disease to lice and human fleas. In addition, recent experimental studies have demonstrated that body lice can transmit the bacteria to naive rabbits (4, 17 ⇓–19). However, the transmission from body lice and human fleas to humans has not yet been documented, and thus the importance of human ectoparasite transmission in current and historic settings remains an open question. Our theoretical analysis demonstrates that human ectoparasites may indeed play such a role.

Mathematical modeling can provide strong insight into mechanisms of plague transmission for past epidemics. Previous epidemiological models of plague during the Second Pandemic are focused mainly on modeling the spread of the disease by commensal rats during a single outbreak (20, 21). In this study, we developed a susceptible–infectious–recovered (SIR) model for plague transmission with a human ectoparasite vector and compared it to models for pneumonic and rat–flea transmission. We applied these models to nine outbreaks during the Second Pandemic, to gain a broad understanding of the transmission dynamics of plague in European epidemics. We identified the best-fitting model for each outbreak and estimated the basic reproduction number, R₀.

Methods

Historical Data.

We used data on the daily and weekly disease-induced mortality for nine plague outbreaks during the Second Pandemic (Table 1). These data were publicly available in secondary sources including published articles, books, and government reports. We digitized the epidemic data from printed tables and graphs, using the entire duration of each outbreak, apart from Eyam, which had two mortality peaks. The deterministic models we used cannot account for the stochasticity of infectious disease processes during the early phase of an epidemic; thus, for the outbreak in Eyam, we removed the first 279 data points and considered only the second, larger epidemic peak. To validate the models for pneumonic and rat-associated plague epidemics, we used three additional mortality curves from epidemics with known transmission routes during the Third Pandemic (Table S1).

View this table:

Table 1.

Summary of the Second Pandemic mortality data

Parameters.

The parameter values and initial conditions used in the models are shown in Table 2 and Table S2. Fixed values were taken from field, experimental, or epidemiological case studies when available. Unobservable parameters were estimated using Bayesian inference.

View this table:

Table 2.

Parameters for three SIR models of plague transmission

Human–Ectoparasite Model.

The transmission of bubonic plague by a human ectoparasite vector, such as human fleas or body lice, is modeled by seven differential equations:dShdt=−βlShIlNh,dIlowdt=βlShIlNh−σbIlow,dIhighdt=(1−gh)σbIlow−γbIhigh,dRhdt=ghσbIlow,dDhdt=γbIhigh,dSldt=rlSl(1−NlKl)−[(βlowIlow+βhighIhigh)SlNh],dIldt=[(βlowIlow+βhighIhigh)SlNh]−γlIl.

The five compartments for humans that are functions of time t: susceptible (Sh), infectious with mild bacteremia (Ilow), infectious with high bacteremia (Ihigh), recovered (Rh), and dead (Dh). The total living population is given by Nh=Sh+Ilow+Ihigh+Rh. The transmission of plague from vectors to humans occurs at rate βl. The model assumes that humans are mildly infectious for an average of 8 d (σb−1), and transmission is unlikely at rate βlow. Humans with mild bacteremia may recover at rate gh, which is around 40% for untreated bubonic plague. Experimental studies have shown that fleas must feed on hosts with high levels of bacteremia for reliable transmission (40). Therefore, the model assumes that moribund humans transmit plague at a high rate to vectors βhigh for an average of 2 d (γb−1). Given the short duration of the outbreaks, we did not model natural births and deaths in the human population.

Human ectoparasite vectors are modeled in two compartments (Sl, Il). The susceptible vector population grows at the intrinsic growth rate rl. The growth of the vector population is limited by the carrying capacity Kl, which is the product of the parasite index and the number of human hosts Nh. Modern studies show that the rate of body louse infestation and abundance in affected human populations ranges from 10.5 to 67.7 lice on average per person (33, 41).

There are a limited number of studies that evaluate human fleas and body lice as vectors for plague (17 ⇓–19). These studies have shown both vectors have similar transmission cycles for Y. pestis, and this makes it difficult to distinguish between the two species with either model structure or parameter values (17 ⇓–19). Our model uses parameters specific to body lice; however, the ranges for the lice and flea parameters overlap. The duration of infection γl−1 has been shown experimentally for both species, and is on average 4.5 d for human fleas and 3 d for body lice (17 ⇓–19). The model assumes that infected human fleas and body lice do not recover. The transmission of plague by human fleas is hypothesized to occur through early phase transmission, an alternative to blocked transmission observed in rat fleas (Xenopsylla cheopis) that does not require a lengthy extrinsic incubation period (42).

Pneumonic Plague Model.

The direct human-to-human transmission of plague is modeled by three differential equations:dShdt=−βpShIhNh,dIhdt=βpShIhNh−γpIh,dDhdt=γpIh.

There are three compartments for humans (Sh, Ih, Dh) and the total human population is Nh=Sh+Ih. There is no compartment for recovered individuals because the case fatality rate of untreated pneumonic plague is close to 100% (43). The human-to-human transmission of pneumonic plague occurs at rate βp. The disease-induced mortality occurs at rate γp per day and is equal to the inverse of the infectious period, which is a mean of 2.5 d for pneumonic plague (5).

Rat–Flea Model.

Based on a metapopulation model for bubonic plague by Keeling and Gilligan (35, 36), the transmission of plague in a rodent epizootic, and the spillover to humans is modeled by 10 differential equations:dSrdt=−βrSrFNr[1−e−aNr],dIrdt=βrSrFNr[1−e−aNr]−γrIr,dRrdt=grγrIr,dDrdt=(1−gr)γrIr,dHdt=rfH(1−HKf),dFdt=(1−gr)γrIrH−dfF,dShdt=−βhShFNh[e−aNr],dIhdt=βhShFNh[e−aNr]−γhIh,dRhdt=ghγhIh,dDhdt=(1−gh)γhIh.

There are four compartments for rats (Sr, Ir,Rr, Dr) and the total rat population is Nr=Sr+Ir+Rr. As epidemics within the rat population can only occur when a large proportion of the rats are susceptible to the disease, we assumed an initial black rat (Rattus rattus) population that was entirely susceptible. Although the expected ratio of urban rats to humans is about 1 rat to every 5 people (44), we allowed the prior in the model to have a maximum ratio of 1:1 rats to humans. Increasing the rat population in medieval cities allowed the simulated rat-borne plague outbreaks to more easily reach the mortality levels observed in humans during the Second Pandemic.

Rat fleas (X. cheopis) are modeled as the average number of fleas per rat, H, and the number of free infectious fleas, F. The flea population has a natural growth rate, rf, that is limited by the carrying capacity Kf. We assumed that the flea population is limited by the number of rat hosts, because X. cheopis does not reproduce on humans (45). Plague is transmitted to rats at rate βr by free infectious fleas searching for a host with searching efficiency a. We further assumed that fleas can transmit plague in the early phase (42). Rats die at a rate equal to the inverse of the infectious period γr−1, or recover with probability gr. When an infected rat dies, a number of free infectious fleas are released into the environment, depending on the average number of fleas per rat. Free infectious fleas die at rate df. The model assumes that plague is a rodent disease and that human cases are a consequence of mortality in the rat population. Therefore, susceptible humans Sh become infected by free infectious fleas at rate βh. Humans remain infected for an average of 10 d (γh−1), at which point they either recover at rate gh or die.

In the model by Keeling and Gilligan (35, 36), it is assumed that the force of infection from free infectious fleas is divided exclusively between rats and humans. However, the authors note that the true force of infection to humans is less because not every flea will find and infect a human (35). For our model, we sought to establish a range for βh that would accurately lower the force of infection to humans. To establish this range, we fitted the model to observed mortality for both rats and humans in Hong Kong in 1903 (Fig. S1) and found that the mean estimate for βh was 0.1 (Table S3). Using simulations, we found that βh should be less than 0.2 to preserve the characteristic delay and higher peak mortality of the rat epizootic compared with the human epidemic. Based on these observations, we constrained the prior for the transmission rate to humans βh to 0.0–0.2, which enabled us to use this model for outbreaks where only human mortality was available.

Bayesian Inference and Markov Chain Monte Carlo.

We fitted the deterministic models to the observed data using Bayesian inference and estimated unobservable parameters of interest. The models had a time-step of 1 d and were fitted to daily mortality or weekly mortality. Denoting the set of model parameters as Θ={S0,β,…}, the probability p of the observed data D1…m given Θ is calculated as the product of a series of Poisson random variables with mean λT equal to the human mortality in the model at times T1…m:p(D|Θ)=∏T=1me−λT(λT)DTDT!.

The parameters that we fitted were the transmission rates for each model (βlow, βhigh, βp, βr, βh) and the size of the initial primary host population that was at risk [S(0)h, S(0)r] or infected [I(0)h, I(0)r]. We assumed uniformly distributed priors and obtained posterior distributions using Markov chain Monte Carlo (MCMC) simulations with an adaptive Metropolis–Hastings algorithm implemented in PyMC2 (46) (for examples of the implementation, see https://zenodo.org/record/1043924). We ran the MCMC simulations for 180,000 iterations with a burn-in of 80,000 iterations and a thinning of 10. We assessed convergence for each model by running three independent MCMC chains and verifying that the Gelman–Rubin statistic (47) was <1.05 for each parameter. We performed model comparison using the Bayesian information criterion (BIC) from the maximum-likelihood estimates of the model parameters (48). The model with the lowest BIC value was the unique preferred model if the second-best model had a BIC value of at least 10 larger (49).

Estimation of the Basic Reproduction Number.

We estimated the basic reproduction number in each model for the primary host using the next-generation matrix method (50).

Reporting Error.

We conducted the analysis again considering different levels of underreporting (10%, 25%, and 50%) for each outbreak. To do so, we incorporated a constant probability of reporting into the likelihood function.

Results

Model Fit and Selection.

We used Bayesian MCMC and the mortality data to fit the three transmission models: human ectoparasite plague (EP), pneumonic plague (PP), and rat-borne plague (RP) (Fig. 1). The posterior means and 95% credible intervals for the fitted parameters in each model are given in Table S3. Fig. 1 shows the fit of each model to the observed mortality. For the smallest outbreaks, Eyam and Givry, it is difficult to visually distinguish between the models because the credible intervals are overlapping. In general, the human ectoparasite model fit the pattern of the observed data for the Second Pandemic outbreaks. However, the model could not account for irregularities in the observed mortality from Malta and Moscow, which have two peaks. For the pneumonic plague model, the mortality curve is right skewed compared with the observed mortality. Mortality in the rat model tended to grow slowly while plague spread through the rat population, and peaked higher than the observed mortality.

Fig. 1.

Fit of three models of plague transmission to mortality during Second Pandemic outbreaks. The observed human mortality data (black dots) and the fit (mean and 95% credible interval) of three models for plague transmission [human ectoparasite (red), pneumonic (blue), and rat–flea (green)] for nine plague outbreaks: (A) Givry, France (1348), (B) Florence, Italy (1400), (C) Barcelona, Spain (1490), (D) London, England (1563), (E) Eyam, England (1665), (F) Gdansk, Poland (1709), (G) Stockholm, Sweden (1710), (H) Moscow, Russia (1772), and (I) Island of Malta, Malta (1813).

We compared the three competing models using the BIC. Our results (Table 3) show that the human ectoparasite model had the lowest BIC value for all outbreaks, except Eyam and Givry. For the remaining outbreaks, the difference in BIC for the human ectoparasite model and the other candidate models was greater than 10, which provides strong evidence against the pneumonic and rat–flea models (50). For Eyam and Givry, the difference between the human ectoparasite model and another model was less than 10; therefore, neither model could be excluded.

View this table:

Table 3.

Comparison of transmission models and posterior estimates for the basic reproduction number for different plague models and outbreaks

To verify our model comparison method, we fitted the models to three additional Third Pandemic outbreaks with known transmission routes (Fig. S2). We found that the model with the lowest BIC matched the known modes of transmission for the outbreaks in Hong Kong (rats) and Harbin (pneumonic) (Table S5). However, we could not distinguish between two of the models for a small outbreak of rat-associated plague in Sydney, suggesting together with the results from Eyam and Givry, that our model comparison method is better suited for sufficiently large outbreaks (>750 deaths).

Basic Reproduction Number R₀.

By definition, the basic reproduction number, R₀, is the average number of secondary cases produced by a primary case, given an entirely susceptible population. In practice, R₀ is an important threshold for disease invasion. For each of the three models, we calculated R₀ from the posterior estimates of the fitted parameters (Table 3). For all of the models, R₀ was greater than 1, which is above the threshold for disease invasion. Using the human ectoparasite model, the estimated R₀ was 1.48–1.91 for all pre-Industrial outbreaks.

Reporting Error.

We considered the impact of different levels of constant underreporting of deaths throughout the epidemics on model selection (Table S6). We found that underreporting of 10% and 25% did not change the results of the model selection; under these conditions, the human ectoparasite model was the best fit for all outbreaks in Europe except Eyam and Givry. Underreporting of 50% changed the best-fitting models of Gdansk and Givry to pneumonic plague. For these cities, 50% underreporting resulted in the death of more than 90% of the population, giving preference to a pneumonic plague model where all infected individuals die from plague.

Discussion

Our study supports human ectoparasite transmission of plague during the Second Pandemic, including the Black Death. Using recent experimental data on human fleas and body lice as plague vectors, we have developed a compartmental model that captures the dynamics of human ectoparasite transmission. We have shown that, in seven out of nine localities, the human ectoparasite model was the preferred model to explain the pattern of plague mortality during an outbreak, rather than models of pneumonic and rat–flea plague transmission (Table 3). The small size of the plague outbreaks in Eyam and Givry made it difficult to distinguish transmission routes based on mortality data. For Eyam, both the human ectoparasite model and the pneumonic model produced a similar quality fit for the observed mortality. This agrees with a previous modeling study of Eyam (1665), which found that the dominant mode of transmission was an unspecified route of human-to-human transmission, rather than rodent transmission (11). Overall, our results suggest that plague transmission in European epidemics occurred predominantly through human ectoparasites, rather than commensal rat or pneumonic transmission.

The strength of our study is that we compared three plague transmission models, each representing a known or hypothetical mode of plague transmission, for nine plague outbreaks across the spatial and temporal extent of the Second Pandemic in Europe. Our study thus provides a more general understanding of plague epidemics in Europe than previous modeling studies that focus on single outbreaks, or single transmission routes (11, 20, 35, 36, 51). However, since we considered nine outbreaks over several centuries, we were limited to using simple models that could be applied systematically. Consequently, these models did not account for local conditions that can affect disease transmission, like war, famine, immunity, and public health interventions. Additionally, we did not model mixed transmission routes, and this makes it difficult to fully assess the contribution of pneumonic plague, which commonly occurs during bubonic outbreaks (52). Secondary pneumonic plague develops in an estimated 20% of bubonic cases, and this creates potential for primary pneumonic spread, even if it is not the dominant transmission route (52). Finally, we do not consider events leading up to the introduction of the disease and our results cannot be extended to plague transmission between localities, which may have involved different transmission mechanisms.

Recent studies have found human ectoparasites during plague outbreaks in the Democratic Republic of Congo (41), Tanzania (53), and Madagascar (54), but their role in these outbreaks is not clear. In the absence of modern studies on human ectoparasites as vectors for plague, our results yield inferences about the conditions necessary to produce outbreaks driven by human ectoparasite transmission. Our estimated values for R₀ using the human ectoparasite model were consistently between 1.5 and 1.9 for all nine cities. The main components of R₀ in the human ectoparasite model are the ectoparasite index and the transmission rates (βl,βlow,βhigh). From the fitted models, we obtained estimates for the transmission rates (βlow,βhigh) from humans to ectoparasites during the early and late stages of plague infection. We found that the majority of ectoparasite infections occurred during the period of high infectivity in humans, consistent with experimental evidence (40). Inferences like these not only improve our understanding of human ectoparasites as plague vectors in the past but also have important implications for limiting the size of plague outbreaks today.

Many studies have sought to clarify the mechanisms underlying the spread and maintenance of plague during the Second Pandemic. Mathematical modeling is an important tool to examine the role of different transmission mechanisms, particularly in the absence of definitive experimental, historical, and archaeological information. Here, we demonstrate that human ectoparasites appear to have been the dominant transmission mode for plague during the Second Pandemic. This alternative mode of transmission could account for many of the epidemiological differences between the Second Pandemic and those caused by rats during the Third Pandemic. Plague is undeniably a disease of significant scientific, historic, and public interest, and is still present in many parts of the world today. It is therefore crucial that we understand the full spectrum of capabilities that this versatile, pandemic disease has exhibited in the past.

Acknowledgments

We thank W. Ryan Easterday and Jukka Corander for their valuable comments. We acknowledge funding from the European Research Council under the FP7-IDEAS-ERC Program (Grant 324249). We also acknowledge funding from the Centre for Ecological and Evolutionary Synthesis.

Footnotes

↵¹To whom correspondence may be addressed. Email: k.r.dean{at}ibv.uio.no, n.c.stenseth{at}ibv.uio.no, or boris.schmid{at}gmail.com.

Author contributions: K.R.D., N.C.S., and B.V.S. designed research; K.R.D. performed research; K.R.D, F.K., and B.V.S. analyzed data; and K.R.D., F.K., L.W., O.C.L., B.B., N.C.S., and B.V.S. wrote the paper.
Reviewers: X.D., Imperial College London; and K.L.G., Centers for Disease Control and Prevention.
The authors declare no conflict of interest.
This article contains supporting information online at www.pnas.org/lookup/suppl/doi:10.1073/pnas.1715640115/-/DCSupplemental.

This open access article is distributed under Creative Commons Attribution-NonCommercial-NoDerivatives License 4.0 (CC BY-NC-ND).

References

↵
1. Bramanti B,
2. Stenseth NC,
3. Walløe L,
4. Lei X
(2016) Plague: A disease which changed the path of human civilization. Adv Exp Med Biol 918:1–26.
OpenUrl
↵
1. Dennis DT, et al.
(1999) Plague Manual: Epidemiology, Distribution, Surveillance and Control (WHO, Geneva).
↵
1. Kugeler KJ,
2. Staples JE,
3. Hinckley AF,
4. Gage KL,
5. Mead PS
(2015) Epidemiology of human plague in the United States, 1900–2012. Emerg Infect Dis 21:16–22.
OpenUrl CrossRef PubMed
↵
1. Drancourt M,
2. Houhamdi L,
3. Raoult D
(2006) Yersinia pestis as a telluric, human ectoparasite-borne organism. Lancet Infect Dis 6:234–241.
OpenUrl CrossRef PubMed
↵
1. Gani R,
2. Leach S
(2004) Epidemiologic determinants for modeling pneumonic plague outbreaks. Emerg Infect Dis 10:608–614.
OpenUrl CrossRef PubMed
↵
1. Hufthammer AK,
2. Walløe L
(2013) Rats cannot have been intermediate hosts for Yersinia pestis during medieval plague epidemics in Northern Europe. J Archaeol Sci 40:1752–1759.
OpenUrl
↵
1. Walløe L
(2008) Medieval and modern bubonic plague: Some clinical continuities. Med Hist Suppl 27:59–73.
OpenUrl
↵
1. Cohn SK
(2002) The Black Death Transformed: Disease and Culture in Early Renaissance Europe (Oxford Univ Press, London).
↵
1. Ell SR
(1979) Some evidence for interhuman transmission of medieval plague. Rev Infect Dis 1:563–566.
OpenUrl PubMed
↵
1. Davis DE
(1986) The scarcity of rats and the Black Death: An ecological history. J Interdiscip Hist 16:455–470.
OpenUrl CrossRef
↵
1. Whittles LK,
2. Didelot X
(2016) Epidemiological analysis of the Eyam plague outbreak of 1665–1666. Proc Biol Sci 283:20160618.
OpenUrl CrossRef PubMed
↵
1. The Advisory Committee appointed by the Secretary of State for India
(1907) Reports on the plague investigations in India. XXIV. General considerations regarding the spread of infection, infectivity of houses etc. in Bombay City and Island. J Hyg (Lond) 7:874–894.
OpenUrl PubMed
↵
1. Brouqui P
(2011) Arthropod-borne diseases associated with political and social disorder. Annu Rev Entomol 56:357–374.
OpenUrl CrossRef PubMed
↵
1. Smallman-Raynor M,
2. Cliff AD
(2004) War Epidemics: An Historical Geography of Infectious Diseases in Military Conflict and Civil Strife, 1850–2000 (Oxford Univ Press, Oxford).
↵
1. Bryceson AD, et al.
(1970) Louse-borne relapsing fever. Q J Med 39:129–170.
OpenUrl PubMed
↵
1. Blanc G,
2. Baltazard M
(1942) Rôle des ectoparasites humains dans la transmission de la peste. Bull Acad Natl Med 126:446.
OpenUrl
↵
1. Houhamdi L,
2. Lepidi H,
3. Drancourt M,
4. Raoult D
(2006) Experimental model to evaluate the human body louse as a vector of plague. J Infect Dis 194:1589–1596.
OpenUrl CrossRef PubMed
↵
1. Zhao WH,
2. Guo M,
3. Duan B,
4. Su LQ
(2016) Study on carrier time in Pulex irritans after infection of Yersinia pestis. China Trop Med 16:28–30.
OpenUrl
↵
1. Ayyadurai S,
2. Sebbane F,
3. Raoult D,
4. Drancourt M
(2010) Body lice, Yersinia pestis orientalis, and Black Death. Emerg Infect Dis 16:892–893.
OpenUrl CrossRef PubMed
↵
1. Massad E,
2. Coutinho FA,
3. Burattini MN,
4. Lopez LF
(2004) The Eyam plague revisited: Did the village isolation change transmission from fleas to pulmonary? Med Hypotheses 63:911–915.
OpenUrl CrossRef PubMed
↵
1. Monecke S,
2. Monecke H,
3. Monecke J
(2009) Modelling the black death. A historical case study and implications for the epidemiology of bubonic plague. Int J Med Microbiol 299:582–593.
OpenUrl PubMed
1. Biraben JN
(1975) Les Hommes et la Peste en France et dans les Pays Européens et Méditerranéens. Tome I. La Peste dans l'Histoire (Mouton, Paris). French.
1. Carmichael AG
(1986) Plague and the Poor in Renaissance Florence (Cambridge Univ Press, Cambridge, UK).
1. Smith RS
(1936) Barcelona “Bills of Mortality” and population, 1457–1590. J Polit Econ 44:84–93.
OpenUrl
1. Ferrán y Clua J,
2. Viñas y Cusí F,
3. Grau Rd
(1907) La Peste Bubónica: Memoria sobre la Epidemia Ocurrida en Porto en 1899 (Tipografia Sucesor de Sánchez, Barcelona). Spanish.
1. Creighton C
(1891) A History of Epidemics in Britain: From A.D. 664 to the Extinction of Plague (Cambridge Univ Press, Cambridge, UK).
1. Frandsen KE
(2010) The Last Plague in the Baltic Region 1709–1713 (Museum Tusculanum Press, Copenhagen).
1. Alexander JT
(2003) Bubonic Plague in Early Modern Russia: Public Health and Urban Disaster (Oxford Univ Press, Oxford).
1. Burrell WH
(1854) Appendix V. To the Second Report on Quarantine: Report of Dr. W. H. Burrell on the Plague of Malta 1813 (George E. Eyre and William Spottiswoode for Her Majesty’s Sationery Office, London).
1. The Advisory Committee appointed by the Secretary of State for India
(1908) Reports on the plague investigations in India XXVIII. Additional experiments on the septicaemia in human plague, with an account of experiments on the infectivity of the excreta. J Hyg (Lond) 8:221–235.
OpenUrl PubMed
1. Evans FC,
2. Smith FE
(1952) The intrinsic rate of natural increase for the human louse, Pediculus humanus L. Am Nat 86:299–310.
OpenUrl CrossRef
1. Peacock AD
(1916) The louse problem at the Western Front. BMJ 1:784–788.
OpenUrl FREE Full Text
↵
1. Foucault C, et al.
(2006) Oral ivermectin in the treatment of body lice. J Infect Dis 193:474–476.
OpenUrl CrossRef PubMed
1. Tollenaere C, et al.
(2010) Susceptibility to Yersinia pestis experimental infection in wild Rattus rattus, reservoir of plague in Madagascar. EcoHealth 7:242–247.
OpenUrl CrossRef PubMed
↵
1. Keeling MJ,
2. Gilligan CA
(2000) Bubonic plague: A metapopulation model of a zoonosis. Proc Biol Sci 267:2219–2230.
OpenUrl CrossRef PubMed
↵
1. Keeling MJ,
2. Gilligan CA
(2000) Metapopulation dynamics of bubonic plague. Nature 407:903–906.
OpenUrl CrossRef PubMed
1. Guernier V, et al.
(2014) Fleas of small mammals on Reunion Island: Diversity, distribution and epidemiological consequences. PLoS Negl Trop Dis 8:e3129.
OpenUrl CrossRef PubMed
1. Carrion AL
(1932) Final report on a rat-flea survey of the city of San Juan, Porto Rico. Public Health Rep 47:193–201.
OpenUrl
1. Bacot AW,
2. Martin CJ
(1924) The respective influences of temperature and moisture upon the survival of the rat flea (Xenopsylla cheopis) away from its host. J Hyg (Lond) 23:98–105.
OpenUrl PubMed
↵
1. Boegler KA,
2. Graham CB,
3. Johnson TL,
4. Montenieri JA,
5. Eisen RJ
(2016) Infection prevalence, bacterial loads, and transmission efficiency in Oropsylla montana (Siphonaptera: Ceratophyllidae) one day after exposure to varying concentrations of Yersinia pestis in blood. J Med Entomol 53:674–680.
OpenUrl CrossRef PubMed
↵
1. Piarroux R, et al.
(2013) Plague epidemics and lice, Democratic Republic of the Congo. Emerg Infect Dis 19:505–506.
OpenUrl PubMed
↵
1. Eisen RJ, et al.
(2006) Early-phase transmission of Yersinia pestis by unblocked fleas as a mechanism explaining rapidly spreading plague epizootics. Proc Natl Acad Sci USA 103:15380–15385.
OpenUrl Abstract/FREE Full Text
↵
1. Kool JL
(2005) Risk of person-to-person transmission of pneumonic plague. Clin Infect Dis 40:1166–1172.
OpenUrl CrossRef PubMed
↵
1. Davis DE,
2. Fales WT
(1949) The distribution of rats in Baltimore, Maryland. Am J Hyg 49:247–254.
OpenUrl PubMed
↵
1. Seal SC,
2. Bhatacharji LM
(1961) Epidemiological studies on plague in Calcutta. I. Bionomics of two species of ratfleas and distribution, densities and resistance of rodents in relation to the epidemiology of plague in Calcutta. Indian J Med Res 49:974–1007.
OpenUrl PubMed
↵
1. Patil A,
2. Huard D,
3. Fonnesbeck CJ
(2010) PyMC: Bayesian stochastic modelling in python. J Stat Softw 35:1–81.
OpenUrl CrossRef PubMed
↵
1. Bernardo JM, et al.
1. Gelman A,
2. Rubin DR
(1992) A single series from the Gibbs sampler provides a false sense of security. Bayesian Statistics 4, ed Bernardo JM, et al. (Oxford Univ Press, Oxford), pp 625–631.
↵
1. Schwartz G
(1978) Estimating the dimension of a model. Ann Stat 6:461–464.
OpenUrl CrossRef
↵
1. Kass RE,
2. Raftery AE
(1995) Bayes factors. J Am Stat Assoc 90:773–795.
OpenUrl CrossRef PubMed
↵
1. Diekmann O,
2. Heesterbeek JA,
3. Roberts MG
(2010) The construction of next-generation matrices for compartmental epidemic models. J R Soc Interface 7:873–885.
OpenUrl CrossRef PubMed
↵
1. Didelot X,
2. Whittles LK,
3. Hall I
(2017) Model-based analysis of an outbreak of bubonic plague in Cairo in 1801. J R Soc Interface 14:20170160.
OpenUrl CrossRef PubMed
↵
1. Alsofrom DJ,
2. Mettler FA Jr,
3. Mann JM
(1981) Radiographic manifestations of plaque in New Mexico, 1975–1980. A review of 42 proved cases. Radiology 139:561–565.
OpenUrl PubMed
↵
1. Laudisoit A, et al.
(2007) Plague and the human flea, Tanzania. Emerg Infect Dis 13:687–693.
OpenUrl CrossRef PubMed
↵
1. Ratovonjato J,
2. Rajerison M,
3. Rahelinirina S,
4. Boyer S
(2014) Yersinia pestis in Pulex irritans fleas during plague outbreak, Madagascar. Emerg Infect Dis 20:1414–1415.
OpenUrl CrossRef PubMed

Article Alerts

Email Article

Citation Tools

Request Permissions

Article Classifications

Biological Sciences
Ecology

This article has a Letter. Please see:

Relationship between Research Article and Letter - August 03, 2018

See related content:

Reply to Park et al.: Human ectoparasite transmission of plague during the Second Pandemic is still plausible
- Aug 03, 2018

Proceedings of the National Academy of Sciences: 115 (6)

Table of Contents

Submit

[1] ↵
Bramanti B,
Stenseth NC,
Walløe L,
Lei X
(2016) Plague: A disease which changed the path of human civilization. Adv Exp Med Biol 918:1–26.
OpenUrl

[2] Bramanti B,

[3] Stenseth NC,

[4] Walløe L,

[5] Lei X

[6] ↵
Dennis DT, et al.
(1999) Plague Manual: Epidemiology, Distribution, Surveillance and Control (WHO, Geneva).

[7] Dennis DT, et al.

[8] ↵
Kugeler KJ,
Staples JE,
Hinckley AF,
Gage KL,
Mead PS
(2015) Epidemiology of human plague in the United States, 1900–2012. Emerg Infect Dis 21:16–22.
OpenUrl CrossRef PubMed

[9] Kugeler KJ,

[10] Staples JE,

[11] Hinckley AF,

[12] Gage KL,

[13] Mead PS

[14] ↵
Drancourt M,
Houhamdi L,
Raoult D
(2006) Yersinia pestis as a telluric, human ectoparasite-borne organism. Lancet Infect Dis 6:234–241.
OpenUrl CrossRef PubMed

[15] Drancourt M,

[16] Houhamdi L,

[17] Raoult D

[18] ↵
Gani R,
Leach S
(2004) Epidemiologic determinants for modeling pneumonic plague outbreaks. Emerg Infect Dis 10:608–614.
OpenUrl CrossRef PubMed

[19] Gani R,

[20] Leach S

[21] ↵
Hufthammer AK,
Walløe L
(2013) Rats cannot have been intermediate hosts for Yersinia pestis during medieval plague epidemics in Northern Europe. J Archaeol Sci 40:1752–1759.
OpenUrl

[22] Hufthammer AK,

[23] Walløe L

[24] ↵
Walløe L
(2008) Medieval and modern bubonic plague: Some clinical continuities. Med Hist Suppl 27:59–73.
OpenUrl

[25] Walløe L

[26] ↵
Cohn SK
(2002) The Black Death Transformed: Disease and Culture in Early Renaissance Europe (Oxford Univ Press, London).

[27] Cohn SK

[28] ↵
Ell SR
(1979) Some evidence for interhuman transmission of medieval plague. Rev Infect Dis 1:563–566.
OpenUrl PubMed

[29] Ell SR

[30] ↵
Davis DE
(1986) The scarcity of rats and the Black Death: An ecological history. J Interdiscip Hist 16:455–470.
OpenUrl CrossRef

[31] Davis DE

[32] ↵
Whittles LK,
Didelot X
(2016) Epidemiological analysis of the Eyam plague outbreak of 1665–1666. Proc Biol Sci 283:20160618.
OpenUrl CrossRef PubMed

[33] Whittles LK,

[34] Didelot X

[35] ↵
The Advisory Committee appointed by the Secretary of State for India
(1907) Reports on the plague investigations in India. XXIV. General considerations regarding the spread of infection, infectivity of houses etc. in Bombay City and Island. J Hyg (Lond) 7:874–894.
OpenUrl PubMed

[36] The Advisory Committee appointed by the Secretary of State for India

[37] ↵
Brouqui P
(2011) Arthropod-borne diseases associated with political and social disorder. Annu Rev Entomol 56:357–374.
OpenUrl CrossRef PubMed

[38] Brouqui P

[39] ↵
Smallman-Raynor M,
Cliff AD
(2004) War Epidemics: An Historical Geography of Infectious Diseases in Military Conflict and Civil Strife, 1850–2000 (Oxford Univ Press, Oxford).

[40] Smallman-Raynor M,

[41] Cliff AD

[42] ↵
Bryceson AD, et al.
(1970) Louse-borne relapsing fever. Q J Med 39:129–170.
OpenUrl PubMed

[43] Bryceson AD, et al.

[44] ↵
Blanc G,
Baltazard M
(1942) Rôle des ectoparasites humains dans la transmission de la peste. Bull Acad Natl Med 126:446.
OpenUrl

[45] Blanc G,

[46] Baltazard M

[47] ↵
Houhamdi L,
Lepidi H,
Drancourt M,
Raoult D
(2006) Experimental model to evaluate the human body louse as a vector of plague. J Infect Dis 194:1589–1596.
OpenUrl CrossRef PubMed

[48] Houhamdi L,

[49] Lepidi H,

[50] Drancourt M,

[51] Raoult D

[52] ↵
Zhao WH,
Guo M,
Duan B,
Su LQ
(2016) Study on carrier time in Pulex irritans after infection of Yersinia pestis. China Trop Med 16:28–30.
OpenUrl

[53] Zhao WH,

[54] Guo M,

[55] Duan B,

[56] Su LQ

[57] ↵
Ayyadurai S,
Sebbane F,
Raoult D,
Drancourt M
(2010) Body lice, Yersinia pestis orientalis, and Black Death. Emerg Infect Dis 16:892–893.
OpenUrl CrossRef PubMed

[58] Ayyadurai S,

[59] Sebbane F,

[60] Raoult D,

[61] Drancourt M

[62] ↵
Massad E,
Coutinho FA,
Burattini MN,
Lopez LF
(2004) The Eyam plague revisited: Did the village isolation change transmission from fleas to pulmonary? Med Hypotheses 63:911–915.
OpenUrl CrossRef PubMed

[63] Massad E,

[64] Coutinho FA,

[65] Burattini MN,

[66] Lopez LF

[67] ↵
Monecke S,
Monecke H,
Monecke J
(2009) Modelling the black death. A historical case study and implications for the epidemiology of bubonic plague. Int J Med Microbiol 299:582–593.
OpenUrl PubMed

[68] Monecke S,

[69] Monecke H,

[70] Monecke J

[71] Biraben JN
(1975) Les Hommes et la Peste en France et dans les Pays Européens et Méditerranéens. Tome I. La Peste dans l'Histoire (Mouton, Paris). French.

[72] Biraben JN

[73] Carmichael AG
(1986) Plague and the Poor in Renaissance Florence (Cambridge Univ Press, Cambridge, UK).

[74] Carmichael AG

[75] Smith RS
(1936) Barcelona “Bills of Mortality” and population, 1457–1590. J Polit Econ 44:84–93.
OpenUrl

[76] Smith RS

[77] Ferrán y Clua J,
Viñas y Cusí F,
Grau Rd
(1907) La Peste Bubónica: Memoria sobre la Epidemia Ocurrida en Porto en 1899 (Tipografia Sucesor de Sánchez, Barcelona). Spanish.

[78] Ferrán y Clua J,

[79] Viñas y Cusí F,

[80] Grau Rd

[81] Creighton C
(1891) A History of Epidemics in Britain: From A.D. 664 to the Extinction of Plague (Cambridge Univ Press, Cambridge, UK).

[82] Creighton C

[83] Frandsen KE
(2010) The Last Plague in the Baltic Region 1709–1713 (Museum Tusculanum Press, Copenhagen).

[84] Frandsen KE

[85] Alexander JT
(2003) Bubonic Plague in Early Modern Russia: Public Health and Urban Disaster (Oxford Univ Press, Oxford).

[86] Alexander JT

[87] Burrell WH
(1854) Appendix V. To the Second Report on Quarantine: Report of Dr. W. H. Burrell on the Plague of Malta 1813 (George E. Eyre and William Spottiswoode for Her Majesty’s Sationery Office, London).

[88] Burrell WH

[89] The Advisory Committee appointed by the Secretary of State for India
(1908) Reports on the plague investigations in India XXVIII. Additional experiments on the septicaemia in human plague, with an account of experiments on the infectivity of the excreta. J Hyg (Lond) 8:221–235.
OpenUrl PubMed

[90] The Advisory Committee appointed by the Secretary of State for India

[91] Evans FC,
Smith FE
(1952) The intrinsic rate of natural increase for the human louse, Pediculus humanus L. Am Nat 86:299–310.
OpenUrl CrossRef

[92] Evans FC,

[93] Smith FE

[94] Peacock AD
(1916) The louse problem at the Western Front. BMJ 1:784–788.
OpenUrl FREE Full Text

[95] Peacock AD

[96] ↵
Foucault C, et al.
(2006) Oral ivermectin in the treatment of body lice. J Infect Dis 193:474–476.
OpenUrl CrossRef PubMed

[97] Foucault C, et al.

[98] Tollenaere C, et al.
(2010) Susceptibility to Yersinia pestis experimental infection in wild Rattus rattus, reservoir of plague in Madagascar. EcoHealth 7:242–247.
OpenUrl CrossRef PubMed

[99] Tollenaere C, et al.

[100] ↵
Keeling MJ,
Gilligan CA
(2000) Bubonic plague: A metapopulation model of a zoonosis. Proc Biol Sci 267:2219–2230.
OpenUrl CrossRef PubMed

[101] Keeling MJ,

[102] Gilligan CA

[103] ↵
Keeling MJ,
Gilligan CA
(2000) Metapopulation dynamics of bubonic plague. Nature 407:903–906.
OpenUrl CrossRef PubMed

[104] Keeling MJ,

[105] Gilligan CA

[106] Guernier V, et al.
(2014) Fleas of small mammals on Reunion Island: Diversity, distribution and epidemiological consequences. PLoS Negl Trop Dis 8:e3129.
OpenUrl CrossRef PubMed

[107] Guernier V, et al.

[108] Carrion AL
(1932) Final report on a rat-flea survey of the city of San Juan, Porto Rico. Public Health Rep 47:193–201.
OpenUrl

[109] Carrion AL

[110] Bacot AW,
Martin CJ
(1924) The respective influences of temperature and moisture upon the survival of the rat flea (Xenopsylla cheopis) away from its host. J Hyg (Lond) 23:98–105.
OpenUrl PubMed

[111] Bacot AW,

[112] Martin CJ

[113] ↵
Boegler KA,
Graham CB,
Johnson TL,
Montenieri JA,
Eisen RJ
(2016) Infection prevalence, bacterial loads, and transmission efficiency in Oropsylla montana (Siphonaptera: Ceratophyllidae) one day after exposure to varying concentrations of Yersinia pestis in blood. J Med Entomol 53:674–680.
OpenUrl CrossRef PubMed

[114] Boegler KA,

[115] Graham CB,

[116] Johnson TL,

[117] Montenieri JA,

[118] Eisen RJ

[119] ↵
Piarroux R, et al.
(2013) Plague epidemics and lice, Democratic Republic of the Congo. Emerg Infect Dis 19:505–506.
OpenUrl PubMed

[120] Piarroux R, et al.

[121] ↵
Eisen RJ, et al.
(2006) Early-phase transmission of Yersinia pestis by unblocked fleas as a mechanism explaining rapidly spreading plague epizootics. Proc Natl Acad Sci USA 103:15380–15385.
OpenUrl Abstract/FREE Full Text

[122] Eisen RJ, et al.

[123] ↵
Kool JL
(2005) Risk of person-to-person transmission of pneumonic plague. Clin Infect Dis 40:1166–1172.
OpenUrl CrossRef PubMed

[124] Kool JL

[125] ↵
Davis DE,
Fales WT
(1949) The distribution of rats in Baltimore, Maryland. Am J Hyg 49:247–254.
OpenUrl PubMed

[126] Davis DE,

[127] Fales WT

[128] ↵
Seal SC,
Bhatacharji LM
(1961) Epidemiological studies on plague in Calcutta. I. Bionomics of two species of ratfleas and distribution, densities and resistance of rodents in relation to the epidemiology of plague in Calcutta. Indian J Med Res 49:974–1007.
OpenUrl PubMed

[129] Seal SC,

[130] Bhatacharji LM

[131] ↵
Patil A,
Huard D,
Fonnesbeck CJ
(2010) PyMC: Bayesian stochastic modelling in python. J Stat Softw 35:1–81.
OpenUrl CrossRef PubMed

[132] Patil A,

[133] Huard D,

[134] Fonnesbeck CJ

[135] ↵
Bernardo JM, et al.
Gelman A,
Rubin DR
(1992) A single series from the Gibbs sampler provides a false sense of security. Bayesian Statistics 4, ed Bernardo JM, et al. (Oxford Univ Press, Oxford), pp 625–631.

[136] Bernardo JM, et al.

[137] Gelman A,

[138] Rubin DR

[139] ↵
Schwartz G
(1978) Estimating the dimension of a model. Ann Stat 6:461–464.
OpenUrl CrossRef

[140] Schwartz G

[141] ↵
Kass RE,
Raftery AE
(1995) Bayes factors. J Am Stat Assoc 90:773–795.
OpenUrl CrossRef PubMed

[142] Kass RE,

[143] Raftery AE

[144] ↵
Diekmann O,
Heesterbeek JA,
Roberts MG
(2010) The construction of next-generation matrices for compartmental epidemic models. J R Soc Interface 7:873–885.
OpenUrl CrossRef PubMed

[145] Diekmann O,

[146] Heesterbeek JA,

[147] Roberts MG

[148] ↵
Didelot X,
Whittles LK,
Hall I
(2017) Model-based analysis of an outbreak of bubonic plague in Cairo in 1801. J R Soc Interface 14:20170160.
OpenUrl CrossRef PubMed

[149] Didelot X,

[150] Whittles LK,

[151] Hall I

[152] ↵
Alsofrom DJ,
Mettler FA Jr,
Mann JM
(1981) Radiographic manifestations of plaque in New Mexico, 1975–1980. A review of 42 proved cases. Radiology 139:561–565.
OpenUrl PubMed

[153] Alsofrom DJ,

[154] Mettler FA Jr,

[155] Mann JM

[156] ↵
Laudisoit A, et al.
(2007) Plague and the human flea, Tanzania. Emerg Infect Dis 13:687–693.
OpenUrl CrossRef PubMed

[157] Laudisoit A, et al.

[158] ↵
Ratovonjato J,
Rajerison M,
Rahelinirina S,
Boyer S
(2014) Yersinia pestis in Pulex irritans fleas during plague outbreak, Madagascar. Emerg Infect Dis 20:1414–1415.
OpenUrl CrossRef PubMed

[159] Ratovonjato J,

[160] Rajerison M,

[161] Rahelinirina S,

[162] Boyer S

Main menu

User menu

Search

Human ectoparasites and the spread of plague in Europe during the Second Pandemic

Significance

Abstract

Methods

Historical Data.

Parameters.

Human–Ectoparasite Model.

Pneumonic Plague Model.

Rat–Flea Model.

Bayesian Inference and Markov Chain Monte Carlo.

Estimation of the Basic Reproduction Number.

Reporting Error.

Results

Model Fit and Selection.

Basic Reproduction Number R₀.

Reporting Error.

Discussion

Acknowledgments

Footnotes

References

Citation Manager Formats

Article Classifications

This article has a Letter. Please see:

See related content:

You May Also be Interested in

Similar Articles

Articles

PNAS Portals

Information

Main menu

User menu

Search

Significance

Abstract

Methods

Historical Data.

Parameters.

Human–Ectoparasite Model.

Pneumonic Plague Model.

Rat–Flea Model.

Bayesian Inference and Markov Chain Monte Carlo.

Estimation of the Basic Reproduction Number.

Reporting Error.

Results

Model Fit and Selection.

Basic Reproduction Number R0.

Reporting Error.

Discussion

Acknowledgments

Footnotes

References

Citation Manager Formats

Article Classifications

This article has a Letter. Please see:

See related content:

Sign up for Article Alerts

Jump to section

You May Also be Interested in

Similar Articles

Articles

PNAS Portals

Information

Basic Reproduction Number R₀.