Human genomes from the Green Sahara reveal ancestral North African lineage
#1
https://www.ebi.ac.uk/ena/browser/view/PRJEB84057

Screening data for libraries created in this study (TKH001.A0101 from TKH001, and XXX001.A18276 and XXX001.A18277 from TKH009) were obtained through shallow shotgun sequencing. SNP enrichment was performed using three capture sets: the 1240k and Twist panels for population genetic analyses, and the Archaic Admixture panel for investigating archaic ancestry. All capture data presented here were aligned to the human reference genome (hg19, GRCh37), filtered for PCR duplicates, filtered for mapping quality of 25 and a minimum length of 30 base pairs. The shallow shotgun data was processed similarly, except with a mapping quality threshold of 30. Additionally, Archaic Admixture capture data were generated from 22 previously existing libraries (van de Loosdrecht et al., 2018) derived from seven Taforalt individuals (TAF09, TAF010, TAF011, TAF012, TAF13, TAF014, TAF015). We provide BAM files from these libraries, aligned to the human reference genome (hg19, GRCh37) and filtered for PCR duplicates. Merged data for each individual, processed analogously to the Takarkori data, are also included.

Additional abstract from last year:

Known as one of the most arid areas on the planet today, the Sahara Desert was, in fact, a green savannah in the Holocene, dotted by forests and water bodies that promoted human occupation and fostered pastoralism. Due to the present-day climatic conditions, ancient DNA does not preserve well in the region, resulting in limited knowledge of the Sahara’s demographic past. Here, we report the first ancient human genome-wide data from the Saharan Pastoral Neolithic. We obtained genomic data from two ca. 7000-yearold female pastoralists buried in the Takarkori rock shelter at the heart of the Tadrart Acacus massif in southwestern Libya, which was used as a burial ground by pastoral communities. We find that the majority of the Takarkori individuals’ ancestry stems from a previously unknown lineage that appears to have remained isolated for most of its existence. Both individuals are most closely related to the preceding 15,000-year-old foragers from Morocco associated with the Iberomaurusian techno-complex, whereas both Takarkori and Iberomaurusian individuals are distantly related to sub-Saharan African lineages. The quality of one of the genomes from Takarkori is sufficient to detect prospective Neandertal ancestry and we find evidence for few segments of ancestry that sum to a total comparable to that detected in the genomes of sub-Saharan Africans. Our results therefore support a model of cultural diffusion, rather than human migration, for the emergence of pastoralist subsistence in the Sahara region.

Related Articles while we wait for the actual article:

Ancestral mitochondrial N lineage from the Neolithic ‘green’ Sahara
https://www.nature.com/articles/s41598-019-39802-1

Aquatic fauna from the Takarkori rock shelter reveals the Holocene central Saharan climate and palaeohydrography
https://journals.plos.org/plosone/articl...0228588%20
AimSmall, alchemist223, Inquirer And 13 others like this post
Reply
#2
To start the discussion, here are some f3 stats, an MDS plot based on them and a qpgraph with the European side based on a simplified version of one of Kale's plus Neandertal and Africans added.

Code:
Source1    Source2    Target    f_3    std.err    Z    SNPs
TKH001    Dinka    Mbuti    0.107739    0.001468    73.384    744819
TKH001    TKH009    Mbuti    0.267794    0.009219    29.047    10881
TKH001    Mende    Mbuti    0.098476    0.001543    63.823    709880
TKH001    Taforalt    Mbuti    0.181195    0.002569    70.543    738796
TKH001    KotiasKide25K    Mbuti    0.140973    0.002207    63.882    730575
TKH001    Natufian    Mbuti    0.147734    0.003824    38.633    73804
TKH001    ZlatyKun    Mbuti    0.130789    0.002199    59.485    717060
TKH001    Iran_N    Mbuti    0.136850    0.001882    72.721    765641
TKH001    ShumLaka    Mbuti    0.066634    0.001315    50.676    720247
TKH001    Yoruba    Mbuti    0.100861    0.001306    77.228    805308
TKH001    Turkey_Boncuklu    Mbuti    0.143429    0.002039    70.331    692302
TKH001    S_Afr_2000BP    Mbuti    0.040332    0.001153    34.993    718738
TKH001    Mota    Mbuti    0.093934    0.001690    55.573    708701
TKH009    Dinka    Mbuti    0.109282    0.003860    28.312    14301
TKH009    TKH001    Mbuti    0.267794    0.009219    29.047    10881
TKH009    Mende    Mbuti    0.097108    0.004302    22.572    13667
TKH009    Taforalt    Mbuti    0.185472    0.006455    28.734    14197
TKH009    KotiasKide25K    Mbuti    0.140198    0.006506    21.550    13988
TKH009    Natufian    Mbuti    0.159012    0.017967    8.850    2277
TKH009    ZlatyKun    Mbuti    0.139927    0.006734    20.779    13788
TKH009    Iran_N    Mbuti    0.147472    0.005543    26.607    14750
TKH009    ShumLaka    Mbuti    0.068599    0.004335    15.826    13814
TKH009    Yoruba    Mbuti    0.102412    0.003220    31.805    15454
TKH009    Turkey_Boncuklu    Mbuti    0.142637    0.006434    22.169    14019
TKH009    S_Afr_2000BP    Mbuti    0.046160    0.004499    10.261    13827
TKH009    Mota    Mbuti    0.096996    0.005752    16.864    13612


.png   Takarkori MDS.png (Size: 20.99 KB / Downloads: 636)

   
JMcB, Megalophias, Piquerobi And 12 others like this post
Reply
#3
For those interested while we wait for more official sources, here are bed/bin/fam files for TKH001 and TKH009. The link will work for a month.

In the zip file I provide 2 versions: the usual without filtering, except for quality of course and those I prefer to use with trimming to remove most of ancient dna damage.  I also provide a version of TKH001 with only Shotgun and Twist Capture merged so as to have as little capture bias as possible.

Here is the link: https://ufile.io/j2wh0t12

Here are the stats:

Code:
Full Versions:
SampleName    TotalSites    NonMissingCalls    avgRawReads    avgDamageCleanedReads    avgSampledFrom
TKH001    1233013    1082262    3.590629620287864    3.590629620287864    3.5755770620423304
TKH009    1233013    23950    2.0442606850049432e-2    2.0442606850049432e-2    2.0382591262217024e-2

Trimmed Versions:

SampleName    TotalSites    NonMissingCalls    avgRawReads    avgDamageCleanedReads    avgSampledFrom
TKH001_tr.SG.TW    1233013    967964    2.3013447546781745    2.3013447546781745    2.294842795655845
TKH001_tr_all    1233013    1019665    2.752770652053141    2.752770652053141    2.7453781914708117
TKH009_tr    1233013    18506    1.5771123256608e-2    1.5771123256608e-2    1.573462729103424e-2
gamalnasser, AimSmall, JMcB And 4 others like this post
Reply
#4
Super interesting stuff! The Takarkori samples seem to really help resolve the phylogeny. Looks like a big impact of ANA on West and East Africa in the qpgraph. Connected to the spread of Y-DNA E and mtDNA L3?
Tatsuya, billh, Psynome And 2 others like this post
Reply
#5
Dang this TKH001 chap looks quite interesting, thanks for the data and that workup Crashdoc.
JMcB, Tatsuya, Necrontyr like this post
Reply
#6
(02-05-2025, 04:55 AM)Lank Wrote: Super interesting stuff! The Takarkori samples seem to really help resolve the phylogeny. Looks like a big impact of ANA on West and East Africa in the qpgraph. Connected to the spread of Y-DNA E and mtDNA L3?

This is in agreement with Lipson et al (2020), as can be seen in Figure 5. It models the bulk of the ancestry in Mota (71%) and Niger-Congo peoples (69%) as having derived from the same source that produced the Aboriginal African ancestry in the Taforalt hunter-gatherers. 

Figure 5:

Show Content

Figure 4(a) models 88% of the ancestry in Niger-Congo peoples as having derived from this source (though it excludes the Taforalt hunter-gatherers).

Show Content
Ronalawe, Psynome, Piquerobi And 3 others like this post
Reply
#7
qpadm(data, left = c("Dinka.DG", "Morocco_Iberomaurusian.AG", "Yoruba.DG"), right = c("Chimp.REF", "Mbuti.DG", "South_Africa_1900BP.SG", "Ethiopia_4500BP.SG", "Russia_Kostenki14_UP.SG", "Russia_UstIshim_IUP.DG", "Czechia_ZlatyKun_IUP.AG.SG", "Russia_MA1_UP.SG", "Cameroon_ShumLaka_SMA.AG", "Belgium_GoyetQ116_1_UP.AG", "Georgia_Kotias_UP.SG", "Italy_Epigravettian.SG", "Turkey_Marmara_Barcin_N.AG"), target = c("Libya_Takarkori_N"))
ℹ Computing f4 stats...
ℹ Computing admixture weights...
ℹ Computing standard errors...
ℹ Computing number of admixture waves...

warning: solve(): system is singular (rcond: 3.80927e-17); attempting approx solution
$weights
# A tibble: 3 × 5
target left weight se z
<chr> <chr> <dbl> <dbl> <dbl>
1 Libya_Takarkori_N Dinka.DG 0.389 0.0878 4.43
2 Libya_Takarkori_N Morocco_Iberomaurusian.AG 0.375 0.0142 26.3
3 Libya_Takarkori_N Yoruba.DG 0.237 0.0804 2.94

$f4
# A tibble: 624 × 9
pop1 pop2 pop3 pop4 est se z p weight
<chr> <chr> <chr> <chr> <dbl> <dbl> <dbl> <dbl> <dbl>
1 Libya_Takarkori_N Dinka.DG Belgium_GoyetQ116_1_UP.AG Cameroon_ShumLaka_SMA.AG 0.00749 0.000506 14.8 1.67e- 49 0.389
2 Libya_Takarkori_N Morocco_Iberomaurusian.AG Belgium_GoyetQ116_1_UP.AG Cameroon_ShumLaka_SMA.AG -0.0158 0.000566 -28.0 3.79e-172 0.375
3 Libya_Takarkori_N Yoruba.DG Belgium_GoyetQ116_1_UP.AG Cameroon_ShumLaka_SMA.AG 0.00942 0.000490 19.2 2.35e- 82 0.237
4 Libya_Takarkori_N fit Belgium_GoyetQ116_1_UP.AG Cameroon_ShumLaka_SMA.AG -0.000796 0.000466 -1.71 8.75e- 2 NA
5 Libya_Takarkori_N Dinka.DG Belgium_GoyetQ116_1_UP.AG Chimp.REF 0.00780 0.000542 14.4 5.98e- 47 0.389
6 Libya_Takarkori_N Morocco_Iberomaurusian.AG Belgium_GoyetQ116_1_UP.AG Chimp.REF -0.0173 0.000604 -28.6 3.41e-180 0.375
7 Libya_Takarkori_N Yoruba.DG Belgium_GoyetQ116_1_UP.AG Chimp.REF 0.0116 0.000513 22.7 9.22e-114 0.237
8 Libya_Takarkori_N fit Belgium_GoyetQ116_1_UP.AG Chimp.REF -0.000694 0.000500 -1.39 1.65e- 1 NA
9 Libya_Takarkori_N Dinka.DG Belgium_GoyetQ116_1_UP.AG Czechia_ZlatyKun_IUP.AG.SG 0.00150 0.000472 3.18 1.47e- 3 0.389
10 Libya_Takarkori_N Morocco_Iberomaurusian.AG Belgium_GoyetQ116_1_UP.AG Czechia_ZlatyKun_IUP.AG.SG -0.00367 0.000584 -6.28 3.39e- 10 0.375
# ℹ 614 more rows
# ℹ Use `print(n = ...)` to see more rows

$rankdrop
# A tibble: 3 × 7
f4rank dof chisq p dofdiff chisqdiff p_nested
<int> <int> <dbl> <dbl> <int> <dbl> <dbl>
1 2 10 14.1 1.69e- 1 12 492. 8.91e-98
2 1 22 507. 3.14e-93 14 8016. 0
3 0 36 8523. 0 NA NA NA

$popdrop
# A tibble: 7 × 14
pat wt dof chisq p f4rank Dinka.DG Morocco_Iberomaurusian.AG Yoruba.DG feasible best dofdiff chisqdiff p_nested
<chr> <dbl> <dbl> <dbl> <dbl> <dbl> <dbl> <dbl> <dbl> <lgl> <lgl> <dbl> <dbl> <dbl>
1 000 0 10 14.1 1.69e- 1 2 0.389 0.375 0.237 TRUE NA NA NA NA
2 001 1 11 73.2 3.00e-11 1 0.649 0.351 NA TRUE TRUE 0 -333. 1
3 010 1 11 406. 3.45e-80 1 3.41 NA -2.41 FALSE TRUE 0 314. 0
4 100 1 11 92.0 6.75e-15 1 NA 0.416 0.584 TRUE TRUE NA NA NA
5 011 2 12 4445. 0 0 1 NA NA TRUE NA NA NA NA
6 101 2 12 3013. 0 0 NA 1 NA TRUE NA NA NA NA
7 110 2 12 8081. 0 0 NA NA 1 TRUE NA NA NA NA
kolompar, Kale, Necrontyr And 1 others like this post
Reply
#8
They apparently have a unique mtDNA N subclade derived from the root. With mostly ANA ancestry, about 23% East African, and 15% Levantine ancestry (according to the qpgraph by crashdoc). I wonder if this indicates an origin of mtDNA N in North Africa?
Ronalawe and alchemist223 like this post
Reply
#9
(02-05-2025, 04:22 AM)crashdoc Wrote: To start the discussion, here are some f3 stats, an MDS plot based on them and a qpgraph with the European side based on a simplified version of one of Kale's plus Neandertal and Africans added.

Code:
Source1    Source2    Target    f_3    std.err    Z    SNPs
TKH001    Dinka    Mbuti    0.107739    0.001468    73.384    744819
TKH001    TKH009    Mbuti    0.267794    0.009219    29.047    10881
TKH001    Mende    Mbuti    0.098476    0.001543    63.823    709880
TKH001    Taforalt    Mbuti    0.181195    0.002569    70.543    738796
TKH001    KotiasKide25K    Mbuti    0.140973    0.002207    63.882    730575
TKH001    Natufian    Mbuti    0.147734    0.003824    38.633    73804
TKH001    ZlatyKun    Mbuti    0.130789    0.002199    59.485    717060
TKH001    Iran_N    Mbuti    0.136850    0.001882    72.721    765641
TKH001    ShumLaka    Mbuti    0.066634    0.001315    50.676    720247
TKH001    Yoruba    Mbuti    0.100861    0.001306    77.228    805308
TKH001    Turkey_Boncuklu    Mbuti    0.143429    0.002039    70.331    692302
TKH001    S_Afr_2000BP    Mbuti    0.040332    0.001153    34.993    718738
TKH001    Mota    Mbuti    0.093934    0.001690    55.573    708701
TKH009    Dinka    Mbuti    0.109282    0.003860    28.312    14301
TKH009    TKH001    Mbuti    0.267794    0.009219    29.047    10881
TKH009    Mende    Mbuti    0.097108    0.004302    22.572    13667
TKH009    Taforalt    Mbuti    0.185472    0.006455    28.734    14197
TKH009    KotiasKide25K    Mbuti    0.140198    0.006506    21.550    13988
TKH009    Natufian    Mbuti    0.159012    0.017967    8.850    2277
TKH009    ZlatyKun    Mbuti    0.139927    0.006734    20.779    13788
TKH009    Iran_N    Mbuti    0.147472    0.005543    26.607    14750
TKH009    ShumLaka    Mbuti    0.068599    0.004335    15.826    13814
TKH009    Yoruba    Mbuti    0.102412    0.003220    31.805    15454
TKH009    Turkey_Boncuklu    Mbuti    0.142637    0.006434    22.169    14019
TKH009    S_Afr_2000BP    Mbuti    0.046160    0.004499    10.261    13827
TKH009    Mota    Mbuti    0.096996    0.005752    16.864    13612
is this qpgraph tree automated ?
Reply
#10
Natufians should be included in the graph to test whether Takarkori has affinity with Natufians independent of Taforalt-like ancestry
Necrontyr, billh, Tatsuya like this post
Reply
#11
Takarkori does not seem to need any Natufian
Tatsuya and Necrontyr like this post


Attached Files Thumbnail(s)
   
Reply
#12
(02-05-2025, 10:23 AM)Norfern-Ostrobothnian Wrote: Takarkori does not seem to need any Natufian

I looked at it looks like Saharan HG lineage , i don't think  they have much Neolithic ancestry, 
they show affinity to Dinka and East african pastoral neolithic ..probably can fit as dinka + some iberomaurusian admixture..not sure tho.
Tatsuya and Norfern-Ostrobothnian like this post
Reply
#13
(02-05-2025, 10:46 AM)masta Wrote:
(02-05-2025, 10:23 AM)Norfern-Ostrobothnian Wrote: Takarkori does not seem to need any Natufian

I looked at it looks like Saharan HG lineage , i don't think  they have much Neolithic ancestry, 
they show affinity to Dinka and East african pastoral neolithic ..probably can fit as dinka + some iberomaurusian admixture..not sure tho.

Both my qpAdm and qpGraph also has some West African ancestry.
I think it can be modelled in G25 once we get the coordinates for it.
masta and Tatsuya like this post
Reply
#14
(02-05-2025, 09:23 AM)Lank Wrote: They apparently have a unique mtDNA N subclade derived from the root. With mostly ANA ancestry, about 23% East African, and 15% Levantine ancestry (according to the qpgraph by crashdoc). I wonder if this indicates an origin of mtDNA N in North Africa?

Not at all it is a very derived variety and the main special thing about it is it went extinct. With the low resolution of mtdna all branches of N are unique.
I looked at the mutations of this N branch it literally has more differences from Basal N either in theory or compared to ancient N samples like Oase than most modern people do. It seems for some weird reason this type of N had a higher rate of mutation. Iirc and afaik
TanTin and Tatsuya like this post
Reply
#15
(02-05-2025, 01:49 PM)venustas Wrote:
(02-05-2025, 09:23 AM)Lank Wrote: They apparently have a unique mtDNA N subclade derived from the root. With mostly ANA ancestry, about 23% East African, and 15% Levantine ancestry (according to the qpgraph by crashdoc). I wonder if this indicates an origin of mtDNA N in North Africa?

Not at all it is a very derived variety and the main special thing about it is it went extinct. With the low resolution of mtdna all branches of N are unique.
I looked at the mutations of this N branch it literally has more differences from Basal N either in theory or compared to ancient N samples like Oase than most modern people do. It seems for some weird reason this type of N had a higher rate of mutation. Iirc and afaik

Obviously, I didn't mean these samples literally carry the root mtDNA N lineage (they lived 7 kya, after all, not 50-60 kya). But their mtDNA N diverged from the root, just like many other N lineage, but in this case they are apparently of mostly non-Eurasian ancestry.

This is obviously not enough to say that mtDNA N originated in North Africa, but it is certainly a possibility. And not such a crazy one, when the Neanderthal admixtures in proto-Eurasians apparently dates to just around 45-50 kya.

I wonder if mtDNA M1 descends from a similar group in the Paleolithic, whether that was in North Africa or just outside in West Asia.
alchemist223 likes this post
Reply


Forum Jump:


Users browsing this thread: 1 Guest(s)