معايرة نماذج PDE على مستوى النسيج لديناميكيات اللِّيجند باستخدام بيانات التنسخ أحادي الخلية والمكانية

كيف تتواصل الخلايا عبر الجرح أثناء الشفاء

عندما يتعرض جلدنا لإصابة، يجب على عدد هائل من الخلايا التنسيق لإصلاح الضرر دون ترك نسيج ندبي مفرط. تقوم الخلايا بذلك عن طريق إرسال واستشعار رسائل كيميائية تسمى اللِّيجندات، التي تنتشر عبر النسيج وتثير استجابات في الخلايا المجاورة. تُظهر هذه الدراسة كيف يمكن تحويل القياسات الجينية الحديثة على مستوى الخلية المفردة والمكانية إلى طريقة قوية جديدة لمعرفة معايرة النماذج الرياضية لتلك الرسائل الكيميائية، مما يساعد الباحثين على فهم أفضل للشفاء الطبيعي والأمراض المرتبطة بالتندب.

لماذا تهم الرسائل الكيميائية داخل الأنسجة

نادراً ما تتصرف الخلايا في أنسجتنا بشكل منعزل. فهي تتبادل باستمرار إشارات تخبرها متى تنقسم أو تتحرك أو تهدأ. العديد من هذه الإشارات هي مواد كيميائية قصيرة المدى تفرزها خلية وتلتقطها مستقبلات على خلية أخرى. هذه الرسائل مركزية في التئام الجروح والالتهاب وأمراض مثل التليف والسرطان. لطالما استخدم الباحثون معادلات رياضية تعرف بنماذج التفاعل–الانتشار لوصف كيفية إنتاج هذه المواد الكيميائية وانتشارها وتحللها. ومع ذلك، ظلت مشكلة قديمة تتمثل في كيفية اختيار قيم معلمية واقعية — مثل معدلات الإنتاج وقوة الارتباط — بحيث تعكس النماذج ما يحدث في الأنسجة الحية فعلًا بدلاً من ظروف المختبر المثالية.

نظرة جديدة على الأنسجة الحية

مؤخرًا، حولت تقنيتان تجريبيتان الطريقة التي يمكن للعلماء بها مراقبة الأنسجة. يقيس التسلسل الجيني للـRNA على مستوى الخلية المفردة أي الجينات النشطة في آلاف الخلايا الفردية، مكشوفًا هوياتها وقدراتها الإشارية. بينما يقيس التنسخ المكاني نشاط الجينات عبر العديد من البقع الصغيرة المرتبة كشبكة على شريحة نسيجية، محافظًا على مواضع الخلايا المختلفة. في هذا العمل، يجمع المؤلفون كلا نوعي البيانات من عينات جلد بشري مأخوذة بعد 30 يومًا من الجرح، عندما يكون النسيج في مرحلة «إعادة التشكيل» من الشفاء. يركزون على ثلاثة أشكال من جزيء إشاري رئيسي، عامل النمو المحول بيتا (TGFβ)، المعروف بتأثيره على تشكيل الندبات وإصلاح النسيج.

خط سير خطوة بخطوة من الجينات إلى معلمات النموذج

يبني الباحثون خط سير حسابي يربط هذه القياسات الجينية الغنية بنموذج على مستوى النسيج لحركة وتأثير TGFβ. أولًا، يستخدمون أدوات معلوماتية حيوية معروفة لتحديد الأنواع الخلوية الرئيسية والأنماط الفرعية — مثل ثلاث فئات مميزة من الخلايا الليفية، والبلاعم، والخلايا البطانية — ولتقدير عدد كل نوع في كل بقعة من شبكة التنسخ المكاني. بعد ذلك يستنتجون مدى قوة احتمالية تواصل كل نوع خلية مع الآخر عبر TGFβ، استنادًا إلى نشاط جينات اللِّيجندات في الخلايا المرسلة والمستقبلات في الخلايا المستقبِلة. تعمل هذه «قوى التفاعل» المستنتجة كهدف تجريبي يجب على النموذج الرياضي مضاهاةه.

تعليم النموذج ليطابق النسيج

بعد ذلك، يمثل الفريق انتشار TGFβ عبر النسيج بمجموعة من المعادلات التفاضلية الجزئية تُحل على شبكة تعكس تخطيط التنسخ المكاني. لأي مجموعة مختارة من قيم المعلمات — مثل سرعة انتشار TGFβ، وسرعة تحلله، وكمية ما تنتجه أو تمتصه كل فئة من الخلايا — يتنبأ النموذج بحقول التركيز للأشكال الثلاثة من TGFβ وبقوى التفاعل الناتجة بين أنواع الخلايا. لضبط تلك المعلمات، يستخدم المؤلفون استراتيجية معايرة من ثلاث مراحل. يطبقون أولًا طرقًا بايزية تقريبية لمسح نطاق واسع من المعلمات المعقولة مستندة إلى تجارب سابقة، مقلصين تدريجيًا إلى تراكيبات تتوافق جيدًا مع قوى التفاعل المستنتجة من البيانات الجينية. ثم يمررون هذه المرشحات الواعدة إلى طريقة تحسين تعتمد على التدرج تقوم بتحسين المعلمات لزيادة الاتفاق إلى أقصى حد. يحقق هذا النهج الهجين تطابقًا خطيًا شبه كامل (معامل ارتباط 0.99) بين قوى التفاعل المتوقعة بواسطة النموذج والمستنتجة من البيانات.

ما يكشفه النموذج المعاير عن عملية الشفاء

مع وجود نموذج مطابق جيدًا، يستكشف الباحثون ما تعنيه المعلمات المستنتجة بيولوجيًا. يؤكدون أن الأنماط الثلاثة للخلايا الليفية المحددة من بيانات التعبير الجيني مختلفة وظيفيًا، حيث يظهر نمط واحد ملفًا التهابيًا قويًا وآخر سمات خلايا ليفية-عضلية تساعد على تقلص الجرح. تشير معدلات الإنتاج المعايرة إلى أن الخلايا الليفية ذات الطابع المؤيد للالتهاب تنتج أقل قدر من TGFβ خلال مرحلة إعادة التشكيل، متناسقًا مع تراجع الالتهاب بحلول اليوم الثلاثين. ومن المثير للدهشة أن نمطًا فرعيًا معروفًا بالخلايا الليفية الحليمية يبدو أنه ينتج كميات نسبية مرتفعة من TGFβ3، وهو شكل غالبًا ما يرتبط بحدوث ندبات أقل ونتائج أكثر تجددًا. تتوافق خرائط تراكيز TGFβ المتوقعة من النموذج مع المناطق التي تتجمع فيها الخلايا الليفية والبلاعم، مما يشير إلى بقع ساخنة للتواصل النشط قد توجه كيفية تشكيل الندبة.

إطار مرن للأنسجة والأمراض المستقبلية

بعيدًا عن حالة التئام الجروح المحددة هذه، يتمثل الإسهام الرئيسي للدراسة في خط سير عام لتحويل بيانات جينية حديثة إلى نماذج نسيجية معايرة جيدًا. من خلال الجمع بين التنسخ أحادي الخلية والمكاني مع أدوات إحصائية وتحسين متقدمة، يمكن للباحثين تقدير معدل حركة اللِّيجندات، وقوة ارتباطها، وأي الخلايا تسهم أكثر في الإشارة بصورة أكثر صرامة. يمكن تكييف هذا الإطار لأنسجة أخرى، وجزيئات إشارية أخرى، وحالات مرضية مثل التندب الليفي أو السرطان، مما يساعد في نهاية المطاف على ترجمة لقطات جزيئية معقدة إلى نماذج تنبؤية لسلوك الأنسجة كاملة.

الاستشهاد: Daher, A., Trucu, D. & Eftimie, R. Calibrating tissue level PDE models of ligand dynamics using single cell and spatial transcriptomics data. npj Syst Biol Appl 12, 44 (2026). https://doi.org/10.1038/s41540-026-00657-8

الكلمات المفتاحية: التئام الجروح, إشارة خلوية, التنسخ المكاني, عامل النمو المحول بيتا, النمذجة الرياضية