Interdisciplinary Analyses of Bronze Age Communities from Western Hungary Reveal Complex Population Histories

Abstract In this study, we report 21 ancient shotgun genomes from present-day Western Hungary, from previously understudied Late Copper Age Baden, and Bronze Age Somogyvár–Vinkovci, Kisapostag, and Encrusted Pottery archeological cultures (3,530–1,620 cal Bce). Our results indicate the presence of high steppe ancestry in the Somogyvár–Vinkovci culture. They were then replaced by the Kisapostag group, who exhibit an outstandingly high (up to ∼47%) Mesolithic hunter–gatherer ancestry, despite this component being thought to be highly diluted by the time of the Early Bronze Age. The Kisapostag population contributed the genetic basis for the succeeding community of the Encrusted Pottery culture. We also found an elevated hunter–gatherer component in a local Baden culture–associated individual, but no connections were proven to the Bronze Age individuals. The hunter–gatherer ancestry in Kisapostag is likely derived from two main sources, one from a Funnelbeaker or Globular Amphora culture–related population and one from a previously unrecognized source in Eastern Europe. We show that this ancestry not only appeared in various groups in Bronze Age Central Europe but also made contributions to Baltic populations. The social structure of Kisapostag and Encrusted Pottery cultures is patrilocal, similarly to most contemporaneous groups. Furthermore, we developed new methods and method standards for computational analyses of ancient DNA, implemented to our newly developed and freely available bioinformatic package. By analyzing clinical traits, we found carriers of aneuploidy and inheritable genetic diseases. Finally, based on genetic and anthropological data, we present here the first female facial reconstruction from the Bronze Age Carpathian Basin.


Introduction
Several studies have addressed major population historical events in prehistoric Europe regarding pre-Neolithic hunter-gatherers (HGs) (Fu et al. 2016;Mittnik et al. 2018;Rivollat et al. 2020), their assimilation to early European farmers (EEF) during the Neolithic era (Haak et al. 2015;Lipson et al. 2017;Mathieson et al. 2018;Mittnik et al. 2018), and the appearance, expansion, and admixture of steppe ancestry between the Eneolithic/Late Copper Age and the dawn of Early Bronze Age (Allentoft et al. 2015;Haak et al. 2015;Olalde et al. 2018;Papac et al. 2021).Though these large studies are vital for understanding the roots of the European gene pool, there are currently only a few studies that uncover regional interaction or social stratification using biological relatedness analyses (Schroeder et al. 2019;Freilich et al. 2021;Žegarac et al. 2021), especially in Hungary during the Copper (∼4,500-2,700 BCE) and Bronze (∼2,700-800 BCE) Ages.A number of cultural transformations occurred in the Carpathian Basin, often attributed to population changes and genetic influxes.The observed shifts in ancestry compositions sparsely covered intensive or direct European HG introgressions into EEF or steppe ancestry groups (Lipson et al. 2017) (besides one such case from today's Romania; González-Fortes et al. 2017), despite the known HG presence in the region at the beginning of the Neolithics (Gamba et al. 2014;Mathieson et al. 2018).In contrast, ancient populations from other parts of Europe, such as Scandinavia (Skoglund et al. 2014;Günther et al. 2018;Mittnik et al. 2018;Malmström et al. 2019), today's Poland (Fernandes et al. 2018) or Iberia (Olalde et al. 2019), show much higher and much later introgression of HG ancestry.Later on, at the beginning of the third millennium BCE, the appearance of steppe-related ancestry shaped the regional genetic landscape extensively, founding the modern-day European genetic makeup (Gamba et al. 2014;Lazaridis et al. 2014;Allentoft et al. 2015;Brandt et al. 2015;Haak et al. 2015;Mathieson et al. 2018;Olalde et al. 2018;Linderholm et al. 2020).
Besides monitoring population events, archeogenetics opens a new window to study health qualities of ancient populations that may lead to a better understanding of the background of recent genetic-related diseases.Studies aimed at uncovering variants under selective pressure in Homo sapiens populations, such as lactase persistence (LCT) or human leukocyte antigen (HLA) genes and pigmentation markers, are beginning to thrive (Mathieson et al. 2015;Childebayeva et al. 2022;Evershed et al. 2022;Lazaridis et al. 2022).However, variants for rare genetic diseases or aneuploidies are sparsely checked on ancient data sets, except for a few cases, such as the study of the Suontaka grave (Moilanen et al. 2022).
Our study aimed to make a transect analysis on a single site from understudied archeological assemblages.This was combined with population genetic analyses, isotope analyses, and phenotype (pigmentation) and clinical variant analyses.Moreover, we present a series of bioinformatic tools for biological relatedness, ploidy, and variant analyses implemented in a new bioinformatic package for ancient DNA analysis.We analyzed archeological finds from the Balatonkeresztúr-Réti-dűlő site in Western Hungary (Transdanubia), where-among others-Bronze Age features and human remains were found during roadwork in 2003.Three Bronze Age archeological horizons were distinguished based on 14 C dates: the Somogyvár-Vinkovci culture (∼2,500-2,200 BCE, n = 1), Kisapostag culture (∼2,200-1,900 BCE, n = 11), and the Encrusted Pottery culture (∼1,900-1,450 BCE, n = 8) that are referred to as Bk-I, II, and III phases in this study, respectively (table 1; supplementary section 1, Supplementary Material online).All three cultural horizons have only a limited number of inhumation remains: this study presents the first validated Somogyvár-Vinkovci culture-associated individual from Hungary, whereas the Kisapostag and Encrusted Pottery cultures have been mainly characterized by cremation burials so far.The archeological origin of the Kisapostag culture is enigmatic, and multiple theories have been proposed to explain its possible connections.One theory is that the pottery decoration technique observed in Kisapostag originated either from Corded Ware in the Middle Dnieper region (Ukraine) or epi-Corded Ware groups (from the northern Carpathians), for example, Chłopice-Veselé (Slovakia).The latter option is also supported by inhumation practices and burial positions (Bóna 1961(Bóna , 1975;;Bándi 1984;Hajdu 2010;Szabó 2010;Kiss 2020).However, local development of communities with eastern (Makó-Kosihy-Čaka) or southern (Somogyvár-Vinkovci) origins, as well as western and southwestern connections (with the Litzenkeramik or Guntramsdorf-Drassburg group in eastern Austria, Slovenia, and western Croatia), have also been raised in the archeological literature (Bóna 1992;Črešnar 2010;Kiss 2015).Additionally, Bell Beaker influence has been proposed based on the craniometry data (so-called Glockenbecher or brachycranic skull type ;Mozsolics 1941;Zoffmann 2007Zoffmann , 2008)).Nevertheless, the archeological connection between Kisapostag and Encrusted Pottery cultures are strong and well established (Kiss 2012).
In order to provide an additional data point to our analysis of population ancestry of the region, a Late Copper Age individual from the Baden culture (∼3,600-2,800 BCE), excavated ∼30 km away from Balatonkeresztúr, was added to the data set.Our data highlight not only detailed population events in a microregion but also reveal hidden processes that formed the genetic landscape of East-Central Europe at the beginning of the Bronze Age.

Results
We shotgun sequenced genomes of 21 individuals yielding between 0.008x and 2.1x average genomic coverage.We also sequenced reads of a capture set consisting 3,000 nuclear single nucleotide polymorphisms (SNPs; see Materials and Methods) and whole mitochondrial DNAs (mtDNAs) of all individuals.The shotgun-and the capture-sequenced samples ultimately resulted in an average ∼144k SNPs/individual using the 1240k SNP panel for genotype calling (Mathieson et al. 2015) (see Materials and Methods; supplementary tables S4 and S7, Supplementary Material online).We utilized short tandem repeat (STR) analysis of the Y chromosome to ascertain direct paternal connections (supplementary table S3, Supplementary Material online).The Y-STR profiles resulted in comparable haplogroup predictions with the NGS (next-generation sequencing) Y-SNP data and were authenticated via repeated reactions (see Materials and Methods).Furthermore, using all known biological and archeological details, we reconstructed the face of individual S13 (from phase Bk-II) (see supplementary section 4, Supplementary Material online).

MBE
The bioarcheological analyses included 14 C dating and 87 Sr/ 86 Sr isotope analyses; the latter is routinely used to trace individual mobility (Alt et al. 2014).

Archeological and Anthropological Evaluation of Samples
We included only one Infans II (late childhood) individual (BAD002) from the site of Balatonlelle belonging to the early phase of the Baden culture 370 cal BCE).A majority of the samples came from the site of Balatonkeresztúr-Réti-dűlő.We sampled and sequenced 20 individuals, labeled S1-S45, skipping numbers that did not belong to the Bronze Age horizon or were not suitable for genetic testing.One male individual (S9) that was associated with Bk-I by 14 C data has a very long (ultradolichocran) skull type, which differentiates him from most individuals in Bk-II and Bk-III, who have very short (brachycranic) skull types (Köhler 2014) (table 1).The male dominance (∼78%) in Bk-II and Bk-III suggests distinctive funeral treatment for males and females.Bk-II phase is represented by an Infans II (7-8-year-old male), three juveniles (16-19-year-old males), and seven adults (30+ years old).They are spatially distributed into grave groups of A and B with two additional, separate inhumations (table 1; supplementary fig.S.1.2.1, Supplementary Material online).Most of the burials contained no grave goods except for small copper jewelry in S10 and S13, and some shell fragments in S45.Radiocarbon dates place these inhumations to ∼2,200-1,770 cal BCE (95.4% CI [confidence interval]); however, with Bayesian analysis using the OxCal software (Bronk Ramsey 2009), the timespan of the Bk-II burials can be reduced to ∼2,050-1,940 cal BCE with a 84.4% CI, whereas only two graves (individuals S10 and S11) were possibly slightly earlier (supplementary section 1.8, Supplementary Material online).The lack of infants at Bk-II is similar to other archeological sites as this phenomenon is common in most periods.It can be associated with different skeletal taphonomy for infants or differential burial treatment as compared with the burial of adults (Alt et al. 2014).The reason for the absence of young adults (∼20-30-year-olds) is unknown.Bk-III is represented by a single mass grave (∼1,870-1,620 cal BCE) with skeletal remains of eight people of various ages.This is an unusual burial treatment as most bodies were cremated or buried as single inhumations in that culture. Fo detailed descriptions of the sites and burials, see supplementary section 1, Supplementary Material online.

Uniparental Genetics and Relatedness Analyses
Analysis of uniparental (maternal and paternal) lineages provides rough estimates of genetic composition and is essential to assess biological relatedness and social structure for the studied population.Additionally, we performed phylogenetic analysis using MrBayes software (Ronquist and Huelsenbeck 2003) on mtDNA to see the phylogeographic affinities of the studied individuals.Accordingly, the valid phylogenetic trees show diverse connections in both Bk-II and Bk-III to various modern and ancient populations, but the overlap between the two horizons are limited.For detailed mtDNA analyses, see supplementary section 2.1, Supplementary Material online.
Contrary to the diverse mtDNA makeup, most male individuals in Bk-II and Bk-III belong to the Y chromosome haplogroup I2a-L1229, except for two individuals that -mtDNA and ChrY denote mitochondrial haplogroup and Y chromosome haplogroup, respectively.In the column "Biological relatedness," 1st and 2nd mean the degree of relations.For the feature, grave ID, and details on newly reported 14 C dates, see supplementary table S1, Supplementary Material online.
belong to haplogroup R1b-Z2103 (table 1).Similar phylogeographic analysis to the mtDNA can be performed on the paternal lineages as well using STR markers.Network analysis (supplementary section 2.2, Supplementary Material online) narrowed down regional Y-chromosomal affinities to the North European Plain and indicated continuity between Bk-II and Bk-III.Uniparental makeup shows a patrilocal social structure that is similar to previously reported Bronze Age findings (Mittnik et al. 2019;Schroeder et al. 2019;Žegarac et al. 2021).Results are highly similar to previous observations on Encrusted Pottery culture's population at the Jagodnjak site, Croatia (Freilich et al. 2021)

Genetic Disorders and Pigmentation
Investigating genetic disorders in archeological data sets can be useful to improve our knowledge of the history of health and medicine and also highlights the overall genetic health of past populations.Genetic disorders, if accompanied with severe phenotypic anomalies, could also explain unusual burial practices, for example, as it was described in cases of dwarfism in the Byzantine era (Slon et al. 2013) or in the case of the Suontaka burial (Moilanen et al. 2022).Therefore, we analyzed the ploidy of the autosomes not only for genetic sex determination but to recover possible aneuploidies resulting in serious health-related traits, such as Turner or Down syndrome (O'Connor 2008).For inferring pigmentation of the studied individuals, we used the HIrisPlex (Chaitanya et al. 2018) system supplemented with further variants obtained from SNPedia database (Cariaso and Lennon 2012).Finally, we created an in-house disease panel described further below and in supplementary section 3, Supplementary Material online.

Aneuploidies
Abnormal numbers of chromosomes result in a few wellknown diseases that we tested for thoroughly.We developed a new method called Z-score-Adjusted Coverage (ZAC) for ploidy estimates by using a set of reference genomes.This can estimate ploidy for samples as low as 0.01x average genomic coverage, enabling genetic sex determination and aneuploidy assessment for all new samples (supplementary section 3.1, Supplementary Material online).As a result, we found one individual, S10-the only child burial in Bk-II-with an XYY gonosomal genotype, described as Jacob's syndrome.This syndrome is relatively frequent (∼0.1%) in today's populations.In most of the cases, it remains silent but occasionally comes with a wide scale of symptoms, mainly behavioral disorders (Bardsley et al. 2013).
mtDNA Diseases We examined the clinical significances of the polymorphisms that can be found in mtDNA by using the mitopathotool software on the AmtDB database (Ehler et al. 2019) and found that individual S1 (40+-year-old male from Bk-II) had one of the defining mutations (T14484C, 48x coverage) of Leber's hereditary optic neuropathy (LHON) causing complete vision loss in ∼50% of males between the ages of 20 and 40 and in rare circumstances it is accompanied by other neuropathies (Tońska et al. 2010).

Nuclear Variants with Clinical Significance
We also examined the nuclear genomes to find regions with clinical significance.Since a complete panel for determining disease susceptibility only exists in commercial DNA kits with nonavailable descriptions similar to the 1240k panel, we created a new SNP calling panel focused on various conditions including amyotrophic lateral sclerosis, Alzheimer disease, autism, Crohn's disease, diabetes, lactose intolerance, mental disorders, Parkinson disease, schizophrenia, and ulcerative colitis.For this study, we used a ∼5k set of clinically significant SNPs, which were marked as "pathogenic" or "likely pathogenic" in the ClinVar database (Landrum et al. 2018), by ignoring deletion, duplication, and copy number variants, as well as SNPs with questionable (signed as "reported," "conflicting reports," etc.) contributions to diseases.The exact method of calling variants can be found in supplementary section 3.2.2,Supplementary Material online.Both the tool and the ∼5k set are built in the PAPline package.We also created a bioinformatic tool that rolls up variant information from input data, which is available in the PAPline package (see the PAPline section for details).After running the panel, we excluded low-coverage transitional variants from the final evaluation due to the possible presence of DNA damage.We only made exceptions when skeletal features supported the presence of the low-quality variant or when more than one sample possessed the same allele.(Freilich et al. 2021) both provided a single read hit.We excluded the Jagodnjak group from our analyses due to the lack of uracil-DNA glycosylase (UDG) treatment, meaning that both hits could be false positives; however, individual S15 possesses the distinguishing skeletal features of this disease, increasing the possibility of the actual presence of this allele in the Encrusted Pottery population.Another ambiguous but possible hit is rs121434442 in individual S6.This SNP is the causative factor for hereditary spastic paraplegia (Shribman et al. 2019), a disease mostly recognized by muscle stiffness in the lower limbs causing movement restrictions.Individual S11, father of S6, shows signs of a limb condition that may be linked to this disease.Finally, autism 15 susceptibility signature transversional variant (rs7794745) (Koeda et al. 2015) was present in individuals S6 and S45.Severe bruxism on the upper front teeth of S45 (supplementary fig.S.3.2.2.4,Supplementary Material online) suggests compulsive behavior that occurs frequently among people with autism spectrum disorder (Muthu and Prathibha 2008;Koeda et al. 2015).Whereas this condition itself could also be linked to some profession-related abrasion, the physical features combined with genetic data and the distinguished burial treatment (supplementary fig.S.1.5.11,Supplementary Material online) support the possibility of the actual onset of symptoms.

Pigmentation
According to our results based on a final set of 58 SNPs, the pigmentation patterns are highly different between horizons, as Bk-I mostly possesses variants for light pigmentation, blue eyes, and blonde hair, whereas Bk-II is more similar to populations of Neolithic Europe of darker coloration (Childebayeva et al. 2022;Lazaridis et al. 2022) (fig.1c), although some variants for lighter pigmentation exist within this group too.Members of Bk-III on the other hand show a wide range from dark to light pigmentation tones and even the presence of variants for red hair (supplementary table S5 and section 3.2.1,Supplementary Material online).

Whole Genome Composition and Genetic Ancestry
Balatonkeresztúr Site Samples To get a general overview of the autosomal composition of the individuals, we performed principal component analysis (PCA) with smartpca software (Patterson et al. 2017) based on 590k nuclear SNPs (Mallick and Reich 2023) and ADMIXTURE (Alexander et al. 2009) analyses based on the 1240k SNP set (Mallick and Reich 2023).According to PCA (fig.2a), Bk-I is clearly separated from Bk-II and Bk-III, where Bk-II has a strong shift toward European HGs overlapping with only a fraction of known ancient samples (Mallick and Reich 2023) and Bk-III.Admixture analyses (fig.2b) for assessing genetic components revealed ∼17% HG, ∼40% EEF, and ∼43% steppe ancestry for Bk-I, similar to average Bronze Age Europeans (Gamba et al. 2014;Allentoft et al. 2015;Brandt et al. 2015;Haak et al. 2015;Olalde et al. 2018) (supplementary tables S9 and S12-S16 and sections 5.2, 5.5.2, and 5.6, Supplementary Material online).According to qpAdm (Patterson et al. 2012), Bk-I is most likely the ∼1:2 mixture of a population, represented by a Vučedol culture-associated individual (Croatia_EBA_Vucedol_3, ∼38 ± 4%) and a mostly steppe characteristic source.This steppe source can be best modeled as a Srubnaya/ Alakul culture-related population (Russia_Srubnaya_Alakul.SG, ∼62 ± 4%), in line with archeological observations (Kulcsár 2009).However, this high proportion of steppe ancestry is likely derived from a previously unsampled group in Eastern Europe, maybe in the vicinity of the Baltics (for details, see supplementary section 5.6.2.1, Supplementary Material online).Bk-II comprises a unique makeup of ∼42% HG, ∼41% EEF, and ∼17% steppe ancestries.qpAdm analysis revealed most plausible sources as a Sweden_FBC (Funnelbeaker culture)-related population and Ukraine_EBA with almost equal contributions (supplementary section 5.6.2.2, Supplementary Material online); however, both populations are likely only an approximation for the actual ancestry of Bk-II, which we discuss further below.Bk-III shows a slight shift in ancestry composition from Bk-II with ∼29% HG, ∼46% EEF, and ∼25% steppe ancestries.qpAdm analyses uncovered that the main ancestry component for Bk-III is Bk-II (∼60 ± 8%), whereas "dilution" of Bk-II to Bk-III is mostly driven by contact with various local populations, genetically best represented by later Transdanubian Hungary_LBA or Serbia_Mokrin_EBA_-Maros (Maros culture) groups (supplementary section 5.6.

Genetic Outliers from Previous Studies and the Origin of HG Ancestry in Bk-II
We aimed to investigate further the exact composition of HG ancestry in Bk-II.qpAdm analysis of the basic composition resulted in EEF ∼40 ± 2%, Eastern HG (EHG) ∼39 ± 3%, Western HG (WHG) ∼13 ± 2.7%, and Caucasus HG ∼8 ± 2%, P = 0.0917, in par with Admixture analysis, pointing toward a rather EHG characteristic composition.Next, we performed an f4 test in the form of f4(test HG, Serbia_IronGates_Mesolithic, Bk-II, Mbuti.DG) (Patterson et al. 2012), to see which HG population relates the best with the Bk-II samples.The aim of this analysis was to detect different HG ancestry contributions besides Iron Gates HG, which has knowingly contributed to Neolithic European populations in the study region and witnessed an intermediate composition between the WHG and EHG (Mathieson et al. 2018).Contrary to the Admixture and qpAdm results, this test revealed that Bk-II individuals have excess in HG ancestry mainly from WHG groups or other mixed characteristic HGs (Croatia_Mesolithic, Poland_BKG_o2.SG [Brześć Kujawski Group outlier], or KO1 [Körös culture outlier HG]), but only marginal relations with the EHG (Lithuania_Mesolithic) populations (for detailed results, see supplementary section 5.3, Supplementary Material online).Surprisingly, none of these HG populations with mixed characteristics (and neither Iron Gates) have enough EHG component to explain the ancestry of the Bk-II samples.The f4 test also revealed that the Bk-II and Bk-III populations differed significantly from other EHG characteristic populations, such as HGs from today's Russia, Ukraine, Baltics (younger phase), or Scandinavia, although we can see some weak connections to older (down to sixth millennium BCE) Lithuanian HGs.These results may reflect the population turnover in the sixth millennium BCE in the Baltics (Mittnik et al. 2018), suggesting that this EHG ancestry is related to the Lithuania_ Mesolithic.On the other hand, qpAdm provides negative weights for this component when we model Bk-II as a combination of WHG (Loschbour_WHG), EHG (Lithuania_Mesolithic), EEF (Turkey_N), and Yamnaya (Russia_EBA_Samara_Yamnaya), suggesting that Lithuania_ Mesolithic is not a good proxy for the actual EHG component.
To infer the timing of HG admixture, we used the DATES (Chintalapati et al. 2022) analysis.This test revealed that the HG ancestry in Bk-II resulted from three independent admixture events: one from Iron Gates HG at the beginning of the Neolithic (similar to other populations at that time), one from a WHG characteristic source around the turn of the fourth and third millennium BCE, and an EHG characteristic source around the second half of the third millennium BCE (for details, see supplementary section 5.4, Supplementary Material online).Summarizing these results, we conclude that the EHG characteristic source of the Bk-II individuals does not exist in the current published database.
We were interested in whether other populations carried this peculiar HG ancestry, to see which region it might originate from.To achieve this, we did a literature search to select individuals with high levels of HG ancestry, who Gerber et al. • https://doi.org/10.1093/molbev/msad182MBE were genetic outliers in their cultural or geographical or temporal context, in order to assess whether they are related to the Bk-II group.Selection was based on previous observations and HG ancestry differences within groups using the results of the Admixture analysis (supplementary section 5.2.2, Supplementary Material online).Then, to reveal similar patterns of HG ancestry, we ran f3 statistics in the form of f3(test HG, test population, Mbuti.DG) on all of the groups (obtained from AADR [Mallick and Reich 2023] database, listed in supplementary table S8, Supplementary Material online).Subsequent Euclidean distance-based clustering of f3 values revealed a number of outliers and even whole populations belonging to the same subcluster as Bk-II (fig.3; supplementary section 5. 5, Supplementary Material online).Accordingly, the earliest signs of such HG ancestry appeared among various Neolithic groups from Western Europe (in line with characteristically high WHG ancestry among Megalithic, Globular Amphora, or FBCs' population) and from Eastern Europe (Bulgaria and Ukraine).Individuals with this ancestry predating Bk-II by only a few generations appeared in Czechia, Northern Hungary, Eastern Germany, and Western Poland, indicating that the Kisapostag-associated population probably came to Transdanubia via a northern route, in line with the observations of Freilich et al. (2021).
Many contemporaneous populations to Bk-II and Bk-III from the British Isles to today's Poland, down to today's Serbia have outliers with a Bk-II-like genomic composition, mostly overlapping with known Kisapostag and Encrusted Pottery contact regions (fig.3).Interestingly, at the end of the second millennium BCE, many Baltic groups appear to be highly similar to Bk-II, indicating long-term success of this ancestry outside the Carpathian Basin.Notably, in the vicinity of Prague, many pre-and post-Bk-II outliers appear along with the archeological presence of the Kisapostag culture, including the Tollense group, which also originates from the region of Bohemia according to isotopic evidence (Price et al. 2019), suggesting a local reservoir of the population.Whereas the appearance of Bk-II ancestry in the Baltics could be connected to this reservoir, especially in the light of the mobility of Tollense group, the 14 C date of Lithuania_LN_o around 2,000 BCE suggests that the population was likely prevalent in nearby unsampled regions of Eastern Europe.
Taking into consideration all of the genetic parallels, their dates, and geographic locations, one plausible scenario is that the EHG characteristic core of Bk-II (which ultimately could be best modeled as Ukraine_EBA by composition) moved northward from the region of today's Eastern Romania, Moldavia, or Western Ukraine, subsequently mixed with FBC-or Globular Amphora culture (GAC)-related populations and then split into two groups: one taking a route to Transdanubia and the other moving further North.These results are highly in par with Mittnik et al. (2018), who suggested population replacement at the end of the second millennium BCE in the Baltic region from a nearby, unsampled region by a population of considerably higher steppe, EEF, and WHG ancestry than the prevailing ones; however, further data are needed from Eastern Europe to affirm this hypothesis.

Isotope Analyses
We sampled molars from the burials and measured the ratio of 87 Sr/ 86 Sr isotopes to evaluate whether individuals were born in the area of their burial.The results (fig.4) show that burials from Balatonkeresztúr have strontium isotope values that match local plants and water, which suggests that none of the studied individuals are firstgeneration occupants.It is, however, interesting how the molar M3 values for individuals S15, S16, and S17 (all from the mass grave) differ from the others.Whereas these values are not out of the estimated local isotope range for the area, they could indicate movement within the region during adolescence when those tissues are forming.This movement could have occurred at the same time for these individuals, as the stronger the divergence from the majority, the younger the individual was at the time of death.It is particularly interesting how individual S15 shows the highest divergence from the others, as this individual had severe complications for walking due to hip dysplasia  (Kiss 2012).Red and purple circles represent individuals that are connected to Bk-II individuals by HG ancestry.Also, red circles are preceding, whereas pinks are succeeding or contemporaneous to the Bk-II horizon.For the detailed list, see supplementary table S5.5.1, Supplementary Material online, and for method description, see supplementary section 5.5, Supplementary Material online.(Lipson et al. 2017), BAD002 has a higher HG component (∼34%), and he also lacked steppe-related ancestry.Therefore, on the genomic PCA, BAD002 relates with Iberian and French Neolithic individuals.According to our ancestry estimates, France_MontAime_MLN.SG describes best the BAD002 individual but other Western European sources, such as Spain_EN, are also plausible (for details, see supplementary section 5.6.1,Supplementary Material online).The pigmentation pattern of BAD002 shows resemblance to average Neolithic Europeans.The foreign cultural traits of the boy's jewelry is in line with his outlier genetic composition in the study region (Bondár and Szécsényi-Nagy 2020); therefore, we conclude that this individual exemplifies large-distance migration in the Copper Age, providing research questions for future studies.Finally, further tests (outgroup f3 statistics and qpAdm) excluded the contribution of BAD002 to Bk-II (supplementary section 5, Supplementary Material online).

PAPline
We introduce a newly developed, freely available bioinformatic package, named PAPline (Performing Archaeogenetic Pipeline), written in linux bash, R, and python v3.8.10 programming languages.One can use this package primarily to analyze next-generation sequencing data of archeogenomic samples, supplemented by tools, including ploidy test, MPMR relatedness analysis, and clinical variant test.The standalone tools and the core workflow of PAPline are available at https://github.com/ArchGenIn/papline.In the future, PAPline is aimed to be compared with the EAGER (Peltzer et al. 2016) and the Paleomix (Schubert et al. 2014) pipelines; for a detailed description, visit the github page.

Discussion
The Carpathian Basin was inhabited by the Baden population at the end of the Copper Age.Their genetic composition was represented by an EEF and-compared with the previous Neolithic populations of the region (Lipson et al. 2017)-a slightly increased HG genetic component.Here, we demonstrated that in the early phase of this culture, a Western European group appeared in Transdanubia, diversifying our previous knowledge about the region's Late Copper Age.
The Carpathian Basin experienced the influx of steppe-related genetic ancestry during the Late Copper Age (Mathieson et al. 2018;Olalde et al. 2018).This transformation was already detectable in the Bronze Age genetic picture of the Balatonkeresztúr-Réti-dűlő, where we could examine multiple populations.The earliest Bronze Age horizon Bk-I (representative of the Somogyvár-Vinkovci culture) is best described by the mixture of local (Vučedol) and a high steppe ancestry population from Eastern Europe; however, since this group has data from only one individual, future studies on this archeological culture and region may diversify the composition of the prevailing population.Bk-I was replaced by the Kisapostag culture-associated group of Bk-II, likely around the 23rd-22nd centuries BCE.According to our results, the Bk-II population had an outstandingly high HG genetic ancestry level compared with other Bronze Age groups of the region.This can be traced back to two main sources, one to a WHG and one to an EHG characteristic population, best modeled as FBC/GAC and Ukraine_EBA; however, likely both are only proximate to the actual source, which are yet to be described.The Y chromosome haplogroup I2a-L1229 can be linked to the FBC/GAC component, although this subgroup only appears first in Bk-II and related groups.The calculated admixture dates suggest the presence of a highly EHG characteristic population in Eastern Europe as late as the beginning of the Bronze Age.This EHG component shows the closest resemblance to Lithuanian Mesolithic individuals, but the best proxy for this population is probably missing from the published database, opening research questions for future studies.
Following the formation of the population represented here by Bk-II, it contributed to various populations in Central-Eastern Europe, whose genetic legacy persisted mostly in the region of today's Hungary and Czechia at least until the end of the Bronze Age and even to the end of the first millennium BCE in the Baltic region.This study does not disclaim any of the archeological theories regarding the origin of the Kisapostag culture (Bóna 1975(Bóna , 1992;;Zoffmann 2008;Črešnar 2010;Kiss 2015), as the EHG core of the Kisapostag-associated group fits really well with the Middle Dnieper origin, whereas further adaptation of cultural elements during their arrival and during their occupation in Transdanubia is plausible.The latter idea is further supported by the 87 Sr/ 86 Sr isotope ratio data that show local isotope ratios for both sexes in both Bk-II and Bk-III.These results place back the time of their arrival a few generations, meaning that both local and southern cultural traits could explain the culture's archeological heterogeneity.
Bk-III was the direct descendant of Bk-II, forming cultural (Encrusted Pottery) and genetic continuity for hundreds of years at the studied site.Observable dilution of HG ancestry in Bk-III compared with Bk-II can be Interdisciplinary Analyses of Bronze Age Communities • https://doi.org/10.1093/molbev/msad182MBE connected to continuous female-biased admixture with nearby communities according to our and previous genetic (Freilich et al. 2021) and archeological (Bóna 1975;Vicze 2011) evidence.The observed shift in pigmentation patterns is likely the result of this admixture.
In both periods, the homogeneity of paternal lineages suggests a patrilocal residence system, similarly to previously described social organizations (Schroeder et al. 2019;Freilich et al. 2021).However, 87 Sr/ 86 Sr isotope data show 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.The overlap between outlier parallels of Bk-II/III and archeological contact regions is also noteworthy, as it suggests smaller scale migrations of Kisapostag/Encrusted Pottery individuals or groups along trading networks, such as mobility possibly connected to wandering merchants.
Notably, none but one (mtDNA haplogroup U5a1g) of the uniparental lineages are the same at the haplogroup level with the individuals from the Croatian Encrusted Pottery culture Jagodnjak site, despite high similarities in cultural traits, social structure, and genomic composition of the communities (Freilich et al. 2021).This supports a regional patrilocal, clanlike superfamily structure of Kisapostag and Encrusted Pottery groups.This finding is particularly interesting in light of a strikingly different social structure observed among the 2,100-1,800 BCE Maros culture individuals from Mokrin, Serbia (Žegarac et al. 2021), that shows extensive amounts of admixture related to the Kisapostag/Encrusted Pottery culture.
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.However, in other cemeteries of the culture, for example, Ordacsehi and Bonyhád in Hungary, males, females, and children were buried close to each other, suggesting high variance in burial practices (Somogyi 2004;Hajdu 2010;Szabó 2010).
Although low genomic coverages did not allow fine SNP recovery, we did find evidence for malignant variants within all of the tested groups and undoubtedly showed the presence of LHON and Jacob's syndrome within Bk-II.Whereas it only remains a possibility, the presence of autism risk factor in the CNTNAP2 gene, signs for severe bruxism, and the distinctive burial treatment of individual S45 suggest the actual onset of symptoms.Additionally, the disease panel we created and made freely available can be extended and used in future studies, providing insight into past population health qualities.
Considering the unstructured age and biological relatedness distribution in the mass grave of Bk-III compared with Bk-II, the coetaneous death of eight people at least, the absence of traumatic or ritual events on bones, and the noncremated nature of the burial all signal a sudden tragic event in the Middle Bronze Age period (Encrusted Pottery population), most likely an epidemic, as first suggested based on the anthropological analyses (Köhler 2006).Careful burial positions also suggest that the deceased were buried by their own community.Interestingly, comparative 87 Sr/ 86 Sr isotope analyses on the first and third molars 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 organization.

Isotope Analyses
Radiocarbon dating was performed at the HEKAL AMS C-14 facility of the Institute for Nuclear Research in Debrecen, Hungary (see supplementary section 1.8, FIG. 4. 87 Sr/ 86 Sr isotope data from the Balatonkeresztúr site.Samples were taken from dental enamel (first, second, and third molars, labeled as M1, M2 and M3) 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 87 Sr/ 86 Sr ratio (green diamonds) data collected from the southern portion of Lake Balaton (Alt et al. 2014).For further data, see supplementary section 1.9, Supplementary Material online.

FIG. 1 .
FIG. 1. Biological relatedness of Bk-II and Bk-III, reconstruction of individual S13.(a) Distribution of graves of Bk-II: individuals S1, S2, S5, S6, and S7 belong to grave group "A", whereas individuals S4, S8, S11, and S13 belong to grave group "B".Individuals S10 and S45 (not shown) placed separately from grave groups.(b) Relatedness network of Bk-II individuals, where colors denote the corresponding grave groups in (a).The dashed line between individuals S1 and S2 represents an undirected second-degree relationship.(c) Relatedness of Bk-III projected onto the photo of the mass grave (obj.B-938).Squares denote males, circles females, solid lines first-degree, and dashed lines second-degree relationships.(d ) Forensic facial and burial reconstruction of individual S13.

FIG. 2 .
FIG. 2. Basic genetic composition of the investigated samples.(a) PCA based on 590k SNPs calculated by the smartpca software (Patterson et al. 2017), where Bk-II (marked with green triangles) clearly separated from all other ancient Central-Eastern European populations.(b) The admixture proportions of the BAD002, Bk-I, Bk-II, and Bk-III samples, where the percentage of steppe ancestry is shown with red, EEF with light green, EHG with light blue, and WHG with dark blue color (supervised Admixture analyses).

FIG. 3 .
FIG. 3. Map of East-CentralEurope with sites and genetic parallels of Kisapostag/Encrusted Pottery culture.The map shows the site of Balatonkeresztúr (red star), Kisapostag and Encrusted Pottery culture archeological sites (dark blue circles), and their archeological connections based on pottery and metal finds (light blue circles) after(Kiss 2012).Red and purple circles represent individuals that are connected to Bk-II individuals by HG ancestry.Also, red circles are preceding, whereas pinks are succeeding or contemporaneous to the Bk-II horizon.For the detailed list, see supplementary tableS5.5.1, Supplementary Material online, and for method description, see supplementary section 5.5, Supplementary Material online.

Table 1 .
Summary of the Investigated Samples.

ID Grave group cal BCE date (95.4% CI) Age Sex mtDNA ChrY Biological relatedness Baden BAD002 3,530-3,370 8-9 M K1a4a1 I-M170 Bk-I: Somogyvár-Vinkovci culture S9 2,560-2,290 35-40 M K1a3a R1a-V2670 Bk-II: Kisapostag or Early Encrusted Pottery culture S1 A 2,120-1,880 40+ M V I2a-L1229 2nd to S2 S2 A 2,120-1,880 30-35
. Inference on biological relatedness was based on READ (Monroy Kuhn et al. 2018) and a newly developed method called Modified Pairwise Mismatch Rate (MPMR; supplementary section 2.3, Supplementary Material online).The relatedness network (fig.1b) of Bk-II approximately follows the distribution of individuals in A-B grave groups (fig. 1a), which were likely established along family relationships and chronology.Individuals buried in the Bk-III mass grave only show a few blood relations (fig.1c): a half-brother or uncle-nephew, a father-son, and a dizygotic twin; to our knowledge, the latter is the oldest detection of such relatedness.None of the distant inhumations (S10, S45) show biological relationship to any other individuals up to second degree.For further details, see supplementary section 2 and tables S1-S3, Supplementary Material online.
Gerber et al. • https://doi.org/10.1093/molbev/msad182MBE the Jagodnjak site of the Encrusted Pottery culture for firm conclusions; however, the aim was more of a technical description of such analyses.Here, we summarize only a few mentionable results of the run, but for the detailed discussion, see supplementary section 3.2, Supplementary Material online.Lig4 syndrome is a transitional mutation (rs104894421) induced disease with skeletal abnormalities (Altmann and Gennery 2016), for which individuals S15 and JAG93 from Gerber et al. • https://doi.org/10.1093/molbev/msad182MBE(seealso supplementary section 3.2.2,Supplementary Material online).Moreover, since M3 values show divergence, but their first molar (M1) do not, these individuals likely grew up in the vicinity of the site, spent some years away from it, and then returned to the same place where they died and were buried.