Skip to main content

Main menu

  • Home
  • Articles
    • Current
    • Special Feature Articles - Most Recent
    • Special Features
    • Colloquia
    • Collected Articles
    • PNAS Classics
    • List of Issues
  • Front Matter
    • Front Matter Portal
    • Journal Club
  • News
    • For the Press
    • This Week In PNAS
    • PNAS in the News
  • Podcasts
  • Authors
    • Information for Authors
    • Editorial and Journal Policies
    • Submission Procedures
    • Fees and Licenses
  • Submit
  • Submit
  • About
    • Editorial Board
    • PNAS Staff
    • FAQ
    • Accessibility Statement
    • Rights and Permissions
    • Site Map
  • Contact
  • Journal Club
  • Subscribe
    • Subscription Rates
    • Subscriptions FAQ
    • Open Access
    • Recommend PNAS to Your Librarian

User menu

  • Log in
  • My Cart

Search

  • Advanced search
Home
Home
  • Log in
  • My Cart

Advanced Search

  • Home
  • Articles
    • Current
    • Special Feature Articles - Most Recent
    • Special Features
    • Colloquia
    • Collected Articles
    • PNAS Classics
    • List of Issues
  • Front Matter
    • Front Matter Portal
    • Journal Club
  • News
    • For the Press
    • This Week In PNAS
    • PNAS in the News
  • Podcasts
  • Authors
    • Information for Authors
    • Editorial and Journal Policies
    • Submission Procedures
    • Fees and Licenses
  • Submit
Research Article

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

View ORCID ProfileKatharine R. Dean, Fabienne Krauer, Lars Walløe, View ORCID ProfileOle Christian Lingjærde, Barbara Bramanti, Nils Chr. Stenseth, and Boris V. Schmid
  1. aCentre for Ecological and Evolutionary Synthesis (CEES), Department of Biosciences, University of Oslo, N-0316 Oslo, Norway;
  2. bDepartment of Physiology, Institute of Basic Medical Sciences, University of Oslo, N-0317 Oslo, Norway;
  3. cDepartment of Computer Science, University of Oslo, N-0316 Oslo, Norway;
  4. dDepartment of Biomedical and Specialty Surgical Sciences, Faculty of Medicine, Pharmacy and Prevention, University of Ferrara, 35-441221 Ferrara, Italy

See allHide authors and affiliations

PNAS February 6, 2018 115 (6) 1304-1309; first published January 16, 2018; https://doi.org/10.1073/pnas.1715640115
Katharine R. Dean
aCentre for Ecological and Evolutionary Synthesis (CEES), Department of Biosciences, University of Oslo, N-0316 Oslo, Norway;
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • ORCID record for Katharine R. Dean
  • For correspondence: k.r.dean@ibv.uio.no n.c.stenseth@ibv.uio.no boris.schmid@gmail.com
Fabienne Krauer
aCentre for Ecological and Evolutionary Synthesis (CEES), Department of Biosciences, University of Oslo, N-0316 Oslo, Norway;
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Lars Walløe
bDepartment of Physiology, Institute of Basic Medical Sciences, University of Oslo, N-0317 Oslo, Norway;
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Ole Christian Lingjærde
cDepartment of Computer Science, University of Oslo, N-0316 Oslo, Norway;
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • ORCID record for Ole Christian Lingjærde
Barbara Bramanti
aCentre for Ecological and Evolutionary Synthesis (CEES), Department of Biosciences, University of Oslo, N-0316 Oslo, Norway;
dDepartment of Biomedical and Specialty Surgical Sciences, Faculty of Medicine, Pharmacy and Prevention, University of Ferrara, 35-441221 Ferrara, Italy
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Nils Chr. Stenseth
aCentre for Ecological and Evolutionary Synthesis (CEES), Department of Biosciences, University of Oslo, N-0316 Oslo, Norway;
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • For correspondence: k.r.dean@ibv.uio.no n.c.stenseth@ibv.uio.no boris.schmid@gmail.com
Boris V. Schmid
aCentre for Ecological and Evolutionary Synthesis (CEES), Department of Biosciences, University of Oslo, N-0316 Oslo, Norway;
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • For correspondence: k.r.dean@ibv.uio.no n.c.stenseth@ibv.uio.no boris.schmid@gmail.com
  1. Contributed by Nils Chr. Stenseth, December 4, 2017 (sent for review September 4, 2017; reviewed by Xavier Didelot and Kenneth L. Gage)

  • Article
  • Figures & SI
  • Info & Metrics
  • PDF
Loading

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.

  • Yersinia pestis
  • Black Death
  • SIR modeling
  • Bayesian analysis
  • Monte Carlo Markov chain

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, R0.

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:
  • View inline
  • View popup
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:
  • View inline
  • View popup
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.
  • Download figure
  • Open in new tab
  • Download powerpoint
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:
  • View inline
  • View popup
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 R0.

By definition, the basic reproduction number, R0, is the average number of secondary cases produced by a primary case, given an entirely susceptible population. In practice, R0 is an important threshold for disease invasion. For each of the three models, we calculated R0 from the posterior estimates of the fitted parameters (Table 3). For all of the models, R0 was greater than 1, which is above the threshold for disease invasion. Using the human ectoparasite model, the estimated R0 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 R0 using the human ectoparasite model were consistently between 1.5 and 1.9 for all nine cities. The main components of R0 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

  • ↵1To 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.

  • Copyright © 2018 the Author(s). Published by PNAS.

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

References

  1. ↵
    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
  2. ↵
    1. Dennis DT, et al.
    (1999) Plague Manual: Epidemiology, Distribution, Surveillance and Control (WHO, Geneva).
  3. ↵
    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.
    OpenUrlCrossRefPubMed
  4. ↵
    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.
    OpenUrlCrossRefPubMed
  5. ↵
    1. Gani R,
    2. Leach S
    (2004) Epidemiologic determinants for modeling pneumonic plague outbreaks. Emerg Infect Dis 10:608–614.
    OpenUrlCrossRefPubMed
  6. ↵
    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
  7. ↵
    1. Walløe L
    (2008) Medieval and modern bubonic plague: Some clinical continuities. Med Hist Suppl 27:59–73.
    OpenUrl
  8. ↵
    1. Cohn SK
    (2002) The Black Death Transformed: Disease and Culture in Early Renaissance Europe (Oxford Univ Press, London).
  9. ↵
    1. Ell SR
    (1979) Some evidence for interhuman transmission of medieval plague. Rev Infect Dis 1:563–566.
    OpenUrlPubMed
  10. ↵
    1. Davis DE
    (1986) The scarcity of rats and the Black Death: An ecological history. J Interdiscip Hist 16:455–470.
    OpenUrlCrossRef
  11. ↵
    1. Whittles LK,
    2. Didelot X
    (2016) Epidemiological analysis of the Eyam plague outbreak of 1665–1666. Proc Biol Sci 283:20160618.
    OpenUrlCrossRefPubMed
  12. ↵
    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.
    OpenUrlPubMed
  13. ↵
    1. Brouqui P
    (2011) Arthropod-borne diseases associated with political and social disorder. Annu Rev Entomol 56:357–374.
    OpenUrlCrossRefPubMed
  14. ↵
    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).
  15. ↵
    1. Bryceson AD, et al.
    (1970) Louse-borne relapsing fever. Q J Med 39:129–170.
    OpenUrlPubMed
  16. ↵
    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
  17. ↵
    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.
    OpenUrlCrossRefPubMed
  18. ↵
    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
  19. ↵
    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.
    OpenUrlCrossRefPubMed
  20. ↵
    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.
    OpenUrlCrossRefPubMed
  21. ↵
    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.
    OpenUrlPubMed
    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.
    OpenUrlPubMed
    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.
    OpenUrlCrossRef
    1. Peacock AD
    (1916) The louse problem at the Western Front. BMJ 1:784–788.
    OpenUrlFREE Full Text
  22. ↵
    1. Foucault C, et al.
    (2006) Oral ivermectin in the treatment of body lice. J Infect Dis 193:474–476.
    OpenUrlCrossRefPubMed
    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.
    OpenUrlCrossRefPubMed
  23. ↵
    1. Keeling MJ,
    2. Gilligan CA
    (2000) Bubonic plague: A metapopulation model of a zoonosis. Proc Biol Sci 267:2219–2230.
    OpenUrlCrossRefPubMed
  24. ↵
    1. Keeling MJ,
    2. Gilligan CA
    (2000) Metapopulation dynamics of bubonic plague. Nature 407:903–906.
    OpenUrlCrossRefPubMed
    1. Guernier V, et al.
    (2014) Fleas of small mammals on Reunion Island: Diversity, distribution and epidemiological consequences. PLoS Negl Trop Dis 8:e3129.
    OpenUrlCrossRefPubMed
    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.
    OpenUrlPubMed
  25. ↵
    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.
    OpenUrlCrossRefPubMed
  26. ↵
    1. Piarroux R, et al.
    (2013) Plague epidemics and lice, Democratic Republic of the Congo. Emerg Infect Dis 19:505–506.
    OpenUrlPubMed
  27. ↵
    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.
    OpenUrlAbstract/FREE Full Text
  28. ↵
    1. Kool JL
    (2005) Risk of person-to-person transmission of pneumonic plague. Clin Infect Dis 40:1166–1172.
    OpenUrlCrossRefPubMed
  29. ↵
    1. Davis DE,
    2. Fales WT
    (1949) The distribution of rats in Baltimore, Maryland. Am J Hyg 49:247–254.
    OpenUrlPubMed
  30. ↵
    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.
    OpenUrlPubMed
  31. ↵
    1. Patil A,
    2. Huard D,
    3. Fonnesbeck CJ
    (2010) PyMC: Bayesian stochastic modelling in python. J Stat Softw 35:1–81.
    OpenUrlCrossRefPubMed
  32. ↵
    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.
  33. ↵
    1. Schwartz G
    (1978) Estimating the dimension of a model. Ann Stat 6:461–464.
    OpenUrlCrossRef
  34. ↵
    1. Kass RE,
    2. Raftery AE
    (1995) Bayes factors. J Am Stat Assoc 90:773–795.
    OpenUrlCrossRefPubMed
  35. ↵
    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.
    OpenUrlCrossRefPubMed
  36. ↵
    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.
    OpenUrlCrossRefPubMed
  37. ↵
    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.
    OpenUrlPubMed
  38. ↵
    1. Laudisoit A, et al.
    (2007) Plague and the human flea, Tanzania. Emerg Infect Dis 13:687–693.
    OpenUrlCrossRefPubMed
  39. ↵
    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.
    OpenUrlCrossRefPubMed
PreviousNext
Back to top
Article Alerts
Email Article

Thank you for your interest in spreading the word on PNAS.

NOTE: We only request your email address so that the person you are recommending the page to knows that you wanted them to see it, and that it is not junk mail. We do not capture any email address.

Enter multiple addresses on separate lines or separate them with commas.
Human ectoparasites and the spread of plague in Europe during the Second Pandemic
(Your Name) has sent you a message from PNAS
(Your Name) thought you would like to see the PNAS web site.
CAPTCHA
This question is for testing whether or not you are a human visitor and to prevent automated spam submissions.
Citation Tools
Human ectoparasites and spread of plague in Europe
Katharine R. Dean, Fabienne Krauer, Lars Walløe, Ole Christian Lingjærde, Barbara Bramanti, Nils Chr. Stenseth, Boris V. Schmid
Proceedings of the National Academy of Sciences Feb 2018, 115 (6) 1304-1309; DOI: 10.1073/pnas.1715640115

Citation Manager Formats

  • BibTeX
  • Bookends
  • EasyBib
  • EndNote (tagged)
  • EndNote 8 (xml)
  • Medlars
  • Mendeley
  • Papers
  • RefWorks Tagged
  • Ref Manager
  • RIS
  • Zotero
Request Permissions
Share
Human ectoparasites and spread of plague in Europe
Katharine R. Dean, Fabienne Krauer, Lars Walløe, Ole Christian Lingjærde, Barbara Bramanti, Nils Chr. Stenseth, Boris V. Schmid
Proceedings of the National Academy of Sciences Feb 2018, 115 (6) 1304-1309; DOI: 10.1073/pnas.1715640115
del.icio.us logo Digg logo Reddit logo Twitter logo CiteULike logo Facebook logo Google logo Mendeley logo
  • Tweet Widget
  • Facebook Like
  • Mendeley logo Mendeley

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

Sign up for Article Alerts

Jump to section

  • Article
    • Abstract
    • Methods
    • Results
    • Discussion
    • Acknowledgments
    • Footnotes
    • References
  • Figures & SI
  • Info & Metrics
  • PDF

You May Also be Interested in

Setting sun over a sun-baked dirt landscape
Core Concept: Popular integrated assessment climate policy models have key caveats
Better explicating the strengths and shortcomings of these models will help refine projections and improve transparency in the years ahead.
Image credit: Witsawat.S.
Model of the Amazon forest
News Feature: A sea in the Amazon
Did the Caribbean sweep into the western Amazon millions of years ago, shaping the region’s rich biodiversity?
Image credit: Tacio Cordeiro Bicudo (University of São Paulo, São Paulo, Brazil), Victor Sacek (University of São Paulo, São Paulo, Brazil), and Lucy Reading-Ikkanda (artist).
Syrian archaeological site
Journal Club: In Mesopotamia, early cities may have faltered before climate-driven collapse
Settlements 4,200 years ago may have suffered from overpopulation before drought and lower temperatures ultimately made them unsustainable.
Image credit: Andrea Ricci.
Steamboat Geyser eruption.
Eruption of Steamboat Geyser
Mara Reed and Michael Manga explore why Yellowstone's Steamboat Geyser resumed erupting in 2018.
Listen
Past PodcastsSubscribe
Birds nestling on tree branches
Parent–offspring conflict in songbird fledging
Some songbird parents might improve their own fitness by manipulating their offspring into leaving the nest early, at the cost of fledgling survival, a study finds.
Image credit: Gil Eckrich (photographer).

Similar Articles

Site Logo
Powered by HighWire
  • Submit Manuscript
  • Twitter
  • Facebook
  • RSS Feeds
  • Email Alerts

Articles

  • Current Issue
  • Special Feature Articles – Most Recent
  • List of Issues

PNAS Portals

  • Anthropology
  • Chemistry
  • Classics
  • Front Matter
  • Physics
  • Sustainability Science
  • Teaching Resources

Information

  • Authors
  • Editorial Board
  • Reviewers
  • Subscribers
  • Librarians
  • Press
  • Site Map
  • PNAS Updates
  • FAQs
  • Accessibility Statement
  • Rights & Permissions
  • About
  • Contact

Feedback    Privacy/Legal

Copyright © 2021 National Academy of Sciences. Online ISSN 1091-6490