Aquisição e pré-processamento de dados
Todos os dados analisados neste estudo foram obtidos do banco de dados público da TCGA listado na Tabela de Materiais. Este estudo seguiu rigorosamente a Política de Compartilhamento de Dados Genômicos do Instituto Nacional de Saúde (NIH) e as diretrizes de publicação fornecidas pela TCGA. Como os conjuntos de dados consistem em dados clínicos e multi-ômicos desidentificados, disponíveis publicamente, este estudo foi isento de aprovação adicional do Conselho de Revisão Institucional (IRB). Os dados do transcriptoma e as informações clínicas correspondentes para o LUSC foram baixados do banco de dados Cancer Genome Atlas (TCGA). Inicialmente, a coorte compreendia 551 amostras (502 tumorais e 49 normais). Para garantir a robustez desse modelo prognóstico, pacientes sem dados completos de sobrevivência geral (OS), status de sobrevivência ou características clínicas essenciais foram excluídos da análise de sobrevivência a jusante. Após esse rigoroso processo de filtragem, foi estabelecida uma coorte analítica final de 470 pacientes com LUSC. Esses 470 pacientes foram posteriormente divididos aleatoriamente em um conjunto de treinamento (n = 235) e um conjunto de testes (n = 235) para construção e validação do modelo de risco. Um total de 792 genes relacionados ao estresse oxidativo foram coletados por meio do banco de dados GeneCards, listado na Tabela de Materiais, utilizando um critério de inclusão com pontuação de relevância ≥ 7. O pacote limma (versão 3.66.0) foi utilizado para rastrear genes de estresse oxidativo diferencialmente expressos. Os plots de vulcão foram produzidos usando o pacote ggplot2 (versão 4.0.2). Entre o grupo de treinamento e o grupo de testes, não houve variação apreciável nas características clínicas (p > 0,05). O conjunto de treinamento era usado para construção de modelos, e o conjunto de testes era usado para validação interna.
Construção e verificação do modelo de risco
Para identificar lncRNAs associados ao estresse oxidativo, foi realizada uma análise de coexpressão de Pearson. Os autores avaliaram a correlação entre os perfis de expressão dos genes relacionados ao estresse oxidativo diferencialmente expressos identificados e todos os lncRNAs anotados dentro do conjunto de dados TCGA-LUSC. Os critérios de correlação para selecionar lncRNAs coexpressos foram definidos em um coeficiente de correlação |R| > 0,4 e um valor-p < 0,001. Com base nesses limiares, um total de 6.088 candidatos a lncRNAs relacionados ao estresse oxidativo foram derivados e extraídos. Posteriormente, esses 6.088 lncRNAs candidatos foram avaliados por meio de regressão LASSO e abordagens analíticas complementares para construir o modelo prognóstico de risco. Utilizando os dados processados, os participantes dos pacientes foram categorizados em camadas de alto e baixo risco.
Fatores independentes e curvas ROC
Modelos de regressão de Cox univariada e multivariada foram utilizados para identificar variáveis prognósticas independentes. Especificamente, para controlar a taxa de erro familiar durante as múltiplas análises univariadas de regressão de Cox (n = 236 testes), foi aplicada a correção de Bonferroni, estabelecendo o limiar de significância rigorosamente ajustado em p < 0,00021 (calculado como 0,05 / 236).
Análise de sobrevivência e análise de componentes principais
O pacote de sobrevivência (versão 3.8-6) foi usado para calcular as taxas de sobrevivência geral (OS), e a análise de componentes principais (PCA) foi usada para avaliar a robustez do modelo de risco.
Nomograma
O índice de ajuste foi avaliado usando o pacote rms em R para avaliar a confiabilidade do nomograma. A carga mutacional tumoral (TMB) foi avaliada e visualizada usando o pacote maftools (versão 2.22.0). Algoritmos, incluindo o ssGSEA, foram usados para detectar infiltração imunológica. Além disso, o pacote ggpubr (versão 0.6.3) foi usado para comparar pontos de verificação imunológico entre grupos de risco. O algoritmo ESTIMATE foi utilizado para caracterizar o microambiente tumoral (TME).
Análise funcional
O pacote clusterProfiler (versão 4.14.6) em R era usado para análise de enriquecimento. Para explorar mais a fundo as vias funcionais, foi realizada uma análise GSEA. Usando o Cytoscape, foi criada uma rede de coexpressões para observação.
Métodos experimentais
Condições e agrupamento da cultura celular
A linhagem celular de carcinoma escamoso pulmonar humano NCI-H520 (RRID: CVCL_1566; listada na Tabela de Materiais) foi cultivada em meio RPMI-1640 suplementado com 10% (v/v) soro fetal bovino e 1% de penicilina-estreptomicina. Células em fase logarítmica foram semeadas em placas de cultura de 6 poços. Após adesão bem-sucedida, as células foram alocadas aleatoriamente em três grupos experimentais distintos: um grupo controle em branco mantido sob condições padrão de cultura, um grupo controle negativo (si-NC) transfectado com siRNA não direcionado, e um grupo experimental de knockdown (si-LINC01615) transfectado com siRNA específico para LINC01615.
Transfecção celular
Os reagentes si-LINC01615 e si-NC foram descongelados no gelo. Para a preparação do complexo de transfecção, 95 μL de meio RPMI-1640 sem soro foram pipetados em tubos de centrífuga estéreis, seguidos pela adição sequencial de 3 μg do respectivo siRNA e 5 μL do reagente de transfecção Lipo3000. O procedimento idêntico foi aplicado tanto para o gruposi LINC01615 quanto para o siNC. As soluções foram suavemente misturadas e incubadas em temperatura ambiente por 5 minutos para facilitar a formação complexa, resultando em um volume total de aproximadamente 200 μL por tubo. Posteriormente, essa mistura foi distribuída uniformemente nos poços de cultivo designados contendo 800 μL de meio basal RPMI-1640. Após um período de incubação de 6 horas a 37 °C, o meio de transfecção foi cuidadosamente aspirado e substituído por meio de cultivo fresco e completo.
Ensaio de cicatrização de feridas (scratch)
Para o ensaio de cicatrização de feridas, células NCI-H520 (RRID: CVCL_1566) que crescem logarítmicamente foram semeadas em placas de cultura de 6 poços com uma densidade de 5,6 × 10 célulasde 5 por poço. Ao atingir 95–100% de confluência, uma ferida linear foi criada artificialmente ao riscar a monocamada celular verticalmente pelo centro do poço usando uma ponta estéril de pipeta de 200 μL. Os poços foram posteriormente lavados três vezes com PBS para remover células destacadas e detritos celulares, e depois incubados em meio RPMI-1640 sem soro. Imagens das áreas idênticas feridas foram capturadas às 0, 24 e 48 horas utilizando um microscópio invertido. A área de cicatrização da ferida foi quantificada medindo a distância da lacuna usando o ImageJ (RRID: SCR_003070; listada na Tabela de Materiais), e a taxa de migração foi calculada em relação à área inicial da ferida em 0 h.
Migração transwell e ensaios de invasão
As capacidades de migração e invasão celular foram avaliadas utilizando câmaras Transwell de 24 poços equipadas com inserções de membrana de policarbonato de 8,0 μm do tamanho de um poro listadas na Tabela de Materiais. Para o ensaio de invasão, as câmaras apicais foram pré-revestidas com 50 μL de Matrigel (diluído 1:8 em meio sem soro; listado na Tabela de Materiais) e incubadas a 37 °C por 2 horas para polimerização, enquanto o ensaio de migração utilizou inserts não revestidos. Após transfecção celular e um período de fome sérica de 24 horas, as células NCI-H520 foram colhidas, ressuspensas em meio RPMI-1640 sem soro e semeadas nas câmaras apicais com densidade de 1,8 × 105 células/poço em 200 μL. As câmaras basais foram preenchidas com 600 μL de meio RPMI-1640 suplementado com 10% de FBS como quimioatrativo. Após uma incubação de 48 horas a 37 °C com 5% de CO₂, os inserts foram removidos e lavados três vezes com PBS. As células remanescentes na superfície apical da membrana foram cuidadosamente removidas com um cotonete úmido. As células que migraram ou invadiram a superfície basolateral foram fixadas com paraformaldeído a 4% por 20 minutos e coradas com solução de Giemsa listada na Tabela de Materiais por 30 minutos em temperatura ambiente. As células migradas e invasas foram visualizadas e enumeradas em três campos ópticos selecionados aleatoriamente por inserto usando um microscópio invertido com ampliação de 40x.
Análise estatística
Todas as análises estatísticas e visualizações de dados foram realizadas usando o software R (versão 4.1.2; RRID: SCR_001905; listados na Tabela de Materiais). Os números de versão específicos de todos os pacotes R utilizados neste estudo foram explicitamente detalhados em suas respectivas subseções metodológicas. Um valor significativo era convencionalmente definido como p < 0,05, a menos que fosse informado de outra forma. Quando múltiplas comparações foram realizadas (por exemplo, as múltiplas regressões univariadas de Cox), o método de Bonferroni foi aplicado rigorosamente para manter a taxa de erro familiar (FWER), com o limiar ajustado exato (por exemplo, p < 0,00021) detalhado na respectiva subseção metodológica. Os dados quantitativos in vitro foram expressos como média ± desvio padrão e analisados usando o teste t de Student ou ANOVA unidirecional, conforme apropriado.