Author response: Eighteenth century Yersinia pestis genomes reveal the long-term persistence of an historical plague focus
2016-01-13

Article Figures and data Abstract eLife digest Introduction Results Discussion Material and methods References Decision letter Author response Article and author information Metrics Abstract The 14th–18th century pandemic of Yersinia pestis caused devastating disease outbreaks in Europe for almost 400 years. The reasons for plague's persistence and abrupt disappearance in Europe are poorly understood, but could have been due to either the presence of now-extinct plague foci in Europe itself, or successive disease introductions from other locations. Here we present five Y. pestis genomes from one of the last European outbreaks of plague, from 1722 in Marseille, France. The lineage identified has not been found in any extant Y. pestis foci sampled to date, and has its ancestry in strains obtained from victims of the 14th century Black Death. These data suggest the existence of a previously uncharacterized historical plague focus that persisted for at least three centuries. We propose that this disease source may have been responsible for the many resurgences of plague in Europe following the Black Death. https://doi.org/10.7554/eLife.12994.001 eLife digest A bacterium called Yersina pestis is responsible for numerous human outbreaks of plague throughout history. It is carried by rats and other rodents and can spread to humans causing what we conventionally refer to as plague. The most notorious of these plague outbreaks – the Black Death – claimed millions of lives in Europe in the mid-14th century. Several other plague outbreaks emerged in Europe over the next 400 years. Then, there was a large gap before the plague re-emerged as threat in the 19th century and it continues to infect humans today, though on a smaller scale. Scientists have extensively studied Y. pestis to understand its origin and how it evolved to become such a deadly threat. These studies led to the assumption that the plague outbreaks of the 14–18th centuries likely originated in rodents in Asia and spread along trade routes to other parts of the world. However, it is not clear why the plague persisted in Europe for 400 years after the Black Death. Could the bacteria have gained a foothold in local rodents instead of being reintroduced from Asia each time? If it did, why did it then disappear for such a long period from the end of the 18th century? To help answer these questions, Bos, Herbig et al. sequenced the DNA of Y. pestis samples collected from the teeth of five individuals who died of plague during the last major European outbreak of plague in 1722 in Marseille, France. The DNA sequences of these bacterial samples were then compared with the DNA sequences of modern day Y. pestis and other historical samples of the bacteria. The results showed the bacteria in the Marseille outbreak likely evolved from the strain that caused the Black Death back in the 14th century. The comparisons showed that the strain isolated from the teeth is not found today, and may be extinct. This suggests that a historical reservoir for plague existed somewhere, perhaps in Asia, or perhaps in Europe itself, and was able to cause outbreaks up until the 18th century.Bos, Herbig et al.'s findings may help researchers trying to control the current outbreaks of the plague in Madagascar and other places. https://doi.org/10.7554/eLife.12994.002 Introduction The bacterium Yersinia pestis is among the most virulent pathogens known to cause disease in humans. As the agent of plague it is an existing threat to public health as the cause of both emerging and re-emerging rodent-derived epidemics in many regions of the world (Duplantier et al., 2005; Vogler et al., 2011; Gage and Kosoy, 2005). This, and its confirmed involvement in three major historical pandemics, have made it the subject of intense study. The first pandemic, also known as the Justinian Plague, occurred from the 6th through the 8th centuries; the second pandemic spanned the 14th to the 18th centuries; and the third pandemic started in the 19th century and persists to the present day. Attempts to date the evolutionary history of Y. pestis using molecular clocks have been compromised by extensive variation in nucleotide substitution rates among lineages (Cui et al., 2013; Wagner et al., 2014), such that there is considerable uncertainty over how long this pathogen has caused epidemic disease in human populations. In addition, there has been lively debate as to whether or not it was the principal cause of the three historical pandemics (Cohn, 2008; Scott and Duncan, 2005). It is well established that extant lineages of Y. pestis circulated during the third pandemic (Achtman et al., 2004; Morelli et al., 2010), and various ancient DNA studies have now unequivocally demonstrated its involvement in the early phase of the first pandemic in the 6th century (Harbeck et al., 2013; Wagner et al., 2014) and the Black Death (1347–1351), which marks the beginning of the second plague pandemic (Bos et al., 2011; Haensch et al., 2010; Schuenemann et al., 2011). All three pandemics likely arose from natural rodent foci in Asia and spread along trade routes to Europe and other parts of the world (Morelli et al., 2010; Wagner et al., 2014). The strains associated with the first and second pandemics represent independent emergence events from these rodent reservoirs in Asia. The on-going third pandemic also originated in Asia, although genetic evidence suggests that it may be derived from strains that descended from those associated with the second wave that spread back to Asia and became re-established in rodent populations there (Wagner et al., 2014). The impact of the second pandemic was extraordinary. During this time period there were hundreds and perhaps thousands of local plague outbreaks in human populations throughout Europe (Schmid et al., 2015). It is very likely that some of these outbreaks were caused by the spread of plague via the maritime transport of humans and cargo, as was undoubtedly the case during the global spread of plague during the third pandemic (Morelli et al., 2010). However, the processes responsible for the potential long-term persistence of plague in Europe during the second pandemic are still subject to debate. It is possible that once introduced to Europe around the time of the Black Death, Y. pestis persisted there for centuries, cycling in and between rodent and human populations and being introduced or reintroduced to various regions throughout Europe (Carmichael, 2014). Another possibility is that plague did not persist long-term in European rodent populations, but rather was continually reintroduced from rodent plague foci in Asia (Schmid et al., 2015). To address this key issue in Y. pestis evolution and epidemiology, we investigated plague-associated skeletal material from one of the last well-documented European epidemics, the Great Plague of Marseille (1720–1722), that occurred in the Provence region of France at the end of the second pandemic. In particular, we compared the evolutionary relationship of these historical Y. pestis strains to those sampled from other time periods, both modern and ancient. Our results demonstrate that the strains responsible for the Black Death left descendants that persisted for several centuries in an as yet unidentified host reservoir population, accumulated genetic variation, and eventually contributed to the Great Plague of Marseille in mid-eighteenth century France. Although these lineages no longer seem to be represented within the genetic diversity of extant sampled Y. pestis, they may have been involved in additional, earlier second-pandemics in the Mediterranean region and beyond. Results Skeletal material used in this investigation was sampled from the Observance (OBS) collection housed in the osteoarcheological library of the Regional Department of Archaeology, French Ministry of Culture, medical faculty of Aix-Marseille Université. Historical records indicate that this collection represents a catastrophic burial for victims of the Great Plague of Marseille relapse in 1722 (Signoli et al., 2002). Rescue excavations carried out in 1994 unearthed the remains of 261 individuals, mostly consisting of adult males (Figure 1). Material from this collection was also used in the first PCR-based investigation of ancient Y. pestis DNA from archaeological tissue (Drancourt et al., 1998). In our study 20 in situ teeth were sampled and screened for Y. pestis DNA through a quantitative PCR assay. Amplification products were detectable for five of the 20 teeth (Table 1). Products were not sequenced. Negative controls were free of amplification products. The five putatively Y. pestis-positive samples were subject to high-throughput DNA sequencing after array-based enrichment for Y. pestis DNA. Mapping to the Y. pestis CO92 reference genome revealed a minimum of 12-fold average genomic coverage for each of the five genomes (Table 2 and Figure 2). Raw sequencing data as well as mapped reads have been deposited at the European Nucleotide Archive under accession PRJEB12163. Comparative analysis together with 133 previously published Y. pestis genomes resulted in 3109 core genome SNPs for phylogenetic analysis. No homoplasies were identified. Figure 1 Download asset Open asset Photograph of rescue excavations at l'Observance in Marseille, France. LAPMO Université de Provence – URA 164 CNRS/AFAN. https://doi.org/10.7554/eLife.12994.003 Table 1 Quantitative PCR data for the pla gene obtained for the 19 Observance teeth extracted. Values below std curve detection are highlighted in grey. https://doi.org/10.7554/eLife.12994.004 PCR Blks:Copies per μlPCR Blk 10PCR Blk 20PCR Blk 30PCR Blk 40Extraction Blks:Copies per μlBLk 10Blk 20Blk 30Blk 40Extracts:SampleCopies per μl101a0101b010201030107391080109011050116105117a0117b1118012301241421250126013001341137301440 Table 2 Mapping statistics. https://doi.org/10.7554/eLife.12994.005 Sample IDRaw read pairsPreprocessed readsMapped reads% mappedDeduppedDuplication factorFold coverage% coveredOBS10714,646,71015,395,2871,850,73412.02%851,0962.1712.2084.79OBS1109,879,25710,458,5711,581,15815.12%855,9881.8514.6190.63OBS11616,858,27317,876,8832,931,55416.40%1,363,4152.1520.0492.56OBS12448,797,60252,345,8844,972,3479.50%973,7635.1113.3485.37OBS137120,018,142129,095,43412,669,6259.81%1,583,5968.0024.4092.76 Figure 2 Download asset Open asset Coverage plots for all reconstructed Observance core genomes. Inner ring: GC content from low (brown, <30%) to high (orange, >55%). Outer concentric rings: from inner to outer, coverage plots of OBS107, 110, 116, 124, and 137 on a logarithmic scale. Axes are at 30x and 100x. Dots: SNPs (black = non-synonymous, grey = synonymous, white = intergenic). Outer ring (larger dots) is SNPs shared by all five strains. Inner rings (smaller dots) are associated with the strain immediately outside and are SNPs that are not shared by all five strains (but may be shared by 2–4 strains). Krzywinski, M. et al. Circos: an Information Aesthetic for Comparative Genomics. Genome Res (2009) 19:1639–1645 https://doi.org/10.7554/eLife.12994.006 The phylogenetic position of the five OBS sequences on a unique branch of the Y. pestis phylogeny (lineages shown in red in Figure 3), yet seemingly derived from those strains associated with the Black Death, both confirms its authenticity and reveals that it is likely a direct descendent of strains that were present in Europe during the Black Death. Although the topological uncertainty (i.e. trichotomy) among the second pandemic strains (London and Observance) suggests that there was a radiation of Y. pestis lineages at the time, it is notable that the OBS sequences do not possess the additional derived position that is present in both the 6330 London strain and the SNP-typed individual from Bergen Op Zoom (defined as the 's12' SNP by Haensch et al., 2010). Indeed, the tree is striking in that there is clear phylogenetic support (88% bootstrap values, Figure 3A) for 6330 clustering with the Branch 1 strains, including those that gave rise to the third (modern) plague pandemic. Finally, although evolutionary rates in Y. pestis are notoriously variable (Cui et al., 2013; Wagner et al., 2014), the comparatively late time period from which our samples derive and the long branch leading to these OBS sequences suggests that this lineage co-existed for an extended period with the Branch 1 strains responsible for all later human plague outbreaks outside of China. Figure 3 Download asset Open asset Phylogenetic tree for Yersinia pestis. (A) Maximum likelihood phylogeny of Y. pestis genomic SNPs showing the position of the Observance (OBS) lineage (red) relative to those of 130 modern (black) and three ancient strains (Black Death in green and Justinian Plague in blue). Modern strains from the third pandemic are shown in purple to highlight their close genetic relatedness. Monophyletic groups of sequences have been collapsed to improve clarity and are shown as triangles. The tree is rooted using single strain of Y. pseudotuberculosis (IP32953), with all derived SNPs removed to assist scaling, with branch lengths reflecting the number of nucleotide substitutions/SNP site. The length of the branch leading to the 0.PE3_Angola (AFR) lineage was reduced because its excessive length adversely affected the scaling of the tree. Location abbreviations are as follows: CNH (China), GEO (Georgia), FSU (Former Soviet Union), MNG (Mongolia), NPL (Nepal), IRN (Iran), AFR (Africa), USA (United States of America), MDG (Madagascar), IND (India). (B) Expanded phylogeny schematic to show the relative positions of the Black Death and the Observance lineages. Numbers on branches correspond to SNPs. https://doi.org/10.7554/eLife.12994.007 Identification of DFR4 in OBS strains Using the pan-array design, we were able to identify a ~15 kb Genomic Island (GI) in the OBS strains as previously observed in the Justinian and Black Death strains (Wagner et al., 2014) and referred to as DFR 4 (difference region 4). Notably, this region has been lost in some Y. pestis strains such as CO92. This island is a striking example of the decay common to the highly plastic Y. pestis genome. Discussion The Observance lineage of Y. pestis identified by our analysis of material from the Great Plague of Marseille (1720–1722) is clearly phylogenetically distinct and has not been identified in any extant plague foci for which full genome data are available. That this lineage is seemingly confined to these ancient samples suggests that it is now extinct, or has yet to be sampled from extant reservoirs. Based on available data, the Y. pestis strains most closely related to the Observance lineage are those obtained from other ancient DNA studies of victims of the Black Death in London and other areas (Haensch et al., 2010; Bos et al., 2011) that occurred some 350 years earlier. The most important conclusion from the current study is that Y. pestis persisted and diversified in at least two lineages throughout the second pandemic, one represented by our Observance strain, and the other by London St. Mary Graces individual 6330. The absence of currently sequenced modern Y. pestis genomes that share an immediate phylogenetic relationship with the Observance lineage suggests that it may no longer exist. The length of the branch leading to the Observance lineage is consistent with a long history of circulation. Of note, this new lineage does not share the single derived mutation that is common to the London St. Mary Graces genome and the Bergen op Zoom plague victim from the Netherlands (Haensch et al., 2010); hence, multiple variants of Y. pestis evidently circulated in London and elsewhere in the 14th century, and at least one of these variants followed a path that differs from the lineage that later gave rise to the 18th century plague epidemic in Provence. Since the novel lineage of Y. pestis identified here was obtained from an active Mediterranean port city that served as a main commercial hub and entry point into Western Europe from various origins, the precise location of the disease's source population cannot be easily determined. Several scenarios can explain these data. First, all Y. pestis diversity documented in the Black Death and later plague outbreaks could have stemmed from foci in Asia, where European epidemics were the result of successive introductions from this distant, albeit prolific, source population (Schmid et al., 2015). This model would not require the back migration of plague into Asia for the third pandemic proposed elsewhere (Wagner et al, 2014), since all preexisting diversity would already be present. A multiple wave theory has previously been proposed to explain the presence of a reportedly different 14th century lineage in the Low Countries (defined by their single 's12' SNP) compared to those circulating in England and France (Haensch et al., 2010). However, the presence of this branch 1 position in both the ancestral and subsequently the derived state in two temporally close London outbreaks (Bos et al., 2011; Figure 3B) suggests that it became fixed during the Black Death or its immediate aftermath, rather than having been introduced via a separate pulse from an Asian focus. To date, all of the Y. pestis material examined from the second pandemic has been assigned to Branch 1 in the Y. pestis phylogeny, indicating a close evolutionary relationship. As the large "'big bang' radiation of Y. pestis at the beginning of the second pandemic (Cui et al., 2013) established at least four new lineages, it seems unlikely that independent waves of plague from Asia into Europe during the Black Death would have involved strains that group on the same branch and differ by so few positions. As China harbours many branch 1 descendent lineages that also share this 's12'" SNP, our observation adds weight to the notion that the branch 1 lineage of Y. pestis spread east sometime during the second pandemic after it entered Europe, and subsequently became established in one or more East Asian reservoir species (Wagner et al., 2014). Here it remained before radiating to other locations in the late 19th century, giving rise to the current worldwide 'third wave' Y. pestis pandemic. In our view, a far more plausible option to account for the distinct position of the Observance lineage is that a Y. pestis strain responsible for the Black Death became established in a natural host reservoir population within Europe or western Asia. Once established, this Y. pestis lineage evolved locally for hundreds of years and contributed to repeated human epidemics in Europe. A similar scenario is possible for the lineage(s) identified in the London St. Mary Graces and Bergen op Zoom material: like Observance, both descend from the Black Death and together they represent at least one source of European plague distinct from that which later caused the Great Plague of Marseille. Plague's temporary persistence somewhere on the European mainland is a compelling possibility (Carmichael, 2014). Alternatively, reservoir populations located in geographically adjacent regions, such as in western Asia or the Caucasus, may have acted as regular sources of disease through successive westerly pulses during the three century-long second pandemic (Schmid et al., 2015). Since our study used material exclusively from a highly active port with trade connections to many areas in the Mediterranean and beyond, it is impossible to identify a source population for the Great Plague of Marseille at our current resolution. Our analysis reveals a previously uncharacterized plague source that could have fuelled the European epidemics from the time of the Black Death until the mid-eighteenth century. This parallels aspects of the third pandemic, that rapidly spread globally and established novel, long-lived endemic rodent foci that periodically emerge to cause human disease (Morelli et al., 2010). The current resolution of Y. pestis phylogeography suggests that the historical reservoir identified here may no longer exist. More extensive sampling of both modern rodent populations and ancient human, as well as rodent, skeletal remains from various regions of Asia, the Caucasus, and Europe may reveal additional clues regarding past ecological niches for plague. The reasons for the apparent disappearance of this historic natural plague focus are, however, unknown, and any single model proposed for plague's disappearance is likely to be overly simplistic. Pre-industrial Europe shared many of the same conditions that are correlated with plague dissemination today including pronounced social inequality, and poor sanitation coupled with high population density in urban centers (Vogler et al., 2011; Barnes, 2014). This constellation of anthropogenic factors, along with the significant social and environmental changes that occurred during the Industrial Revolution must be considered alongside models of climate and vector-driven dynamics as contributing to the rapid decline in historical European plague outbreaks, where genomic data provide but one piece of critical information in this consilient Material and methods this study 20 in situ teeth were at the in a ancient DNA and to of from the inner of the or the DNA were carried out following elsewhere et al., and et al., and of Y. pestis DNA was through an established quantitative PCR for the pla gene et al., 2011). DNA for the five putatively Y. pestis-positive each from a separate were into and 2010; et al., following an DNA et al., to the most common of ancient DNA were with to a of 20 and in over two 1 et al., were with three the CO92 Y. pestis reference with additional regions from Y. pestis Y. pseudotuberculosis and Y. pseudotuberculosis that are in CO92. including those from were sequenced on one of an Raw reads were reads were as elsewhere et al., and reads were subsequently for a length of reads were mapped to the CO92 reference genome using the and with and a of in SNP reads were at three and the was used as our SNP 1). the of reads were as and SNP was to a by Bos et al. using the of the Genome et al., 2011). for 133 previously published Y. pestis strains were 2). A was used for the of the SNPs were called with a coverage of and a of for the SNP as identified by Morelli et al. genomic regions as by et al. as well as regions, and were from the SNP analysis. SNPs were using et al., with The and region was to reads were and using a minimum of to than were out using reads and reads were and mapped to the CO92 reference using the and with the in Wagner et al. were collapsed using a which reads with and were into and SNPs with at least coverage and were SNPs were for reads were with et al., 2014) and the CO92 reference using SNPs were called with the in and were by minimum coverage and SNPs in genome were identified by a direct the reference using et al., 2002). These methods were by the SNP Phylogenetic analysis of the Y. pestis A phylogenetic tree of Y. pestis strains was using the likelihood available in et al., 2010), the model of nucleotide substitution available from the on and a of and The of individual was using bootstrap using the same substitution model and as The phylogeny data from the five Observance (OBS) strains, one from Black Death victims who died in London between and of strains and one from a outbreak in – one from the Plague of and 130 modern Y. pestis strains. A single of Y. pseudotuberculosis was used as an to the although with all derived SNPs removed to assist Identification of DFR4 in OBS strains a and reads were mapped to the Y. pestis at the DFR 4 region position and the unique were also mapped back to the to of the region using read length data. All were compared using and to identify of the that were unique to compared to CO92 and were by reads from our ancient a Open The current on this is being References Morelli Vogler Wagner and history of the plague Yersinia pestis of the of and the 1994 outbreak in a for data Bos Herbig A Schuenemann genomes reveal as a source of new world human Bos Schuenemann A genome of pestis from victims of the of and detection of in in ancient DNA Plague persistence in western a in the the Black