$$\rightleftharpoonup{xx}$$
$$\longleftharp{xx}$$,
$$\longrightharp{xx}$$,
Em um estudo separado11, o DNA foi extraído do pó ósseo gerado de cada local de amostragem anatômica em 11 indivíduos, utilizando um protocolo padrão de extração de DNA otimizado para fragmentos curtos de tecido calcificado2. Bibliotecas de fita simples foram então produzidas28 e sequenciadas em um HiSeq 4000 (75 pb emparelhado-end) a uma profundidade de ~20.000.000 leituras por amostra. Os dados da sequência resultante foram então avaliados quanto ao conteúdo endógeno de DNA humano usando o pipeline EAGER29 (configurações BWA: comprimento da semente de 32, penalidade de incompatibilidade de 0,1, filtro de qualidade de mapeamento de 37). Todos os resultados representativos são relatados usando as mesmas métricas de Parker et al. 202011 para consistência. As bibliotecas das porções em pó da pars petrosa produziram, em média, DNA endógeno mais alto do que qualquer um dos outros 23 locais de amostragem anatômica pesquisados (Figura 6A-B). Os sete locais adicionais de amostragem anatômica apresentados neste protocolo (cemento, primeira passagem da câmara da polpa dentária e dentina dos molares permanentes; osso cortical do corpo vertebral e arco vertebral superior da vértebra torácica; osso cortical do tufo apical da falange distal; e osso cortical do colo do tálus) produziram os próximos maiores rendimentos (sem significância estatística entre esses locais de amostragem anatômica; Figura 6A-B; Arquivo Suplementar 1: EndogenousDNAPreCap). Todos esses locais alternativos produziram consistentemente rendimentos de DNA adequados para análises genéticas populacionais padrão, como análises mitocondriais e análises de polimorfismo de nucleotídeo único (SNP). As taxas de duplicação em bibliotecas decorrentes de todos os locais de amostragem anatômica foram baixas (fatores de conglomerados < 1,2 em média, calculados como a razão entre todas as leituras de mapeamento e leituras únicas de mapeamento, Tabela 2; Arquivo Suplementar 1: ClusterFactor), indicando que todas as bibliotecas examinadas eram de altíssima complexidade. Da mesma forma, as estimativas médias de contaminação exógena do DNA humano foram baixas, com média < 2% (contaminação do cromossomo X em machos, n = 7, conforme relatado pelo pipeline ANGSD30) em todos os locais de amostragem anatômica, exceto no arco vertebral superior (contaminação média estimada: 2,11%, com uma amostra removida como outlier; KRA005: 19,52%, ver Tabela 2; Arquivo Suplementar 1: Xcontaminação). O comprimento médio do fragmento (após filtragem para remoção de todas as leituras < 30 pb) foi menor no material coletado da câmara da polpa dentária e dentina, sem variação significativa entre os demais locais de amostragem anatômica (55,14 pb e 60,22 pb, respectivamente, em comparação com uma mediana média de 62,87, valores de p pareados, < 0,019, Tabela 2; Arquivo suplementar 1: AvgFragLength). Além disso, os dentes e as vértebras torácicas contêm múltiplos locais de amostragem anatômica onde foi observada alta recuperação endógena de DNA, tornando-os particularmente adequados como alternativas à pars petrosa.

Figura 6: Conteúdo de ADN humano para todas as amostras triadas. As linhas pretas representam a média geral, enquanto as linhas vermelhas representam a mediana (sólido: proporção do DNA humano, tracejado: leituras humanas mapeadas por milhão de leituras geradas). Locais de amostragem anatômica individuais com uma proporção média de DNA humano maior do que a média geral (8,16%) são coloridos em todas as análises. (A) A proporção de leituras mapeadas para o genoma de referência hg19. A linha tracejada azul representa o máximo teórico dado os parâmetros de mapeamento da tubulação (gerados usando Gargammel31 para simular uma distribuição aleatória de 5.000.000 de leituras do genoma de referência hg19 com danos simulados). Médias individuais (X preto) e medianas (círculo vermelho) são relatadas para as amostras com uma proporção média de DNA humano maior do que a média geral. Os intervalos de confiança indicam limites superiores e inferiores, excluindo valores atípicos estatísticos. (B) O número de leituras únicas mapeadas para o genoma de referência hg19 por milhão de leituras de esforço de sequenciamento (extremidade emparelhada de 75 pb). Os intervalos de confiança indicam limites superiores e inferiores, excluindo valores atípicos estatísticos. Este número foi adaptado de Parker, C. et al. 202011. Por favor, clique aqui para ver uma versão maior desta figura.
Tabela 2: Níveis médios de duplicação (leituras de mapeamento/leituras únicas), comprimentos médios e medianos de fragmentos e estimativas de contaminação do cromossomo X para todos os locais de amostragem anatômica. Erro relatado como o erro padrão da média. Esta tabela foi adaptada de Parker, C. et al. 202011.
| Local de amostragem | Fator de duplicação médio (# leituras mapeadas /# leituras mapeadas exclusivas) | Comprimento médio do fragmento em pb | Proporção média estimada de contaminação do cromossomo X |
| Pirâmide petrosa | 1,188 ± 0,006 | 65,40 ± 1,36 | 0,000 ± 0,003 |
| Cimento | 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 distal | 1,191 ± 0,049 | 65,95 ± 1,08 | 0,013 ± 0,005 |
| Corpo 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* |
| Tálus | 1,198 ± 0,010 | 68,20 ± 1,24 | 0,011 ± 0,003 |
| *Amostra KRA005 removida como um outlier em 0,1952 | | | |
Disponibilidade de código
Todos os programas de análise e módulos R utilizados nas análises deste manuscrito estão disponíveis gratuitamente junto aos seus respectivos autores. Todo o código R personalizado está disponível mediante solicitação.
Disponibilidade dos dados
Todos os dados brutos utilizados no cálculo dos resultados representativos estão disponíveis gratuitamente no repositório de dados ENA do Arquivo Europeu de Nucleotídeos (número de adesão PRJ-EB36983) ou materiais suplementares de Parker, C. et al.11.
Arquivo Suplementar 1. Clique aqui para baixar este arquivo.