$$\rightleftharpoonup{xx}$$
$$\longleftharp{xx}$$,
$$\longrightharp{xx}$$,
In uno studio separato 11, il DNA è stato estratto dalla polvere ossea generata da ciascuna posizione di campionamento anatomico in11 individui, utilizzando un protocollo di estrazione del DNA standard ottimizzato per brevi frammenti di tessuto calcificato2. Le librerie a singolo filamento sono state quindi prodotte28 e sequenziate su un HiSeq 4000 (75 bp paired-end) ad una profondità di ~ 20.000.000 di letture per campione. I dati di sequenza risultanti sono stati quindi valutati per il contenuto di DNA umano endogeno utilizzando la pipeline EAGER29 (impostazioni BWA: lunghezza del seme di 32, penalità di mancata corrispondenza di 0,1, filtro di qualità della mappatura di 37). Tutti i risultati rappresentativi sono riportati utilizzando le stesse metriche di Parker et al. 202011 per coerenza. Le librerie delle porzioni in polvere della pars petrosa hanno prodotto, in media, DNA endogeno superiore rispetto a qualsiasi altra delle altre 23 posizioni di campionamento anatomico esaminate (Figura 6A-B). Le sette ulteriori sedi di campionamento anatomico presentate in questo protocollo (il cemento, il primo passaggio della camera pulpare dentale e la dentina dei molari permanenti; l'osso corticale dal corpo vertebrale e l'arco vertebrale superiore della vertebra toracica; l'osso corticale dal ciuffo apicale della falange distale; e l'osso corticale dal collo dell'astragalo) hanno prodotto le rese successive più elevate (senza alcuna significatività statistica tra queste posizioni di campionamento anatomico; Figura 6A-B; File supplementare 1: DNAPreCap endogeno). Queste posizioni alternative producono DNA in modo coerente adeguate per le analisi standard di genetica delle popolazioni, come le analisi mitocondriali e le analisi del polimorfismo a singolo nucleotide (SNP). I tassi di duplicazione nelle librerie derivanti da tutte le posizioni di campionamento anatomico erano bassi (fattori di cluster < 1,2 in media, calcolati come rapporto tra tutte le letture di mappatura e letture di mappatura univoche, Tabella 2; Supplemental File 1: ClusterFactor), indicando che tutte le librerie sottoposte a screening erano di complessità molto elevata. Allo stesso modo, le stime medie di contaminazione del DNA umano esogeno erano basse, con una media < 2% (contaminazione del cromosoma X nei maschi, n = 7, come riportato dalla pipeline ANGSD30) in tutte le sedi di campionamento anatomico ad eccezione dell'arco vertebrale superiore (contaminazione media stimata: 2,11%, con un campione rimosso come anomalia; KRA005: 19,52%, cfr. tabella 2; File supplementare 1: Xcontamination). La lunghezza media del frammento (dopo il filtraggio per rimuovere tutte le letture < 30 bp) era più bassa nel materiale raccolto dalla camera pulpare dentale e dalla dentina, senza variazioni significative tra le altre posizioni di campionamento anatomico (55,14 bp e 60,22 bp, rispettivamente rispetto a una mediana media di 62,87, valori p a coppie < 0,019, Tabella 2; File supplementare 1: AvgFragLength). Inoltre, i denti e le vertebre toraciche contengono ciascuno più punti di campionamento anatomico in cui è stato osservato un elevato recupero endogeno del DNA, rendendoli particolarmente adatti come alternative alla pars petrosa.

Figura 6: Contenuto di DNA umano per tutti i campioni esaminati. Le linee nere rappresentano la media complessiva, mentre le linee rosse rappresentano la mediana (solido: proporzione del DNA umano, tratteggiato: letture umane mappate per milione di letture generate). Le singole posizioni di campionamento anatomico con una proporzione media di DNA umano superiore alla media complessiva (8,16%) sono colorate in tutte le analisi. (A) La proporzione di letture mappate al genoma di riferimento hg19. La linea tratteggiata blu rappresenta il massimo teorico dati i parametri di mappatura della pipeline (generati utilizzando Gargammel31 per simulare una distribuzione casuale di 5.000.000 di letture dal genoma di riferimento hg19 con danno simulato). Le medie individuali (X nera) e le mediane (cerchio rosso) sono riportate per quei campioni con una proporzione media di DNA umano più elevata rispetto alla media complessiva. Gli intervalli di confidenza indicano i limiti superiore e inferiore escludendo i valori anomali statistici. (B) Il numero di letture uniche mappate al genoma di riferimento hg19 per milione di letture dello sforzo di sequenziamento (75 bp accoppiati end). Gli intervalli di confidenza indicano i limiti superiore e inferiore escludendo i valori anomali statistici. Questa figura è stata adattata da Parker, C. et al. 202011. Fare clic qui per visualizzare una versione ingrandita di questa figura.
Tabella 2: Livelli medi di duplicazione (letture di mappatura/letture univoche), lunghezze medie e mediane dei frammenti e stime della contaminazione del cromosoma X per tutte le posizioni di campionamento anatomico. Errore segnalato come errore standard della media. Questa tabella è stata adattata da Parker, C. et al. 202011.
| Luogo di campionamento | Fattore di duplicazione medio (# letture mappate /# letture mappate univoche) | Lunghezza media del frammento in bp | Percentuale media stimata di contaminazione del cromosoma X |
| Piramide petrosa | 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 |
| Dentina | 1,188 ± 0,061 | 60,22 ± 2,37 | 0,002 ± 0,007 |
| Polpa | 1,179 ± 0,024 | 55,14 ± 2,90 | 0,013 ± 0,006 |
| Falange distale | 1,191 ± 0,049 | 65,95 ± 1,08 | 0,013 ± 0,005 |
| Corpo vertebrale | 1,194 ± 0,037 | 66,14 ± 1,03 | 0,008 ± 0,003 |
| Arco vertebrale superiore | 1,19 ± 0,017 | 63,02 ± 1,23 | 0,021 ± 0,009* |
| Astragalo | 1,198 ± 0,010 | 68,20 ± 1,24 | 0,011 ± 0,003 |
| *Campione KRA005 rimosso come valore anomalo a 0,1952 | | | |
Disponibilità del codice
Tutti i programmi di analisi e i moduli R utilizzati nelle analisi di questo manoscritto sono liberamente disponibili presso i rispettivi autori. Tutti i codici R personalizzati sono disponibili su richiesta.
Disponibilità dei dati
Tutti i dati grezzi utilizzati nel calcolo dei risultati rappresentativi sono liberamente disponibili nell'archivio europeo dei nucleotidi ENA (numero di adesione PRJ-EB36983) o nei materiali supplementari di Parker, C. et al.11.
File supplementare 1. Clicca qui per scaricare questo file.