Monitoring DNA contamination in handled vs. directly excavated ancient human skeletal remains

COL1A2 Barcoding: Bone Species Identification via Shotgun Proteomics

Species identification of fragmentary bones remains a challenging task in archeology and forensics. A species identification method for such fragmentary bones that has recently attracted interest is the use of bone collagen proteins. Here, we describe a method similar to DNA barcoding that reads collagen protein sequences in bone and automatically determines the species by performing sequence database searches. The method is almost identical to conventional shotgun proteomics analysis of bone samples, except that the database used by the SEQUEST search engine consisted only of entries for collagen type 1 alpha 2 (COL1A2) proteins from various vertebrates. Accordingly, the COL1A2 peptides that differ in sequence among species act as species marker peptides. In SEQUEST-based shotgun proteomics, the protein entries that contain more marker peptide sequences are assigned higher scores; therefore, the highest-scoring protein entry will be the COL1A2 entry for the species from which the analyzed bone was derived. We tested our method using bone samples from 30 vertebrate species and found that all species were correctly identified. In conclusion, COL1A2 can be used as a bone protein barcode and can be read through shotgun proteomics, allowing for automatic bone species identification. Data are available via ProteomeXchange with the identifier PXD045402.

A melting pot of Roman dogs north of the Alps with high phenotypic and genetic diversity and similar diets

Several dog skeletons were excavated at the Roman town of Augusta Raurica and at the military camp of Vindonissa, located in the northern Alpine region of Switzerland (Germania Superior). The relationships between them and the people, the nature of their lives, and the circumstances of their deaths are unclear. In order to gain insight into this dog population, we collected 31 dogs deposited almost simultaneously in two wells (second half of the third century CE), three dogs from burial contexts (70-200 CE and third to fifth century CE) at Augusta Raurica, and two dogs from burial contexts at Vindonissa (ca. first century CE). We detected a mixed population of young and adult dogs including small, medium and large sized individuals. Three small dogs had conspicuous phenotypes: abnormally short legs, and one with a brachycephalic skull. Stable isotope analysis of a subset of the dogs showed that their diets were omnivorous with a substantial input of animal proteins and little variation, except one with a particularly low δ15N value, indicating a diet low in animal proteins. Partial mitochondrial DNA sequences from 25 dogs revealed eight haplotypes within canine haplogroup A (11 dogs; 44%; 5 haplotypes), C (8 dogs; 32%; 1 haplotype), D (4 dogs, 16%; 1 haplotype) and B (2 dogs, 8%; 1 haplotype). Based on shotgun sequencing, four Roman mitogenomes were assembled, representing sub-haplogroups A1b3, A1b2 and C2. No canine pathogens were identified, weakening the assumption of infectious disease as a cause for dog disposal. The genetic and morphological diversity observed in dogs of Augusta Raurica and Vindonissa is similar to modern dog diversity.

Mitochondrial DNA control region typing from highly degraded skeletal remains by single-multiplex next-generation sequencing

Poor nuclear DNA preservation from highly degraded skeletal remains is the most limiting factor for the genetic identification of individuals. Mitochondrial DNA (mtDNA) typing, and especially of the control region (CR), using next-generation sequencing (NGS), enables retrieval of valuable genetic information in forensic contexts where highly degraded human skeletal remains are the only source of genetic material. Currently, NGS commercial kits can type all mtDNA-CR in fewer steps than the conventional Sanger technique. The PowerSeq CRM Nested System kit (Promega Corporation) employs a nested multiplex-polymerase chain reaction (PCR) strategy to amplify and index all mtDNA-CR in a single reaction. Our study analyzes the success of mtDNA-CR typing of highly degraded human skeletons using the PowerSeq CRM Nested System kit. We used samples from 41 individuals from different time periods to test three protocols (M1, M2, and M3) based on modifications of PCR conditions. To analyze the detected variants, two bioinformatic procedures were compared: an in-house pipeline and the GeneMarker HTS software. The results showed that many samples were not analyzed when the standard protocol (M1) was used. In contrast, the M3 protocol, which includes 35 PCR cycles and longer denaturation and extension steps, successfully recovered the mtDNA-CR from highly degraded skeletal samples. Mixed base profiles and the percentage of damaged reads were both indicators of possible contamination and can provide better results if used together. Furthermore, our freely available in-house pipeline can provide variants concordant with the forensic software.

Ancient Yersinia pestis genomes lack the virulence-associated YpfΦ prophage present in modern pandemic strains

Yersinia pestis is the causative agent of at least three major plague pandemics (Justinianic, Medieval and Modern). Previous studies on ancient Y. pestis genomes revealed that several genomic alterations had occurred approximately 5000-3000 years ago and contributed to the remarkable virulence of this pathogen. How a subset of strains evolved to cause the Modern pandemic is less well-understood. Here, we examined the virulence-associated prophage (YpfΦ), which had been postulated to be exclusively present in the genomes of strains associated with the Modern pandemic. The analysis of two new Y. pestis genomes from medieval/early modern Denmark confirmed that the phage is absent from the genome of strains dating to this time period. An extended comparative genome analysis of over 300 strains spanning more than 5000 years showed that the prophage is found in the genomes of modern strains only and suggests an integration into the genome during recent Y. pestis evolution. The phage-encoded Zot protein showed structural homology to a virulence factor of Vibrio cholerae. Similar to modern Y. pestis, we observed phages with a common origin to YpfΦ in individual strains of other bacterial species. Our findings present an updated view on the prevalence of YpfΦ, which might contribute to our understanding of the host spectrum, geographical spread and virulence of Y. pestis responsible for the Modern pandemic.

Estimation of DNA Degradation in Archaeological Human Remains

The evaluation of the integrity and quantity of DNA extracted from archaeological human remains is a fundamental step before using the latest generation sequencing techniques in the study of evolutionary processes. Ancient DNA is highly fragmented and chemically modified; therefore, the present study aims to identify indices that can allow the identification of potentially amplifiable and sequenceable DNA samples, reducing failures and research costs. Ancient DNA was extracted from five human bone remains from the archaeological site of Amiternum L'Aquila, Italy dating back to the 9th-12th century and was compared with standard DNA fragmented by sonication. Given the different degradation kinetics of mitochondrial DNA compared to nuclear DNA, the mitochondrially encoded 12s RNA and 18s ribosomal RNA genes were taken into consideration; fragments of various sizes were amplified in qPCR and the size distribution was thoroughly investigated. DNA damage degree was evaluated by calculating damage frequency (λ) and the ratio between the amount of the different fragments and that of the smallest fragment (Q). The results demonstrate that both indices were found to be suitable for identifying, among the samples tested, those less damaged and suitable for post-extraction analysis; mitochondrial DNA is more damaged than nuclear, in fact, amplicons up to 152 bp and 253 bp, respectively are obtained.

A mobile lab for ancient DNA extraction in Perug

We report the use of a mobile laboratory set up to extract ancient DNA (aDNA) from 34 human coprolites (fossilized faeces) samples. Our approach enabled the rapid genetic characterization of 5,000 years old archeological samples. It is useful for the on-site screening of museums and freshly excavated samples for DNA. This approach is accessible to other investigators as the mobile laboratory was set up using commercially available instruments.

A Case Study for the Recovery of Authentic Microbial Ancient DNA from Soil Samples

High Throughput DNA Sequencing (HTS) revolutionized the field of paleomicrobiology, leading to an explosive growth of microbial ancient DNA (aDNA) studies, especially from environmental samples. However, aDNA studies that examine environmental microbes routinely fail to authenticate aDNA, examine laboratory and environmental contamination, and control for biases introduced during sample processing. Here, we surveyed the available literature for environmental aDNA projects—from sample collection to data analysis—and assessed previous methodologies and approaches used in the published microbial aDNA studies. We then integrated these concepts into a case study, using shotgun metagenomics to examine methodological, technical, and analytical biases during an environmental aDNA study of soil microbes. Specifically, we compared the impact of five DNA extraction methods and eight bioinformatic pipelines on the recovery of microbial aDNA information in soil cores from extreme environments. Our results show that silica-based methods optimized for aDNA research recovered significantly more damaged and shorter reads (<100 bp) than a commercial kit or a phenol−chloroform method. Additionally, we described a stringent pipeline for data preprocessing, efficiently decreasing the representation of low-complexity and duplicated reads in our datasets and downstream analyses, reducing analytical biases in taxonomic classification.

A qPCR-duplex assay for sex determination in ancient DNA

Molecular biology techniques are increasingly being used in sex identification of skeletal remains when traditional anthropometric analyzes are not successful in identifying sex of remains that are incomplete, fragmented and /or of immature individuals. In the present work, we investigated the possibility of determining sex by using the qPCR-duplex method for both ancient and modern DNA samples. This method involves the co-amplification of two genes in a single reaction system and the subsequent analysis of the fusion curves; the gene sequences used for the construction of suitable primers are those of steroid sulfatase (STS) and testis specific protein Y-linked 1 (TSPY) genes which turned out to be two sensitive markers as they have a detection limit of 60 pg and 20 pg respectively on modern DNA. The validity of the method was verified on modern DNA in which gender was identified in all the samples with 100% accuracy; thus, allowing for the same results as the classic method with amelogenin, but in a faster and more immediate way, as it allows for sex determination solely by analyzing the denaturation curves without having to perform an electrophoretic run. The proposed molecular technique proves to be sensitive and precise even on degraded DNA, in fact on 9 archaeological finds dating from the VII-XII century in which sex had been identified through anthropometric analysis, it confirmed the sex of 8 out of 9 finds correctly.

Nanobioengineering: A promising approach for early detection of COVID-19

Unique pneumonia due to an unknown source emerged in December 2019 in the city of Wuhan, China. Consequently, the World Health Organization (WHO) declared this condition as a new coronavirus disease-19 also known as COVID-19 on February 11, 2020, which on March 13, 2020 was declared as a pandemic. The virus that causes COVID-19 was found to have a similar genome (80% similarity) with the previously known acute respiratory syndrome also known as SARS-CoV. The novel virus was later named Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). SARS-CoV-2 falls in the family of Coronaviridae which is further divided into Nidovirales and another subfamily called Orthocoronavirinae. The four generations of the coronaviruses belongs to the Orthocoronavirinae family that consists of alpha, beta, gamma and delta coronavirus which are denoted as α-CoV, β-CoV, γ-CoV, δ-CoV respectively. The α-CoV and β-CoVs are mainly known to infect mammals whereas γ-CoV and δ-CoV are generally found in birds. The β-CoVs also comprise of SARS-CoV and also include another virus that was found in the Middle East called the Middle East respiratory syndrome virus (MERS-CoV) and the cause of current pandemic SARS-CoV-2. These viruses initially cause the development of pneumonia in the patients and further development of a severe case of acute respiratory distress syndrome (ARDS) and other related symptoms that can be fatal leading to death.

Refining the evolutionary time machine: An assessment of whole genome amplification using single historical Daphnia eggs

Whole genome sequencing is instrumental for the study of genome variation in natural populations, delivering important knowledge on genomic modifications and potential targets of natural selection at the population level. Large dormant eggbanks of aquatic invertebrates such as the keystone herbivore Daphnia, a microcrustacean widespread in freshwater ecosystems, provide detailed sedimentary archives to study genomic processes over centuries. To overcome the problem of limited DNA amounts in single Daphnia dormant eggs, we developed an optimized workflow for whole genome amplification (WGA), yielding sufficient amounts of DNA for downstream whole genome sequencing of individual historical eggs, including polyploid lineages. We compare two WGA kits, applied to recently produced Daphnia magna dormant eggs from laboratory cultures, and to historical dormant eggs of Daphnia pulicaria collected from Arctic lake sediment between 10 and 300 years old. Resulting genome coverage breadth in most samples was ~70%, including those from >100-year-old isolates. Sequence read distribution was highly correlated among samples amplified with the same kit, but less correlated between kits. Despite this, a high percentage of genomic positions with single nucleotide polymorphisms in one or more samples (maximum of 74% between kits, and 97% within kits) were recovered at a depth required for genotyping. As a by-product of sequencing we obtained 100% coverage of the mitochondrial genomes even from the oldest isolates (~300 years). The mitochondrial DNA provides an additional source for evolutionary studies of these populations. We provide an optimized workflow for WGA followed by whole genome sequencing including steps to minimize exogenous DNA.

Insights into Punic genetic signatures in the southern necropolis of Tharros (Sardinia)

BACKGROUND

Phoenician and Punic expansions have been protagonists of intense trade networks and settlements in the Mediterranean Sea.

AIMS

The maternal genetic variability of ancient Punic samples from the Sardinian necropolis of Tharros was analysed, with the aim to explore genetic interactions and signatures of past population events.

SUBJECTS AND METHODS

The mtDNA HVS-I and coding region SNPs were analysed in 14 Punic samples and 74 modern individuals from Cabras and Belvì (for which the HVS-II region was also analysed). The results were compared with 5,590 modern Euro-Mediterranean sequences and 127 ancient samples.

RESULTS

While contemporary groups fall within the genetic variability of other modern Sardinians, our Punic samples reveal proximity to present-day North-African and Iberian populations. Furthermore, Cabras and Belvì cluster mainly with pre-Phoenician groups, while samples from Tharros project with other Punic Sardinian individuals.

CONCLUSION

This study provides the first preliminary insights into the population dynamics of the Punic site of Tharros. While the number of currently available samples does not allow definitive investigation of the connection with indigenous Sardinian groups, our results seem to confirm internal migratory phenomena in the central-western Mediterranean and female participation in the Punic mobility.

Present-Day DNA Contamination in Ancient DNA Datasets

Present-day contamination can lead to false conclusions in ancient DNA studies. A number of methods are available to estimate contamination, which use a variety of signals and are appropriate for different types of data. Here an overview of currently available methods highlighting their strengths and weaknesses is provided, and a classification based on the signals used to estimate contamination is proposed. This overview aims at enabling researchers to choose the most appropriate methods for their dataset. Based on this classification, potential avenues for the further development of methods are discussed.

Targeted analysis of polymorphic loci from low-coverage shotgun sequence data allows accurate genotyping of HLA genes in historical human populations

The highly polymorphic human leukocyte antigen (HLA) plays a crucial role in adaptive immunity and is associated with various complex diseases. Accurate analysis of HLA genes using ancient DNA (aDNA) data is crucial for understanding their role in human adaptation to pathogens. Here, we describe the TARGT pipeline for targeted analysis of polymorphic loci from low-coverage shotgun sequence data. The pipeline was successfully applied to medieval aDNA samples and validated using both simulated aDNA and modern empirical sequence data from the 1000 Genomes Project. Thus the TARGT pipeline enables accurate analysis of HLA polymorphisms in historical (and modern) human populations.

A mobile laboratory for ancient DNA analysis

Mobile devices for on-field DNA analysis have been used for medical diagnostics at the point-of-care, forensic investigations and environmental surveys, but still have to be validated for ancient DNA studies. We report here on a mobile laboratory that we setup using commercially available devices, including a compact real-time PCR machine, and describe procedures to perform DNA extraction and analysis from a variety of archeological samples within 4 hours. The process is carried out on 50 mg samples that are identified at the species level using custom TaqMan real-time PCR assays for mitochondrial DNA fragments. We evaluated the potential of this approach in museums lacking facilities for DNA studies by analyzing samples from the Enlène (MIS 2 layer) and the Portel-Ouest cave (MIS 3 deposits), and also performed experiments during an excavation campaign at the Roc-en-Pail (MIS 5) open-air site. Enlène Bovinae bone samples only yielded DNA for the extinct steppe bison (Bison priscus), whereas Portel-Ouest cave coprolites contained cave hyena (Crocuta crocuta spelaea) DNA together, for some of them, with DNA for the European bison sister species/subspecies (Bison schoetensacki/Bb1-X), thus highlighting the cave hyena diet. Roc-en-Pail Bovinae bone and tooth samples also contained DNA for the Bison schoetensacki/Bb1-X clade, and Cervidae bone samples only yielded reindeer (Rangifer tarandus) DNA. Subsequent DNA sequencing analyses confirmed that correct species identification had been achieved using our TaqMan assays, hence validating these assays for future studies. We conclude that our approach enables the rapid genetic characterization of tens of millennia-old archeological samples and is expected to be useful for the on-site screening of museums and freshly excavated samples for DNA content. Because our mobile laboratory is made up of commercially available instruments, this approach is easily accessible to other investigators.

Distinguishing extant elephants ivory from mammoth ivory using a short sequence of cytochrome b gene

Trade in ivory from extant elephant species namely Asian elephant (Elephas maximus), African savanna elephant (Loxodonta africana) and African forest elephant (Loxodonta cyclotis) is regulated internationally, while the trade in ivory from extinct species of Elephantidae, including woolly mammoth, is unregulated. This distinction creates opportunity for laundering and trading elephant ivory as mammoth ivory. The existing morphological and molecular genetics methods do not reliably distinguish the source of ivory items that lack clear identification characteristics or for which the quality of extracted DNA cannot support amplification of large gene fragments. We present a PCR-sequencing method based on 116 bp target sequence of the cytochrome b gene to specifically amplify elephantid DNA while simultaneously excluding non-elephantid species and ivory substitutes, and while avoiding contamination by human DNA. The partial Cytochrome b gene sequence enabled accurate association of ivory samples with their species of origin for all three extant elephants and from mammoth. The detection limit of the PCR system was as low as 10 copy numbers of target DNA. The amplification and sequencing success reached 96.7% for woolly mammoth ivory and 100% for African savanna elephant and African forest elephant ivory. This is the first validated method for distinguishing elephant from mammoth ivory and it provides forensic support for investigation of ivory laundering cases.

Analysis of ancient mtDNA from the medieval archeological site of Amiternum (L'Aquila), central Italy

Study of ancient DNA makes it possible to analyze genetic relationships between individuals and populations of past and present. In this paper we have analyzed remains of human bones, dating back to the 8th-10th century AD, from the burials found in the Cathedral of Santa Maria in Civitate, archaeological site of Amiternum, L'Aquila, Italy. As a genetic marker, the hypervariable region 1 of mitochondrial DNA (HVR1) was selected. To obtain reliable sequences from the hypervariable region 1 of mtDNA (HVR1) were performed: multiple extractions, template quantification and cloning of PCR products. The sequences obtained were compared with Anderson's sequence for the identification of polymorphisms (SNP) and haplogroups. The data obtained were analyzed with various software and phylogenetic methods. For the comparison between populations, ancient and modern sequences found in databases and literature have been used. This work provides preliminary information on the correlation between the population of Amiternum, the migrant populations transited and/or established in the territory of Amiternum such as Byzantines, Longobards (Lombards), which dominated the Italian peninsula between 568 and 774 AD, and the current populations of Italy.

The study of haplogroups, the analysis of genetic variability and phylogenesis studies on the sequences considered show a genetic closeness between the individuals of Amiternum, the current population of central-northern Italy and the Germanic tribe of Longobards, however, also highlights genetic traits of Byzantines in some samples of Amiternum. Using the analysis of amelogenin gene fragments, we successfully determined the sex of the bone remains on all samples.

Impact of DNA degradation on massively parallel sequencing-based autosomal STR, iiSNP, and mitochondrial DNA typing systems

Biological samples, including skeletal remains exposed to environmental insults for extended periods of time, exhibit increasing levels of DNA damage and fragmentation. Human forensic identification methods typically use a combination of mitochondrial (mt) DNA sequencing and short tandem repeat (STR) analysis, which target segments of DNA ranging from 80 to 500 base pairs (bps). Larger templates are often unavailable as skeletal samples age and the associated DNA degrades. Single-nucleotide polymorphism (SNP) loci target shorter templates and may serve as a solution to the problem. Recently developed assays for STR and SNP analysis using a massively parallel sequencing approach, such as the ForenSeq kit (Verogen, San Diego, CA), offer a means for generating results from degraded samples as they target templates down to 60 to 170 bps. We performed a modeling study that demonstrates that SNPs can increase the significance of an identification when analyzing DNA down to an average size of 100 bps for input amounts between 0.375 and 1 ng of nuclear DNA. Observations from this study were then compared with human skeletal material results (n = 14, ninth to eighteenth centuries), which further demonstrated the utility of the ForenSeq kit for degraded samples. The robustness of the Promega PowerSeq™ Mito System was also tested with human skeletal remains (n = 70, ninth to eighteenth centuries), resulting in successful coverage of 99.29% of the mtDNA control region at 50× coverage or more. This was accompanied by modifications to a mainstream DNA extraction technique for skeletal remains that improved recovery of shorter templates.

Defending Our Public Biological Databases as a Global Critical Infrastructure

Progress in modern biology is being driven, in part, by the large amounts of freely available data in public resources such as the International Nucleotide Sequence Database Collaboration (INSDC), the world's primary database of biological sequence (and related) information. INSDC and similar databases have dramatically increased the pace of fundamental biological discovery and enabled a host of innovative therapeutic, diagnostic, and forensic applications. However, as high-value, openly shared resources with a high degree of assumed trust, these repositories share compelling similarities to the early days of the Internet. Consequently, as public biological databases continue to increase in size and importance, we expect that they will face the same threats as undefended cyberspace. There is a unique opportunity, before a significant breach and loss of trust occurs, to ensure they evolve with quality and security as a design philosophy rather than costly "retrofitted" mitigations. This Perspective surveys some potential quality assurance and security weaknesses in existing open genomic and proteomic repositories, describes methods to mitigate the likelihood of both intentional and unintentional errors, and offers recommendations for risk mitigation based on lessons learned from cybersecurity.

Plasmid DNA contaminant in molecular reagents

Background noise in metagenomic studies is often of high importance and its removal requires extensive post-analytic, bioinformatics filtering. This is relevant as significant signals may be lost due to a low signal-to-noise ratio. The presence of plasmid residues, that are frequently present in reagents as contaminants, has not been investigated so far, but may pose a substantial bias. Here we show that plasmid sequences from different sources are omnipresent in molecular biology reagents. Using a metagenomic approach, we identified the presence of the (pol) of equine infectious anemia virus in human samples and traced it back to the expression plasmid used for generation of a commercial reverse transcriptase. We found fragments of multiple other expression plasmids in human samples as well as commercial polymerase preparations. Plasmid contamination sources included production chain of molecular biology reagents as well as contamination of reagents from environment or human handling of samples and reagents. Retrospective analyses of published metagenomic studies revealed an inaccurate signal-to-noise differentiation. Hence, the plasmid sequences that seem to be omnipresent in molecular biology reagents may misguide conclusions derived from genomic/metagenomics datasets and thus also clinical interpretations. Critical appraisal of metagenomic data sets for the possibility of plasmid background noise is required to identify reliable and significant signals.

Authentication and Assessment of Contamination in Ancient DNA

Contamination from both present-day humans and postmortem microbial sources is a common challenge in ancient DNA studies. Here we present a suite of tools to assist in the assessment of contamination in ancient DNA data sets. These tools perform standard tests of authenticity of ancient DNA data including detecting the presence of postmortem damage signatures in sequence alignments and quantifying the amount of present-day human contamination.

Ancient DNA study reveals HLA susceptibility locus for leprosy in medieval Europeans

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.

Leprosy, caused by infection with Mycobacterium leprae, was common in Europe in the Middle Ages. Here, Krause-Kyora et al. analyze ancient DNA from a medieval Danish leprosarium to assemble 10 complete bacterial genomes and perform association analysis of the DRB1*15:01 allele with risk of leprosy infection.

Maternal DNA lineages at the gate of Europe in the 10th century AD

Given the paucity of archaeogenetic data available for medieval European populations in comparison to other historical periods, the genetic landscape of this age appears as a puzzle of dispersed, small, known pieces. In particular, Southeastern Europe has been scarcely investigated to date. In this paper, we report the study of mitochondrial DNA in 10^th century AD human samples from Capidava necropolis, located in Dobruja (Southeastern Romania, Southeastern Europe). This geographical region is particularly interesting because of the extensive population flux following diverse migration routes, and the complex interactions between distinct population groups during the medieval period. We successfully amplified and typed the mitochondrial control region of 10 individuals. For five of them, we also reconstructed the complete mitochondrial genomes using hybridization-based DNA capture combined with Next Generation Sequencing. We have portrayed the genetic structure of the Capidava medieval population, represented by 10 individuals displaying 8 haplotypes (U5a1c2a, V1a, R0a2’3, H1, U3a, N9a9, H5e1a1, and H13a1a3). Remarkable for this site is the presence of both Central Asiatic (N9a) and common European mtDNA haplotypes, establishing Capidava as a point of convergence between East and West. The distribution of mtDNA lineages in the necropolis highlighted the existence of two groups of two individuals with close maternal relationships as they share the same haplotypes. We also sketch, using comparative statistical and population genetic analyses, the genetic relationships between the investigated dataset and other medieval and modern Eurasian populations.

Improved recovery of ancient DNA from subfossil wood - application to the world's oldest Late Glacial pine forest

Ancient DNA from historical and subfossil wood has a great potential to provide new insights into the history of tree populations. However, its extraction and analysis have not become routine, mainly because contamination of the wood with modern plant material can complicate the verification of genetic information. Here, we used sapwood tissue from 22 subfossil pines that were growing c. 13 000 yr bp in Zurich, Switzerland. We developed and evaluated protocols to eliminate surface contamination, and we tested ancient DNA authenticity based on plastid DNA metabarcoding and the assessment of post-mortem DNA damage. A novel approach using laser irradiation coupled with bleaching and surface removal was most efficient in eliminating contaminating DNA. DNA metabarcoding confirmed which ancient DNA samples repeatedly amplified pine DNA and were free of exogenous plant taxa. Pine DNA sequences of these samples showed a high degree of cytosine to thymine mismatches, typical of post-mortem damage. Stringent decontamination of wood surfaces combined with DNA metabarcoding and assessment of post-mortem DNA damage allowed us to authenticate ancient DNA retrieved from the oldest Late Glacial pine forest. These techniques can be applied to any subfossil wood and are likely to improve the accessibility of relict wood for genome-scale ancient DNA studies.

Cheek tooth morphology and ancient mitochondrial DNA of late Pleistocene horses from the western interior of North America: Implications for the taxonomy of North American Late Pleistocene Equus

Horses were a dominant component of North American Pleistocene land mammal communities and their remains are well represented in the fossil record. Despite the abundant material available for study, there is still considerable disagreement over the number of species of Equus that inhabited the different regions of the continent and on their taxonomic nomenclature. In this study, we investigated cheek tooth morphology and ancient mtDNA of late Pleistocene Equus specimens from the Western Interior of North America, with the objective of clarifying the species that lived in this region prior to the end-Pleistocene extinction. Based on the morphological and molecular data analyzed, a caballine (Equus ferus) and a non-caballine (E. conversidens) species were identified from different localities across most of the Western Interior. A second non-caballine species (E. cedralensis) was recognized from southern localities based exclusively on the morphological analyses of the cheek teeth. Notably the separation into caballine and non-caballine species was observed in the Bayesian phylogenetic analysis of ancient mtDNA as well as in the geometric morphometric analyses of the upper and lower premolars. Teeth morphologically identified as E. conversidens that yielded ancient mtDNA fall within the New World stilt-legged clade recognized in previous studies and this is the name we apply to this group. Geographic variation in morphology in the caballine species is indicated by statistically different occlusal enamel patterns in the specimens from Bluefish Caves, Yukon Territory, relative to the specimens from the other geographic regions. Whether this represents ecomorphological variation and/or a certain degree of geographic and genetic isolation of these Arctic populations requires further study.

Revising mtDNA haplotypes of the ancient Hungarian conquerors with next generation sequencing

As part of the effort to create a high resolution representative sequence database of the medieval Hungarian conquerors we have resequenced the entire mtDNA genome of 24 published ancient samples with Next Generation Sequencing, whose haplotypes had been previously determined with traditional PCR based methods. We show that PCR based methods are prone to erroneous haplotype or haplogroup determination due to ambiguous sequence reads, and many of the resequenced samples had been classified inaccurately. The SNaPshot method applied with published ancient DNA authenticity criteria is the most straightforward and cheapest PCR based approach for testing a large number of coding region SNP-s, which greatly facilitates correct haplogroup determination.

Dental pulp as a source of low-contaminated DNA

The in-laboratory contamination of the ancient samples hinders the result interpretation of the investigations in the field of paleomicrobiology. We had promoted the dental pulp as a sample that limits the risks of in-laboratory contamination of the ancient material. In this work, we measured the contamination of the dental pulp manipulated according to paleomicrobiology protocol, used as a source of a total DNA for metagenomics. First, total DNA extracted from two dog canines was sequenced using next generation sequencing. This yielded a total of 487,828 trimmed reads with a length of 227 ± 35 bp. Sequence analysis of the final dataset using Blast algorithm search and stringent thresholds for sequence identity and coverage against a database including both Canis lupus familiaris and Homo sapiens complete genomes showed that 95% of reads were assigned to C. familiaris whereas 0.03% was assigned to H. sapiens. In a second step, two teeth collected from two 12th century mammals were manipulated following the same protocol. A total of 13,890 trimmed reads with a 157 ± 67 bp length yielded 0-0.35% reads assigned to H. sapiens. This study indicates that the dental pulp is a useful for detecting the significant nucleic sequences in both modern and ancient samples.

Molecular Genealogy of a Mongol Queen's Family and Her Possible Kinship with Genghis Khan

Members of the Mongol imperial family (designated the Golden family) are buried in a secret necropolis; therefore, none of their burial grounds have been found. In 2004, we first discovered 5 graves belonging to the Golden family in Tavan Tolgoi, Eastern Mongolia. To define the genealogy of the 5 bodies and the kinship among them, SNP and/or STR profiles of mitochondria, autosomes, and Y chromosomes were analyzed. Four of the 5 bodies were determined to carry the mitochondrial DNA haplogroup D4, while the fifth carried haplogroup CZ, indicating that this individual had no kinship with the others. Meanwhile, Y-SNP and Y-STR profiles indicate that the males examined belonged to the R1b-M343 haplogroup. Thus, their East Asian D4 or CZ matrilineal and West Eurasian R1b-M343 patrilineal origins reveal genealogical admixture between Caucasoid and Mongoloid ethnic groups, despite a Mongoloid physical appearance. In addition, Y chromosomal and autosomal STR profiles revealed that the four D4-carrying bodies bore the relationship of either mother and three sons or four full siblings with almost the same probability. Moreover, the geographical distribution of R1b-M343-carrying modern-day individuals demonstrates that descendants of Tavan Tolgoi bodies today live mainly in Western Eurasia, with a high frequency in the territories of the past Mongol khanates. Here, we propose that Genghis Khan and his family carried Y-haplogroup R1b-M343, which is prevalent in West Eurasia, rather than the Y-haplogroup C3c-M48, which is prevalent in Asia and which is widely accepted to be present in the family members of Genghis Khan. Additionally, Tavan Tolgoi bodies may have been the product of marriages between the lineage of Genghis Khan’s Borjigin clan and the lineage of either the Ongud or Hongirad clans, indicating that these individuals were members of Genghis Khan’s immediate family or his close relatives.

Human Contamination in Public Genome Assemblies

Contamination in genome assembly can lead to wrong or confusing results when using such genome as reference in sequence comparison. Although bacterial contamination is well known, the problem of human-originated contamination received little attention. In this study we surveyed 45,735 available genome assemblies for evidence of human contamination. We used lineage specificity to distinguish between contamination and conservation. We found that 154 genome assemblies contain fragments that with high confidence originate as contamination from human DNA. Majority of contaminating human sequences were present in the reference human genome assembly for over a decade. We recommend that existing contaminated genomes should be revised to remove contaminated sequence, and that new assemblies should be thoroughly checked for presence of human DNA before submitting them to public databases.

Mitochondrial DNA Suggests a Western Eurasian Origin for Ancient (Proto-) Bulgarians

Ancient (proto-) Bulgarians have long been thought of as a Turkic population. However, evidence found in the past three decades shows that this is not the case. Until now, this evidence has not included ancient mitochondrial DNA (mtDNA) analysis. To fill this void, we collected human remains from the 8th to the 10th century AD located in three necropolises in Bulgaria: Nojarevo (Silistra region) and Monastery of Mostich (Shumen region), both in northeastern Bulgaria, and Tuhovishte (Satovcha region) in southwestern Bulgaria. The phylogenetic analysis of 13 ancient DNA samples (extracted from teeth) identified 12 independent haplotypes, which we further classified into mtDNA haplogroups found in present-day European and western Eurasian populations. Our results suggest a western Eurasian matrilineal origin for proto-Bulgarians, as well as a genetic similarity between proto- and modern Bulgarians. Our future work will provide additional data that will further clarify proto-Bulgarian origins, thereby adding new clues to the current understanding of European genetic evolution.

Inferring chronological age from DNA methylation patterns of human teeth

OBJECTIVE

Current methods to determine chronological age from modern and ancient remains rely on both morphological and molecular approaches. However, low accuracy and the lack of standardized protocols make the development of alternative methods for the estimation of individual's age even more urgent for several research fields, such as biological anthropology, biodemography, forensics, evolutionary genetics, and ancient DNA studies. Therefore, the aim of this study is to identify genomic regions whose DNA methylation level correlates with age in modern teeth.

METHODS

We used MALDI-TOF mass spectrometry to analyze DNA methylation levels of specific CpGs located in the ELOVL2, FHL2, and PENK genes. We considered methylation data from cementum, dentin and pulp of 21 modern teeth (from 17 to 77 years old) to construct a mathematical model able to exploit DNA methylation values to predict age of the individuals.

RESULTS

The median difference between the real age and that estimated using DNA methylation values is 1.20 years (SD = 1.9) if DNA is recovered from both cementum and pulp of the same modern teeth, 2.25 years (SD = 2.5) if DNA is recovered from dental pulp, 2.45 years (SD = 3.3) if DNA is extracted from cementum and 7.07 years (SD = 7.0) when DNA is recovered from dentin only.

DISCUSSION

We propose for the first time the evaluation of DNA methylation at ELOVL2, FHL2, and PENK genes as a powerful tool to predict age in modern teeth for anthropological applications. Future studies are needed to apply this method also to historical and relatively ancient human teeth.

Genealogical relationships between early medieval and modern inhabitants of Piedmont

In the period between 400 to 800 AD, also known as the period of the Barbarian invasions, intense migration is documented in the historical record of Europe. However, little is known about the demographic impact of these historical movements, potentially ranging from negligible to substantial. As a pilot study in a broader project on Medieval Europe, we sampled 102 specimens from 5 burial sites in Northwestern Italy, archaeologically classified as belonging to Lombards or Longobards, a Germanic people ruling over a vast section of the Italian peninsula from 568 to 774. We successfully amplified and typed the mitochondrial hypervariable region I (HVR-I) of 28 individuals. Comparisons of genetic diversity with other ancient populations and haplotype networks did not suggest that these samples are heterogeneous, and hence allowed us to jointly compare them with three isolated contemporary populations, and with a modern sample of a large city, representing a control for the effects of recent immigration. We then generated by serial coalescent simulations 16 millions of genealogies, contrasting a model of genealogical continuity with one in which the contemporary samples are genealogically independent from the medieval sample. Analyses by Approximate Bayesian Computation showed that the latter model fits the data in most cases, with one exception, Trino Vercellese, in which the evidence was compatible with persistence up to the present time of genetic features observed among this early medieval population. We conclude that it is possible, in general, to detect evidence of genealogical ties between medieval and specific modern populations. However, only seldom did mitochondrial DNA data allow us to reject with confidence either model tested, which indicates that broader analyses, based on larger assemblages of samples and genetic markers, are needed to understand in detail the effects of medieval migration.

Aspects of ancient mitochondrial DNA analysis in different populations for understanding human evolution

The evolution of modern humans is a long and difficult process which started from their first appearance and continues to the present day. The study of the genetic origin of populations can help to determine population kinship and to better understand the gradual changes of the gene pool in space and time. Mitochondrial DNA (mtDNA) is a proper tool for the determination of the origin of populations due to its high evolutionary importance. Ancient mitochondrial DNA retrieved from museum specimens, archaeological finds and fossil remains can provide direct evidence for population origins and migration processes. Despite the problems with contaminations and authenticity of ancient mitochondrial DNA, there is a developed set of criteria and platforms for obtaining authentic ancient DNA. During the last two decades, the application of different methods and techniques for analysis of ancient mitochondrial DNA gave promising results. Still, the literature is relatively poor with information for the origin of human populations. Using comprehensive phylogeographic and population analyses we can observe the development and formation of the contemporary populations. The aim of this study was to shed light on human migratory processes and the formation of populations based on available ancient mtDNA data.

Ancient DNA reveals matrilineal continuity in present-day Poland over the last two millennia

While numerous ancient human DNA datasets from across Europe have been published till date, modern-day Poland in particular, remains uninvestigated. Besides application in the reconstruction of continent-wide human history, data from this region would also contribute towards our understanding of the history of the Slavs, whose origin is hypothesized to be in East or Central Europe. Here, we present the first population-scale ancient human DNA study from the region of modern-day Poland by establishing mitochondrial DNA profiles for 23 samples dated to 200 BC – 500 AD (Roman Iron Age) and for 20 samples dated to 1000–1400 AD (Medieval Age). Our results show that mitochondrial DNA sequences from both periods belong to haplogroups that are characteristic of contemporary West Eurasia. Haplotype sharing analysis indicates that majority of the ancient haplotypes are widespread in some modern Europeans, including Poles. Notably, the Roman Iron Age samples share more rare haplotypes with Central and Northeast Europeans, whereas the Medieval Age samples share more rare haplotypes with East-Central and South-East Europeans, primarily Slavic populations. Our data demonstrates genetic continuity of certain matrilineages (H5a1 and N1a1a2) in the area of present-day Poland from at least the Roman Iron Age until present. As such, the maternal gene pool of present-day Poles, Czechs and Slovaks, categorized as Western Slavs, is likely to have descended from inhabitants of East-Central Europe during the Roman Iron Age.

Sex determination in highly fragmented human DNA by high-resolution melting (HRM) analysis

Sex identification in ancient human remains is a common problem especially if the skeletons are sub-adult, incomplete or damaged. In this paper we propose a new method to identify sex, based on real-time PCR amplification of small fragments (61 and 64 bp) of the third exon within the amelogenin gene covering a 3-bp deletion on the AMELX-allele, followed by a High Resolution Melting analysis (HRM). HRM is based on the melting curves of amplified fragments. The amelogenin gene is located on both chromosomes X and Y, showing dimorphism in length. This molecular tool is rapid, sensitive and reduces the risk of contamination from exogenous genetic material when used for ancient DNA studies. The accuracy of the new method described here has been corroborated by using control samples of known sex and by contrasting our results with those obtained with other methods. Our method has proven to be useful even in heavily degraded samples, where other previously published methods failed. Stochastic problems such as the random allele drop-out phenomenon are expected to occur in a less severe form, due to the smaller fragment size to be amplified. Thus, their negative effect could be easier to overcome by a proper experimental design.

Targeted sampling of cementum for recovery of nuclear DNA from human teeth and the impact of common decontamination measures

Background

Teeth are a valuable source of DNA for identification of fragmented and degraded human remains. While the value of dental pulp as a source of DNA is well established, the quantity and presentation of DNA in the hard dental tissues has not been extensively studied. Without this knowledge common decontamination, sampling and DNA extraction techniques may be suboptimal. Targeted sampling of specific dental tissues could maximise DNA profiling success, while minimising the need for laborious sampling protocols and DNA extraction techniques, thus improving workflows and efficiencies. We aimed to determine the location of cellular DNA in non-degraded human teeth to quantify the yield of nuclear DNA from cementum, the most accessible and easily sampled dental tissue, and to investigate the effect of a common decontamination method, treatment with sodium hypochlorite (bleach).

We examined teeth histologically and subsequently quantified the yield of nuclear DNA from the cementum of 66 human third molar teeth. We also explored the effects of bleach (at varying concentrations and exposure times) on nuclear DNA within teeth, using histological and quantitative PCR methods.

Results

Histology confirmed the presence of nucleated cells within pulp and cementum, but not in dentine. Nuclear DNA yields from cementum varied substantially between individuals but all samples gave sufficient DNA (from as little as 20 mg of tissue) to produce full short tandem repeat (STR) profiles. Variation in yield between individuals was not influenced by chronological age or sex of the donor. Bleach treatment with solutions as dilute as 2.5% for as little as 1 min damaged the visible nuclear material and reduced DNA yields from cementum by an order of magnitude.

Conclusions

Cementum is a valuable, and easily accessible, source of nuclear DNA from teeth, and may be a preferred source where large numbers of individuals need to be sampled quickly (for example, mass disaster victim identification) without the need for specialist equipment or from diseased and degraded teeth, where pulp is absent. Indiscriminant sampling and decontamination protocols applied to the outer surface of teeth can destroy this DNA, reducing the likelihood of successful STR typing results.

mtDNA from the early Bronze Age to the Roman period suggests a genetic link between the Indian subcontinent and Mesopotamian cradle of civilization

Ancient DNA methodology was applied to analyse sequences extracted from freshly unearthed remains (teeth) of 4 individuals deeply deposited in slightly alkaline soil of the Tell Ashara (ancient Terqa) and Tell Masaikh (ancient Kar-Assurnasirpal) Syrian archaeological sites, both in the middle Euphrates valley. Dated to the period between 2.5 Kyrs BC and 0.5 Kyrs AD the studied individuals carried mtDNA haplotypes corresponding to the M4b1, M49 and/or M61 haplogroups, which are believed to have arisen in the area of the Indian subcontinent during the Upper Paleolithic and are absent in people living today in Syria. However, they are present in people inhabiting today’s Tibet, Himalayas, India and Pakistan. We anticipate that the analysed remains from Mesopotamia belonged to people with genetic affinity to the Indian subcontinent since the distribution of identified ancient haplotypes indicates solid link with populations from the region of South Asia-Tibet (Trans-Himalaya). They may have been descendants of migrants from much earlier times, spreading the clades of the macrohaplogroup M throughout Eurasia and founding regional Mesopotamian groups like that of Terqa or just merchants moving along trade routes passing near or through the region. None of the successfully identified nuclear alleles turned out to be ΔF508 CFTR, LCT-13910T or Δ32 CCR5.