Datenerfassung und Vorverarbeitung
Alle in dieser Studie analysierten Daten stammen aus der öffentlich zugänglichen TCGA-Datenbank , die in der Materialtabelle aufgeführt ist. Diese Studie hielt sich strikt an die Genomic Data Sharing Policy des National Institutes of Health (NIH) und die von TCGA vorgegebenen Publikationsrichtlinien. Da die Datensätze aus nicht identifizierten, öffentlich zugänglichen klinischen und Multi-Omics-Daten bestehen, war diese Studie von einer weiteren Genehmigung durch das Institutional Review Board (IRB) ausgenommen. Transkriptomdaten und entsprechende klinische Informationen für LUSC wurden aus der Cancer Genome Atlas (TCGA)-Datenbank heruntergeladen. Anfangs umfasste die Kohorte 551 Proben (502 Tumor und 49 normal). Um die Robustheit dieses prognostischen Modells zu gewährleisten, wurden Patienten ohne vollständige Gesamtüberlebensdaten (OS), Überlebensstatus oder wesentliche klinische Merkmale aus der nachgelagerten Überlebensanalyse ausgeschlossen. Nach diesem strengen Filterprozess wurde eine abschließende analytische Kohorte von 470 LUSC-Patienten eingerichtet. Diese 470 Patienten wurden anschließend zufällig in ein Trainingsset (n = 235) und ein Testset (n = 235) für den Aufbau und die Validierung von Risikomodellen eingeteilt. Insgesamt wurden 792 Gene mit oxidativem Stress über die GeneCards-Datenbank gesammelt, die in der Materialübersicht aufgeführt ist, unter Verwendung eines Einschlusskriteriums mit einem Relevanzscore ≥ 7. Das Limma-Paket (Version 3.66.0) wurde verwendet, um auf unterschiedlich exprimierte oxidative Stressgene zu screenen. Vulkandiagramme wurden mit dem ggplot2-Paket (Version 4.0.2) erstellt. Zwischen der Ausbildungsgruppe und der Testgruppe gab es keine nennenswerten Unterschiede in den klinischen Merkmalen (p > 0,05). Der Trainingssatz wurde für den Modellbau verwendet, der Testsatz für die interne Validierung.
Entwicklung und Verifikation des Risikomodells
Zur Identifizierung von lncRNAs, die mit oxidativem Stress assoziiert sind, wurde eine Pearson-Koexpressionsanalyse durchgeführt. Die Autoren bewerteten die Korrelation zwischen den Expressionsprofilen der identifizierten differenziell exprimierten oxidativ-stressbezogenen Gene und allen annotierten lncRNAs im TCGA-LUSC-Datensatz. Die Korrelationskriterien zur Auswahl ko-exprimierter lncRNAs wurden auf einen Korrelationskoeffizienten |R| > 0,4 und ein p-Wert < 0,001. Basierend auf diesen Schwellenwerten wurden insgesamt 6.088 Kandidaten für oxidative, stressbezogene lncRNAs gewonnen und extrahiert. Anschließend wurden diese 6.088 Kandidaten-lncRNAs mittels LASSO-Regression und komplementären analytischen Ansätzen bewertet, um das prognostische Risikomodell zu erstellen. Mit Hilfe der verarbeiteten Daten wurden die Patiententeilnehmer in Hoch- und Niedrigrisiko-Schichten eingeteilt.
Unabhängige Faktoren und ROC-Kurven
Univariate und multivariate Cox-Regressionsmodelle wurden verwendet, um unabhängige prognostische Variablen zu identifizieren. Konkret wurde zur Kontrolle der familienweisen Fehlerrate während der multiplen univariaten Cox-Regressionsanalysen (n = 236 Tests) die Bonferroni-Korrektur angewendet, wobei die streng angepasste Signifikanzschwelle auf p < 0,00021 (berechnet als 0,05 / 236) gesetzt wurde.
Überlebensanalyse und Hauptkomponentenanalyse
Das Überlebenspaket (Version 3.8-6) wurde zur Berechnung der Gesamtüberlebensraten (OS) verwendet, und die Hauptkomponentenanalyse (PCA) wurde verwendet, um die Robustheit des Risikomodells zu bewerten.
Nomogramm
Der Fit-Index wurde mit dem RMS-Paket in R bewertet, um die Zuverlässigkeit des Nomogramms zu bewerten. Die Tumormutationslast (TMB) wurde mit dem maftools-Paket (Version 2.22.0) bewertet und visualisiert. Algorithmen, darunter ssGSEA, wurden verwendet, um Immuninfiltration zu erkennen. Zusätzlich wurde das ggpubr-Paket (Version 0.6.3) verwendet, um Immunkontrollpunkte zwischen Risikogruppen zu vergleichen. Der ESTIMATE-Algorithmus wurde verwendet, um die Tumor-Mikroumgebung (TME) zu charakterisieren.
Funktionalanalyse
Das clusterProfiler-Paket (Version 4.14.6) in R wurde für die Anreicherungsanalyse verwendet. Um funktionale Wege weiter zu erforschen, wurde eine GSEA-Analyse durchgeführt. Mit Cytoscape wurde ein Koexpressionsnetzwerk zur Beobachtung erstellt.
Experimentelle Methoden
Zellkulturbedingungen und Gruppierungen
Die NCI-H520-Zelllinie des menschlichen Lungenplattenepithelkarzinoms (RRID: CVCL_1566; in der Materialtabelle aufgeführt) wurde in RPMI-1640-Medium kultiviert, ergänzt mit 10 % (v/v) fetalem Rinderserum und 1 % Penicillin-Streptomycin. Die Logaritam-Phasenzellen wurden in 6-Well-Kulturplatten eingesät. Nach erfolgreicher Adhärenz wurden die Zellen zufällig in drei verschiedene experimentelle Gruppen eingeteilt: eine blanke Kontrollgruppe, die unter Standardkulturbedingungen aufrechterhalten wurde, eine negative Kontrollgruppe (si-NC), die mit nicht-zielgerichteter siRNA transfiziert wurde, und eine experimentelle Knockdown-Gruppe (si-LINC01615), die mit LINC01615-spezifischer siRNA transfektiert wurde.
Zelltransfektion
Die si-LINC01615 - und si-NC-Reagenzien wurden auf Eis aufgetaut. Für die Präparation des Transfektionskomplexes wurden 95 μL serumfreies RPMI-1640-Medium in sterile Zentrifugenröhren pipettiert, gefolgt von sequentieller Zugabe von 3 μg der jeweiligen siRNA und 5 μL des Transfektionsreagens Lipo3000. Das gleiche Verfahren wurde sowohl für diesi-LINC01615 - als auch für siNC-Gruppen angewandt. Die Lösungen wurden vorsichtig gemischt und bei Raumtemperatur 5 Minuten inkubiert, um die Komplexbildung zu erleichtern, was zu einem Gesamtvolumen von etwa 200 μL pro Röhre führte. Anschließend wurde dieses Gemisch gleichmäßig in die vorgesehenen Kulturbrunnen verteilt, die 800 μL RPMI-1640 Basalmedium enthalten. Nach einer 6-stündigen Brutzeit bei 37 °C wurde das Transfektionsmedium sorgfältig aspiriert und durch frisches, vollständiges Kulturmedium ersetzt.
Wundheilungstest (Scratch-)Assay
Für den Wundheilungstest wurden logarithmisch wachsende NCI-H520-Zellen (RRID: CVCL_1566) in 6-Well-Kulturplatten mit einer Dichte von 5,6 × 105-Zellen pro Bohrloch ausgesaat. Nach Erreichen von 95–100 % Konfluenz wurde eine lineare Wunde künstlich erzeugt, indem die Zellmonoschicht vertikal über die Mitte des Brunnens mit einer sterilen 200-μL-Pipettenspitze gekratzt wurde. Die Brunnen wurden anschließend dreimal mit PBS gewaschen, um abgetrennte Zellen und Zellreste zu entfernen, und anschließend in serumfreiem RPMI-1640-Medium inkubiert. Bilder der identischen verletzten Bereiche wurden bei 0, 24 und 48 Stunden mit einem umgekehrten Mikroskop aufgenommen. Der Wundheilungsbereich wurde durch Messung der Lückenentfernung mit ImageJ (RRID: SCR_003070; in der Materialtabelle aufgeführt) quantifiziert, und die Migrationsrate wurde relativ zur anfänglichen Wundfläche bei 0 Stunden berechnet.
Transwell-Migration und Invasionstests
Zelluläre Migrations- und Invasionskapazitäten wurden unter Verwendung von 24-Well-Transwell-Kammern bewertet, die mit 8,0 μm porengroßen Polycarbonatmembran-Einsätzen ausgestattet sind, die in der Materialtabelle aufgeführt sind. Für den Invasionstest wurden die apikalen Kammern mit 50 μL Matrigel (1:8 im serumfreien Medium verdünnt; in der Materialtabelle aufgeführt) vorbeschichtet und bei 37 °C für 2 Stunden inkubiert, während der Migrationstest unbeschichtete Einsätze verwendete. Nach der Zelltransfektion und einer 24-stündigen Serum-Hungerphase wurden NCI-H520-Zellen entnommen, in serumfreiem RPMI-1640-Medium wieder suspendiert und mit einer Dichte von 1,8 × 105 Zellen /Brunnen in 200 μL in die apikalen Kammern eingesät. Die Basalkammern wurden mit 600 μL RPMI-1640-Medium gefüllt, ergänzt mit 10 % FBS als Chemomagnet. Nach einer 48-stündigen Inkubation bei 37 °C mit 5 % CO₂ wurden die Einsätze entfernt und dreimal mit PBS gewaschen. Zellen auf der apikalen Oberfläche der Membran wurden vorsichtig mit einem feuchten Wattestäbchen entfernt. Die Zellen, die die basolaterale Oberfläche gewandert oder invadiert hatten, wurden 20 Minuten lang mit 4 % Paraformaldehyd fixiert und bei Raumtemperatur 30 Minuten mit Giemsa-Lösung aus der Materialtabelle gefärbt. Die migrierten und invadierten Zellen wurden visualisiert und über drei zufällig ausgewählte optische Felder pro Insert mit einem umgekehrten Mikroskop mit 40-facher Vergrößerung aufgezählt.
Statistische Analyse
Alle statistischen Analysen und Datenvisualisierungen wurden mit R-Software (Version 4.1.2; RRID: SCR_001905; aufgeführt in der Materialtabelle). Die spezifischen Versionsnummern aller in dieser Studie verwendeten R-Pakete wurden in ihren jeweiligen methodischen Unterabschnitten explizit detailliert dargestellt. Ein signifikanter Wert wurde konventionell als p < 0,05 definiert, sofern nicht anders angegeben. Wo mehrere Vergleiche durchgeführt wurden (z. B. die multiplen univariaten Cox-Regressionen), wurde die Bonferroni-Methode streng angewandt, um die familienweise Fehlerquote (FWER) aufrechtzuerhalten, wobei die exakt angepasste Schwelle (z. B. p < 0,00021) im jeweiligen methodologischen Unterabschnitt detailliert dargestellt wurde. Quantitative In-vitro-Daten wurden als Mittelwert ± Standardabweichung angegeben und je nach Bedarf mit Student's t-Test oder Einweg-ANOVA analysiert.