$$\rightleftharpoonup{xx}$$
$$\longleftharp{xx}$$,
$$\longrightharp{xx}$$,
In een afzonderlijke studie11 werd DNA geëxtraheerd uit botpoeder gegenereerd uit elke anatomische bemonsteringslocatie bij 11 personen, met behulp van een standaard DNA-extractieprotocol dat is geoptimaliseerd voor korte fragmenten uit verkalkt weefsel2. Enkelstrengsbibliotheken werden vervolgensgeproduceerd 28 en gesequenced op een HiSeq 4000 (75 bp paired-end) tot een diepte van ~ 20.000.000 reads per monster. De resulterende sequentiegegevens werden vervolgens geëvalueerd op endogene menselijke DNA-inhoud met behulp van de EAGER-pijplijn29 (BWA-instellingen: zaadlengte van 32, 0,1 mismatch-straf, mappingkwaliteitsfilter van 37). Alle representatieve resultaten worden gerapporteerd met behulp van dezelfde statistieken als Parker et al. 202011 voor consistentie. Bibliotheken van de poedervormige delen van de pars petrosa leverden gemiddeld een hoger endogene DNA op dan alle andere 23 onderzochte anatomische bemonsteringslocaties (figuur 6A-B). De zeven aanvullende anatomische bemonsteringslocaties die in dit protocol worden gepresenteerd (het cementum, de eerste doorgang van de tandpulpkamer en dentine van permanente kiezen; corticale bot van het wervellichaam en de superieure wervelboog van de thoracale wervel; corticale bot uit het apicale plukje van de distale falanx; en corticale bot uit de nek van de talus) produceerden de op een na hoogste opbrengsten (zonder statistische significantie tussen deze anatomische bemonsteringslocaties; Figuur 6A-B; Aanvullend bestand 1: EndogenousDNAPreCap). Deze alternatieve locaties hebben allemaal consistent geproduceerde DNA-opbrengsten opgeleverd die voldoende zijn voor standaard populatiegenetische analyses zoals mitochondriale analyses en single nucleotide polymorphism (SNP) analyses. Duplicatiepercentages in bibliotheken afkomstig van alle anatomische bemonsteringslocaties waren laag (clusterfactoren < gemiddeld 1,2, berekend als de verhouding tussen alle mapping reads en unieke mapping reads, tabel 2; Supplemental File 1: ClusterFactor), wat aangeeft dat alle gescreende bibliotheken van zeer hoge complexiteit waren. Evenzo waren de gemiddelde schattingen van exogene menselijke DNA-besmetting laag, gemiddeld < 2% (X-chromosoomverontreiniging bij mannen, n = 7, zoals gerapporteerd door de ANGSD30-pijplijn) in alle anatomische bemonsteringslocaties behalve de superieure wervelboog (gemiddelde geschatte besmetting: 2,11%, waarbij één monster als uitschieter werd verwijderd; KRA005: 19,52%, zie tabel 2; Aanvullend bestand 1: Xcontaminatie). De gemiddelde fragmentlengte (na filtering om alle afgelezen waarden < 30 bp te verwijderen) was het laagst in het materiaal verzameld uit de tandpulpkamer en dentine, zonder significante variatie tussen andere anatomische bemonsteringslocaties (respectievelijk 55,14 bp en 60,22 bp, in vergelijking met een gemiddelde mediaan van 62,87, paarsgewijze p-waarden < 0,019, tabel 2; Aanvullend bestand 1: AvgFragLength). Bovendien bevatten de tanden en borstwervels elk meerdere anatomische bemonsteringslocaties waar hoog endogene DNA-herstel werd waargenomen, waardoor ze bijzonder geschikt zijn als alternatieven voor de pars petrosa.

Figuur 6: Menselijk DNA-gehalte voor alle gescreende monsters. Zwarte lijnen vertegenwoordigen het totale gemiddelde, terwijl rode lijnen de mediaan vertegenwoordigen (solide: menselijk DNA-aandeel, onderbroken: in kaart gebrachte menselijke metingen per miljoen gegenereerde metingen). Individuele anatomische bemonsteringslocaties met een gemiddeld menselijk DNA-aandeel dat hoger is dan het totale gemiddelde (8,16%) worden in alle analyses ingekleurd. (A) Het aandeel van reads mapping tot het hg19-referentiegenoom. De blauwe stippellijn vertegenwoordigt het theoretische maximum gezien de mappingparameters van de pijplijn (gegenereerd met Gargammel31 om een willekeurige verdeling van 5.000.000 aflezingen van het hg19-referentiegenoom met gesimuleerde schade te simuleren). Individuele gemiddelden (zwarte X) en medianen (rode cirkel) worden gerapporteerd voor die monsters met een hoger gemiddeld menselijk DNA-aandeel dan het totale gemiddelde. Betrouwbaarheidsintervallen geven boven- en ondergrenzen aan, met uitzondering van statistische uitschieters. (B) Het aantal unieke reads mapping naar het hg19-referentiegenoom per miljoen reads van sequencing-inspanning (75 bp gepaard uiteinde). Betrouwbaarheidsintervallen geven boven- en ondergrenzen aan, met uitzondering van statistische uitschieters. Dit cijfer is overgenomen uit Parker, C. et al. 202011. Klik hier om een grotere versie van deze figuur te bekijken.
Tabel 2: Gemiddelde duplicatieniveaus (mapping reads/unique reads), gemiddelde en mediane fragmentlengtes en X-chromosoombesmettingsschattingen voor alle anatomische bemonsteringslocaties. Fout gerapporteerd als de standaardfout van het gemiddelde. Deze tabel is aangepast van Parker, C. et al. 202011.
| Bemonsteringslocatie | Gemiddelde duplicatiefactor (# in kaart gebrachte leeswaarden /# unieke toegewezen leeswaarden) | Gemiddelde fragmentlengte in bp | Gemiddeld geschat aandeel X-chromosoombesmetting |
| Petrous piramide | 1.188 ± 0.006 | 65,40 ± 1,36 | 0,000 ± 0,003 |
| Cementum | 1.197 ± 0.028 | 67,28 ± 1,76 | 0,011 ± 0,003 |
| Dentine | 1.188 ± 0.061 | 60,22 ± 2,37 | 0,002 ± 0,007 |
| Pulp | 1.179 ± 0.024 | 55,14 ± 2,90 | 0,013 ± 0,006 |
| Distale falanx | 1.191 ± 0.049 | 65,95 ± 1,08 | 0,013 ± 0,005 |
| Wervellichaam | 1.194 ± 0.037 | 66,14 ± 1,03 | 0,008 ± 0,003 |
| Superieure wervelboog | 1,19 ± 0,017 | 63,02 ± 1,23 | 0,021 ± 0,009* |
| Talus | 1.198 ± 0.010 | 68,20 ± 1,24 | 0,011 ± 0,003 |
| *Monster KRA005 verwijderd als uitschieter op 0,1952 | | | |
Beschikbaarheid van code
Alle analyseprogramma's en R-modules die in de analyses van dit manuscript worden gebruikt, zijn vrij beschikbaar bij hun respectieve auteurs. Alle aangepaste R-code is op aanvraag beschikbaar.
Beschikbaarheid van gegevens
Alle ruwe gegevens die worden gebruikt bij de berekening van representatieve resultaten zijn vrij beschikbaar in de gegevensopslagplaats van het European Nucleotide Archive ENA (toetredingsnummer PRJ-EB36983) of aanvullende materialen van Parker, C. et al.11.
Aanvullend dossier 1. Klik hier om dit bestand te downloaden.