نسبة Ka/Ks مفصلة

خطوات الاستخدام

الخطوة 1: إدخال تسلسلين من الكودونات (CDS) بصيغة FASTA أو نص عادي. يقوم الحاسب بمحاذاة التسلسلين باستخدام خوارزمية مثل Needleman-Wunsch لضمان مقارنة دقيقة لكل موقع.

الخطوة 2: تصنيف كل موقع كـ 'متزامن' إذا كان التغيير في النوكليوتيد لا يغير الحمض الأميني، أو 'غير متزامن' إذا غيّر الحمض الأميني. يُستخدم الجدول الوراثي القياسي لتحديد ذلك.

الخطوة 3: حساب عدد المواقع المتزامنة (S) وغير المتزامنة (N) المتاحة للطفرة، مع مراعاة الشفرة الوراثية المتنوعة (مثل الشفرة العالمية أو الخاصة بالميتوكوندريا).

الخطوة 4: حساب معدل الطفرات المتزامنة (Ks = عدد الطفرات المتزامنة / عدد المواقع المتزامنة) ومعدل الطفرات غير المتزامنة (Ka = عدد الطفرات غير المتزامنة / عدد المواقع غير المتزامنة). يتم استخدام تصحيح بواسون أو نموذج Jukes-Cantor لتعويض الطفرات المتعددة.

الأسئلة الشائعة

نسبة Ka/Ks هي مقياس لقوة ونوع الانتقاء الطبيعي على مستوى تطور البروتين. Ka هو معدل التبديلات النوكليوتيدية غير المتزامنة (التي تغير الحمض الأميني) لكل موقع غير متزامن، وKs هو معدل التبديلات المتزامنة (التي لا تغير الحمض الأميني) لكل موقع متزامن. بما أن التبديلات المتزامنة تُعتبر محايدة غالباً، فإن Ks يعمل كمرجع لمعدل الطفرة المحايد. إذا كانت النسبة > 1، فهذا يعني أن التبديلات غير المتزامنة تتراكم بمعدل أعلى من المتوقع تحت التحييد، مما يشير إلى انتقاء إيجابي (تطور تكيفي). إذا كانت النسبة < 1، فهذا يعني أن التبديلات غير المتزامنة يتم إزالتها بواسطة الانتقاء التطهيري (السالب). إذا كانت النسبة ≈ 1، فهذا يشير إلى تطور محايد. ومع ذلك، النسبة المتوسطة على طول الجين بالكامل قد تخفي انتقاءً إيجابياً في مواقع محددة، لذا يُفضل استخدام تحليل النافذة المنزلقة أو نماذج الكودون. من المهم أيضاً استخدام تصحيحات الطفرات المتعددة (مثل Jukes-Cantor أو Kimura) لتجنب تقدير أقل من الواقع للمسافة الحقيقية.

هناك عدة طرق: (1) طريقة Nei-Gojobori: وهي الأبسط وتعتمد على حساب متوسط عدد المواقع المتزامنة وغير المتزامنة عبر جميع الكودونات، ثم تقدير المسافة باستخدام تصحيح بواسون أو Jukes-Cantor. (2) طريقة Li-Wu-Luo: تحسن دقة تقدير المواقع المتزامنة وغير المتزامنة عن طريق تصنيف كل موضع في الكودون (أول، ثان، ثالث) نظراً لاختلاف معدلات الطفرة. (3) طريقة Yang-Nielsen: تستخدم نماذج تطور الكودون (مثل نموذج Goldman-Yang) وتقدر المعاملات عبر الاحتمال الأقصى (ML). هذه الطريقة الأكثر دقة وتعقيداً، وتسمح بتقدير النسبة لكل موقع. (4) طرق بايزية: مثل التي تستخدم في برنامج PAML (yn00)، والتي تدمج معلومات مسبقة عن المعدلات. لكل طريقة مزايا وعيوب؛ الطرق البسيطة سريعة لكنها قد تكون متحيزة في ظل اختلاف كبير في التركيب النوكليوتيدي أو اختلاف في نسب الانتقال/التحول.

تحليل النافذة المنزلقة هو تقنية لفحص التباين المحلي لنسبة Ka/Ks على طول الجين. بدلاً من حساب نسبة واحدة للجين بأكمله، نقسم الجين إلى نوافذ صغيرة (مثلاً 10-50 كودوناً) ونحسب النسبة لكل نافذة، ثم نزلق النافذة بخطوة (مثلاً 1-10 كودونات) للحصول على أطياف متصلة. هذا يساعد في تحديد مناطق الجين التي قد تخضع لانتقاء إيجابي (نسبة >1 في نافذة) حتى لو كان متوسط الجين أقل من 1. على سبيل المثال، في جين ليزوزيم القردة العليا، وجدت دراسات أن بعض النوافذ في الموقع النشط أظهرت نسبة >1، مما يشير إلى تطور تكيفي في تلك المنطقة. لتنفيذ ذلك في الحاسبة، يتم إدخال حجم النافذة والخطوة، ويقوم الحاسب بحساب Ka و Ks لكل نافذة باستخدام نفس الطريقة المختارة (مثل ML أو Nei-Gojobori). النتائج تعرض كجدول أو رسم بياني يوضح النسبة عبر المواقع.

الشفرة الوراثية تختلف بين النوى والميتوكوندريا والبعض الآخر. في الميتوكوندريا، بعض الكودونات ترمز لأحماض أمينية مختلفة (مثل AUA ترمز للميثيونين بدلاً من الأيزوليوسين). هذا يؤثر على تصنيف المواقع المتزامنة وغير المتزامنة. على سبيل المثال، طفرة بين AUA وAUG: في الشفرة النووية، AUA ترمز للأيزوليوسين وAUG للميثيونين، لذا التغيير غير متزامن. لكن في الشفرة الميتوكوندرية البشرية، كلاهما يرمزان للميثيونين، لذا تكون متزامنة. لذلك، يجب تحديد الشفرة الوراثية المستخدمة قبل الحساب. في الحاسبة التفصيلية، نسمح للمستخدم باختيار الشفرة الوراثية (قياسية، ميتوكوندرية حيوانات، ميتوكوندرية خميرة، إلخ). كذلك، بعض الشفرات الوراثية لديها كودونات توقف مختلفة، مما قد يؤثر على عد المواقع الصالحة. تجاهل هذه الاختلافات قد يؤدي إلى نتائج خاطئة، خاصة في دراسات تطور الجينات الميتوكوندرية.

نسبة Ka/Ks مفصلة

أدخل القيم

خطوات الاستخدام

أمثلة محلولة

أمثلة محسوبة

الأسئلة الشائعة

أدخل القيم

خطوات الاستخدام

1مثال: مقارنة تسلسلين لجين P53 بين نوعين من الثدييات. التسلسل الأول: ATG... (طول 10 كودونات) والتسلسل الثاني: ATG... مع بعض الاختلافات. احسب نسبة Ka/Ks وفسر النتيجة.

2مثال آخر: مقارنة جين lysozyme في الإنسان والقردة العليا. طول التسلسل 50 كودوناً. وجد 10 تغيرات غير متزامنة و5 تغيرات متزامنة. احسب النسبة وفسر.

أمثلة محسوبة

الأسئلة الشائعة