Zbieranie danych i przetwarzanie wstępne
Wszystkie dane przeanalizowane w niniejszym badaniu zostały uzyskane z publicznie dostępnej bazy danych TCGA wymienionej w Tabeli Materiałów. Niniejsze badanie ściśle przestrzegło Polityki Dotyczącej Udostępniania Danych Genomowych Narodowych Instytutów Zdrowia (NIH) oraz wytycznych publikacji dostarczonych przez TCGA. Ponieważ zbiory danych składają się z anonimizowanych, dostępnych publicznie danych klinicznych i wieloomowych, niniejsze badanie było zwolnione z dalszej aprobaty Instytucjonalnej Rady Etycznej (IRB). Dane transkryptomowe i odpowiadające im informacje kliniczne dotyczące LUSC zostały pobrane z bazy danych Cancer Genome Atlas (TCGA). Początkowo kohortę stanowiło 551 próbek (502 guzów i 49 normalnych). Aby zapewnić solidność tego modelu prognostycznego, pacjenci brakujący pełne dane dotyczące przeżycia ogólnego (OS), statusu przeżycia lub istotnych cech klinicznych zostali wykluczeni z dalszej analizy przeżycia. Po tej surowej filtracji kohorty, ustanowiono kohorty analitycznej 470 pacjentów z LUSC. Ci 470 pacjentów zostali następnie losowo podzieleni na zbiór treningowy (n = 235) i zbiór testowy (n = 235) w celu konstrukcji i walidacji modelu ryzyka. Łącznie zebrano 792 geny związane ze stresem oksydacyjnym za pośrednictwem bazy danych GeneCards, wymienionej w Tabeli Materiałów, z wykorzystaniem kryterium inkluzji wyniku trafności ≥ 7. Do przesiewania różnic w ekspresji genów związanych ze stresem oksydacyjnym zastosowano pakiet limma (wersja 3.66.0). Wykresy wulkaniczne stworzono za pomocą pakietu ggplot2 (wersja 4.0.2). Pomiędzy grupą treningową a grupą testową nie było zauważalnej różnicy w cechach klinicznych (p > 0.05). Zestaw treningowy został użyty do konstrukcji modelu, a zestaw testowy do wewnętrznej walidacji.
Konstrukcja i weryfikacja modelu ryzyka
Aby zidentyfikować lncRNA związane ze stresem oksydacyjnym, przeprowadzono analizę koekspresji Pearsona. Autorzy ocenili korelacje między profilami ekspresji zidentyfikowanych genów związanych ze stresem oksydacyjnym a wszystkimi zaanonsowanymi lncRNA w zbiorze TCGA-LUSC. Kryteria korelacji dla wyboru koekspresowanych lncRNA zostały ustawione na współczynniku korelacji |R| > 0.4 i wartości p < 0.001. Na podstawie tych progów wyprowadzono i wyodrębniono łącznie 6088 kandydatów lncRNA związanych ze stresem oksydacyjnym. Następnie te 6088 kandydatów lncRNA zostało ocenione za pomocą regresji LASSO i uzupełniających podejść analitycznych do konstrukcji prognostycznego modelu ryzyka. Wykorzystując przetworzone dane, uczestnicy pacjentów zostali sklasyfikowani na wysokie i niskie warstwy ryzyka.
Niezależne czynniki i krzywe ROC
Wykorzystano uniwariantne i multiwariantne modele regresji Coxa do zidentyfikowania niezależnych zmiennych prognostycznych. W szczególności, aby kontrolować rodzinne wskaźniki błędów podczas wielu uniwariantnych analiz regresji Coxa (n = 236 testów), zastosowano korekcję Bonferroniego, ustawiając surowy dostosowany próg istotności na p < 0.00021 (obliczony jako 0.05 / 236).
Analiza przeżycia i analiza głównych składowych
Pakiet survival (wersja 3.8-6) został użyty do obliczenia wskaźników przeżycia ogólnego (OS), a analiza głównych składowych (PCA) została zastosowana do oceny solidności modelu ryzyka.
Nomogramma
Wskaźnik dopasowania został oceniony za pomocą pakietu rms w R w celu oceny wiarygodności nomogramma. Obciążenie mutacyjne guza (TMB) zostało ocenione i wizualizowane za pomocą pakietu maftools (wersja 2.22.0). Algorytmy, w tym ssGSEA, zostały użyte do wykrycia infiltracji immunologicznej. Ponadto pakiet ggpubr (wersja 0.6.3) został użyty do porównania punktów kontrolnych odpowiedzi immunologicznej między grupami ryzyka. Algorytm ESTIMATE został użyty do charakteryzacji mikrośrodowiska guza (TME).
Analiza funkcjonalna
Pakiet clusterProfiler (wersja 4.14.6) w R został użyty do analizy wzbogacania. Aby dalej zbadać funkcjonalne ścieżki, przeprowadzono analizę GSEA. Za pomocą Cytoscape stworzono sieć koekspresji dla obserwacji.
Metody eksperymentalne
Warunki hodowli komórek i grupowanie
Ludzka linia komórek raka płaskonabłonkowego płuc NCI-H520 (RRID: CVCL_1566; wymieniona w Tabeli Materiałów) była hodowana w podłożu RPMI-1640 uzupełnionym o 10% (v/v) surowicy wołowej oraz 1% penicyliny-streptomycyny. Komórki w fazie logarytmicznej zostały posadzone w 6-dołkowych płytkach hodowlanych. Po pomyślnym przyczepianiu, komórki zostały losowo przydzielone do tr