A expressão de FOXD1 está elevada em tecidos humanos de CCEP
Uma consulta ao Gene Expression Profiling Interactive Analysis (GEPIA) mostrou maior abundância de transcritos de FOXD1 no LUSC em comparação com amostras normais de pulmão, enquanto a comparação com adenocarcinoma de pulmão não alcançou significância estatística (Figura 1A). A análise de 12 pares de tecidos de LUSC e tecidos não tumorais adjacentes revelou a mesma tendência de alteração por qRT-PCR e por western blotting (Figura 1B,C). A coloração por imunohistoquímica também mostrou um sinal mais intenso de FOXD1 em seções tumorais do que em tecido adjacente pareado (Figura 1D).
A diminuição de FOXD1 reduz a viabilidade e a migração e aumenta a apoptose em células de carcinoma de pulmão de células escamosas (LUSC)
A abundância do transcrito FOXD1 foi maior nas células SK-MES-1 e NCI-H226 do que nas células BEAS-2B (Figura 2A). Três RNAs interferentes pequenos direcionados ao FOXD1 foram testados; o siRNA-1 produziu a maior redução detectada por qRT-PCR e por western blotting e foi selecionado para ensaios subsequentes (Figura 2B,C). Em relação ao controle negativo, o silenciamento do FOXD1 reduziu o sinal do Cell Counting Kit-8 (Figura 2D), aumentou o número de células positivas para Annexina V (Figura 2E) e diminuiu a migração em Transwell (Figura 2F) em ambas as linhagens celulares de LUSC. Um único reagente selecionado foi utilizado nas etapas seguintes, de modo que os resultados funcionais não constituem exclusão independente de sequência de efeitos fora do alvo.
A expressão de METTL3 correlaciona-se com FOXD1 em amostras de LUSC
A metilação m6A é uma modificação química dinâmica e reversível amplamente presente em eucariotos, regulada por metiltransferases e desmetilases, e reconhecida e processada por proteínas ligadoras de m6A18. O sinal global de m6A foi maior em 12 tecidos de LUSC do que em tecidos adjacentes pareados e também foi maior em ambas as linhagens celulares de LUSC do que nas células BEAS-2B (Figura 3A). Os transcritos que codificam METTL3, METTL14, WTAP e RBM15 estavam aumentados nas amostras tumorais pareadas (Figura 3B). Nos dados do GEPIA e no conjunto de 12 tumores, FOXD1 mostrou sua associação positiva mais clara com METTL3; a associação com RBM15 foi mais fraca (Figura 3C,D). A análise por Western blot confirmou ainda maior abundância de METTL3 nos tecidos tumorais (Figura 3E). Essas análises estabelecem co-variação, não regulação direta.
A diminuição de METTL3 altera a abundância de FOXD1 e os fenótipos das células de LUSC
Três RNAs interferentes pequenos direcionados ao METTL3 foram testados em células SK-MES-1 e NCI-H226 por qRT-PCR e western blotting (Figura 4A,B). A diminuição do METTL3 reduziu o sinal global de m6A (Figura 4C). A superexpressão de FOXD1 aumentou os níveis de FOXD1 sem elevar o METTL3 (Figura 4D), e o experimento de resgate com quatro grupos confirmou a menor abundância de FOXD1 após a diminuição do METTL3 e sua restauração pelo construto de FOXD1 (Figura 4E,F). A diminuição do METTL3 reduziu a viabilidade e a migração e aumentou a apoptose; a reexpressão de FOXD1 contrapôs parcialmente essas alterações (Figura 4G–I). Assim, os dados de resgate são consistentes com a atuação de FOXD1 downstream do METTL3 em células cultivadas.
A perturbação de METTL3 está associada à estabilidade do RNA e da proteína FOXD1
Evidências crescentes sugerem que a modificação de mRNA m6A modula diversos aspectos do ciclo do mRNA, incluindo o processamento de pré-mRNA, tradução, exportação nucleocitoplasmática, degradação e splicing19,20. O SRAMP identificou motivos candidatos de m6A no RNA de FOXD1, e o m6A2Target listou METTL3 como uma enzima escritora predita associada a FOXD1 em células H1299 (Figura 5A,B). Essas saídas de banco de dados têm caráter gerador de hipóteses e não demonstram ligação ou metilação em sítios específicos em células de LUSC. Após o bloqueio transcricional com actinomicina D, o RNA de FOXD1 diminuiu mais rapidamente após o silenciamento de METTL3 do que em células controle (Figura 5C). A abundância de FOXD1 foi maior em células com superexpressão de METTL3 tratadas com MG132 do que em controles vetoriais tratados com MG132, enquanto esse padrão não foi observado após tratamento com CHX (Figura 5D,E). O curso temporal com CHX mostrou perda mais rápida da proteína FOXD1 após o silenciamento de METTL3 (Figura 5F). Em conjunto, esses ensaios são compatíveis com alterações associadas a METTL3 na estabilidade do RNA e da proteína FOXD1, mas não identificam a adenosina modificada, não demonstram ligação direta nem definem um leitor de m6A.
A supressão do FOXD1 inibe o crescimento de xenografias de LUSC
Para testar a função da FOXD1 in vivo, células SK-MES-1 controle e com silenciamento de FOXD1 foram implantadas subcutaneamente em camundongos nu (n = 8 por grupo). A partir de uma semana após o implante, os tumores foram medidos a cada 4 dias durante 20 dias e coletados no dia 27 após a inoculação. Os tumores ressecados foram menores no grupo com silenciamento de FOXD1 (Figura 6A), suas curvas de crescimento aumentaram mais lentamente (Figura 6B) e seus pesos finais foram menores (Figura 6C). A coloração para FOXD1 e Ki-67 também foi reduzida nos tumores com silenciamento (Figura 6D). Esses experimentos avaliam apenas o silenciamento de FOXD1; eles não testam a inibição de METTL3 in vivo.
Em conjunto, FOXD1 foi superexpresso em tecidos e células de carcinoma de pulmão de células escamosas e promoveu a proliferação, migração e crescimento de xenografte. A diminuição de METTL3 reduziu os níveis globais de m6A, diminuiu a expressão de FOXD1 e prejudicou fenótipos celulares malignos, enquanto a reexpressão de FOXD1 reverteu parcialmente esses efeitos. Os ensaios de estabilidade apoiam a regulação pós-transcricional de FOXD1 associada a METTL3, mas não estabelecem metilação direta em sítios específicos ou o resgate in vivo de toda a via. O delineamento do estudo é resumido na Figura 7.
DISPONIBILIDADE DE DADOS:
Todos os dados analisados neste estudo estão incluídos no artigo e nas suas figuras. Os dados numéricos com identificação removida que sustentam os resumos relatados estão disponíveis mediante solicitação razoável ao autor correspondente, sujeitos às exigências institucionais de governança e de proteção aos participantes. Nenhum arquivo de simulação ou código personalizado de análise foi utilizado.

Figura 1: Expressão de FOXD1 em amostras de câncer de pulmão humano. (A) Comparação pelo GEPIA da abundância de transcritos de FOXD1 em adenocarcinoma de pulmão (LUAD) e carcinoma de células escamosas de pulmão (LUSC) versus amostras normais de pulmão. (B) Abundância de transcritos de FOXD1 em 12 pares de tecidos de LUSC e tecidos não tumorais adjacentes. (C) Abundância de proteína FOXD1 nos mesmos tecidos pareados. (D) Imunohistoquímica representativa de FOXD1. Número de pares de tecidos humanos = 12; repetições técnicas no estilo in vitro não alteram o valor de n ao nível do paciente. Barra de escala no painel D = 50 µm. ***P < 0,001. Clique aqui para visualizar uma versão maior desta figura.

Figura 2: Efeitos do silenciamento de FOXD1 em células SK-MES-1 e NCI-H226. (A) Abundância do transcrito de FOXD1 em células BEAS-2B, SK-MES-1 e NCI-H226. (B) Triagem por qRT-PCR de três RNAs interferentes pequenos direcionados a FOXD1, 48 h após a transfeção. (C) Triagem por Western blot dos mesmos três RNAs interferentes pequenos direcionados a FOXD1. (D) Ensaio com o kit de contagem celular-8 após o silenciamento de FOXD1. (E) Análise de apoptose com anexina V–isotiocianato de fluoresceína/iodeto de propídio após 48 h. (F) Migração em Transwell medida 24 h após a semeadura das células transfectadas; as imagens arquivadas do campo de visão são apresentadas sem barra de escala retrospectiva, pois os metadados de calibração estavam indisponíveis. Três experimentos biológicos independentes foram realizados. *P < 0,05; **P < 0,01; ***P < 0,001. Clique aqui para visualizar uma versão maior desta figura.

Figura 3: N6-metiladenosina global e a associação entre METTL3 e FOXD1. (A) Sinal global de m6A em tecidos pareados e nas células BEAS-2B, SK-MES-1 e NCI-H226. (B) Transcritos de METTL3, METTL14, WTAP e RBM15 em 12 pares de tecidos de CECP e tecidos adjacentes. (C) Correlações GEPIA entre FOXD1 e cada enzima metiltransferase. (D) Correlações entre FOXD1 e os transcritos das enzimas metiltransferases em 12 tumores de CECP. (E) Proteína METTL3 em tecidos de CECP pareados e tecidos adjacentes. Número de pares de tecidos humanos = 12; os ensaios celulares representam três experimentos biológicos independentes. *P < 0,05; **P < 0,01; ***P < 0,001. Clique aqui para visualizar uma versão ampliada desta figura.

Figura 4: Silenciamento de METTL3, resgate de FOXD1 e fenótipos de células de LUSC. (A) Triagem por qRT-PCR de três RNAs interferentes pequenos direcionados ao METTL3, 48 h após a transfeção. (B) Triagem por western blot dos mesmos três RNAs interferentes pequenos direcionados ao METTL3. (C) Sinal global de m6A após o silenciamento de METTL3. (D) Abundância do transcrito e da proteína FOXD1 após superexpressão de FOXD1. (E) Abundância dos transcritos de METTL3 e FOXD1 nas condições controle, si-METTL3, ov-FOXD1 e si-METTL3 + ov-FOXD1 combinadas, medida por qRT-PCR. (F) Abundância proteica de METTL3 e FOXD1 nos mesmos grupos de resgate, medida por western blot. (G) Ensaio com o Kit Cell Counting Kit-8. (H) Apoptose por citometria de fluxo. (I) Migração em Transwell. Três experimentos biológicos independentes foram realizados. Os números abaixo das bandas do western blot indicam as razões normalizadas entre a proteína-alvo e a GAPDH. ns, não significativo; *P < 0,05; **P < 0,01; ***P < 0,001. Clique aqui para visualizar uma versão ampliada desta figura.

Figura 5: METTL3 está associado à regulação do transcripto e da proteína FOXD1. (A) Sítios de N6-metiladenosina previstos pelo SRAMP dentro do transcripto FOXD1. (B) Previsão por m6A2Target identificando METTL3 como um escritor candidato associado ao FOXD1. (C) Abundância relativa de mRNA de FOXD1 após tratamento com actinomicina D a 5 µg/mL em células controle negativo e com redução de METTL3 nos tempos de 0 h, 2 h, 4 h e 6 h. (D) Abundância da proteína FOXD1 em células controle com vetor não tratadas, células controle com vetor tratadas com MG132 e células com superexpressão de METTL3 tratadas com MG132 (10 µM MG132). (E) Abundância da proteína FOXD1 em células controle com vetor não tratadas, células controle com vetor tratadas com CHX e células com superexpressão de METTL3 tratadas com CHX (20 µM CHX). (F) Expressão da proteína FOXD1 após tratamento com CHX a 20 µM por 0 h, 3 h, 6 h e 9 h em células controle negativo e com redução de METTL3. GAPDH foi utilizada como controle de carregamento. Os experimentos celulares foram realizados em três réplicas biológicas independentes. Os ensaios de blot mostrados são representativos; os números abaixo das bandas indicam as razões normalizadas de FOXD1/GAPDH. *P < 0,05; **P < 0,01; ***P < 0,001. Clique aqui para visualizar uma versão maior desta figura.

Figura 6: A supressão de FOXD1 inibe o crescimento de xenografias de SK-MES-1. (A) Tumores removidos dos grupos controle negativo e com supressão de FOXD1; barra de escala = 1 cm. (B) Crescimento tumoral medido a cada 4 dias durante 20 dias, iniciando 1 semana após a implantação, com coleta no dia 27 após a inoculação. (C) Peso tumoral na coleta. (D) Imunohistoquímica representativa de FOXD1 e Ki-67 em tecido de xenograftia; barra de escala = 50 µm. n = 8 camundongos por grupo; **P < 0,01. Clique aqui para visualizar uma versão maior desta figura.

Figura 7: Desenho geral do estudo e limites das evidências. Tecidos LUSC pareados (n = 12 pares) e linhagens celulares foram utilizados para perfis de FOXD1, METTL3 e m6A global; após perturbação genética e experimentos de resgate, foram realizados ensaios de viabilidade, apoptose, migração e estabilidade. Um modelo de xenotransplante com silenciamento de FOXD1 (n = 8 camundongos por grupo) foi monitorado por 20 dias. O diagrama distingue as evidências baseadas em células para METTL3 das evidências específicas de FOXD1 in vivo e identifica a ligação direta, a metilação m6A em sítios específicos, a identidade do leitor e a atividade de METTL3 in vivo como não testadas. Clique aqui para visualizar uma versão maior desta figura.
| Gene | Sentido | Sequência (5′–3′) |
| FOXD1 | Sentido | TCTGAGTAACGCGTTATGCTT |
| FOXD1 | Antissentido | CGGAATCGGACGCAAGTCTT |
| GAPDH | Sentido | GGAGCGAGATCCCTCCAAAAT |
| GAPDH | Antissentido | GGCTGTTGTCATACTTCTCATGG |
| METTL14 | Sentido | CTGAAAGTGCCGACAGCATTGG |
| METTL14 | Antissentido | CTCTCCTTCATCCAGATACTTACG |
| METTL3 | Sentido | CCAGCACAGCTTCAGCAGTTCC |
| METTL3 | Antissentido | GCGTGGAGATGGCAAGACAGATG |
| RBM15 | Sentido | GTTGTGGCTTATGTGGAGTTTAC |
| RBM15 | Antissentido | CACTTAAAACACCGGCATTGG |
| WTAP | Sentido | GCAACAACAGCAGGAGTCTGCA |
| WTAP | Antissentido | CTGCTGGACTTGCTTGAGGTAC |
Tabela 1: Sequências de iniciadores para reação de polimerase em cadeia com transcrição reversa quantitativa. As sequências são apresentadas na direção 5′–3′.