$$\rightleftharpoonup{xx}$$
$$\longleftharp{xx}$$,
$$\longrightharp{xx}$$,
במחקר נפרד מס' 11, דנ"א הופק מאבקת עצם שנוצרה מכל מיקום דגימה אנטומיב-11 פרטים, תוך שימוש בפרוטוקול מיצוי דנ"א סטנדרטי הממוטב למקטעים קצרים מרקמה2. ספריות חד-גדיליות יוצרואז 28 ורוצפו על HiSeq 4000 (75 bp זוגי-קצה) לעומק של ~20,000,000 קריאות לכל דגימה. נתוני הרצף שהתקבלו הוערכו לאחר מכן עבור תוכן דנ"א אנושי אנדוגני באמצעות צינור EAGER29 (הגדרות BWA: אורך זרע של 32, עונש אי התאמה 0.1, מסנן איכות מיפוי של 37). כל התוצאות המייצגות מדווחות באמצעות אותם מדדים כמו Parker et al. 202011 לצורך עקביות. ספריות מהחלקים האבקתיים של ה-pars petrosa הניבו, בממוצע, דנ"א אנדוגני גבוה יותר מכל אחד מ-23 אתרי הדגימה האנטומיים האחרים שנסקרו (איור 6A-B). שבעת מיקומי הדגימה האנטומיים הנוספים המוצגים בפרוטוקול זה (הצמנטום, המעבר הראשון של תא מוך השן ודנטין של טוחנות קבועות; עצם קליפת המוח מגוף החוליה וקשת החוליות העליונה של חוליית החזה; עצם קליפת המוח מהציצית האפית של הפלנקס הדיסטלי; ועצם קליפת המוח מצוואר הטלוס) הניבו את התשואות הגבוהות הבאות (ללא מובהקות סטטיסטית בין מיקומי דגימה אנטומיים אלה; איור 6A-B; קובץ משלים 1: אנדוגניDNAPreCap). כל המיקומים החלופיים האלה מייצרים באופן עקבי תפוקות דנ"א המתאימות לניתוחים גנטיים סטנדרטיים של אוכלוסיות, כגון ניתוחים מיטוכונדריאליים וניתוחי פולימורפיזם של נוקלאוטידים בודדים (SNP). שיעורי הכפילויות בספריות שנבעו מכל מיקומי הדגימה האנטומיים היו נמוכים (גורמי אשכול < 1.2 בממוצע, המחושבים כיחס בין כל קריאות המיפוי לקריאות מיפוי ייחודיות, טבלה 2; קובץ משלים 1: ClusterFactor), המציין כי כל הספריות שנבדקו היו בעלות מורכבות גבוהה מאוד. באופן דומה, הערכות הזיהום הממוצעות של דנ"א אנושי אקסוגני היו נמוכות, בממוצע < 2% (זיהום כרומוזום X אצל זכרים, n = 7, כפי שדווח על ידי צינור ANGSD30) בכל אתרי הדגימה האנטומיים למעט קשת החוליות העליונה (זיהום משוער ממוצע: 2.11%, כאשר דגימה אחת הוסרה כחריגה; KRA005: 19.52%, ראה לוח 2; קובץ משלים 1: Xcontamination). אורך השבר הממוצע (לאחר סינון להסרת כל הקריאות < 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: רמות שכפול ממוצעות (מיפוי קריאות/קריאות ייחודיות), אורכי מקטעים ממוצעים וחציוניים ואומדני זיהום כרומוזום X עבור כל מיקומי הדגימה האנטומית. שגיאה שדווחה כשגיאת התקן של הממוצע. טבלה זו הותאמה מתוך Parker, C. et al. 202011.
| מיקום הדגימה | מקדם שכפול ממוצע (# קריאות ממופות /# קריאות ממופות ייחודיות) | אורך מקטע ממוצע ב- bp | שיעור משוער ממוצע של זיהום כרומוזום X |
| פירמידת פטרוס | 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 מותאם אישית זמין על פי בקשה.
זמינות נתונים
כל הנתונים הגולמיים המשמשים לחישוב התוצאות המייצגות זמינים באופן חופשי במאגר הנתונים האירופי Nucleotide Archive ENA (מספר הצטרפות PRJ-EB36983) או בחומרים משלימים של Parker, C. et al.11.
קובץ משלים 1. אנא לחץ כאן כדי להוריד קובץ זה.