A comprehensive RNA-seq dataset of embryonic and larval development in Senegalese sole (Solea senegalensis) across three thermal regimes
تقدم هذه الدراسة مجموعة بيانات عالية الجودة لتسلسل الحمض النووي الريبي (RNA-seq) تضم 48 عينة من التطور الجنيني واليرقي لسمك موسى السنغالي عبر ثماني مراحل تحت ثلاثة أنظمة حرارية، مما يوفر مورداً أساسياً لاستقصاء اللدونة الحرارية، والتنظيم التطوري، وتطبيقات الاستزراع المائي في الأسماك المفلطحة.
المؤلفون الأصليون:Lourdes Vega, Daniel Ramirez, Ismael Cross, Oscar Aramburu, Dorinda Torres-Sabino, Maialen Carballeda, María Esther Rodríguez, Marco Mendizabal, Paulino Martínez, Laureana Rebordinos
تخيل المحيط كموقع بناء ضخم وصاخب، حيث يعمل مهندسون صغار غير مرئيين على بناء آلات معقدة تسمى الأسماك. في هذا العالم، ليست درجة الحرارة مجرد رقم على ميزان حرارة؛ بل هي "المراقب" الذي يصرخ بالأوامر، ويخبر المهندسين بمدى سرعة العمل ونوع المواد التي يجب استخدامها. عندما تكون المياه دافئة، قد يهرع الطاقم في عمله؛ وعندما تكون باردة، قد يبطئون سرعتهم لمراجعة مخططاتهم بدقة. يطلق العلماء على هذه القدرة اسم "اللدونة الحرارية" (thermal plasticity)—وهي قدرة الكائن الحي على تغيير سلوكه أو جسده بناءً على الطقس. ولكن هنا يكمن اللغز: بينما يمكننا رؤية الأسماك المكتملة وهي تسبح حولنا، لا يمكننا بسهًا رؤية التعليمات الصغيرة غير المرئية (الجينات) التي تخبر الخلايا بكيفية بناء الجسد في المقام الأول. إن فهم هذه التعليمات أمر بالغ الأهمية لأن صغار الأسماك تشبه دمى الزجاج الهشة؛ فإذا أصبحت الحرارة مرتفعة جدًا أو منخفضة جدًا أثناء عملية البناء، فقد تنتهي بقطع مكسورة أو قد لا تكتمل عملية بنائها على الإطلاق. وهذا أمر مهم جدًا للأشخاص الذين يربون الأسماك من أجل الغذاء، حيث يحتاجون إلى معرفة كيفية الحفاظ على الماء في الظروف المثالية تمامًا لضمان بقاء الصغار ونموهم بقوة.
وهنا، يظهر سمك "سول السينغالي" (Senegalese sole)، وهو سمك مسطح يشبه الفطيرة المسطحة ذات عينين على جانب واحد. قرر فريق من العلماء استراق النظر خلف الستار في المراحل الأولى من حياة هذا السمك لمعرفة كيف تغير درجات الحرارة خطط البناء الخاصة به. لم يكتفوا بمراقبة السمكة فحسب، بل نظروا إلى "كتاب الوصفات" داخل كل خلية، المعروف باسم الحمض النووي الريبوزي (RNA). فكر في الـ RNA كنسخة نشطة من وصفة يقرأها طاقم المطبخ (الخلايا) حاليًا لخبز كعكة ما. إذا تغيرت درجة الحرارة، فقد يختار الطهاة وصفات مختلفة أو يقرؤونها بسرعة أكبر.
أجرى الباحثون تجربة ضخمة شملت 48 مجموعة مختلفة من صغار هذه الأسماك. وقد قاموا بتربيتها في ثلاثة "مطابخ" مختلفة: مطبخ بارد عند 14 درجة مئوية، ومطبخ مريح عند 18 درجة مئوية، ومطبخ دافئ عند 21 درجة مئوية. وتفقدوا الأسماك في ثماني لحظات مختلفة من حياتها، بدءًا من أول انقسام لخلية واحدة (المرحلة البلستولية - blastula) وصولًا إلى المرحلة التي أنهت فيها تحولها الكبير إلى سمكة مسطحة تعيش في القاع (مرحلة ما بعد التحول - post-metamorphic). وفي كل لحظة ومع كل درجة حرارة، أخذوا لقطة للنشاط الجيني للسمكة باستخدام ماسح ضوئي عالي التقنية يسمى (RNA-seq). أنتج هذا عن مكتبة ضخمة تضم حوالي 3-30 مليون "صفحة" من النصوص لكل عينة، مما خلق مجموعة بيانات هائلة تعمل مثل فيلم "التايم لابس" (الفاصل الزمني) لحياة السمكة الجينية.
ماذا وجدوا؟ أظهرت البيانات أن كتاب الوصفات الجينية للسمكة يتغير بشكل كبير مع نموها. حدثت أكبر التغييرات عندما كانت الأسماك تبدأ للتو في بناء أجسامها، وكانت الأنماط مختلفة تمامًا اعتمادًا على ما إذا كانت في مياه باردة، أو دافئة، أو مثالية. لاحظ العلماء أن الأسماك في درجات الحرارة المختلفة لم تكتفِ بالنمو بسرعات متفاوتة فحسب، بل استخدمت بالفعل مجموعات مختلفة من التعليمات. على سبيل المثال، في المراحل المبكرة جدًا، تسببت درجة الحرارة في تشغيل أو إيقاف عدد هائل من الجينات، وكأن طاقم البناء كان يتجادل حول المخططات. ومع تقدم الأسماك في العمر واقترابها من الفقس، بدت تلك الجدالات وكأنها تهدأ، مع تغير أقل في سلوك الجينات بناءً على الحرارة.
وجد الفريق أيضًا بعض "اللاعبين النجوم" في النص الجيني الذين ظهروا مرارًا وتكرارًا، بغض النظر عن المرحلة أو درجة الحرارة. شمل ذلك جينات مثل hsp90aa1.2 و hsc70، والتي تعمل مثل فرق الإصلاح في حالات الطوارئ التي تتدخل عندما تشتد الضغوط (مثل الحرارة العالية). كما رأوا أن مسارات بناء معينة، مثل تلك المخصصة لبناء جلد وعضلات السمكة، كانت متأثرة بشدة بدرجة الحرارة. تؤكد الدراسة أن درجة الحرارة هي "المخرج" القوي في الحياة المبكرة للسمكة، حيث تشكل كيفية نموها وتطورها. ورغم أن الورقة البحثية لا تدعي أنها حلت كل الأسرار المتعلقة بتربية الأسماك، إلا أنها توفر خريطة عالية الجودة ودقيقة للتغيرات الجينية التي تحدث عندما تواجه هذه الأسماك درجات حرارة مختلفة. وهذه الخريطة متاحة الآن للعلماء الآخرين لاستخدامها، مما يساعدهم في معرفة كيفية تربية أسماك أكثر صحة وفهم كيفية تكيف هذه الكائنات المذهلة مع عالم متغير.
ملخص تقني: مجموعة بيانات شاملة لتسلسل الحمض النووي الريبوزي (RNA-seq) للنمو الجنيني واليرقي في سمك موسى السنغالي (Solea senegalensis) عبر ثلاثة أنظمة حرارية
المشكلة والسياق يعد سمك الموسى السنغالي (Solea senegalensis) نوعاً ذا قيمة تجارية عالية، إلا أن تربيته في المزارع السمكية تواجه تحديات كبيرة خلال مراحل الحياة المبكرة، بما في ذلك ارتفاع معدلات النفوق، وتفاوت النمو، والتشوهات الهيكلية. وترتبط هذه المشكلات غالباً بالعوامل البيئية، حيث تم تحديد درجة الحرارة كمحرك رئيسي لمعدل النمو، وأداء اليرقات، وأنماة التعبير الجيني. وبينما توسعت الموارد الجينومية لسمك S. senegalensis في السنوات الأخيرة — بما في ذلك التجميعات على مستوى الكروموسومات ودراسات تحديد الجنس — لا تزال هناك فجوة حرجة في البيانات المتاحة علنًا لنسخ الحمض النووي الريبوزي (transcriptomic data) التي تغطي بشكل شامل مراحل تطورية متعددة تحت ظروف حرارية متباينة ومتحكم بها. ركزت الدراسات الحالية إلى حد كبير على جوانب محددة مثل التكيف الحراري أو التحول الشكلي بمعزل عن غيرها، مما يفتقر إلى وجود مجموعة بيانات موحدة تمتد عبر الانتقال الكامل من التكوين الجنيني إلى اليرقات ما بعد التحول.
المنهجية لمعالجة هذه الفجوة، أنتج المؤلفون مجموعة بيانات عالية الجودة لـ RNA-seq توضح التطور الجنيني واليرقي تحت ثلاثة أنظمة درجات حرارة ثابتة: 14 درجة مئوية، و18 درجة مئوية، و21 درجة مئوية. تضمن التصميم التجريبي أخذ عينات عبر ثماني مراحل تطورية متميزة: ست مراحل جنينية (المرحلة البلستية، والمرحلة العتامية، ومراحل التجزئة المبكرة والمتوسطة والمتأخرة، ومرحلة ما قبل الفقس) ومرحلتين يرقيتين (ما قبل التحول وما بعد التحول).
جمع العينات: تم إنشاء ما مجموعه 48 مكتبة RNA-seq، تضم مكررات بيولوجية اثنين لكل مرحلة من المراحل الثماني عبر ظروف درجات الحرارة الثلاثة. جُمعت العينات في مركز "إسيمات" (ECIMAT) للبحوث البحرية، وحُفظت في محلول RNAlater، وخُزنت في درجة حرارة -80 مئوية.
التسلسل: استُخلص الحمض النووي الريبوزي الكلي (Total RNA)، وجُهزت المكتبات، وسُلسلت باستخدام منصة Illumina NovaSeq 6000، مما أنتج حوالي 30 مليون قراءة مزدوج النهاية بطول 150 زوجاً قاعدياً لكل مكتبة.
معالجة البيانات: تمت معالجة القراءات الخام باستخدام مسار nf-core/rnaseq (الإصدار 3.10.1). أُجري ضبط الجودة باستخدام FastQC وTrimmomatic. تمت محاذاة القراءات إلى المرجع الجينومي لـ S. senegalensis باستخدام STAR، وقُدر وفرة النسخ باستخدام RSEM لإنتاج العد الخام وقيم "النسخ لكل مليون" (TPM).
التحليل الإحصائي: أُجريت التحليلات اللاحقة في بيئة R (الإصدار 4.3.0). تم تحليل التعبير التفاضلي باستخدام DESeq2، مع اتخاذ حالة 18 درجة مئوية كمجموعة مرجعية. اعتُبرت الجينات ذات التعبير التفاضلي (DEGs) معتبرة إذا كانت نسبة الخطأ في الاكتشاف (FDR) أقل من 0.05. كما أُجريت تحليلات الإثراء الوظيفي (مصطلحات الوجود الجيني - Gene Ontology - ومسارات KEGG) باستخدام clusterProfiler، مع ربط نسخ S. senegalensis بنظائرها في سمك الدانيو ريريو (Danio rerio).
النتائج الرئيسية أنتجت الدراسة مجموعة بيانات قوية تتميز بجودة تسلسل عالية وتكرار بيولوجي متسق.
جودة البيانات والتجميع: أكد تحليل المكونات الرئيسية (PCA) وتحليل ارتباط سبيرمان أن العينات تجمعت أساساً حسب مرحلة التطور، مع تجمع المكررات البيولوجية بشكل وثيق. يشير هذا إلى قابلية عالية للتكرار ويجسد عمليات إعادة تشكيل النسخ الواسعة التي تحدث طوال فترة التطور. لوحظ تجميع خاص بدرجة الحرارة في المراحل المبكرة (البلستية، والعتامية، والتجزئة المبكرة) ويرقات ما قبل التحول، بينما أظهرت المراحل اللاحقة تداخلاً أكبر بين المعالجات الحرارية.
ديناميكيات التعبير التفاضلي: تباين عدد الجينات ذات التعبير التفاضلي (DEGs) بشكل كبير عبر مراحل التطور. أظهرت المرحلة البلستية أعلى عدد من الجينات المستجيبة للحرارة (6,983 جين تفاضلي لـ 14 درجة مئوية مقابل 18 درجة مئوية؛ و6,923 لـ 21 درجة مئوية مقابل 18 درجة مئوية). وفي المقابل، أظهرت مرحلة ما قبل الفقس أقل عدد من الجينات التفاضلية. ومن الجدير بالذكر أنه بينما سجلت المرحلة البلستية أعلى عدد من الجينات، أظهرت مرحلتا التجزئة المتأخرة وما قبل الفقس أعلى متوسط لقيم لوغاريتم (log₂) التغير المطلق.
الجينات المرشحة المتكررة: أدت عملية تصفية الجينات بناءً على الأهمية الإحصائية وحجم التغير إلى تحديد مجموعة من الجينات المتكررة المستجيبة للحرارة. كان الجين الأكثر اكتشافاً هو cyp8b1 (عضو في عائلة سيتوكروم P450 تحت العائلة الفرعية B رقم 1)، والذي ظهر في خمس مقارنات مستقلة. وشملت الجينات الأخرى التي تم تحديدها بشكل متكرر بروتينات الصدمة الحرارية (hsc70, hsp90aa1.2, hspa8b, hsp71) والجين rbp2a.
الإثراء الوظيفي:
الوجود الجيني (GO): تم تحديد 278 مصطلحاً مثرياً بشكل معنوي عبر جميع المقارنات، تغطي مجالات العمليات البيولوجية، والوظيفة الجزيئية، والمكون الخلوي.
مسارات KEGG: تم اكتشاف 129 مساراً مثرياً بشكل معنوي. شملت المسارات المثيرة بشكل متكرر إشارات الكاد هيرين (Cadherin signaling)، والوصلات الالتصاقية (Adherens junction)، والوصلات المحكمة (Tight junction)، وإشارات IgSF CAM، وإشارات TGF-beta، والفسفرة التأكسدية (Oxidative phosphorylation).
الأهمية والادعاءات يضع المؤلفون هذه المجموعة كمرجع أساسي يسد فجوة حرجة في بيانات النسخ المتاحة علناً لتطور الأسماك المفلطحة في مراحلها المبكرة. ويدعي البحث أن هذا المورد يتيح:
التوصيف الشامل: رؤية تفصيلية لديناميكيات التعبير الجيني الخاصة بالمرحلة ودرجة الحرارة من المرحلة البلستية إلى يرقات ما بعد التحول.
استقصاء اللدونة الحرارية: أساس لدراسة كيفية دفع درجة الحرارة للتباين المظهري والتنظيم التطوري في سمك S. senegalensis.
دعم الاستزراع السمكي: بيانات يمكن أن تدعم الأبحاث المتعلقة بالتباين المظهري المدفوع بالحرارة، مما قد يساعد في تطوير مؤشرات حيوية للإجهاد الحراري أو جودة التطور في بيئات الاستزراع السمكي.
البيولوجيا المقارنة: مورد للبيولوجيا المقارنة للأسماك المفلطحة ولدمج مجموعات البيانات "الأوميكية" الأخرى (مثل الإبيجينوميات، والبروتيوميات، والميتابولوميات).
يصرح المؤلفون صراحةً بأن مجموعة البيانات مناسبة لتوصيف ديناميكيات التعبير الجيني، وتحديد جينات الاستجابة الحرارية، وتطوير المؤشرات الحيوية، مع الامتناع عن اقتراح بروتوكولات تجريبية مستقبلية محددة بخلاف هذه التطبيقات المذكورة. وتتوفر البيانات، بما في ذلك القراءات الخام، ومصفوفات العد، ونتائج الإثراء، عبر مستودع Figshare لتسهيل هذه المسارات البحثية.