Cases of trisomy 21 and trisomy 18 among historic and prehistoric individuals discovered from ancient DNA.

Aneuploidies, and in particular, trisomies represent the most common genetic aberrations observed in human genetics today. To explore the presence of trisomies in historic and prehistoric populations we screen nearly 10,000 ancient human individuals for the presence of three copies of any of the target autosomes. We find clear genetic evidence for six cases of trisomy 21 (Down syndrome) and one case of trisomy 18 (Edwards syndrome), and all cases are present in infant or perinatal burials. We perform comparative osteological examinations of the skeletal remains and find overlapping skeletal markers, many of which are consistent with these syndromes. Interestingly, three cases of trisomy 21, and the case of trisomy 18 were detected in two contemporaneous sites in early Iron Age Spain (800-400 BCE), potentially suggesting a higher frequency of burials of trisomy carriers in those societies. Notably, the care with which the burials were conducted, and the items found with these individuals indicate that ancient societies likely acknowledged these individuals with trisomy 18 and 21 as members of their communities, from the perspective of burial practice.

Polygenic risk scores as a marker for epilepsy risk across lifetime and after unspecified seizure events.

A diagnosis of epilepsy has significant consequences for an individual but is often challenging in clinical practice. Novel biomarkers are thus greatly needed. Here, we investigated how common genetic factors (epilepsy polygenic risk scores, [PRSs]) influence epilepsy risk in detailed longitudinal electronic health records (EHRs) of > 360k Finns spanning up to 50 years of individuals' lifetimes. Individuals with a high genetic generalized epilepsy PRS (PRSGGE) in FinnGen had an increased risk for genetic generalized epilepsy (GGE) (hazard ratio [HR] 1.55 per PRSGGE standard deviation [SD]) across their lifetime and after unspecified seizure events. Effect sizes of epilepsy PRSs were comparable to effect sizes in clinically curated data supporting our EHR-derived epilepsy diagnoses. Within 10 years after an unspecified seizure, the GGE rate was 37% when PRSGGE > 2 SD compared to 5.6% when PRSGGE < -2 SD. The effect of PRSGGE was even larger on GGE subtypes of idiopathic generalized epilepsy (IGE) (HR 2.1 per SD PRSGGE). We further report significantly larger effects of PRSGGE on epilepsy in females and in younger age groups. Analogously, we found significant but more modest focal epilepsy PRS burden associated with non-acquired focal epilepsy (NAFE). We found PRSGGE specifically associated with GGE in comparison with >2000 independent diseases while PRSNAFE was also associated with other diseases than NAFE such as back pain. Here, we show that epilepsy specific PRSs have good discriminative ability after a first seizure event i.e. in circumstances where the prior probability of epilepsy is high outlining a potential to serve as biomarkers for an epilepsy diagnosis.

Epilepsies of presumed genetic etiology show enrichment of rare variants that occur in the general population.

Previous studies suggested that severe epilepsies, e.g., developmental and epileptic encephalopathies (DEEs), are mainly caused by ultra-rare de novo genetic variants. For milder disease, rare genetic variants could contribute to the phenotype. To determine the importance of rare variants for different epilepsy types, we analyzed a whole-exome sequencing cohort of 9,170 epilepsy-affected individuals and 8,436 control individuals. Here, we separately analyzed three different groups of epilepsies: severe DEEs, genetic generalized epilepsy (GGE), and non-acquired focal epilepsy (NAFE). We required qualifying rare variants (QRVs) to occur in control individuals with an allele count ≥ 1 and a minor allele frequency ≤ 1:1,000, to be predicted as deleterious (CADD ≥ 20), and to have an odds ratio in individuals with epilepsy ≥ 2. We identified genes enriched with QRVs primarily in NAFE (n = 72), followed by GGE (n = 32) and DEE (n = 21). This suggests that rare variants may play a more important role for causality of NAFE than for DEE. Moreover, we found that genes harboring QRVs, e.g., HSGP2, FLNA, or TNC, encode proteins that are involved in structuring the brain extracellular matrix. The present study confirms an involvement of rare variants for NAFE that occur also in the general population, while in DEE and GGE, the contribution of such variants appears more limited.

Implementation and evaluation of personal genetic testing as part of genomics analysis courses in German universities.

PURPOSE: Due to the increasing application of genome analysis and interpretation in medical disciplines, professionals require adequate education. Here, we present the implementation of personal genotyping as an educational tool in two genomics courses targeting Digital Health students at the Hasso Plattner Institute (HPI) and medical students at the Technical University of Munich (TUM). METHODS: We compared and evaluated the courses and the students' perceptions on the course setup using questionnaires. RESULTS: During the course, students changed their attitudes towards genotyping (HPI: 79% [15 of 19], TUM: 47% [25 of 53]). Predominantly, students became more critical of personal genotyping (HPI: 73% [11 of 15], TUM: 72% [18 of 25]) and most students stated that genetic analyses should not be allowed without genetic counseling (HPI: 79% [15 of 19], TUM: 70% [37 of 53]). Students found the personal genotyping component useful (HPI: 89% [17 of 19], TUM: 92% [49 of 53]) and recommended its inclusion in future courses (HPI: 95% [18 of 19], TUM: 98% [52 of 53]). CONCLUSION: Students perceived the personal genotyping component as valuable in the described genomics courses. The implementation described here can serve as an example for future courses in Europe.

The role of rare genetic variants enrichment in epilepsies of presumed genetic etiology.

Previous studies suggested that severe epilepsies e.g., developmental and epileptic encephalopathies (DEE) are mainly caused by ultra-rare de novo genetic variants. For milder phenotypes, rare genetic variants could contribute to the phenotype. To determine the importance of rare variants for different epilepsy types, we analyzed a whole-exome sequencing cohort of 9,170 epilepsy-affected individuals and 8,436 controls. Here, we separately analyzed three different groups of epilepsies : severe DEEs, genetic generalized epilepsy (GGE), and non-acquired focal epilepsy (NAFE). We required qualifying rare variants (QRVs) to occur in controls at a minor allele frequency ≤ 1:1,000, to be predicted as deleterious (CADD≥20), and to have an odds ratio in epilepsy cases ≥2. We identified genes enriched with QRVs in DEE (n=21), NAFE (n=72), and GGE (n=32) - the number of enriched genes are found greatest in NAFE and least in DEE. This suggests that rare variants may play a more important role for causality of NAFE than in DEE. Moreover, we found that QRV-carrying genes e.g., HSGP2, FLNA or TNC are involved in structuring the brain extracellular matrix. The present study confirms an involvement of rare variants for NAFE, while in DEE and GGE, the contribution of such variants appears more limited.

FinnGen provides genetic insights from a well-phenotyped isolated population.

Population isolates such as those in Finland benefit genetic research because deleterious alleles are often concentrated on a small number of low-frequency variants (0.1% ≤ minor allele frequency < 5%). These variants survived the founding bottleneck rather than being distributed over a large number of ultrarare variants. Although this effect is well established in Mendelian genetics, its value in common disease genetics is less explored1,2. FinnGen aims to study the genome and national health register data of 500,000 Finnish individuals. Given the relatively high median age of participants (63 years) and the substantial fraction of hospital-based recruitment, FinnGen is enriched for disease end points. Here we analyse data from 224,737 participants from FinnGen and study 15 diseases that have previously been investigated in large genome-wide association studies (GWASs). We also include meta-analyses of biobank data from Estonia and the United Kingdom. We identified 30 new associations, primarily low-frequency variants, enriched in the Finnish population. A GWAS of 1,932 diseases also identified 2,733 genome-wide significant associations (893 phenome-wide significant (PWS), P < 2.6 × 10-11) at 2,496 (771 PWS) independent loci with 807 (247 PWS) end points. Among these, fine-mapping implicated 148 (73 PWS) coding variants associated with 83 (42 PWS) end points. Moreover, 91 (47 PWS) had an allele frequency of <5% in non-Finnish European individuals, of which 62 (32 PWS) were enriched by more than twofold in Finland. These findings demonstrate the power of bottlenecked populations to find entry points into the biology of common diseases through low-frequency, high impact variants.

Rare coding variants in ten genes confer substantial risk for schizophrenia.

Rare coding variation has historically provided the most direct connections between gene function and disease pathogenesis. By meta-analysing the whole exomes of 24,248 schizophrenia cases and 97,322 controls, we implicate ultra-rare coding variants (URVs) in 10 genes as conferring substantial risk for schizophrenia (odds ratios of 3-50, P < 2.14 × 10-6) and 32 genes at a false discovery rate of <5%. These genes have the greatest expression in central nervous system neurons and have diverse molecular functions that include the formation, structure and function of the synapse. The associations of the NMDA (N-methyl-D-aspartate) receptor subunit GRIN2A and AMPA (α-amino-3-hydroxy-5-methyl-4-isoxazole propionic acid) receptor subunit GRIA3 provide support for dysfunction of the glutamatergic system as a mechanistic hypothesis in the pathogenesis of schizophrenia. We observe an overlap of rare variant risk among schizophrenia, autism spectrum disorders1, epilepsy and severe neurodevelopmental disorders2, although different mutation types are implicated in some shared genes. Most genes described here, however, are not implicated in neurodevelopment. We demonstrate that genes prioritized from common variant analyses of schizophrenia are enriched in rare variant risk3, suggesting that common and rare genetic risk factors converge at least partially on the same underlying pathogenic biological processes. Even after excluding significantly associated genes, schizophrenia cases still carry a substantial excess of URVs, which indicates that more risk genes await discovery using this approach.

The origin and legacy of the Etruscans through a 2000-year archeogenomic time transect.

The origin, development, and legacy of the enigmatic Etruscan civilization from the central region of the Italian peninsula known as Etruria have been debated for centuries. Here we report a genomic time transect of 82 individuals spanning almost two millennia (800 BCE to 1000 CE) across Etruria and southern Italy. During the Iron Age, we detect a component of Indo-European­associated steppe ancestry and the lack of recent Anatolian-related admixture among the putative non­Indo-European­speaking Etruscans. Despite comprising diverse individuals of central European, northern African, and Near Eastern ancestry, the local gene pool is largely maintained across the first millennium BCE. This drastically changes during the Roman Imperial period where we report an abrupt population-wide shift to ~50% admixture with eastern Mediterranean ancestry. Last, we identify northern European components appearing in central Italy during the Early Middle Ages, which thus formed the genetic landscape of present-day Italian populations.

A novel compound heterozygous leptin receptor mutation causes more severe obesity than in Leprdb/db mice.

The leptin receptor (Lepr) pathway is important for food intake regulation, energy expenditure, and body weight. Mutations in leptin and the Lepr have been shown to cause early-onset severe obesity in mice and humans. In studies with C57BL/6NCrl mice, we found a mouse with extreme obesity. To identify a putative spontaneous new form of monogenic obesity, we performed backcross studies with this mouse followed by a quantitative trait locus (QTL) analysis and sequencing of the selected chromosomal QTL region. We thereby identified a novel Lepr mutation (C57BL/6N-LeprL536Hfs*6-1NKB), which is located at chromosome 4, exon 11 within the CRH2-leptin-binding site. Compared with C57BL/6N mice, LeprL536Hfs*6 develop early onset obesity and their body weight exceeds that of Leprdb/db mice at an age of 30 weeks. Similar to Leprdb/db mice, the LeprL536Hfs*6 model is characterized by hyperphagia, obesity, lower energy expenditure and activity, hyperglycemia, and hyperinsulinemia compared with C57BL/6N mice. Crossing Leprdb/wt with LeprL536Hfs*6/wt mice results in compound heterozygous LeprL536Hfs*6/db mice, which develop even higher body weight and fat mass than both homozygous Leprdb/db and LeprL536Hfs*6 mice. Compound heterozygous Lepr deficiency affecting functionally different regions of the Lepr causes more severe obesity than the parental homozygous mutations.

KCNT1-related epilepsies and epileptic encephalopathies: phenotypic and mutational spectrum.

Variants in KCNT1, encoding a sodium-gated potassium channel (subfamily T member 1), have been associated with a spectrum of epilepsies and neurodevelopmental disorders. These range from familial autosomal dominant or sporadic sleep-related hypermotor epilepsy to epilepsy of infancy with migrating focal seizures (EIMFS) and include developmental and epileptic encephalopathies. This study aims to provide a comprehensive overview of the phenotypic and genotypic spectrum of KCNT1 mutation-related epileptic disorders in 248 individuals, including 66 previously unpublished and 182 published cases, the largest cohort reported so far. Four phenotypic groups emerged from our analysis: (i) EIMFS (152 individuals, 33 previously unpublished); (ii) developmental and epileptic encephalopathies other than EIMFS (non-EIMFS developmental and epileptic encephalopathies) (37 individuals, 17 unpublished); (iii) autosomal dominant or sporadic sleep-related hypermotor epilepsy (53 patients, 14 unpublished); and (iv) other phenotypes (six individuals, two unpublished). In our cohort of 66 new cases, the most common phenotypic features were: (i) in EIMFS, heterogeneity of seizure types, including epileptic spasms, epilepsy improvement over time, no epilepsy-related deaths; (ii) in non-EIMFS developmental and epileptic encephalopathies, possible onset with West syndrome, occurrence of atypical absences, possible evolution to developmental and epileptic encephalopathies with sleep-related hypermotor epilepsy features; one case of sudden unexplained death in epilepsy; (iii) in autosomal dominant or sporadic sleep-related hypermotor epilepsy, we observed a high prevalence of drug-resistance, although seizure frequency improved with age in some individuals, appearance of cognitive regression after seizure onset in all patients, no reported severe psychiatric disorders, although behavioural/psychiatric comorbidities were reported in ∼50% of the patients, sudden unexplained death in epilepsy in one individual; and (iv) other phenotypes in individuals with mutation of KCNT1 included temporal lobe epilepsy, and epilepsy with tonic-clonic seizures and cognitive regression. Genotypic analysis of the whole cohort of 248 individuals showed only missense mutations and one inframe deletion in KCNT1. Although the KCNT1 mutations in affected individuals were seen to be distributed among the different domains of the KCNT1 protein, genotype-phenotype considerations showed many of the autosomal dominant or sporadic sleep-related hypermotor epilepsy-associated mutations to be clustered around the RCK2 domain in the C terminus, distal to the NADP domain. Mutations associated with EIMFS/non-EIMFS developmental and epileptic encephalopathies did not show a particular pattern of distribution in the KCNT1 protein. Recurrent KCNT1 mutations were seen to be associated with both severe and less severe phenotypes. Our study further defines and broadens the phenotypic and genotypic spectrums of KCNT1-related epileptic conditions and emphasizes the increasingly important role of this gene in the pathogenesis of early onset developmental and epileptic encephalopathies as well as of focal epilepsies, namely autosomal dominant or sporadic sleep-related hypermotor epilepsy.

Comprehensive characterization of amino acid positions in protein structures reveals molecular effect of missense variants.

Interpretation of the colossal number of genetic variants identified from sequencing applications is one of the major bottlenecks in clinical genetics, with the inference of the effect of amino acid-substituting missense variations on protein structure and function being especially challenging. Here we characterize the three-dimensional (3D) amino acid positions affected in pathogenic and population variants from 1,330 disease-associated genes using over 14,000 experimentally solved human protein structures. By measuring the statistical burden of variations (i.e., point mutations) from all genes on 40 3D protein features, accounting for the structural, chemical, and functional context of the variations' positions, we identify features that are generally associated with pathogenic and population missense variants. We then perform the same amino acid-level analysis individually for 24 protein functional classes, which reveals unique characteristics of the positions of the altered amino acids: We observe up to 46% divergence of the class-specific features from the general characteristics obtained by the analysis on all genes, which is consistent with the structural diversity of essential regions across different protein classes. We demonstrate that the function-specific 3D features of the variants match the readouts of mutagenesis experiments for BRCA1 and PTEN, and positively correlate with an independent set of clinically interpreted pathogenic and benign missense variants. Finally, we make our results available through a web server to foster accessibility and downstream research. Our findings represent a crucial step toward translational genetics, from highlighting the impact of mutations on protein structure to rationalizing the variants' pathogenicity in terms of the perturbed molecular mechanisms.

Predicting functional effects of missense variants in voltage-gated sodium and calcium channels.

Malfunctions of voltage-gated sodium and calcium channels (encoded by SCNxA and CACNA1x family genes, respectively) have been associated with severe neurologic, psychiatric, cardiac, and other diseases. Altered channel activity is frequently grouped into gain or loss of ion channel function (GOF or LOF, respectively) that often corresponds not only to clinical disease manifestations but also to differences in drug response. Experimental studies of channel function are therefore important, but laborious and usually focus only on a few variants at a time. On the basis of known gene-disease mechanisms of 19 different diseases, we inferred LOF (n = 518) and GOF (n = 309) likely pathogenic variants from the disease phenotypes of variant carriers. By training a machine learning model on sequence- and structure-based features, we predicted LOF or GOF effects [area under the receiver operating characteristics curve (ROC) = 0.85] of likely pathogenic missense variants. Our LOF versus GOF prediction corresponded to molecular LOF versus GOF effects for 87 functionally tested variants in SCN1/2/8A and CACNA1I (ROC = 0.73) and was validated in exome-wide data from 21,703 cases and 128,957 controls. We showed respective regional clustering of inferred LOF and GOF nucleotide variants across the alignment of the entire gene family, suggesting shared pathomechanisms in the SCNxA/CACNA1x family genes.

MISCAST: MIssense variant to protein StruCture Analysis web SuiTe.

Human genome sequencing efforts have greatly expanded, and a plethora of missense variants identified both in patients and in the general population is now publicly accessible. Interpretation of the molecular-level effect of missense variants, however, remains challenging and requires a particular investigation of amino acid substitutions in the context of protein structure and function. Answers to questions like 'Is a variant perturbing a site involved in key macromolecular interactions and/or cellular signaling?', or 'Is a variant changing an amino acid located at the protein core or part of a cluster of known pathogenic mutations in 3D?' are crucial. Motivated by these needs, we developed MISCAST (missense variant to protein structure analysis web suite; http://miscast.broadinstitute.org/). MISCAST is an interactive and user-friendly web server to visualize and analyze missense variants in protein sequence and structure space. Additionally, a comprehensive set of protein structural and functional features have been aggregated in MISCAST from multiple databases, and displayed on structures alongside the variants to provide users with the biological context of the variant location in an integrated platform. We further made the annotated data and protein structures readily downloadable from MISCAST to foster advanced offline analysis of missense variants by a wide biological community.

Identification of pathogenic variant enriched regions across genes and gene families.

Missense variant interpretation is challenging. Essential regions for protein function are conserved among gene-family members, and genetic variants within these regions are potentially more likely to confer risk to disease. Here, we generated 2871 gene-family protein sequence alignments involving 9990 genes and performed missense variant burden analyses to identify novel essential protein regions. We mapped 2,219,811 variants from the general population into these alignments and compared their distribution with 76,153 missense variants from patients. With this gene-family approach, we identified 465 regions enriched for patient variants spanning 41,463 amino acids in 1252 genes. As a comparison, by testing the same genes individually, we identified fewer patient variant enriched regions, involving only 2639 amino acids and 215 genes. Next, we selected de novo variants from 6753 patients with neurodevelopmental disorders and 1911 unaffected siblings and observed an 8.33-fold enrichment of patient variants in our identified regions (95% C.I. = 3.90-Inf, P-value = 2.72 × 10-11). By using the complete ClinVar variant set, we found that missense variants inside the identified regions are 106-fold more likely to be classified as pathogenic in comparison to benign classification (OR = 106.15, 95% C.I = 70.66-Inf, P-value < 2.2 × 10-16). All pathogenic variant enriched regions (PERs) identified are available online through "PER viewer," a user-friendly online platform for interactive data mining, visualization, and download. In summary, our gene-family burden analysis approach identified novel PERs in protein sequences. This annotation can empower variant interpretation.

Paternal-age-related de novo mutations and risk for five disorders.

There are established associations between advanced paternal age and offspring risk for psychiatric and developmental disorders. These are commonly attributed to genetic mutations, especially de novo single nucleotide variants (dnSNVs), that accumulate with increasing paternal age. However, the actual magnitude of risk from such mutations in the male germline is unknown. Quantifying this risk would clarify the clinical significance of delayed paternity. Using parent-child trio whole-exome-sequencing data, we estimate the relationship between paternal-age-related dnSNVs and risk for five disorders: autism spectrum disorder (ASD), congenital heart disease, neurodevelopmental disorders with epilepsy, intellectual disability and schizophrenia (SCZ). Using Danish registry data, we investigate whether epidemiologic associations between each disorder and older fatherhood are consistent with the estimated role of dnSNVs. We find that paternal-age-related dnSNVs confer a small amount of risk for these disorders. For ASD and SCZ, epidemiologic associations with delayed paternity reflect factors that may not increase with age.

Targeted gene sequencing in 6994 individuals with neurodevelopmental disorder with epilepsy.

PURPOSE: We aimed to gain insight into frequencies of genetic variants in genes implicated in neurodevelopmental disorder with epilepsy (NDD+E) by investigating large cohorts of patients in a diagnostic setting. METHODS: We analyzed variants in NDD+E using epilepsy gene panel sequencing performed between 2013 and 2017 by two large diagnostic companies. We compared variant frequencies in 6994 panels with another 8588 recently published panels as well as exome-wide de novo variants in 1942 individuals with NDD+E and 10,937 controls. RESULTS: Genes with highest frequencies of ultrarare variants in NDD+E comprised SCN1A, KCNQ2, SCN2A, CDKL5, SCN8A, and STXBP1, concordant with the two other epilepsy cohorts we investigated. In only 46% of the analyzed 262 dominant and X-linked panel genes ultrarare variants in patients were reported. Among genes with contradictory evidence of association with epilepsy, CACNB4, CLCN2, EFHC1, GABRD, MAGI2, and SRPX2 showed equal frequencies in cases and controls. CONCLUSION: We show that improvement of panel design increased diagnostic yield over time, but panels still display genes with low or no diagnostic yield. With our data, we hope to improve current diagnostic NDD+E panel design and provide a resource of ultrarare variants in individuals with NDD+E to the community.

GRIN2A-related disorders: genotype and functional consequence predict phenotype.

Alterations of the N-methyl-d-aspartate receptor (NMDAR) subunit GluN2A, encoded by GRIN2A, have been associated with a spectrum of neurodevelopmental disorders with prominent speech-related features, and epilepsy. We performed a comprehensive assessment of phenotypes with a standardized questionnaire in 92 previously unreported individuals with GRIN2A-related disorders. Applying the criteria of the American College of Medical Genetics and Genomics to all published variants yielded 156 additional cases with pathogenic or likely pathogenic variants in GRIN2A, resulting in a total of 248 individuals. The phenotypic spectrum ranged from normal or near-normal development with mild epilepsy and speech delay/apraxia to severe developmental and epileptic encephalopathy, often within the epilepsy-aphasia spectrum. We found that pathogenic missense variants in transmembrane and linker domains (misTMD+Linker) were associated with severe developmental phenotypes, whereas missense variants within amino terminal or ligand-binding domains (misATD+LBD) and null variants led to less severe developmental phenotypes, which we confirmed in a discovery (P = 10-6) as well as validation cohort (P = 0.0003). Other phenotypes such as MRI abnormalities and epilepsy types were also significantly different between the two groups. Notably, this was paralleled by electrophysiology data, where misTMD+Linker predominantly led to NMDAR gain-of-function, while misATD+LBD exclusively caused NMDAR loss-of-function. With respect to null variants, we show that Grin2a+/- cortical rat neurons also had reduced NMDAR function and there was no evidence of previously postulated compensatory overexpression of GluN2B. We demonstrate that null variants and misATD+LBD of GRIN2A do not only share the same clinical spectrum (i.e. milder phenotypes), but also result in similar electrophysiological consequences (loss-of-function) opposing those of misTMD+Linker (severe phenotypes; predominantly gain-of-function). This new pathomechanistic model may ultimately help in predicting phenotype severity as well as eligibility for potential precision medicine approaches in GRIN2A-related disorders.

De novo variants in neurodevelopmental disorders with epilepsy.

Epilepsy is a frequent feature of neurodevelopmental disorders (NDDs), but little is known about genetic differences between NDDs with and without epilepsy. We analyzed de novo variants (DNVs) in 6,753 parent-offspring trios ascertained to have different NDDs. In the subset of 1,942 individuals with NDDs with epilepsy, we identified 33 genes with a significant excess of DNVs, of which SNAP25 and GABRB2 had previously only limited evidence of disease association. Joint analysis of all individuals with NDDs also implicated CACNA1E as a novel disease-associated gene. Comparing NDDs with and without epilepsy, we found missense DNVs, DNVs in specific genes, age of recruitment, and severity of intellectual disability to be associated with epilepsy. We further demonstrate the extent to which our results affect current genetic testing as well as treatment, emphasizing the benefit of accurate genetic diagnosis in NDDs with epilepsy.

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.