في علم الوراثة السكانية، تُدرس إحصاءات ترددات الأليلات عادةً في ظل انتقاء ضعيف أو حيادية. لكن في كثير من السيناريوهات البيولوجية، مثل تكيّف الفيروسات والبكتيريا أو الانتقاء الاصطناعي، يكون الانتقاء قوياً بالنسبة إلى الانحراف العشوائي. وتقع إلى جانب هذا العمل دراسات إما تركز على النتائج بعيدة المدى تحت انتقاء قوي، وإما تستخدم مساراً حتمياً واحداً. وتسدّ هذه الورقة تلك الفجوة: فهي تتناول إحصاءات انتقالية (معتمدة على الزمن) في نظام الانتقاء القوي، حيث تعطي التقريبات المعتادة أخطاء كبيرة بسبب تجاهل المسارات المؤدية إلى فقدان الأليل.
يعمل المؤلف في إطار تقريب الانتشار لنموذج رايت-فيشر. يُعرَّف المعامل R = 2Ne|s|، حيث Ne هو الحجم الفعّال للسكان، وs معامل الانتقاء. وعندما يكون R >> 1 يكون الانتقاء قوياً. ويُمثَّل إحصاء تردد الأليل كمجموع مرجَّح لتوقعين مشروطين: بشرط التثبيت، وبشرط الفقدان. ثم يُعاد كتابة كل توقع مشروط عبر عملية معدَّلة بانتفاء مختلف. بعد ذلك يُطبَّق تفكيك وفق المعامل الصغير ε = 1/R، حيث √ε هو سعة الضجيج (الانحراف). ويحتفظ التقريب بحدود حتى الرتبة ε، ما يعطي إسهاماً من مسارين حتميين ومن التقلبات حولهما. وتُحل معادلات المسارات الحتمية وتصحيحات التقلبات عددياً.
تم الحصول على تعبيرات صريحة للمتوسط والتباين في تردد الأليل كدالتين للزمن. ويُقرَّب المتوسط بمجموع مسارات حتمية مرجَّحة وتصحيح من الرتبة 1/R. ويحتوي التباين على حد ينشأ من اختلاف المسارين الحتميين، وحد من التقلبات. وتُظهر المقارنة العددية مع النتائج الدقيقة لنموذج رايت-فيشر أن التقريب الحتمي المحض قد يعطي خطأً يصل إلى 100% في التباين، خصوصاً عندما يكون التردد الابتدائي صغيراً. أما إدخال تصحيح التقلبات فيخفض الخطأ إلى بضعة في المئة. فعلى سبيل المثال، عند R=50 وتردد ابتدائي y=0.001 ينخفض خطأ المتوسط من 10.6% إلى 3.2%، وخطأ التباين من 12.6% إلى 4.9%. ويكون إسهام التقلبات مهماً بشكل خاص في الأزمنة المبكرة، حين يكون التباين الحتمي شبه معدوم.
يُظهر العمل أنه حتى في ظل انتقاء قوي لا يجوز إغفال المسارات المؤدية إلى فقدان الأليل. وخلافاً للدراسات المجاورة التي تستخدم مساراً حتمياً واحداً أو تتناول الحدود بعيدة المدى فقط، يُقدَّم هنا أسلوب منهجي للإحصاءات الانتقالية. وهذا يسدّ فجوة مهمة: فبات من الممكن الآن حساب الديناميكا الزمنية لترددات الأليلات حساباً سليماً في نظام الانتقاء القوي، وهو ما كان متعذراً بالأساليب المعتادة. ويكشف الأسلوب أيضاً أن التقلبات ليست دائماً تصحيحات صغيرة: فقد تهيمن على التباين في فترات زمنية معينة.
يمكن توسيع الأسلوب ليشمل حالة أليلات متعددة، حيث يلزم الجمع على تثبيتات أليلات مختلفة. كما يمكن إدخال الطفرات، وهو ما يتطلب مقاربة أخرى، لأن التوزيع المستقر لا يُوصف بمسارين حتميين. ويظل التعميم على حجم سكان متغير وانتقاء متغير موضع اهتمام أيضاً. وإلى جانب ذلك، يمكن دراسة السلوك في الطبقات الحدودية قرب الترددين 0 و1، حيث قد تتطلب التفكيكات المقاربة معالجة مستقلة.
الأسلوب قابل للتطبيق في تحليل السلاسل الزمنية لترددات الأليلات في التطور التجريبي، وفي مجتمعات الفيروسات والبكتيريا، وكذلك في برامج الانتقاء الاصطناعي. وقد يحسّن تقدير معاملات الانتقاء من البيانات، ويساعد في التنبؤ بمصير الأليلات المفيدة في المجموعات السكانية الصغيرة، وهو أمر مهم لعلم الوراثة الخاص بحفظ الطبيعة.
قد تكون الخطوة التالية تطبيق الأسلوب على بيانات حقيقية لسلاسل زمنية من ترددات الأليلات، للتحقق من فعاليته في تقدير معاملات الانتقاء. ويجدر أيضاً بحث التوسيع ليشمل حالة أليلات متعددة وإدخال الطفرات.
لا تزال مشكلة الوصف الدقيق للديناميكا الانتقالية لترددات الأليلات تحت انتقاء قوي مفتوحة في علم الوراثة السكانية، خصوصاً بالنسبة إلى التباين والعزوم الأعلى. ويعالج هذا العمل هذه المشكلة مباشرةً في حالة أليلين دون طفرات، مقترحاً أسلوباً فعّالاً حسابياً. كما يطرح مسألة قابلية تطبيق التفكيكات المقاربة قرب حدود نطاق التردد.
حتى في ظل انتقاء قوي، حين تكون القوة الحتمية كبيرة، قد يتحدد تباين تردد الأليل تقريباً بالكامل بتقلبات عشوائية على مدى زمني طويل. فعلى سبيل المثال، عند R=50 وتردد ابتدائي y=0.001 يكون الإسهام الحتمي في التباين شبه معدوم حتى الزمن τ≈5، بينما يصف إسهام التقلبات التباين الحقيقي وصفاً دقيقاً.