Investigating kinship of Neolithic post-LBK human remains from Krusza Zamkowa, Poland using ancient DNA

Development and validation of YARN: A novel SE-400 MPS kit for East Asian paternal lineage analysis

Y-chromosomal short tandem repeat polymorphisms (Y-STRs) and Y-chromosomal single nucleotide polymorphisms (Y-SNPs) are valuable genetic markers used in paternal lineage identification and population genetics. Currently, there is a lack of an effective panel that integrates Y-STRs and Y-SNPs for studying paternal lineages, particularly in East Asian populations. Hence, we developed a novel Y-chromosomal targeted panel called YARN (Y-chromosome Ancestry and Region Network) based on multiplex PCR and a single-end 400 massive parallel sequencing (MPS) strategy, consisting of 44 patrilineage Y-STRs and 260 evolutionary Y-SNPs. A total of 386 reactions were validated for the effectiveness and applicability of YARN according to SWGDAM validation guidelines, including sensitivity (with a minimum input gDNA of 0.125 ng), mixture identification (ranging from 1:1-1:10), PCR inhibitor testing (using substances such as 50 μM hematin, 100 μM hemoglobin, 100 μM humic acid, and 2.5 mM indigo dye), species specificity (successfully distinguishing humans from other animals), repeatability study (achieved 100% accuracy), and concordance study (with 99.91% accuracy for 1121 Y-STR alleles). Furthermore, we conducted a pilot study using YARN in a cohort of 484 Han Chinese males from Huaiji County, Zhaoqing City, Guangdong, China (GDZQHJ cohort). In this cohort, we identified 52 different Y-haplogroups and 73 different surnames. We found weak to moderate correlations between the Y-haplogroups, Chinese surnames, and geographical locations of the GDZQHJ cohort (with λ values ranging from 0.050 to 0.340). However, when we combined two different categories into a new independent variable, we observed stronger correlations (with λ values ranging from 0.617 to 0.754). Overall, the YARN panel, which combines Y-STR and Y-SNP genetic markers, meets forensic DNA quality assurance guidelines and holds potential for East Asian geographical origin inference and paternal lineage analysis.

Genetic continuity, isolation, and gene flow in Stone Age Central and Eastern Europe

The genomic landscape of Stone Age Europe was shaped by multiple migratory waves and population replacements, but different regions do not all show similar patterns. To refine our understanding of the population dynamics before and after the dawn of the Neolithic, we generated and analyzed genomic sequence data from human remains of 56 individuals from the Mesolithic, Neolithic, and Eneolithic across Central and Eastern Europe. We found that Mesolithic European populations formed a geographically widespread isolation-by-distance zone ranging from Central Europe to Siberia, which was already established 10,000 years ago. We found contrasting patterns of population continuity during the Neolithic transition: people around the lower Dnipro Valley region, Ukraine, showed continuity over 4000 years, from the Mesolithic to the end of the Neolithic, in contrast to almost all other parts of Europe where population turnover drove this cultural change, including vast areas of Central Europe and around the Danube River.

Patrilocality and hunter-gatherer-related ancestry of populations in East-Central Europe during the Middle Bronze Age

The demographic history of East-Central Europe after the Neolithic period remains poorly explored, despite this region being on the confluence of various ecological zones and cultural entities. Here, the descendants of societies associated with steppe pastoralists form Early Bronze Age were followed by Middle Bronze Age populations displaying unique characteristics. Particularly, the predominance of collective burials, the scale of which, was previously seen only in the Neolithic. The extent to which this re-emergence of older traditions is a result of genetic shift or social changes in the MBA is a subject of debate. Here by analysing 91 newly generated genomes from Bronze Age individuals from present Poland and Ukraine, we discovered that Middle Bronze Age populations were formed by an additional admixture event involving a population with relatively high proportions of genetic component associated with European hunter-gatherers and that their social structure was based on, primarily patrilocal, multigenerational kin-groups.

Kinship analysis of 5th- to 6th-century skeletons of Romanized indigenous people from the Bled-Pristava archaeological site

The familial relationship between skeletons buried together in a shared grave is important for understanding the burial practices of past human populations. Four skeletons were excavated from the Late Antiquity part of the Bled-Pristava burial site in Slovenia, dated to the 5th to 6th century. They were anthropologically characterized as two adults (a middle-aged man and a young woman) and two non-adults (of unknown sex). Based on stratigraphy, the skeletons were considered to be buried simultaneously in one grave. Our aim was to determine whether the skeletons were related. Petrous bones and teeth were used for genetic analysis. Specific precautions were followed to prevent contamination of ancient DNA with contemporary DNA, and an elimination database was established. Bone powder was obtained using a MillMix tissue homogenizer. Prior to extracting the DNA using Biorobot EZ1, 0.5 g of powder was decalcified. The PowerQuant System was used for quantification, various autosomal kits for autosomal short tandem repeat (STR) typing, and the PowerPlex Y23 kit for Y-STR typing. All analyses were performed in duplicate. Up to 28 ng DNA/g of powder was extracted from the samples analyzed. Almost full autosomal STR profiles obtained from all four skeletons and almost full Y-STR haplotypes obtained from two male skeletons were compared, and the possibility of a familial relationship was evaluated. No amplification was obtained in the negative controls, and no match was found in the elimination database. Autosomal STR statistical calculations confirmed that the adult male was the father of two non-adult individuals and one young adult individual from the grave. The relationship between the males (father and son) was additionally confirmed by an identical Y-STR haplotype that belonged to the E1b1b haplogroup, and a combined likelihood ratio for autosomal and Y-STRs was calculated. Kinship analysis confirmed with high confidence (kinship probability greater than 99.9% was calculated for all three children) that all four skeletons belonged to the same family (a father, two daughters, and a son). Through genetic analysis, the burial of members of the same family in a shared grave was confirmed as a burial practice of the population living in the Bled area in Late Antiquity.

Maternal genetic origin of the late and final Neolithic human populations from present-day Poland

OBJECTIVE

We aim to identify maternal genetic affinities between the Middle to Final Neolithic (3850-2300 BC) populations from present-day Poland and possible genetic influences from the Pontic steppe.

MATERIALS AND METHODS

We conducted ancient DNA studies from populations associated with Złota, Globular Amphora, Funnel Beaker, and Corded Ware cultures (CWC). We sequenced genomic libraries on Illumina platform to generate 86 complete ancient mitochondrial genomes. Some of the samples were enriched for mitochondrial DNA using hybridization capture.

RESULTS

The maternal genetic composition found in Złota-associated individuals resembled that found in people associated with the Globular Amphora culture which indicates that both groups likely originated from the same maternal genetic background. Further, these two groups were closely related to the Funnel Beaker culture-associated population. None of these groups shared a close affinity to CWC-associated people. Haplogroup U4 was present only in the CWC group and absent in Złota group, Globular Amphora, and Funnel Beaker cultures.

DISCUSSION

The prevalence of mitochondrial haplogroups of Neolithic farmer origin identified in Early, Middle and Late Neolithic populations suggests a genetic continuity of these maternal lineages in the studied area. Although overlapping in time - and to some extent - in cultural expressions, none of the studied groups (Złota, Globular Amphora, Funnel Beaker), shared a close genetic affinity to CWC-associated people, indicating a larger extent of cultural influence from the Pontic steppe than genetic exchange. The higher frequency of haplogroup U5b found in populations associated with Funnel Beaker, Globular Amphora, and Złota cultures suggest a gradual maternal genetic influx from Mesolithic hunter-gatherers. Moreover, presence of haplogroup U4 in Corded Ware groups is most likely associated with the migrations from the Pontic steppe at the end of the Neolithic and supports the observed genetic distances.

Biological kinship in 750 year old human remains from Central Argentina with signs of interpersonal violence

Human skeletal remains of an adult male (20-24 years old) and a juvenile (4-8 years old), dated to 750 ± 85 ¹⁴C years BP, were found on the southern margin of Mar Chiquita Lagoon (Córdoba, Argentina). Both individuals show signs of being victims of interpersonal violence, with arrowheads associated with the remains and perimortem lesions on the juvenile, as well as an unusual form of burial, with the juvenile partially overlapped with the adult. The aim of this work is to study a possible kin relationship between these two individuals through ancient DNA analysis. Biological kinship was evaluated by autosomal and Y-chromosome STR (short tandem repeat) typing, PCR-APLP for SNP determination and hypervariable region I sequencing of the mitochondrial DNA. Genetic analyses indicated that these individuals shared the same Y-chromosomal haplotype but different mitochondrial lineages. The likelihood ratio based on autosomal loci indicates that the genetic profiles of the human remains would be more likely to be that indicating a father-son bond. The paleogenetic approach combined with forensic genetic methods applied to this study allowed us to confirm a hypothesis that originated in bioarchaeological evidence. This study constitutes a unique case in Argentina of kinship determination based on DNA profiles of human remains in an archaeological context of interpersonal violence. It is important to highlight the contribution made by these studies to address topics usually hidden in bioarchaeological studies, such as community organization, cultural customs and mortuary practices.

AmtDB: a database of ancient human mitochondrial genomes

Ancient mitochondrial DNA is used for tracing human past demographic events due to its population-level variability. The number of published ancient mitochondrial genomes has increased in recent years, alongside with the development of high-throughput sequencing and capture enrichment methods. Here, we present AmtDB, the first database of ancient human mitochondrial genomes. Release version contains 1107 hand-curated ancient samples, freely accessible for download, together with the individual descriptors, including geographic location, radiocarbon dating, and archaeological culture affiliation. The database also features an interactive map for sample location visualization. AmtDB is a key platform for ancient population genetic studies and is available at https://amtdb.org.

Ancient Mitochondrial Genomes Reveal the Absence of Maternal Kinship in the Burials of Çatalhöyük People and Their Genetic Affinities

Çatalhöyük is one of the most widely recognized and extensively researched Neolithic settlements. The site has been used to discuss a wide range of aspects associated with the spread of the Neolithic lifestyle and the social organization of Neolithic societies. Here, we address both topics using newly generated mitochondrial genomes, obtained by direct sequencing and capture-based enrichment of genomic libraries, for a group of individuals buried under a cluster of neighboring houses from the classical layer of the site’s occupation. Our data suggests a lack of maternal kinship between individuals interred under the floors of Çatalhöyük buildings. The findings could potentially be explained either by a high variability of maternal lineages within a larger kin group, or alternatively, an intentional selection of individuals for burial based on factors other than biological kinship. Our population analyses shows that Neolithic Central Anatolian groups, including Çatalhöyük, share the closest affinity with the population from the Marmara Region and are, in contrast, set further apart from the Levantine populations. Our findings support the hypothesis about the emergence and the direction of spread of the Neolithic within Anatolian Peninsula and beyond, emphasizing a significant role of Central Anatolia in this process.

Mitochondrial genomes reveal an east to west cline of steppe ancestry in Corded Ware populations

From around 4,000 to 2,000 BC the forest-steppe north-western Pontic region was occupied by people who shared a nomadic lifestyle, pastoral economy and barrow burial rituals. It has been shown that these groups, especially those associated with the Yamnaya culture, played an important role in shaping the gene pool of Bronze Age Europeans, which extends into present-day patterns of genetic variation in Europe. Although the genetic impact of these migrations from the forest-steppe Pontic region into central Europe have previously been addressed in several studies, the contribution of mitochondrial lineages to the people associated with the Corded Ware culture in the eastern part of the North European Plain remains contentious. In this study, we present mitochondrial genomes from 23 Late Eneolithic and Bronze Age individuals, including representatives of the north-western Pontic region and the Corded Ware culture from the eastern part of the North European Plain. We identified, for the first time in ancient populations, the rare mitochondrial haplogroup X4 in two Bronze Age Catacomb culture-associated individuals. Genetic similarity analyses show close maternal genetic affinities between populations associated with both eastern and Baltic Corded Ware culture, and the Yamnaya horizon, in contrast to larger genetic differentiation between populations associated with western Corded Ware culture and the Yamnaya horizon. This indicates that females with steppe ancestry contributed to the formation of populations associated with the eastern Corded Ware culture while more local people, likely of Neolithic farmer ancestry, contributed to the formation of populations associated with western Corded Ware culture.

Whole mitochondrial genome diversity in two Hungarian populations

Complete mitochondrial genomics is an effective tool for studying the demographic history of human populations, but there is still a deficit of mitogenomic data in European populations. In this paper, we present results of study of variability of 80 complete mitochondrial genomes in two Hungarian populations from eastern part of Hungary (Szeged and Debrecen areas). The genetic diversity of Hungarian mitogenomes is remarkably high, reaching 99.9% in a combined sample. According to the analysis of molecular variance (AMOVA), European populations showed a low, but statistically significant level of between-population differentiation (Fst = 0.61%, p = 0), and two Hungarian populations demonstrate lack of between-population differences. Phylogeographic analysis allowed us to identify 71 different mtDNA sub-clades in Hungarians, sixteen of which are novel. Analysis of ancestry-informative mtDNA sub-clades revealed a complex genetic structure associated with the genetic impact of populations from different parts of Eurasia, though the contribution from European populations is the most pronounced. At least 8% of ancestry-informative haplotypes found in Hungarians demonstrate similarity with East and West Slavic populations (sub-clades H1c23a, H2a1c1, J2b1a6, T2b25a1, U4a2e, K1c1j, and I1a1c), while the influence of Siberian populations is not so noticeable (sub-clades A12a, C4a1a, and probably U4b1a4).

Analysis of LCT-13910 genotypes and bone mineral density in ancient skeletal materials

The relation of LCT-13910 genotypes and bone mineral density (BMD) has been the subject of modern-day human population studies, giving inconsistent results. In the present study we analyze for the first time a relation of LCT-13910 genotypes and BMD in historical skeletal individuals. Ancient population might be a model for testing this association due to elimination of non-natural factors affecting bone density. Among 22 medieval individuals from Sanok churchyard (South-Eastern Poland; dated from XIV to XVII c. AD) we identified 4 individuals with osteoporosis (mean BMD = 0.468 g/cm2, SD = 0.090), 10 individuals with osteopenia (mean BMD = 0.531 g/cm2, SD = 0.066) and 8 individuals with normal BMD values (mean BMD = 0,642 g/cm2, SD = 0.060). Analyses of BMD and LCT-13910 genotypes revealed that mean BMD was the highest (0.583 g/cm2, SD = 0.065) in the individuals with lactose tolerance genotypes (TT and CT). We also found possible association of lower BMD at the radius and CC genotypes due to higher but not statistically significant frequency of osteoporosis in the lactose intolerant group (p = 0.60). Statistically significant correlation was found between BMD and females aged 20-35 years, with tendency to reduce BMD with age (p = 0.02).

Biological Sexing of a 4000-Year-Old Egyptian Mummy Head to Assess the Potential of Nuclear DNA Recovery from the Most Damaged and Limited Forensic Specimens

High throughput sequencing (HTS) has been used for a number of years in the field of paleogenomics to facilitate the recovery of small DNA fragments from ancient specimens. Recently, these techniques have also been applied in forensics, where they have been used for the recovery of mitochondrial DNA sequences from samples where traditional PCR-based assays fail because of the very short length of endogenous DNA molecules. Here, we describe the biological sexing of a ~4000-year-old Egyptian mummy using shotgun sequencing and two established methods of biological sex determination (R_X and R_Y), by way of mitochondrial genome analysis as a means of sequence data authentication. This particular case of historical interest increases the potential utility of HTS techniques for forensic purposes by demonstrating that data from the more discriminatory nuclear genome can be recovered from the most damaged specimens, even in cases where mitochondrial DNA cannot be recovered with current PCR-based forensic technologies. Although additional work remains to be done before nuclear DNA recovered via these methods can be used routinely in operational casework for individual identification purposes, these results indicate substantial promise for the retrieval of probative individually identifying DNA data from the most limited and degraded forensic specimens.

A review of bioinformatic methods for forensic DNA analyses

Short tandem repeats, single nucleotide polymorphisms, and whole mitochondrial analyses are three classes of markers which will play an important role in the future of forensic DNA typing. The arrival of massively parallel sequencing platforms in forensic science reveals new information such as insights into the complexity and variability of the markers that were previously unseen, along with amounts of data too immense for analyses by manual means. Along with the sequencing chemistries employed, bioinformatic methods are required to process and interpret this new and extensive data. As more is learnt about the use of these new technologies for forensic applications, development and standardization of efficient, favourable tools for each stage of data processing is being carried out, and faster, more accurate methods that improve on the original approaches have been developed. As forensic laboratories search for the optimal pipeline of tools, sequencer manufacturers have incorporated pipelines into sequencer software to make analyses convenient. This review explores the current state of bioinformatic methods and tools used for the analyses of forensic markers sequenced on the massively parallel sequencing (MPS) platforms currently most widely used.

Comprehensive analysis of microorganisms accompanying human archaeological remains

Metagenome analysis has become a common source of information about microbial communities that occupy a wide range of niches, including archaeological specimens. It has been shown that the vast majority of DNA extracted from ancient samples come from bacteria (presumably modern contaminants). However, characterization of microbial DNA accompanying human remains has never been done systematically for a wide range of different samples. We used metagenomic approaches to perform comparative analyses of microorganism communities present in 161 archaeological human remains. DNA samples were isolated from the teeth of human skeletons dated from 100 AD to 1200 AD. The skeletons were collected from 7 archaeological sites in Central Europe and stored under different conditions. The majority of identified microbes were ubiquitous environmental bacteria that most likely contaminated the host remains not long ago. We observed that the composition of microbial communities was sample-specific and not correlated with its temporal or geographical origin. Additionally, traces of bacteria and archaea typical for human oral/gut flora, as well as potential pathogens, were identified in two-thirds of the samples. The genetic material of human-related species, in contrast to the environmental species that accounted for the majority of identified bacteria, displayed DNA damage patterns comparable with endogenous human ancient DNA, which suggested that these microbes might have accompanied the individual before death. Our study showed that the microbiome observed in an individual sample is not reliant on the method or duration of sample storage. Moreover, shallow sequencing of DNA extracted from ancient specimens and subsequent bioinformatics analysis allowed both the identification of ancient microbial species, including potential pathogens, and their differentiation from contemporary species that colonized human remains more recently.

Late Danubian mitochondrial genomes shed light into the Neolithisation of Central Europe in the 5th millennium BC

Background

Recent aDNA studies are progressively focusing on various Neolithic and Hunter - Gatherer (HG) populations, providing arguments in favor of major migrations accompanying European Neolithisation. The major focus was so far on the Linear Pottery Culture (LBK), which introduced the Neolithic way of life in Central Europe in the second half of 6th millennium BC. It is widely agreed that people of this culture were genetically different from local HGs and no genetic exchange is seen between the two groups. From the other hand some degree of resurgence of HGs genetic component is seen in late Neolithic groups belonging to the complex of the Funnel Beaker Cultures (TRB). Less attention is brought to various middle Neolithic cultures belonging to Late Danubian sequence which chronologically fall in between those two abovementioned groups. We suspected that genetic influx from HG to farming communities might have happened in Late Danubian cultures since archaeologists see extensive contacts between those two communities.

Results

Here we address this issue by presenting 5 complete mitochondrial genomes of various late Danubian individuals from modern-day Poland and combining it with available published data. Our data show that Late Danubian cultures are maternally closely related to Funnel Beaker groups instead of culturally similar LBK.

Conclusions

We assume that it is an effect of the presence of individuals belonging to U5 haplogroup both in Late Danubians and the TRB. The U5 haplogroup is thought to be a typical for HGs of Europe and therefore we argue that it is an additional evidence of genetic exchange between farming and HG groups taking place at least as far back as in middle Neolithic, in the Late Danubian communities.

Electronic supplementary material

The online version of this article (doi:10.1186/s12862-017-0924-0) contains supplementary material, which is available to authorized users.

Diverse origin of mitochondrial lineages in Iron Age Black Sea Scythians

Scythians were nomadic and semi-nomadic people that ruled the Eurasian steppe during much of the first millennium BCE. While having been extensively studied by archaeology, very little is known about their genetic identity. To fill this gap, we analyzed ancient mitochondrial DNA (mtDNA) from Scythians of the North Pontic Region (NPR) and successfully retrieved 19 whole mtDNA genomes. We have identified three potential mtDNA lineage ancestries of the NPR Scythians tracing back to hunter-gatherer and nomadic populations of east and west Eurasia as well as the Neolithic farming expansion into Europe. One third of all mt lineages in our dataset belonged to subdivisions of mt haplogroup U5. A comparison of NPR Scythian mtDNA linages with other contemporaneous Scythian groups, the Saka and the Pazyryks, reveals a common mtDNA package comprised of haplogroups H/H5, U5a, A, D/D4, and F1/F2. Of these, west Eurasian lineages show a downward cline in the west-east direction while east Eurasian haplogroups display the opposite trajectory. An overall similarity in mtDNA lineages of the NPR Scythians was found with the late Bronze Age Srubnaya population of the Northern Black Sea region which supports the archaeological hypothesis suggesting Srubnaya people as ancestors of the NPR Scythians.