Genetic factors have been estimated to account for about 25% of the variation in an adult's life span. The complement component C4 with the isotypes C4A and C4B is an effector protein of the immune system, and differences in the overall C4 copy number or gene size (long C4L; short C4S) may influence the strength of the immune response and disease susceptibilities. Previously, an association between C4B copy number and life span was reported for Hungarians and Icelanders, where the C4B*Q0 genotype, which is defined by C4B gene deficiency, showed a decrease in frequency with age. Additionally, one of the studies indicated that a low C4B copy number might be a genetic trait that is manifested only in the presence of the environmental risk factor "smoking". These observations prompted us to investigate the role of the C4 alleles in our large German longevity sample (∼ 700 cases; 94-110 years and ∼ 900 younger controls). No significant differences in the number of C4A, C4B and C4S were detected. Besides, the C4B*Q0 carrier state did not decrease with age, irrespective of smoking as an interacting variable. However, for C4L*Q0 a significantly different carrier frequency was observed in the cases compared with controls (cases: 5.08%; controls: 9.12%; p = 0.003). In a replication sample of 714 German cases (91-108 years) and 890 controls this result was not replicated (p = 0.14) although a similar trend of decreased C4L*Q0 carrier frequency in cases was visible (cases: 7.84%; controls: 10.00%).

Mesolithic populations throughout Europe used diverse resource exploitation strategies that focused heavily on collecting and hunting wild prey. Between 5500 and 4200 cal BC, agriculturalists migrated into northwestern Europe bringing a suite of Neolithic technologies including domesticated animals. Here we investigate to what extent Mesolithic Ertebølle communities in northern Germany had access to domestic pigs, possibly through contact with neighbouring Neolithic agricultural groups. We employ a multidisciplinary approach, applying sequencing of ancient mitochondrial and nuclear DNA (coat colour-coding gene MC1R) as well as traditional and geometric morphometric (molar size and shape) analyses in Sus specimens from 17 Neolithic and Ertebølle sites. Our data from 63 ancient pig specimens show that Ertebølle hunter-gatherers acquired domestic pigs of varying size and coat colour that had both Near Eastern and European mitochondrial DNA ancestry. Our results also reveal that domestic pigs were present in the region ~500 years earlier than previously demonstrated.

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.

Studying ancient DNA allows us to retrace the evolutionary history of human pathogens, such as Mycobacterium leprae, the main causative agent of leprosy. Leprosy is one of the oldest recorded and most stigmatizing diseases in human history. The disease was prevalent in Europe until the 16th century and is still endemic in many countries with over 200,000 new cases reported annually. Previous worldwide studies on modern and European medieval M. leprae genomes revealed that they cluster into several distinct branches of which two were present in medieval Northwestern Europe. In this study, we analyzed 10 new medieval M. leprae genomes including the so far oldest M. leprae genome from one of the earliest known cases of leprosy in the United Kingdom-a skeleton from the Great Chesterford cemetery with a calibrated age of 415-545 C.E. This dataset provides a genetic time transect of M. leprae diversity in Europe over the past 1500 years. We find M. leprae strains from four distinct branches to be present in the Early Medieval Period, and strains from three different branches were detected within a single cemetery from the High Medieval Period. Altogether these findings suggest a higher genetic diversity of M. leprae strains in medieval Europe at various time points than previously assumed. The resulting more complex picture of the past phylogeography of leprosy in Europe impacts current phylogeographical models of M. leprae dissemination. It suggests alternative models for the past spread of leprosy such as a wide spread prevalence of strains from different branches in Eurasia already in Antiquity or maybe even an origin in Western Eurasia. Furthermore, these results highlight how studying ancient M. leprae strains improves understanding the history of leprosy worldwide.

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.

Human longevity is heritable, but genome-wide association (GWA) studies have had limited success. Here, we perform two meta-analyses of GWA studies of a rigorous longevity phenotype definition including 11,262/3484 cases surviving at or beyond the age corresponding to the 90th/99th survival percentile, respectively, and 25,483 controls whose age at death or at last contact was at or below the age corresponding to the 60th survival percentile. Consistent with previous reports, rs429358 (apolipoprotein E (ApoE) ε4) is associated with lower odds of surviving to the 90th and 99th percentile age, while rs7412 (ApoE ε2) shows the opposite. Moreover, rs7676745, located near GPR78, associates with lower odds of surviving to the 90th percentile age. Gene-level association analysis reveals a role for tissue-specific expression of multiple genes in longevity. Finally, genetic correlation of the longevity GWA results with that of several disease-related phenotypes points to a shared genetic architecture between health and 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.

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.

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).

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.

Apolipoprotein E (APOE) is a multifunctional plasma protein mainly acting in lipid metabolism. Human APOE is polymorphic with three major isoforms (APOE2, APOE3 and APOE4). Up to 75% of the body's APOE is produced by the liver. There is increasing evidence from studies in brain-derived cells that APOE4 affects mitochondrial function and biogenesis as well as stress and inflammatory responses - processes, whose disturbances are considered hallmarks of the ageing process. However, although the liver is the main production site of APOE, knowledge about the impact of the APOE genotype on hepatic stress response-related processes is rather limited. Therefore, we studied biomarkers of oxidative status (glutathione levels, 3-nitrotyrosine adducts, protein carbonyl concentration), ER stress (XBP1(S), BiP, DDIT3), proteasome activity, mitochondrial function (respiratory complexes, ATP levels and mitochondrial membrane potential as well as biomarkers of mitochondrial biogenesis, fission and fusion), autophagy (LC3, LAMP2A), apoptosis (BCL2, BAX, CYCS) and DNA damage in the liver of APOE targeted replacement (TR) mice and in Huh7 hepatocytes overexpressing the APOE3 and the APOE4 isoform, respectively. APOE4 mice exhibited a lower chymotrypsin-like and a higher trypsin-like proteasome activity. Levels of protein carbonyls were moderately higher in liver tissue of APOE4 vs. APOE3 mice. Other biomarkers of oxidative stress were similar between the two genotypes. Under basal conditions, the stress-response pathways investigated appeared largely unaffected by the APOE genotype. However, upon stress induction, APOE4 expressing cells showed lower levels of adenosine triphosphate (ATP) and lower mRNA levels of the ATP-generating complex V of the mitochondrial respiratory chain. Overall, our findings provide evidence for a rather low influence of the APOE genotype on the hepatic stress response processes investigated in this study.

Interferon gamma is a major macrophage-activating cytokine during infection with Mycobacterium tuberculosis, the causative pathogen of tuberculosis, and its role has been well established in animal models and in humans. This cytokine is produced by activated T helper 1 cells, which can best deal with intracellular pathogens such as M. tuberculosis. Based on the hypothesis that genes which regulate interferon gamma may influence tuberculosis susceptibility, we investigated polymorphisms in eight candidate genes.

The analysis of structural variants, in particular of copy-number variations (CNVs), has proven valuable in unraveling the genetic basis of human diseases. Hence, a large number of algorithms have been developed for the detection of CNVs in SNP array signal intensity data. Using the European and African HapMap trio data, we undertook a comparative evaluation of six commonly used CNV detection software tools, namely Affymetrix Power Tools (APT), QuantiSNP, PennCNV, GLAD, R-gada and VEGA, and assessed their level of pair-wise prediction concordance. The tool-specific CNV prediction accuracy was assessed in silico by way of intra-familial validation. Software tools differed greatly in terms of the number and length of the CNVs predicted as well as the number of markers included in a CNV. All software tools predicted substantially more deletions than duplications. Intra-familial validation revealed consistently low levels of prediction accuracy as measured by the proportion of validated CNVs (34-60%). Moreover, up to 20% of apparent family-based validations were found to be due to chance alone. Software using Hidden Markov models (HMM) showed a trend to predict fewer CNVs than segmentation-based algorithms albeit with greater validity. PennCNV yielded the highest prediction accuracy (60.9%). Finally, the pairwise concordance of CNV prediction was found to vary widely with the software tools involved. We recommend HMM-based software, in particular PennCNV, rather than segmentation-based algorithms when validity is the primary concern of CNV detection. QuantiSNP may be used as an additional tool to detect sets of CNVs not detectable by the other tools. Our study also reemphasizes the need for laboratory-based validation, such as qPCR, of CNVs predicted in silico.

Triosephosphate isomerase (TPI) is a central and conserved glycolytic enzyme. In humans, TPI is encoded by a single gene on 12p13, and associated with a rare genetic disorder, TPI deficiency. Reduced TPI activity can increase specific oxidant resistances of model organisms and TPI null-alleles have been hypothesized to promote a heterozygote advantage in man. However, comprehensive genetic information about the TPI1 locus is still lacking.

Human longevity is a complex phenotype influenced by genetic and environmental components. Unraveling the contribution of genetic vs. nongenetic factors to longevity is a challenging task. Here, we conducted a large-scale RNA-sequencing-based expression quantitative trait loci study (eQTL) with subsequent heritability analysis. The investigation was performed on blood samples from 244 individuals from Germany and Denmark, representing various age groups including long-lived subjects up to the age of 104 years. Our eQTL-based approach revealed for the first time that human longevity is associated with a depletion of metabolic pathways in a genotype-dependent and independent manner. Further analyses indicated that 20% of the differentially expressed genes are influenced by genetic variants in cis. The subsequent study of twins showed that the transcriptional activity of a third of the differentially regulated genes is heritable. These findings suggest that longevity-associated biological processes such as altered metabolism are, to a certain extent, also the driving force of longevity rather than just a consequence of old age.

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 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.

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.

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.

A 5,000-year-old Yersinia pestis genome (RV 2039) is reconstructed from a hunter-fisher-gatherer (5300-5050 cal BP) buried at Riņņukalns, 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.