🔬 المنهجية العلمية

شرح مفصّل وصادق لكل خطوة علمية نستخدمها. لا صندوق أسود.

خط المعالجة تحديد السلالة T2T-CHM13 جودة التغطية تصنيف SNP عمق كل SNP TMRCA الحدود الصادقة
1خط المعالجة الكامل

عند رفع ملف BAM أو VCF، يمر عبر هذه المراحل بالترتيب:

Yleaf v4.02تحديد الطفرات المعروفة ومقارنتها بشجرة YFull v10
Connected-path algorithmخوارزميتنا لتحديد الفرع الطرفي (القسم 2)
HipSTRاستخراج Y-STR markers من ملف BAM
pyliftoverتحويل كل إحداثية من hg38 إلى T2T-CHM13v2.0
samtoolsحساب جودة التغطية والعمق لكل SNP
NIST/GIAB masksتصنيف كل SNP بالموثوقية
2كيف نحدد سلالتك

لا نختار ببساطة أعمق فرع في الشجرة. نستخدم خوارزمية من خطوتين:

الخطوة 1: تحديد العائلة الرئيسية بعدد الطفرات الأكبر
الخطوة 2: المشي من ROOT عبر فروع مثبتة فقط، حتى ننقطع عند أول فرع بلا دليل
⚠️ لماذا هذا مهم
نظام أبسط قد يختار فرعاً عميقاً بدليل ضعيف (SNP واحد) بدل فرع أعلى بدليل قوي. خوارزميتنا تعطي الأولوية للدليل الأقوى.

→ شاهد مثالاً حياً: أدلة تصنيف السلالة

3T2T-CHM13 — الميزة الحصرية

مرجع hg38 يحتوي فجوات كبيرة في كروموسوم Y — أكثر من 50% كان مفقوداً حتى 2022. T2T-CHM13v2.0 أكمل هذه الفجوات.

نحوّل كل SNP تلقائياً بين المرجعَين عبر UCSC liftOver. الموقع الحقيقي قد يتغيّر بمئات الآلاف من القواعد.

مصدر البيانات
Nurk et al. 2023, "The complete sequence of a human Y chromosome", Nature
4جودة التغطية

لملفات BAM، نحسب إحصائيات تغطية chrY الكامل باستخدام samtools:

samtools coverage -r chrY sample.bam samtools depth -r chrY -q 20 -Q 30 sample.bam | awk '$3>=5'
⚠️ فرق مهم عن YFull
نحسب على chrY الكامل (57.2Mb). YFull يقتصر على MSY (~10.4Mb) بعد استبعاد المناطق المتكررة. الأرقام غير قابلة للمقارنة المباشرة.
5تصنيف الطفرات: Clean مقابل H

كروموسوم Y يحتوي مناطق متكررة (ampliconic) تظهر عدة مرات. طفرة فيها قد تكون قراءة من نسخة مختلفة، وليست طفرة فريدة.

نستخدم ملفات BED رسمية من NIST/GIAB Genome Stratifications (v3.5):

التصنيفالمعنىالمصدر
Cleanمنطقة عالية الموثوقيةخارج المناطق المقنّعة
Hمنطقة متكررة (Homologous)GRCh38_chrY_ampliconic.bed.gz + GRCh38_chrY_XTR.bed.gz
⚠️ ملاحظة صادقة
طفراتنا «النظيفة» ليست كلها بالضرورة جديدة فعلاً. YFull يقارن مع شجرة أضخم بكثير فيصل لعدد أدق.

→ مصدر البيانات: NIST/GIAB Genome Stratifications (GitHub)

6عمق القراءة لكل طفرة

لكل private SNP في ملف BAM، نحسب عدد القراءات الداعمة عند ذلك الموضع:

samtools depth -b positions.bed -q 20 sample.bam

نستخدم فلتر جودة المحاذاة (≥20) فقط. فلتر جودة القاعدة الإضافي (≥30) كان يُخفي قراءات حقيقية موثوقة، فقرّرنا عدم استخدامه.

العمقالتصنيف اللوني
≥ 20x أخضر — موثوق جداً
10–19x أصفر — مقبول
< 10x أحمر — ضعيف
7حساب TMRCA

عندما يوجد مستخدمان أو أكثر بنفس السلالة، نحسب TMRCA بمعادلة YFull القياسية:

TMRCA = (private_SNPs_shared × 144.41) + 60
⚠️ متى لا نحسب TMRCA
إذا كنت المستخدم الوحيد في سلالتك، نعرض «غير متاح» بصراحة. الدقة تزداد مع رفع عيّنات أكثر.
8حدودنا الصادقة مقابل YFull
الجانبDeepAncestry™YFull
حجم قاعدة SNPs المرجعيةصغيرة (محلية)ضخمة (آلاف العيّنات، 10+ سنوات)
مراجعة يدوية للشجرةآلية بالكاملفريق علمي متخصص
T2T-CHM13حصري✗ hg38 فقط
شفافية المنهجيةهذه الصفحةمحدودة

نحن لا ندّعي أننا أدق من YFull في الفروع النادرة أو TMRCA المعقّد. نقدّم شفافية كاملة وميزات مختلفة، وليس بديلاً كاملاً.

جرّب بنفسك