تحليل التعبير في مجموعات البيانات العامة والأنسجة السريرية المقترنة
أعيد تقييم نمط تعبير miR-192 باستخدام بيانات miRNA-seq من TCGA-COAD وTCGA-READ المستمدة من GDC (الشكل 1). في هذه الملفات، تم تعيين الإدخال على مستوى قاعدة البيانات على أنه hsa-mir-192 بدلاً من تعيينه مباشرة كنوع hsa-miR-192-5p الناضج؛ وبناءً عليه، فُسرت نتيجة قاعدة البيانات العامة هذه على أنها إشارة TCGA لـ miR-192. في مجموعة COADREAD المجمعة، تم تحليل 616 ورماً أولياً و11 عينة طبيعية من الأنسجة الصلبة. وأظهرت عينات الأورام تعبيراً أعلى لـ miR-192 مقارنة بالعينات الطبيعية (الوسيط 16.05 مقابل 13.46 log2[RPM + 1]، وقيمة Mann-Whitney p = 9.96 x 10-8). ولوحظ نفس الاتجاه المتمثل في زيادة التعبير في الأورام في COAD وحده (455 ورماً و8 عينات طبيعية، p = 9.89 x 10-8) وفي READ وحده (161 ورماً و3 عينات طبيعية، p = 1.94 x 10-5). كما أظهر تحليل TCGA المقترن تعبيراً أعلى في الأورام في 11 حالة متطابقة (الوسيط لفرق الورم-الطبيعي = 3.13 log2[RPM + 1]، وقيمة Wilcoxon p = 0.0029). وظلت هذه النتائج متضاربة مع نتيجة qRT-PCR المقترنة المحلية، والتي دعمت انخفاض تعبير miR-192-5p الناضج في أنسجة CRC (الشكل 2).
الأداء التشخيصي والارتباطات المرضية السريرية
تكرر تحليل ROC بناءً على بيانات TCGA باستخدام تعريف التمييز بين العينات الورمية والعينات الطبيعية الموضح أعلاه (الشكل 3). بلغت قيمة AUC نحو 0.968 لمجموعة COADREAD المدمجة، و0.947 لـ COAD وحده، و0.994 لـ READ وحده. وأسفرت عملية إعادة أخذ العينات المتوازنة الإضافية عن قيم وسيطة لـ AUC بلغت 0.983 لـ COADREAD، و0.984 لـ COAD، و1.000 لـ READ؛ ومع ذلك، كانت فترات الثقة 95% واسعة لكل من COADREAD (0.868–1.000) وCOAD (0.781–1.000)، كما اشتملت مجموعة READ على ثلاث عينات طبيعية فقط. وأظهر تحليل استبعاد عينة واحدة (Leave-one-out) للعينات الطبيعية أيضاً أن تقدير ROC كان حساساً لصغر حجم المجموعة الطبيعية. وبناءً على ذلك، تم الاحتفاظ بنتيجة ROC كدليل استكشافي فقط على الفصل بين الورم والأنسجة الطبيعية، وليس كدليل كافٍ ليكون مؤشراً حيوياً للتشخيص السريري. وأظهر تحليل الارتباط العيادي المرضي بناءً على وسيط تعبير miR-192 في TCGA ارتباطات معنوية مع مرحلة N (p = 0.019) والعمر (p = 0.004)، بينما لم تكن مرحلة T (p = 0.110) ومرحلة M (p = 0.932) معنوية (الجدول 1).
تحليل البقاء
لم يدعم تحليل البقاء الإجمالي وجود ارتباط إنذاري ملحوظ للتعبير عن miR-192 في بيانات TCGA لسرطان القولون والمستقيم CRC (الشكل 3). وفي مجموعة COADREAD المدمجة، لم تكن مقارنة Kaplan-Meier المعتمدة على القطع الوسيط ذات دلالة إحصائية (log-rank p = 0.857)، كما أظهر تحليل انحدار Cox أحادي المتغير عدم وجود ارتباط ملحوظ (HR = 0.934 لكل وحدة تعبير log2، p = 0.374). ولوحظ الاستنتاج نفسه في COAD (log-rank p = 0.889؛ Cox p = 0.493) وREAD (log-rank p = 0.807؛ Cox p = 0.716). وبناءً عليه، فإن البيانات لا تدعم الادعاء بأن miR-192 يمثل مؤشراً حيوياً إنذارياً قوياً في سرطان القولون والمستقيم CRC.
الجينات المستهدفة المرشحة والإثراء الوظيفي
تم إجراء تحليل الجينات المستهدفة المرشحة من خلال دمج تنبؤات miRWalk مع الجينات ذات التعبير التفاضلي من مجموعات بيانات GEO وهي GSE89076 وGSE156355. وقد حدد هذا الفرز 93 جيناً مرشحاً لتحليلات الشبكة والإثراء اللاحقة (الشكل 4). وأبقى تحليل التفاعل البروتيني على KIF20A، وTPX2، وCDCA5، وCCNB1، وCDK1، وPLP1، وNRXN1، وGRIK3، وKIF5C كجينات مركزية استكشافية (الشكل 5). وأظهر تحليل حساسية RNA-seq من TCGA أن KIF20A، وTPX2، وCDCA5، وCCNB1، وCDK1 كانت أعلى في الأنسجة الورمية منها في الأنسجة الطبيعية، بينما كانت PLP1، وNRXN1، وGRIK3، وKIF5C أقل في الأنسجة الورمية. وكان تحليل الارتباط مع تعبير miR-192 في TCGA مختلطاً؛ حيث أظهرت CDCA5 وCCNB1 وCDK1 ارتباطات إيجابية، بينما أظهرت PLP1 وNRXN1 وKIF5C ارتباطات سلبية. وسلط تحليل الإثراء الضوء على مصطلحات ممثلة تتعلق باستقلاب الأدوية بواسطة سيتوكروم P450، واستقلاب المواد الغريبة بواسطة سيتوكروم P450، واستقلاب الريتينول، والمصطلحات المتعلقة بالتشابك العصبي الجلوتاماتي، وتنظيم ما قبل التشابك (الشكل 6). تدعم هذه النتائج الأهمية البيولوجية لشبكة الجينات المركزية ولكنها لا تثبت استهدافاً مباشراً لـ miR-192-5p.
مقايسة هجرة الخلايا
تم إجراء عفن (transfection) لخلايا HT29 باستخدام محاكيات miR-192 لتقييم التأثيرات الوظيفية لزيادة تنظيم miR-192. وأكد الفحص بالمجهر الفلوري نجاح عملية العفن في مجموعات محاكي miR-192-5p ومجموعات العفن الضابطة السلبية (الشكل 7). وفي اختبار التئام الجروح، كان معدل هجرة خلايا HT29 التي تم عفنها بـ miR-192 أقل من معدل هجرة خلايا HT29 غير المعفنة والخلايا المعفنة بالضابط السلبي عند كل من 24 h و 48 h (الشكل 8). ففي غضون 24 h، بلغ معدل الهجرة 9.75 +/- 2.43% في المجموعة المعفنة بـ miR-192 مقارنة بـ 15.69 +/- 3.47% في المجموعة غير المعفنة و 15.45 +/- 3.92% في المجموعة الضابطة السلبية (P = 0.010). وعند 48 h، بلغ معدل الهجرة 20.04 +/- 2.54% في المجموعة المعفنة بـ miR-192 مقارنة بـ 26.84 +/- 7.65% و 27.98 +/- 6.78% في المجموعتين الضابطتين، على التوالي (P = 0.021). وتشير هذه النتائج إلى أن زيادة تنظيم miR-192 تثبط هجرة خلايا HT29 في المختبر (in vitro).
تحليل ارتباط البصمة المناعية
أُعيد تقييم الارتباط بين تعبير miR-192 في TCGA والبصمات المرتبطة بالمناعة باستخدام بيانات تسلسل الحمض النووي الريبوزي (RNA-seq) للأورام الأولية (الشكل 9). شمل تحليل ارتباط درجة البصمة ما مجموعه 425 عينة ورمية ذات بيانات متطابقة من miRNA وRNA-seq. أظهر تعبير miR-192 ارتباطات سلبية ضعيفة مع بصمات خلايا مناعية متعددة، بما في ذلك خلايا Tfh (rho = -0.205, p = 2.15e-05)، والخلايا الشجيرية المنشطة (rho = -0.186, p = 0.00011)، وخلايا Treg (rho = -0.185, p = 0.00013)، والخلايا الشجيرية البلازمية (rho = -0.169, p = 0.00046)، والبلاعم (rho = -0.149, p = 0.0020)، وخلايا Tcm (rho = -0.136, p = 0.0050)، وخلايا NK CD56bright (rho = -0.132, p = 0.0065)، وخلايا Tem (rho = -0.129, p = 0.0078). كما لوحظت ارتباطات إيجابية ضعيفة لبصمات الخلايا الحمضية (rho = 0.139, p = 0.0040) وخلايا Th17 (rho = 0.121, p = 0.0129). وبناءً على ذلك، يبدو الارتباط المناعي لـ miR-192 في سرطان القولون والمستقيم (CRC) متباينًا، ولا ينبغي وصفه بأنه إيجابي بشكل موحد عبر جميع الرشاحات المناعية.
توافر البيانات
تتوفر مجموعات البيانات العامة التي تم تحليلها في هذه الدراسة من بوابة GDC/TCGA وقاعدة بيانات Gene Expression Omnibus (GEO) تحت أرقام الوصول GSE89076 و GSE156355. تتوفر تحليلات المختبر الجاف وعمليات إعادة التحليل التي تم إجراؤها خلال هذه الدراسة في الجدول التكميلي 1. أما البيانات التجريبية الأساسية التي تدعم تحليلات qRT-PCR والتئام الجروح، بما في ذلك قيم Ct الخام لـ qRT-PCR، وقياسات عرض الجرح، وصور المجهر الأصلية، والبيانات التجريبية ذات الصلة، فهي متوفرة في المجلد التكميلي 1.

الشكل 1: نظرة عامة على تعبير miR-192 في مختلف أنواع السرطانات في قاعدة بيانات TCGA. يمثل المحور السيني اختصارات الأورام، ويمثل المحور الصادي التعبير النسبي لـ miR-192. يشير اللون الأزرق إلى الأنسجة الطبيعية، واللون الأحمر إلى أنسجة الورم. ويتم الإشارة إلى الدلالة الإحصائية كما يلي: * = p < 0.05، ** = p < 0.01، *** = p < 0.001؛ NS = غير دال إحصائياً. يرجى النقر هنا لعرض نسخة أكبر من هذا الشكل.

الشكل 2: تعبير miR-192/miR-192-5p في مجموعات بيانات سرطان القولون والمستقيم (CRC) والأنسجة المقترنة. (أ) تحليل miRNA-seq من TCGA/GDC يظهر إشارة hsa-mir-192 على مستوى قاعدة البيانات أعلى في عينات أورام CRC منها في العينات الطبيعية. (ب) تحليل qRT-PCR لعينات أنسجة محلية مقترنة يظهر تعبيراً أقل لـ miR-192-5p الناضج في أنسجة CRC مقارنة بالأنسجة غير السرطانية المجاورة. يرجى النقر هنا لعرض نسخة أكبر من هذا الشكل.

الشكل 3تحليلات تشخيصية وبقائية استكشافية بناءً على تعبير hsa-mir-192 في أطلس جينوم السرطان (TCGA). (أ–جمنحنيات كابلان-ماير للبقاء على قيد الحياة. لم يدعم تحليل البقاء وجود ارتباط إنذاري ذي دلالة إحصائية.دمنحنى خصائص تشغيل المستقبل (ROC) للتمييز بين الورم والأنسجة الطبيعية في مجموعة TCGA. يجب تفسير المساحة تحت المنحنى (AUC) بحذر نظرًا لقلة عدد العينات الطبيعية. يرجى النقر هنا لعرض نسخة أكبر من هذا الشكل.

الشكل 4: فحص الجينات المستهدفة المرشحة وتحليل التفاعل البروتيني. (A) مخطط فين يوضح التداخل بين توقعات أهداف miRNA والجينات ذات التعبير المتباين من GSE89076 وGSE156355. (B) شبكة التفاعل بين البروتينات لـ 93 جيناً مرشحاً. تمثل العقد الجينات، بينما تمثل الحواف الارتباطات الجينية المتوقعة أو المنسقة. يرجى النقر هنا لعرض نسخة أكبر من هذا الشكل.

الشكل 5: فحص الجينات المركزية من شبكة الأهداف المرشحة. (A,B) وحدات الشبكة التي تم تحديدها من خلال فحص الوحدات. (C) ترتيب الجينات المركزية بناءً على اتصال الشبكة. (D) الأهداف المركزية المرشحة التي تم الاحتفاظ بها بعد دمج الفحص القائم على الوحدات والدرجة. يرجى النقر هنا لعرض نسخة أكبر من هذا الشكل.

الشكل 6: مصطلحات إثراء GO وKEGG المختارة للجينات المرشحة ووحدات الشبكة ذات الدلالة. يلخص المخطط الفقاعي الوظائف والمسارات البيولوجية المثرية والممثلة. وتعد نتائج الإثراء استكشافية وتتطلب تحققاً تجريبياً. يرجى النقر هنا لعرض نسخة أكبر من هذا الشكل.

الشكل 7: صور ترانسفكشن خلايا HT29. (أ) خلايا HT29 نموذجية. (ب) صورة فلورية لخلايا HT29 تم ترانسفكشنها بمحاكي miR-192-5p. (ج) صورة فلورية لخلايا HT29 تم ترانسفكشنها بضابط سلبي con238. شريط المقياس = 100 µm. يرجى النقر هنا لعرض نسخة أكبر من هذا الشكل.

الشكل 8: مقايسة التئام الجروح في خلايا HT29. (أ) صور تمثيلية للخدش عند 0 h و 24 h و 48 h في خلايا HT29 غير منقولة، وخلايا HT29 منقولة بـ miR-192-5p mimic، وخلايا HT29 منقولة بـ con238. (ب) تقدير كمي لمعدلات التئام الجروح عند 24 h و 48 h، موضحة كنسب مئوية. شريط المقياس = 100 µm. يرجى النقر هنا لعرض نسخة أكبر من هذا الشكل.

الشكل 9: الارتباط بين تعبير miR-192 في قاعدة بيانات TCGA ودرجات البصمة المناعية في عينات أورام CRC. أظهر التحليل ارتباطات مناعية غير متجانسة، مع وجود العديد من الارتباطات السلبية الضعيفة، وارتباطات إيجابية، خاصة بالنسبة لبصمات الخلايا الحمضية (eosinophil) وخلايا Th17. يرجى النقر هنا لعرض نسخة أكبر من هذا الشكل.
| الخصائص | تعبير منخفض لـ TCGA miR-192 | تعبير مرتفع لـ TCGA miR-192 | القيمة الاحتمالية p value | الطريقة |
| n | 308 | 308 | | |
| مرحلة T، العدد (٪) | 0.11 | اختبار مربع كاي Chisq.test |
| T1 | 9 (1.5%) | 11 (1.8%) | | |
| T2 | 44 (7.2%) | 59 (9.6%) | | |
| T3 | 209 (34.1%) | 211 (34.4%) | | |
| T4 | 43 (7%) | 27 (4.4%) | | |
| مرحلة N، العدد (٪) | 0.019 | اختبار مربع كاي Chisq.test |
| N0 | 156 (25.5%) | 190 (31%) | | |
| N1 | 88 (14.4%) | 63 (10.3%) | | |
| N2 | 61 (10%) | 54 (8.8%) | | |
| مرحلة M، العدد (٪) | 0.932 | اختبار مربع كاي Chisq.test |
| M0 | 214 (39.9%) | 234 (43.7%) | | |
| M1 | 41 (7.6%) | 47 (8.8%) | | |
| العمر، الوسيط (IQR) | 66 (56, 75) | 69 (60, 77) | 0.004 | اختبار ويلكوكسون Wilcoxon |
جدول 1: الارتباط بين مجموعات تعبير miR-192 في TCGA والخصائص السريرية المرضية في CRC. تم تحديد مجموعات التعبير العالي والمنخفض بناءً على متوسط مستوى تعبير miR-192 في TCGA. في الجدول الملخص، كانت مرحلة N والعمر ذوي دلالة إحصائية، بينما لم تكن مرحلة T ومرحلة M ذات دلالة إحصائية.
الجدول التكميلي 1: نتائج إعادة التحليل المختبري الجاف التي تم إنتاجها لهذه الدراسة. يتضمن دفتر العمل مقارنات التعبير الجيني لـ TCGA/GDC، وتحليلات حساسية ROC المتوازنة، وتحليلات استبعاد عينة طبيعية واحدة، وتحليلات الارتباط العيادي المرضي، وتحليلات البقاء، وملخصات تعبير وارتباط الجينات المركزية، ونتائج ارتباط التوقيع المناعي.يرجى النقر هنا لتنزيل هذا الملف.
المجلد التكميلي 1: البيانات التجريبية الخام الداعمة لتجارب qRT-PCR، وترصيف الخلايا، والتئام الجروح.تتضمن هذه الملفات قيم Ct الأصلية لـ qRT-PCR، وبيانات قياس التدفق الخلوي/الترصيف، وقياسات عرض الجرح، وصور مجهرية ممثلة استُخدمت في التحليلات المعروضة في المخطوطة.يرجى النقر هنا لتحميل هذا الملف.