Exploration of key genes associated with pyruvate metabolism in nonalcoholic fatty liver disease based on single-cell and bulk RNA-seq data
تدمج هذه الدراسة بين تسلسل الحمض النووي الريبي (RNA) للكتلة وتسلسل الحمض النووي الريبي للخلية الواحدة مع التحقق التجريبي لتحديد جينَي ACOT12 وSLC16A1 كجينات محورية مرتبطة باستقلاب البيروفات في مرض الكبد الدهني غير الكحولي، مما يوضح أدوارهما في الارتشاح المناعي وتطور المرض مع اقتراح جزيئات صغيرة علاجية محتملة.
يعتبر الكبد بمثابة محطة المعالجة الكيميائية في الجسم، حيث يقوم باستمرار بتصفية الدم، وتخزين الطاقة، وتفكيك السموم. وللحفاظ على سير هذه العمليات بسلاسة، فإنه يعتمد على تدفق مستمر للوقود، تماماً كما تحتاج المصانع إلى إمداد موثوق من المواد الخام. وأحد أهم أنواع الوقود لهذه العملية هو "البيروفات"، وهو جزيء صغير ينتج عندما يفكك الجسم السكريات والكربوهيدرات. وفي الكبد السليم، يتم تحويل البيروفات بكفاءة إلى طاقة أو يُستخدم لبناء جزيئات جديدة. ومع ذلك، في حالة تُعرف باسم مرض الكبد الدهني غير الكحولي، يتعطل خط الوقود هذا؛ حيث يغرق الكبد بالدهون، مما يؤدي إلى الالتهاب والتندب، وهو تطور يمكن أن يتسبب في النهاية في أضرار جسيمة للعضو. وبينما يعلم الأطباء أن الاختلالات الأيضية هي التي تقود هذا المرض، إلا أن المفاتيح الجزيئية المحددة التي تحول الكبد السليم إلى كبد مريض ظلت صعبة التحديد. ويعد فهم هذه المفاتيح أمراً ضرورياً لأن العلاجات الحالية تركز في الغالب على تغيير نمط الحياة، مما يترك فجوة للعلاجات المستهدفة التي يمكنها إيقاف المرض في مساره.
لقد وضع فريق من الباحثين في جامعة كونمينغ الطبية هدفاً يتمثل في العثور على هذه المفاتيح المفقودة من خلال البحث في التعليمات الجينية للكبد. وقد جمعوا بين نهجين قويين: تحليل مجموعات كبيرة من البيانات الجينية من العديد من المرضى، وفحص نشاط الخلايا الفردية داخل الكبد. ومن خلال التركيز تحديداً على الجينات المشاركة في عملية أيض البيروفات، بحثوا عن التعليمات المحددة التي تسير بشكل خاطئ عندما يبدأ الكبد في تراكم الدهون. وقد قادتهم تحقيقاتهم إلى جينين محددين، هما ACOT12 وSLC16A1، واللذان بدوا كلاعبين مركزيين في المرض. ففي المرضى الذين يعانون من مرض الكبد الدهني، كانت التعليمات الخاصة بجين ACOT12 مرتفعة للغاية، بينما كانت التعليمات الخاصة بجين SLC16A1 منخفضة. هذا السلوك المتعاكس يشير إلى أن الكبد يحاول التعويض عن أزمة أيضية، ربما عن طريق تفكيك الدهون الزائدة بينما يكافح لنقل الوقود من وإلى الخلايا بكفاءة.
لم يتوقف الباحثون عند تحديد هذه الجينات، بل أرادوا فهم كيفية تفاعلها مع بيئة الكبد. فقد اكتشفوا أن هذين الجينين يبدو أنهما يؤثران على استجابة الجهاز المناعي للمرض. فجين SLC16A1 الذي انخفض نشاطه، بدا مرتبطاً بوجود أعلى لبعض الخلايا المناعية التي تجوب الكبد، بينما ارتبط الجين المرتفع ACOT12 بوجود عدد أقل من هذه الخلاوة. وهذا يشير إلى أن الحالة الأيضية لخلايا الكبد تتواصل مباشرة مع الجهاز المناعي، مما قد يؤدي إلى الالتهاب الذي يحول مجرد تراكم الدهون إلى كبد أكثر خطورة ومصاب بالتندب. ولرؤية أين يحدث هذا بالضبط، فحص الفريق الكبد خلية تلو الأخرى، ووجدوا أن نوعاً معيناً من الخلايا، يُعرف باسم "الخلية النجمية الكبدية"، هو الموقع الرئيسي الذي تنشط فيه هذه التغييرات الجينية بشكل أكبر. هذه الخلايا تكون هادئة في العادة، ولكن عندما تشعر بالخطر، تستيقظ وتبدأ في بناء الأنسجة الندبية، وهي عملية تؤدي إلى تليف الكبد. وقد أظهرت الدراسة أنه مع زيادة نشاط هذه الخلايا وتحركها نحو حالة تكوين الندبات، تتغير مستويات هذين الجينين الرئيسيين وفق نمط يمكن التنبؤ به، مما يعزز دورهما في تطور المرض.
ولاختبار ما إذا كانت هذه النتائج ستصمد في نظام حي، استخدم الباحثون أسماك الزيبرا وخلايا كبد بشرية في طبق مخبري. فقد عرضوا هذه النماذج لمستويات عالية من الدهون، محاكيةً بذلك ظروف المرض. وكما توقعت النماذج الحاسوبية، أظهرت الأسماك والخلايا زيادة في تخزين الدهون ونفس التحولات في النشاط الجيني: ارتفع ACOT12 وانخفض SLC16A1. وقد أكد هذا أن الأنماط الجينية الملاحظة في المرضى البشريين لم تكن مجرد ضجيج إحصائي، بل تعكس تغيرات بيولوجية حقيقية. ثم طرح الفريق سؤالاً جديداً: هل يمكن لأي شيء أن يوقف هذه العملية؟ وباستخدام المحاكاة الحاسوبية، قاموا بفحص آلاف المركبات الطبيعية الموجودة في المحيط لمعرفة ما إذا كان أي منها يمكن أن يرتبط بالبروتينات التي تنتجها هذه الجينات. وقد حددوا جزيئين محددين من السلاسل الطويلة اللذين يتناسبان بإحكام مع هياكل البروتينات المستهدفة، تماماً مثل مفتاح يتناسب مع قفل. وأظهرت هذه الجزيئات، المستمدة من مصادر بحرية، وعوداً في النماذج الحاسوبية كأدوات محتملة للتأثير على نشاط هذه الجينات، مما يفتح مساراً جديداً لتطوير الأدوية.
تخلص الدراسة إلى أن اختلال أيض البيروفات ليس مجرد عرض جانبي لمرض الكبد الدهني، بل هو جزء جوهري من آليته، مدفوعاً بالعمليات المتضادة لـ ACOT12 وSLC16A1. ويبدو أن هذين الجينين يعملان كجسر بين كيفية تعامل الكبد مع الطاقة وكيفية استجابته للإشارات المناعية والتندب. وبينما لا تزال الأدوية المحتملة التي تم تحديدها في المراحل الأولى من الاكتشاف وتتطلب مزيداً من الاختبارات لإثبات فعاليتها في البشر، فإن البحث يوفر خارطة طريق واضحة للاعبين الجزيئيين المعنيين. ومن خلال تحديد هذه الجينات المحددة والخلايا التي تعتمد عليها، تقدم الدراسة هدفاً جديداً للعلاجات المستقبلية، مما يقرب المجال من الوصول إلى علاجات يمكنها معالجة الأسباب الجذرية لتلف الكبد بدلاً من مجرد إدارة الأعراض.
ملخص تقني: استكشاف الجينات الرئيسية المرتبطة باستقلاب البيروفات في مرض الكبد الدهني غير الكحولي
بيان المشكلة يُعد مرض الكبد الدهني غير الكحولي (NAFLD)، والذي أُعيد تصنيفه مؤخرًا تحت إطار مرض الكبد الدهني المرتبط بالخلل الوظيفي الأيضي (MASALD)، حالة كبدية مزمنة منتشرة تتميز بترسب الدهون، والالتهاب، والقدرة على التطور إلى تليف. وبينما تُعتبر مقاومة الأنسولين والاختلال الأيضي من الدوافع الراسخة، فإن المنظمات الجزيئية المحددة التي تربط إعادة البرمجة الأيضية — وخاصة استقلاب البيروفات — بالتليف وتنشيط الجهاز المناسي لا تزال غير محددة بشكل كامل. تركز العلاجات الحالية على عوامل الخطر الأيضية ولكنها تفتقر إلى الفعالية الكافية في منع التليف. تهدف هذه الدراسة إلى تحديد الجينات الرئيسية المرتبطة باستقلاب البيروفات في مرض NAFLD، وتوضيح سياقاتها الخلوية، واستكشاف الأهداف العلاجية المحتملة.
المنهجية استخدمت الدراسة نهجًا متكاملًا متعدد الأوميكس يجمع بين بيانات تسلسل الحمض النووي الريبي للكتلة (bulk RNA-seq) وتسلسل الحمض النووي الريبي أحادي الخلية (scRNA-seq)، متبوعًا بالتحقق التجريبي والمسح الحاسوبي للأدوية.
الحصول على البيانات ومعالجتها:
تسلسل الحمض النووي الريبي للكتلة (Bulk RNA-seq): تم استرداد الملفات التعريفية النسخية من قاعدة بيانات GEO (المجموعة التدريبية GSE126848، ومجموعة التحقق GSE89632).
تسلسل الحمض النووي الريبي أحادي الخلية (scRNA-seq): تم تحليل البيانات من GSE212837 (11 عينة كبد لمرضى NAFLD و3 عينات ضابطة) لتوضيح التغاير الخلوي.
مجموعات الجينات: تم تجميع الجينات المرتبطة باستقلاب البيروفات (PMRGs) من مسار KEGG في MSigDB (97 جينًا).
التحليل المعلوماتي الحيوي:
اختيار المرشحين: تم تحديد الجينات المعبر عنها تفاضليًا (DEGs) باستخدام DESeq2. وتقاطع جينات استقلاب البيروفات (PMRGs) مع الـ DEGs لتشكيل مجموعة مرشحة.
تعلم الآلة: استُخدمت ثلاث خوارزميات — انحدار LASSO، وشجرة القرار المعززة بالتدرج (GBDT)، وتعزيز التدرج الشديد (XGBoost) — لتنقية الجينات المميزة. وأسفر تقاطع مخرجات هذه الخوارزميات عن مجموعة الجينات النهائية. �����* الأداء التشخيصي: تم تقييم قدرة التشخيص باستخدام تحليل منحنى خصائص تشغيل المستقبِل (ROC) ومنطقة تحت المنحنى (AUC) للتحقق من الإمكانات التشخيصية.
التحليل الوظيفي والمناعي: تم إجراء تحليل إثراء مجموعة الجينات (GSEA)، وبناء شبكة التفاعل بين البروتينات (PPI) باستخدام STRING، وتوصيف ارتشاح الخلايا المناعية (GSVA). كما ربطت تحليلات الارتباط بين الجينات الرئيسية والخلايا المناعية المرتشحة تفاضليًا (DICs).
الشبكات التنظيمية: تم بناء شبكات التعبير المشترك، وmiRNA-mRNA، وlncRNA-miRNA-mRNA باستخدام قواعد بيانات مثل GeneMANIA وmultiMiR وStarBase.
تحليل الخلية الواحدة (Single-Cell Analysis): تمت معالجة بيانات scRNA-seq باستخدام Selet لتقسيم الخلايا وتصنيفها. كما طُبق تحليل المسار الزمني الكاذب (Monocle) واستدلال التواصل بين الخلايا (CellChat) على الخلايا النجمية الكبدية (HSCs). وتم تقييم نشاط العوامل النسخية (TF) ونشاط المسارات الأيضية باستخدام DoRothEA وscMetabolism وReactomeGSA.
التحقق التجريبي:
النماذج: استُخدم نموذج سمك الزيبرا المحفز بنظام غذائي عالي الدهون، ونموذج خلايا LX-2 (خلايا نجمية كبدية بشرية) المعالجة بحمض البالمتيك (PA).
الاختبارات: استُخدم صبغ ألوان الزيت الحمراء (Oil Red O)، وصبغة الهيماتوكسيلين والإيوسين (H&E)، وصبغة ماسون (Masson trichrome)، وتفاعل البوليميراز المتسلسل الكمي (qPCR) للتحقق من تراكم الدهون، والتغيرات النسيجية، ومستويات التعبير عن mRNA لكل من ACOT12 وSLC16A1.
اكتشاف الأدوية:
المسح: قام نموذج GraphBAN، المدرب على مجموعات بيانات التفاعل بين المركب والبروتين، بمسح قاعدة بيانات المنتجات الطبيعية البحرية الشاملة (CMNPD) بحثًا عن روابط محتملة.
التصفية: تمت تصفية المرشحين بناءً على خصائص ADME (امتصاص معوي عالٍ، الامتثال لقاعدة ليبينسكي، وانخ้อย من تنبيهات PAINS).
تحديد الجينات الرئيسية: من خلال دمج تعلم الآلة والتحقق من ROC، تم تحديد ACOT12 وSLC16A1 كجينات محورية مرتبطة باستقلاب البيروفات في NAFLD.
أنماط التعبير: لوحظ ارتفاع معنوي في تعبير ACOT12 في عينات NAFLD، بينما انخفض تعبير SLC16A1 بشكل ملحوظ. وأظهر كلاهما اتجاهات متسقة عبر مجموعات التدريب والتحقق مع AUC > 0.7.
الرؤى الوظيفية:
المسارات: كشف تحليل GSEA أن كلا الجينين يلتقيان في مسار البروتيزوم. كما ارتبط ACOT12 إضافيًا بمسارات التوازن البروتيني (مثل الترجمة، والاستهداف المشترك للترجمة المعتمد على SRP).
الارتباط المناعي: أظهر SLC16A1 ارتباطات إيجابية مع معظم الخلايا المناعية المرتشحة تفاضليًا (DICs)، بما في ذلك خلايا B الذاكرة والخلايا التائية التنظيمية. وفي المقابل، أظهر ACOT12 ارتباطات سلبية مع هذه المجموعات المناعية.
التمركز داخل الخلية: يُتوقع أن يكون ACOT12 متمركزًا في السيتوبلازم/الغشاء السيتوبلازمي، بينما يتمركز SLC16A1 بشكل أساسي في الغشاء البلازمي والشبكة الإندوبلازمية.
السياق الخلوي (scRNA-seq):
المجتمع الخلوي الرئيسي: تم تحديد الخلايا النجمية الكبدية (HSCs) كعنصر خلوي رئيسي يساهم في تقدم NAFLD، حيث أظهرت تغيرات كبيرة في النسب وتعبير الجينات.
ديناميكيات HSC: أشار تحليل المسار الزمني الكاذب إلى أن تعبير ACOT12 يزداد أثناء تمايز HSC (من الحالة الساكنة إلى الحالة النشطة)، بينما يظل SLC16A1 منخفضًا ومستقرًا.
التواصل: كشف تحليل التواصل بين الخلايا عن تعزيز التفاعلات بين HSC والخلايا الكبدية، وانخفاض التفاعلات بين HSC والخلايا اللمفاوية التائية في عينات NAFLD مقارنة بالعينات الضابطة.
التحقق التجريبي:
في يرقات سمك الزيبرا وخلايا LX-2، تأكد ارتفاع تعبير ACOT12 وانخفاض SLC16A1 في نماذج NAFLD، وهو ما يتوافق مع التنبؤات المعلوماتية الحيوية.
أكدت الصبغات النسيجية تراكم الدهون وتورمًا طفيفًا في الخلايا الكبدية في هذه النماذج.
المرشحات العلاجية:
حدد التنبؤ الدوائي مركبين طبيعيين بحريين: [(1Z,13Z)-1,14-diisothiocyanatotetradeca-1,13-diene] و [(1Z,12Z)-1,13-diisothiocyanatotrideca-1,12-diene].
أشارت عمليات الإرساء الجزيئي ومحاكاة MD إلى أن هذه المركبات تتفاعل مع ACOT12 وSLC16A1 على التوالي، مع طاقات ارتباط مواتية وتكوين معقدات مستقرة.
الأهمية والادعاءات يفترض المؤلفون أن ACOT12 وSLC16A1 يعملان كسمات جزيئية رئيسية مرتبطة باستقلاب البيروفات في NAFLD. تقترح الدراسة نموذجًا تشارك فيه هذه الجينات في تقدم المرض من خلال:
توازن البروتيزوم: ربط الإجهاد الأيضي بدوران البروتين وإجهاد الشبكة الإندوبلازمية.
إعادة تشكيل البيئة المناعية الدقيقة: إظهار ارتباطات متضادة مع ارتشاح الخلايا المناعية.
تنشيط HSC: العمل كمنظمات تعتمد على السياق في الانتقال من الخلايا النجمية الكبدية الساكنة إلى النشطة، مما يؤثر على إعادة التشكيل التليفي.
تخلص الورقة إلى أن هذه النتائج توفر أدلة ميكانيكية تربط إعادة البرمجة الأيضية بالتليف، وتسلط الضوء على ACOT12 وSLC16A1 كأهداف جزيئية محتملة. ويُقدم المؤلفون المركبات الطبيعية البحرية المحددة كمرشحات أولية للاستكشاف العلاجي. ويذكر المؤلفون صراحة أنه على الرغم من أن البيانات الحاسوبية والتجريبية الأولية تدعم هذه النتائج، إلا أن الوظائف الدقيقة لهذه الجينات في الشبكة الأوسع للاستقلاب والمناعة والتليف، وفعالية المركبات المتوقعة، تتطلب مزيدًا من التحقق من خلال تجارب وظيفية مفصلة في المختبر (in vitro) وفي الكائن الحي (in vivo).