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.

Current allele distribution of the human longevity gene APOE in Europe can mainly be explained by ancient admixture.

Variation in apolipoprotein E (APOE) has been shown to have the strongest genetic effect on human longevity. The aim of this study was to unravel the evolutionary history of the three major APOE alleles in Europe by analysing ancient samples up to 12,000 years old. We detected significant allele frequency shifts between populations and over time. Our analyses indicated that selection led to large frequency differences between the earliest European populations (i.e., hunter-gatherers vs. first farmers), possibly due to changes in diet/lifestyle. In contrast, the allele distributions in populations from ~4000 BCE onward can mainly be explained by admixture, suggesting that it also played an important role in shaping current APOE variation. In any case, the resulting allele frequencies strongly influence the predisposition for longevity today, likely as a consequence of past adaptations and demographic processes.

Genome-wide association study of population-standardised cognitive performance phenotypes in a rural South African community.

Cognitive function is an indicator for global physical and mental health, and cognitive impairment has been associated with poorer life outcomes and earlier mortality. A standard cognition test, adapted to a rural-dwelling African community, and the Oxford Cognition Screen-Plus were used to capture cognitive performance as five continuous traits (total cognition score, verbal episodic memory, executive function, language, and visuospatial ability) for 2,246 adults in this population of South Africans. A novel common variant, rs73485231, reached genome-wide significance for association with episodic memory using data for ~14 million markers imputed from the H3Africa genotyping array data. Window-based replication of previously implicated variants and regions of interest support the discovery of African-specific associated variants despite the small population size and low allele frequency. This African genome-wide association study identifies suggestive associations with general cognition and domain-specific cognitive pathways and lays the groundwork for further genomic studies on cognition in Africa.

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.

Long-Lived Individuals Show a Lower Burden of Variants Predisposing to Age-Related Diseases and a Higher Polygenic Longevity Score.

Longevity is a complex phenotype influenced by both environmental and genetic factors. The genetic contribution is estimated at about 25%. Despite extensive research efforts, only a few longevity genes have been validated across populations. Long-lived individuals (LLI) reach extreme ages with a relative low prevalence of chronic disability and major age-related diseases (ARDs). We tested whether the protection from ARDs in LLI can partly be attributed to genetic factors by calculating polygenic risk scores (PRSs) for seven common late-life diseases (Alzheimer's disease (AD), atrial fibrillation (AF), coronary artery disease (CAD), colorectal cancer (CRC), ischemic stroke (ISS), Parkinson's disease (PD) and type 2 diabetes (T2D)). The examined sample comprised 1351 German LLI (≥94 years, including 643 centenarians) and 4680 German younger controls. For all ARD-PRSs tested, the LLI had significantly lower scores than the younger control individuals (areas under the curve (AUCs): ISS = 0.59, p = 2.84 × 10-35; AD = 0.59, p = 3.16 × 10-25; AF = 0.57, p = 1.07 × 10-16; CAD = 0.56, p = 1.88 × 10-12; CRC = 0.52, p = 5.85 × 10-3; PD = 0.52, p = 1.91 × 10-3; T2D = 0.51, p = 2.61 × 10-3). We combined the individual ARD-PRSs into a meta-PRS (AUC = 0.64, p = 6.45 × 10-15). We also generated two genome-wide polygenic scores for longevity, one with and one without the TOMM40/APOE/APOC1 gene region (AUC (incl. TOMM40/APOE/APOC1) = 0.56, p = 1.45 × 10-5, seven variants; AUC (excl. TOMM40/APOE/APOC1) = 0.55, p = 9.85 × 10-3, 10,361 variants). Furthermore, the inclusion of nine markers from the excluded region (not in LD with each other) plus the APOE haplotype into the model raised the AUC from 0.55 to 0.61. Thus, our results highlight the importance of TOMM40/APOE/APOC1 as a longevity hub.

Detailed stratified GWAS analysis for severe COVID-19 in four European populations.

Given the highly variable clinical phenotype of Coronavirus disease 2019 (COVID-19), a deeper analysis of the host genetic contribution to severe COVID-19 is important to improve our understanding of underlying disease mechanisms. Here, we describe an extended genome-wide association meta-analysis of a well-characterized cohort of 3255 COVID-19 patients with respiratory failure and 12 488 population controls from Italy, Spain, Norway and Germany/Austria, including stratified analyses based on age, sex and disease severity, as well as targeted analyses of chromosome Y haplotypes, the human leukocyte antigen region and the SARS-CoV-2 peptidome. By inversion imputation, we traced a reported association at 17q21.31 to a ~0.9-Mb inversion polymorphism that creates two highly differentiated haplotypes and characterized the potential effects of the inversion in detail. Our data, together with the 5th release of summary statistics from the COVID-19 Host Genetics Initiative including non-Caucasian individuals, also identified a new locus at 19q13.33, including NAPSA, a gene which is expressed primarily in alveolar cells responsible for gas exchange in the lung.

Corrigendum: Apolipoprotein E Genetic Variation and its Association With Cognitive Function in Rural-Dwelling Older South Africans.

[This corrects the article DOI: 10.3389/fgene.2021.689756.].

Apolipoprotein E Genetic Variation and Its Association With Cognitive Function in Rural-Dwelling Older South Africans.

Apolipoprotein E (APOE) 𝜀4 allele carrier status is well known for its association with an increased likelihood of developing Alzheimer's disease, but its independent role in cognitive function is unclear. APOE genetic variation is understudied in African populations; hence, this cross-sectional study in a rural South African community examined allele and genotype frequencies, and their associations with cognitive function. Cognitive function was assessed using two different screening methods to produce a total cognition score and four domain-specific cognition scores for verbal episodic memory, executive function, language, and visuospatial ability. Cognitive phenotype and APOE genotype data were used to determine whether APOE variation was significantly associated with cognitive function in this population. Observed allele frequencies for 1776 participants from the HAALSI study [age 40-80years (mean=56.19); 58.2% female] were 58.1% (𝜀3), 25.4% (𝜀4) and 16.5% (𝜀2). Allele distributions were similar to the African super population, but different from all non-African super populations from the 1,000 Genomes Project. The 𝜀3 homozygous genotype was most common (34.9%) and used as the base genotype for comparison in regression models. Four models were tested for each of the five cognitive phenotypes to explore association of APOE variation with cognitive function. In the first model assessing association with all genotypes for all individuals, marginally significant associations were observed for 𝜀2 homozygotes where executive function scored higher by ~0.5 standard deviations (p=0.037, SE=0.23), and for 𝜀3/𝜀4 heterozygotes where visuospatial ability scores were lower (p=0.046, SE=0.14). These did not survive correction for multiple testing. Regional African population differences were observed at the APOE locus. Marginally, significant associations between APOE genotype, and executive function and visuospatial ability indicate the need for larger studies to better examine these associations in African populations. Furthermore, longitudinal data could shed light on APOE genetic association with rate of change, or decline, in cognitive function.

Ancient DNA Study in Medieval Europeans Shows an Association Between HLA-DRB1*03 and Paratyphoid Fever.

Outbreaks of infectious diseases repeatedly affected medieval Europe, leaving behind a large number of dead often inhumed in mass graves. Human remains interred in two burial pits from 14th century CE Germany exhibited molecular evidence of Salmonella enterica Paratyphi C (S. Paratyphi C) infection. The pathogen is responsible for paratyphoid fever, which was likely the cause of death for the buried individuals. This finding presented the unique opportunity to conduct a paratyphoid fever association study in a European population. We focused on HLA-DRB1*03:01 that is a known risk allele for enteric fever in present-day South Asians. We generated HLA profiles for 29 medieval S. Paratyphi C cases and 24 contemporaneous controls and compared these to a modern German population. The frequency of the risk allele was higher in the medieval cases (29.6%) compared to the contemporaneous controls (13%; p = 0.189), albeit not significantly so, possibly because of small sample sizes. Indeed, in comparison with the modern controls (n = 39,689; 10.2%; p = 0.005) the frequency difference became statistically significant. This comparison also suggested a slight decrease in the allele's prevalence between the medieval and modern controls. Up to now, this is the first study on the genetic predisposition to Salmonella infection in Europeans and the only association analysis on paratyphoid fever C. Functional investigation using computational binding prediction between HLA variants and S. Paratyphi and S. Typhi peptides supported a reduced recognition capacity of bacterial proteins by DRB1*03:01 relative to other common DRB1 variants. This pattern could potentially explain the disease association. Our results suggest a slightly reduced predisposition to paratyphoid fever in modern Europeans. The causative allele, however, is still common today, which can be explained by a trade-off, as DRB1*03:01 is protective against infectious respiratory diseases such as severe respiratory syndrome (SARS). It is thus possible that the allele also provided resistance to corona-like viruses in the past.

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.

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.

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.

Exome-Wide Association Study Identifies FN3KRP and PGP as New Candidate Longevity Genes.

Despite enormous research efforts, the genetic component of longevity has remained largely elusive. The investigation of common variants, mainly located in intronic or regulatory regions, has yielded only little new information on the heritability of the phenotype. Here, we performed a chip-based exome-wide association study investigating 62 488 common and rare coding variants in 1248 German long-lived individuals, including 599 centenarians and 6941 younger controls (age < 60 years). In a single-variant analysis, we observed an exome-wide significant association between rs1046896 in the gene fructosamine-3-kinase-related-protein (FN3KRP) and longevity. Noteworthy, we found the longevity allele C of rs1046896 to be associated with an increased FN3KRP expression in whole blood; a database look-up confirmed this effect for various other human tissues. A gene-based analysis, in which potential cumulative effects of common and rare variants were considered, yielded the gene phosphoglycolate phosphatase (PGP) as another potential longevity gene, though no single variant in PGP reached the discovery p-value (1 × 10E-04). Furthermore, we validated the previously reported longevity locus cyclin-dependent kinase inhibitor 2B antisense RNA 1 (CDKN2B-AS1). Replication of our results in a French longevity cohort was only successful for rs1063192 in CDKN2B-AS1. In conclusion, we identified 2 new potential candidate longevity genes, FN3KRP and PGP which may influence the phenotype through their role in metabolic processes, that is, the reverse glycation of proteins (FN3KRP) and the control of glycerol-3-phosphate levels (PGP).

Rare genetic coding variants associated with human longevity and protection against age-related diseases.

Extreme longevity in humans has a strong genetic component, but whether this involves genetic variation in the same longevity pathways as found in model organisms is unclear. Using whole-exome sequences of a large cohort of Ashkenazi Jewish centenarians to examine enrichment for rare coding variants, we found most longevity-associated rare coding variants converge upon conserved insulin/insulin-like growth factor 1 signaling and AMP-activating protein kinase signaling pathways. Centenarians have a number of pathogenic rare coding variants similar to control individuals, suggesting that rare variants detected in the conserved longevity pathways are protective against age-related pathology. Indeed, we detected a pro-longevity effect of rare coding variants in the Wnt signaling pathway on individuals harboring the known common risk allele APOE4. The genetic component of extreme human longevity constitutes, at least in part, rare coding variants in pathways that protect against aging, including those that control longevity in model organisms.

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.

High mitochondrial diversity of domesticated goats persisted among Bronze and Iron Age pastoralists in the Inner Asian Mountain Corridor.

Goats were initially managed in the Near East approximately 10,000 years ago and spread across Eurasia as economically productive and environmentally resilient herd animals. While the geographic origins of domesticated goats (Capra hircus) in the Near East have been long-established in the zooarchaeological record and, more recently, further revealed in ancient genomes, the precise pathways by which goats spread across Asia during the early Bronze Age (ca. 3000 to 2500 cal BC) and later remain unclear. We analyzed sequences of hypervariable region 1 and cytochrome b gene in the mitochondrial genome (mtDNA) of goats from archaeological sites along two proposed transmission pathways as well as geographically intermediary sites. Unexpectedly high genetic diversity was present in the Inner Asian Mountain Corridor (IAMC), indicated by mtDNA haplotypes representing common A lineages and rarer C and D lineages. High mtDNA diversity was also present in central Kazakhstan, while only mtDNA haplotypes of lineage A were observed from sites in the Northern Eurasian Steppe (NES). These findings suggest that herding communities living in montane ecosystems were drawing from genetically diverse goat populations, likely sourced from communities in the Iranian Plateau, that were sustained by repeated interaction and exchange. Notably, the mitochondrial genetic diversity associated with goats of the IAMC also extended into the semi-arid region of central Kazakhstan, while NES communities had goats reflecting an isolated founder population, possibly sourced via eastern Europe or the Caucasus region.

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.

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.

DNA methylation QTL analysis identifies new regulators of human longevity.

Human longevity is a complex trait influenced by both genetic and environmental factors, whose interaction is mediated by epigenetic mechanisms like DNA methylation. Here, we generated genome-wide whole-blood methylome data from 267 individuals, of which 71 were long-lived (90-104 years), by applying reduced representation bisulfite sequencing. We followed a stringent two-stage analysis procedure using discovery and replication samples to detect differentially methylated sites (DMSs) between young and long-lived study participants. Additionally, we performed a DNA methylation quantitative trait loci analysis to identify DMSs that underlie the longevity phenotype. We combined the DMSs results with gene expression data as an indicator of functional relevance. This approach yielded 21 new candidate genes, the majority of which are involved in neurophysiological processes or cancer. Notably, two candidates (PVRL2, ERCC1) are located on chromosome 19q, in close proximity to the well-known longevity- and Alzheimer's disease-associated loci APOE and TOMM40. We propose this region as a longevity hub, operating on both a genetic (APOE, TOMM40) and an epigenetic (PVRL2, ERCC1) level. We hypothesize that the heritable methylation and associated gene expression changes reported here are overall advantageous for the LLI and may prevent/postpone age-related diseases and facilitate survival into very old age.

Quantitative comparison of within-sample heterogeneity scores for DNA methylation data.

DNA methylation is an epigenetic mark with important regulatory roles in cellular identity and can be quantified at base resolution using bisulfite sequencing. Most studies are limited to the average DNA methylation levels of individual CpGs and thus neglect heterogeneity within the profiled cell populations. To assess this within-sample heterogeneity (WSH) several window-based scores that quantify variability in DNA methylation in sequencing reads have been proposed. We performed the first systematic comparison of four published WSH scores based on simulated and publicly available datasets. Moreover, we propose two new scores and provide guidelines for selecting appropriate scores to address cell-type heterogeneity, cellular contamination and allele-specific methylation. Most of the measures were sensitive in detecting DNA methylation heterogeneity in these scenarios, while we detected differences in susceptibility to technical bias. Using recently published DNA methylation profiles of Ewing sarcoma samples, we show that DNA methylation heterogeneity provides information complementary to the DNA methylation level. WSH scores are powerful tools for estimating variance in DNA methylation patterns and have the potential for detecting novel disease-associated genomic loci not captured by established statistics. We provide an R-package implementing the WSH scores for integration into analysis workflows.