Home / News / Genome


May 21, 2023May 21, 2023

Communications Biology volume 6, Article number: 319 (2023) Cite this article

9271 Accesses

66 Altmetric

Metrics details

Despite the localisation of the southern Caucasus at the outskirt of the Fertile Crescent, the Neolithisation process started there only at the beginning of the sixth millennium with the Shomutepe-Shulaveri culture of yet unclear origins. We present here genomic data for three new individuals from Mentesh Tepe in Azerbaijan, dating back to the beginnings of the Shomutepe-Shulaveri culture. We evidence that two juveniles, buried embracing each other, were brothers. We show that the Mentesh Tepe Neolithic population is the product of a recent gene flow between the Anatolian farmer-related population and the Caucasus/Iranian population, demonstrating that population admixture was at the core of the development of agriculture in the South Caucasus. By comparing Bronze Age individuals from the South Caucasus with Neolithic individuals from the same region, including Mentesh Tepe, we evidence that gene flows between Pontic Steppe populations and Mentesh Tepe-related groups contributed to the makeup of the Late Bronze Age and modern Caucasian populations. Our results show that the high cultural diversity during the Neolithic period of the South Caucasus deserves close genetic analysis.

In the Near East, the Neolithic way of life emerged between 9000 and 7000 BCE. Several centres of Neolithisation have been identified, such as the Levant or Southern China, from which the agropastoral way of life diffused to other regions. The mechanisms of this diffusion have attracted tremendous attention for the last few decades. In some places, the Neolithic gained ground through the acculturation of local hunter-gatherers (for instance, in Anatolia or Iran);1 but in most regions (Europe, South-East Asia), farmer populations spread, and assimilation processes took place with a degree of admixture (for a review, see ref. 2).

The mechanism of Neolithisation in the South Caucasus, a region located between the Black and Caspian Seas on the southern slope of the Greater Caucasus Mountains, remains poorly understood. Mesolithic sites are known at Damjili Cave, unit 5 (Western Azerbaijan)3, Kmlo-2 Rock Shelter (Western Armenia) and Kotias Klde Cave (Western Georgia). Paleogenetic analyses of human bones excavated from the Kotias Klde Cave showed a genetic continuity with earlier Upper Palaeolithic (post-Last Glacial Maximum/LGM) sites4 but a discontinuity with pre-LGM individuals5. Their genetic ancestry shares, to a certain extent, a common origin with ancient Iranian populations6 and differs from that of Anatolian and Levant hunter-gatherer groups1,7, demonstrating a high genetic differentiation at this time between geographically close populations8.

The first settlements attributed to the Early Neolithic period belong to an aceramic culture, evidenced in several places in Central Georgia, as at Nagutni, in Western Georgia, as at Paluri9 and in Western Azerbaijan at Damjili Cave, unit 43. However, evidence of agriculture and herding remains scarce, suggesting that these sites represent a transitional phase between the Mesolithic and the Neolithic.

In this context, the Shomutepe-Shulaveri culture (SSC) is the most ancient Caucasus culture with a complete Neolithic package9,10,11. Found in several clusters of settlements in the northern foothills of the lesser Caucasus, the SSC is characterised by circular mud-brick houses, domestic animals and cereals, handmade pottery, sometimes with incised and relief decoration, and obsidian and bone industries. Variants, such as the Aratashen/Aknashen culture9,12 (Ararat Plain) and other Neolithic contemporaneous cultures like the Kültepe13 Culture (Nakhchivan region) are also found in the South Caucasus. The slightly later Kamiltepe culture (Mil steppe culture), which probably includes the site of Polutepe, differs by its architecture, the use of flint tools instead of obsidian, and pottery-painted patterns that are rather related to Northern Iran and the Zagros14,15.

The origins of the SSC are still discussed10,11. Due to the rapid transition from the aceramic stage to the SSC, population continuity during the Neolithisation process is possible3,12. However, several cultural and biological features are nonlocal. Domesticated animals, such as cattle, pigs or goats, originate from Eastern Anatolia and from the Zagros moutains16,17,18. Similarly, the glume wheat and barley recovered in the SSC sites have been domesticated elsewhere in the Middle East19, even if the ancestors of the naked wheat Aegilops tauschii are found in the Caucasus and may have been local20. The material culture and architecture evidence technical transfers with neighbouring regions such as Southern Anatolia, Pre-Halafian and Halafian culture of Northern Mesopotamia and Zagros11,21. Taken together, these data suggest a strong cultural connection, and maybe a degree of admixture, with other groups from the Fertile Crescent. Indeed, the genome-wide data for one individual coming from the same collective grave at Mentesh Tepe as the samples we analyse here, and confirmed by that of other individuals from Polutepe in the Mughan steppe or from Aknashen and Masis Blur in Armenia22 already showed that southern Caucasian groups are part of a cline that connected Eastern Anatolian and Zagros populations and evidence a gene flow that began around 6500 years BC23.

One of the oldest sites of the SSC, Mentesh Tepe, is located in the Tovuz district of western Azerbaijan and has been excavated between 2007 and 2015 (Fig. 1a). Several occupations were revealed, the earliest dating back to the Neolithic SSC period24,25. At this time, the botanical assemblage is dominated by cereals, especially barley, naked wheat and emmer26, a common association during the Neolithic Southern Caucasus. Animal remains consist largely of domesticated ones (ovicaprines, cattle, pigs, dogs), and wild animals are rare27. The diet of the Neolithic individuals from Mentesh Tepe relied mainly on C3-plants, such as wheat, barley and lentils, with some evidence of freshwater fishes; the consumption of animal proteins varies between individuals28. The pottery differs from classical SSC sites by being vegetal-tempered, a characteristic shared with Kamiltepe or the first occupation of the Nakhchivan site of Kültepe29. As observed in many SSC sites, houses are circular and made of mudbricks, with or without the addition of straw or other organic material21. Two SSC occupation phases are represented, separated in some places by a thick layer of ashes.

a Localisation of Mentesh Tepe and other Neolithic and Protohistoric places that provided ancient genomes. Map tiles by Stamen Design, under CC BY 3.0. Data by OpenStreetMap, under ODbL. b Structure 342 and localisation of each individual. c Photograph of MT23 and MT26 (credit L. Pecqueur).

In a context where Neolithic burials are rare, Mentesh Tepe is exceptional for the discovery of a collective burial containing around 30 individuals (Fig. 1b), which is associated with the end of the first phase of frequentation of the site. Archaeoanthropological analysis has shown that it was a complex funerary gesture with mostly simultaneous deposits and, in contrast, some successive deposits which permitted manipulations on not completely decomposed bodies30. The number of individuals in the burial, as well as their sex and age bias, suggest a dramatic event such as an epidemic, a famine, or a sudden episode, but no trace of violence has been evidenced on the bones.

There is no specific orientation or position of the corpses, but some intentional arrangements are visible. The most striking is formed by two juveniles embracing each other (Fig. 1c). Such an arrangement is rare, but other examples have been found in Neolithic and Protohistoric times, such as in Diyarbakir (Turkey, 6100 BC) or Valdaro (Italy, 3000 BC). Double burials are often considered as a lover's embrace, but arguments for this explanation are often elusive31.

To better understand the origin of the Shomu-Shulaveri population and the structuration of this community, we performed paleogenetic studies of some of the individuals found in the collective burial. The genetic data obtained are then compared to those of another individual from the structure already published and to contemporaneous southwestern Asian genomes.

From the 30 individuals from Structure 342 of the site, we sampled 23 petrous bones but could only obtain genome-wide data for one female (Individual 7, later called MT7) and two males (Individual 23 and Individual 26, MT23 and MT26, respectively), with coverage ranging from 0.1 to 0.3X, a number of SNPs hit on the 1240k from 61,151 to 205,055 and ancient DNA damages (Supplementary Fig. 1). For four other individuals, shallow sequencing and estimation of Rx and Ry only allowed determining biological sex (Supplementary Data 1): two were males and two were females. There is no discrepancy compared to the already published sex determinations made from adult pelvic measurement30, but genetic data allowed a robust determination of the sex for some juveniles.

We merged our genome-wide data with the Human Origins dataset (HO-dataset)1, as well as with 3529 previously unrelated published ancient genomes (Supplementary Data 2). To decipher the genetic relations between the new Mentesh Tepe individuals and other ancient populations from the Caucasus, Anatolia, the Near East, and the Middle East, we performed: (1) a PCA32 on the modern dataset on which the ancient genomes have been projected and (2) an unsupervised ADMIXTURE33 analysis with the HO and the ancient dataset (Fig. 2a, b, Supplementary Fig. 2). The PCA shows that the Mentesh individuals overlap with some other previously published Neolithic or Chalcolithic individuals from the South Caucasus22,23, but the individual from Aknashen falls a bit closer to CHG than the main neolithic cluster (Fig. 2a), and fall intermediate between the Iran Neolithic cluster and the Neolithic Anatolian Farmer group. The ADMIXTURE analysis suggests that the Mentesh individuals carry three main components: i.e., ca. 30% Neolithic Iran (Iran_N; green), 15% Levant Neolithic (PPN; pale rose) and 55% blue and pink components shared with Anatolian or European Neolithic populations (Fig. 2b). Both analyses show that the new Mentesh Tepe individuals present similar profile as that already published from individual MTT001 (despite different sequencing strategy) and from other Neolithic (Polutepe, Azerbaijan, Masis Blur and Aknashen, Armenia) or Chalcolithic (Alkhantepe, Azerbaijan) sites in the South Caucasus without true significant variations in percentage (Kruskal–Wallis test, p-value = 0.25), even though the Aknashen individual has the highest IranN/CHG percentage. They are also very similar to the Chalcolithic and Bronze Age Anatolian populations from Arslantepe and with other Chalcolithic and Bronze Age Anatolian populations but are quite distant from the Late Neolithic Anatolian Tell Kurdu individuals (five individuals identified as Tell_Kurdu_EC by the original publication, whom we refer as TellKurdu_LN) who instead clusters with the Neolithic Anatolian populations. Due to the lack of genomic data from South-eastern Mesopotamia, only the ancestry of groups that lived in North-eastern Mesopotamia during the PPN period and the Late Bronze Age could be considered. We observe that the South Caucasus population displays a different profile from the PPN Mesopotamian one, as it has more Anatolian affinity, but that it presents a profile similar to the LBA individual from North-eastern Mesopotamia (Nemrik9_LBA).

a Principal Component Analysis calculated on 1390 published present-day individuals (in grey) with ancient individuals projected onto it. The new Mentesh Tepe samples are indicated by filled dark pink circles and are labelled. Zoom panel allows to label all the Caucasus Neolithic individuals. b Unsupervised ADMIXTURE analysis, K = 9 for a selection of ancient populations. The Mentesh Tepe individuals appear like other Neolithic and Chalcolithic individuals from the South Caucasus (lowlands_N and LateC). CHG, Caucasus Hunter-Gatherer; Meso, Mesolithic; Levant_EP, Levant Epipaleolithic; C or Chalco, Chalcolithic; EN, Early Neolithic; Neo, Neolithic; BA, Bronze Age. c D-statistics D(Mbuti, Y; Z, Mentesh Tepe). D-statistics tests whether either Z (or Mentesh Tepe) has excess affinity with Y and becomes negative (or positive). Values that deviate from 0 by the ±2 Standard Error (SE)—displayed by the error bars—are represented with a filled circle. d Ancestry proportions of the Neolithic group on autosomes calculated with qpAdm. Ancestry proportions are plotted with error bars representing ±1 SE. POT = Polutepe, MB = Masis Blur, Akn = Aknashen.

To formally test for the genetic affinity with earlier populations in Western Eurasia observed in the PCA in our samples, we performed D-statistics of the form D(Mbuti, Y; Z, MT). The statistic deviates significantly from zero if the pair of Anatolian/Caucasian/Mesopotamian groups (Z) and Mentesh (MT) do not have the same genetic relation to Western European populations (Y). We confirm that the new Mentesh samples are very similar to MTT001 as almost no D-statistics of the form D(Mbuti, Z; MTT001, MT23/MT7/MT26) significantly deviate from 0, showing a high level of homogeneity in this site (Supplementary Fig. 3).

When performing the same D-statistic (Fig. 2c) with MT7, MT23 and MTT001 reunited within a single Mentesh Tepe group (but excluding MT26, as he is related to MT23—cf. infra), we observe that they share more alleles with the Early European Farmers (Greece_N, Serbia_EN), Anatolian farmers, Anatolian Epipaleolithic individual (Anatolia EP) and Neolithic Levant (PPN) than with the Caucasus hunter-gatherers (CHG from the Satsurblia and Kotias Klde caves). They also deviate from the Neolithic Anatolian populations (Barcın site in North Anatolia, Tell Kurdu in South-eastern Anatolia and Tepecik Çiftlik in Central Anatolia) by sharing more alleles with CHG and Iran_N (Ganj Dareh in Zagros mountains) but sharing fewer alleles with Western European Hunter-Gatherers (Iberia_HG, Loschbourg, Goyet_Q2, Latvia_HG), Early European Farmers, Anatolia EP, and with the Neolithic Levant. We also note that they do not differ from the Neolithic individuals from Polutepe and from a Late Chalcolithic individual from Alkhantepe, both sites being further east in the Mughan plain of Azerbaijan (Supplementary Fig. 3). They also do not diverge from the Neolithic Armenian groups (Aknashen and Masis Blur). In the end, in regard to the Eurasian diversity, they all form a homogeneous cluster representing a Neolithic ancestry in the South Caucasus, from east to west, sharing a recent ancestry with Neolithic groups found in South-eastern Anatolia.

Using qpAdm34, we could successfully model the Mentesh Tepe group and each of the individuals as a two-way admixture of TellKurdu_LN, representing South-eastern Anatolian ancestry, and Iran_N (p = 0.20, 37% Iran_N source) (Fig. 2d). Replacing TellKurdu_LN by Barcın_N (p = 0.03–0.75, 33–45% Iran_N source) and replacing Iran_N by CHG (p = 0.01–0.45, 36–49% CHG source) also provide acceptable models for Mentesh as a group but not always for all the samples when modelled individually. We also modelled Mentesh as a mixture between a Neolithic Anatolian group (Tell Kurdu or Barcın) and a North Mesopotamian group (Nemrik9 or Mardin). Eventually, we tested 3-populations models and obtained acceptable models for Mentesh, involving CHG (15–25%), an Anatolian population (Barcın 40% or Tell Kurdu 45%) and a Neolithic North-eastern Mesopotamian population (Nemrik9 29–45%). The same analysis was performed on the other Neolithic South-Caucasus groups (Polutepe, Aknashen and Masis Blur), and we obtained models with a p > 0.05, involving the same kind of combinations as for Mentesh. The amount of Anatolian ancestry in these three populations is not significatively different than that in Mentesh Tepe (Wilcoxon test, p-value > 0.05 for all three), even though the Aknashen individual presents the highest CHG percentage. (Supplementary Fig. 4).

Eventually, we used DATES35 to estimate the date of the admixture event. To account for the low number of individuals representing the Neolithic Ganj Dareh in Iran and the poor coverage, we added other individuals from Neolithic sites in Iran (Seh Gabi, Tepe Abdul Hossein, Wezmeh). We found that the admixture between the Anatolian source and the Iranian source only took place 15 ± 5 generations before Mentesh Tepe occupancy (Z = 2.5 and nrmsd = 0.348). With 28 years per generation, this dates the admixture event around 6300 BC (Phase 1 occupancy: 5880 BC). However, the decay of ancestry covariance estimated by DATES for Mentesh poorly fits with the data (Supplementary Fig. 5).

Having established a general profile for the Neolithic South-Caucasus, we explored the transitions in the genetic structure of Bronze Age populations from the South-Caucasus, including Kura-Araxes individuals from Kalavan-1, the Talin necropolis and the tombs of Kaps in Armenia1,36 as well as all the individuals recently published by Lazaridis et al.22. Sadly, no ancient DNA could be retrieved from the Mentesh Tepe Chalcolithic levels. With almost all D-stats of the form D(Mbuti, Mentesh; Caucasus BA, Anatolia BA/Caucasus BA) being null, we do not see any preferential gene flow from Mentesh Tepe into one Bronze Age population from the South Caucasus or Anatolia (Supplementary Fig. 6).

The PCA shows that all the Bronze Age individuals from Armenia plot together and are shifted toward the Steppe cluster. In the ADMIXTURE analysis, they all exhibit a red component, absent in the Neolithic Mentesh Tepe individuals but maximised in Steppe populations and present, also, in CHG individuals. Interestingly, individuals from Chalcolithic Armenia (from Areni-1 cave, four of whom are directly dated by C14) do carry this Steppe/CHG component, whereas a Chalcolithic individual from Alkhantepe in Azerbaijan does not. D-statistics of the form D(Mbuti, Steppe Eneolithic; Mentesh Tepe, South Caucasus Bronze Age) are almost all significantly positive (Z-score: +2.1 to +5.8), highlighting a gene flow to the South Caucasus from the Steppes or from a population linked to CHG after the Neolithic period. This result can be interpreted with two different hypotheses: either a Neolithic population from North-Caucasus or an ancestral population from the South Caucasus carrying a small proportion of Steppe/CHG ancestry replaced the local Mentesh-like Neolithic population in South-Caucasus, or a population from the steppe north of the Caucasus migrated south and admixed with the local population.

To test for these hypotheses, we used qpAdm with the rotating method37 and modelled the Chalcolithic and Bronze Age populations found in Armenia. The only fitting model for Areni-1 cave (Chalcolithic Armenia) is an admixture between 25% Steppe and 75% Mentesh (p-value = 0.02) (Supplementary Fig. 7). During the Bronze Age, we observe an increase in Steppe contribution from the Early Bronze Age Kura-Araxes (0–10% Steppe contribution) to the Middle and Late Bronze Age individuals (around 40% Steppe contribution). This increase could be linked to a wave of migration from the north during the Bronze Age, or to a continuous admixture between Steppe and South Caucasus populations, maybe through North-Caucasus groups, as the latter are genetically close to the South Caucasus population during the Maikop period/Bronze Age36. For the Early Bronze Age populations of Armenia, we also note that the best models (p = 0.08 for Talin and Karavan and p = 0.51 for Kaps) are involving CHGs instead of a Steppe population. This suggests that the admixture at the base of the Kura-Araxes ancestry occurred between an unsampled population from the Caucasus with a profile more similar to the Caucasus Hunter-Gatherers and a Mentesh-like population from the Late Chalcolithic period. Thus, it is unlikely that Kura-Araxes populations had yet received a significant gene flow from the Steppe at this period. Models involving the only Middle Chalcolithic individual from the North Caucasus (from the Unakozovskaya cave)36 available to date did fit but not as well as the other models (0.01 < p-value < 0.05). Though we note that for Kura-Araxes the models involving the North Caucasus individuals were more successful than for the LBA populations, no nested model involving 100% of North Caucasus Neolithic was detected. In other words, the Bronze Age or Chalcolithic individuals found in Armenia cannot be modelled as 100% North-Caucasus Neolithic. This result suggests that the admixture scenario is more likely than the migration one, given the individuals sequenced for now.

To establish the kinship structure and social organisation of the four sequenced Neolithic Mentesh Tepe individuals, we first look at uniparental markers. Three different mitochondrial haplogroups were determined: U7, K1 and U1a1. Both MT23 and MT26 shared the same mitochondrial haplogroup, as well as the same Y chromosome haplogroup (J2b).

Next, genetic relatedness analyses were performed to further investigate potential close familial ties within individuals of structure 342. To this intent and given the low coverage of our samples, kinship analyses were carried out using READ38.

Although READ can infer relatedness from extremely low-coverage data, a limitation of this method is that it generally requires a cohort of input individuals. It is based on the assumption that most of them are genetically unrelated to compute a sound estimate of the genetic diversity found within our population of interest. To overcome this issue, we extracted the genotype of the Neolithic MTT001, along with 21 previously published Anatolian individuals from the Late Chalcolithic/Early Bronze-Age Arslantepe (ART) and Tell Kurdu from both Neolithic and Early Chalcolithic period23 (Fig. 2a). Here, the intent is to use these individuals as a proxy population, to calculate the median of the average pairwise distances between individuals (P̅0), while minimising the risk of introducing biases arising from a mismatch between the genetic diversity found within burial 342 and that of our surrogate population.

Using this approach, we detected elevated levels of genetic relatedness between individuals MT23–MT26, with a non-normalised average pairwise SNP mismatch ratio of P0 0.239. This value equates to P̅0 734 ± 0.024 when normalised over the median of ART individuals, which can be confidently translated as first-degree individuals (Fig. 3, Supplementary Fig. 8 and Supplementary Data 5), and with a high number of overlapping SNPs (n = 17864) which is well over READ's satisfactory threshold of around 1500 SNPs. These results were subsequently confirmed using TKGWV239, an alternative kinship estimation method which provided us with the same findings (Supplementary Fig. 9).

Sorted normalised average P0 values for all pairwise comparisons between Mentesh Tepe individuals (blue markers) and Arslan Tepe Late Chalcolithic/Early Bronze-Age Arslantepe individuals (orange boxplot), using the READ estimation method. Marker error bars indicate 95% confidence intervals for the observed average P0 estimation. Coloured areas represent the thresholds at which relatedness orders are considered by READ. Boxplot highlights the distribution of the observed normalised P0 of the Arslan Tepe individuals used for the normalisation step of READ, with each point corresponding to a given pair of proxy individuals (jittered along the horizontal axis). Lower and upper fences respectively represent the minimum and maximum sample point, excluding observations located below or above 1.5 the interquartile range. Boxspan: 25th–75th percentile.

However, close genetic ties could not be detected regarding the other pairs involving individuals MT7 and MTT001. The pairs MT23–MT7 and MT26–MT7 were both estimated with a normalised value neighbouring the observed median of ART individuals (P̅0 0.952 ± 0.037 and P̅0 0.963 ± 0.045, respectively, Supplementary Data 5). However, READ's inability to detect any order of kinship greater than the second degree and the high uncertainty in these latter results should be considered when interpreting them, as is demonstrated by the high amplitude of the associated confidence intervals and given the low number of overlapping SNPs between these individuals and MT7 (n = 8674 and n = 5561 for MT23–MT7 and MT26–MT7, respectively).

Our genetic kinship analysis results regarding pair MT23–MT26 unequivocally indicate that these two individuals share a first order of relatedness. This observation, combined with their close age of death, their shared mitochondrial haplogroups, and the context in which the remains of these two individuals were discovered within structure 342, is strongly consistent with the hypothesis that they are, in fact, siblings. It should be noted that the degree of relationship estimated for this pair remained consistent throughout subsequent attempts involving different surrogate populations, including Late Neolithic individuals of Tell Kurdu (Supplementary Fig. 8).

Even after removing MT26, as he is the brother of MT23, the other samples, MT23, MTT001, and MT7, show strong genetic homogeneity despite the different sequencing approaches between MTT001 and the others. However, no nuclear or mitochondrial genetic relationships have been observed between these individuals. This suggests that Mentesh Tepe represents a highly homogeneous population or an extended family uniting people related by more than 3° or more. This observation is not surprising considering the high number of individuals in the structure and the fact that they do not represent the whole community: indeed, adult males are underrepresented, and the bodies have been buried by people aware of the family relationships between them (as exemplified by the staging of MT23 and MT26).

On a regional scale, a certain degree of homogeneity is observed. The Polutepe group, though belonging to a different Neolithic culture, is highly similar to that of Mentesh Tepe, as is the Masis Blur Neolithic group, belonging to the Aknashen culture, a variant of the Shomu-Shulaveri culture. While the Aknashen individual is more Caucasian-like than the other South Caucasus Neolithic individuals, this does not constitute a significant difference, showing that a coherent culture displays a small genetic heterogeneity. Moreover, this period is characterised by several other archaeological cultures, such as Kültepe or Kamiltepe, that have not been studied yet through paleogenetic studies. As the Shomutepe-Shulaveri is already not completely homogenous, it would not be surprising that the population linked to these cultures are genetically different. In addition, the lack of archaeological background for Polutepe prevents us from drawing any conclusions about the associated culture. However, we do show here that the Aknashen culture presents a good degree of genetic homogeneity with the Shomu-Shulaveri individuals but is not fully parallel to the overall cultural unity Aratashen-Shomu-Shulaveri culture seen by the archaeologists.

The genetic analyses performed on the Mentesh Tepe samples allow a better understanding of the diffusion of the Neolithic in this area. Based on the material culture, two models were considered: after a very short transitional phase of an aceramic culture, the Neolithic package appeared with the SSC, either from acculturation or from the migration of human groups from Western Asia9. The genetic data from Mentesh Tepe shows that this SSC site shared one ancestry probably inherited from the Mesolithic Caucasian Hunter-Gatherers (or the Neolithic Iranian Farmers), one from the early Anatolian farmers and one from North-eastern Mesopotamia. The genetic proximity to southern populations could be an indication of the role that Halaf and Hassuna communities played in the formation of the SSC, as already identified through the material culture. Due to the poor DNA preservation, genetic data from the core of Mesopotamia remains scarce, and the genetic affinity of Halaf groups can only be determined by their proximity to South-eastern Anatolia and North-eastern Mesopotamia. Interestingly, Tell Kurdu, which has known a strong influence on the Halafian culture, appears as the best proxy for the ancestors of Mentesh Tepe individuals in the qpAdm analyses. The other part of Mentesh Tepe ancestry is best modelled by Iranian Neolithic groups from Ganj Dareh40. This genetic component is also shared with Caucasian and Iranian Hunter-Gatherer, suggesting that a large population inhabited the southern fringes of the Caucasus Mountains and coasts of the Caspian Sea at the end of the Neolithic period. The currently available data could not discriminate if this ancestry was local (for instance, from the Early Neolithic groups evidenced at Chokh, Kmlo-2 or Damjili Cave) or came from Iranian farmers, thus leading to a complete replacement of the more ancient groups.

Like Skourtanioti et al.23. analyses, which date the admixture event to 6500 BC, we obtain an admixture event around the middle of the 7th millennium, even if their work used two individuals from two different sites with diverse backgrounds, dates and possibly genetic origin, while we focused our analyses on a homogeneous population. However, due to the high number of missing data in the populations used to calculate the date of admixture, the fit is not perfect. Mentesh Tepe appears as a site that closely followed the process of migration and admixture of Anatolian farmers and of North-eastern Mesopotamian populations with the South Caucasus earlier Hunter-Gatherers’ descendants, which was contemporaneous with the establishment of the Shomu-Shulaveri culture and is coherent with the ancient dates obtained for this site.

The genetic data from Mentesh Tepe also helps understanding better the subsequent periods, as previous studies exploring the Bronze Age in the South Caucasus only had available the Chalcolithic genome-wide data from the North Caucasus1,36.

Indeed, we find that the North Caucasus Chalcolithic individual used previously36 as a proxy for the Chalcolithic and Neolithic Caucasus shows a genetic profile that differs from the Mentesh Tepe one, with a slight amount of Eastern Hunter-Gatherer component likely due to a gene flow from the adjacent steppe. The comparison was not directly feasible in the former study due to the lack of Neolithic data from the southern slope of the Great Caucasus, which emphasises the need for a thorough sampling in ancient DNA studies, both in space and time.

Thanks to the new data, we show that, after the end of the Neolithic, a genetic component related to the CHG reappeared in the South Caucasus and that, later in the Bronze Age, a genetic component related to the Pontic steppe groups arrived. Chalcolithic and Late Bronze Age groups in the South Caucasus can be modelled as a mixture of Mentesh Tepe-related populations and the Steppe population. Interestingly, Steppe ancestry is not spread homogeneously during the Chalcolithic period, whereas Areni-1 individuals display a high amount of steppe ancestry; none is found in the Alkhantepe genome. Such diversification in the population's structures could mirror the one observed in cultural features. The funerary practices were also more diversified during the Chalcolithic, and while the dead at Alkhantepe are still buried in pits according to ancient practices, the first kurgans are observed in Kavtiskhevi (Georgia) and Soyuq Bulaq (Azerbaijan)41,42. Part of the Areni-1 pottery as well as the funerary architecture and pottery from Soyuq Bulaq connect them with the Maikop tradition43. The pre-Kura-Araxes South-Caucasus thus appears as a mosaic of populations of different ancestry, with the intrusion of more northern populations. Nevertheless, the precise dating of one individual from Areni 1 could be uncertain (cf. supplementary text).

Compared to the Mentesh Tepe Neolithic individuals, an archaic ancestry associated with CHG re-emerges during the Early Bronze Age/Kura-Araxes period. We also observed some Steppe ancestry in two Kura-Araxes groups (Kaps, Shengavit), but at the lowest proportion than during the Late Bronze Age. At the end of the Kura-Araxes period, more frequent connections with the Steppes are also evidenced in the burial rituals44. This gene flow from the north mirrors the earlier (Maikop period) increase of the South Caucasus ancestries in the populations living on the northern slopes of the Great Caucasus after the end of the Neolithic, thus showing that the mountains act as a bridge rather than a frontier.

The analysis of ancient DNA from the Mentesh Tepe Neolithic funerary pit reveals an uneven genetic homogeneity inside this Shomu-Shulaveri culture archaeological site. After the dramatic event that caused their death, the two brothers were buried, embracing each other, by people knowing their close relationship. The Shomu-Shulaveri culture took place between two intense population events: a population expansion of Neolithic groups from South-eastern Anatolia and Northern Mesopotamia at the beginning of the Neolithic and diffusion of steppe populations starting with the Middle of the Bronze Age period and following the resurgence of the CHG ancestry during the Kura-Araxes period. Even if the genetic analyses presented here help to better understand the timing and direction of these events, a more detailed analysis of other samples from the South Caucasus is called for, as this area is characterised by a high cultural diversity during the Neolithic period.

We obtained ethical authorisation from the Institute of Archaeology and Ethnography of Azerbaijan to study samples from Mentesh Tepe. Extraction was attempted on the 30 samples found in Structure 342 of the site, except MTT001 (Individual 1), which has already been published elsewhere23. Petrous bone fragments from 23 among 30 individuals of the SSC structure S342 have been selected (Supplementary Data 1) for ancient DNA analysis based on overall bone preservation.

All pre-amplification steps were carried out in the clean room dedicated to ancient DNA of the Paleogenomic and Molecular genetics platform, set in the Musée de l’Homme (Paris). As in ref. 45, ancient DNA extraction was performed using a protocol adapted from ref. 46. Briefly, 50–200 mg of petrous bone was powdered by drilling and incubated in 1 ml of lysis buffer (0.45 M EDTA, 10 mM Tris-HCl (pH 8.0), 0.1% SDS, 65 mM DTT, and 0.5 mg/ml proteinase K), at 37 °C, for 14 h. After centrifugation, 1 ml of supernatant was recovered and purified with 13 ml of binding buffer (5 M GuHCl, 40% 2-propanol 0.05% Tween20, 90 mM sodium acetate 2 M, 1x Phenol Red). The mixture was then transferred to a High Pure Extender Assembly column (Roche High Pure Viral Nucleic Acid Large Volume Kit) and centrifuged. Then, the column was washed using manufacturer's recommendations (briefly, 500 µl of inhibitor removal buffer, centrifugation, twice 450 µl of wash buffer). DNA was eluted in 100 µl of elution buffer. First, we tested the presence of human ancient DNA through PCR by targeting a portion of the mitochondrial genome using the same protocol as in ref. 45. Only samples with positive mitochondrial DNA amplification were used for library preparation.

The ancient DNA extract was converted to a TruSeq Nano Illumina library using the manufacturer's protocol with slight modifications that account for the ancient DNA damage. First, DNA was not fragmented for obvious reasons; only 25 µl of ancient DNA extract was used; after end-repair, the libraries were purified using a MinElute column (Qiagen ©); the adaptor mix was gently pre-heated before adding ligase enzyme; libraries were amplified using 10–12 PCR cycles and purified on a MinElute column (Qiagen ©). Analyses on a LabChip ® GX provided an estimated size distribution of fragments with a peak length of 150–250 bp. Due to the limited preservation of genetic material in the site, reliable ancient DNA was found only in seven individuals, and genome-wide data were generated only for three of them.

For the three selected samples (endogenous DNA content >1%), genomic capture was performed using the myBaits Expert Whole Genome Enrichment (WGE) kit (Arbor Biosciences), and the manufacturer's instructions were followed. Baits were formed from the genomic DNA of three individuals of different (African, European and Asian) ancestries.

Aliquots of the amplified libraries were first pooled and sequenced on MiSeq instrument (2 × 75 bp) on the IGenSeq platform. Captured libraries were sequenced on a NextSeq 500 (2 × 75 bp) instrument.

Considering that our kinship analysis protocol makes use of previously published ancient individuals as a surrogate population, a separate data-preprocessing workflow mirroring that of ref. 23 was devised to make our results as comparable as possible with these previous analyses. All of our newly sequenced samples were thus specifically preprocessed through EAGER-v1.92.3747. From the raw shotgun sequencing data, sequencing adaptors were thus clipped, and paired-end reads were merged using AdapterRemoval-v2.2.1a48, all the while discarding any read shorter than 30 bp. Reads were then aligned against the 1000Genomes-phase2 reference genome (hs37d5) using BWA-v0.7.17-r118849, with a quality threshold for read trimming of q30. Samples were then merged across sequencing runs for each individual, using samtools merge v1.6.050. PCR duplicates were marked and removed using DeDup v0.12.147 prior to and following the previously described merging. PMD-damage statistics and frequency of nucleotide misincorporation were then estimated until the 50th nucleotide (starting from the 3′ and 5′ ends) using MapDamage v2.0.651 (Supplementary Fig. S1). Quality scores within the bam files were then rescaled accordingly with the latter software, using default parameters. We extracted reads overlapping known variants from the 1240K dataset from the rescaled bam file using samtools mpileup, while disabling per-Base Alignment Quality (-B flag) and skipping any alignment carrying a base and/or mapping Phred-quality score which was lower than 30 (−q and −Q flags). Random pseudo-haploid variant calling was then performed using SequenceTools PileupCaller-v1.5.0 (Download:

To evaluate contamination, we used AuthentiCT52, and for the male samples, we calculated contamination using the X chromosome53 with ANGSD54. (Supplementary Data 1).

We selected 3529 published ancient human genomes from Eurasia (Supplementary Data 2) from the Palaeolithic to Middle Age, whom DNA sequencing data were generated with whole genome shotgun or hybridisation capture technics from the merge dataset v42.4 available at and in two recent paper about the South Caucasus and Anatolia22,23. We retained the nonrelated individuals with more than 10,000 SNPs hit on the 1240k panel. We analysed ancient data with 1587 Eurasian individuals from the Human Origin dataset. We merged our genome-wide data to this dataset using mergeit from the EIGENSOFT package suite55, and we haploidised them by randomly selecting one allele per position. The final merge includes 597,573 SNPs for 5116 individuals, and we call it the HO-dataset. For the analysis requiring more SNPs, we use a merge without the modern individuals from the HO dataset, covering 1,233,013 SNPs, which we call the 1240k dataset.

We ran PCA with smartpca55 using the HO-dataset, on 1390 European and Middle Eastern individuals, and we projected all the ancient samples. We used the default parameters with lsqproject: YES, and numoutlieriter: 0 settings.

We computed ADMIXTURE33 analysis using the HO dataset, where all populations were downsampled to a maximum of 20 individuals. We performed the analysis on 1266 Eurasian modern individuals, including East Asian, and 2526 ancient samples from the first dataset on a subset of 365,075 SNPs pruned for linkage disequilibrium (by using PLINK --indep-pairwise 200 25 0.4 function). We ran 20 replicate ADMIXTURE analyses for K between 2 and 14, from which we kept the most likely.

We performed D-statistics on the 1240k dataset using the qpDstat program of the ADMIXTOOLS package56 with Mbuti as outgroup for both statistics. We computed all the D-statistics for each Mentesh Tepe individual and for the group Mentesh Tepe formed by MT23, MT7 and MTT001. We computed D-statistics of the form D(Mbuti, Y; Z, MT), with Z being Palaeolithic or Neolithic groups from Anatolia (Barcın site in North Anatolia, TellKurdu site from South-eastern Anatolia and Tepecik site from Central Anatolia), Northern Mesopotamia (Nemrik 9 and Shanidar sites from Irak, Mardin site from Turkey) or Caucasus (CHG, Polutepe, Armenia Neolithic, Armenia Chalcolithic, North Caucasus Chalcolithic) and Y ancient Palaeolithic, Mesolithic or Neolithic populations from Western Eurasia to test the affinity of Mentesh Tepe to them when compared to an older contemporaneous group from the same region. We also performed D-stat of the form D(Mbuti, MT; Anatolia or Caucasus BA, Caucasus BA) to test the affinity of Mentesh Tepe with the populations that succeeded it in the region. These last results are shown in the supplementary information (Supplementary Fig. 3). All the results of the D-statistics analyses, including the data used to generate the figures, are in Supplementary Data 3.

We performed rotating qpAdm analysis from ADMIXTOOLS package56, using the 1240k dataset, to model the ancestry of Mentesh Tepe and Polutepe populations. We used Mbuti.DG, Ami.DG, Mixe.DG, Russia_Kostenki14.SG, Russia_MA1_HG.SG, EEHG, Italy_North_Villabruna_HG, Natufian as reference populations. Prior to the analysis, we checked if the reference populations could discriminate well between the source populations by computing qpWave.

We tested a rotating group with Iran_GanjDareh_N, CHG, Russia_Caucasus_Eneolithic, PPN, Barcın_N, Anatolia_TellKurdu_LN IRQ_Nemrik9_PPN and TUR_SE_Mardin_PPN to assert the best model of a two-way admixture for each of our individuals. We then run the same qpAdm analysis for the Mentesh group form by MT23, MT7 and MTT001 and all the South-Caucasus Neolithic (SCN) individuals in one group. We also tested for 3-way admixture models; we could obtain model without nested 2-way models for Mentesh. All the results of the qpAdm modelling, including the data used to generate the figures, are in Supplementary Data 4.

We used DATES v7531835 to estimate the time of admixture events between Iran_N and Anatolia_N in the Mentesh group. To convert the estimated admixture date in a generation into years, we assumed 29 years per generation13. The standard errors of DATES estimates come from the weighted block jackknife with "binsize: 0.001", "maxdis: 1", "runmode: 1", "mincount: 1", "lovalfit: 0.45" as parameters as in the example file at We tested three different configurations to account for the poor coverage of the ancient individuals: the first test was conducted with only ancient Iranian Neolithic individuals as a source and Barcın_N or Tell Kurdu_LN as the second source. The second, with modern populations from Armenia, Georgia and the Russian Caucasus added in the Iranian source. A third and final test was made with Barcın_N and Tell Kurdu_LN individuals put together as a source. We performed these three configurations for the Mentesh individuals grouped together and for all the Neolithic South-Caucasus individuals as a group.

Y chromosome and mitochondrial haplogroups were determined by comparison of vcf files to databases. For mitochondrial DNA, a comparison was performed on Phylotree 17, using Haplogrep v2.2;57 for Y chromosome, the comparison on ISOGG version July 11, 2020, was done manually.

Using PLINK-v1.90b6.1658,59, we first filtered out from our 1240K dataset anything but biallelic autosomal SNP positions carrying a minor allele frequency greater or equal to 5% and subsequently extracted the genotypes of individuals MT23, MT26, MT7 and MTT001, along with 21 unrelated Late Chalcolithic and Early Bronze-Age individuals from Arslantepe (ART). The resulting dataset was then transposed to the tped/tfam format and used as input for READ. Next, we devised a two-pass estimation approach by running the READ method twice: A first run was therefore performed using only the Arslantepe individuals and with default parameters to compute a normalisation value (normalisation method: median; nonoverlapping window-size: 1000000). A second run was then carried out against the Mentesh Tepe individuals alone, using the previously obtained median P̄0 as a predefined normalisation value. We re-applied this protocol a second time using six individuals from the Late Neolithic and Early Chalcolithic site of Tell Kurdu to substantiate our initial results (Supplementary Fig. 8).

We also used TKGWV2 to substantiate our initial results (Supplementary Fig. 9). As per recommendations from the authors, the final BAM files from all Mentesh individuals were downsampled to around 1.8M reads beforehand39. We then performed the analysis using the 22M nonfixed biallelic variant set provided by Fernandes et al.39. To assess confidence in the raw relatedness coefficients given by the method, we obtained posterior probabilities by simulating distribution curves for each tested pair of individuals (n = 6000 simulation replicates).

The dataset used for PCA and ADMIXURE includes 597,573 SNPs for 5116 individuals, and the dataset used for Dstats, qpAdm and DATES is covering 1,233,013 SNPs for 5,116 individuals. When an analysis targets Mentesh as a group, it means three individuals grouped together (MT7, MT23 and MT001). ADMIXTURE analysis has been replicated 20 times. For Kinship analysis, READ analysis has been replicated twice: one with 21 individuals from Arslantepe and one with six individuals from Tell Kurdu.

Further information on research design is available in the Nature Portfolio Reporting Summary linked to this article.

The dataset generated for this study (bam file) can be found in the European Nucleotide Archive under accession number PRJEB54894. The data used to generate Fig. 2c are in Supplementary Data 3, for Fig. 2d in Supplementary Data 4 and for Fig. 3 in Supplementary Data 5.

Code to reproduce the analysis and the figure of this paper can be obtained here:, or here:

Lazaridis, I. P. et al. Genomic insights into the origin of farming in the ancient Near East. Nature 536, 419–424 (2016).

Article CAS PubMed PubMed Central Google Scholar

Racimo, F., Sikora, M. Vander Linden, M., Schroeder, H. & Lalueza-Fox, C. Beyond broad strokes: sociocultural insights from the study of ancient genomes. Nat. Rev. Genet. 21, 355–366 (2020).

Article CAS PubMed Google Scholar

Nishiaki, Y. et al. The Mesolithic-Neolithic interface in the Southern Caucasus: 2016–2017 excavations at Damjili Cave, West Azerbaijan. Archaeol. Res. Asia 19, 100140 (2019).

Article Google Scholar

Jones, E. R. et al. Upper Palaeolithic genomes reveal deep roots of modern Eurasians. Nat. Commun. 6, 1–8 (2015).

Article CAS Google Scholar

Gelabert, P. et al. Genome-scale sequencing and analysis of human, wolf, and bison DNA from 25,000-year-old sediment. Curr. Biol. 31, 3564–3574.e9 (2021).

Article CAS PubMed PubMed Central Google Scholar

Marchi, N. et al. The genomic origins of the world's first farmers. Cell 185, 1842–1859.e18 (2022).

Article CAS PubMed PubMed Central Google Scholar

Feldman, M. et al. Late Pleistocene human genome suggests a local origin for the first farmers of central Anatolia. Nat. Commun. 10, 1218 (2019).

Article PubMed PubMed Central Google Scholar

Lazaridis, I. et al. Ancient DNA from Mesopotamia suggests distinct Pre-Pottery and Pottery Neolithic migrations into Anatolia. Science 377, 982–987 (2022).

Article CAS PubMed PubMed Central Google Scholar

Chataigner, C., Badalyan, R. & Arimura, M. The Neolithic of the Caucasus. Oxford Handbooks Online. (Oxford University Press, 2014).

Nishiaki, Y., Guliyev, F. & Kadowaki, S. (eds). Hacı Elamxanlı Tepe. The Archaeological Investigations of an Early Neolithic Settlement in West Azerbaijan. (Bibliotheca neolithica Asiae meridionalis et occidentalis ex oriente, 2021).

Badalyan, R., Chataigner, C. & Harutyunyan, A. (eds) The Neolithic Settlement of Aknashen (Ararat valley, Armenia). (Archeopress, 2022).

Nishiaki, Y., Guliyev, F. & Kadowaki, S. Chronological contexts of the earliest Pottery Neolithic in the South Caucasus: radiocarbon dates for Göytepe and Hacı Elamxanlı Tepe, Azerbaijan. Am. J. Archaeol. 119, 279 (2015).

Article Google Scholar

Marro, C., Bakhshaliyev, V., Berthon, R. & Thomalsky, J. New light on the Late Prehistory of the South Caucasus: data from the recent excavation campaigns at Kültepe I in Nakhchivan, Azerbaijan (2012–2018). Paléorient 45, 81–113 (2019).

Article Google Scholar

Lyonnet, B. et al. Ancient Kura 2010–2011: The first two seasons of joint field work in the Southern Caucasus. Archaeologische Mitteilungen aus Iran und Turan. 44, 1–190 (Deutsches Archäologisches Institut, Eurasien-Abteilung, 2012).

Helwing, B. et al. The Kura Projects New Research on the Later Prehistory of the Southern Caucasus (Dietrich Reimer, 2017).

MacHugh, D. E., Larson, G. & Orlando, L. Taming the past: ancient DNA and the study of animal domestication. Annu. Rev. Anim. Biosci. 5, 329–351 (2017).

Article CAS PubMed Google Scholar

Verdugo, M. P. et al. Ancient cattle genomics, origins, and rapid turnover in the Fertile Crescent. Science 365, 173–176 (2019).

Article CAS PubMed Google Scholar

Kadowaki, S., Ohnishi, K., Arai, S., Guliyev, F. & Nishiaki, Y. Mitochondrial DNA analysis of ancient domestic goats in the Southern Caucasus: a preliminary result from Neolithic settlements at Göytepe and Hacı Elamxanlı Tepe. Int. J. Osteoarchaeol. 27, 245–260 (2017).

Article Google Scholar

Akashi, C., Tanno, K., Guliyev, F. & Nishiaki, Y. Neolithisation processes of the South Caucasus: as viewed from macro-botanical analyses at Hacı Elamxanlı Tepe, West Azerbaijan. Paléorient 44, 75–89 (2018).

Google Scholar

Mizuno, N., Yamasaki, M., Matsuoka, Y., Kawahara, T. & Takumi, S. Population structure of wild wheat D-genome progenitor Aegilops tauschii Coss.: implications for intraspecific lineage diversification and evolution of common wheat. Mol. Ecol. 19, 999–1013 (2010).

Article PubMed Google Scholar

Baudouin, E. Rethinking architectural techniques of the Southern Caucasus in the 6th millennium BC: a re-examination of former data and new insights. Paléorient 45, 115–150 (2019).

Article Google Scholar

Lazaridis, I. et al. The genetic history of the Southern Arc: a bridge between West Asia and Europe. Science 377, eabm4247 (2022).

Article CAS PubMed Google Scholar

Skourtanioti, E. et al. Genomic history of Neolithic to Bronze Age Anatolia, Northern Levant, and Southern Caucasus. Cell 181, 1158–1175.e28 (2020).

Article CAS PubMed Google Scholar

Lyonnet, B. et al. Mentesh Tepe, an early settlement of the Shomu-Shulaveri Culture in Azerbaijan. Quat. Int. 395, 170–183 (2016).

Article Google Scholar

Lyonnet, B. et al. In The Kura Projects. New Research on the Later Prehistory of the Southern Caucasus. Mentesh Tepe (Azerbaijan), a preliminary report on the 2012–2014 excavations, 125–140 (D. Reimer, 2017).

Decaix, A. et al. Vegetation and plant exploitation at Mentesh Tepe (Azerbaijan), 6th–3rd millennium BC initial results of the archaeobotanical study. Quat. Int. 395, 19–30 (2016).

Article Google Scholar

Benecke, N. In The Kura Projects. New Research on the Later Prehistory of the Southern Caucasus (eds. Helwing, B. et al.) Exploitation of animal resources in neolithic settlements of the Kura region (South Caucasia), 357–369 (D. Reimer, 2017).

Herrscher, E. et al. Dietary inferences through stable isotope analysis at the Neolithic and Bronze Age in the southern Caucasus (sixth to first millenium BC, Azerbaijan): from environmental adaptation to social impacts. Am. J. Phys. Anthropol. 167, 856–875 (2018).

Article PubMed Google Scholar

Lyonnet, B. Mentesh Tepe 2012–2014. In The Kura Projects. New Research on the Later Prehistory of the Southern Caucasus. (eds. Helwing, B. et al.) The Pottery, 141–151 (D. Reimer, 2017).

Pecqueur, L. & Jovenet, E. In The Kura Projects New Research on the Later Prehistory of the Southern Caucasus. (eds. Helwing, B. et al.) La sépulture 342 de Mentesh Tepe (Azerbaïdjan): un exemple de chaîne opératoire funéraire complexe au Néolithique. Etude préliminaire, 163–178 (D. Reimer, 2017).

Geller, P. L. The Bioarchaeology of Socio-Sexual Lives. Queering Common Sense About Sex, Gender, and Sexuality. (Springer, 2017).

Price, A. L. et al. Principal components analysis corrects for stratification in genome-wide association studies. Nat. Genet. 38, 904–909 (2006).

Article CAS PubMed Google Scholar

Alexander, D. H., Novembre, J. & Lange, K. Fast model-based estimation of ancestry in unrelated individuals. Genome Res. 19, 1655–1664 (2009).

Article CAS PubMed PubMed Central Google Scholar

Haak, W. et al. Massive migration from the steppe was a source for Indo-European languages in Europe. Nature 522, 207–211 (2015).

Article CAS PubMed PubMed Central Google Scholar

Narasimhan, V. M. et al. The formation of human populations in South and Central Asia. Science 365, eaat7487 (2019).

Article CAS PubMed PubMed Central Google Scholar

Wang, C.-C. et al. Ancient human genome-wide data from a 3000-year interval in the Caucasus corresponds with eco-geographic regions. Nat. Commun. 10, 590 (2019).

Article CAS PubMed PubMed Central Google Scholar

Harney, É., Patterson, N., Reich, D. & Wakeley, J. Assessing the performance of qpAdm: a statistical tool for studying population admixture. Genetics 217, iyaa045 (2021).

Article PubMed PubMed Central Google Scholar

Kuhn, J. M. M., Jakobsson, M. & Günther, T. Estimating genetic kin relationships in prehistoric populations. PLoS ONE 13, e0195491 (2018).

Article Google Scholar

Fernandes, D. M., Cheronet, O., Gelabert, P. & Pinhasi, R. TKGWV2: an ancient DNA relatedness pipeline for ultra-low coverage whole genome shotgun data. Sci. Rep. 11, 21262 (2021).

Article CAS PubMed PubMed Central Google Scholar

Gallego-Llorente, M. et al. The genetics of an early Neolithic pastoralist from the Zagros, Iran. Sci. Rep. 6, 31326 (2016).

Article CAS PubMed PubMed Central Google Scholar

Poulmarc’h, M. & Le. Mort, F. Diversification of the funerary practices in the Southern Caucasus from the Neolithic to the Chalcolithic. Quat. Int. 395, 184–193 (2016).

Article Google Scholar

Lyonnet, B. et al. Late Chalcolithic Kurgans in Transcaucasia. The cemetery of Soyuq Bulaq (Azerbaijan). Archäologische Mitteilungen aus Iran und Turan 40, 27–44 (2008).

Google Scholar

Areshian, G. E. et al. The chalcolithic of the Near East and south-eastern Europe: discoveries and new perspectives from the cave complex Areni-1, Armenia. Antiquity 86, 115–130 (2012).

Article Google Scholar

Lyonnet, B. In V. Sarianidi Memorial Volume, Transaction of Margiana Archaeological Expedition. (eds. Dubova, N. et al.) A grave with a wooden wagon in Transcaucasia (Azerbaijan). Its relations with Central Asia,191–198 (2016).

Guarino-Vignon, P. et al. Genetic analysis of a Bronze Age individual from Ulug-depe (Turkmenistan). Front. Genet. 13, 884612 (2022).

Article CAS PubMed PubMed Central Google Scholar

Dabney, J. et al. Complete mitochondrial genome sequence of a Middle Pleistocene cave bear reconstructed from ultrashort DNA fragments. Proc. Natl Acad. Sci. USA 110, 15758–15763 (2013).

Article CAS PubMed PubMed Central Google Scholar

Peltzer, A. et al. EAGER: efficient ancient genome reconstruction. Genome Biol. 17, 60 (2016).

Article PubMed PubMed Central Google Scholar

Schubert, M., Lindgreen, S. & Orlando, L. AdapterRemoval v2: rapid adapter trimming, identification, and read merging. BMC Res. Notes 9, 88 (2016).

Article PubMed PubMed Central Google Scholar

Li, H. & Durbin, R. Fast and accurate short read alignment with Burrows-Wheeler transform. Bioinformatics 25, 1754–1760 (2009).

Article CAS PubMed PubMed Central Google Scholar

Li, H. et al. The Sequence Alignment/Map format and SAMtools. Bioinformatics 25, 2078–2079 (2009).

Article PubMed PubMed Central Google Scholar

Jónsson, H., Ginolhac, A., Schubert, M., Johnson, P. L. F. & Orlando, L. mapDamage2.0: fast approximate Bayesian estimates of ancient DNA damage parameters. Bioinformatics 29, 1682–1684 (2013).

Article PubMed PubMed Central Google Scholar

Peyrégne, S. & Peter, B. M. AuthentiCT: a model of ancient DNA damage to estimate the proportion of present-day DNA contamination. Genome Biol. 21, 246 (2020).

Article PubMed PubMed Central Google Scholar

Rasmussen, M. et al. An Aboriginal Australian genome reveals separate human dispersals into Asia. Science 334, 94–98 (2011).

Article CAS PubMed PubMed Central Google Scholar

Korneliussen, T. S., Albrechtsen, A. & Nielsen, R. ANGSD: analysis of next generation sequencing data. BMC Bioinformatics 15, 356 (2014).

Article PubMed PubMed Central Google Scholar

Patterson, N., Price, A. L. & Reich, D. Population structure and Eigenanalysis. PLoS Genet. 2, e190 (2006).

Article PubMed PubMed Central Google Scholar

Patterson, N. et al. Ancient admixture in human history. Genetics 192, 1065–1093 (2012).

Article PubMed PubMed Central Google Scholar

Kloss-Brandstätter, A. et al. HaploGrep: a fast and reliable algorithm for automatic classification of mitochondrial DNA haplogroups. Hum. Mutat. 32, 25–32 (2011).

Article PubMed Google Scholar

Purcell, S. et al. PLINK: a tool set for whole-genome association and population-based linkage analyses. Am. J. Hum. Genet. 81, 559–575 (2007).

Article CAS PubMed PubMed Central Google Scholar

Chang, C. C. et al. Second-generation PLINK: rising to the challenge of larger and richer datasets. Gigascience 4, 7 (2015).

Article PubMed PubMed Central Google Scholar

Download references

The authors express their gratitude to the French-Azerbaijani excavation team, who contributed to the sampling. The excavation in Mentesh Tepe is partly supported by the Ministry of Foreign and European Affairs and the ANR projects "Ancient Kura" and "Kura in Motion", directed by B. Lyonnet and B. Helwing. This project has been funded by the ANR project "Kura in Motion" and the Sorbonne University Emergence grant (PaleOxus). P.G.V. was funded by a PhD grant given by ENS de Lyon. We thank the Paleogenomic and molecular genetics platform (P2GM) of the MNHN (Paris) as well as the IGenSeq platform of the Institut du Cerveau et de la Moelle (Paris). We thank Nina Marchi for the comments and critics and Lasse Vinner for help in setting up a first version of the DNA extraction protocole.

UMR7206 Éco-Anthropologie (EA), CNRS, Muséum National d’Histoire Naturelle, Université Paris-Cité, Paris, France

Perle Guarino-Vignon, Maël Lefeuvre, Amélie Chimènes, Laure Pecqueur & Céline Bon

UMR5288 CAGT, CNRS, Université Paul Sabatier, Toulouse, France

Perle Guarino-Vignon

Department of Archaeology, University of York, BioArCh, Environment Building Wentworth Way Heslington, York, YO10 5NG, UK

Aurore Monnereau

Head of the Science Fund and the Museum Department of the Institute of Archeology, Ethnography and Anthropology, Azerbaijan National Academy of Science, Baku, Azerbaijan

Farhad Guliyev

Institut National de Recherches Archéologiques Préventives, Centre-île de France, Paris, France

Laure Pecqueur & Elsa Jovenet

UMR7192 PROCLAC « Proche-Orient – Caucase: langues, archéologie, cultures », CNRS, Paris, France

Bertille Lyonnet

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

C.B. and B.L. designed the study. B.L., F.G., E.J. and L.P. provided archaeological material and advised on the archaeological background. A.C. and A.M. performed the laboratory work. P.G.V. and M.L. performed the data analyses, with C.B. providing guidance. P.G.V., M.L., B.L. and C.B. wrote the manuscript with input from all co-authors.

Correspondence to Perle Guarino-Vignon or Céline Bon.

The authors declare no competing interests.

Communications Biology thanks the anonymous reviewers for their contribution to the peer review of this work. Primary handling editor: George Inglis. Peer reviewer reports are available.

Publisher's note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit

Reprints and Permissions

Guarino-Vignon, P., Lefeuvre, M., Chimènes, A. et al. Genome-wide analysis of a collective grave from Mentesh Tepe provides insight into the population structure of early neolithic population in the South Caucasus. Commun Biol 6, 319 (2023).

Download citation

Received: 07 September 2022

Accepted: 08 March 2023

Published: 25 March 2023


Anyone you share the following link with will be able to read this content:

Sorry, a shareable link is not currently available for this article.

Provided by the Springer Nature SharedIt content-sharing initiative

By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.