Human ectoparasites and the spread of plague in Europe during the Second Pandemic
- aCentre for Ecological and Evolutionary Synthesis (CEES), Department of Biosciences, University of Oslo, N-0316 Oslo, Norway;
- bDepartment of Physiology, Institute of Basic Medical Sciences, University of Oslo, N-0317 Oslo, Norway;
- cDepartment of Computer Science, 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
See allHide authors and affiliations
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, 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).
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.
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:
The five compartments for humans that are functions of time t: susceptible (
Human ectoparasite vectors are modeled in two compartments (
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
Pneumonic Plague Model.
The direct human-to-human transmission of plague is modeled by three differential equations:
There are three compartments for humans (
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:
There are four compartments for rats (
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,
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
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
The parameters that we fitted were the transmission rates for each model (
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.
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.
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 (
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
- ↵
- Bramanti B,
- Stenseth NC,
- Walløe L,
- Lei X
- ↵
- Dennis DT, et al.
- ↵
- ↵
- ↵
- ↵
- Hufthammer AK,
- Walløe L
- ↵
- Walløe L
- ↵
- Cohn SK
- ↵
- ↵
- ↵
- ↵
- ↵
- ↵
- Smallman-Raynor M,
- Cliff AD
- ↵
- ↵
- Blanc G,
- Baltazard M
- ↵
- ↵
- Zhao WH,
- Guo M,
- Duan B,
- Su LQ
- ↵
- ↵
- ↵
- Biraben JN
- Carmichael AG
- Smith RS
- Ferrán y Clua J,
- Viñas y Cusí F,
- Grau Rd
- Creighton C
- Frandsen KE
- Alexander JT
- Burrell WH
- Peacock AD
- ↵
- ↵
- ↵
- Carrion AL
- ↵
- Boegler KA,
- Graham CB,
- Johnson TL,
- Montenieri JA,
- Eisen RJ
- ↵
- ↵
- Eisen RJ, et al.
- ↵
- ↵
- ↵
- ↵
- ↵
- Bernardo JM, et al.
- Gelman A,
- Rubin DR
- ↵
- ↵
- ↵
- ↵
- ↵
- ↵
- ↵
Citation Manager Formats
Article Classifications
- Biological Sciences
- Ecology
This article has a Letter. Please see:
- Relationship between Research Article and Letter - August 03, 2018