Y-chromosomal lineage determination
We determined Y-chromosome haplogroups in ancient male samples using clean_tree(85). This software implements SAMtools (86) mpileup to call alleles at given genomic coordinates (2,710 SNPs), which we obtained from ISOGG 2013 (International Society of Genetic Geneology). The number of SNPs identified in our samples are shown in Supplementary Table S10. It was not possible to determine the paternal lineage of AH2 due to insufficient number of reads. Derived alleles for F38 and WC1 can be found in Supplementary Table S11 and S12, respectively. WC1 belongs to haplogroup G2b. G2 Y-chromosomal haplogroups are more frequent in Turkey, southern Caucasus and Iran (87, 88). Furthermore, some populations from Pakistan and Southern Asia such as Kalash, Brahui and Pashtun also have the G haplogroup at a higher frequency than other populations (89, 90). The very unique G2b subclade is rarely represented around the world. G2b is mostly found in Pashtuns (91) and among Europeans mainly in Ashkenazi Jews (92). However, it is also found at very low numbers in some other regions in the Near East (87, 90, 93). G2b lineages have not yet been identified in ancient samples. The dominating G-derived lineage in the European Neolithic so far is G2a (3, 4). F38 belongs to sub-haplogroup R1b1a2a2-CTS1078/Z2103. This lineage can be included in the L23(xM412) clade, which is characterized by frequencies higher than 10% in the Caucasus, Turkey, Southeastern Europe, and Circum-Uralic populations, and is mostly found at very low frequencies in Western Europe. This pattern sharply contrasts with the distribution of M412-derived Y-chromosomes, which are very common in Western Europe but rare in the East (94). In ancient DNA studies, derived alleles that define the R1b1a2a2 lineage were found in five Yamnaya individuals and in two samples from the Poltkava culture (3, 17).
Although we observed two reads covering marker J2a1h1-M158, this is a C→T mutation and in addition, there are no upstream markers (ancestral alleles at J-M304, J2-M172 and J2a-M410) that would support including F38 in this lineage. In modern-day Iranian population, the most frequent Y-chromosome haplogroups are J2-M172 (22.5%) and J1-M267 (~10%), and the less frequent ones are R1a-M198 (0-25%) and R1b-M269 (95).