Construction of a High-Density SNP Genetic Map and QTL Mapping of Cold Tolerance with Candidate Gene Screening in Peach
أنشأت هذه الدراسة خريطة جينية عالية الكثافة لتعدد أشكال النوكليوتيدات المفردة (SNP) وحددت 12 موقعاً للسمات الكمية (QTLs) مرتبطة بتحمل البرودة في الخوخ، مما أدى إلى تحديد الأولوية لتسعة جينات مرشحة من خلال دمج رسم خرائط الـ QTL والتحليل النسخي لتسهيل التربية الجزيئية.
المؤلفون الأصليون:Qihang Tian, Jie Li, Zhendong Tian, Yu Wang, Zhaoyuan Wang, Guojian Liu, Ruifeng Chang, Fei Wang, Yonghong Li
تُعد أشجار الخوخ محبوبة لثمارها الحلوة، لكنها أيضاً ناجات رقيقات في ظل مناخ متغير. هذه الأشجار، الموطن الأصلي لها هو الصين والتي تُزرع الآن في جميع المناطق المعتدلة، تواجه معركة مستمرة ضد البرد. فعندما تنخفض درجات حرارة الشتاء إلى مستويات منخفضة جداً، يمكن للمياه داخل خلايا الشجرة أن تتجمد، فتتمدد وتؤدي إلى تمزق الجدران الخلوية، تماماً مثل انفجار بالون مائي في مجمد (فريزر). يؤدي هذا الضرر إلى موت الأغصان، وفي الحالات الشديدة، موت الشجرة بأكملة. وللبقاء على قيد الحياة، طورت أشجار الخوخ دفاعات داخلية معقدة: فهي تضخ مركبات طبيعية مضادة للتجمد مثل السكريات والأحماض الأمينية، وتقوي أغشيتها الخلوية لإبقاء البرد في الخارج. ومع ذلك، فإن القدرة على تحمل ظروف التجمد هذه تختلف بشكل هائل من شجرة إلى أخرى، وهي سمة تتحكم فيها شبكة خفية من الجينات التي ظل العلماء يكافحون طويْلًا لرسم خريطتها. ويعد فهم هذه الشفرة الوراثية هو المفتاح لاستنباط أصناف جديدة يمكنها الازدهار في المناطق الأكثر برودة دون أن تموت في فصل الشتاء.
في دراسة حديثة، شرع الباحثون في فك هذا اللغز الجيني عبر بناء خريطة تفصيلية لجينوم الخوخ، وبالتحديد البحث عن التعليمات التي تحكم البقاء في البرد. بدأوا بتهجين بين نوعين مختلفين تماماً من الآباء: 'دونغهي رقم 1' (Donghe No. 1)، وهو صنف قوي معروف بقدرته على تحمل درجات حرية تصل إلى 30 درجة مئوية تحت الصفر، و'القرن الحادي والعشرين' (21st Century)، وهو صنف شائع ولكنه هش يعاني في البرد. ومن خلال هذا التهجين بين الاثنين، أنشأ الفريق عائلة مكونة من 150 سليلًا، حيث ورث كل منها مزيجاً فريداً من الجينات من والديه. ولكي يحددوا أي الجينات كانت مسؤولة عن الاختلاف في تحمل البرد، احتاج الباحثون إلى وسيلة لتتبع وراثة كل قطعة صغيرة من الحمض النووي. وقد استخدموا تقنية تسلسل حديثة لمسح الحمض النووي للآباء وجميع السلال الـ 150، مما أدى إلى تحديد أكثر من 2600 علامة جينية محددة. تعمل هذه العلامات كلوحات إرشادية على طول كروموسومات الشجرة، مما يسمح للعلماء برؤية أي أجزاء من الجينوم جاءت من الأب القوي وأيها جاء من الأب الهش.
وكانت النتيجة خريطة جينية عالية الدقة، وهي ترتيب خطي لهذه اللوحات الإرشادية التي تغطي جينوم الخوخ بأكمله. وتتميز هذه الخريطة بدقة متناهية، حيث تفصل بين العلامات مسافات ضيقة جداً بحيث يقل متوسط المسافة بينها عن ثلث سنتيمورغان (centimorgan)، وهي وحدة تستخدم لقياس المسافة الجينية. ثم وضع الباحثون هذه الأشجار تحت الاختبار في الحقل، وراقبوها على مدار ثلاثة فصول شتاء متتالية. وقاموا بقياس خمس علامات مختلفة لصحة الشجرة تجاه البرد في الأغصان التي يبلغ عمرها عاماً واحداً: مقدار التيار الكهربائي الذي يمكن أن يمر عبر الأنسجة التالفة، وكمية السكر والنشا التي خزنتها الأشجار، ومستويات الأحماض الأمينية الواقية، ودرجة الحرارة التي بدأت عندها أنسجة النبات في التجمد. ومن خلال مقارنة هذه القياسات الفيزيائية بالخرائط الجينية، استطاع الفريق تحديد أي أقسام من الجينوم مرتبطة بتحسن البقاء على قيد الحياة بدقة.
وكشف التحليل عن اثني عشر منطقة متميزة على الكروموسومات تؤثر بشكل كبير على تحمل البرد. وُجدت معظم هذه المناطق الهامة في كروموسومين محددين، حيث برزت منطقة واحدة على الكروموسوم الثامن كمركز رئيسي لمقاومة البرد. وقد فسرت هذه المنطقة المحددة ما يقرب من 15 بالمائة من الاختلافات المشاهدة في مدى قدرة الأشجار على البقاء في البرد، وهي مساهمة كبيرة بالنسبة لسمة معقدة كهذه. كما وجد الباحثون أن بعض مناطق الجينات المفيدة هذه ظهرت باستمرار عبر جميع سنوات الدراسة الثلاث، مما يشير إلى أنها أهداف مستقرة وموثوقة للتربية. ولتضييق نطاق البحث أكثر، نظر الفريق في نشاط الجينات في شجرتي الأبوين عندما تعرضتا لبرد شديد في بيئة مختبرية. وقارنوا الجينات التي تم تفعيلها أو إيقافها أثناء الإجهاد بمواقع مناطق تحمل البرد في خريطتهم الجينية.
سمحت عملية التدقيق المزدوج هذه للعلماء بتصفية آلاف الجينات غير ذات الصلة والتركيز على قائمة قصيرة تضم تسعة مرشحين رئيسيين. وتشارك هذه الجينات في مهام بقاء حاسمة، مثل إرسال الإشارات عند انخفاض درجة الحرارة، وإعادة بناء الغشاء الخلوي للحفاظ على مرونته، وإنتاج البروتينات التي تحمي الشجرة من أضرار التجمد. فعلى سبيل المثال، يساعد أحد الجينات المحددة في إعادة تشكيل الدهون في الغشاء الخلوي، بينما يعمل آخر كمستشعر يحفز نظام الدفاع لدى الشجرة. وتؤكد الدراسة أن تحمل البرد في الخوخ لا يتحكم فيه جين "سحري" واحد، بل بمجموعة من العوامل الجينية الصغيرة التي تعمل معاً. ورغم أن الباحثين لم يثبتوا بعد بدقة كيفية عمل كل جين من هذه الجينات التسعة في شجرة حية، إلا أن اكتشافهم يوفر أساساً متيناً للمستقبل. فمن خلال تحديد هذه المناطق الجينية المحددة والجينات الموجودة بداخلها، أصبح أمام المربين الآن مسار واضح للمضي قدماً لتطوير أصناف جديدة من الخوخ يمكنها تحمل فصول الشتاء الأكثر قسوة، مما يضمن بقاء المحصول في عالم يزداد دفئاً وتقلبات.
ملخص تقني: بناء خريطة وراثية عالية الكثافة لتعدد أشكال النوكليوتيدات المفردة (SNP) ورسم خرائط الصفات الكمية (QTL) لتحمل البرودة مع فحص الجينات المرشحة في الخوخ
بيان المشكلة يعد تحمل البرودة سمة كمية بالغة الأهمية ومعقدة في زراعة الخوخ (Prunus persica)، حيث تؤثر بشكل مباشر على البقاء خلال فصل الشتاء، والتكيف الإقليمي، واستقرار المحصول. تسبب درجات الحرارة المنخفضة تلفاً في الأغشية، وتسرب الأيونات، وتراكم أنواع الأكسجين التفاعلية (ROS)، مما يؤدي غالباً إلى موت النبات. وبينما ركزت الأبحاث السابقة على التقييمات الفسيولوجية وفحص الأصول الوراثية، لا تزال الآليات الوراثية الجزيئية الكامنة وراء تحمل البرودة في الخوخ غير مفهومة بشكل جيد. كما أن الخرائط الوراثية التقليدية التي تستخدم علامات مثل RAPD وAFLP وSSR تفتقر غالباً إلى الكثافة والدقة المطلوبة لرسم خرائط الصفات الكمية (QTL) بدقة. علاوة على ذلك، فإن التحديد الدقيق لخرائط QTL الخاصة بتحمل البرودة في الخوخ نادر بسبب الطبيعة متعددة الجينات والارتباط القوي بين النمط الجيني والبيئة.
المنهجية استخدمت هذه الدراسة نهجاً متكاملاً يجمع بين التسلسل عالي الإنتاجية، والخرائط الوراثية عالية الكثافة، والنمذجة الظاهرية لعدة سنوات، والتحليل الترانسكريبتومي:
المواد النباتية: تم إنتاج مجتمع F1 مكون من 150 شتلة عن طريق تهجين الصنف المتحمل للبرودة 'Donghe No. 1' (أنثى) مع الصنف الحساس للبرودة '21st Century' (ذكر).
النمذجة الظاهرية: تم تقييم تحمل البرودة على مدار ثلاث سنوات متتالية (2022–2024). تعرضت الفروع الساكنة ذات العمر سنة واحدة لإجهاد البرودة عند درجة حرارة -25 مئوية لمدة 12 ساعة. وتم قياس خمسة مؤشرات فسيولوجية وكيميائية حيوية: الموصلية الكهربائية النسبية (REC)، والسكريات الذائبة (SS)، والبرولين (Pro)، والنشا (St)، والتحليل الحراري التفاضلي (DTA).
التنميط الجيني وبناء الخريطة: استُخدم تسلسل تضخيم الموقع النوعي (SLAF-seq). تم هضم الحمض النووي الريبوزي منقوص الأكسجين (DNA) للآباء والأبناء مزدوجاً باستخدام الإنزيمين RsaI وHaeIII. وبعد استخدام تقنية Illumina PE150 وتصفية صارمة (العمق، السلامة، تشوه الانعزال)، تم بناء خريطة وراثية عالية الكثافة باستخدام طريقة الاختبار الكاذب المزدوج (double pseudo-testcross).
رسم خرائط QTL: تم إجراء رسم خرائط الفترات المركبة باستخدام برنامج MapQTL 6.0 على البيانات الظاهرية من 2022–2024. وتم تحديد عتبات المعنوية عبر 1,000 اختبار تبديل (permutation tests).
فحص الجينات المرشحة:
الجينومي: تمت محاذاة فترات QTL مع المرجع الجينومي للخوخ (v2.0.a1) لاستخراج تعليقات الجينات.
الترانسكريبتومي: أُجري تسلسل RNA-seq على الخطوط الأبوية تحت درجات الحرارة العادية وتحت إجهاد -25 درجة مئوية. وتم تحديد الجينات ذات التعبير التفاضلي (DEGs) بمعيار (|Log2FC| ≥ 1، FDR < 0.05).
التكامل: تم إعطاء الأولوية للجينات الموجودة ضمن فترات QTL والتي كانت أيضاً ذات تعبير تفاضلي (سواء بشكل مستمر أو بشكل محدد تحت إجهاد البرودة).
التحقق: تم التحقق من الجينات المرشحة المختارة عبر تقنية qRT-PCR.
النتائج الرئيسية
الخريطة الوراثية عالية الكثافة: أنتج نهج SLAF-seq ما مجموعه 2,601 علامة SNP عالية الجودة تم رسم خرائطها على ثماني مجموعات ربط (LGs). بلغ الطول الإجمالي للخريطة 761.55 سنتيمورغان (cM)، بمتوسط فاصل بين العلامات قدره 0.29 cM، وبسلامة علامات بلغت 99.76%. أظهرت الخريطة توافقاً عالياً مع الجينوم الفيزيائي (معاملات ارتباط سبيرمان تراوحت بين 0.86 و1.00).
تحديد QTL: تم تحديد اثني عشر QTL معنوياً مرتبطة بتحمل البرودة عبر ست مجموعات ربط (Pp03, Pp04, Pp05, Pp06, Pp07, Pp08).
تراوح التباين الظاهري المفسر (PVE) بواسطة QTLs الفردية من 6.68% إلى 16.38%.
تم تحديد عنقود QTL رئيسي على Pp08، يحتوي على مواقع للسكريات الذائبة (qSS08-3)، والموصلية الكهربائية النسبية (qREC08-1, qREC08-2)، وسمات أخرى. فسر الموقع qREC08-1 نسبة 14.80% من التباين الظاهري.
كان الـ QTL ذو التأثير الأقوى هو qDTA06-1 على Pp06، حيث فسر 16.38% من التباين.
لوحظ وجود تموضع مشترك لـ QTLs الخاصة بالبرولين والنشا على Pp03.
فحص الجينات المرشحة:
أسفر التوصيف الأولي لفترات QTL عن 1,290 جيناً فريداً.
أدى التكامل مع البيانات الترانسكريبتومية (DEGs تحت ظروف الحرارة العادية وإجهاد البرودة) إلى تضييق النطاق إلى 61 جيناً من الـ DEGs الواقعة ضمن نطاق QTL.
حددت عملية التصفية الوظيفية الإضافية تسعة جينات مرشحة أساسية. وتشمل هذه الجينات:
عوامل النسخ: DREB1F، وبروتين DREB1F، وبروتين NAC (NAC domain-containing protein 78).
الإشارات والتمثيل الغذائي: GDSL esterase/lipase (المشارك في إعادة تشكيل الدهون)، وLRR receptor-like kinase PSKR1 (إشارات تعتمد على الكالسيوم)، ومكونات مسار MAPK.
أكد التحقق باستخدام qRT-PCR اتجاهات التعبير للجينات المختارة، مما يدعم موثوقية البيانات الترانسكريبوتومية.
الأهمية والادعاءات يزعم المؤلفون أن هذه الدراسة توفر مورداً وراثياً قيماً للتحقيق في الآليات الجزيئية لتحمل البرودة في الخوخ. ومن خلال بناء خريطة SNP عالية الكثافة ذات دقة فائقة مقارنة بالخرائط التقليدية، تتغلب الدراسة على القيود السابقة في اكتشاف QTLs تحمل البرودة. إن تحديد عنقود QTL رئيسي على Pp08 وما تلاه من فحص لتسعة جينات مرشحة يقدم أهدافاً محددة للتربية الجزيئية.
تضع الورقة هذا العمل كخطوة نحو كسر قيود الزراعة الإقليمية لأصناف الخوخ عالية الجودة. وتؤكد أن دمج رسم خرائط QTL مع التوصيف الترانسكريبتومي يوفر نهجاً فعالاً لفحص الجينات المرتبطة بالإجهاد غير الحيوي في الأشجار المثمرة المعمرة، مما يتجاوز الحاجة إلى عمليات فحص واسعة النطاق ومجهدة للجهد البشري من أجل التحديد الدقيق (fine-mapping). وتعمل الجينات المرشحة و QTLs المحددة كقاعدة نظرية لتطوير علامات جزيئية وظيفية وتعزيز الانتخاب بمساعدة العلامات (marker-assisted selection) للأصول الوراثية من الخوخ المقاوم للبرد.
القيود الملحوظة يقر المؤلفون بأن وظائف الجينات المرشحة تم استنتاجها بشكل أساسي من بيانات التوصيف والتعبير دون تحقق وظيفي في الجسم الحي أو في المختبر (مثل الإفراط في التعبير أو التعطيل الجيني). وبناءً على ذلك، لا تزال الآليات التنظيمية الدقيقة بحاجة إلى توضيح في العمل المستقبلي. بالإضافة إلى ذلك، تشير الدراسة إلى أن اكتشاف QTL قد تباين عبر السنوات بسبب التأثيرات البيئية، وهو أمر نموذجي للصفات الكمية المعقدة في الأشجار المثمرة المزروعة في الحقول.