Hybrid assembly of the African savanna elephant whole blood transcriptome using short- and long-read data
تُنتج هذه الدراسة نسخة شاملة وعالية الجودة من النسخ المتماثل للترانسكريبتوم في الدم الكامل لفيل السافانا الأفريقي عبر مقارنة طرق تجميع مختلفة وتحديد أداة StringTie2 مع التجميع الموجه بالتوصيف كنهج أمثل، مما يؤسس لمورد حيوي لأبحاث الأمراض، والجينوم الوظيفي، وأبحاث الحفظ المستقبلية لهذا النوع المهدد بالانقراض.
المؤلفون الأصليون:Stacey Caitlin Engel, Tanya Jane Kerr, Gian van der Spuy, Peter Erik Buss, Léanie Kleynhans, Michele Ann Miller, Tracey Jooste
في المناظر الطبيعية الشاسعة والمشمسة للسافانا الأفريقية، يقف الفيل كعملاق لطيف، ومع ذلك فهو يواجه تهديدات غير مرئية يمكن أن تكون قاتلة مثل أي مفترس. تشكل أمراض مثل السل وحالة فيروسية نزفية شديدة مخاطر كبيرة على هذه المجموعات المهددة بالانقراض، لكن العلماء جاهدوا لفهم كيفية مقاومة الأجهزة المناعية للفيلة لهذه الأمراض. ويعود ذلك إلى حد كبير لافتقارهم إلى خريطة حاسمة: قائمة مفصلة بالجينات النشطة في دم الفيل. في علم الأحياء، تسمى هذه القائمة "الترانسكريبتوم" (transcriptome). فكر في الأمر كقائمة فورية بالتعليمات التي تستخدمها الخلية حاليًا لبناء البروتينات، مما يكشف عن الدفاعات التي في حالة استنفار وتلك التي في حالة راحة. وبدون هذه القائمة لفيل السافانا الأفريقي، كان الباحثون يحاولون التنقل في لغز طبي معقد باستخدام رسم تخطيطي غامض فقط، غير قادرين على رؤية الآليات الجزيئية المحددة التي تحدد ما إذا كان الحيوان سينجو من العدوى أم لا.
ولحل هذه المشكلة، شرع فريق من الباحثين في بناء أول ترانسكريبتوم عالي الجودة مستخلص خصيصًا من الدم الكامل لفيلة السافانا الأفريقية التي تعيش في البرية. لقد أدركوا أنه للحصول على صورة واضحة، يحتاجون إلى النظر إلى المادة الوراثية من زاويتين مختلفتين. إحدى الطرق، باستخدام تسلسل القراءة القصيرة (short-read sequencing)، تشبه قراءة كتاب عبر أخذ لقطات صغيرة ومثالية للكلمات الفردية؛ فهي دقيقة للغاية ولكن من الصعب من خلالها رؤية كيفية تدفق الجمل معًا. أما الطريقة الأخرى، باستخدام تسلسل القراءة الطويلة (long-read sequencing)، فهي تشبه قراءة فقرات طويلة ومتصلة؛ فهي تلتقط البنية الكاملة للقصة ولكنها قد تحتوي أحيانًا على بعض الأخطاء المطبعية. ومن خلال الجمع بين هذين النهجين، هدف العلماء إلى إنشاء خريطة جينية كاملة ودقيقة يمكن أن تعمل كمرجع للدراسات المستقبلية حول صحة الفيلة.
جمع الفريق عينات دم من أربعة ذكور بالغين في حديقة كروجر الوطنية بجنوب أفريقيا. واستخرجوا بعناية الحمض النووي الريبوزي (RNA)، وهو الجزيء الذي يحمل التعليمات الوراثية، وجهزوه للتسلسل باستخدام تقنيات القراءة القصيرة والطويلة. وبمجرد توليد البيانات، واجه الباحثون تحديًا حاسوبيًا هائلًا: تجميع ملايين الشظايا الجينية في وحدة متماسكة. لقد اختبروا عدة برامج حاسوبية مختلفة لمعرفة أي منها يمكنه إعادة بناء القصة الجينية للفيل بشكل أفضل. حاولت بعض البرامج بناء الخريطة من الصفر دون أي مساعدة، بينما استخدمت أخرى الجينوم الموجود للفيل، رغم كونه غير مكتمل، كدليل. قام الباحثون بتشغيل هذه البرامج في تركيبات متنوعة، مع تغذيتها بمزيج مختلف من البيانات القصيرة والطويلة لمعرفة أي نهج ينتج النتيجة الأكثر موثوقية.
بعد معالجة الأرقام، برزت طريقة واحدة كفائزة واضحة. تم إنشاء أفضل تجميع بواسطة برنامج يسمى StringTie2، ولكن مع لمسة خاصة: فقد كان موجهًا بالتعليقات التوضيحية الجينية المعروفة للفيل. أنتج هذا النهج ترانسكريبتوم كان مكتملًا بشكل ملحوظ، حيث التقط 98.4 بالمائة من الجينات الأساسية المتوقع وجودها في الثدييات. في المقابل، أنتجت الطرق التي حاولت بناء الخريطة من الصفر دون دليل نتائج أكثر ضجيجًا، مع وجود العديد من القطع المفقودة والأجزاء المكررة. أظهر التجميع الفائز أن الباحثين نجحوا في إعادة بناء خيوط طويلة ومتصلة من التعليمات الجينية، بدلاً من مجرد كومة من الشظايا القصيرة. وكشفت هذه الخريطة عالية الجودة أن عينات الدم احتوت على مجموعة متنوعة من الجينات النشطة المشاركة في الاستجابات المناعية، والتنظيم الخلوي، والعمليات الأيضية، وهو بالضبط ما يتوقع وجوده في حيوان حي يتنفس يواجه عدوى محتملة.
لم تقتصر الدراسة على مجرد إنشاء خريطة؛ بل حددت أيضًا آلاف التسلسلات الجينية الجديدة التي لم تكن موجودة في البيانات المرجعية السابقة. تشير هذه النسخ الجديدة إلى أن الفهم الحالي لجينوم الفيل غير مكتمل، وأن هذا المورد الجديد القائم على الدم يقدم صورة أغنى لبيولوجيا الحيوان. وبينما كانت الدراسة محدودة بأربعة ذكور بالغين، وقد يكون لضغط عملية الإمساك بالحيوان تأثير على بعض نشاط الجينات، فإن عمق البيانات ومتانة التجميع يوفران أساسًا متينًا للمستقبل. يقدم هذا العمل أداة حيوية للأطباء البيطريين وعلماء الحفاظ على البيئة، مما يسم يسمح لهم أخيرًا بالتحقيق في كيفية استجابة هذه الحيوانات المهيبة للأمراض على المستوى الجزيئي. ومع وجود هذا المرجع الجديد في أيديهم، يمكن للعلماء الآن البدء في كشف أسرار مناعة الفيلة، مما قد يؤدي إلى طرق أفضل لحماية هذه العمالقة المهددة بالانقراض من الأمراض التي تهدد بقاءها.
ملخص تقني: التجميع الهجين لترانسكريبتوم الدم الكامل للفيل الأفريقي فيل السافانا
بيان المشكلة لقد تعرقلت الأبحاث المتعلقة بالاستجابات المناعية واكتشاف المؤشرات الحيوية في الفيل الأفريقي فيل السافانا (Loxodonta africana) بشكل كبير بسبب ندرة الموارد الترانسكريبتومية. وبينما يعد التوصيف الترانسكريبتومي أداة قياسية لفهم الآليات المناعية للمضيف والنشوء المرضي للأمراض في البشر والنماذج البيطرية الأخرى، إلا أن الفيل يفتقر إلى ترانسكريتوم مرجعي شامل وعالي الجودة. تقتصر البيانات الموجودة حاليًا على مجموعة بيانات واحدة مستمدة من الأنسجة المشيمية، مما يحد من فائدتها في دراسة الاستجابات المناعية الجهازية للتهديدات الحرجة مثل السل (TB) ومرض هستيروسيس (EEHV) النزفي. علاوة على ذلك، فإن الاستراتيجية المثلى لتجميع الترانسكريبتومات في الأنواع غير النموذجية — وتحديدًا ما إذا كان ينبغي استخدام نهج موجه بالمرجع أو نهج التجميع من الصفر (de novo)، وكيفية دمج بيانات القراءة القصيرة والطويلة بأفضل طريقة — تتطلب التحقق من صحتها في هذا السياق البيولوجي المحدد.
المنهجية أنتجت الدراسة ترانسكريبتوم عالي الجودة للدم الكامل باستخدام نهج تجميع هجين يجمع بين تقنيات تسلسل القراءات القصيرة والطويلة.
جمع العينات: تم استخراج الحمض النووي الريبوزي (RNA) من الدم الكامل لأربعة فيلة ذكور بالغة حرة في متنزه كروجر الوطني بجنوب أفريقيا. خضعت العينات لعملية استنفاد لـ mRNA الخاص بالجلوبين و الـ rRNA لإثراء الـ mRNA.
التسلسل:
القراءات القصيرة: تم إنتاجها باستخدام منصة MGI Tech (DNBSEQ-G400RS)، مما أنتج قراءات مزدوجة النهاية بطول 100 زوج من القواعد.
القراءات الطويلة: تم إنتاجها باستخدام Oxford Nanopore Technologies (ONT) PromethION 2 Solo مع مجموعة أدوات تسلسل cDNA-PCR V14 Barcoding Kit. تم إجراء استدعاء القواعد (Base calling) باستخدام نموذجي الدقة العالية (HAC) والدقة الفائقة (SUP).
سير العمل المعلوماتي الحيوي:
المعالجة المسبقة: تم فحص جودة القراءات (FastQC)، وتقليمها (Trimmomatic لـ NGS؛ وPorechop/NanoFilt لـ TGS)، ومحاذاتها مع الجينوم المرجعي لفيل السافانا الأفريقي (mLoxAfr1.hap2) باستخدام STAR (لـ NGS) وMinimap2 (لـ TGS).
استراتيجيات التجميع: تم تقييم مسارات تجميع متعددة:
التجميع من الصفر (De novo): rnaSPAdes (باستخدام كل من الأنماط الهجينة والقراءة القصيرة فقط)، وRNA-Bloom2، وTrinity (قراءة قصيرة فقط).
التقييم: تم تقييم التجميعات باستخدام rnaQUAST للإحصاءات الهيكلية، وBUSCO لاكتمال المحتوى مقابل سلالة mammalia_odb12، والتوصيف الوظيفي عبر eggNOG-mapper. تم التنبؤ بالمناطق المشفرة باستخدام TransDecoder، وتم تقييم التكرار عبر CD-HIT.
النتائج الرئيسية
التجميع الأمثل: برز التجميع الموجه بالتوصيف لـ StringTie2 كأعلى ترانسكريبتوم جودة. حيث حقق درجة اكتمال BUSCO بلغت 98.4%، وهي الأعلى بين جميع الطرق المختبرة.
الجودة الهيكلية: أظهر التجميع الموجه بالتوصភាព لـ StringTie2 استمرارية فائقة في الترانسكريبت، مع N50 قدره 13,206 زوج قاعدي، و8 نسبة مئوية 80.3% من الترانسكريبتات التي يتجاوز طولها 1,000 زوج قاعدي. كما أظهر أقل نسبة من الترانسكريبتات غير الموصوفة (11.7%) وأعلى نسبة من الترانسكريبتات ذات الأطر القرائية المفتوحة المتوقعة (83.3%).
مقارنة النهج:
الموجه بالمرجع مقابل التجميع من الصفر: أنتجت النهج الموجه بالمرجع عمومًا ترانسكريبتات أقل عددًا، وأطول، وأكثر استمرارية مع نسب تكرار أقل مقارنة بمجمعات التجميع من الصفر (Trinity، rnaSPAdes)، والتي أنتجت تجميعات أكبر مع معدلات تكرار أعلى وأطوال متوسطة أقصر للترانسكريبتات.
فائدة البيانات الهجينة: بينما لم يؤدِ إضافة بيانات القراءة الطويلة إلى تجميع rnaSPAdes من الصفر إلى تغيير جذري في إجمالي حجم التجميع، إلا أنه زاد من نسبة الترانسكريبتات التي تزيد عن 1 كيلو قاعدة وحسّن من استعادة الجينات والأنماط المتعددة (isoforms) المجمعة بنسبة 95% مقارنة بالتجميع المعتمد على القراءات القصيرة فقط.
RNA-Bloom2: أنتجت هذه الأداة تجميعًا مدمجًا، لكنها عانت من انخفاض الاكتمال (13.1% BUSCO) ومعدل مرتفع من الجينات المفقودة.
التوصيف الوظيفي: احتوى الترانسكريبتوم الناتج على عمليات بيولوجية متنوعة، ووظائف جزيئية، ومسارات ذات صلة بالتمثيل الغذائي، والسرطان، والعدوى الفيروسية، والتنظيم الخلوي والمناعي. سلط تحليل KEGG الضوء على المسارات المتعلقة بالتمثيل الغذائي، والسرطان، والعدوى الفيروسية، والإشارات (مثل MAPK، وPI3K-Akt).
الجدة (Novelty): حددت الدراسة 2,231 جينًا جديدًا محتملاً و2,342 ترانسكريبتًا جديدًا محتملاً لم يتداخلوا مع التوصيفات المرجعية الموجودة، مما يشير إلى أن هذا المورد يلتقط ملفات تعبير فريدة من الدم الكامل غير الموجودة في مرجع المشيمة.
الأهمية والادعاءات يزعم البحث أن الترانسكريبتوم الموجه بالتوصيف لـ StringTie2 الذي تم إنتاجه يعمل كمورد مرجعي قوي للتحقيقات المستقبلية في آليات المرض (تحديدًا السل وEEHV) والاستجابات المناعية في فيلة السافانا الأفريقية. ويؤكد المؤلفون أن هذا المورد يوفر أساسًا ضروريًا لدراسات التعبير الجيني، والجينوميات الوظيفية، وتوصيف الجينوم في هذا النوع المهدد بالانقراض.
بالإضافة إلى ذلك، تصدق الدراسة على سير عمل محدد لتجميع الترانسكربتوم في الأنواع غير النموذجية حيث يوجد جينوم مرجعي ولكن ليس مثاليًا. تشير النتائج إلى أنه عند توفر جينوم مرجعي، فإن النهج الموجه بالتوصيف يعطي استمرارية واكتمالًا فائقين مقارنة بطرق التجميع من الصفر. كما يوضح العمل قيمة استراتيجيات التسلسل الهجين، مشيرًا إلى أنه بينما قد لا تغير بيانات القراءة الطويلة مقاييس التجميع بشكل كبير في سياق الموجه بالمرجع، إلا أنها تساهم في استعادة ترانسكريبتات وأنماط متعددة (isoforms) أطول عند استخدامها في سياقات التجميع من الصفر. ويخلص المؤلفون إلى أن هذا المورد، المستمد من الدم الكامل، يقدم صورة أكثر شمولاً للنشاط المناعي الجهازي مقارنة بمجموعات البيانات السابقة الخاصة بأنسجة معينة.