The population genomic legacy of the second plague pandemic

Human populations have been shaped by catastrophes that may have left long-lasting signatures in their genomes. One notable example is the second plague pandemic that entered Europe in ca. 1,347 CE and repeatedly returned for over 300 years, with typical village and town mortality estimated at 10%-40%.¹ It is assumed that this high mortality affected the gene pools of these populations. First, local population crashes reduced genetic diversity. Second, a change in frequency is expected for sequence variants that may have affected survival or susceptibility to the etiologic agent (Yersinia pestis).² Third, mass mortality might alter the local gene pools through its impact on subsequent migration patterns. We explored these factors using the Norwegian city of Trondheim as a model, by sequencing 54 genomes spanning three time periods: (1) prior to the plague striking Trondheim in 1,349 CE, (2) the 17^th-19^th century, and (3) the present. We find that the pandemic period shaped the gene pool by reducing long distance immigration, in particular from the British Isles, and inducing a bottleneck that reduced genetic diversity. Although we also observe an excess of large F_ST values at multiple loci in the genome, these are shaped by reference biases introduced by mapping our relatively low genome coverage degraded DNA to the reference genome. This implies that attempts to detect selection using ancient DNA (aDNA) datasets that vary by read length and depth of sequencing coverage may be particularly challenging until methods have been developed to account for the impact of differential reference bias on test statistics.

The influence of temperature on the seasonality of historical plague outbreaks

Modern plague outbreaks exhibit a distinct seasonal pattern. By contrast, the seasonality of historical outbreaks and its drivers has not been studied systematically. Here, we investigate the seasonal pattern, the epidemic peak timing and growth rates, and the association with latitude, temperature, and precipitation using a large, novel dataset of plague- and all-cause mortality during the Second Pandemic in Europe and the Mediterranean. We show that epidemic peak timing followed a latitudinal gradient, with mean annual temperature negatively associated with peak timing. Based on modern temperature data, the predicted epidemic growth of all outbreaks was positive between 11.7°C and 21.5°C with a maximum around 17.3°C. Hence, our study provides evidence that the growth of plague epidemics across the whole study region depended on similar absolute temperature thresholds. Here, we present a systematic analysis of the seasonality of historical plague in the Northern Hemisphere, and we show consistent evidence for a temperature-related process influencing the epidemic peak timing and growth rates of plague epidemics.