W oddzielnym badaniu11 DNA wyekstrahowano z proszku kostnego uzyskanego z każdego anatomicznego miejsca pobrania u 11 osobników, stosując standardowy protokół ekstrakcji DNA zoptymalizowany dla krótkich fragmentów z tkanki zwapniałej2. Następnie przygotowano biblioteki jednoniciowe28 i zsekwencjonowano je na platformie HiSeq 4000 (75 bp paired-end) do głębokości ok. 20 000 000 odczytów na próbkę. Uzyskane dane sekwencyjne oceniono pod kątem zawartości endogennego ludzkiego DNA przy użyciu potoku EAGER29 (ustawienia BWA: długość ziarna 32, kara za niedopasowanie 0,1, filtr jakości mapowania 37). Dla zachowania spójności wszystkie reprezentatywne wyniki przedstawiono przy użyciu tych samych wskaźników, co w pracy Parker et al. 202011. Biblioteki przygotowane z rozdrobnionych części pars petrosa przyniosły średnio wyższą zawartość endogennego DNA niż każde z pozostałych 23 badanych anatomicznych miejsc pobrania (Rycina 6A-B). Siedem dodatkowych anatomicznych miejsc pobrania przedstawionych w niniejszym protokole (cementum, pierwsze przejście przez jamę miazgi zębowej oraz zębinę zębów trzonowych stałych; kość korowa z trzonu kręgu i górnego łuku kręgu kręgów piersiowych; kość korowa z kępki apikalnej paliczka dalszego oraz kość korowa z szyjki kości skokowej) dało następne najwyższe wydajności (bez istotności statystycznej między tymi anatomicznymi miejscami pobrania; Rycina 6A-B; Plik uzupełniający 1: EndogenousDNAPreCap). Wszystkie te alternatywne miejsca konsekwentnie dostarczały wydajności DNA wystarczającej do standardowych analiz genetyki populacyjnej, takich jak analizy mitochondrialne i analizy polimorfizmu pojedynczego nukleotydu (SNP). Wskaźniki duplikacji w bibliotekach pochodzących ze wszystkich anatomicznych miejsc pobrania były niskie (średni współczynnik klastrowania < 1,2, obliczony jako stosunek wszystkich zmapowanych odczytów do unikalnych odczytów zmapowanych, Tabela 2; Plik uzupełniający 1: ClusterFactor), co wskazuje, że wszystkie przesiewowe biblioteki charakteryzowały się bardzo wysoką złożonością. Podobnie, średnie szacunki kontaminacji egzogennym ludzkim DNA były niskie, wynosząc średnio < 2% (kontaminacja chromosomem X u mężczyzn, n = 7, zgodnie z raportem potoku ANGSD30) we wszystkich anatomicznych miejscach pobrania, z wyjątkiem górnego łuku kręgu (średnia szacowana kontaminacja: 2,11%, z usunięciem jednej próbki jako wartości odstającej; KRA005: 19,52%, patrz Tabela 2; Plik uzupełniający 1: Xcontamination). Średnia długość fragmentów (po odfiltrowaniu wszystkich odczytów < 30 bp) była najniższa w materiale pobranym z jamy miazgi zębowej i zębiny, przy braku istotnych różnic między pozostałymi anatomicznymi miejscami pobrania (odpowiednio 55,14 bp i 60,22 bp w porównaniu do średniej mediany 62,87, wartości p dla porównań parami < 0,019, Tabela 2; Plik uzupełniający 1: AvgFragLength). Dodatkowo, zęby i kręgi piersiowe zawierają po kilka anatomicznych miejsc pobrania, w których zaobserwowano wysoką odzyskiwalność endogennego DNA, co czyni je szczególnie odpowiednimi alternatywami dla pars petrosa.

Rycina 6: Zawartość ludzkiego DNA dla wszystkich badanych próbek. Czarne linie reprezentują średnią ogólną, natomiast czerwone linie reprezentują medianę (linia ciągła: proporcja ludzkiego DNA, linia przerywana: zmapowane ludzkie odczyty na milion wygenerowanych odczytów). Poszczególne anatomiczne miejsca pobrania próbek ze średnią proporcją ludzkiego DNA wyższą niż średnia ogólna (8,16%) są wyróżnione kolorami we wszystkich analizach. (A) Proporcja odczytów mapujących do genomu referencyjnego hg19. Niebieska przerywana linia reprezentuje teoretyczne maksimum przy danych parametrach mapowania potoku przetwarzania (wygenerowane przy użyciu Gargammel31 w celu symulacji losowego rozkładu 5 000 000 odczytów z genomu referencyjnego hg19 z symulowanym uszkodzeniem). Średnie poszczególnych próbek (czarny znak X) i mediany (czerwone koło) są podane dla próbek o średniej proporcji ludzkiego DNA wyższej niż średnia ogólna. Przedziały ufności wskazują górne i dolne granice z wyłączeniem wartości odstających statystycznie. (B) Liczba unikalnych odczytów mapujących do genomu referencyjnego hg19 na milion odczytów z sekwencjonowania (75 bp paired end). Przedziały ufności wskazują górne i dolne granice z wyłączeniem wartości odstających statystycznie. Rycina ta została zaadaptowana z Parker, C. et al. 202011. Kliknij tutaj, aby wyświetlić powiększoną wersję tej ryciny.
Tabela 2: Średnie poziomy duplikacji (odczyty zmapowane/odczyty unikalne), średnie i medianowe długości fragmentów oraz szacunki zanieczyszczenia chromosomem X dla wszystkich anatomicznych miejsc pobrania próbek. Błąd podano jako błąd standardowy średniej. Tabela została opracowana na podstawie pracy Parker, C. et al. 202011.
| Miejsce poboru próbek | Średni współczynnik duplikacji (liczba zmapowanych odczytów / liczba unikalnych zmapowanych odczytów) | Średnia długość fragmentu w bp | Średni szacowany stopień zanieczyszczenia chromosomem X |
| Piramida skroniowa | 1.188 ± 0.006 | 65.40 ± 1.36 | 0.000 ± 0.003 |
| Cement** | 1.197 ± 0.028 | 67.28 ± 1.76 | 0.011 ± 0.003 |
| Dentyna | 1.188 ± 0.061 | 60.22 ± 2.37 | 0.002 ± 0.007 |
| Miazga zęba | 1.179 ± 0.024 | 55.14 ± 2.90 | 0.013 ± 0.006 |
| Paliczek dalszy | 1.191 ± 0.049 | 65.95 ± 1.08 | 0.013 ± 0.005 |
| Trzon kręgu | 1.194 ± 0.037 | 66.14 ± 1.03 | 0.008 ± 0.003 |
| Górny łuk kręgowy | 1.19 ± 0.017 | 63.02 ± 1.23 | 0.021 ± 0.009* |
| Kość skokowa | 1.198 ± 0.010 | 68.20 ± 1.24 | 0.011 ± 0.003 |
| *Próbkę KRA005 usunięto jako wartość odstającą przy 0,1952 | | | |
Dostępność kodu
Wszystkie programy analityczne oraz moduły R wykorzystane w analizach w niniejszym manuskrypcie są bezpłatnie dostępne u ich odpowiednich autorów. Wszystkie niestandardowe kody R są dostępne na żądanie.
Dostępność danych
Wszystkie surowe dane wykorzystane do obliczenia reprezentatywnych wyników są swobodnie dostępne w repozytorium danych European Nucleotide Archive ENA (numer dostępu PRJ-EB36983) lub w materiałach uzupełniających pracy Parker, C. et al.11.
Plik uzupełniający 1. Kliknij tutaj, aby pobrać ten plik.