Genetic Dissection of Spot Blotch Resistance in Wheat: Novel Loci, Candidate Genes, and Genomic Prediction
توضح هذه الدراسة البنية الجينية لمقاومة تبقع الأوراق في القمح من خلال تحديد 185 ارتباطاً بين العلامات والسمات، واكتشاف موقع جيني ثلاثي التكرار (3D locus) جديد، وتوصيف الجينات المرشحة ضمن مناطق سمات الكمية (QTL) الغنية، وإثبات فعالية التنبؤ الجينومي في تسريع استنباط الأصناف المقاومة.
يعد القمح المحصول الغذائي الأكثر حيوية في العالم، وهو غذاء أساسي يغذي المليارات. ومع ذلك، فإن نموه مهدد باستمرار بمرض فطري يسمى "لفحة البقع" (spot blotch). يزدهر هذا الممرض في المناخات الدافئة والرطبة، لا سيما في أجزاء من جنوب آسيا، وأمريكا الجنوبية، وأفريقيا، حيث يمكنه التهام مساحات شاسعة من المحاصيل. وعندما يضرب المرض بقوة، يمكنه تدمير ما يصل إلى سبعين بالمائة من الحصاد، مما يترك المزارعين بصوامع فارغة والمجتمعات تواجه الجوع. ولعقود من الزمن، كان الدفاع الأساسي ضد هذا العدو هو الرش الكيميائي، لكن هذه المواد مكلفة وغير متوفرة غالبًا لصغار المزارعين الذين هم في أمس الحاجة إليها. ويكمن الحل الأكثر استدامة داخل القمح نفسه: وهو استنباط أصناف تقاوم الفطر بشكل طبيعي. ومع ذلك، فإن المقاومة في القمح ليست مجرد مفتاح بسيط؛ بل هي سمة معقدة تتحكم فيها أجزاء مختلفة عديدة من الجينوم تعمل معًا، مما يجعل من الصعب تحديد أي الجينات توفر الحماية بدقة.
ولحل هذا اللغز، وضع فريق من الباحثين من الصين والمكسيك خطة لرسم الخريطة الجينية لمقاومة لفحة البقع. قاموا بجمع مجموعة متنوعة تضم 255 صنفًا من القمح، تمثل أصولاً وتاريخًا مختلفًا من التربية، وزرعوها في تجارب ميدانية في المكسيك على مدار ثلاث سنوات متتالية. تعمد الفريق إصابة المحاصيل بالفطر ليروا كيف ستصمد كل سلالة. وقاموا بقياس الضرر بعناية، مع تعديل بياناتهم لتفسير كيفية نضج كل نبات مبكرًا أو متأخرًا، نظرًا لأن التوقيت يمكن أن يؤثر على كمية المرض التي يتعرض لها النبات. كانت النتائج مذهلة: كان التركيب الجيني للقمح هو العامل المهيمن في تحديد البقاء، بينما لعبت البيئة دورًا ثانويًا. والأهم من ذلك، كشفت الدراسة عن نمط واضح من التفوق؛ حيث أظهرت أصناف القمح التي تنحدر من الصين باستمرار مستويات أقل بكثير من أضرار المرض مقارنة بالأصناف القادمة من أجزاء أخرى من العالم، بما في ذلك تلك القادمة من المركز الدولي لتحسين الذرة والقمح في المكسيك.
وبالتعمق أكثر في الحمض النووي لهذه النباتات، استخدم الباحثون طريقة إحصائية قوية لمسح جينوم القمح بأكمله بحثًا عن علامات محددة مرتبطة بالمقاومة. وقد حددوا 185 إشارة جينية متميزة مرتبطة بالقدرة على محاربة الفطر. كانت هذه الإشارات منتشرة عبر كروموسومات القمح، لكنها لم تكن موزعة بالتساوي. برز الجينوم الفرعي B، وهو أحد المجموعات الجينية الثلاث الرئيسية التي يتكون منها القمح، كبؤرة ساخنة للمقاومة، حيث احتوى على أكبر عدد من الإشارات الواقية. وكان الاكتشاف الهام بشكل خاص هو منطقة مقاومة جديدة وقوية وجدت في الكروموسوم 3D. لم تكن هذه المنطقة قد سُلط الضوء عليها كمنطقة دفاع رئيسية في الدراسات السابقة، مما يمثل توسعًا حقيقيًا في الأدوات الجينية المعروفة المتاحة للمربين.
ثم نظر الباحثون عن كثب في النسخ المحددة من الجينات، أو الأنماط الفردانية (haplotypes)، الموجودة في هذه المناطق المقاومة. ووجدوا أن الأصناف الصينية تحمل مزيجًا فريدًا من هذه الكتل الجينية الواقية التي كانت نادرة أو غائبة في الأصناف غير الصينية. على سبيل المثال، وُجد نمط واقٍ محدد في الكروموسوم 3D في أربعة وثلاثين خطًا صينيًا ولكن في أربعة خطوط أخرى فقط. لم يكن هذا "التفوق الصيني" ناتجًا عن جين خارق واحد، بل عن تراكم العديد من الاختلافات الجينية الصغيرة والمواتية التي تعمل بتضافر. كما حددت الدراسة 575 جينًا تقع ضمن هذه المناطق المقاومة. ومن خلال تحليل وظائف هذه الجينات، أعاد الفريق بناء نظام دفاع متعدد الطبقات. ويبدو أن نباتات القمح تقاوم من خلال تنظيم كيفية تفكيك البروتينات، وتعزيز جدران خلاياها بحواجز كيميائية، وإنتاج مواد تثبط قدرة الفطر على الهجوم.
أخيرًا، اختبر الفريق ما إذا كان بإمكانهم استخدام هذه المعلومات الجينية للتنبؤ بأداء أصناف القمح المستقبلية قبل زراعتها في الحقل. قاموا ببناء نموذج حاسوبي يستخدم البيانات الجينية للتنبؤ بمدى مقاومة النبات لمرض لفحة البقع. أثبت النموذج دقة عالية، حيث تنبأ بالنتائج بدرجة عالية من الموثوقية. ومن الأهمية بمكان أن النموذج عمل بشكل أفضل عندما استخدم المجموعة الكاملة من العلامات الجينية المرتبطة بالمرض، بدلاً من مجرد العلامات القليلة التي بدت مستقرة عبر جميع السنوات. وهذا يشير إلى أن حتى التأثيرات الجينية الصغيرة والدقيقة تساهم في القوة الإجمالية لدفاع النبات. توفر هذه النتائج خارطة طريق مفصلة للمربين؛ فمن خلال التركيز على المناطق الجينية المحددة، وخاصة تلك المتوفرة بكثرة في السلالات الصينية، واستخدام التنبؤ الجينومي لتسريع عملية الاختيار، يمكن للعلماء الآن تطوير أصناف جديدة من القمح مجهزة طبيعيًا لمواجهة لفحة البقع، مما يوفر مستقبلًا أكثر أمانًا واستدامة لإنتاج الغذاء العالمي.
بيان المشكلة يعد مرض تبقع البقع (SB) في القمح، والذي يسببه الفطر نصف الحيوي (hemibiotrophic) Bipolaris sorokiniana، عائقاً حرجاً أمام إنتاج القمح العالمي، لا سيما في المناطق الدافئة والرطبة في جنوب آسيا، وأمريكا الجنوبية، وأفريقيا. يتسبب هذا المرض في خسائر كبيرة في المحاصيل (بمتوسط 15-20%، وتصل إلى 70% في حالات الأوبئة)، ويتم مكافحته أساساً باستخدام مبيدات الفطريات، وهي مكلفة وغير مستدامة بيئياً. وبينما تعتبر مقاومة العائل هي الاستراتيجية المفضلة للإدارة، فإن مقاومة تبقع البقع هي صفة كمية تتأثر بالظروف البيئية. وعلى الرغم من أن الدراسات السابقة قد حددت عدة مواقع سمات كمية (QTLs) للمقاومة، مثل Sb1 (المتمركز مع Lr34) و tsn1، إلا أن غالبية المواقع المحددة لها تأثيرات صغيرة، ولا تزال العلامات التشخيصية للاختيار بمساعدة العلامات (MAS) محدودة. علاوة على ذلك، فإن الأساس الجيني للمقاومة المتفوقة الملحوظة في بعض الأصناف، وخاصة الأصناف الصينية، يتطلب مزيداً من التوضيح.
المنهجية استخدمت الدراسة نهجاً شاملاً يدمج التقييم المظهري، والنمط الجيني عالي الكثافة، والنمذجة الإحصائية المتقدمة:
المواد النباتية: تم تقييم مجموعة متنوعة تضم 255 من أصناف القمح، تتكون من 138 صنفاً صينياً و117 صنفاً من CIMMYT، وأمريكا الجنوبية، وأمريكا الشمالية، وأوروبا، وآسيا، وأوقيانوسيا، وأفريقيا.
التقييم المظهري: أُجريت تجارب ميدانية في محطة "Agua Fria" التابعة لـ CIMMYT في المكسيك على ثلاثة مواسم نمو متتالية (2020-21 إلى 2022-23). تم تلقيح الأصناف اصطناعياً بمزيج من عزلات B. sorokiniana شديدة الإمراض. تم تسجيل شدة المرض باستخدام مقياس "Saari and Prescott" ذي الرقمين، وحُسبت المساحة تحت منحنى تقدم المرض (AUDPC). وللتحكم في الارتباط السلبي المعروف بين تاريخ الإزهار وشدة المرض، تم تعديل بيانات الشدة الخام (SBadj) باستخدام نموذج انحدار.
النمط الجيني: تم إجراء النمط الجيني باستخدام تقنية DArTSeq، مما أسفر عن 18,436 علامة عالية الجودة. وبعد تصفية القيم المفقودة وتكرار الأليل الأدنى، استُخدمت 14,195 من تعدد أشكال النوكليوتيدات المفردة (SNPs) ذات المواقع الفيزيائية المعروفة للتحليل.
التحليل الإحصائي:
دراسة الارتباط الجينومي الكامل (GWAS): أُجريت دراسات GWAS متعددة النماذج باستخدام نماذج BLINK وFarmCPU وMLM لتحديد الارتباطات بين العلامات والسمات (MTAs).
تحليل الهابلوتايب (النمط الفرداني): أُجري تحليل الهابلوتايب القائم على عدم التوازن الارتباطي (LD) لتحديد مناطق الـ QTL وتقييم ترددات الأليلات عبر أصول الأنماط الجينية المختلفة.
تنقيب الجينات المرشحة: تم تحديد الجينات داخل مناطق QTL الهامة باستخدام المرجع IWGSC RefSeq v1.1 وقواعد البيانات الوظيفية (مثل EggNOG وUniProt وغيرها).
التنبؤ الجينومي (GP): استُخدم انحدار ريج (Ridge regression) البايزي للتنبؤ بمقاومة تبقع البقع. تمت مقارنة مصفوفتين من الـ SNPs: المصفوفة 1 (جميع الـ SNPs المرتبطة بالسمة) والمصفوفة 2 (الـ SNPs الهامة في بيئتين على الأقل).
النتائج الرئيسية
التباين المظهري والوراثية: أظهرت شدة تبقع البقع وراثة واسعة النطاق عالية (H2=0.92)، مما يشير إلى تحكم جيني قوي. أظهرت الأصناف الصينية باستمرار شدة أقل بكثير في مرض تبقع البقع مقارنة بالأصناف غير الصينية عبر جميع المواسم.
نتائج GWAS: حددت الدراسة 185 من الـ MTAs الهامة الموزعة عبر جميع الكروموسومات باستثناء 6D. احتوى الجينوم الفرعي B على أعلى نسبة من الارتباطات (69 MTA)، وبرز الكروموسوم 2B كمنطقة ساخنة رئيسية (18 MTA).
المواقع الجديدة والمستقرة:
تم تحديد تجمع مقاومة قوي وجديد على الكروموسوم 3D (14 MTA)، وهي منطقة لم يتم تسليط الضوء عليها سابقاً كمنطقة ساخنة رئيسية في دراسات مماثلة.
كشف تحليل الهابلوتايب لتجمع 3D عن أنماط متميزة: كان الهابلوتايب المفضل لـ QTL-SB-9 غنياً جداً في الأصناف الصينية، بينما كان الهلوتايب المفضل لـ QTL-SB-10 أكثر شيوعاً في الأصناف غير الصينية.
كما أظهرت المواقع المستقرة مثل QTL-SB-1 (2A) و QTL-SB-19 (5D) إثراءً كبيراً للهابلوتايبات المقاومة في الأصناف الصينية.
الجينات المرشحة: ضمن 26 منطقة QTL المحددة، تم تحديد 575 جيناً. سلط تحليل الإثراء الوظيفي الضوء على بنية دفاعية متعددة الطبقات تشمل:
نظام اليوبيكويتين-بروتيازوم (UPS): إثراء إنزيمات E3 ligases وإنزيمات ربط اليوبيكويتين، مما يشير إلى تنظيم الاستجابات المناعية عبر دوران البروتين.
الدفاع الكيميائي والهيكلي: مسارات التخليق الحيوي للألكالويدات الأيزوكوينولينية، والبنزوكسازينويدات، والكيوتين/الشمع، إلى جانب وجود مثبطات البروتياز وبروتينات نقل الدهون (LTPs) التي تعمل في الأبوبلاست (apoplast).
عوامل النسخ: تحديد عوامل الاستجابة للإجهاد (WRKY, MYB) وجينات إشارات الهرمونات (SA, ABA).
تم تحديد ما مجموعه 82 من الجينات المرشحة عالية الثقة.
التنبؤ الجينومي: حققت نماذج التنبؤ الجينومي (GP) دقة تنبؤية متوسطة إلى عالية (معاملات ارتباط بيرسون تصل إلى 0.833). تفوقت النماذج التي استخدمت المجموعة الكاملة من الـ SNPs المرتبطة بالسمة (المصفوفة 1) باستمرار على تلك التي استخدمت فقط الـ SNPs المستقرة ومتعددة البيئات (المصفوفة 2)، مما يدل على استقرار تنبؤي فائق.
الأهمية والادعاءات يزعم المؤلفون أن هذه الدراسة توفر تشريحاً جينياً شاملاً لمقاومة تبقع البقع، وتقدم عدة مساهمات رئيسية لتربية القمح:
التحقق من قيمة الأصناف الصينية: يؤكد البحث أن أصناف القمح الصينية تعمل كمستودع حيوي للأليلات الممتازة لمقاومة تبقع البقع، ويرجع ذلك إلى تجميع (pyramiding) الهابلوتايبات المفضلة في مواقع جديدة (مثل 3D) ومستقرة.
توسيع نطاق الترسانة الدفاعية: إن تحديد تجمع مقاومة جديد وقوي على الكروموسوم 3D يوسع المشهد الجيني المعروف لمقاومة تبقع البقع، ويوفر أهدافاً جديدة للتربية.
الرؤى الوظيفية: من خلال ربط إشارات GWAS بجينات مرشحة محددة، توضح الدراسة آلية دفاع "التعرف-التنظيم-التنفيذ" التي تتضمن دوران البروتين بوساطة نظام UPS والدفاعات الكيميائية في الأبوبلاست.
الفائدة في التربية: يوضح إثبات أن استخدام المجموعة الكاملة من الـ SNPs المرتبطة بالسمة يعطي دقة تنبؤ جينومي أعلى، دعم استخدام الاختيار الجينومي (GS) لتسريع تطوير أصناف القمح المقاومة لتبقع البقع. توفر الـ QTLs والجينات المرشحة المحددة موارد قيمة لاستراتيجيات الاختيار بمساعدة العلامات (MAS) وتجميع الجينات (gene pyramiding).
تخلص الورقة إلى أن دمج هذه الموارد الجينية مع التنبؤ الجينومي يمكن أن يعزز بفعالية تربية أصناف القمح ذات المقاومة المستدامة لمرض تبقع البقع.