Disentangling mitochondrial copy number variation and PCR amplification bias in DNA metabarcoding

إليك ملخص تقني مفصل للورقة البحثية بعنوان: "فك الارتباط بين تباين عدد نسخ الحمض النووي الميتوكوندري وانحياز التضخيم عبر تفاعل البوليميراز المتسلسل (PCR) في استقصاء الميتاباركودينج للحمض النووي (DNA metabarcoding)."

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

يعد استقصاء الميتاباركودج للحمض النووي أداة قوية لمراقبة التنوع البيولوجي، ولكنه يواجه قيداً حرجاً: وهو عدم القدرة على تقدير الكتلة الحيوية أو الوفرة بشكل موثوق من خلال عدد القراءات التسلسلية. حدد المؤلفون مصدرين رئيسيين للانحياز اللذين يفصلان بين عدد القراءات والواقع البيولوجي:

1. تباين عدد نسخ الحمض النووي الميتوكوندري (mtDNA): يختلف عدد نسخ الجينات العلامة الميتوكوندرية (مثل COI) لكل وحدة كتلة حيوية بشكل كبير بين الأنواع وداخلها، وكذلك بين الأنسجة والمراحل الحياتية. يمكن أن يصل هذا التباين إلى 200 ضعف، مما يعني أن أعداد التسلسلات لا تعكس الكتلة الحيوية للكائن الحي بشكل خطي.
2. انحياز التضخيم عبر الـ PCR: تؤدي عدمات التطابق في البادئات (Primer mismatches) والاختلافات النوعية في التسلسل إلى تباين في كفاءة التضخيم ($E$) أثناء عملية الـ PCR. يتسبب هذا في تمثيل مفرط أو ناقص لبعض الأصناف بشكل غير متناسب في بيانات التسلسل النهائية.

بينما اقتُرحت استراتيجيات مثل معايرة دورات الـ PCR (تغيير عدد الدورات لتقدير تركيز القالب الأولي) لتصحيح انحياز التضخيم، إلا أن فعاليتها لا تزال مثيرة للجدل. علاوة على على ذلك، فإن التفاعل بين تباين عدد نسخ الـ mtDNA وانحياز الـ PCR في المجتمعات المعقدة متعددة الأنواع ليس مفهوماً بشكل جيد.

2. المنهجية

استخدمت الدراسة تصميماً تجريبياً صارماً باستخدام مجتمعات وهمية (mock communities) لعزل وقياس هذه المتغيرات:

• تجميع المجتمعات الوهمية:

  • تم إنشاء 81 مجتمعاً وهمياً باستخدام خمسة أصناف من المفصليات تغطي رتباً مختلفة: Acromyrmex (Hymenoptera)، و Blaptica dubia (Blattodea)، و Galleria mellonella (Lepidoptera)، و Gammarus (Amphipoda)، و Hermetia illucens (Diptera).
  • غطت العينات تدرجات منهجية لنسب الكتلة الحيوية (10%–90%) ومجموعات متنوعة من الأنواع (من 2 إلى 5 أنواع).
  • تم تجميع إجمالي 84 مجتمعاً؛ واستُبعد 3 منها بسبب أخطاء في التعامل.

• تقنيات القياس الكمي:

  • تفاعل البوليميراز المتسلسل الرقمي القطيري (ddPCR): استُخدم لتحديد عدد نسخ الـ mtDNA المطلقة لكل نوع في كل عينة، ليكون بمثابة "الحقيقة الأرضية" لوفرة القالب.
  • الميتاباركودينج (Metabarcoding): نُفذ باستخدام بروتوكول PCR ذي خطوتين مع مجموعتين مختلفتين من البادئات تستهدف جين COI:
    • Fwh2/Fwh2Rn: مجموعة بادئات متحللة (degenerate) معروفة بتغطيتها الواسعة.
    • BF3/BR2: مجموعة بادئات ثانية للتحقق من الصحة.
  • معايرة دورات الـ PCR: بالنسبة لبادئ Fwh2، تم تشغيل خطوة الـ PCR الأولى عند 8 أعداد مختلفة من الدورات (6، 8، 10، 12، 14، 16، 18، و20 دورة) لاختبار ما إذا كانت نسب القراءات ستتغير بتغير عدد الدورات.

• تحليل البيانات:

  • شملت المعالجة المعلوماتية الحيوية فك التجميع (demultiplexing)، وإزالة الضجيج (denoising)، والتصنيف (الذي تم تجميعه على مستوى الرتبة).
  • تم حساب درجة التمثيل (Representation Score) لقياس الانحراف بين القراءات النسبية ونسخ الـ mtDNA النسبية.
  • تم تطوير إطار رياضي لنمذجة انحياز التضخيم غير الأسي واشتقاق كفاءات التضخيم النسبية الخاصة بكل صنف.

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

تقدم الورقة ثلاث مساهمات نظرية ومنهجية رئيسية:

1. رفض معايرة دورات الـ PCR: أظهرت الدراسة أن تغيير عدد دورات الـ PCR لا يغير النسب النسبية لقراءات الأصناف في المجتمع. وهذا يناقض الافتراض بأن الانحياز يتراكم بشكل أسي مع الدورات. يجادل المؤلفون بأن عدمات تطابق البادئات يتم حلها بعد أول دورتين (عندما يتم استبدال القالب الأصلي بمنتجات تضخيم ذات مواقع ارتباط مثالية للبادئات)، مما يجعل معدل التضخيم ثابتاً بعد ذلك.
2. الاشتقاق الرياضي للانحياز غير الأسي: يقترح المؤلفون نموذجاً ميكانيكياً حيث لا يكون التضخيم أسياً بحتًا فيما يتعلق بالانحياز. لقد اشتقوا معادلات لحساب كفاءة التضخيم النسبية ($E_{rel}$) بناءً على صنف مرجعي، مما يسمح بتقدير أعداد نسخ القالب الأولية من بيانات القراءات دون الحاجة إلى تغيير عدد الدورات.
3. القياس الكمي لتباين الـ mtDNA: تقدم الدراسة دليلاً تجريبياً على التباين الهائل في أعداد نسخ الـ mtDNA لكل وحدة كتلة حيوية عبر مختلف رتب المفصليات، مما يسلط الضوء على هذا كعائق أساسي أمام قياس الكتلة الحيوية.

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

• الـ mtDNA مقابل الكتلة الحيوية: بينما تدرجت أعداد نسخ الـ mtDNA طردياً مع الكتلة الحيوية المدخلة ($\rho=0.65$)، كان هناك تشتت كبير. تباينت أعداد النسخ بشكل ملحوظ داخل وبين الأنواع، مما يؤكد أنه لا يمكن استنتاج الكتلة الحيوية مباشرة من أعداد الـ mtDNA دون عوامل تصحيح خاصة بكل نوع.
• القراءات مقابل أعداد نسخ الـ mtDNA: أظهرت قراءات الميتاباركودينج غير المصححة ارتباطاً معنوياً ولكنه مشتت للغاية مع نسخ الـ mtDNA ($\rho=0.75$).
  • تم تمثيل Blattodea و Diptera و Lepidoptera بشكل جيد نسبياً.
  • تم تمثيل Amphipoda و Hymenoptera بشكل ناقص بشدة (وصل نقص التمثيل في Amphipoda إلى 2000 ضعف).
• فشل معايرة دورات الـ PCR: لم تؤدِ زيادة دورات الـ PCR من 6 إلى 20 إلى تغيير النسب النسبية للقراءات. وهذا يؤكد أن الانحياز يتحدد في وقت مبكر ويظل مستقراً، مما يجعل معايرة الدورات استراتيجية تصحيح غير فعالة.
• نجاح التصحيح:
  • باستخدام النموذج الرياضي المشتق، حسب المؤلفون كفاءات التضخيم النسبية لكل صنف بالنسبة لمرجع (Blattica dubia).
  • أدى تطبيق تصحيحات الكفاءة هذه إلى تحسين دقة أعداد نسخ الـ mtDNA المتوقعة من بيانات القراءات بشكل كبير.
  • بالنسبة لبادئ Fwh2، كان التصحيح فعالاً للغاية (Spearman $\rho=0.95$). أما بالنسبة لـ BF3، فقد كان التحسن كبيراً أيضاً ولكنه أقل وضوحاً ($\rho=0.96$)، لأن BF3 يمتلك بالفعل تمثيلاً ذاتياً أفضل.

5. الأهمية والقيود

الأهمية:

• الرؤية الميكانيكية: توضح الدراسة أن انحياز الـ PCR في الميتاباركودينج هو على الأرجح غير أسي ومدفوع بتفاعلات البادئ والقالب في الدورات المبكرة، مما يتحدى الافتراضات السابقة بأن الانحياز يتراكم عبر الدورات.
• إطار التصحيح: توفر نهجاً رياضياً قابلاً للتطبيق لتصحيح انحياز التضخيم باستخدام عوامل الكفاءة النسبية، بشرما توفر صنف مرجعي.
• التوقعات الواقعية: تؤكد الدراسة أنه حتى مع تصحيح الانحياز، فإن تقدير الكتلة الحيوية أو الوفرة الفعلية يظل غير دقيق للغاية بسبب التباين البيولوجي المتأصل وغير القابل للتصحيح في عدد نسخ الـ mtDNA لكل خلية/نسيج.

القيود والتوجهات المستقبلية:

• القابلية للتطبيق العملي: تتطلب طريقة التصحيح معرفة كفاءة التضخيم النسبية لكل نوع في العينة بالنسبة لمرجع ما، وهو أمر مستحيل حالياً في العينات البيئية عالية التنوع (مثل مصائد Malaise التي تحتوي على أنواع غير معروفة).
• تقدير الكتلة الحيوية: حتى مع التصحيح المثالي للانحياز، فإن التباين العالي في أعداد نسخ الـ mtDNA لكل وحدة كتلة حيوية يمنع التقدير الدقيق للكتلة الحيوية.
• الاحتياجات المستقبلية: يقترح المؤلفون أن الميتاباركودينج الكمي المستقبلي قد يتطلب:
  • استخدام علامات نووية ذات نسخة واحدة (أقل تبايناً في عدد النسخ).
  • تطوير معايير "spike-in" عالمية.
  • التطبيق على العينات منخفضة التنوع (مثل مراقبة اللافقاريات الكبيرة) حيث تكون قوائم الأنواع معروفة.

في الختام، تقدم الدراسة رؤى مفاهيمية أساسية حول آليات انحياز الـ PCR وتباين الـ mtDNA، وتوفر أداة رياضية للتصحيح، بينما تسلط الضوء في الوقت نفسه على القيود البيولوجية الأساسية التي تحد حالياً من التطبيق الكمي للميتاباركودينج للحمض النووي.