$$\rightleftharpoonup{xx}$$
$$\longleftharp{xx}$$,
$$\longrightharp{xx}$$,
В отдельном исследовании11 ДНК была извлечена из костного порошка, полученного из каждого анатомического места отбора проб у 11 человек, с использованием стандартного протокола экстракции ДНК, оптимизированного для коротких фрагментов из кальцинированной ткани2. Затем было создано28 одноцепочечных библиотек и секвенировано на HiSeq 4000 (75 bp paired-end) на глубину ~ 20 000 000 считываний на образец. Полученные данные последовательности затем оценивали на предмет эндогенного содержания ДНК человека с использованием конвейера EAGER29 (настройки BWA: длина семени 32, штраф за несоответствие 0,1, фильтр качества отображения 37). Все репрезентативные результаты сообщаются с использованием тех же показателей, что и Parker et al. 202011 для согласованности. Библиотеки из порошкообразных частей pars petrosa дали, в среднем, более высокую эндогенную ДНК, чем любое из других 23 исследованных анатомических мест отбора проб (рисунок 6A-B). Семь дополнительных анатомических мест отбора проб, представленных в настоящем протоколе (цемент, первый проход камеры пульпы зуба и дентин постоянных моляров; кортикальная кость из тела позвонка и верхняя позвоночная дуга грудного позвонка; кортикальная кость из апикального пучка дистальной фаланги; и кортикальная кость из шейки таранной кости) дали следующий по величине урожай (без статистической значимости между этими анатомическими местами отбора проб; Рисунок 6А-В; Дополнительный файл 1: EndogenousDNAPreCap). Все эти альтернативные места последовательно продуцируют ДНК, достаточные для стандартных популяционных генетических анализов, таких как митохондриальный анализ и анализ однонуклеотидного полиморфизма (SNP). Частота дублирования в библиотеках, обусловленных всеми анатомическими местами выборки, была низкой (кластерные коэффициенты < в среднем 1,2, рассчитанные как отношение всех сопоставлений считываний к уникальным показаниям сопоставления, таблица 2; Дополнительный файл 1: ClusterFactor), указывающий на то, что все экранируемые библиотеки имеют очень высокую сложность. Аналогичным образом, средние оценки экзогенного загрязнения ДНК человека были низкими, составляя в среднем < 2% (загрязнение Х-хромосом у мужчин, n = 7, как сообщается в конвейере ANGSD30) во всех анатомических местах отбора проб, за исключением верхней позвоночной дуги (среднее оценочное загрязнение: 2,11%, при этом один образец был удален в качестве выброса; KRA005: 19,52%, см. таблицу 2; Дополнительный файл 1: Загрязнение). Средняя длина фрагмента (после фильтрации для удаления всех показаний < 30 bp) была самой низкой в материале, собранном из камеры пульпы зуба и дентина, без существенных различий между другими анатомическими местами отбора проб (55,14 bp и 60,22 bp, соответственно, по сравнению со средней медианой 62,87, попарные p-значения < 0,019, таблица 2; Дополнительный файл 1: AvgFragLength). Кроме того, зубы и грудные позвонки содержат несколько анатомических мест отбора проб, где наблюдалось высокое эндогенное восстановление ДНК, что делает их особенно подходящими в качестве альтернативы pars petrosa.

Рисунок 6: Содержание ДНК человека для всех проверенных образцов. Черные линии представляют общее среднее значение, в то время как красные линии представляют медиану (сплошная: пропорция ДНК человека, пунктирная: нанесенное на карту человеческое чтение на миллион сгенерированных чтений). Отдельные анатомические места отбора проб со средней долей ДНК человека выше общей средней (8,16%) окрашиваются во всех анализах. (A) Доля считываемых карт с эталонным геномом hg19. Синяя пунктирная линия представляет собой теоретический максимум с учетом параметров отображения конвейера (сгенерированный с использованием Gargammel31 для моделирования случайного распределения 5 000 000 считываний из эталонного генома hg19 с смоделированным повреждением). Индивидуальные средства (черный X) и медианы (красный круг) сообщаются для образцов с более высокой средней долей ДНК человека, чем общее среднее значение. Доверительные интервалы обозначают верхнюю и нижнюю границы, исключая статистические выбросы. (B) Количество уникальных считываний, сопоставленных с эталонным геномом hg19 на миллион считываний усилий секвенирования (75 bp парный конец). Доверительные интервалы обозначают верхнюю и нижнюю границы, исключая статистические выбросы. Эта цифра была адаптирована из Parker, C. et al. 202011. Пожалуйста, нажмите здесь, чтобы просмотреть увеличенную версию этого рисунка.
Таблица 2: Средние уровни дупликации (картирование считываний/уникальных считываний), средняя и медианная длина фрагментов и оценки загрязнения Х-хромосом для всех анатомических мест выборки. Ошибка, сообщаемая как стандартная погрешность среднего значения. Эта таблица была адаптирована из Parker, C. et al. 202011.
| Место отбора проб | Средний коэффициент дублирования (# сопоставленных чтений /# уникальных сопоставленных чтений) | Средняя длина фрагмента в bp | Средняя оценочная доля загрязнения Х-хромосом |
| Пирамида Петруса | 1.188 ± 0.006 | 65.40 ± 1.36 | 0.000 ± 0.003 |
| Цемент | 1.197 ± 0.028 | 67.28 ± 1.76 | 0.011 ± 0.003 |
| Дентин | 1.188 ± 0.061 | 60.22 ± 2.37 | 0.002 ± 0.007 |
| Пульпа | 1.179 ± 0.024 | 55.14 ± 2.90 | 0.013 ± 0.006 |
| Дистальная фаланга | 1.191 ± 0.049 | 65.95 ± 1.08 | 0.013 ± 0.005 |
| Тело позвонка | 1.194 ± 0.037 | 66.14 ± 1.03 | 0,008 ± 0,003 |
| Верхняя позвоночная дуга | 1.19 ± 0.017 | 63.02 ± 1.23 | 0.021 ± 0.009* |
| Осыпь | 1.198 ± 0.010 | 68.20 ± 1.24 | 0.011 ± 0.003 |
| *Образец KRA005 удален как выброс при 0,1952 | | | |
Доступность кода
Все аналитические программы и модули R, используемые при анализе этой рукописи, находятся в свободном доступе у их соответствующих авторов. Весь пользовательский код R доступен по запросу.
Доступность данных
Все исходные данные, используемые при расчете репрезентативных результатов, находятся в свободном доступе в репозитории данных ENA Европейского нуклеотидного архива (номер присоединения PRJ-EB36983) или дополнительных материалах Parker, C. et al.11.
Дополнительный файл 1. Пожалуйста, нажмите здесь, чтобы загрузить этот файл.