$$\rightleftharpoonup{xx}$$
$$\longleftharp{xx}$$,
$$\longrightharp{xx}$$,
In einer separaten Studie 11 wurde DNA aus Knochenpulver extrahiert, das von jedem anatomischen Probenahmeort bei11 Personen erzeugt wurde, wobei ein Standard-DNA-Extraktionsprotokoll verwendet wurde, das für kurze Fragmente aus verkalktem Gewebe optimiert war2. Einzelsträngige Bibliotheken wurden dann28 produziert und auf einem HiSeq 4000 (75 bp paired-end) bis zu einer Tiefe von ~20.000.000 Lesevorgängen pro Sample sequenziert. Die resultierenden Sequenzdaten wurden dann mit Hilfe der EAGER-Pipeline29 auf endogenen menschlichen DNA-Gehalt ausgewertet (BWA-Einstellungen: Seedlänge von 32, 0,1 Mismatch-Strafe, Mapping-Qualitätsfilter von 37). Alle repräsentativen Ergebnisse werden aus Gründen der Konsistenz unter Verwendung derselben Metriken wie Parker et al. 202011 berichtet. Bibliotheken aus den pulverisierten Teilen der Pars petrosa ergaben im Durchschnitt eine höhere endogene DNA als alle anderen 23 untersuchten anatomischen Probenahmestellen (Abbildung 6A-B). Die sieben zusätzlichen anatomischen Probenahmestellen, die in diesem Protokoll vorgestellt werden (Zement, erster Durchgang der Zahnpulpakammer und Dentin der permanenten Molaren; kortikaler Knochen aus dem Wirbelkörper und oberer Wirbelbogen des Brustwirbels; kortikaler Knochen aus dem apikalen Büschel der distalen Phalanx; und kortikaler Knochen aus dem Hals des Talus) ergaben die nächsthöheren Erträge (ohne statistische Signifikanz zwischen diesen anatomischen Probenahmestellen; Abbildung 6A-B; Ergänzende Datei 1: EndogenousDNAPreCap). Diese alternativen Standorte produzierten alle konsistent DNA-Ausbeuten, die für populationsgenetische Standardanalysen wie mitochondriale Analysen und Einzelnukleotid-Polymorphismus-Analysen (SNP) geeignet sind. Die Duplikationsraten in Bibliotheken, die von allen anatomischen Stichprobenorten stammten, waren niedrig (Clusterfaktoren < 1,2 im Durchschnitt, berechnet als Verhältnis aller Mapping-Lesevorgänge zu eindeutigen Mapping-Lesevorgängen, Tabelle 2; Supplemental File 1: ClusterFactor), was darauf hinweist, dass alle untersuchten Bibliotheken von sehr hoher Komplexität waren. In ähnlicher Weise waren die durchschnittlichen Schätzungen der exogenen DNA-Kontamination beim Menschen niedrig und betrugen durchschnittlich < 2% (X-Chromosomenkontamination bei Männern, n = 7, wie von der ANGSD30-Pipeline berichtet) an allen anatomischen Probenahmestellen mit Ausnahme des oberen Wirbelbogens (durchschnittliche geschätzte Kontamination: 2,11%, wobei eine Probe als Ausreißer entfernt wurde; KRA005: 19,52%, siehe Tabelle 2; Ergänzende Datei 1: Xkontamination). Die durchschnittliche Fragmentlänge (nach Filterung, um alle Messwerte < 30 bp zu entfernen) war am niedrigsten in dem Material, das aus der Zahnpulpakammer und dem Dentin gesammelt wurde, ohne signifikante Unterschiede zwischen anderen anatomischen Probenahmestellen (55,14 bp bzw. 60,22 bp im Vergleich zu einem durchschnittlichen Median von 62,87, paarweise p-Werte < 0,019, Tabelle 2; Ergänzende Datei 1: AvgFragLength). Darüber hinaus enthalten die Zähne und Brustwirbel jeweils mehrere anatomische Probenahmestellen, an denen eine hohe endogene DNA-Erholung beobachtet wurde, was sie besonders als Alternative zur Pars petrosa geeignet macht.

Abbildung 6: Menschlicher DNA-Gehalt für alle untersuchten Proben. Schwarze Linien stellen den Gesamtmittelwert dar, während rote Linien den Median darstellen (durchgehend: menschlicher DNA-Anteil, gestrichelt: abgebildete menschliche Lesevorgänge pro Million generierter Lesevorgänge). Einzelne anatomische Probenahmestellen mit einem durchschnittlichen menschlichen DNA-Anteil über dem Gesamtmittelwert (8,16%) werden in allen Analysen eingefärbt. (A) Der Anteil der Lesevorgänge, die auf das hg19-Referenzgenom abgebildet sind. Die blaue gestrichelte Linie stellt das theoretische Maximum unter Berücksichtigung der Kartierungsparameter der Pipeline dar (generiert mit Gargammel31 , um eine zufällige Verteilung von 5.000.000 Lesevorgängen aus dem hg19-Referenzgenom mit simulierten Schäden zu simulieren). Individuelle Mittelwerte (schwarzes X) und Mediane (roter Kreis) werden für diejenigen Proben mit einem höheren durchschnittlichen menschlichen DNA-Anteil als der Gesamtmittelwert angegeben. Konfidenzintervalle geben obere und untere Grenzen ohne statistische Ausreißer an. (B) Die Anzahl der eindeutigen Lesevorgänge, die dem hg19-Referenzgenom zugeordnet werden, pro Million Lesungen des Sequenzierungsaufwands (75 bp gepaartes Ende). Konfidenzintervalle geben obere und untere Grenzen ohne statistische Ausreißer an. Diese Abbildung wurde von Parker, C. et al. 202011 angepasst. Bitte klicken Sie hier, um eine größere Version dieser Abbildung zu sehen.
Tabelle 2: Durchschnittliche Duplikationsgrade (Mapping-Lesevorgänge/eindeutige Lesevorgänge), durchschnittliche und mediane Fragmentlängen und Schätzungen der X-Chromosomenkontamination für alle anatomischen Probenahmestellen. Fehler, der als Standardfehler des Mittelwerts gemeldet wird. Diese Tabelle wurde angepasst von Parker, C. et al. 202011.
| Probenahmeort | Durchschnittlicher Duplizierungsfaktor (# zugeordnete Lesevorgänge /# eindeutige zugeordnete Lesevorgänge) | Durchschnittliche Fragmentlänge in bp | Durchschnittlicher geschätzter Anteil der X-Chromosom-Kontamination |
| Petrous-Pyramide | 1,188 ± 0,006 | 65,40 ± 1,36 | 0,000 ± 0,003 |
| Zement | 1.197 ± 0.028 | 67,28 ± 1,76 | 0,011 ± 0,003 |
| Dentin | 1.188 ± 0.061 | 60,22 ± 2,37 | 0,002 ± 0,007 |
| Fruchtfleisch | 1,179 ± 0,024 | 55,14 ± 2,90 | 0,013 ± 0,006 |
| Endglied | 1,191 ± 0,049 | 65,95 ± 1,08 | 0,013 ± 0,005 |
| Wirbelkörper | 1,194 ± 0,037 | 66,14 ± 1,03 | 0,008 ± 0,003 |
| Oberer Wirbelbogen | 1,19 ± 0,017 | 63,02 ± 1,23 | 0,021 ± 0,009* |
| Sprungbein | 1.198 ± 0.010 | 68,20 ± 1,24 | 0,011 ± 0,003 |
| *Probe KRA005 als Ausreißer bei 0,1952 entfernt | | | |
Code-Verfügbarkeit
Alle Analyseprogramme und R-Module, die in den Analysen dieses Manuskripts verwendet werden, sind bei ihren jeweiligen Autoren frei erhältlich. Alle benutzerdefinierten R-Codes sind auf Anfrage erhältlich.
Datenverfügbarkeit
Alle Rohdaten, die für die Berechnung repräsentativer Ergebnisse verwendet werden, sind im European Nucleotide Archive ENA Data Repository (Zugangsnummer PRJ-EB36983) oder ergänzenden Materialien von Parker, C. et al.11 frei verfügbar.
Ergänzende Datei 1. Bitte klicken Sie hier, um diese Datei herunterzuladen.