Seasonal fluctuations in fitness result in severe reductions in effective population size

1. بيان المشكلة

يُعد حجم الجماعة الفعال ($N_e$) معلماً أساسياً في علم الوراثة الجماعية، حيث يحدد قوة الانحراف الوراثي والحفاظ على التنوع الجيني. وبينما يُنظر إلى $N_e$ غالباً من منظور ديموغرافي، فإنه يتأثر بشدة بالانتخاب الطبيعي. وقد أثبتت الأبحاث السابقة أن الانتخاب الإيجابي (المسحات الانتخابية) والانتخاب المنقي (الانتخاب الخلفي) يقللان من $N_e$، بينما يمكن للانتخاب المتوازن أن يزيده.

ومع ذلك، فإن تأثير الانتخاب المتقلب — حيث تتغير ضغوط الانتخاب في الاتجاه أو الكثافة بمر over الزمن (على سبيل المثال، موسمياً) — على $N_e$ على مستوى الجينوم لا يزال غير مفهوم بشكل جيد. لقد حددت الدراسات التجريبية في ذبابة الفاكهة (Drosophila melanogaster) مئات المواقع التي تخضع لتذبذبات تكيفية عبر المواسم، ومع ذلك اعتمدت النماذج النظرية إلى حد كبير على افتراضات مجردة أو ضعيفة المعلمات. ومن غير المعروف كيف يجمع التأثير التراكمي لمئات المواقع المتقلبة موسمياً، وتفاعلاتها فوق الجينية (epistatic interactions)، والعوامل الديموغرافية، لتغيير $N_e$ في الجماعات الطبيعية.

2. المنهجية

استخدم المؤلفون محاكاة التطور الزمني الأمامي باستخدام برنامج SLiM (v4.0.1) لنمذجة التكيف متعدد المواقع تحت تأثير التحولات البيئية الموسمية.

• نموذج الرفع الانعزالي (Segregation Lift - SL): استخدمت الدراسة نموذجاً تم تطويره مؤخراً (Wittmann et al., 2017) حيث يتم تحديد اللياقة الفردية من خلال درجة موسمية ($z$). هذه الدرجة هي مجموع المساهمات من مواقع تكيفية متعددة. ويتم حساب اللياقة ($w$) عبر تحويل القوة ($w = (1+z)^y$)، حيث يتحكم الأس ($y$) في شدة التفاعل فوق الجيني (تحديداً، تفاعل العوائد المتناقصة).
• معلمات المحاكاة:
  • الديموغرافيا: أُجريت عمليات المحاكاة مع حجم جماعة تعداد ثابت ($N_c$) يبلغ مليون فرد ثنائي الصيغة الكروموسومية (تم خفضه إلى 500 ألف لأغراض الكفاءة الحسابية في بعض التجارب) على مدار 36,000 جيل (ما يعادل 18,000 جيل في النموذج المخفض).
  • الموسمية: تتبادل البيئة بين الصيف والشتاء، بمعدل 10 أجيال لكل موسم.
  • المواقع: اختبرت المحاكاة الأولية من 100 إلى 500 موقع موسمي مع تباين في شدة التفاعل فوق الجيني ($y$ تتراوح من 0.5 إلى 20). بدأت الترددات الأليلية من 50%.
  • المعايرة التجريبية: لمحاكاة حالة D. melanogaster الطبيعية، استخلص المؤلفون معلمات محددة (عدد المواقع وشدة التفاعل فوق الجيني) باستخدام نماذج الانحدار الخطي المتعدد المطبقة على بيانات من دراسات تجريبية حديثة (مثل Bergland et al. 2014; Machado et al. 2021; Rudman et al. 2022).
• تقدير $N_e$: قامت الدراسة بحساب حجم الجماعة الفعال للتباين ($N_e$) عن طريق تتبع التباين الزمني في الترددات الأليلية لـ العلامات المحايدة (التي تم إدراجها على فترات 100 كيلوبايت على أذرع كروموسومية غير مرتبطة) عبر دورة موسمية كاملة. سمح هذا للمؤلفين بقياس التأثير الجينومي للانتخاب على التنوع المحايد المرتبط.
• اختبارات المتانة: اختبر المؤلفون سيناريوهات مختلفة، بما في ذلك:
  • تغيير أحجام تعداد الجماعة ($N_c$ من 20 ألف إلى مليون).
  • ديموغرافيا "الازدهار والانهيار" الموسمية (تذبذب حجم الجماعة بين مليون في الصيف و100 ألف في الشتاء).
  • وضع حد لعدد النسل (تحديد التباين التكاثري إلى 20 نسلاً لكل والد) لمعالجة العيوب المحتملة في نموذج SL.

3. المساهمات الرئيسية

• تكميم تأثير الانتخاب المتقلب: تعد هذه أول دراسة تقدم تقديرات مباشرة ومستندة إلى بيانات تجريبية حول كيفية تقليل الانتخاب الموسمي متعدد المواقع لـ $N_e$ في سياق جينومي واقعي.
• النمذجة التنبؤية: طور المؤلفون نماذج انحدار تتنبأ بدقة بسعة التردد الأليلي بناءً على عدد المواقع المتعددة وشدة التفاعل فوق الجيني ($R^2 > 0.96$).
• تحديد المحرك الأساسي: أظهروا أن النقص في $N_e$ لا ينتج عن متوسط تقلب جميع المواقع، بل يرتبط ارتباطاً قوًا بـ أقصى سعة تردد أليلي للموقع الواحد الأكثر تقلباً.
• الفصل عن حجم التعداد: تُظهر الدراسة أن النقص النسبي في $N_e$ مستقل إلى حد كبير عن حجم تعداد الجماعة ($N_c$)، مما يتحدى الفكرة القائلة بأن الانتخاب المتقلب وحده يمكن أن يفسر "مفارقة ليوونتن" (الارتباط الضعيف بين $N_e$ و $N_c$ عبر الأنواع).

4. النتائج الرئيسية

• انخفاض حاد في $N_e$: في السيناريوهات الواقعية تجريبياً لـ D. melanogaster، يقلل الانتخاب الموسمي المتقلب من حجم الجماعة الفعال بنسبة تقارب 50% مقارنة بحجم التعداد ($N_e/N_c \approx 0.5$). تراوحت الانخفاضات من 7% إلى 76% اعتماداً على مجموعة المعلمات المحددة.
• دور السعة القصوى: وُجد ارتباط خطي قوي ($R^2 = 0.702$) بين أقصى سعة تردد أليلي موسمية وحجم النقص في $N_e$. المواقع ذات التقلبات الأكثر تطرفاً (على سبيل المثال، سعة ~35-45% الملاحظة في الطبيعة) هي التي تقود تأثير عنق الزجاجة.
• الآلية: ينتج النقص عن عنق زجاجة تكاثري متسلسل. مع تغير المواسم، يصبح الأفراد ذوو الأنماط الجينية المفضلة في الموسم السابق غير متكيفين، بينما يكتسب الأفراد الذين كانوا غير متكيفين سابقاً ميزة تكاثرية. هذا يخلق تبايناً عالياً في النجاح التكاثري (عدد النسل) عند نقاط الانتقال بين المواسم، مما يؤدي فعلياً إلى تقليص $N_e$.
• الاستقرار: تستقر تذبذبات التردد الأليلي بسرعة (خلال ~1,200-2,000 جيل)، ويستقر النقص الناتج في $N_e$ بعد ذلك بوقت قصير.
• الاستقلال الديموغرافي: ظل النقص النسبي في $N_e$ ثابتاً عبر مختلف أحجام التعداد الثابتة. وحتى في سيناريوهات "الازدهار والانهيار"، كان النقص يتناسب مع المتوسط التوافقي لحجم الجماعة، محافظاً على نسبة النقص (~50%) بالنسبة للحجم الفعال.
• تحديد عدد النسل: أدى وضع حد أقصى لعدد النسل لكل والد (لمنع التباين غير الواقعي) إلى تقليل طفيف في حجم فقدان $N_e$ في السيناريوهات المتطرفة، لكنه لم يلغِ النقص بنسبة ~50% الملاحظ في السيناريوهات المتوسطة.

5. الأهمية

• إعادة تقييم التنوع الجيني: تشير النتائج إلى أن النماذج القياسية للتنوع الجيني قد تبالغ في تقدير $N_e$ بشكل كبير في الأنواع التي تخضع لانتخاب موسمي قوي (مثل الحشرات وغيرها من ذوات الدم البارد). إن التباين "الخفي" في النجاح التكاثري الناتج عن التكيف الموسمي هو قيد رئيسي، لم يكن مقدراً حق قدره من قبل، ويؤثر على التنوع الجيني.
• مفارقة ليوونتن: بينما يسبب الانتخاب المتقلب انخفاضات كبيرة في $N_e$، تخلص الدراسة إلى أنه غير كافٍ لحل مفارقة ليوونتن بالكامل بمفرده لأن النقص مستقل إلى حد كبير عن $N_c$. بدلاً من ذلك، فإنه يعمل بالتزامن مع الانتخاب المرتبط والتاريخ الديموغرافي.
• الديناميكيات التطورية: تسلط النتائج الضوء على أن الانتخاب المتقلب يخلق "سحباً" مستمراً على التنوع الجيني، وهو أمر متميز عن الآثار العابرة للمسحات الانتخابية. وهذا له تداعيات على علم الوراثة الحفظي وعلى التنبؤ بالقدرة التكيفية في ظل التغيرات المناخية السريعة.
• تقدم منهجي: من خلال النجاح في إعادة إنتاج تذبذبات التردد الأليلي التجريبية وربطها بـ $N_e$ عبر المحاكاة الأمامية، توفر هذه الدراسة إطاراً قوياً للتحقيقات المستقبلية حول الانتخاب الزمني في مجموعات ضريبية متنوعة.