|
Post by Admin on Mar 28, 2022 20:58:23 GMT
Genetic outliers from previous studies Many samples were defined as genetic outliers in their genetic context by previous studies from Bronze Age Europe. We selected such outlier individuals with high HG ancestry components to assess whether they are related to Bk-II. Selection was based on previous observations and also by using Dixon’s Q-test25 at 90% confidence interval on HG component’s upper deviation using results of the Admixture analysis. First, we ran f4-statistics in the form of f4(W=test outlier, X=corresponding population, Y=Bk-II, Z=Yoruba)26. This test resulted in positive values for some outliers (W) meaning that these are genetically closer to Bk-II (Y) than its presumed population (X). However, Z-scores are low in many cases, and false positives may appear solely by high HG component, not by true relationship (see Supplementary Table S10 and Supplementary Information section 5.3). To check true relationship between Bk-II and groups/samples with high HG ancestry, we performed an outgroup f3-statistics in the form of f3(X=Bk-II, Y=test HG-s, Z=Yoruba) for all relevant archaic populations and outliers23 (Supplementary Table S11) that resulted a table of allele frequency based distances between test pairs Xs and Ys. Euclidean clustering based on the results of f3-statistics revealed that a number of samples and even three populations from the Baltic (Fig. 3 b, Supplementary Information section 5.4) from AADR23 form a cluster with Bk-II and Bk-III, suggesting actual genetic connection via a common HG ancestry source. Origin of high HG ancestry component in BK-II and the Kisapostag population f4-statistics in the form of f4(W=Sweden HG, X=Serbia IronGates HG, Y=Bk-II, Z=Yoruba) revealed that while Pitted Ware culture associated individuals have similarly high HG levels to Bk-II from Scandinavian HG-s, in Bk-II this type of HG component is significantly negative meaning no or minimal shared ancestry (Supplementary Information section 5.3). Results of outgroup in f3-statistics in the form of f3(X=Bk-II, Y=HG-s, Z=Yoruba) showed that HGs fit the best to Bk-II are geographically widely distributed (Supplementary Table S11), thus we conclude that the true source for this component is yet to be described. We can take into consideration the outgroup f3-statistics results, chronology, timing of HG admixture according to Freilich 202110 between ~3400-2400 BCE, qpAdm results and the geographical distribution of groups and outliers of similar HG makeup. These suggest a dated migration pattern for this undescribed population with dominant HG genetic ancestry from what is present-day Bulgaria to the Baltic through the Eastern borders of the Carpathians (Fig. 3 b). Ancestors of Bk-II likely branched off from this migration route and started to move towards West, by at least around ~2500 BCE, subsequently intermixing with various groups. Interestingly, the phylogeography of mitochondrial haplogroup U4b1b1 (individual S15) perfectly fits this scenario. Notably, strontium (87Sr/86Sr ratio) isotope results from molars (Fig. 4, Supplementary Information section 1.9) indicate that every individual of this study was raised and lived at least close proximity to Balatonkeresztúr site in their childhood and early adolescence, suggesting that Bk-II group does not represent the first generation of the Kisapostag culture associated population in Transdanubia (Western Hungary). Correspondingly, the high number of supported models in qpAdm results could be the result of subsequent local admixtures. Fig. 4. Sr isotope data from the Balatonkeresztúr site. Samples were taken from dental enamel (first, second and third molars) to evaluate whether individuals were born in the area, or grew up in a geologically distinct region. All of the samples are consistent with previously published plant and water 87Sr/86Sr ratio (green diamonds) data collected from the southern portion of Lake Balaton14. For further data see Supplementary Information section 1.9.
|
|
|
Post by Admin on Mar 29, 2022 0:37:43 GMT
A Late Copper Age outlier individual from Balatonlelle site We included a Late Copper Age individual from Balatonlelle site to this study for its high HG genomic ancestry component. Mitogenome of BAD002 (K1a4a1) shows affinity to Iberian BBC associated individuals (Supplementary Information Fig. S.2.1.1), and his Y chromosomal haplogroup belongs to I-M170. On genomic PCA this sample groups with Iberian and French Neolithic individuals, due to higher HG component (~34%) than known Neolithic and Copper Age populations in the Carpathian Basin have, and due to the lack of steppe related ancestry. qpAdm estimates pointed to a source of Neolithic communities from present-day (Northwestern) France (87±8%)27 with a minor extra HG component (13±8%) with p=0.948 (Supplementary Information section 5.5, Supplementary Tables S12 and S13). Pigmentation pattern of BAD002 shows resemblance to average Neolithic Europeans. The foreign cultural traits of the boy’s jewellery is in line with his outlier genetic composition in the study region28. Notably, further tests (outgroup f3-statistics, qpAdm) excluded contribution nor any connection of BAD002 to Bk-II (Supplementary Tables S11, S12 and S13, Supplementary Information section 5). Therefore we conclude that this individual testifies large-scale migration in the Copper Age, providing research questions for future studies.
PAPline We introduce our newly developed, freely available bioinformatic pipeline, named PAPline (Performing Archaeogenetic Pipeline), written in linux bash, R, and python v3.8.10 programming languages. One can use this software primarily to analyse next generation sequencing data of archaeogenomic samples, supplemented by non-pipeline scripts. The main distinguishing feature of the PAPline compared to the EAGER29 and the Paleomix30 pipelines is its user friendliness, as its installation process is less complicated and it provides a graphical interface with almost complete automatisation with practical details. For detailed description visit the github page or see Supplementary Information section 6.
|
|
|
Post by Admin on Mar 29, 2022 14:45:27 GMT
Discussion The Carpathian Basin was inhabited by the Baden cultures’ population at the end of the Copper Age, and their genetic composition was represented by an early farmer and a slightly increased HG genetic component, compared to the previous Neolithic populations of the region6. Here we demonstrated that in the early phase of this culture, a group of Western European origin appeared in Transdanubia, diversifying what we knew about the region’s Late Copper Age substrate up to now.
The Carpathian Basin experienced the influx of steppe-related genetic ancestry at the dawn of the Bronze Age5,8, and this transformation was already detectable at Balatonkeresztúr-Réti-dűlő site as well, where we could examine multiple populations. The earliest Bronze Age horizon Bk-I (representative of the Somogyvár-Vinkovci culture) shows similarities to Poland Southeast BBC associated population with high steppe ancestry that was replaced by the Kisapostag culture associated Bk-II around the 23-22th century BCE, while at least some sort of genetic ancestry of Bk-I in this population can not be excluded. According to our results, the Bk-II population had outstandingly high HG genetic ancestry levels, compared to other Bronze Age groups of the region, which can be traced back to today’s Ukraine, Belarus, Moldavia or Romania, pointing to a long standing previously unsampled population with dominant HG ancestry. Calculated admixture dates10 suggest the presence of a genetically pure or at least highly HG specific population in Eastern Europe as late as the end of the Copper Age. Part of this group subsequently admixed with populations of mainly steppe (likely Poland Southeast BBC) and early farmer (most likely a Globular Amphora culture related) ancestry during their westward migration on a Northern route, leaving genetic traces in Corded Ware culture, BBC, and other Bronze Age populations. The paternal lineage of BK-II was likely linked to the farmer component, as I2a-M223 (upstream to I2a-L1229) was a frequent paternal lineage among Globular Amphora culture and related populations. Looking for the possible source areas of the Kisapostag culture, a number of archaeological theories need to be considered. The pottery decoration technique originated either from Corded Ware in the Middle Dnieper region (Ukraine), epi-Corded Ware groups (northern Carpathians), e.g. Chłopice-Veselé or Nitra groups (Slovakia), the latter two is also supported by inhumation practises31–36. However, connections with the Litzenkeramik or Guntramsdorf-Drassburg group (eastern Austria, Slovenia, western Croatia) were also raised37,38. Pottery forms were connected to local development of communities with eastern (Makó–Kosihy–Čaka) or southern (Somogyvár–Vinkovci) origins, too39. BBC influence was also mentioned based on connections of pottery and craniometry data (so called Glockenbecher or brachycran skull type40–42). The results of this study fit best with the Middle Dnieper area origin of BK-II, especially when we consider individual I4110 from Dereivka I (Ukraine Eneolithic) as one of the earliest representatives of their genomic makeup.
Strontium isotope (87Sr/86Sr ratio) data, representing through nutrition the bioavailable Sr in the area where people lived in a certain age interval, shows local values for both sexes in both Bk-II and Bk-III. These results push back the timing of their arrival a few generations, meaning that local or southern impact of cultural traits and maybe even genetic admixtures likely occurred during this short period as well, which also could explain the culture’s archaeological heterogeneity.
The population of Bk-III was the direct descendant of Bk-II, forming not just cultural (Encrusted pottery) but also genetic continuity for at least ~500 years, even if the radiocarbon sequences allow a few decades of hiatus at the studied site. Continuous female-biased admixture with various groups occurred during this period according to our and previous genetic10 and archaeological31,43 evidence, diluting the BK-II genetic ancestry.
In both periods, the homogeneity of paternal lineages suggest a similar social organisation described in9,10 of a patrilocal residence system. However, strontium isotope data shows local values for both sexes, which along with similar genomic makeup of females and males suggest exogamy most probably between villages of the same population. Two pairs of half-sibling graves in the two periods may indicate polygamy, although remarriage for high female mortality is more plausible. Notably, almost none of the uniparental markers are identical even at the haplogroup level with individuals from the Croatian Encrusted Pottery culture Jagodnjak site, despite high similarities in cultural traits, social structure and genomic composition of the communities. This points to clan-like or patriarchal superfamily structure of Kisapostag and Encrusted Pottery groups. The relatively limited presence of female and children burials in both Bk-II and Bk-III periods may suggest distinctive treatment or another (here undiscovered) burial group for women and children at the same site. Although, in other Bronze Age cemeteries, e.g. Ordacsehi and Bonyhád in Hungary, males, females and children were buried close to each other, suggesting high variance of burial practises34,35,44.
While low genomic coverage did not allow fine SNP recovery, we did find evidence for malignant variants within all of our tested groups, and undoubtedly showed the presence of LHON and Jacob’s syndrome within Bk-II. Additionally, the disease panel we created can be extended and used in future studies, providing insight into past population health qualities.
Considering the unstructured age and kinship distribution in the mass grave Bk-III compared to Bk-II, the coetaneous death of eight people at least, the absence of traumatic or ritual events on bones, and non-cremated nature of the burial all signals a sudden tragic event in the Encrusted pottery period, most likely an epidemic, as first suggested based on the anthropological analyses45. Interestingly, comparative strontium isotope analyses on the first and third molar of the individuals in the BK-III mass grave indicate that subadult males - including a severely disabled individual (S15) with hip dysplasia - left their community for a while and then returned to their birthplace prior to their death, raising further questions for future studies on prehistoric lifeways and social organisations.
|
|
|
Post by Admin on Mar 29, 2022 17:45:55 GMT
Materials and Methods Isotope analyses Radiocarbon dating was performed at the HEKAL AMS C-14 facility of the Institute for Nuclear Research in Debrecen, Hungary (see Supplementary Information section 1.8). 87Sr/86Sr isotope measurements were performed in the ICER Centre, Institute for Nuclear Research Debrecen, Hungary and at Quinnipiac and Yale University, Connecticut, USA (see Supplementary Information section 1.9). Ancient DNA laboratory work Petrous bones and teeth were taken from skulls for genetic investigation (Supplementary Table S1). Laboratory work was performed in a dedicated ancient DNA laboratory facility (Institute of Archaeogenomics, Research Centre for the Humanities, Eötvös Loránd Research Network, Budapest, Hungary). Each step was carried out in separate rooms under sterile conditions, during work protective clothing was used. Irradiated UV-C light, DNA-ExitusPlus™ (AppliChem) and/or bleach were applied for cleaning after and between work stages, and also, blank controls were utilised at all times. Sample surfaces were cleaned by sandblasting and mechanically milled to powder. DNA extraction was performed according to Dabney et al. 201346 with minor changes according to Lipson et al. 20176. DNA extraction success was verified by PCR using mtDNA primer pairs (F16209-R06348; F16045-R06240). Half-UDG treated libraries were used according to Rohland et al. 201547 with minor changes. Unique double internal barcode combinations were used for each library (Supplementary Table S1). Libraries were amplified with TwistAmp Basic (Twist DX Ltd) and purified with AMPure XP beads (Agilent). Then, concentration measurements were taken on Qubit 2.0 fluorometer, fragment sizes were checked on Agilent 4200 TapeStation System (Agilent High Sensitivity D1000 ScreenTape Assay). Hybridisation capture method for mtDNA and 3k nuclear SNP was used besides whole genome shotgun, as described by Haak et al. 2015, Lipson et al. 2017 and Csáky et al. 20204,6,48. Bait production was based on Fu et al. 20161 and N. Rohland’s personal communication, the oligos as a pool was ordered from CustomArray Inc. Both for shotgun and capture libraries, universal iP5 and unique iP7 indexes were used. Sequencing was done on Illumina MiSeq and NovaSeq platforms with custom setup and 150, 200 and 300 cycles, respectively. Additionally, we investigated Y chromosome STR profiles (17 markers) with AmpFISTR® Yfiler® PCR Amplification Kit (Applied Biosystems), having one blank and one positive control at each reaction preparation. The workflow followed the recommended protocol except the PCR cycles were increased from 30 to 34 and reactions were halved in volume. Two repeats were done where at least 4 markers yielded results. Data analyses were carried out in GeneMapper® ID Software v3.2.1 (Applied Biosystems), results are summarised in Supplementary Table S3. Bioinformatic analyses Illumina sequencing reads were processed by the PAPline, for details, see Supplementary Information section 6. We used the GRCH37.p13 reference sequence for calling pseudohaploid genomes. For kinship inferences we applied the READ software49 and a custom script (named MPMR, see Supplementary Information section 2.3 and Supplementary Table S2). MtDNA analyses included phylogenetic analyses using the MrBayes v3.2.650 and the BEAST v1.10.451 softwares and diversity tests using the Popgenome52 R package, see Supplementary Information section 2.1. For Y chromosome haplogroup determination the Yleaf v153 software was applied. For network analysis of STR data we used Network v10.1.0.0 and Network publisher v2.1.2.554,55, see Supplementary Information section 2.2. Due to low genomic coverages (<10,000 SNPs) we discarded individuals S2, S5 and S17 from the population genetic analyses. PCA was made by the Eigensoft smartpca software56 using the Human Origins Panel SNP set26, for other analyses the 1240k array12 was used for SNP call, for results, see Supplementary Table S7. Individuals S4, S5, S6 and S20 were discarded from further tests for them being first degree relatives of other samples. For investigating ancestry estimates we used supervised admixture analysis calculated by the ADMIXTURE v1.3.0 software57. The results were visualised by custom R scripts. f-statistics and qpAdm were performed using the admixr v0.9.158 and the admixtools v2.0.026 R packages. Data availability All other study data are included in the article and/or Supplementary Information and tables. Sequencing data are deposited in the European Nucleotide Archive (ENA) and are available to download under accession number PRJEB49524. Our new software ‘PAPline’ is freely available at www.github.com/gerberd-workshop/papline. Author contribution D.G, V.K., A.Sz-N. conceived and designed the experiments. D.G. processed the sequencing data, created the PAPline, and performed downstream bioinformatic analyses. B.Sz. did all molecular laboratory work. O.Sz. created the mtDNA database for phylogenetic analyses. B.E. obtained Y chromosome STR data. B.Gy. and E.A. optimised genetic analyses. J.I.G., A.H., L.P. performed Sr isotope analysis. G.K., Sz.F., V. K. and M.B. evaluated the archaeological context. B.G.M., Á.K. and K.K. did the anthropological examination of the remains. Á.K. made the facial reconstruction. V.Sz. performed radiocarbon calibrations and modelling. B.G.M. sampled the remains. E.A., V.K. and A.Sz-N. jointly supervised the research and wrote the paper with D.G. All authors provide critical feedback for this study and contribute to the final manuscript. Ethics declarations The authors declare that they had requested and got permission for the destructive bioarchaeological analyses of the archaeological material in the study from the stakeholders, excavator and processor archaeologists. Supplementary Information[supplements/478968_file02.pdf] www.biorxiv.org/content/10.1101/2022.02.03.478968v1.supplementary-material
|
|
|
Post by Admin on Mar 29, 2022 20:14:16 GMT
Genetic Continuity of Bronze Age Ancestry with Increased Steppe-Related Ancestry in Late Iron Age Uzbekistan Vikas Kumar, E Andrew Bennett, Dongyue Zhao, Yun Liang, Yunpeng Tang, Meng Ren, Qinyan Dai, Xiaotian Feng, Peng Cao, Ruowei Yang Molecular Biology and Evolution, Volume 38, Issue 11, November 2021, Pages 4908–4917, doi.org/10.1093/molbev/msab216Abstract Although Uzbekistan and Central Asia are known for the well-studied Bronze Age civilization of the Bactria–Margiana Archaeological Complex (BMAC), the lesser-known Iron Age was also a dynamic period that resulted in increased interaction and admixture among different cultures from this region. To broaden our understanding of events that impacted the demography and population structure of this region, we generated 27 genome-wide single-nucleotide polymorphism capture data sets of Late Iron Age individuals around the Historical Kushan time period (∼2100–1500 BP) from three sites in South Uzbekistan. Overall, Bronze Age ancestry persists into the Iron Age in Uzbekistan, with no major replacements of populations with Steppe-related ancestry. However, these individuals suggest diverse ancestries related to Iranian farmers, Anatolian farmers, and Steppe herders, with a small amount of West European Hunter Gatherer, East Asian, and South Asian Hunter Gatherer ancestry as well. Genetic affinity toward the Late Bronze Age Steppe herders and a higher Steppe-related ancestry than that found in BMAC populations suggest an increased mobility and interaction of individuals from the Northern Steppe in a Southward direction. In addition, a decrease of Iranian and an increase of Anatolian farmer-like ancestry in Uzbekistan Iron Age individuals were observed compared with the BMAC populations from Uzbekistan. Thus, despite continuity from the Bronze Age, increased admixture played a major role in the shift from the Bronze to the Iron Age in southern Uzbekistan. This mixed ancestry is also observed in other parts of the Steppe and Central Asia, suggesting more widespread admixture among local populations. Introduction Uzbekistan, in Central Asia, includes diverse populations with both East and West Eurasian ancestries (Irwin et al. 2010). During the Bronze Age (BA), the emergence of pastoral economies led to an increased mobility and the development of many settlements across Central Asia (Kohl 2007; Frachetti 2012; Spengler 2015). Many such pastoral settlements from the Mid-BA (∼2100–1800 BCE), particularly in and around South Uzbekistan, possess a material culture associated with the Bactria–Margiana Archaeological Complex (BMAC) (Kohl 2007; Cunliffe 2015). In particular, BA Uzbekistan populations show similar ancestry profiles to those found in other BMAC settlements around the Amu Dariya River, suggesting interactions and connections between them (Narasimhan et al. 2019). The BMAC populations were previously shown to be primarily a mixture of Iranian (∼60–65%) and Anatolian (∼20–25%) farmer ancestries (Narasimhan et al. 2019). Some BMAC individuals were found to have high Yamnaya/Steppe-related ancestry, suggesting this ancestry began appearing in Central Asia by around ∼4100 BP (Narasimhan et al. 2019). Later, in the Mid and Late BA, communities residing in the Bactrian region of Uzbekistan showed higher Steppe-related ancestry compared with the Early BA populations, suggesting an increased influence from the Steppe herders in Uzbekistan (Narasimhan et al. 2019). In Central Asia, the transition from the BA to the Iron Age (IA) toward the middle of the second millennium BCE is characterized by major shifts in material culture with increased mobility and interaction (Cunliffe 2015). During the Early IA, South-Central Asia contained many important settlements, such as Bactria, located in the Amu-Darya basin, Margiana, located near the delta of the Murghab River, Sgodia, around the Zeravshan River basin, and Fergana valley, near the Tianshan mountain range (Lhuillier and Mashkour 2017). During the beginning of the IA (∼1500–1400 BCE) a new culture, “Yaz,” appeared in South-Central Asia characterized by a shift in material culture that included a distinct hand-made pottery and different funerary practices, which lacked prominent graves like the Kurgan-style burial graves observed in the Steppe region (Sarmiento and Lhuillier 2011; Lhuillier and Mashkour 2017). Later, during the middle IA, which extended from ∼1000 BCE until the conquest of Central Asia by the Achaemenids, the Persian Empire saw the development of large settlements, advancement in Iron metallurgy technology and with a new wheel-made pottery style by the “Yaz II” culture, followed by the late IA establishment of the Achaemenid Empire (∼550–330 BCE), also sometimes referred to as “Yaz III” (Kuhrt 2001; Wu et al. 2015; Lhuillier and Mashkour 2017). Later, this region was under the control of the Greco-Bactrian-Kingdom (∼250–125 BCE) (Abazov 2008), during which major upheavals occurred, with migrations and settlements by Indo-European tribes such as the Sakas from North Asia and later by the Yuezhi people from East Asia, who went on to establish the Kushan Empire (∼1st century CE) (Abdullaev 2007; Abazov 2008; Cunliffe 2015; Lhuillier and Mashkour 2017). In general, in archaeological and genetic studies, the impact of Steppe-related culture and ancestry increased from the BA to IA in Central Asia (Cunliffe 2015; Narasimhan et al. 2019), thus genetically, a question remains of how Steppe-related migration into this region from the BA impacted the late IA populations in Central Asia. To address this question, and the general lack of post-BA genomic data from Uzbekistan, we sequenced 27 ancient samples from the Late IA encompassing the Kushan time period (Uz_IA) (∼2100–1500 BP), from three sites in South Uzbekistan: Rabat, Serkharakat, and Dehkan (fig. 1A and supplementary table S1, Supplementary Material online).
|
|