$$\rightleftharpoonup{xx}$$
$$\longleftharp{xx}$$,
$$\longrightharp{xx}$$,
Dans une étude distincte 11, l’ADN a été extrait de la poudre d’os générée à partir de chaque emplacement d’échantillonnage anatomique chez11 individus, en utilisant un protocole d’extraction d’ADN standard optimisé pour de courts fragments de tissu calcifié2. Des bibliothèques simple brin ont ensuite été produites28 et séquencées sur un HiSeq 4000 (extrémité appariée de 75 pb) à une profondeur de ~20 000 000 lectures par échantillon. Les données de séquence résultantes ont ensuite été évaluées pour le contenu endogène de l’ADN humain à l’aide du pipeline EAGER29 (paramètres BWA: longueur de la graine de 32, pénalité de non-concordance de 0,1, filtre de qualité de cartographie de 37). Tous les résultats représentatifs sont présentés en utilisant les mêmes paramètres que Parker et al. 202011 pour la cohérence. Les bibliothèques des parties en poudre de la pars petrosa ont produit, en moyenne, un ADN endogène plus élevé que n’importe lequel des 23 autres sites d’échantillonnage anatomique étudiés (figure 6A-B). Les sept autres sites d’échantillonnage anatomique présentés dans ce protocole (le cément, le premier passage de la chambre pulpaire dentaire et la dentine des molaires permanentes; l’os cortical du corps vertébral et l’arc vertébral supérieur de la vertèbre thoracique; l’os cortical de la touffe apicale de la phalange distale; et l’os cortical du col du talus) ont produit les rendements les plus élevés suivants (sans signification statistique entre ces emplacements d’échantillonnage anatomique; Figure 6A-B; Fichier supplémentaire 1 : EndogenousDNAPreCap). Ces emplacements alternatifs ont tous produit de manière cohérente des rendements d’ADN adéquats pour les analyses de génétique des populations standard telles que les analyses mitochondriales et les analyses de polymorphisme mononucléotidique (SNP). Les taux de duplication dans les bibliothèques provenant de tous les emplacements d’échantillonnage anatomique étaient faibles (facteurs de grappe < 1,2 en moyenne, calculés comme le rapport entre toutes les lectures cartographiques et les lectures cartographiques uniques, tableau 2; Fichier supplémentaire 1 : ClusterFactor), indiquant que toutes les bibliothèques examinées étaient d’une très grande complexité. De même, les estimations moyennes de la contamination exogène de l’ADN humain étaient faibles, atteignant en moyenne < 2 % (contamination par le chromosome X chez les mâles, n = 7, telle que rapportée par le pipeline ANGSD30) dans tous les sites d’échantillonnage anatomique, à l’exception de l’arc vertébral supérieur (contamination moyenne estimée : 2,11 %, avec un échantillon prélevé comme valeur aberrante; KRA005: 19,52%, voir tableau 2; Fichier supplémentaire 1 : Xcontamination). La longueur moyenne des fragments (après filtrage pour éliminer toutes les lectures < 30 pb) était la plus faible dans le matériau prélevé dans la chambre pulpaire dentaire et la dentine, sans variation significative entre les autres emplacements d’échantillonnage anatomique (55,14 pb et 60,22 pb, respectivement, par rapport à une médiane moyenne de 62,87, valeurs p par paires < 0,019, tableau 2; Fichier supplémentaire 1 : AvgFragLength). De plus, les dents et les vertèbres thoraciques contiennent chacune plusieurs emplacements d’échantillonnage anatomique où une récupération élevée de l’ADN endogène a été observée, ce qui les rend particulièrement appropriées comme alternatives à la pars petrosa.

Figure 6 : Contenu de l’ADN humain pour tous les échantillons sélectionnés. Les lignes noires représentent la moyenne globale, tandis que les lignes rouges représentent la médiane (solide: proportion d’ADN humain, en pointillés: lectures humaines cartographiées par million de lectures générées). Les emplacements d’échantillonnage anatomique individuels avec une proportion moyenne d’ADN humain supérieure à la moyenne globale (8,16%) sont colorisés dans toutes les analyses. (A) La proportion de lectures cartographiées avec le génome de référence hg19. La ligne pointillée bleue représente le maximum théorique compte tenu des paramètres de cartographie du pipeline (générés à l’aide de Gargammel31 pour simuler une distribution aléatoire de 5 000 000 de lectures du génome de référence hg19 avec des dommages simulés). Les moyennes individuelles (X noir) et les médianes (cercle rouge) sont rapportées pour les échantillons dont la proportion moyenne d’ADN humain est supérieure à la moyenne globale. Les intervalles de confiance indiquent les limites supérieure et inférieure à l’exclusion des valeurs statistiques aberrantes. (B) Le nombre de lectures uniques correspondant au génome de référence hg19 par million de lectures d’effort de séquençage (extrémité appariée de 75 pb). Les intervalles de confiance indiquent les limites supérieure et inférieure à l’exclusion des valeurs statistiques aberrantes. Cette figure a été adaptée de Parker, C. et al. 202011. Veuillez cliquer ici pour voir une version agrandie de cette figure.
Tableau 2 : Niveaux moyens de duplication (lectures cartographiques/lectures uniques), longueurs moyennes et médianes des fragments et estimations de la contamination du chromosome X pour tous les emplacements d’échantillonnage anatomique. Erreur signalée comme erreur-type de la moyenne. Ce tableau a été adapté de Parker, C. et al. 202011.
| Lieu d’échantillonnage | Facteur de duplication moyen (# lectures mappées /# lectures mappées uniques) | Longueur moyenne des fragments en pb | Proportion moyenne estimée de contamination par le chromosome X |
| Pyramide pétreuse | 1,188 ± 0,006 | 65,40 ± 1,36 | 0,000 ± 0,003 |
| Cément | 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 |
| Pulpe | 1,179 ± 0,024 | 55,14 ± 2,90 | 0,013 ± 0,006 |
| phalange distale | 1,191 ± 0,049 | 65,95 ± 1,08 | 0,013 ± 0,005 |
| Corps vertébral | 1,194 ± 0,037 | 66,14 ± 1,03 | 0,008 ± 0,003 |
| Arc vertébral supérieur | 1,19 ± 0,017 | 63,02 ± 1,23 | 0,021 ± 0,009* |
| Astragale | 1,198 ± 0,010 | 68,20 ± 1,24 | 0,011 ± 0,003 |
| *Échantillon KRA005 retiré comme valeur aberrante à 0,1952 | | | |
Disponibilité du code
Tous les programmes d’analyse et modules R utilisés dans les analyses de ce manuscrit sont disponibles gratuitement auprès de leurs auteurs respectifs. Tous les codes R personnalisés sont disponibles sur demande.
Disponibilité des données
Toutes les données brutes utilisées dans le calcul des résultats représentatifs sont disponibles gratuitement dans le dépôt de données ENA European Nucleotide Archive (numéro d’acquisition PRJ-EB36983) ou dans des documents supplémentaires de Parker, C. et al.11.
Fichier supplémentaire 1. Veuillez cliquer ici pour télécharger ce fichier.