מידע חדש - השערה שייכת לא"ח

בפוסט הקודם שלי ניתחתי את המשמעות הסטטיסטית של ממצא השערה של א”ח. מצאתי שיש לו מספר מוטציות נדירות שבהערכה שמרנית מביאות את עוצמת ההתאמה ל-1:100,000.

מאז קיבלתי מידע נוסף, ולאחר התייעצות עם מומחה הגנטיקה המאוד חביב דביר דהרי, אפשר להגיד שהממצא חזק אפילו יותר, והשערה בוודאות שייכת לא”ח.

המידע החדש:

1. מאגרים חדשים - אני והמכון לרפואה משפטית התבססנו על מאגר EMPOP, שכולל כ-5000 דגימות. הוא נותן את אפשרויות החיפוש הכי מתקדמות להתאמת רצפים, אבל אם רוצים לאמוד נדירות של מוטציות ספציפיות, יש מאגרים הרבה יותר גדולים: MITOMAP עם כ-57,000 דגימות, ו-Helix עם כ-200,000 דגימות.

2. מיפוי פילוגנטי של המיטוכנודריה האנושית שמאפשר להבין אילו מוטציות התקבצו יחד לאורך האבולוציה.

הממצאים החדשים:

1. מתוך 37 המוטציות של א”ח, 30 הן בדיוק המוטציות של Haplogroup I, שנוצרה באיזור אירן לפני 20,000 שנה.
   לפני שבוע, כשעוד היה לנו רק את המיפוי החלקי של השערה והערכת התפוצה הייתה 1 ל-1000 בעולם (ע”פ EMPOP), ו-1 ל-222 בארץ (שאח”כ עודכן ל-1 ל-700 עם הגדלת המאגר) - אז זה בעצם התייחס לתפוצה של השושלת הזו.

2. עוד 2 מוטציות (309.1C 315.1C) מגיעות בקורלציה מאד גבוהה לקודמות, ואפשר להתייחס אליהן כחלק מאותה שושלת I.

3. נשארות 5 מוטציות, שלא מופיעות בשושלת I, או באף שושלת ידועה שנגזרת ממנה.

4. 2 הראשונות, שהן בעצם מוטציה אחת (מחיקת שני מיקומים צמודים 523del 524del), מאד נפוצות בכלל האוכלוסיה, אבל נדירות בקבוצה I. ב-EMPOP יש 25 משושלת I, ולאף אחד מהם אין את המוטציה הזו.

5. השלוש שנשארו הן המוטציות הנדירות מהפוסט הקודם, אבל עכשיו יש לנו מידע יותר מדוייק על התפוצה שלהן, שמסוכם בטבלה הבאה:

   

6. אז 2120A נשארת בערך באותה תפוצה. קודם הייתה 1 ל-1000, עכשיו זה 1 ל-2000 או 1 ל-3000.

7. 14484G שלא הופיעה ב-EMPOP, לא מופיעה גם ב-MITOMAP, אבל מופיעה פעם אחת ב-Helix. אז אפשר לתת לה הערכה יותר מדוייקת של כ-1:250,000 (!).
   ההופעה הזו לא קשורה לשושלת I (הדגימה שהתאימה היא משושלת H), אז גם פה אין התאמה מלאה לא”ח.

8. 5444T שהייתה קודם בתפוצה של 1 ל-2000, הרבה יותר נדירה במאגרים האחרים: 1 ל-50,000 בשניהם יחד. כשבודקים את מקור הדגימות ב-EMPOP מגלים ששתיהן מגיעות מאותה שושלת ספציפית במזרח טימור, וזה מה שגרם שם להערכת היתר של התפוצה.

אז בעצם יש לנו ארבעה גורמים מאוד נדירים שונים בדנ”א המיטוכונדרי של א”ח: כל אחת משלוש המוטציות הנדירות לבדן, והשילוב של שושלת I עם המוטציות 523del 524del.

למי שלא רוצה להסתבך, הממצא הכי פשוט להבנה הוא שיש מוטציה אחת (14484G) שהיא לבדה כל כך נדירה שזה יותר ממספיק כדי להכריע שהשערה שייכת לא”ח בוודאות, בהינתן שהוא כבר נקשר לזירה לפני-כן.

להמחשה אינטואיטיבית: זה כמו למצוא בזירה קרע של כרטיס ביקור ורואים בו 8-1648 ואז מתגלה שהטלפון של אדם שכבר היה לגביו חשד מסויים הוא 052-538-1648.

אבל בשביל הספורט, נבדוק אם אפשר להגיד יותר מזה. השאלה המעניינת היא מה הקורלציה בין ארבעת הגורמים השונים. אם הם בלתי תלויים אז מכפלת הסיכויים שלהם היא במיליארדים, וברור שהם יימצאו רק אצל א”ח ועוד מספר קטן של קרובי משפחה שלו מהשושלת האמהית.

בפוסט הקודם מצאתי ש-2120A ו-5444T לא מופיעים יחד במאגר EMPOP, אז הן כנראה לא מאוד קשורות. ב-Helix לא ראיתי אופציה לבדוק מופעים משותפים, אבל ב-57,000 הדגימות של MITOMAP אפשר, ואכן גם שם אין אף מופע משותף לזוג הזה.

עוד בדקתי את התפוצה של המוטציות הבודדות של שושלת I בשלוש אופנים: א) לבד, ב) עם 2120A, ג) ועם 5444T, וכמו שרואים בגרף, אין קשר חזק - כלומר, אם למוטציה מסויימת יש תפוצה X באוכלוסיה הכללית, אז יש לה תפוצה X גם בקרב בעלי 2120A או 5444T. זה רחוק מחד-משמעי, אבל נותן אינדיקציה מסויימת שהמוטציות האלה לא קורלטיביות לשושלת I (בדומה ל-324del 325del, שהייתה אפילו בקורלציה הפוכה לשושלת I).

המסקנה הסבירה מכל זה: א”ח הוא משושלת I שהתפתחה לפני כ-20,000 שנה. לאורך אלפי השנים האחרונות נוספו בשושלת האמהית שלו, באופן אקראי ובלתי-תלוי, עוד ארבע מוטציות. סטטיסטית, סביר שהאחרונה שבהם נוצרה במאות השנים האחרונות, שזה אומר שהרצף שלו יכול להתאים רק לכמה עשרות אנשים בכל העולם (לא להתבלבל עם מספר הצאצאים שעולה מעריכית - השושלת האמהית גדלה בקצב איטי כי רק הבנות מעבירות את הגנים).

מי שרוצה להמשיך לטעון שזו ראייה עם “משקל מוגבל”, כלומר שיש עוד עשרות אלפי ישראלים עם הרצף הזה, חייב לטעון ש:

1. רוב השושלת הזו בארץ - כי אנחנו יודעים ש-14484G לבדה היא מקסימום עשרות אלפי אנשים בעולם.

2. השושלת הזו מרכיבה את רוב שושלת I בארץ - כי אנחנו יודעים שכל שושלת I בארץ היא אלפים או עשרות אלפים.

מעבר לכך שזה דורש כמה צירופי מקרים נדירים, זה גם מנוגד למה שאנחנו יודעים על שושלות אמהיות. שושלות כאלה נוצרות לאורך אלפי ועשרות-אלפי שנים. אלה קבועי זמן שלא רלבנטים לשושלות יהודיות, בטח לא ישראליות. בנוסף, אנחנו יודעים על עירוב בין יהודים למקומיים בכל מקום בו שהו. לכן אנחנו רואים פיזור גאוגרפי די רחב של כל השושלות האמהיות (יש כמובן הטיות גאוגרפיות, אבל לא באופן מוחלט). יוצאים מן הכלל אלה אוכלוסיות מבודדות באיים, כפי שראינו למעלה במזרח טימור, וגם שם רואים הגירה בין האיים.

כלומר, אין שום אפשרות שיש בישראל שושלת גדולה שעברה מתחת לרדאר של כל המחקרים העולמיים.

האופציה היחידה שבה השערה היא לא של א”ח זה במקרה של טעות מעבדה חמורה, או הונאה. מישהו רואה אפשרות נוספת?