רכישת נתונים ועיבוד מוקדם
כל הנתונים שנותחו במחקר זה התקבלו ממאגר הנתונים הנגיש לציבור של TCGA המופיע בטבלת החומרים. מחקר זה דבק בקפדנות במדיניות שיתוף הנתונים הגנומיים של המכון הלאומי לבריאות (NIH) ובהנחיות הפרסום שסיפקה TCGA. מכיוון שמאגרי הנתונים מורכבים מנתונים קליניים ומולטי-אומיקס שאינם מזוהים, זמינים לציבור, מחקר זה היה פטור מאישור נוסף של ועדת הביקורת המוסדית (IRB). נתוני טרנסקריפטום ומידע קליני מקביל עבור LUSC הורדו ממאגר אטלס הגנום של סרטן (TCGA). בתחילה, הקבוצה כללה 551 דגימות (502 גידולים ו-49 תקינות). כדי להבטיח את עמידות המודל הפרוגנוסטי הזה, מטופלים שחסרו נתוני הישרדות כוללים מלאים, מצב הישרדות או מאפיינים קליניים חיוניים הוצאו מניתוח ההישרדות בהמשך. לאחר תהליך סינון קפדני זה, הוקם קבוצה אנליטית סופית של 470 מטופלי LUSC. 470 המטופלים הללו חולקו באקראי לקבוצת אימון (n = 235) וקבוצת בדיקות (n = 235) לבניית מודל סיכון ואימותו. בסך הכל נאספו 792 גנים הקשורים ללחץ חמצוני באמצעות מאגר GeneCards, המופיע בטבלת החומרים, תוך שימוש בקריטריון הכללה של ציון רלוונטיות ≥ 7. חבילת הלימה (גרסה 3.66.0) שימשה לסינון גנים של לחץ חמצוני המובע באופן שונה. חלקות הר געש יוצרו באמצעות חבילת ggplot2 (גרסה 4.0.2). בין קבוצת האימון לקבוצת הבדיקה, לא נראתה שונות ניכרת במאפיינים הקליניים (p > 0.05). סט האימון שימש לבניית מודלים, וסט הבדיקה שימש לאימות פנימי.
בניית ואימות מודל הסיכון
כדי לזהות lncRNA הקשורים ללחץ חמצוני, בוצע ניתוח ביטוי משותף של פירסון. המחברים העריכו את המתאם בין פרופילי הביטוי של גנים הקשורים ללחץ חמצוני שזוהו במונחים שונים לבין כל ה-lncRNA המוערכים במאגר הנתונים של TCGA-LUSC. קריטריוני המתאם לבחירת lncRNA מבוטאים במשותף נקבעו במקדם קורלציה |R| > 0.4 וערך p < 0.001. בהתבסס על ספים אלו, נגזרו והופקו בסך הכל 6,088 מועמדים ל-lncRNA הקשורים ללחץ חמצוני. בהמשך, 6,088 מועמדים אלו הוערכו באמצעות רגרסיית LASSO ושיטות אנליטיות משלימות, לבניית מודל הסיכון הפרוגנוסטי. באמצעות הנתונים המעובדים, משתתפי המטופלים סווגו לשכבות סיכון גבוה ונמוך.
גורמים בלתי תלויים ועקומות ROC
מודלים של רגרסיה חד-משתנית ורב-משתנית של קוקס שימשו לזיהוי משתנים פרוגנוסטיים עצמאיים. בפרט, כדי לשלוט בשיעור השגיאה לפי משפחה במהלך ניתוחי רגרסיה חד-משתניים של קוקס (n = 236 מבחנים), יושם תיקון בונפרוני, וקבע את סף המובהקות המותאם בקפדנות ל-p < 0.00021 (מחושב כ-0.05 / 236).
ניתוח הישרדות וניתוח רכיבים עיקריים
חבילת ההישרדות (גרסה 3.8-6) שימשה לחישוב שיעורי ההישרדות הכוללים (OS), וניתוח רכיבים עיקריים (PCA) שימש להערכת עמידות מודל הסיכון.
נומוגרמה
מדד ההתאמה הוערך באמצעות חבילת rms ב-R כדי להעריך את אמינות הנומוגרמה. נטל המוטציה של הגידול (TMB) הוערך והוצג באמצעות חבילת maftools (גרסה 2.22.0). אלגוריתמים, כולל ssGSEA, שימשו לזיהוי חדירה חיסונית. בנוסף, חבילת ggpubr (גרסה 0.6.3) שימשה להשוואת נקודות בדיקה חיסוניות בין קבוצות סיכון. אלגוריתם ESTIMATE שימש לאפיון מיקרו-סביבת הגידול (TME).
ניתוח פונקציונלי
חבילת clusterProfiler (גרסה 4.14.6) ב-R שימשה לניתוח העשרה. כדי לחקור מסלולים פונקציונליים לעומק, בוצע ניתוח GSEA. באמצעות Cytoscape, נוצרה רשת קו-ביטוי לתצפית.
שיטות ניסוי
תנאי תרבית תאים וקיבוץ
קו תאי קרצינומה של תאי קשקש ריאות אנושיים NCI-H520 (RRID: CVCL_1566; רשום בטבלת החומרים) תורתח במדיום RPMI-1640 שהושלף עם סרום בקר עוברי של 10% (v/v) ו-1% פניצילין-סטרפטומיצין. תאי לוג-פאזה הוזרעו ללוחות תרבית של 6 בארות. לאחר הקפדה מוצלחת, התאים חולקו באקראי לשלוש קבוצות ניסוי מובחנות: קבוצת ביקורת ריקה שנשמרה בתנאי תרבית סטנדרטיים, קבוצת ביקורת שלילית (si-NC) שהועברה ל-siRNA לא מכוון, וקבוצת הנקה ניסיונית (si-LINC01615) שהועברה עם siRNA ספציפי ל-LINC01615.
טרנספקציה תאית
החומרים si-LINC01615 ו-si-NC הופשרו על קרח. להכנת קומפלקס הטרנספקציה, 95 מיקרוליטר של מדיום RPMI-1640 ללא סרום הועברו לפיפטה לצינורות צנטריפוגות סטריליות, ולאחר מכן הוספנו ברצף 3 מיקרוגרם של siRNA ו-5 מיקרו-ליטר של תגובת טרנספקציה Lipo3000. אותו הליך הוחל הן לקבוצות siLINC01615 והן לקבוצות siNC. התמיסות עורבבו בעדינות ודגרו בטמפרטורת החדר במשך 5 דקות כדי להקל על היווצרות מורכבת, מה שהביא לנפח כולל של כ-200 מיקרו-ליטר לכל צינור. לאחר מכן, תערובת זו חולקה באופן אחיד לבארות התרבות המיועדות המכילות 800 מיקרוליטר של מדיום בסיס RPMI-1640. לאחר תקופת דגירה של 6 שעות בטמפרטורה של 37 מעלות צלזיוס, מדיום הטרנספקציה נשאף בזהירות והוחלף במדיום תרבית טרי ומלא.
בדיקת ריפוי פצעים (שריטות)
למבחן ריפוי הפצע, תאי NCI-H520 שגדלו לוגריתמית (RRID: CVCL_1566) זרעו בלוחות תרבית של 6 בארות בצפיפות של 5.6 ×-105 תאים לבאר. בהגעה ל-95–100% מפגש מגוון, נוצר פצע ליניארי באופן מלאכותי על ידי שריטה אנכית של שכבת התא במרכז הבאר באמצעות קצה פיפטה סטרילי בנפח 200 מיקרוליטר. הבארות נשטפו לאחר מכן שלוש פעמים באמצעות PBS להסרת תאים מנותקים ופסולת תאית, ולאחר מכן דגרו במדיום RPMI-1640 ללא סרום. תמונות של אותם אזורים פצועים צולמו ב-0, 24 ו-48 שעות באמצעות מיקרוסקופ הפוך. אזור ההחלמה של הפצע נמדד על ידי מדידת מרחק הפער באמצעות ImageJ (RRID: SCR_003070; רשום בטבלת החומרים), וקצב הנדידה חושב ביחס לאזור הפצע ההתחלתי ב-0 שעות.
בדיקות הגירה ופלישה של טרנסוול
קיבולת הנדידה והפלישה התאית הוערכו באמצעות תאי טרנסוול בעלי 24 בארות המצוידים בכניסות ממברנה פוליקרבונטית בגודל נקבוביות של 8.0 מיקרון , המפורטות בטבלת החומרים. למבחן הפלישה, תאי האפיקל צופו מראש ב-50 מיקרו ליטר של Matrigel (מדולל 1:8 במדיום ללא סרום; מופיע בטבלת החומרים) ודגרו בטמפרטורה של 37°C במשך שעתיים לפולימריזציה, בעוד שבדיקת ההגירה השתמשה בהכנסות לא מצופות. לאחר טרנספקציה של תאים ותקופת רעב סרום של 24 שעות, תאי NCI-H520 נאספו, הושמו מחדש במדיום RPMI-1640 ללא סרום, והוזרעו לתאי האפיקלי בצפיפות של 1.8 ×-105 תאים לבאר ב-200 מיקרוליטר. תאי הבסיס מולאו ב-600 מיקרוליטר של מדיום RPMI-1640 בתוספת 10% FBS ככימואטרטנט. לאחר דגירה של 48 שעות בטמפרטורה של 37 מעלות צלזיוס עם 5% CO₂, הוסרו ושטפו שלוש פעמים עם PBS. התאים שנותרו על פני השטח האפיקליים של הממברנה נלקחו בעדינות באמצעות מטוש כותנה רטוב. התאים שהיגרו או פלשו למשטח הבזולטרלי הוקבעו עם 4% פרפורמלדהיד למשך 20 דקות ונצבעו בתמיסת גימסה המופיעה בטבלת החומרים למשך 30 דקות בטמפרטורת החדר. התאים המהגרים והפולשים הוצגו ונספרו בשלושה שדות אופטיים שנבחרו באקראי בכל אינסרט באמצעות מיקרוסקופ הפוך בהגדלה של פי 40.
ניתוח סטטיסטי
כל הניתוחים הסטטיסטיים והוויזואליזציות של הנתונים בוצעו באמצעות תוכנת R (גרסה 4.1.2; RRID: SCR_001905; מופיע בטבלת החומרים). מספרי הגרסאות הספציפיים לכל חבילות ה-R ששימשו במחקר זה פורטו במפורש בתת-הסעיפים המתודולוגיים שלהן. ערך משמעותי הוגדר באופן קונבנציונלי כ-p < 0.05 אלא אם כן דווח אחרת. כאשר בוצעו השוואות מרובות (למשל, רגרסיות קוקס חד-משתניות מרובות), שיטת בונפרוני יושמה בקפדנות לשמירה על שיעור השגיאה לפי משפחה (FWER), כאשר הסף המדויק המותאם (למשל, p < 0.00021) מפורט בסעיף המתודולוגי המתאימה. נתונים כמותיים במבחנה הובאו לידי ביטוי כסטיית תקן ממוצעת ± ונותחו באמצעות מבחן t של סטודנט או ANOVA חד-כיווני לפי הצורך.