DOI: https://doi.org/10.1038/s41698-023-00488-4
PMID: https://pubmed.ncbi.nlm.nih.gov/38200223
تاريخ النشر: 2024-01-10
المؤلف: Alberto Valdeolivas وآخرون
الموضوع الرئيسي: علم النسخ الجيني أحادي الخلية والمكاني
نظرة عامة
تناقش الورقة البحثية قيود الأنماط الجزيئية التوافقية (CMS) لسرطان القولون والمستقيم (CRC) بسبب التباين داخل الورم، مما يعيق فائدتها السريرية. لمعالجة هذه المشكلة، استخدم المؤلفون تقنيات النسخ الجزيئي المكاني (ST) جنبًا إلى جنب مع تسلسل RNA على مستوى الخلية الواحدة (scRNA-seq) لتحليل التركيب الخلوي والجزيئي المكاني لـ CRC. سمح هذا النهج برسم خرائط مفصلة للتباين داخل الورم ضمن CMS والبيئة الدقيقة المحيطة بها.
من المهم أن الدراسة حددت أحداث التواصل الخلوي الرئيسية عند واجهة الورم والستروما في سرطانات CMS2، كاشفة عن إشارات تعيق نمو الورم وأخرى تنشطه. تشمل النتائج الملحوظة التنظيم المحتمل لنشاط النسخ لـ ETV4 بواسطة الديكورين (DCN) والتفاعل بين الليغاندات والمستقبلات PLAU و PLAUR. تؤكد النتائج على فائدة ST في توضيح التباين الجزيئي لـ CRC، مما يسهل التقدم في استراتيجيات العلاج الشخصي.
مقدمة
يظل سرطان القولون والمستقيم (CRC) تحديًا صحيًا عالميًا كبيرًا، مع أكثر من 1.85 مليون حالة و850,000 وفاة سنويًا. على الرغم من أن العلاجات الشخصية ساهمت في انخفاض معدلات الوفيات، إلا أن التباين الداخلي بين المرضى وداخل الورم يعقد العلاج الفعال. يتضح هذا التباين على المستويات الجينومية، والإبيجينومية، والنسخية، وكذلك في البيئة الدقيقة للورم (TME). تم اعتماد تصنيف ائتلاف تصنيف سرطان القولون والمستقيم إلى أربعة أنماط جزيئية توافقية (CMS) على نطاق واسع؛ ومع ذلك، فإن فائدته السريرية محدودة بسبب الاعتماد على التسلسل الكمي، الذي يفشل في التقاط الأنماط الظاهرية المختلطة لـ CMS والتركيب الخلوي المعقد للأورام بدقة.
بدأت التطورات الأخيرة في تسلسل RNA على مستوى الخلية الواحدة (scRNA-seq) في توضيح خصائص CMS على المستوى الخلوي، كاشفة عن وجود عدة CMS داخل المرضى الفرديين. ومع ذلك، لا تزال التوزيعات المكانية والتفاعلات لهذه CMS مع بيئاتها الدقيقة غير مفهومة بشكل كافٍ. تقدم النسخ الجزيئي المكاني (ST) نهجًا واعدًا للتغلب على هذه القيود من خلال تمكين قياس التعبير الجيني عبر المساحات النسيجية، مما يدمج البيانات الشكلية والنسخية. في هذه الدراسة، استخدمنا ST لتحليل 14 عينة من سبعة مرضى مصابين بـ CRC، بهدف تعزيز فهمنا للخصائص المكانية والتباين لـ CMS. تشمل نتائجنا رسم خرائط التركيب الخلوي، وارتباط الميزات الجزيئية والشكلية المتميزة مع CMS المختلفة، ورؤى حول التفاعلات بين الخلايا في سرطانات CMS2، مما يبرز إمكانيات ST في إبلاغ استراتيجيات العلاج الشخصية لـ CRC.
الطرق
تحدد قسم الطرق تصميم التجارب والتقنيات التحليلية المستخدمة في الدراسة. استخدم الباحثون نهجًا كميًا، حيث قاموا بإجراء تحليلات إحصائية لتقييم البيانات التي تم جمعها من تجارب مختلفة. تضمنت المنهجيات المحددة تجارب مختبرية محكومة، حيث تم التلاعب بالمتغيرات بشكل منهجي لملاحظة آثارها على النتائج المعنية.
شمل جمع البيانات استخدام أدوات وبروتوكولات موحدة لضمان الموثوقية والصلاحية. تم إجراء التحليل باستخدام أدوات برمجية سهلت تطبيق الاختبارات الإحصائية المناسبة، مثل ANOVA وتحليل الانحدار، لتحديد أهمية النتائج. يؤكد القسم على أهمية القابلية للتكرار والشفافية في عملية البحث، موضحًا الخطوات المتخذة للتخفيف من التحيزات وتعزيز قوة النتائج.
النتائج
في هذا القسم، قام المؤلفون بتقييم اتساق نتائج التفكيك عبر النسخ التقنية من خلال استخدام تصحيح الدفعة على الملفات النسخية باستخدام Harmony. بعد ذلك، استخدموا خوارزمية لوفين، المطبقة في وظيفة FindClusters لحزمة Seurat، لتجميع تمثيلات Harmony. تم اختيار مجموعة من معلمات الدقة الكبيرة (من 1 إلى 2، بزيادة 0.1) لإنشاء مجموعات دقيقة تتوافق مع المناطق التشريحية ذات التوزيعات الخلوية المماثلة عبر النسخ.
لتحديد اتساق نتائج التفكيك، قام المؤلفون بحساب متوسط عدد المعرفات الجزيئية الفريدة (UMIs) المقدرة بواسطة Cell2Location لكل نوع خلية ومجموعة. ثم طبقوا معامل ارتباط بيرسون لتقييم تشابه هذه التقديرات بين النسخ التقنية، مما يوفر مقياسًا إحصائيًا لإعادة إنتاج نتائجهم.
المناقشة
استخدم البحث تقنيات النسخ الجزيئي المكاني (ST) وتقنيات التفكيك لتحليل التركيب الخلوي والتوزيع المكاني لأنواع الخلايا في أورام سرطان القولون والمستقيم (CRC) من سبعة مرضى بخصائص متنوعة. تم معالجة ما مجموعه 20,733 نقطة Visium، كاشفة عن تباين كبير بين المرضى في نسب أنواع الخلايا، لا سيما في توزيعات الخلايا الورمية والمناعية. من الجدير بالذكر أن الخلايا الورمية شكلت حوالي 5% فقط في عينة أحد المرضى، والتي كانت في الغالب نسيجًا غير ورمي. أظهرت الدراسة قابلية عالية للتكرار عبر النسخ التقنية، حيث تجاوزت معاملات ارتباط بيرسون 0.9 لمعظم العينات. تم تأكيد النتائج من خلال تعليقات أخصائي الأمراض، مما يؤكد أن التوزيع المكاني لأنواع الخلايا يتماشى مع الأقسام التشريحية المتوقعة.
علاوة على ذلك، وصف البحث البيئة الدقيقة للورم (TME) والأنماط الجزيئية التوافقية (CMS) لـ CRC، محددًا هيمنة خلايا الورم CMS2 في عدة عينات. كشفت التحليلات عن تركيبات خلوية مناعية وستروما متميزة مرتبطة بأنماط CMS المختلفة، حيث أظهرت أورام CMS1-CMS2 المختلطة نسبًا أعلى من خلايا T وB. كما سلطت الدراسة الضوء على دور عوامل النسخ ومسارات الإشارات المحددة، مثل WNT وVEGF، في تقدم الورم. من المهم أن البحث أسس صلة بين التفاعلات الخلوية المحددة والميزات الجزيئية لـ CRC، مما يوفر رؤى حول الآليات الكامنة وراء التباين وتقدم الورم. بشكل عام، تؤكد هذه الدراسة على فائدة ST في توضيح الهيكل الخلوي المعقد وعمليات التواصل داخل أورام CRC.
DOI: https://doi.org/10.1038/s41698-023-00488-4
PMID: https://pubmed.ncbi.nlm.nih.gov/38200223
Publication Date: 2024-01-10
Author(s): Alberto Valdeolivas et al.
Primary Topic: Single-cell and spatial transcriptomics
Overview
The research paper discusses the limitations of the consensus molecular subtypes (CMS) of colorectal cancer (CRC) due to intratumoral heterogeneity, which hampers its clinical utility. To address this issue, the authors employed Spatial Transcriptomics (ST) alongside single-cell RNA sequencing (scRNA-seq) to analyze the spatial cellular and molecular composition of CRC. This approach allowed for a detailed mapping of the intratumoral heterogeneity within CMS and their surrounding microenvironment.
Significantly, the study identified key cell communication events at the tumor-stroma interface in CMS2 carcinomas, revealing both tumor growth-inhibiting and -activating signals. Notable findings include the potential regulation of ETV4 transcriptional activity by decorin (DCN) and the interaction between the PLAU and PLAUR ligands and receptors. The results underscore the utility of ST in elucidating the molecular heterogeneity of CRC, thereby facilitating advancements in personalized therapeutic strategies.
Introduction
Colorectal cancer (CRC) remains a significant global health challenge, with over 1.85 million cases and 850,000 deaths annually. Although personalized treatments have contributed to a decline in mortality rates, the inherent inter-patient and intra-tumor heterogeneity complicates effective treatment. This heterogeneity is evident at genomic, epigenomic, and transcriptomic levels, as well as in the tumor microenvironment (TME). The CRC subtyping consortium’s classification into four Consensus Molecular Subtypes (CMS) has been widely adopted; however, its clinical utility is limited by the reliance on bulk sequencing, which fails to accurately capture mixed CMS phenotypes and the intricate cellular composition of tumors.
Recent advancements in single-cell RNA sequencing (scRNA-seq) have begun to elucidate CMS characteristics at the cellular level, revealing the coexistence of multiple CMS within individual patients. Nonetheless, the spatial distribution and interactions of these CMS with their respective TMEs remain inadequately understood. Spatial transcriptomics (ST) offers a promising approach to overcome these limitations by enabling the measurement of gene expression across tissue spaces, thereby integrating morphological and transcriptomic data. In this study, we employed ST to analyze 14 samples from seven CRC patients, aiming to enhance our understanding of the spatial properties and heterogeneity of CMS. Our findings include the spatial mapping of cell type composition, the association of distinct molecular and morphological features with different CMS, and insights into intercellular interactions in CMS2 carcinomas, underscoring ST’s potential to inform personalized treatment strategies for CRC.
Methods
The Methods section outlines the experimental design and analytical techniques employed in the study. The researchers utilized a quantitative approach, employing statistical analyses to evaluate the data collected from various experiments. Specific methodologies included controlled laboratory experiments, where variables were systematically manipulated to observe their effects on the outcomes of interest.
Data collection involved the use of standardized instruments and protocols to ensure reliability and validity. The analysis was conducted using software tools that facilitated the application of appropriate statistical tests, such as ANOVA and regression analysis, to determine the significance of the findings. The section emphasizes the importance of replicability and transparency in the research process, detailing the steps taken to mitigate biases and enhance the robustness of the results.
Results
In this section, the authors assessed the consistency of deconvolution results across technical replicates by employing batch correction on transcriptomic profiles using Harmony. Following this, they utilized the Louvain algorithm, implemented in the Seurat package’s FindClusters function, to cluster the Harmony embeddings. A range of large resolution parameters (from 1 to 2, incremented by 0.1) was selected to generate fine-grained clusters that correspond to anatomical regions with similar cell type distributions across the replicates.
To quantify the consistency of the deconvolution results, the authors calculated the mean number of unique molecular identifiers (UMIs) estimated by Cell2Location for each cell type and cluster. They then applied Pearson’s correlation coefficient to evaluate the similarity of these estimates between the technical replicates, providing a statistical measure of the reproducibility of their findings.
Discussion
The research utilized spatial transcriptomics (ST) and deconvolution techniques to analyze the cellular composition and spatial distribution of cell types in colorectal cancer (CRC) tumors from seven patients with varying characteristics. A total of 20,733 Visium spots were processed, revealing significant inter-patient heterogeneity in cell type proportions, particularly in tumor and immune cell distributions. Notably, tumor cells constituted only about 5% in one patient’s sample, which was predominantly non-neoplastic tissue. The study demonstrated high reproducibility across technical replicates, with Pearson’s correlation coefficients exceeding 0.9 for most samples. The findings were corroborated by a pathologist’s annotations, confirming that the spatial distribution of cell types aligned with expected anatomical compartments.
Furthermore, the research characterized the tumor microenvironment (TME) and consensus molecular subtypes (CMS) of CRC, identifying a predominance of CMS2 tumor cells in several samples. The analysis revealed distinct immune and stromal cell compositions associated with different CMS phenotypes, with mixed CMS1-CMS2 tumors exhibiting higher T and B cell proportions. The study also highlighted the role of specific transcription factors and signaling pathways, such as WNT and VEGF, in tumor progression. Importantly, the research established a connection between the identified cellular interactions and the molecular features of CRC, providing insights into the mechanisms underlying tumor heterogeneity and progression. Overall, this work underscores the utility of ST in elucidating the complex cellular architecture and communication processes within CRC tumors.