تحليل التعبير في مجموعات البيانات العامة والأنسجة السريرية المقترنة
أُعيد تقييم نمط تعبير 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 ثلاث عينات طبيعية فقط. وأظهر تحليل استبعاد عينة طبيعية واحدة أيضاً أن تقدير 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 المشتركة، لم تكن مقارنة كابلان-ماير القائمة على القطع المتوسطي معنوية (log-rank p = 0.857)، كما لم يظهر تحليل انحدار كوكس أحادي المتغير أي ارتباط معنوي (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). وأظهر تحليل الحساسية لبيانات TCGA RNA-seq أن تعبير KIF20A وTPX2 وCDCA5 وCCNB1 وCDK1 كان أعلى في الأنسجة الورمية مقارنة بالأنسجة الطبيعية، بينما كان تعبير PLP1 وNRXN1 وGRIK3 وKIF5C أقل في الأنسجة الورمية. وكان تحليل الارتباط مع تعبير miR-192 في TCGA مختلطاً: حيث أظهرت جينات CDCA5 وCCNB1 وCDK1 ارتباطات إيجابية، بينما أظهرت PLP1 وNRXN1 وKIF5C ارتباطات سلبية. وسلط تحليل الإثراء الضوء على مصطلحات ممثلة تتعلق بتمثيل الأدوية بواسطة سيتوكروم P450، وتمثيل المواد الغريبة عن الجسم بواسطة سيتوكروم P450، وتمثيل الريتينول، والمصطلحات المتعلقة بالمشبك الغلوتاماتي، وتنظيم ما قبل المشبك (الشكل 6). وتدعم هذه النتائج الصلة البيولوجية لشبكة الجينات المركزية ولكنها لا تثبت استهدافاً مباشراً لـ miR-192-5p.
مقايسة هجرة الخلايا
تم تحويل خلايا 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 في المختبر.
تحليل ارتباط البصمة المناعية
أُعيد تقييم الارتباط بين تعبير 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) والأنسجة المزدوجة. (A) تحليل miRNA-seq من TCGA/GDC يظهر إشارة hsa-mir-192 أعلى على مستوى قاعدة البيانات في عينات أورام CRC مقارنة بالعينات الطبيعية. (B) تحليل qRT-PCR لعينات أنسجة محلية مزدوجة يظهر تعبيراً أقل لـ miR-192-5p الناضج في أنسجة CRC مقارنة بالأنسجة المجاورة غير السرطانية. يرجى النقر هنا لعرض نسخة أكبر من هذا الشكل.

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

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

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

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

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

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

الشكل 9: الارتباط بين تعبير miR-192 في TCGA ودرجات البصمة المرتبطة بالمناعة في عينات أورام CRC. أظهر التحليل ارتباطات مناعية غير متجانسة، مع وجود عدة ارتباطات سلبية ضعيفة، وارتباطات إيجابية بشكل أساسي لبصمات الخلايا الحمضية وخلايا Th17. يرجى النقر هنا لعرض نسخة أكبر من هذا الشكل.
| الخصائص | تعبير منخفض لـ TCGA miR-192 | تعبير مرتفع لـ TCGA miR-192 | القيمة الاحتمالية (p value) | الطريقة |
| n | 308 | 308 | | |
| مرحلة T، العدد (n) (%) | 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، العدد (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، العدد (n) (%) | 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 المتوازنة، وتحليلات استبعاد عينة طبيعية واحدة (leave-one-normal-out)، وتحليلات الارتباط السريري المرضي، وتحليلات البقاء على قيد الحياة، وملخصات تعبير وارتباط الجينات المركزية (hub-gene)، ونتائج ارتباط البصمة المناعية.يرجى النقر هنا لتحميل هذا الملف.
المجلد التكميلي 1: البيانات التجريبية الأولية التي تدعم تجارب qRT-PCR، وترصيع الخلايا، والتئام الجروح.تتضمن هذه الملفات قيم Ct الأصلية لـ qRT-PCR، وبيانات قياس التدفق الخلوي/الترصيع، وقياسات عرض الجرح، وصور مجهرية ممثلة استُخدمت في التحليلات المعروضة في المخطوطة.يرجى النقر هنا لتحميل هذا الملف.