SpaCoDA: fast compositional deconvolution maps recurrent cellular niches across hundreds of spatial transcriptomics sections
يُعد SpaCoDA إطار عمل حوسبيًا سريعًا وفعالًا من حيث استهلاك الذاكرة، يقوم بفك التشفين بدقة لبيانات النسخ المكاني القائمة على البقع لرسم خرائط للمنافذ الخلوية المتكررة عبر مجموعات كبيرة، مما يتيح اكتشاف البيئات الدقيقة للأنسجة ذات الصلة سريريًا في سرطان الثدي والورم الأرومي الدبقي.
الأنسجة ليست مجرد أكياس من الخلايا؛ بل هي مجتمعات منظمة تعيش فيها أنواع مختلفة من الخلايا معاً في أحياء محددة، تماماً مثل الأحياء المتميزة في مدينة صاخبة. في العضو السليم، تعمل هذه الأحياء بسلاسة، ولكن في أمراض مثل السرطان، يصبح المخطط فوضوًا، حيث تكتظ الخلايا العدوانية وتزيح الخلايا الطبيعية، أو تفشل الخلايا المناعية في الوصول إلى الأماكن الصحيحة. ولفهم كيفية عمل هذه الأنسجة وأسباب فشلها، يستخدم العلماء تقنية تسمى "علم النسخ المكاني" (spatial transcriptomics). تعمل هذه الطريقة كخريطة عالية الدقة، حيث تقرأ النشاط الجيني للخلايا مع تتبع موقعها بدقة. ومع ذلك، يوجد عائق رئيسي: فالعديد من أدوات التخطيط الأكثر قوة تلتقط رقعة صغيرة من النسيج تحتوي على خليط من عدة خلايا مختلفة ممتزجة معاً. ولرؤية الصورة الحقيقية، يجب على الباحثين فصل هذا الخليط رياضياً لمعرفة أنواع الخلايا الموجودة بالضبط وأعدادها، وهي عملية تُعرف باسم "تفكيك المزيج" (deconvolution). وحتى الآن، كان القيام بذلك لمجموعات كبيرة من المرضى مهمة بطيئة وتستهلك الذاكرة، وغالباً ما كانت تتطلب أياماً من وقت الكمبيوتر، مما حد من عدد الأشخاص الذين يمكن دراستهم في آن واحد.
قدم فريق من الباحثين الآن أداة جديدة تسمى "SpaCoDA" تحل هذه العقبة، مما يسمح للعلماء برسم خرائط للأحياء الخلوية عبر مئات المقاطع النسيجية في جزء ضئيل من الوقت الذي كان مطلوباً سابقاً. قام الباحثون ببناء نظام يقوم أولاً بتقدير مزيج الخلايا في كل رقعة صغيرة من النسيج عن طريق مقارنة البيانات الجينية بمكتبة مرجعية لأنواع الخلايا المعروفة. وبدلاً من الاعتماد على نماذج إحصائية بطيئة ومعقدة، تستخدم طريقتهم نهجاً رياضياً أسرع يعدل الفروق بين المكتبة المرجعية والنسيج الفعلي الذي تتم دراسته. بعد ذلك، تقوم بتجميع هذه الرقع في أحياء محلية بناءً على تركيبها الجيني وموقعها الفيزيائي. وأخيراً، تجمع هذه الأحياء المحلية من جميع المرضى المختلفين في مجموعة واحدة موحدة من "الأحياء الميتاوية" (meta-neighborhoods). وهذا يسم يسمح للباحثين برؤية أي البيئات الخلوية المحددة تتكرر لدى الكثير من الناس وكيف ترتبط تلك البيئات بنتائج المرض.
اختبر الباحثون هذا النهج الجديد الجديد بصرامة قبل تطبيقه على بيانات طبية حقيقية. فقد أنشأوا آلاف خرائط الأنسجة المحاكات حيث كان تكوين الخلايا الحقيقي معروفاً، ووضعوا أداة جديدة في مواجهة ثماني طرق أخرى قائمة. في هذه المحاكاة، أدت الأداة الجديدة أداءً بدقة أفضل الطرق الموجودة، حيث جاءت في المرتبة الثانية فقط بعد المنافس الأكثر دقة في كل سيناريو، لكنها فعلت ذلك أسرع بمئات المرات. وفي إحدى المقارنات، أنهت مهمة استغرقت حوالي أحد عشر ثانية في حين استغرقت الطريقة التالية الأفضل أكثر من أربعين دقيقة. كما أثبتت الأداة دقة عالية عند اختبارها على بيانات نسيج حقيقية من أجهزة التصوير التي يمكنها رؤية الخلايا الفردية، حيث حددت أنواع الخلايا بدقة حتى عندما كانت البيانات المرجعية من تقنية مختلفة. وتعني هذه السرعة والدقة أن الباحثين يمكنهم الآن تحليل مجموعات ضخمة من عينات الأنسجة، وهو أمر كان غير عملي سابقاً.
باستخدام هذه القدرة الجديدة، طبق الفريق "SpaCoDA" على مجموعة ضخمة من عينات سرطان الثدي وأورام الدماغ، حيث حلل أكثر من مليون نقطة نسيجية متميزة من 254 مقطعاً لسرطان الثدي و181 مقطعاً للورم الأرومي الدبقي (glioblastoma). ومن هذا الكم الهائل من البيانات، حددت الأداة 29 حياً خلوياً متكرراً في سرطان الثدي و31 في الورم الأرومي الدبقي. لم تكن هذه الأحياء عشوائية؛ بل مثلت بيئات بيولوجية متميزة، مثل المناطق الغنية بالخلايا المناعية، أو المناطق التي تهيمن عليها الخلايا الورمية العدوانية، أو المناطق التي تتفاعل فيها الأوعية الدموية والأنسجة الداعمة. ثم بحث الباحثون فيما إذا كان وجود هذه الأحياء المحددة يتنبأ بكيفية تدهور حالة المرضى. في سرطان الثدي، وجدوا أن الحي المليء بخلايا B والخلايا البلازمية كان مرتبطاً بفترات بقاء أطول واستجابة أفضل للعلاج. وعلى العكس من ذلك، في الورم الأرومي الدبقي، ارتبط وجود حي معين يتميز بخلايا تعاني من نقص الأكسجين وتنقسم بسرعة، وآخر مليء بالخلايا المناعية التي تعرض المستضدات بالقرب من الأوعية الدموية، بفترات بقاء أقصر.
تشير هذه النتائج إلى أن الترتيب المحدد للخلايا داخل الورم هو مؤشر قوي على سلوكه. لا تدعي الدراسة علاج هذه الأمراض، لكنها توفر إطاراً عملياً لتحديد البيئات الخلوية الأكثر أهمية. ومن خلال جعل تحليل مئات المرضى ممكناً بسرعة، تفتح الأداة الباب لاكتشاف أنماط جديدة يمكن أن تساعد الأطباء في تصنيف الأمراض بدقة أكبر أو تحديد أهداف جديدة للعلاج. ويؤكد الباحثون أنه بينما الروابط بين هذه الأحياء ونتائج المرضى قوية، إلا أنها استكشافية وتحتاج إلى مزيد من التحقق. ومع ذلك، فإن القدرة على رسم خرائط لهذه المناظر الطبيعية الخلوية عبر مجموعات كبيرة من الناس تمثل خطوة هائلة للأمام في فهم الجغرافيا المعقدة للمرض.
ملخص تقني لـ SpaCoDA
بيان المشكلة تقيس تقنيات النسخ المكاني القائمة على البقع (مثل 10x Visium) الحمض النووي الريبوزي (RNA) المجمع من خلايا متعددة داخل كل بقعة، مما يستلزم عملية "إزالة الخلط" (deconvolution) لأنواع الخلايا لرسم خرائط البيئات المجهرية للأنسجة. وبينما توجد طرق عديدة لإزالة الخلط، فإن النهج الأكثر دقة غالبًا ما تكون مكلفة حاسوبيًا وتستهلك ذاكرة مكثفة، مما يجعلها غير عملية للمجموعات الضخمة التي تضم مئات المقاطع وملايين البقع. علاوة على ذلك، هناك نقص في الأطر الراسخة لجعل "البيئات المجهرية" (niches) للأنسجة (وهي جيران خلويون متكررون) قابلة للمقارنة عبر مئات المقاطع ومجموعات البيانات المختلفة، وهو مطلب ضروري للأسئلة ذات الدوافع السريرية المتعلقة بتصنيف الأمراض والاستجابة للعلاج. غالبًا ما تستهدف أدوات تحديد البيئات المجهرية الحالية بيانات ذات دقة الخلية الواحدة، مما يترك فجوة في تحليل تركيبات البقع التي تمت إزالة الخلط منها على نطاق المجموعات الكبيرة.
المنهجية يقدم المؤلفون SpaCoDA، وهو إطار عمل يتكون من مكونين رئيسيين: مُقدّر سريع لإزالة الخلط يعتمد على الأوزان المتباينة، ومسار تحليل بيئات مجهرية على مستويين.
مُقدّر إزالة الخلط:
إعداد المرجع: يتم تحجيم بيانات المرجع أحادية الخلية المُعنونة لتصل إلى حجم مكتبة موحد، وتُجمع حسب نوع الخلية، وتُطبع (normalize) لتلخص الواحد الصحيح لإنشاء مصفوفة التوقيع (signature matrix). يتم ترتيب الجينات العلامة (marker genes) حسب النوعية والتعبير.
المربعات الصغرى غير السالبة الموزونة بالتباين (Variance-Weighted NNLS): يتم مطابقة كل بقعة مع توقيع المرجع (المقتصر على الجينات العلامة) باستخدام طريقة المربعات الصغرى غير السالبة الموزونة بالتباين (NNLS). تُشتق الأوزان من علاقة متوسط-تباين ذات توزيع نيغ-بينوميال (negative-binomial).
تصحيح المنصة التكراري: بين التكرارات، يتم إعادة تقدير عامل المنصة على مستوى الجين من خلال مقارنة ملف التعبير المرصود المتوسط مع الملف المتوقع من المتوسط التركيبي. يقوم هذا العامل بإعادة تحجيم التوقيع لمراعاة تأثيرات المنصة المحددة (مثل الاختلافات في كفاءة التقاط الجينات بين scRNA-seq والبيانات المكانية).
التحويل إلى كسور: تُحول حصص الـ RNA المُطابقة إلى كسور خلوية عن طريق القسمة على عوامل حجم الخلية (متوسط محتوى الـ RNA لكل نوع) وجزء الـ RNA الذي تم التقاطه بواسطة الجينات العلامة المختارة.
تحليل البيئات المجهرية (Niche Analysis) على مستويين:
البيئات المجهرية المحلية: داخل كل مقطع، يتم تجميع البقع بناءً على مقياس مسافة مدمج يجمع بين مسافة آيتشيسون (Aitchison distance) (التشابه التركيبي) والمسافة الإقليدية المكانية. يتم دمج هذه المصفوفة عبر تحليل الإحداثيات الرئيسية (PCoA) وتجميعها باستخدام رسم بياني للجار الأقرب المشترك (SNN) وخوارزمية لوفين (Louvain). يعمل متوسط التركيب لكل عنقود محلي كـ "نموذج أولي" (prototype).
البيئات المجهرية العليا (Meta-Niches): يتم تطبيع النماذج الأولية من جميع المقاطع، وتقليل أبعادها عبر PCA، ودمجها عبر المجموعة بأكملها باستخدام تقنية Harmony لتصحيح تأثيرات "مقطع الأصل". يتم تجميع هذه النماذج الأولية المدمجة لتحديد "البيئات المجهرية العليا"، وهي حالات تركيبية متكررة قابلة للمقارنة عبر المجموعة بأكملها.
المساهمات الرئيسية
الكفاءة الخوارزمية: يستبدل SpaCoDA احتمالات العد (مثل توزيع Poisson في RCTD) بطريقة NNLS الموزونة بالتباين، مما يسمح بحلول دقيقة لجميع البقع في وقت واحد. يؤدي هذا إلى تحسينات في السرعة بمقدار عدة مراتب عشرية وتقليل استهلاك الذاكرة مقارنة بالطرق الرائدة.
تكامل البيئات المجهرية على نطاق المجموعة: يقدم البحث سير عمل محدد لتحديد ودمج "البيئات المجهرية العليا" من بيانات البقع التي تمت إزالة الخلط منها، مما يتيح مقارنة البيئات المجهرية للأنسجة عبر مئات المقاطع ودراسات متعددة.
المتانة تجاه عدم التطابق: تتضمن الطريقة تصحيحًا تكراريًا للمنصة للتعامل مع التباينات بين المرجع والبيانات المستهدفة، بما في ذلك تأثيرات المنصة على مستوى الجين وضوضاء أخذ العينات.
النتائج
الاختبار المعياري على البيانات المحاكات: عبر 123 تكوينًا محاكى تحت ثلاثة سيناريوهات لعدم تطابق المرجع (الضوضاء المستقلة، تأثيرات المنصة، وضوضاء أخذ العينات)، صنف SpaCoDA باستمرار ضمن أفضل طريقتين في جميع مقاييس التقييم الأربعة (ارتباط بيرسون/سبيرمان، RMSE، JSD). جاء في المرتبة الثانية فقط بعد RCTD في متوسط التصنيف ولكنه كان أسرع بشكل ملحوظ.
الاختبارات المستمدة من التصوير: في ست مجموعات بيانات من Xenium/CosMx ذات تركيبات معروفة (اختبارات لنفس المقطع وعبر المنصات)، حقق SpaCoDA أفضل متوسط تصنيف عبر جميع المقاييس. كما كان أسرع من الطرق المنافسة بمعدل يتراوح من 4 إلى حوالي 1,80{0} ضعفًا (متوسط 10.8 ثانية لكل مجموعة بيانات مقابل 2,517 ثانية لـ RCTD) واستخدم أقل ذروة ذاكرة (~1.69 جيجابايت).
أطلس سرطان الثدي: عند تطبيقه على 254 مقطعًا لسرطان الثدي (~594 ألف بقعة)، حدد SpaCoDA 29 بيئة مجهرية عليا. ارتبطت البيئة المجهرية العليا للخلايا البائية/الخلايا البلازمية (Meta-niche 23) بطول البقاء الإجمالي في مجموعات TCGA-BRCA و METABRIC، وكانت أعلى بشكل ملحوظ في الأورام التي حققت استجابة مرضية كاملة (pCR) في تجربة I-SPY2.
أطلس الورم الأرومي الدبقي (Glioblastoma): عند تطبيقه على 181 مقطعًا للورم الأرومي الدبقي (~585 ألف بقعة)، حدد SpaCoDA 31 بيئة مجهرية عليا. ارتبطت بصمات البيئة المجهرية العليا للخلايا المقدمة للمستضد حول الأوعية (Meta-niche 1) والبيئة المجهرية التكاثرية ونقص الأكسجين (Meta-niche 29) بقصر مدة البقاء الإجمالي في مجموعات TCGA-GBM وثلاث مجموعات CGGA.
الأهمية والادعاءات يدعي المؤلفون أن SpaCoFA يوفر توازنًا مثاليًا بين الدقة التنافسية والتكلفة الحاسوبية، مما يجعل تقدير التركيب وتحليل البيئات المجهرية أمرًا عمليًا على نطاق الأطالس المكانية الحالية (مئات المقاطع). يؤكد البحث أن SpaCoDA يوفر إطار عمل قابل للتكرار لتحديد البيئات الخلوية المتكررة عبر المرضى، مما يمكن أن يكون نقطة انطلاق لتصنيف الأمراض وتحديد الأهداف العلاجية الجديدة. ويذكر المؤلفون صراحة أن الارتباطات السريرية الموجودة (مثل ارتباطات البقاء على قيد الحياة) هي استكشافية وتولد فرضيات للتحقق المستقبلي. كما يقرون بالقيود، مثل عدم القدرة على تصحيح أنواع الخلايا المفقودة من المرجع أو الخطأ في التسمية، ويشيرون إلى أن البيئات المجهرية العليا تمثل حالات تركيبية تقديرية متكررة وليست دليلًا مباشرًا على الترتيبات الخلوية المشتركة.