Category Archives: bioarchaeology

Evolutionary Clues in 17th-Century Smallpox Genome

By Michelle Ziegler

Smallpox is one of those diseases long believed to have an ancient pedigree, the suspected culprit of legendary epidemics. Yet, so far, smallpox hasn’t made a big impression in ancient DNA surveys. If it was truly endemic throughout the Old World before 1492, so much so that it pops up in the New World almost immediately after contact, it’s odd that it has not been more prominent in ancient DNA surveys. Be ready for a smallpox paradigm shift and reexamination of its reputed history.

In December, Ana Duggan, Maria Perdomo, and the McMaster Ancient DNA Centre team announced the first full ancient smallpox genome isolated from a mummified 17th-century child in Vilnius, Lithuania. Radiocarbon dates of the child place him or her in the mid 17th century (est. c. 1654) in the midst of dated smallpox epidemics from all over Europe.

gr1.jpg
Fig 1 (Duggan et al, 2016): Left: Distribution of smallpox records in Europe. Right upper: Dominican Church of the Holy Spirit, Vilnius. Lithuania. Right lower: Crypt containing the child’s remains.

Their finding was unexpected. They were not looking for smallpox at all; the child had no observable lesions. They were hoping to find JC polyomavirus, of particular interest to one of the co-authors, and so they first enriched the specimen for this virus (McKenna 2016). After detecting variola virus (VARV), smallpox, instead they then enriched for VARV to confirm the initial signal.

duggan_cbsmallpox_finalMore than just confirming the signal, they were able to reconstruct the entire genome producing the entire sequence at an average depth of 18X. The surprising child had more revelations in his or her viral sequence. The sequence is ancestral to all existing reference strains. This is consistent with short stretches of aDNA amplified from 300-year-old frozen Siberian remains (Biagini et al, 2012). Unfortunately, the sequences from 2012 were not distinctive enough from the new Lithuanian sequence to give phylogicial resolution between them. Oddly, the frozen Siberian remains also lacked smallpox skin lesions with one showing signs of pulmonary hemorrhages.

Its ancestral position in the phylogeny suggests that a severe bottleneck occurred before c. 1654. As Duggan et al (2016) remark, vaccination would cause a very strong bottleneck, but this occurs after 1654 and there is new diversity among the descendent reference specimens producing two major clades. Yet to be determined is the evolutionary effect of extensive variolation practices in the early modern period. In contrast to vaccination, variolation is a form of intentional smallpox transmission that sometimes went horribly wrong.

Evolution continues unabated. The molecular clock is consistent among the 20th-century specimens and the latest aDNA from the Lithuanian child. The two clades of smallpox collected from 20th-century specimens diverged from each other sometime around the mid-17th century after vaccination began. Interestingly, the less virulent Variola minor strain is not predicted to have emerged from Variola major clade P-II until the mid-19th century.

screen-shot-2016-12-28-at-2-39-12-pm
Evolutionary history of Variola (Duggan et al, 2016)

It’s not entirely surprising that smallpox, a highly transmissive human-only virus, has a relatively recent last common ancestor; other viruses like measles do as well. Measles last common ancestor is probably in the early 20th century (Furuse, Suzuki & Oshitani, 2010). The dominance of the 1918 influenza strain in recent influenza phylogeny is another example; incomplete because influenza swaps genes with influenza viruses that are circulating primarily in birds, but also in swine and earlier in equines (Taubenberger & Morens, 2005). Improved transmission strains are likely to out-compete strains with a lower transmission rate if they achieve a global spread. For some viruses, though not necessarily all, improved transmissibility and virulence go hand in hand. So, in the end, the relatively recent last common ancestor says more about its global transmissibility than anything else.

The antiquity of the virus needs two components to estimate. The molecular clock must be steady, and it is so far (though this could change with more ancient specimens), and a near relative ‘out group’, related strains outside the Variola clade (a branch of the larger genetic tree). One potential problem here is that as transmissibility improves the clock may speed up. The speed of the clock is determined by the reproduction rate. The relatively steady clock back to this 17th-century specimen suggests that the transmission rate was pretty steady — after the evolutionary/transmission leap that swept aside other Variola strains. The inactivation of several orthopoxvirus genes in smallpox that are functional in vaccinia (used in smallpox vaccines), camelpox, and teterapox  (the ‘out groups’ used) may suggest that one or more of these genes had been protective. When the genes were inactivated, smallpox probably became a much more dangerous virus to humans.

Historical epidemiology suggests that there was once more variation in the virulence of smallpox epidemics.  Securely identifying smallpox epidemics in the historical record is much harder than is generally assumed, and it is harder yet to make a claim for a significant demographic impact prior to the Renaissance (Carmichael & Silverman, 1987). This is the problem with theories that smallpox was the cause of the second century Antonine plague and then failed to cause an epidemic with a major demographic effect for many centuries. I find this very hard to believe. Additionally, the infamous smallpox epidemics in the New World are now also be reevaluated in ways that diminish smallpox’s toll and add in a wide variety of contributing factors to produce a colonization syndemic. This has most recently been summarized in essays collected in Beyond Germs: Native Depopulation in North America (2015).

One other observation from these studies: All ancient smallpox DNA to date has been extracted from mummy tissue, not a tooth or bone. This may point toward one of the limitations of ancient DNA pathogen surveys that currently use primarily teeth. Since neither mummy had visible smallpox lesions, smallpox should be considered a possibility in any mummy.


References

Duggan, A. T., Perdomo, M. F., Piombino-Mascali, D., Marciniak, S., Poinar, D., Emery, M. V., et al. (2016). 17th Century Variola Virus Reveals the Recent History of Smallpox. Current Biology, 1–7. http://doi.org/10.1016/j.cub.2016.10.061

Biagini, P., Thèves, C., Balaresque, P., Géraut, A., Cannet, C., Keyser, C., et al. (2012). Variola virus in a 300-year-old Siberian mummy. The New England Journal of Medicine, 367(21), 2057–2059. http://doi.org/10.1056/NEJMc1208124

McKenna, Maryn (8 Dec 2016) Child Mummy Found with Oldest Known Smallpox Virus. National Geographic. (online)

Carmichael, A. G., & Silverstein, A. M. (1987). Smallpox in Europe before the seventeenth century: virulent killer or benign disease? Journal of the History of Medicine and Allied Sciences, 42(2), 147–168.

Furuse Y, Suzuki A, & Oshitani H (2010). Origin of measles virus: divergence from rinderpest virus between the 11th and 12th centuries. Virology journal, 7 PMID: 20202190

Taubenberger, J. K., & Morens, D. M. (2005). 1918 Influenza: the mother of all pandemics. Emerging Infectious Diseases, 12(1), 15–22. http://doi.org/10.3201/eid1201.050979

Beyond Germs: Native Depopulation in North America. Edited by Catherine Cameron, Paul Kelton, and Alan Swedlund. University of Arizona Press, 2015.

Plague in 6th century Aschheim and Altenerding, Bavaria

Since I last wrote about Bavaria, the aDNA centers have been busy. With the accepted manuscript of the second new paper available this past week, its time for an update. The fourth paper on Aschheim not only confirmed the first three, but it also produced the first full genome of Yersinia pestis for the Plague of Justinian (Wagner et al, 2014). This paper also confirmed the Bavarian strain’s placement in the phylogeny of Y. pestis. The availability of the first full genome will primarily be important for comparison to newly discovered samples from elsewhere. Using newer technology, the newest paper refined some of the Aschheim sequence and produced a full genome of Y. pestis from a woman buried at Altenerding, about 20 km from Aschheim (Feldman et al, 2016). Radiocarbon dating from both sites places the epidemic in the mid-sixth century; it can not differentiate which specific epidemic ‘wave’.  The Altenerding epidemic was from the same Y. pestis lineage as Aschheim proving that this was a regional epidemic, possibly the same epidemic event. The phylogeny for the first pandemic is still based on a single epidemic from one geographic region, so the time is not yet ripe to use the phylogeny to tell inform us on the transmission or route of the pandemic.

6th cent Bavaria
Map of Roman Bavaria showing the Roman roads with Aschheim and Altenerding marked. The half circle/mound mark designates Roman villas. (modified from the Pelagios project)

It is, however,  time to start thinking a little more about the environment of these sites. They are both located on the Munich gravel plain, foreland (foothills) north of the Alps. Aschheim is located closer to the Alps at an elevation of 500 meters with Altenerding 20 km further north at a lower elevation in small valley formed by a tributary of the River Isar. The Roman road running horizontally across the map runs west to Augsburg, the capital of the Roman province of Raetia Secunda and east to the city of Batavia, a colony in the province of Noricum. The road running by Altenerding would take traffic eventually north toward Regensburg (Casta Regina).

Large water feature is Speichersee lake with a man-made 20th century reservoir used to power hydroelectric plants and serve some of the water needs of the Munich region. As far as I can tell, none of this would have been present in the Late Antique period. The River Isar is the green line to the west of both sites. Munich will later be founded where the road crosses the river from monastic land in about 1158. There was nothing special at the river crossing in the sixth century. Although the road crosses the river, there is no indication of a Roman bridge on the map.

Both Aschheim and Altenerding are located in what would have been the province of Raetia II. While they are along Roman roads, this would have been a rural area. Both Aschheim and Altenerding were sites of Roman villas and Dornach near Aschheim was a small settlement. How much of this would have been occupied and further developed (or not) after the Roman army left is unclear. The cemetery at Altenerding is triple the size of Aschheim. Yet, there is reason to think that Aschheim was hit harder by the plague and based on the carbon dates of graves with some molecular plague signal, probably more than once. Michael McCormick (2015:83) suggests that the Aschheim cemetery gathered graves from a dispersed settlement that probably had fewer than 70 people at any one time.

A living history museum in Munich area at Kirchheim has reconstructed typical buildings from the early medieval Merovingian period. Although this area was nominally under Merovingian Frankish hegemony there is little specifically Frankish about the archaeology. They were all wooden construction. Below is a picture of a sunken pit building, an ‘out building’ and a long house.

1280px-Bajuwarenhof_Kirchheim_Übersicht_2012-08-05
Reconstruction of 6th-7th century Bavarian buildings at Kirchheim in the Munich district close to Aschheim. (Photo by Leporollo, Wikipedia CC3.0)

Continue to think of the Plague of Justinian in Constantinople and Pelusium, it was surely there. Just remember that most of its geographic spread may have looked more like this picture.

References:

Feldman, M., Harbeck, M., Keller, M., Spyrou, M. A., Rott, A., Trautmann, B., et al. (2016). A high-coverage Yersinia pestis Genome from a 6th-century Justinianic Plague Victim. Molecular Biology and Evolution, 1–31. [Accepted manuscript]

McCormick, M. (2015). Tracking mass death during the fall of Rome’s empire (I). Journal of Roman Archaeology, 28, 325–357.

Wagner, D. M., Klunk, J., Harbeck, M., Devault, A., Waglechner, N., Sahl, J. W., et al. (2014). Yersinia pestis and the Plague of Justinian 541–543 AD: a genomic analysis. The Lancet Infectious Diseases, 14(4), 1–8. http://doi.org/10.1016/S1473-3099(13)70323-2

cropped-venice-plague.jpg

Plague Dialogues: Monica Green and Boris Schmid on Plague Phylogeny (II)

Monica H. Green (monica.green@asu.edu,@MonicaMedHist) is a historian of medieval medicine. An elected Fellow of the Medieval Academy of America, she teaches both global history and the global history of health. She was the editor in 2014 of Pandemic Disease in the Medieval World: Rethinking the Black Death, the inaugural issue of a new journal, The Medieval Globe.

Boris Schmid (@BorisVSchmid) is a theoretical biologist at the University of Oslo, Norway, and specializes in disease ecology and epidemiology. He recently described a link between climate fluctuations in medieval Central Asia and what looks like repeated introductions of plague into Europe’s harbors, a hypothesis that can be tested by the analysis of ancient DNA samples of Y. pestis. He works in a multidisciplinary team of theoreticians, archeologists, microbiologists and historians, led by Nils Chr. Stenseth.


Boris:

In our previous blog post, Monica and I discussed how different lineages of plague – Yersinia pestis – collected their own genetic signature (SNP profile) as they diversified from a common ancestor. Monica also summarized in broad terms what ancient DNA samples of Y. pestis (extracted from plague victims) are now available from the initial Black Death outbreak and how they are related, using the latest plague studies of Haensch, Bos and Spyrou.

In this blog post Monica will delve into the nitty gritty details of these aDNA plague studies, and give an example of how to transform those details into a new understanding of the past rodent reservoirs and global mobility of plague, one of the deadliest diseases of our collective past. And I close the post by reflecting on the potential of aDNA to connect the fields of history and biology.

Monica:

Thanks, Boris. It might be important to remind readers that we don’t have any aDNA evidence from past rodent populations yet. All the samples to date have been retrieved from human victims. But the SNP study that Seifert et al. published earlier this year from samples in Brandenburg, an inhumation from the time of the 30 Years War (1618-1648); the whole genome study that Bos and Herbig et al. also published this year, reporting on the samples from 18th-century Marseille; and the sample from Ellwangen included in the new study by Spyrou et al., all document that Branch 1A (see tree in our previous post) “focalized,” that is, it set up shop in some rodent population(s) and happily continued to proliferate for another 400 years. But all that happened, it is clear, separately from what was going on with Branch 1B, what I have taken to call the pestis secunda.

In their most recent study, the Tübingen/Jena team headed by Krause give us further insight into the early stages of Branch 1B. The beginning of Branch 1B was first documented in the 2011 London study, though it was only earlier this year that I realized that London sample 6330 likely dates from the 1360s and does not come from the initial Black Death outbreak. (It comes from a different burial ground, St Mary Graces.) In Bos and Herbig et al. 2016, it was reported that sample 6330 differed from the 1348-50 London Y. pestis genome by two SNPs. In the present study, interestingly, Spyrou et al. report something slightly different. Sample 6330 does indeed differ from the London Black Death genome by two SNPs (p3 and p4), but a third SNP in sample 6330 they are reporting here for the first time (p5) seems to be unique, a ‘G’ to ‘T’ switch at position 4,301,295 not found in any other historic genome or in the reference strain, CO92. (Spyrou et al. did not include London 6330 on their Table 1, so we offer a modified version of it here in fig. 3.)

CORRECTED fig03 for Contagions blog, (a) Spyrou et al 2016 fig02B Y pestis phylogenetic tree (detail of origins of Branches 1A and 1B), (b) Table 1 with 6330 SNPs added (06272016)
Fig. 3: (a) Detail of Spyrou et al. 2016, fig 2B: Yersinia pestis phylogeny – SNPs distinguishing Branches 1A and 1B; (b) Spyrou et al. 2016, Table 1, modified with data on London sample 6330 drawn from Spyrou et al. 2016, Table S4, SNP table. The SNPs unique to London 6330 and Bolgar City are highlighted in yellow.

 

First of all, we might say that those three SNPs are significant for the time gap they suggest between the Black Death and the 2nd wave of plague to hit London in 1361-63; for the sake of argument, we’ll say “21 years,” to use your formula (3 x 7 years), Boris. But think about the implications of that: plague arrives in London at the end of 1348 as a new disease, and a new strain (with 3 new SNPs) is causing a major new outbreak in 1361-63, which is when this burial seems to date from. 21 years have not passed since the previous outbreak. So what gives? Obviously, the “time-to-SNP” calculus we’re using is an average, not an absolute. But it does make us stop and wonder: did all this really happen so fast? And did it really happen in western Europe?

Which brings us to the Bolgar City sample. It, clearly, is a “descendant” of the same strain as London 6330: it has the two new SNPs, p3 and p4. (It doesn’t have that unique p5 SNP of London 6330, but that may have arisen as little as three days before this person died. We cannot attach any evolutionary significance to it until we see it documented somewhere else.) But note this: the Bolgar City strain has evolved further. It now has the p6 SNP that will define all the rest of Branch 1B and it has its own unique SNP, p7. And again, we have a problem of time compression: the Bolgar City sample (if we can trust the dating of the coins which were said to have been found with the body) may date as early as the late 1360s.

Remember what we need to have an outbreak in humans from a new lineage of Y. pestis: not simply does that new lineage have to arise from a single change in a single cell of Y. pestis, but that new SNP needs to proliferate enough in a reservoir rodent population to cause a new epidemic in humans. So looking over all these SNPs, p1-p7, we can see that they cluster into two “founder effect” phenomena: one that creates the initial Black Death lineage (Branch 1A) and one that creates the pestis secunda lineage (Branch 1B).

Where did those two lineage foundations happen? Let’s go back to Caffa, the “hurling bodies over the walls” scenario. Clearly, if we can believe that story (and remember, we have only one account of it, and that from a non-eye-witness), it tells of an already proliferating plague outbreak. By October 1346, Y. pestis was multiplying by the millions in rats and mice and rodents of whatever kind that lived in and around Caffa.

One, and only one, of those gazillions of Caffese offspring gave rise to Branch 1B. It, too, needed to find a place to set up shop and proliferate to make gazillions of (nearly) identical copies. And where was that place? Was it (as Krause seemed to imply in his April lecture) in London? Maybe it was in or near Bergen op Zoom (NL), where we find a sample with the same SNP profile as the London 6330 sample (Haensch et al. 2010)? Or was it near the same place where Branch 1A had already established its original home, before it reached the Black Sea? Haensch et al. had already proposed in 2010 a “northern” route for the introduction of the pestis secunda strain that reached the Netherlands. I’ll admit, I was skeptical for the longest time. But now I see that this possibility might bear more analysis. At the very least, the question shifts our focus away from western Europe and back to the areas around the Black and Caspian Seas. And that’s exactly where our Bolgar City sample is from, the one that is already showing two SNPs of further evolution beyond London 6330 but might not be a whole lot younger than it. As we said, jetting out of Heathrow wasn’t yet an option in the 14th century. But there was plenty of activity in these central Eurasian areas dominated by the Mongol Golden Horde to connect lots of rodent reservoirs to a bacterium looking for a new place to call home.

Boris:

Thanks Monica! The amount of information that follows from a few different nucleotides between aDNA samples is quite amazing, and learning how to interpret this data historically is rightly one of the transformative processes now happening in Biology (and if I say so, in Medical History as well).

Monica’s interpretation of plague’s past mobility is based on the same genetic data as the one sketched out in Spyrou 2016, and highlights the challenge of interpreting ancient DNA, given that the ancient DNA sequences of plague are still so sparsely sampled across time and space. One thing that strikes me as especially important is how much the argument of “favor the most simple, parsimonious explanation” changes based on whether you think of plague largely in terms of a human epidemic (which Wagner 2014, and by extension Spyrou 2016 appear to do), or as a disease that spread through human and wildlife both, as Monica and I do. If you include the possibility of new wildlife reservoirs of plague (and plague has created numerous new wildlife reservoirs in time), say near Bolgar City, the logic of how plague moved across Eurasia changes.

As more aDNA data becomes available, it will be very interesting to see the geographic range that a lineage of plague bacteria can spread without collecting changes in its SNP profile. Once we have a good idea of that, and a more complete view of the SNP profiles that existed during the past pandemics, SNP profiles might be used to shed light on the actual source of a historic plague outbreak, and thus offer an independent way of checking the reliability of historic sources that blame particular smugglers, ships, refugees or clothing as the source of a plague outbreak.

Wow, thanks, Monica, for this great discussion. This is an example on how history and biology can intertwine, and while we are all waiting for more revelations from aDNA and historic sources, it seems prudent to start more interactions between historians and biologists. There is an inherent bias to doubt your own data too much, and trust another fields’ data too blindly, leading to mistakes at both sides: we blindly pick some historic report as authoritative, or put too much faith in a report on the (in-) efficiency of plague transmission by different flea species, whilst a single mutation that causes the loss of a gene can have drastic effects on how well the disease transmits (Hinnebusch, 2016). The only practical way to avoid falling into such pitfalls is by investing in cross-talk between scholars of the humanities and natural sciences!

Monica:

Thank you, Boris. This was great. And very special thanks to Michelle Ziegler, for hosting our discussion on her super blog, Contagions.

References

Bos, K. I., Herbig, A., Sahl, J., Waglechner, N., Fourment, M., Forrest, S. A., et al. (2016). Eighteenth century Yersinia pestis genomes reveal the long-term persistence of an historical plague focus. eLife, 5, 17837. http://doi.org/10.7554/eLife.12994

Cui, Y., Yu, C., Yan, Y., Li, D., Li, Y., Jombart, T., et al. (2012). Historical variations in mutation rate in an epidemic pathogen, Yersinia pestis. Proceedings of the National Academy of Sciences, 110(2), 577-582.  http://doi.org/10.1073/pnas.1205750110/-/DCSupplemental/sd01.xls

Haensch, S., Bianucci, R., Signoli, M., Rajerison, M., Schultz, M., Kacki, S., et al. (2010). Distinct Clones of Yersinia pestis Caused the Black Death. PLoS Pathogens, 6(10), e1001134. http://doi.org/10.1371/journal.ppat.1001134.t001

Hinnebusch, B. J., Chouikha, I., & Sun, Y.-C. (2016). Ecological Opportunity, Evolution, and the Emergence of Flea-borne Plague. Infection and Immunity, IAI.00188–16–31. http://doi.org/10.1128/IAI.00188-16

Krause, Johannes (4-12-2016)  Oral Presentation #S577:  Ancient pathogen genomics: what we learn from historic pandemics. European Congress  of Clinical Microbiology and Infectious Diseaseshttp://eccmidlive.org/#resources/ancient-pathogen-genomics-what-we-learn-from-historic-pandemics

Seifert, L., Wiechmann, I., Harbeck, M., Thomas, A., Grupe, G., Projahn, M., et al. (2016). Genotyping Yersinia pestis in Historical Plague: Evidence for Long-Term Persistence of Y. pestis in Europe from the 14th to the 17th Century. PLoS ONE, 11(1), e0145194-8.  http://doi.org/10.1371/journal.pone.0145194

Spyrou, M. A., Tukhbatova, R. I., Feldman, M., Drath, J., Kacki, S., de Heredia, J. B., et al. (2016). Historical Y. pestisGenomes Reveal the European Black Death as the Source of Ancient and Modern Plague Pandemics. Cell Host and Microbe, 19(6), 874–881.  http://doi.org/10.1016/j.chom.2016.05.012

Wagner, D. M., Klunk, J., Harbeck, M., Devault, A., Waglechner, N., Sahl, J. W., et al. (2014). Yersinia pestis and the Plague of Justinian 541–543 AD: a genomic analysis. The Lancet Infectious Diseases, 14(4), 1–8. http://doi.org/10.1016/S1473-3099(13)70323-2