$$\rightleftharpoonup{xx}$$
$$\longleftharp{xx}$$,
$$\longrightharp{xx}$$,
En un estudio separado 11, se extrajo ADN del polvo óseo generado a partir de cada ubicación de muestreo anatómico en11 individuos, utilizando un protocolo estándar de extracción de ADN optimizado para fragmentos cortos de tejido calcificado2. Luego se produjeron bibliotecas monocatenarias28 y se secuenciaron en un HiSeq 4000 (extremo pareado de 75 pb) a una profundidad de ~ 20,000,000 lecturas por muestra. Los datos de secuencia resultantes se evaluaron para determinar el contenido de ADN humano endógeno utilizando la tubería EAGER29 (configuración de BWA: longitud de semilla de 32, penalización de desajuste de 0.1, filtro de calidad de mapeo de 37). Todos los resultados representativos se informan utilizando las mismas métricas que Parker et al. 202011 para mayor coherencia. Las bibliotecas de las porciones en polvo de la pars petrosa produjeron, en promedio, ADN endógeno más alto que cualquiera de los otros 23 lugares de muestreo anatómico estudiados (Figura 6A-B). Los siete lugares de muestreo anatómico adicionales presentados en este protocolo (el cemento, el primer paso de la cámara de la pulpa dental y la dentina de los molares permanentes; el hueso cortical del cuerpo vertebral y el arco vertebral superior de la vértebra torácica; el hueso cortical del mechón apical de la falange distal; y el hueso cortical del cuello del astrágalo) produjeron los siguientes rendimientos más altos (sin significación estadística entre estos lugares de muestreo anatómico; Figura 6A-B; Archivo complementario 1: DNAPreCap endógeno). Todas estas ubicaciones alternativas produjeron consistentemente rendimientos de ADN adecuados para los análisis genéticos de poblaciones estándar, como los análisis mitocondriales y los análisis de polimorfismo de nucleótido único (SNP). Las tasas de duplicación en las bibliotecas derivadas de todas las ubicaciones de muestreo anatómico fueron bajas (factores de conglomerado < 1,2 en promedio, calculado como la relación entre todas las lecturas de mapeo y lecturas de mapeo únicas, Tabla 2; Archivo complementario 1: ClusterFactor), lo que indica que todas las bibliotecas examinadas eran de muy alta complejidad. Del mismo modo, las estimaciones promedio de contaminación por ADN humano exógeno fueron bajas, con un promedio de < 2% (contaminación del cromosoma X en hombres, n = 7, según lo informado por la tubería ANGSD30) en todos los lugares de muestreo anatómico, excepto en el arco vertebral superior (contaminación promedio estimada: 2.11%, con una muestra eliminada como valor atípico; KRA005: 19,52%, véase el cuadro 2; Archivo complementario 1: Xcontaminación). La longitud promedio del fragmento (después de filtrar para eliminar todas las lecturas < 30 pb) fue más baja en el material recolectado de la cámara de la pulpa dental y la dentina, sin variación significativa entre otras ubicaciones de muestreo anatómico (55.14 pb y 60.22 pb, respectivamente, en comparación con una mediana promedio de 62.87, valores de p por pares < 0.019, Tabla 2; Archivo complementario 1: AvgFragLength). Además, los dientes y las vértebras torácicas contienen múltiples ubicaciones de muestreo anatómico donde se observó una alta recuperación endógena de ADN, lo que las hace particularmente adecuadas como alternativas a la pars petrosa.

Figura 6: Contenido de ADN humano para todas las muestras examinadas. Las líneas negras representan la media general, mientras que las líneas rojas representan la mediana (sólido: proporción de ADN humano, discontinuo: lecturas humanas mapeadas por millón de lecturas generadas). Los lugares de muestreo anatómico individuales con una proporción promedio de ADN humano superior a la media general (8,16%) se colorean en todos los análisis. (A) La proporción de lecturas que se asignan al genoma de referencia hg19. La línea discontinua azul representa el máximo teórico dados los parámetros de mapeo de la tubería (generados usando Gargammel31 para simular una distribución aleatoria de 5,000,000 lecturas del genoma de referencia hg19 con daño simulado). Las medias individuales (X negra) y las medianas (círculo rojo) se informan para aquellas muestras con una proporción promedio de ADN humano más alta que la media general. Los intervalos de confianza indican los límites superior e inferior, excluidos los valores atípicos estadísticos. (B) El número de lecturas únicas que se asignan al genoma de referencia hg19 por millón de lecturas de esfuerzo de secuenciación (75 pb de final pareado). Los intervalos de confianza indican los límites superior e inferior, excluidos los valores atípicos estadísticos. Esta figura ha sido adaptada de Parker, C. et al. 202011. Haga clic aquí para ver una versión más grande de esta figura.
Tabla 2: Niveles promedio de duplicación (lecturas de mapeo/lecturas únicas), longitudes promedio y mediana de fragmentos y estimaciones de contaminación del cromosoma X para todas las ubicaciones de muestreo anatómico. Error notificado como el error estándar de la media. Esta tabla ha sido adaptada de Parker, C. et al. 202011.
| Lugar de muestreo | Factor de duplicación promedio (# lecturas asignadas /# lecturas asignadas únicas) | Longitud media del fragmento en pb | Proporción media estimada de contaminación por el cromosoma X |
| Pirámide de Petrous | 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 |
| Pulpa | 1.179 ± 0.024 | 55,14 ± 2,90 | 0,013 ± 0,006 |
| Falange distal | 1.191 ± 0.049 | 65,95 ± 1,08 | 0,013 ± 0,005 |
| Cuerpo vertebral | 1.194 ± 0.037 | 66,14 ± 1,03 | 0,008 ± 0,003 |
| Arco vertebral superior | 1,19 ± 0,017 | 63,02 ± 1,23 | 0,021 ± 0,009* |
| Astrágalo | 1.198 ± 0.010 | 68,20 ± 1,24 | 0,011 ± 0,003 |
| *Muestra KRA005 eliminada como valor atípico en 0.1952 | | | |
Disponibilidad del código
Todos los programas de análisis y módulos R utilizados en los análisis de este manuscrito están disponibles gratuitamente de sus respectivos autores. Todo el código R personalizado está disponible a pedido.
Disponibilidad de datos
Todos los datos brutos utilizados en el cálculo de resultados representativos están disponibles gratuitamente en el repositorio de datos ENA del Archivo Europeo de Nucleótidos (número de acceso PRJ-EB36983) o en materiales complementarios de Parker, C. et al.11.
Archivo complementario 1. Haga clic aquí para descargar este archivo.