New insights into the Tyrolean Iceman's origin and phenotype as inferred by whole-genome sequencing. Nat Commun. 2012;3:698. [PMID: 22426219 DOI: 10.1038/ncomms1701]

Diagnostic metagenomics: potential applications to bacterial, viral and parasitic infections

The term ‘shotgun metagenomics’ is applied to the direct sequencing of DNA extracted from a sample without culture or target-specific amplification or capture. In diagnostic metagenomics, this approach is applied to clinical samples in the hope of detecting and characterizing pathogens. Here, I provide a conceptual overview, before reviewing several recent promising proof-of-principle applications of metagenomics in virus discovery, analysis of outbreaks and detection of pathogens in contemporary and historical samples. I also evaluate future prospects for diagnostic metagenomics in the light of relentless improvements in sequencing technologies.

A paleogenomic perspective on evolution and gene function: new insights from ancient DNA

The publication of partial and complete paleogenomes within the last few years has reinvigorated research in ancient DNA. No longer limited to short fragments of mitochondrial DNA, inference of evolutionary processes through time can now be investigated from genome-wide data sampled as far back as 700,000 years. Tremendous insights have been made, in particular regarding the hominin lineage. With rare exception, however, a paleogenomic perspective has been mired by the quality and quantity of recoverable DNA. Though conceptually simple, extracting ancient DNA remains challenging, and sequencing ancient genomes to high coverage remains prohibitively expensive for most laboratories. Still, with improvements in DNA isolation and declining sequencing costs, the taxonomic and geographic purview of paleogenomics is expanding at a rapid pace. With improved capacity to screen large numbers of samples for those with high proportions of endogenous ancient DNA, paleogenomics is poised to become a key technology to better understand recent evolutionary events.

The impact of whole-genome sequencing on the reconstruction of human population history

Examining patterns of molecular genetic variation in both modern-day and ancient humans has proved to be a powerful approach to learn about our origins. Rapid advances in DNA sequencing technology have allowed us to characterize increasing amounts of genomic information. Although this clearly provides unprecedented power for inference, it also introduces more complexity into the way we use and interpret such data. Here, we review ongoing debates that have been influenced by improvements in our ability to sequence DNA and discuss some of the analytical challenges that need to be overcome in order to fully exploit the rich historical information that is contained in the entirety of the human genome.

Separating endogenous ancient DNA from modern day contamination in a Siberian Neandertal

One of the main impediments for obtaining DNA sequences from ancient human skeletons is the presence of contaminating modern human DNA molecules in many fossil samples and laboratory reagents. However, DNA fragments isolated from ancient specimens show a characteristic DNA damage pattern caused by miscoding lesions that differs from present day DNA sequences. Here, we develop a framework for evaluating the likelihood of a sequence originating from a model with postmortem degradation-summarized in a postmortem degradation score-which allows the identification of DNA fragments that are unlikely to originate from present day sources. We apply this approach to a contaminated Neandertal specimen from Okladnikov Cave in Siberia to isolate its endogenous DNA from modern human contaminants and show that the reconstructed mitochondrial genome sequence is more closely related to the variation of Western Neandertals than what was discernible from previous analyses. Our method opens up the potential for genomic analysis of contaminated fossil material.

Derived immune and ancestral pigmentation alleles in a 7,000-year-old Mesolithic European

Ancient genomic sequences have started to reveal the origin and the demographic impact of farmers from the Neolithic period spreading into Europe. The adoption of farming, stock breeding and sedentary societies during the Neolithic may have resulted in adaptive changes in genes associated with immunity and diet. However, the limited data available from earlier hunter-gatherers preclude an understanding of the selective processes associated with this crucial transition to agriculture in recent human evolution. Here we sequence an approximately 7,000-year-old Mesolithic skeleton discovered at the La Braña-Arintero site in León, Spain, to retrieve a complete pre-agricultural European human genome. Analysis of this genome in the context of other ancient samples suggests the existence of a common ancient genomic signature across western and central Eurasia from the Upper Paleolithic to the Mesolithic. The La Braña individual carries ancestral alleles in several skin pigmentation genes, suggesting that the light skin of modern Europeans was not yet ubiquitous in Mesolithic times. Moreover, we provide evidence that a significant number of derived, putatively adaptive variants associated with pathogen resistance in modern Europeans were already present in this hunter-gatherer.

Genome-wide nucleosome map and cytosine methylation levels of an ancient human genome

Epigenetic information is available from contemporary organisms, but is difficult to track back in evolutionary time. Here, we show that genome-wide epigenetic information can be gathered directly from next-generation sequence reads of DNA isolated from ancient remains. Using the genome sequence data generated from hair shafts of a 4000-yr-old Paleo-Eskimo belonging to the Saqqaq culture, we generate the first ancient nucleosome map coupled with a genome-wide survey of cytosine methylation levels. The validity of both nucleosome map and methylation levels were confirmed by the recovery of the expected signals at promoter regions, exon/intron boundaries, and CTCF sites. The top-scoring nucleosome calls revealed distinct DNA positioning biases, attesting to nucleotide-level accuracy. The ancient methylation levels exhibited high conservation over time, clustering closely with modern hair tissues. Using ancient methylation information, we estimated the age at death of the Saqqaq individual and illustrate how epigenetic information can be used to infer ancient gene expression. Similar epigenetic signatures were found in other fossil material, such as 110,000- to 130,000-yr-old bones, supporting the contention that ancient epigenomic information can be reconstructed from a deep past. Our findings lay the foundation for extracting epigenomic information from ancient samples, allowing shifts in epialleles to be tracked through evolutionary time, as well as providing an original window into modern epigenomics.

Neolithic mitochondrial haplogroup H genomes and the genetic origins of Europeans

Haplogroup H dominates present-day Western European mitochondrial DNA variability (>40%), yet was less common (~19%) among Early Neolithic farmers (~5450 BC) and virtually absent in Mesolithic hunter-gatherers. Here we investigate this major component of the maternal population history of modern Europeans and sequence 39 complete haplogroup H mitochondrial genomes from ancient human remains. We then compare this 'real-time' genetic data with cultural changes taking place between the Early Neolithic (~5450 BC) and Bronze Age (~2200 BC) in Central Europe. Our results reveal that the current diversity and distribution of haplogroup H were largely established by the Mid Neolithic (~4000 BC), but with substantial genetic contributions from subsequent pan-European cultures such as the Bell Beakers expanding out of Iberia in the Late Neolithic (~2800 BC). Dated haplogroup H genomes allow us to reconstruct the recent evolutionary history of haplogroup H and reveal a mutation rate 45% higher than current estimates for human mitochondria.

Genomics and the Evolution of Phenotypic Traits

Evolutionary genetics has entered an unprecedented era of discovery, catalyzed in large part by the development of technologies that provide information about genome sequence and function. An important benefit is the ability to move beyond a handful of model organisms in lab settings to identify the genetic basis for evolutionarily interesting traits in many organisms in natural settings. Other benefits are the abilities to identify causal mutations and validate their phenotypic consequences more readily and in many more species. Genomic technologies have reinvigorated interest in some of the most fundamental and persistent questions in evolutionary genetics, revealed previously unsuspected evolutionary phenomena, and opened the door to a wide range of new questions.

Pulling out the 1%: whole-genome capture for the targeted enrichment of ancient DNA sequencing libraries

Most ancient specimens contain very low levels of endogenous DNA, precluding the shotgun sequencing of many interesting samples because of cost. Ancient DNA (aDNA) libraries often contain <1% endogenous DNA, with the majority of sequencing capacity taken up by environmental DNA. Here we present a capture-based method for enriching the endogenous component of aDNA sequencing libraries. By using biotinylated RNA baits transcribed from genomic DNA libraries, we are able to capture DNA fragments from across the human genome. We demonstrate this method on libraries created from four Iron Age and Bronze Age human teeth from Bulgaria, as well as bone samples from seven Peruvian mummies and a Bronze Age hair sample from Denmark. Prior to capture, shotgun sequencing of these libraries yielded an average of 1.2% of reads mapping to the human genome (including duplicates). After capture, this fraction increased substantially, with up to 59% of reads mapped to human and enrichment ranging from 6- to 159-fold. Furthermore, we maintained coverage of the majority of regions sequenced in the precapture library. Intersection with the 1000 Genomes Project reference panel yielded an average of 50,723 SNPs (range 3,062-147,243) for the postcapture libraries sequenced with 1 million reads, compared with 13,280 SNPs (range 217-73,266) for the precapture libraries, increasing resolution in population genetic analyses. Our whole-genome capture approach makes it less costly to sequence aDNA from specimens containing very low levels of endogenous DNA, enabling the analysis of larger numbers of samples.

Molecular phylogeography of a human autosomal skin color locus under natural selection

Divergent natural selection caused by differences in solar exposure has resulted in distinctive variations in skin color between human populations. The derived light skin color allele of the SLC24A5 gene, A111T, predominates in populations of Western Eurasian ancestry. To gain insight into when and where this mutation arose, we defined common haplotypes in the genomic region around SLC24A5 across diverse human populations and deduced phylogenetic relationships between them. Virtually all chromosomes carrying the A111T allele share a single 78-kb haplotype that we call C11, indicating that all instances of this mutation in human populations share a common origin. The C11 haplotype was most likely created by a crossover between two haplotypes, followed by the A111T mutation. The two parental precursor haplotypes are found from East Asia to the Americas but are nearly absent in Africa. The distributions of C11 and its parental haplotypes make it most likely that these two last steps occurred between the Middle East and the Indian subcontinent, with the A111T mutation occurring after the split between the ancestors of Europeans and East Asians.

Ligation bias in illumina next-generation DNA libraries: implications for sequencing ancient genomes

Ancient DNA extracts consist of a mixture of endogenous molecules and contaminant DNA templates, often originating from environmental microbes. These two populations of templates exhibit different chemical characteristics, with the former showing depurination and cytosine deamination by-products, resulting from post-mortem DNA damage. Such chemical modifications can interfere with the molecular tools used for building second-generation DNA libraries, and limit our ability to fully characterize the true complexity of ancient DNA extracts. In this study, we first use fresh DNA extracts to demonstrate that library preparation based on adapter ligation at AT-overhangs are biased against DNA templates starting with thymine residues, contrarily to blunt-end adapter ligation. We observe the same bias on fresh DNA extracts sheared on Bioruptor, Covaris and nebulizers. This contradicts previous reports suggesting that this bias could originate from the methods used for shearing DNA. This also suggests that AT-overhang adapter ligation efficiency is affected in a sequence-dependent manner and results in an uneven representation of different genomic contexts. We then show how this bias could affect the base composition of ancient DNA libraries prepared following AT-overhang ligation, mainly by limiting the ability to ligate DNA templates starting with thymines and therefore deaminated cytosines. This results in particular nucleotide misincorporation damage patterns, deviating from the signature generally expected for authenticating ancient sequence data. Consequently, we show that models adequate for estimating post-mortem DNA damage levels must be robust to the molecular tools used for building ancient DNA libraries.

Biomolecular Archaeology

Ancient biomolecules including DNA, proteins, and lipids are often preserved in archaeological skeletons or artifacts such as potsherds from cooking vessels. Techniques for analyzing these molecules have improved dramatically in recent years, though challenges remain in ensuring that results are authentic and not confused by the presence of contaminating modern biomolecules. Ancient DNA (aDNA) can be used to identify the sex, kinship relationships, and population affinities of human skeletons, and also to detect the presence of disease-causing organisms such as the plague and tuberculosis bacteria. Stable isotope ratios in collagen and other skeletal proteins enable past diets to be studied, and similar work with lipids from potsherds have revealed that dairying in Europe began 2,000 years earlier than previously thought. Biomolecular archaeology has therefore developed into a mature discipline that is making a significant contribution to different aspects of our understanding of the human past.

Paleoproteomic study of the Iceman's brain tissue

The Tyrolean Iceman, a Copper-age ice mummy, is one of the best-studied human individuals. While the genome of the Iceman has largely been decoded, tissue-specific proteomes have not yet been investigated. We studied the proteome of two distinct brain samples using gel-based and liquid chromatography-mass spectrometry-based proteomics technologies together with a multiple-databases and -search algorithms-driven data-analysis approach. Thereby, we identified a total of 502 different proteins. Of these, 41 proteins are known to be highly abundant in brain tissue and 9 are even specifically expressed in the brain. Furthermore, we found 10 proteins related to blood and coagulation. An enrichment analysis revealed a significant accumulation of proteins related to stress response and wound healing. Together with atomic force microscope scans, indicating clustered blood cells, our data reopens former discussions about a possible injury of the Iceman's head near the site where the tissue samples have been extracted.

Complete mitochondrial genome sequence of a Middle Pleistocene cave bear reconstructed from ultrashort DNA fragments

Although an inverse relationship is expected in ancient DNA samples between the number of surviving DNA fragments and their length, ancient DNA sequencing libraries are strikingly deficient in molecules shorter than 40 bp. We find that a loss of short molecules can occur during DNA extraction and present an improved silica-based extraction protocol that enables their efficient retrieval. In combination with single-stranded DNA library preparation, this method enabled us to reconstruct the mitochondrial genome sequence from a Middle Pleistocene cave bear (Ursus deningeri) bone excavated at Sima de los Huesos in the Sierra de Atapuerca, Spain. Phylogenetic reconstructions indicate that the U. deningeri sequence forms an early diverging sister lineage to all Western European Late Pleistocene cave bears. Our results prove that authentic ancient DNA can be preserved for hundreds of thousand years outside of permafrost. Moreover, the techniques presented enable the retrieval of phylogenetically informative sequences from samples in which virtually all DNA is diminished to fragments shorter than 50 bp.

Low-pass DNA sequencing of 1200 Sardinians reconstructs European Y-chromosome phylogeny

Genetic variation within the male-specific portion of the Y chromosome (MSY) can clarify the origins of contemporary populations, but previous studies were hampered by partial genetic information. Population sequencing of 1204 Sardinian males identified 11,763 MSY single-nucleotide polymorphisms, 6751 of which have not previously been observed. We constructed a MSY phylogenetic tree containing all main haplogroups found in Europe, along with many Sardinian-specific lineage clusters within each haplogroup. The tree was calibrated with archaeological data from the initial expansion of the Sardinian population ~7700 years ago. The ages of nodes highlight different genetic strata in Sardinia and reveal the presumptive timing of coalescence with other human populations. We calculate a putative age for coalescence of ~180,000 to 200,000 years ago, which is consistent with previous mitochondrial DNA-based estimates.

Patterns of population differentiation and natural selection on the celiac disease background risk network

Celiac disease is a common small intestinal inflammatory condition induced by wheat gluten and related proteins from rye and barley. Left untreated, the clinical presentation of CD can include failure to thrive, malnutrition, and distension in juveniles. The disease can additionally lead to vitamin deficiencies, anemia, and osteoporosis. Therefore, CD potentially negatively affected fitness in past populations utilizing wheat, barley, and rye. Previous analyses of CD risk variants have uncovered evidence for positive selection on some of these loci. These studies also suggest the possibility that risk for common autoimmune conditions such as CD may be the result of positive selection on immune related loci in the genome to fight infection. Under this evolutionary scenario, disease phenotypes may be a trade-off from positive selection on immunity. If this hypothesis is generally true, we can expect to find a signal of natural selection when we survey across the network of loci known to influence CD risk. This study examines the non-HLA autosomal network of gene loci associated with CD risk in Europe. We reject the null hypothesis of neutrality on this network of CD risk loci. Additionally, we can localize evidence of selection in time and space by adding information from the genome of the Tyrolean Iceman. While we can show significant differentiation between continental regions across the CD network, the pattern of evidence is not consistent with primarily recent (Holocene) selection across this network in Europe. Further localization of ancient selection on this network may illuminate the ecological pressures acting on the immune system during this critically interesting phase of our evolution.

Emulsion PCR-coupled target enrichment: an effective fishing method for high-throughput sequencing of poorly preserved ancient DNA

Due to the difficulties in deep sequencing, high-throughput sequencing of ancient DNA has been limited to exceptionally well-preserved ancient materials. The primary factor is microbial attack popularly observed in the buried materials, and it causes drastic increase in relative ratio of microbial DNA in the extracted DNA. We present a unified strategy in which emulsion PCR is coupled with target enrichment followed by next-generation sequencing. The method made it possible to obtain efficiently non-duplicated reads mapped to target sequences of interest, and this can achieve deep and reliable sequencing of ancient DNA from typical materials, even though poorly preserved.

Scalable metagenomic taxonomy classification using a reference genome database

Motivation: Deep metagenomic sequencing of biological samples has the potential to recover otherwise difficult-to-detect microorganisms and accurately characterize biological samples with limited prior knowledge of sample contents. Existing metagenomic taxonomic classification algorithms, however, do not scale well to analyze large metagenomic datasets, and balancing classification accuracy with computational efficiency presents a fundamental challenge.

Results: A method is presented to shift computational costs to an off-line computation by creating a taxonomy/genome index that supports scalable metagenomic classification. Scalable performance is demonstrated on real and simulated data to show accurate classification in the presence of novel organisms on samples that include viruses, prokaryotes, fungi and protists. Taxonomic classification of the previously published 150 giga-base Tyrolean Iceman dataset was found to take <20 h on a single node 40 core large memory machine and provide new insights on the metagenomic contents of the sample.

Availability: Software was implemented in C++ and is freely available at http://sourceforge.net/projects/lmat

Contact:allen99@llnl.gov

Supplementary information:Supplementary data are available at Bioinformatics online.

High resolution mapping of Y haplogroup G in Tyrol (Austria)

The distribution of Y-chromosomal haplogroup G2a (G-P15) in present-day paternal lineages in Tyrol (Austria) was analyzed by applying a high-density regional sampling scheme that also covered remote mountain areas. There is evidence from ancient genetic data for a high frequency of Y-chromosomal haplogroup G in prehistoric populations of Central Europe, whilst nowadays levels well below 10% are routinely observed. A population sample comprising ∼3700 specimens was analyzed for Y-chromosomal variation by genotyping Y-SNPs and Y-STRs. The set of binary markers included nine SNPs specific for sub-lineages of haplogroup G. The frequency of haplogroup G in 2379 unrelated men born in Tyrol amounted to 11.3%. Nearly all of these Y chromosomes belonged to haplogroup G2a. The main sub-haplogroup within G2a was defined by the SNP L497 (G2a3b1c) and reached a population frequency of 8.6%. Although this average level is higher than reported for other countries the geographical distribution of haplogroup G-L497 showed a differentiated pattern with a clustered distribution within some alpine valleys, where maxima above 40% were found. Both, the estimation of coalescent times and a principle coordinates analysis based on RST values derived from Y-STR haplotypes from different sub-regions of Tyrol revealed evidence for an old settlement history associated with Y chromosomes belonging to haplogroup G in the Tyrolean Alps.

Ancestry of modern Europeans: contributions of ancient DNA

Understanding the peopling history of Europe is crucial to comprehend the origins of modern populations. Of course, the analysis of current genetic data offers several explanations about human migration patterns which occurred on this continent, but it fails to explain precisely the impact of each demographic event. In this context, direct access to the DNA of ancient specimens allows the overcoming of recent demographic phenomena, which probably highly modified the constitution of the current European gene pool. In recent years, several DNA studies have been successfully conducted from ancient human remains thanks to the improvement of molecular techniques. They have brought new fundamental information on the peopling of Europe and allowed us to refine our understanding of European prehistory. In this review, we will detail all the ancient DNA studies performed to date on ancient European DNA from the Middle Paleolithic to the beginning of the protohistoric period.

A call to order at the spirochaetal host-pathogen interface

As the Lyme disease spirochaete Borrelia burgdorferi shuttles back and forth between arthropod vector and vertebrate host, it encounters vastly different and hostile environments. Major mechanisms contributing to the success of this pathogen throughout this complex transmission cycle are phase and antigenic variation of abundant and serotype-defining surface lipoproteins. These peripherally membrane-anchored virulence factors mediate niche-specific interactions with vector/host factors and protect the spirochaete from the perils of the mammalian immune response. In this issue of Molecular Microbiology, Tilly, Bestor and Rosa redefine the roles of two lipoproteins, OspC and VlsE, during mammalian infection. Using a variety of promoter fusions in combination with a sensitive in vivo 'use it or lose it' gene complementation assay, the authors demonstrate that proper sequential expression of OspC followed by VlsE indeed matters. A previously suggested general functional redundancy between these and other lipoproteins is shown to be limited and dependent on an immunodeficient experimental setting that is arguably of diminished ecological relevance. These data reinforce the notion that OspC plays a unique role during initial infection while the antigenically variant VlsE proteins allow for persistence in the mammalian host.

Efficient moment-based inference of admixture parameters and sources of gene flow

The recent explosion in available genetic data has led to significant advances in understanding the demographic histories of and relationships among human populations. It is still a challenge, however, to infer reliable parameter values for complicated models involving many populations. Here, we present MixMapper, an efficient, interactive method for constructing phylogenetic trees including admixture events using single nucleotide polymorphism (SNP) genotype data. MixMapper implements a novel two-phase approach to admixture inference using moment statistics, first building an unadmixed scaffold tree and then adding admixed populations by solving systems of equations that express allele frequency divergences in terms of mixture parameters. Importantly, all features of the model, including topology, sources of gene flow, branch lengths, and mixture proportions, are optimized automatically from the data and include estimates of statistical uncertainty. MixMapper also uses a new method to express branch lengths in easily interpretable drift units. We apply MixMapper to recently published data for Human Genome Diversity Cell Line Panel individuals genotyped on a SNP array designed especially for use in population genetics studies, obtaining confident results for 30 populations, 20 of them admixed. Notably, we confirm a signal of ancient admixture in European populations-including previously undetected admixture in Sardinians and Basques-involving a proportion of 20-40% ancient northern Eurasian ancestry.

Phenotypes from ancient DNA: approaches, insights and prospects

The great majority of phenotypic characteristics are complex traits, complicating the identification of the genes underlying their expression. However, both methodological and theoretical progress in genome-wide association studies have resulted in a much better understanding of the underlying genetics of many phenotypic traits, including externally visible characteristics (EVCs) such as eye and hair color. Consequently, it has become possible to predict EVCs from human samples lacking phenotypic information. Predicting EVCs from genetic evidence is clearly appealing for forensic applications involving the personal identification of human remains. Now, a recent paper has reported the genetic determination of eye and hair color in samples up to 800 years old. The ability to predict EVCs from ancient human remains opens up promising perspectives for ancient DNA research, as this could allow studies to directly address archaeological and evolutionary questions related to the temporal and geographical origins of the genetic variants underlying phenotypes.

The geography of recent genetic ancestry across Europe

A genomic survey of recent genealogical relatedness reveals the close ties of kinship and the impact of events across the past 3,000 years of European history.

The recent genealogical history of human populations is a complex mosaic formed by individual migration, large-scale population movements, and other demographic events. Population genomics datasets can provide a window into this recent history, as rare traces of recent shared genetic ancestry are detectable due to long segments of shared genomic material. We make use of genomic data for 2,257 Europeans (in the Population Reference Sample [POPRES] dataset) to conduct one of the first surveys of recent genealogical ancestry over the past 3,000 years at a continental scale. We detected 1.9 million shared long genomic segments, and used the lengths of these to infer the distribution of shared ancestors across time and geography. We find that a pair of modern Europeans living in neighboring populations share around 2–12 genetic common ancestors from the last 1,500 years, and upwards of 100 genetic ancestors from the previous 1,000 years. These numbers drop off exponentially with geographic distance, but since these genetic ancestors are a tiny fraction of common genealogical ancestors, individuals from opposite ends of Europe are still expected to share millions of common genealogical ancestors over the last 1,000 years. There is also substantial regional variation in the number of shared genetic ancestors. For example, there are especially high numbers of common ancestors shared between many eastern populations that date roughly to the migration period (which includes the Slavic and Hunnic expansions into that region). Some of the lowest levels of common ancestry are seen in the Italian and Iberian peninsulas, which may indicate different effects of historical population expansions in these areas and/or more stably structured populations. Population genomic datasets have considerable power to uncover recent demographic history, and will allow a much fuller picture of the close genealogical kinship of individuals across the world.

Few of us know our family histories more than a few generations back. It is therefore easy to overlook the fact that we are all distant cousins, related to one another via a vast network of relationships. Here we use genome-wide data from European individuals to investigate these relationships over the past 3,000 years, by looking for long stretches of genome that are shared between pairs of individuals through their inheritance from common genetic ancestors. We quantify this ubiquitous recent common ancestry, showing for instance that even pairs of individuals from opposite ends of Europe share hundreds of genetic common ancestors over this time period. Despite this degree of commonality, there are also striking regional differences. Southeastern Europeans, for example, share large numbers of common ancestors that date roughly to the era of the Slavic and Hunnic expansions around 1,500 years ago, while most common ancestors that Italians share with other populations lived longer than 2,500 years ago. The study of long stretches of shared genetic material promises to uncover rich information about many aspects of recent population history.

No ancient DNA damage in Actinobacteria from the Neanderthal bone

Background

The Neanderthal genome was recently sequenced using DNA extracted from a 38,000-year-old fossil. At the start of the project, the fraction of mammalian and bacterial DNA in the sample was estimated to be <6% and 9%, respectively. Treatment with restriction enzymes prior to sequencing increased the relative proportion of mammalian DNA to 15%, but the large majority of sequences remain uncharacterized.

Principal Findings

Our taxonomic profiling of 3.95 Gb of Neanderthal DNA isolated from the Vindija Neanderthal Vi33.16 fossil showed that 90% of about 50,000 rRNA gene sequence reads were of bacterial origin, of which Actinobacteria accounted for more than 75%. Actinobacteria also represented more than 80% of the PCR-amplified 16S rRNA gene sequences from a cave sediment sample taken from the same G layer as the Neanderthal bone. However, phylogenetic analyses did not identify any sediment clones that were closely related to the bone-derived sequences. We analysed the patterns of nucleotide differences in the individual sequence reads compared to the assembled consensus sequences of the rRNA gene sequences. The typical ancient nucleotide substitution pattern with a majority of C to T changes indicative of DNA damage was observed for the Neanderthal rRNA gene sequences, but not for the Streptomyces-like rRNA gene sequences.

Conclusions/Significance

Our analyses suggest that the Actinobacteria, and especially members of the Streptomycetales, contribute the majority of sequences in the DNA extracted from the Neanderthal fossil Vi33.16. The bacterial DNA showed no signs of damage, and we hypothesize that it was derived from bacteria that have been enriched inside the bone. The bioinformatic approach used here paves the way for future studies of microbial compositions and patterns of DNA damage in bacteria from archaeological bones. Such studies can help identify targeted measures to increase the relative amount of endogenous DNA in the sample.

Immunity, atherosclerosis and cardiovascular disease

Atherosclerosis, the major cause of cardiovascular disease (CVD), is a chronic inflammatory condition with immune competent cells in lesions producing mainly pro-inflammatory cytokines. Dead cells and oxidized forms of low density lipoproteins (oxLDL) are abundant. The major direct cause of CVD appears to be rupture of atherosclerotic plaques. oxLDL has proinflammatory and immune-stimulatory properties, causes cell death at higher concentrations and contains inflammatory phospholipids with phosphorylcholine (PC) as an interesting epitope. Antibodies against PC (anti-PC) may be atheroprotective, one mechanism being anti-inflammatory. Bacteria and virus have been discussed, but it has been difficult to find direct evidence, and antibiotic trials have not been successful. Heat shock proteins could be one major target for atherogenic immune reactions. More direct causes of plaque rupture include pro-inflammatory cytokines, chemokines, and lipid mediators. To prove that inflammation is a cause of atherosclerosis and CVD, clinical studies with anti-inflammatory and/or immune-modulatory treatment are needed. The potential causes of immune reactions and inflammation in atherosclerosis and how inflammation can be targeted therapeutically to provide novel treatments for CVD are reviewed.

Updating the Y-chromosomal phylogenetic tree for forensic applications based on whole genome SNPs

The Y-chromosomal phylogenetic tree has a wide variety of important forensic applications and therefore it needs to be state-of-the-art. Nevertheless, since the last 'official' published tree many publications reported additional Y-chromosomal lineages and other phylogenetic topologies. Therefore, it is difficult for forensic scientists to interpret those reports and use an up-to-date tree and corresponding nomenclature in their daily work. Whole genome sequencing (WGS) data is useful to verify and optimise the current phylogenetic tree for haploid markers. The AMY-tree software is the first open access program which analyses WGS data for Y-chromosomal phylogenetic applications. Here, all published information is collected in a phylogenetic tree and the correctness of this tree is checked based on the first large analysis of 747 WGS samples with AMY-tree. The obtained result is one phylogenetic tree with all peer-reviewed reported Y-SNPs without the observed recurrent and ambiguous mutations. Nevertheless, the results showed that currently only the genomes of a limited set of Y-chromosomal (sub-)haplogroups is available and that many newly reported Y-SNPs based on WGS projects are false positives, even with high sequencing coverage methods. This study demonstrates the usefulness of AMY-tree in the process of checking the quality of the present Y-chromosomal tree and it accentuates the difficulties to enlarge this tree based on only WGS methods.

Oral pathologies of the Neolithic Iceman, c.3,300 BC

The famous Iceman 'Ötzi' (South Tyrol Museum of Archaeology, Bolzano, Italy), a Neolithic human ice mummy, offers a unique opportunity to study evolutionary aspects of oral disease. The aim of this study was to assess, for the very first time, his oral cavity, which surprisingly had never been examined systematically. Based on several computed tomography (CT) scans from 1991 onwards and on macroscopic investigation, only a few findings, such as a central maxillary diastema, heavy abrasions, and missing wisdom teeth, were known. We re-evaluated the latest CT scans from 2005 and found various oral pathologies. In line with the increase of tooth decay in the Neolithic - because of diet change in this historic transition phase - several carious lesions were found, one of which penetrated into the dental pulp. In accordance with the Iceman's troubled life, as several injuries on his body and his violent death attest, mechanical trauma of one of his upper front teeth is evident. Finally, the poor periodontal condition of the Iceman's dentition (e.g. loss of alveolar bone), indicative of periodontitis, was assessed. These oral pathological findings in the Iceman's dentition provide a unique glimpse into the evolutionary history of oral conditions.

Revising a personal genome by comparing and combining data from two different sequencing platforms

For the robust practice of genomic medicine, sequencing results must be compatible, regardless of the sequencing technologies and algorithms used. Presently, genome sequencing is still an imprecise science and is complicated by differences in the chemistry, coverage, alignment, and variant-calling algorithms. We identified ~3.33 million single nucleotide variants (SNVs) and ~3.62 million SNVs in the SJK genome using SOLiD and Illumina data, respectively. Approximately 3 million SNVs were concordant between the two platforms while 68,532 SNVs were discordant; 219,616 SNVs were SOLiD-specific and 516,080 SNVs were Illumina-specific (i.e., platform-specific). Concordant, discordant, and platform-specific SNVs were further analyzed and characterized. Overall, a large portion of heterozygous SNVs that were discordant with genotyping calls of single nucleotide polymorphism chips were highly confident. Approximately 70% of the platform-specific SNVs were located in regions containing repetitive sequences. Such platform-specificity may arise from differences between platforms, with regard to read length (36 bp and 72 bp vs. 50 bp), insert size (~100-300 bp vs. ~1-2 kb), sequencing chemistry (sequencing-by-synthesis using single nucleotides vs. ligation-based sequencing using oligomers), and sequencing quality. When data from the two platforms were merged for variant calling, the proportion of callable regions of the reference genome increased to 99.66%, which was 1.43% higher than the average callability of the two platforms, representing ~40 million bases. In this study, we compared the differences in sequencing results between two sequencing platforms. Approximately 90% of the SNVs were concordant between the two platforms, yet ~10% of the SNVs were either discordant or platform-specific, indicating that each platform had its own strengths and weaknesses. When data from the two platforms were merged, both the overall callability of the reference genome and the overall accuracy of the SNVs improved, demonstrating the likelihood that a re-sequenced genome can be revised using complementary data.

First insights into the metagenome of Egyptian mummies using next-generation sequencing

We applied, for the first time, next-generation sequencing (NGS) technology on Egyptian mummies. Seven NGS datasets obtained from five randomly selected Third Intermediate to Graeco-Roman Egyptian mummies (806 BC-124AD) and two unearthed pre-contact Bolivian lowland skeletons were generated and characterised. The datasets were contrasted to three recently published NGS datasets obtained from cold-climate regions, i.e. the Saqqaq, the Denisova hominid and the Alpine Iceman. Analysis was done using one million reads of each newly generated or published dataset. Blastn and megablast results were analysed using MEGAN software. Distinct NGS results were replicated by specific and sensitive polymerase chain reaction (PCR) protocols in ancient DNA dedicated laboratories. Here, we provide unambiguous identification of authentic DNA in Egyptian mummies. The NGS datasets showed variable contents of endogenous DNA harboured in tissues. Three of five mummies displayed a human DNA proportion comparable to the human read count of the Saqqaq permafrost-preserved specimen. Furthermore, a metagenomic signature unique to mummies was displayed. By applying a "bacterial fingerprint", discrimination among mummies and other remains from warm areas outside Egypt was possible. Due to the absence of an adequate environment monitoring, a bacterial bloom was identified when analysing different biopsies from the same mummies taken after a lapse of time of 1.5 years. Plant kingdom representation in all mummy datasets was unique and could be partially associated with their use in embalming materials. Finally, NGS data showed the presence of Plasmodium falciparum and Toxoplasma gondii DNA sequences, indicating malaria and toxoplasmosis in these mummies. We demonstrate that endogenous ancient DNA can be extracted from mummies and serve as a proper template for the NGS technique, thus, opening new pathways of investigation for future genome sequencing of ancient Egyptian individuals.

AMY-tree: an algorithm to use whole genome SNP calling for Y chromosomal phylogenetic applications

Background

Due to the rapid progress of next-generation sequencing (NGS) facilities, an explosion of human whole genome data will become available in the coming years. These data can be used to optimize and to increase the resolution of the phylogenetic Y chromosomal tree. Moreover, the exponential growth of known Y chromosomal lineages will require an automatic determination of the phylogenetic position of an individual based on whole genome SNP calling data and an up to date Y chromosomal tree.

Results

We present an automated approach, ‘AMY-tree’, which is able to determine the phylogenetic position of a Y chromosome using a whole genome SNP profile, independently from the NGS platform and SNP calling program, whereby mistakes in the SNP calling or phylogenetic Y chromosomal tree are taken into account. Moreover, AMY-tree indicates ambiguities within the present phylogenetic tree and points out new Y-SNPs which may be phylogenetically relevant. The AMY-tree software package was validated successfully on 118 whole genome SNP profiles of 109 males with different origins. Moreover, support was found for an unknown recurrent mutation, wrong reported mutation conversions and a large amount of new interesting Y-SNPs.

Conclusions

Therefore, AMY-tree is a useful tool to determine the Y lineage of a sample based on SNP calling, to identify Y-SNPs with yet unknown phylogenetic position and to optimize the Y chromosomal phylogenetic tree in the future. AMY-tree will not add lineages to the existing phylogenetic tree of the Y-chromosome but it is the first step to analyse whole genome SNP profiles in a phylogenetic framework.

The musculoskeletal abnormalities of the Similaun Iceman ("ÖTZI"): clues to chronic pain and possible treatments

BACKGROUND AND INTRODUCTION

In 1991, a deceased human male was found frozen in a glacier pool in the Italian Alps in north west Italy, and is now carefully preserved in the South Tyrol Museum of Archaeology, in Bolzano, Italy. The bodily tissues of the 5,300 year old male (colloquially referred to as the Iceman or Ötzi) were well preserved despite damage related to freezing, and glacial movement. Associated articles of well-preserved clothing, tools, weapons and other devices were also present and have been studied in detail. Clinical examination and imaging investigations have also shown that the Icemen had experienced possible illnesses in his lifetime and had identifiable areas of arthritis and musculoskeletal injury. This report includes some key observations on the musculoskeletal state of Ötzi and reference to the involvement of tattoo markings. Some aspects about the aetiology of his abnormalities and inflammatory arthritis are considered along with possible treatments that he might have employed.

METHODS AND RESULTS

We (WFK and MK) undertook a clinical musculoskeletal examination of the Iceman, details of which with available photographs and radiographic imaging pertaining to the musculoskeletal findings of the Iceman are reported here. The skin of the Iceman has numerous linear carbon tattoos, which are not of a decorative type. These have been presumed to possibly be "medicinal" tattoos administered for therapeutic reasons and may have been used in acupuncture-like treatment of pain. Spinal imaging identified areas of spinal damage and our observations have provided clues as to possible sites of spinal initiated pain and hence sites for administration of the "medicinal" tattoos. We observed body areas of the Iceman, in which imaging demonstrated arthritis and other forms of long-term musculoskeletal damage, but which do not have adjacent or corresponding "medicinal" tattoos. We contend that the back and leg "medicinal" tattoos correspond directly to sites of chronic right knee and right ankle pain, and left thoracolumbar pain. They also correspond to lower lumbar and sciatic referred radicular pain which may have a contributory cause related to the presence of a transitional lumbar 5 vertebra. Using recent published data (Keller et al. in Nature Commun 3:698, 2012. doi: 10.1038/ncomms1701 ) of the genome structure of the Iceman, we suggest some potential causes of the osteoarthritis or inflammatory joint injury may relate to presence of coronary heart disease (CHD) and Lyme disease (Borrelia burgdorferi) infection. We speculate on possible medical applications of natural products for self-medication.

CONCLUSIONS

These observations highlight several diagnostic features of musculoskeletal conditions in the Iceman with the possibility that tattoos may have been used for diagnosis or location of his painful states. The origins of his musculoskeletal conditions are unclear but there are indications that Lyme disease and CHD may have been factors. The associations or use of natural products may give insights into their applications at the time of the life of the Iceman.

Bona fide colour: DNA prediction of human eye and hair colour from ancient and contemporary skeletal remains

BACKGROUND

DNA analysis of ancient skeletal remains is invaluable in evolutionary biology for exploring the history of species, including humans. Contemporary human bones and teeth, however, are relevant in forensic DNA analyses that deal with the identification of perpetrators, missing persons, disaster victims or family relationships. They may also provide useful information towards unravelling controversies that surround famous historical individuals. Retrieving information about a deceased person's externally visible characteristics can be informative in both types of DNA analyses. Recently, we demonstrated that human eye and hair colour can be reliably predicted from DNA using the HIrisPlex system. Here we test the feasibility of the novel HIrisPlex system at establishing eye and hair colour of deceased individuals from skeletal remains of various post-mortem time ranges and storage conditions.

METHODS

Twenty-one teeth between 1 and approximately 800 years of age and 5 contemporary bones were subjected to DNA extraction using standard organic protocol followed by analysis using the HIrisPlex system.

RESULTS

Twenty-three out of 26 bone DNA extracts yielded the full 24 SNP HIrisPlex profile, therefore successfully allowing model-based eye and hair colour prediction. HIrisPlex analysis of a tooth from the Polish general Władysław Sikorski (1881 to 1943) revealed blue eye colour and blond hair colour, which was positively verified from reliable documentation. The partial profiles collected in the remaining three cases (two contemporary samples and a 14th century sample) were sufficient for eye colour prediction.

CONCLUSIONS

Overall, we demonstrate that the HIrisPlex system is suitable, sufficiently sensitive and robust to successfully predict eye and hair colour from ancient and contemporary skeletal remains. Our findings, therefore, highlight the HIrisPlex system as a promising tool in future routine forensic casework involving skeletal remains, including ancient DNA studies, for the prediction of eye and hair colour of deceased individuals.

The half-life of DNA in bone: measuring decay kinetics in 158 dated fossils

Claims of extreme survival of DNA have emphasized the need for reliable models of DNA degradation through time. By analysing mitochondrial DNA (mtDNA) from 158 radiocarbon-dated bones of the extinct New Zealand moa, we confirm empirically a long-hypothesized exponential decay relationship. The average DNA half-life within this geographically constrained fossil assemblage was estimated to be 521 years for a 242 bp mtDNA sequence, corresponding to a per nucleotide fragmentation rate (k) of 5.50 × 10(-6) per year. With an effective burial temperature of 13.1°C, the rate is almost 400 times slower than predicted from published kinetic data of in vitro DNA depurination at pH 5. Although best described by an exponential model (R(2) = 0.39), considerable sample-to-sample variance in DNA preservation could not be accounted for by geologic age. This variation likely derives from differences in taphonomy and bone diagenesis, which have confounded previous, less spatially constrained attempts to study DNA decay kinetics. Lastly, by calculating DNA fragmentation rates on Illumina HiSeq data, we show that nuclear DNA has degraded at least twice as fast as mtDNA. These results provide a baseline for predicting long-term DNA survival in bone.

Biodiversity and disease: a synthesis of ecological perspectives on Lyme disease transmission

Recent reviews have argued that disease control is among the ecosystem services yielded by biodiversity. Lyme disease (LD) is commonly cited as the best example of the 'diluting' effect of biodiversity on disease transmission, but many studies document the opposite relationship, showing that human LD risk can increase with forestation. Here, we unify these divergent perspectives and find strong evidence for a positive link between biodiversity and LD at broad spatial scales (urban to suburban to rural) and equivocal evidence for a negative link between biodiversity and LD at varying levels of biodiversity within forests. This finding suggests that, across zoonotic disease agents, the biodiversity-disease relationship is scale dependent and complex.

Ancient admixture in human history

Population mixture is an important process in biology. We present a suite of methods for learning about population mixtures, implemented in a software package called ADMIXTOOLS, that support formal tests for whether mixture occurred and make it possible to infer proportions and dates of mixture. We also describe the development of a new single nucleotide polymorphism (SNP) array consisting of 629,433 sites with clearly documented ascertainment that was specifically designed for population genetic analyses and that we genotyped in 934 individuals from 53 diverse populations. To illustrate the methods, we give a number of examples that provide new insights about the history of human admixture. The most striking finding is a clear signal of admixture into northern Europe, with one ancestral population related to present-day Basques and Sardinians and the other related to present-day populations of northeast Asia and the Americas. This likely reflects a history of admixture between Neolithic migrants and the indigenous Mesolithic population of Europe, consistent with recent analyses of ancient bones from Sweden and the sequencing of the genome of the Tyrolean "Iceman."

An overview of the genetic structure within the Italian population from genome-wide data

In spite of the common belief of Europe as reasonably homogeneous at genetic level, advances in high-throughput genotyping technology have resolved several gradients which define different geographical areas with good precision. When Northern and Southern European groups were considered separately, there were clear genetic distinctions. Intra-country genetic differences were also evident, especially in Finland and, to a lesser extent, within other European populations. Here, we present the first analysis using the 125,799 genome-wide Single Nucleotide Polymorphisms (SNPs) data of 1,014 Italians with wide geographical coverage. We showed by using Principal Component analysis and model-based individual ancestry analysis, that the current population of Sardinia can be clearly differentiated genetically from mainland Italy and Sicily, and that a certain degree of genetic differentiation is detectable within the current Italian peninsula population. Pair-wise F(ST) statistics Northern and Southern Italy amounts approximately to 0.001 between, and around 0.002 between Northern Italy and Utah residents with Northern and Western European ancestry (CEU). The Italian population also revealed a fine genetic substructure underscoring by the genomic inflation (Sardinia vs. Northern Italy = 3.040 and Northern Italy vs. CEU = 1.427), warning against confounding effects of hidden relatedness and population substructure in association studies.

From next-generation sequencing alignments to accurate comparison and validation of single-nucleotide variants: the pibase software

Scientists working with single-nucleotide variants (SNVs), inferred by next-generation sequencing software, often need further information regarding true variants, artifacts and sequence coverage gaps. In clinical diagnostics, e.g. SNVs must usually be validated by visual inspection or several independent SNV-callers. We here demonstrate that 0.5-60% of relevant SNVs might not be detected due to coverage gaps, or might be misidentified. Even low error rates can overwhelm the true biological signal, especially in clinical diagnostics, in research comparing healthy with affected cells, in archaeogenetic dating or in forensics. For these reasons, we have developed a package called pibase, which is applicable to diploid and haploid genome, exome or targeted enrichment data. pibase extracts details on nucleotides from alignment files at user-specified coordinates and identifies reproducible genotypes, if present. In test cases pibase identifies genotypes at 99.98% specificity, 10-fold better than other tools. pibase also provides pair-wise comparisons between healthy and affected cells using nucleotide signals (10-fold more accurately than a genotype-based approach, as we show in our case study of monozygotic twins). This comparison tool also solves the problem of detecting allelic imbalance within heterozygous SNVs in copy number variation loci, or in heterogeneous tumor sequences.

The genetic history of Europeans

The evolutionary history of modern humans is characterized by numerous migrations driven by environmental change, population pressures, and cultural innovations. In Europe, the events most widely considered to have had a major impact on patterns of genetic diversity are the initial colonization of the continent by anatomically modern humans (AMH), the last glacial maximum, and the Neolithic transition. For some decades it was assumed that the geographical structuring of genetic diversity within Europe was mainly the result of gene flow during and soon after the Neolithic transition, but recent advances in next-generation sequencing (NGS) technologies, computer simulation modeling, and ancient DNA (aDNA) analyses are challenging this simplistic view. Here we review the current knowledge on the evolutionary history of humans in Europe based on archaeological and genetic data.

Ancient migratory events in the Middle East: new clues from the Y-chromosome variation of modern Iranians

Knowledge of high resolution Y-chromosome haplogroup diversification within Iran provides important geographic context regarding the spread and compartmentalization of male lineages in the Middle East and southwestern Asia. At present, the Iranian population is characterized by an extraordinary mix of different ethnic groups speaking a variety of Indo-Iranian, Semitic and Turkic languages. Despite these features, only few studies have investigated the multiethnic components of the Iranian gene pool. In this survey 938 Iranian male DNAs belonging to 15 ethnic groups from 14 Iranian provinces were analyzed for 84 Y-chromosome biallelic markers and 10 STRs. The results show an autochthonous but non-homogeneous ancient background mainly composed by J2a sub-clades with different external contributions. The phylogeography of the main haplogroups allowed identifying post-glacial and Neolithic expansions toward western Eurasia but also recent movements towards the Iranian region from western Eurasia (R1b-L23), Central Asia (Q-M25), Asia Minor (J2a-M92) and southern Mesopotamia (J1-Page08). In spite of the presence of important geographic barriers (Zagros and Alborz mountain ranges, and the Dasht-e Kavir and Dash-e Lut deserts) which may have limited gene flow, AMOVA analysis revealed that language, in addition to geography, has played an important role in shaping the nowadays Iranian gene pool. Overall, this study provides a portrait of the Y-chromosomal variation in Iran, useful for depicting a more comprehensive history of the peoples of this area as well as for reconstructing ancient migration routes. In addition, our results evidence the important role of the Iranian plateau as source and recipient of gene flow between culturally and genetically distinct populations.

Distinguishing the co-ancestries of haplogroup G Y-chromosomes in the populations of Europe and the Caucasus

Haplogroup G, together with J2 clades, has been associated with the spread of agriculture, especially in the European context. However, interpretations based on simple haplogroup frequency clines do not recognize underlying patterns of genetic diversification. Although progress has been recently made in resolving the haplogroup G phylogeny, a comprehensive survey of the geographic distribution patterns of the significant sub-clades of this haplogroup has not been conducted yet. Here we present the haplogroup frequency distribution and STR variation of 16 informative G sub-clades by evaluating 1472 haplogroup G chromosomes belonging to 98 populations ranging from Europe to Pakistan. Although no basal G-M201* chromosomes were detected in our data set, the homeland of this haplogroup has been estimated to be somewhere nearby eastern Anatolia, Armenia or western Iran, the only areas characterized by the co-presence of deep basal branches as well as the occurrence of high sub-haplogroup diversity. The P303 SNP defines the most frequent and widespread G sub-haplogroup. However, its sub-clades have more localized distribution with the U1-defined branch largely restricted to Near/Middle Eastern and the Caucasus, whereas L497 lineages essentially occur in Europe where they likely originated. In contrast, the only U1 representative in Europe is the G-M527 lineage whose distribution pattern is consistent with regions of Greek colonization. No clinal patterns were detected suggesting that the distributions are rather indicative of isolation by distance and demographic complexities.