Quantifying cell shape and density fluctuations in epithelial tissue in vivo

ملخص تقني: قياس تقلبات شكل الخلايا وكثافتها في الأنسجة الطلائية داخل الكائن الحي

بيان المشكلة
تعد الأنسجة الطلائية أنظمة غير متجانسة، حيث تلعب الانحرافات المحلية عن متوسط شكل وحجم الخلايا (التقلبات) دوراً حاسماً في العمليات البيولوجية مثل التطور والالتئام الجرحي. وبينما تقوم النماذج الحالية إما بمحاكاة الخلايا الفردية (مثل نماذج Vertex أو Cellular Potts) أو باستخدام معادلات متصلة تتجاهل المعلومات على مقياس الخلية، تبرز الحاجة إلى إطار عمل يربط بين هذه المقاييس. وتحديداً، يهدف المؤلفون إلى قياس كيفية تطور التقلبات المحلية في شكل الخلايا وكثافتها بمرور الوقت، وكيف يمكن لهذه التقلبات أن تكشف عن الخصائص الميكانيكية العينية، مثل "صلابة" النسيج الفعالة وإعادة التشكيل النشط للخلايا، داخل جناح ذبابة الدروسوفيلا (Drosophila) في مرحلة العذراء.

المنهجية
يجمع البحث بين التصوير عالي الدقة وبين إطار نظري اختزالي:

1. الحصول على البيانات التجريبية:

   • النظام: يركز البحث على نسيج ظهارة جناح ذبابة الدروسوفيلا في مرحلة العذراء من الساعة 18 إلى 21 بعد تشكل الغمد (APF)، وهي فترة تتميز بانكماش المفصل، وتمدد نصل الجناح، واستطالة الخلايا على طول المحور القريب-البعيد (P/D).
   • التصوير: استُخدم المجهر البؤري (Confocal microscopy) على عذارى تعبر عن بروتين E-cadherin-GFP (لحدود الخلايا) وHistone2-RFP (للنوى). تم التقاط الصور كل دقيقتين على مدار نافذة زمنية مدتها 3 ساعات.
   • التحليل: استخدمت خوارزمية التعلم العميق للكشف عن حدود الخلايا، والتي تم مطابقتها مع مضلعات لاستخراج مساحة الخلية وشكلها. تم ترميز أشكال الخلايا باستخدام موتر $q$ من الرتبة الثانية، والذي يمثل الانحرافات عن الشكل الدائري (الاستطالة والتوجه). تم جمع أكثر من 1,000,000 قياس على مستوى الخلية.

2. الإطار النظري:

   • النموذج المجهري: تُنمذج الخلايا كأشكال بيضاوية لينة متفاعلة ذات أحجام ونسب استطالة متفاوتة. يخضع النظام لمعادلات تفاضلية عشوائية (SDEs) لمواضع الخلايا ($r_i$) وموترات الشكل ($q_i$). تتضمن هذه المعادلات حدوداً حتمية (المحاذاة النشطة، والتفاعلات الزوجية عبر جهد $U$) وحدود ضوضاء عشوائية (تمثل التقلبات الميكانيكية والكيميائية).
   • التجميع الخشن (Coarse-Graining): يتم إجراء عملية التجميع الخشن للمعادلات المجهرية لاشتقاق معادلات تفاضلية عشوائية متصلة للحقل المحلي للكثافة ($\hat{\rho}$) والحقل المحلي للشكل ($\hat{Q}$).
   • الخطية (Linearization): تمت خطية المعادلات حول القيم المتوسطة، مما فصل الديناميكيات إلى تقلبات "مُلدنة" (تعتمد على الزمن) وانحرافات "مُثبتة" (غير معتمدة على الزمن هيكلياً).
   • دوال الارتباط: اشتق المؤلفون تعبيرات تحليلية لدوال الارتباط الزمكاني للتقلبات في الكثافة ($C_{\rho\rho}$) والشكل ($C_{qq}$). تعتمد هذه الدوال على معاملات تمثل الانتشار، وتغيرات الشكل النشطة، وصلابة النسيج الفعالة.

المساهمات والنتائج الرئيسية
تتمثل المساهمة الأساسية في إطار نظري نجح في مطابقة دوال الارتباط التجريبية لاستخراج معاملات ميكانيكية كمية من الأنسجة الحية.

1. استخراج المعاملات من الأنسجة البرية (Wild-Type):
   من خلال مطابقة دوال الارتباط النظرية مع البيانات التجريبية، حدد المؤلفون عدة معاملات رئيسية (الجدول I):

   • الصلابة الفعالة ($L$): $10.6 \pm 0.5 \, \mu\text{m}^4 \text{min}^{-1}$. يرمز هذا المعامل إلى مقاومة النسيج للتشوه.
   • تغير الشكل النشط ($A$): $(2.83 \pm 0.48) \times 10^{-3} \, \text{min}^{-1}$. يمثل معدل الاسترخاء النشط لأشكال الخلايا نحو حالة مفضلة.
   • الانتشار والضوضاء: وُجد أن ثابت الانتشار للكثافة ($D$) هو $7.4 \, \mu\text{m}^2 \text{min}^{-1}$، وهو أكبر بكثير من انتشار مركز الخلية ($D_r$)، مما يشير إلى أن الحركة الجماعية هي السائدة.
   • الانحرافات الهيكلية: حددت الدراسة مكون كثافة هيكلي "مُثبت" بطول اضمحلال ($\lambda$) قدره $23 \, \mu\text{m}$، ويُعزى ذلك إلى عمليات نمو سابقة (مثل تكوين العروق).

2. التحقق عبر تجارب الاضطراب:
   تم اختبار الإطار على أنسجة مضطربة للتحقق من قدرته على التمييز بين الحالات الميكانيكية:

   • جرح الليزر: بعد إغلاق الجرح، أظهر النسيج بالقرب من موقع الجرح انخفاضاً بنسبة 22% في الصلابة الفعالة ($L$)، مما يشير إلى أن النسيج يصبح أكثر قابلية للتشوه لتسهيل الالتئام. وفي المقابل، زاد معامل تغير الشكل النشط ($A$) بمقدار 1.6 مرة بالقرب من الجرح، مما يشير إلى تسارع إعادة التشكيل النشط. أما بعيداً عن الجرح، فقد ظل $A$ مشابهاً للأنسجة السليمة، مما يشير إلى نشاط محصور مكانياً.
   • الاضطراب الجيني (إشارات JNK): في الأنسجة التي تم فيها تثبيط إشارات JNK (عبر التعبير عن بروتين bsk سالب هيمنة)، تأخر التئام الجروح. كشف النموذج أنه بينما انخفضت الصلابة ($L$) (بشكل مشابه للنسيج المجروح في النوع البري)، فإن معامل تغير الشكل النشط ($A$) فشل في الارتفاع. هذا النقص في إعادة التشكيل النشط يرتبط ببطء الالتئام الملحوظ، مما يثبت قدرة النموذج على الفصل بين الصلابة السلبية والدفع البيولوجي النشط.

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

• يمكن لدوال الارتباط الزمكاني أن تعمل كمسبار غير جراحي لقياس "صلابة" النسيج الفعالة ومعدلات إعادة التشكيل النشط.
• يمكن للإطار التمييز بين التغيرات الميكانيكية السلبية (تغير الصلابة) والاستجابات البيولوجية النشطة (تغير معدلات إعادة التشكيل) استجابة للاضطرابات مثل الجرح أو التثبيط الجيني.
• يوفر النموذج جسراً بين التفاعلات المجهرية بين الخلايا والسلوك العيني للنسيج، مما يقدم طريقة لتتبع كيفية تطور خصائص النسيج أثناء التخلق (morphogenesis) أو الإصلاح دون الحاجة إلى تجارب ارتداد الليزر الجراحية.

يشير المؤلفون إلى أنه على الرغم من أن النموذج يوفر "ملاءمة جيدة بشكل معقول" ويلتقط الفيزياء الأساسية، إلا أنه يعتمد على افتراضات تبسيطية (مثل تجانس الصلابة، وتجاهل كسر التماثل النيماتي). وهم يقترحون أن التحسينات المستقبلية قد تتضمن معاملات صلابة متباينة الخواص (anisotropic) لتحقيق مطابقة أفضل مع البيانات، لكن الإطار الحالي كافٍ لاستخراج معاملات فعالة ذات معنى من الأنظمة المعقدة داخل الكائن الحي.