Ancient Yersinia pestis genomes lack the virulence-associated YpfΦ prophage present in modern pandemic strains.

Yersinia pestis is the causative agent of at least three major plague pandemics (Justinianic, Medieval and Modern). Previous studies on ancient Y. pestis genomes revealed that several genomic alterations had occurred approximately 5000-3000 years ago and contributed to the remarkable virulence of this pathogen. How a subset of strains evolved to cause the Modern pandemic is less well-understood. Here, we examined the virulence-associated prophage (YpfΦ), which had been postulated to be exclusively present in the genomes of strains associated with the Modern pandemic. The analysis of two new Y. pestis genomes from medieval/early modern Denmark confirmed that the phage is absent from the genome of strains dating to this time period. An extended comparative genome analysis of over 300 strains spanning more than 5000 years showed that the prophage is found in the genomes of modern strains only and suggests an integration into the genome during recent Y. pestis evolution. The phage-encoded Zot protein showed structural homology to a virulence factor of Vibrio cholerae. Similar to modern Y. pestis, we observed phages with a common origin to YpfΦ in individual strains of other bacterial species. Our findings present an updated view on the prevalence of YpfΦ, which might contribute to our understanding of the host spectrum, geographical spread and virulence of Y. pestis responsible for the Modern pandemic.

Analysis of Genomic DNA from Medieval Plague Victims Suggests Long-Term Effect of Yersinia pestis on Human Immunity Genes.

Pathogens and associated outbreaks of infectious disease exert selective pressure on human populations, and any changes in allele frequencies that result may be especially evident for genes involved in immunity. In this regard, the 1346-1353 Yersinia pestis-caused Black Death pandemic, with continued plague outbreaks spanning several hundred years, is one of the most devastating recorded in human history. To investigate the potential impact of Y. pestis on human immunity genes, we extracted DNA from 36 plague victims buried in a mass grave in Ellwangen, Germany in the 16th century. We targeted 488 immune-related genes, including HLA, using a novel in-solution hybridization capture approach. In comparison with 50 modern native inhabitants of Ellwangen, we find differences in allele frequencies for variants of the innate immunity proteins Ficolin-2 and NLRP14 at sites involved in determining specificity. We also observed that HLA-DRB1*13 is more than twice as frequent in the modern population, whereas HLA-B alleles encoding an isoleucine at position 80 (I-80+), HLA C*06:02 and HLA-DPB1 alleles encoding histidine at position 9 are half as frequent in the modern population. Simulations show that natural selection has likely driven these allele frequency changes. Thus, our data suggest that allele frequencies of HLA genes involved in innate and adaptive immunity responsible for extracellular and intracellular responses to pathogenic bacteria, such as Y. pestis, could have been affected by the historical epidemics that occurred in Europe.

Neolithic Yersinia pestis infections in humans and a dog.

Yersinia pestis has been infecting humans since the Late Neolithic (LN). Whether those early infections were isolated zoonoses or initiators of a pandemic remains unclear. We report Y. pestis infections in two individuals (of 133) from the LN necropolis at Warburg (Germany, 5300-4900 cal BP). Our analyses show that the two genomes belong to distinct strains and reflect independent infection events. All LN genomes known today (n = 4) are basal in the phylogeny and represent separate lineages that probably originated in different animal hosts. In the LN, an opening of the landscape resulted in the introduction of new rodent species, which may have acted as Y. pestis reservoirs. Coincidentally, the number of dogs increased, possibly leading to Y. pestis infections in canines. Indeed, we detect Y. pestis in an LN dog. Collectively, our data suggest that Y. pestis frequently entered human settlements at the time without causing significant outbreaks.

A melting pot of Roman dogs north of the Alps with high phenotypic and genetic diversity and similar diets.

Several dog skeletons were excavated at the Roman town of Augusta Raurica and at the military camp of Vindonissa, located in the northern Alpine region of Switzerland (Germania Superior). The relationships between them and the people, the nature of their lives, and the circumstances of their deaths are unclear. In order to gain insight into this dog population, we collected 31 dogs deposited almost simultaneously in two wells (second half of the third century CE), three dogs from burial contexts (70-200 CE and third to fifth century CE) at Augusta Raurica, and two dogs from burial contexts at Vindonissa (ca. first century CE). We detected a mixed population of young and adult dogs including small, medium and large sized individuals. Three small dogs had conspicuous phenotypes: abnormally short legs, and one with a brachycephalic skull. Stable isotope analysis of a subset of the dogs showed that their diets were omnivorous with a substantial input of animal proteins and little variation, except one with a particularly low δ15N value, indicating a diet low in animal proteins. Partial mitochondrial DNA sequences from 25 dogs revealed eight haplotypes within canine haplogroup A (11 dogs; 44%; 5 haplotypes), C (8 dogs; 32%; 1 haplotype), D (4 dogs, 16%; 1 haplotype) and B (2 dogs, 8%; 1 haplotype). Based on shotgun sequencing, four Roman mitogenomes were assembled, representing sub-haplogroups A1b3, A1b2 and C2. No canine pathogens were identified, weakening the assumption of infectious disease as a cause for dog disposal. The genetic and morphological diversity observed in dogs of Augusta Raurica and Vindonissa is similar to modern dog diversity.

Pathogen genomics study of an early medieval community in Germany reveals extensive co-infections.

BACKGROUND: The pathogen landscape in the Early European Middle Ages remains largely unexplored. Here, we perform a systematic pathogen screening of the rural community Lauchheim "Mittelhofen," in present-day Germany, dated to the Merovingian period, between fifth and eighth century CE. Skeletal remains of individuals were subjected to an ancient DNA metagenomic analysis. Genomes of the detected pathogens were reconstructed and analyzed phylogenetically. RESULTS: Over 30% of the individuals exhibit molecular signs of infection with hepatitis B virus (HBV), parvovirus B19, variola virus (VARV), and Mycobacterium leprae. Seven double and one triple infection were detected. We reconstructed four HBV genomes and one genome each of B19, VARV, and M. leprae. All HBV genomes are of genotype D4 which is rare in Europe today. The VARV strain exhibits a unique pattern of gene loss indicating that viruses with different gene compositions were circulating in the Early Middle Ages. The M. leprae strain clustered in branch 3 together with the oldest to-date genome from the UK. CONCLUSIONS: The high burden of infectious disease, together with osteological markers of physiological stress, reflect a poor health status of the community. This could have been an indirect result of the climate decline in Europe at the time, caused by the Late Antique Little Ice Age (LALIA). Our findings suggest that LALIA may have created an ecological context in which persistent outbreaks set the stage for major epidemics of severe diseases such as leprosy and smallpox hundreds of years later.

Phylogeography in an "oyster" shell provides first insights into the genetic structure of an extinct Ostrea edulis population.

The historical phylogeography of Ostrea edulis was successfully depicted in its native range for the first time using ancient DNA methods on dry shells from museum collections. This research reconstructed the historical population structure of the European flat oyster across Europe in the 1870s-including the now extinct population in the Wadden Sea. In total, four haplogroups were identified with one haplogroup having a patchy distribution from the North Sea to the Atlantic coast of France. This irregular distribution could be the result of translocations. The other three haplogroups are restricted to narrow geographic ranges, which may indicate adaptation to local environmental conditions or geographical barriers to gene flow. The phylogenetic reconstruction of the four haplogroups suggests the signatures of glacial refugia and postglacial expansion. The comparison with present-day O. edulis populations revealed a temporally stable population genetic pattern over the past 150 years despite large-scale translocations. This historical phylogeographic reconstruction was able to discover an autochthonous population in the German and Danish Wadden Sea in the late nineteenth century, where O. edulis is extinct today. The genetic distinctiveness of a now-extinct population hints at a connection between the genetic background of O. edulis in the Wadden Sea and for its absence until today.

Mass burial genomics reveals outbreak of enteric paratyphoid fever in the Late Medieval trade city Lübeck.

Medieval Europe was repeatedly affected by outbreaks of infectious diseases, some of which reached epidemic proportions. A Late Medieval mass burial next to the Heiligen-Geist-Hospital in Lübeck (present-day Germany) contained the skeletal remains of more than 800 individuals who had presumably died from infectious disease. From 92 individuals, we screened the ancient DNA extracts for the presence of pathogens to determine the cause of death. Metagenomic analysis revealed evidence of Salmonella enterica subsp. enterica serovar Paratyphi C, suggesting an outbreak of enteric paratyphoid fever. Three reconstructed S. Paratyphi C genomes showed close similarity to a strain from Norway (1200 CE). Radiocarbon dates placed the disease outbreak in Lübeck between 1270 and 1400 cal CE, with historical records indicating 1367 CE as the most probable year. The deceased were of northern and eastern European descent, confirming Lübeck as an important trading center of the Hanseatic League in the Baltic region.

A 5,000-year-old hunter-gatherer already plagued by Yersinia pestis.

A 5,000-year-old Yersinia pestis genome (RV 2039) is reconstructed from a hunter-fisher-gatherer (5300-5050 cal BP) buried at Rinnukalns, Latvia. RV 2039 is the first in a series of ancient strains that evolved shortly after the split of Y. pestis from its antecessor Y. pseudotuberculosis ∼7,000 years ago. The genomic and phylogenetic characteristics of RV 2039 are consistent with the hypothesis that this very early Y. pestis form was most likely less transmissible and maybe even less virulent than later strains. Our data do not support the scenario of a prehistoric pneumonic plague pandemic, as suggested previously for the Neolithic decline. The geographical and temporal distribution of the few prehistoric Y. pestis cases reported so far is more in agreement with single zoonotic events.

Genome-wide study of a Neolithic Wartberg grave community reveals distinct HLA variation and hunter-gatherer ancestry.

The Wartberg culture (WBC, 3500-2800 BCE) dates to the Late Neolithic period, a time of important demographic and cultural transformations in western Europe. We performed genome-wide analyses of 42 individuals who were interred in a WBC collective burial in Niedertiefenbach, Germany (3300-3200 cal. BCE). The results showed that the farming population of Niedertiefenbach carried a surprisingly large hunter-gatherer ancestry component (34-58%). This component was most likely introduced during the cultural transformation that led to the WBC. In addition, the Niedertiefenbach individuals exhibited a distinct human leukocyte antigen gene pool, possibly reflecting an immune response that was geared towards detecting viral infections.

Yersinia pestis strains from Latvia show depletion of the pla virulence gene at the end of the second plague pandemic.

Ancient genomic studies have identified Yersinia pestis (Y. pestis) as the causative agent of the second plague pandemic (fourteenth-eighteenth century) that started with the Black Death (1,347-1,353). Most of the Y. pestis strains investigated from this pandemic have been isolated from western Europe, and not much is known about the diversity and microevolution of this bacterium in eastern European countries. In this study, we investigated human remains excavated from two cemeteries in Riga (Latvia). Historical evidence suggests that the burials were a consequence of plague outbreaks during the seventeenth century. DNA was extracted from teeth of 16 individuals and subjected to shotgun sequencing. Analysis of the metagenomic data revealed the presence of Y. pestis sequences in four remains, confirming that the buried individuals were victims of plague. In two samples, Y. pestis DNA coverage was sufficient for genome reconstruction. Subsequent phylogenetic analysis showed that the Riga strains fell within the diversity of the already known post-Black Death genomes. Interestingly, the two Latvian isolates did not cluster together. Moreover, we detected a drop in coverage of the pPCP1 plasmid region containing the pla gene. Further analysis indicated the presence of two pPCP1 plasmids, one with and one without the pla gene region, and only one bacterial chromosome, indicating that the same bacterium carried two distinct pPCP1 plasmids. In addition, we found the same pattern in the majority of previously published post-Black Death strains, but not in the Black Death strains. The pla gene is an important virulence factor for the infection of and transmission in humans. Thus, the spread of pla-depleted strains may, among other causes, have contributed to the disappearance of the second plague pandemic in eighteenth century Europe.

Gene-flow from steppe individuals into Cucuteni-Trypillia associated populations indicates long-standing contacts and gradual admixture.

The Cucuteni-Trypillia complex (CTC) flourished in eastern Europe for over two millennia (5100-2800 BCE) from the end of the Neolithic to the Early Bronze Age. Its vast distribution area encompassed modern-day eastern Romania, Moldova and western/central Ukraine. Due to a lack of existing burials throughout most of this time, only little is known about the people associated with this complex and their genetic composition. Here, we present genome-wide data generated from the skeletal remains of four females that were excavated from two Late CTC sites in Moldova (3500-3100 BCE). All individuals carried a large Neolithic-derived ancestry component and were genetically more closely related to Linear Pottery than to Anatolian farmers. Three of the specimens also showed considerable amounts of steppe-related ancestry, suggesting influx into the CTC gene-pool from people affiliated with, for instance, the Ukraine Mesolithic. The latter scenario is supported by archaeological evidence. Taken together, our results confirm that the steppe component arrived in eastern Europe farming communities maybe as early as 3500 BCE. In addition, they are in agreement with the hypothesis of ongoing contacts and gradual admixture between incoming steppe and local western populations.

Targeted analysis of polymorphic loci from low-coverage shotgun sequence data allows accurate genotyping of HLA genes in historical human populations.

The highly polymorphic human leukocyte antigen (HLA) plays a crucial role in adaptive immunity and is associated with various complex diseases. Accurate analysis of HLA genes using ancient DNA (aDNA) data is crucial for understanding their role in human adaptation to pathogens. Here, we describe the TARGT pipeline for targeted analysis of polymorphic loci from low-coverage shotgun sequence data. The pipeline was successfully applied to medieval aDNA samples and validated using both simulated aDNA and modern empirical sequence data from the 1000 Genomes Project. Thus the TARGT pipeline enables accurate analysis of HLA polymorphisms in historical (and modern) human populations.

The Third Annual Meeting of the European Virus Bioinformatics Center.

The Third Annual Meeting of the European Virus Bioinformatics Center (EVBC) took place in Glasgow, United Kingdom, 28-29 March 2019. Virus bioinformatics has become central to virology research, and advances in bioinformatics have led to improved approaches to investigate viral infections and outbreaks, being successfully used to detect, control, and treat infections of humans and animals. This active field of research has attracted approximately 110 experts in virology and bioinformatics/computational biology from Europe and other parts of the world to attend the two-day meeting in Glasgow to increase scientific exchange between laboratory- and computer-based researchers. The meeting was held at the McIntyre Building of the University of Glasgow; a perfect location, as it was originally built to be a place for "rubbing your brains with those of other people", as Rector Stanley Baldwin described it. The goal of the meeting was to provide a meaningful and interactive scientific environment to promote discussion and collaboration and to inspire and suggest new research directions and questions. The meeting featured eight invited and twelve contributed talks, on the four main topics: (1) systems virology, (2) virus-host interactions and the virome, (3) virus classification and evolution and (4) epidemiology, surveillance and evolution. Further, the meeting featured 34 oral poster presentations, all of which focused on specific areas of virus bioinformatics. This report summarizes the main research findings and highlights presented at the meeting.

Neolithic and medieval virus genomes reveal complex evolution of hepatitis B.

The hepatitis B virus (HBV) is one of the most widespread human pathogens known today, yet its origin and evolutionary history are still unclear and controversial. Here, we report the analysis of three ancient HBV genomes recovered from human skeletons found at three different archaeological sites in Germany. We reconstructed two Neolithic and one medieval HBV genome by de novo assembly from shotgun DNA sequencing data. Additionally, we observed HBV-specific peptides using paleo-proteomics. Our results demonstrated that HBV has circulated in the European population for at least 7000 years. The Neolithic HBV genomes show a high genomic similarity to each other. In a phylogenetic network, they do not group with any human-associated HBV genome and are most closely related to those infecting African non-human primates. The ancient viruses appear to represent distinct lineages that have no close relatives today and possibly went extinct. Our results reveal the great potential of ancient DNA from human skeletons in order to study the long-time evolution of blood borne viruses.

Ancient DNA study reveals HLA susceptibility locus for leprosy in medieval Europeans.

Leprosy, a chronic infectious disease caused by Mycobacterium leprae (M. leprae), was very common in Europe till the 16th century. Here, we perform an ancient DNA study on medieval skeletons from Denmark that show lesions specific for lepromatous leprosy (LL). First, we test the remains for M. leprae DNA to confirm the infection status of the individuals and to assess the bacterial diversity. We assemble 10 complete M. leprae genomes that all differ from each other. Second, we evaluate whether the human leukocyte antigen allele DRB1*15:01, a strong LL susceptibility factor in modern populations, also predisposed medieval Europeans to the disease. The comparison of genotype data from 69 M. leprae DNA-positive LL cases with those from contemporary and medieval controls reveals a statistically significant association in both instances. In addition, we observe that DRB1*15:01 co-occurs with DQB1*06:02 on a haplotype that is a strong risk factor for inflammatory diseases today.