La expresión de FOXD1 está elevada en tejidos humanos de CME
Una consulta de análisis interactivo de perfiles de expresión génica (GEPIA) mostró una mayor abundancia de transcritos de FOXD1 en LUSC en comparación con muestras pulmonares normales, mientras que la comparación con adenocarcinoma de pulmón no alcanzó significancia (Figura 1A). El análisis de 12 tejidos emparejados de LUSC y tejidos adyacentes no tumorales mostró la misma dirección del cambio mediante qRT-PCR e inmunotransferencia (Figura 1B,C). La tinción inmunohistoquímica también mostró una señal de FOXD1 más intensa en secciones tumorales que en tejido adyacente emparejado (Figura 1D).
La inhibición de FOXD1 reduce la viabilidad y la migración, y aumenta la apoptosis en células de LUSC
La abundancia del ARN mensajero de FOXD1 fue mayor en las células SK-MES-1 y NCI-H226 que en las células BEAS-2B (Figura 2A). Se evaluaron tres ARN de interferencia pequeños dirigidos contra FOXD1; el siARN-1 produjo la mayor reducción detectada mediante qRT-PCR y western blotting, y se seleccionó para los ensayos posteriores (Figura 2B,C). En comparación con el control negativo, el silenciamiento de FOXD1 redujo la señal del kit de conteo celular CCK-8 (Figura 2D), aumentó el número de células positivas para Annexina V (Figura 2E) y disminuyó la migración en Transwell (Figura 2F) en ambas líneas celulares de carcinoma de células escamosas pulmonares. Se utilizó un único reactivo seleccionado, por lo que los resultados funcionales no constituyen una exclusión independiente de la secuencia respecto a efectos fuera de blanco.
La expresión de METTL3 se correlaciona con FOXD1 en muestras de CUP
La metilación m6A es una modificación química dinámica y reversible ampliamente presente en eucariotas, que es regulada por metiltransferasas y desmetilasas, y es reconocida y procesada por proteínas que se unen a m6A18. La señal global de m6A fue mayor en los 12 tejidos de LUSC que en los tejidos adyacentes apareados, y también fue mayor en ambas líneas celulares de LUSC que en las células BEAS-2B (Figura 3A). Los transcritos que codifican METTL3, METTL14, WTAP y RBM15 aumentaron en las muestras tumorales apareadas (Figura 3B). Tanto en los datos de GEPIA como en el conjunto de 12 tumores, FOXD1 mostró su asociación positiva más clara con METTL3; la asociación con RBM15 fue más débil (Figura 3C,D). El análisis por inmunotransferencia mostró además una mayor abundancia de METTL3 en los tejidos tumorales (Figura 3E). Estos análisis establecen una co-variación, no una regulación directa.
La reducción de METTL3 altera la abundancia de FOXD1 y los fenotipos de células de CUSP
Se ensayaron tres ARN de interferencia pequeños dirigidos contra METTL3 en células SK-MES-1 y NCI-H226 mediante qRT-PCR y transferencia western (Figura 4A,B). La reducción de METTL3 disminuyó la señal global de m6A (Figura 4C). La sobreexpresión de FOXD1 incrementó FOXD1 sin aumentar METTL3 (Figura 4D), y el experimento de rescate de cuatro grupos confirmó una menor abundancia de FOXD1 tras la reducción de METTL3 y su restauración mediante el constructo de FOXD1 (Figura 4E,F). La reducción de METTL3 disminuyó la viabilidad y la migración, y aumentó la apoptosis; la reexpresión de FOXD1 se opuso parcialmente a estos cambios (Figura 4G–I). Por lo tanto, los datos de rescate son consistentes con un papel de FOXD1 actuando aguas abajo de METTL3 en células en cultivo.
La perturbación de METTL3 se asocia con la estabilidad del ARN y de la proteína FOXD1
Creciente evidencia sugiere que la modificación de ARNm por m6A regula diversos aspectos del ciclo del ARNm, incluyendo el empalme de pre-ARNm, la traducción, la exportación nucleocitoplasmática, la degradación y el empalme19,20. SRAMP identificó motivos candidatos de m6A en el ARN de FOXD1, y m6A2Target listó a METTL3 como una enzima escritora predicha asociada a FOXD1 en células H1299 (Figura 5A,B). Estas salidas de bases de datos generan hipótesis, pero no demuestran unión ni metilación específica en sitios en células de LUSC. Tras el bloqueo transcripcional con actinomicina D, el ARN de FOXD1 disminuyó más rápidamente tras el silenciamiento de METTL3 que en células control (Figura 5C). La abundancia de FOXD1 fue mayor en células que sobreexpresan METTL3 tratadas con MG132 que en controles vectoriales tratados con MG132, mientras que este patrón no se observó tras el tratamiento con CHX (Figura 5D,E). El curso temporal con CHX mostró una pérdida más rápida de la proteína FOXD1 tras el silenciamiento de METTL3 (Figura 5F). En conjunto, estos ensayos son consistentes con cambios asociados a METTL3 en la estabilidad del ARN y la proteína FOXD1, pero no identifican la adenosina modificada, no demuestran unión directa ni definen un lector de m6A.
La inhibición de FOXD1 suprime el crecimiento del xenoinjerto de LUSC
Para evaluar la función de FOXD1 in vivo, se implantaron células control y células SK-MES-1 con silenciamiento de FOXD1 subcutáneamente en ratones desnudos (n = 8 por grupo). A partir de la primera semana tras el implante, los tumores se midieron cada 4 días durante 20 días y se recolectaron el día 27 después de la inoculación. Los tumores extirpados fueron más pequeños en el grupo con silenciamiento de FOXD1 (Figura 6A), sus curvas de crecimiento aumentaron más lentamente (Figura 6B) y sus pesos finales fueron menores (Figura 6C). La tinción de FOXD1 y Ki-67 también se redujo en los tumores con silenciamiento (Figura 6D). Estos experimentos evalúan únicamente el silenciamiento de FOXD1; no evalúan la inhibición de METTL3 in vivo.
En conjunto, FOXD1 se sobreexpresó en tejidos y células de carcinoma de células escamosas de pulmón y promovió la proliferación, migración y crecimiento de xenoinjertos. La inhibición de METTL3 redujo los niveles globales de m6A, disminuyó la expresión de FOXD1 y alteró los fenotipos celulares malignos, mientras que la reexpresión de FOXD1 revirtió parcialmente estos efectos. Los ensayos de estabilidad respaldan la regulación postranscripcional de FOXD1 asociada a METTL3, pero no establecen una metilación directa y específica en sitios determinados ni la recuperación in vivo de la vía completa. El diseño del estudio se resume en la Figura 7.
DISPONIBILIDAD DE LOS DATOS:
Todos los datos analizados en este estudio se incluyen en el artículo y sus figuras. Los datos numéricos desidentificados que sustentan los resúmenes informados están disponibles a solicitud razonable del autor correspondiente, sujeto a los requisitos institucionales de gobernanza y de protección de los participantes. No se utilizaron archivos de simulación ni código personalizado de análisis.

Figura 1: Expresión de FOXD1 en muestras de cáncer de pulmón humano. (A) Comparación mediante GEPIA de la abundancia de transcritos de FOXD1 en adenocarcinoma de pulmón (LUAD) y carcinoma de células escamosas de pulmón (LUSC) frente a muestras normales de pulmón. (B) Abundancia de transcritos de FOXD1 en 12 tejidos emparejados de LUSC y tejidos adyacentes no tumorales. (C) Abundancia de proteína FOXD1 en los mismos tejidos emparejados. (D) Inmunohistoquímica representativa de FOXD1. Tejido humano n = 12 pares; la repetición técnica al estilo in vitro no cambia el valor de n a nivel de paciente. Barra de escala en el panel D = 50 µm. ***P < 0,001. Haga clic aquí para ver una versión más grande de esta figura.

Figura 2: Efectos del silenciamiento de FOXD1 en células SK-MES-1 y NCI-H226. (A) Abundancia del transcrito de FOXD1 en células BEAS-2B, SK-MES-1 y NCI-H226. (B) Análisis por qRT-PCR de tres ARN interferentes pequeños dirigidos contra FOXD1 a las 48 h posteriores a la transfección. (C) Análisis por inmunotransferencia de los mismos tres ARN interferentes pequeños dirigidos contra FOXD1. (D) Ensayo con el kit de conteo celular-8 tras el silenciamiento de FOXD1. (E) Análisis de apoptosis mediante V-FITC/ioduro de propidio a las 48 h. (F) Migración en Transwell medida a las 24 h posteriores a la siembra de las células transfectadas; las imágenes archivadas del campo de visión se presentan sin barra de escala retrospectiva debido a la indisponibilidad de metadatos de calibración. Se realizaron tres experimentos biológicos independientes. *P < 0,05; **P < 0,01; ***P < 0,001. Haga clic aquí para ver una versión más grande de esta figura.

Figura 3: N6-metiladenosina global y la asociación entre METTL3 y FOXD1. (A) Señal global de m6A en tejidos apareados y en las células BEAS-2B, SK-MES-1 y NCI-H226. (B) Transcritos de METTL3, METTL14, WTAP y RBM15 en 12 pares de tejidos de carcinoma escamoso de pulmón (LUSC) y tejidos adyacentes. (C) Correlaciones GEPIA entre FOXD1 y cada enzima metiltransferasa. (D) Correlaciones entre FOXD1 y los transcritos de las enzimas metiltransferasas en 12 tumores de LUSC. (E) Proteína METTL3 en tejidos de LUSC y tejidos adyacentes apareados. Número de pares de tejidos humanos: n = 12; los ensayos en células representan tres experimentos biológicos independientes. *P < 0,05; **P < 0,01; ***P < 0,001. Haga clic aquí para ver una versión más grande de esta figura.

Figura 4: Silenciamiento de METTL3, rescate de FOXD1 y fenotipos de células de LUSC. (A) Análisis por qRT-PCR de tres ARN interferentes pequeños dirigidos contra METTL3 a las 48 h tras la transfección. (B) Análisis por inmunotransferencia de los mismos tres ARN interferentes pequeños dirigidos contra METTL3. (C) Señal global de m6A tras el silenciamiento de METTL3. (D) Abundancia del ARN mensajero y de la proteína FOXD1 tras la sobreexpresión de FOXD1. (E) Abundancia de los transcritos de METTL3 y FOXD1 en condiciones de control, si-METTL3, ov-FOXD1 y si-METTL3 + ov-FOXD1 combinadas, medida por qRT-PCR. (F) Abundancia proteica de METTL3 y FOXD1 en los mismos grupos de rescate medida por inmunotransferencia. (G) Ensayo con el kit Cell Counting Kit-8. (H) Apoptosis mediante citometría de flujo. (I) Migración en Transwell. Se realizaron tres experimentos biológicos independientes. Los números bajo las bandas de inmunotransferencia indican las relaciones normalizadas de proteína diana/GAPDH. ns, no significativo; *P < 0,05; **P < 0,01; ***P < 0,001. Haga clic aquí para ver una versión ampliada de esta figura.

Figura 5: METTL3 está asociado con la regulación del ARNm y la proteína FOXD1. (A) Sitios predichos por SRAMP de N6-metiladenosina dentro del ARNm de FOXD1. (B) Predicción mediante m6A2Target que identifica a METTL3 como un escritor candidato asociado con FOXD1. (C) Abundancia relativa de ARNm de FOXD1 tras el tratamiento con actinomicina D a 5 µg/mL en células con control negativo y células con reducción de METTL3 a las 0 h, 2 h, 4 h y 6 h. (D) Abundancia de la proteína FOXD1 en células con vector de control sin tratar, células con vector de control tratadas con MG132 y células que sobreexpresan METTL3 tratadas con MG132 (10 µM MG132). (E) Abundancia de la proteína FOXD1 en células con vector de control sin tratar, células con vector de control tratadas con CHX y células que sobreexpresan METTL3 tratadas con CHX (20 µM CHX). (F) Expresión de la proteína FOXD1 tras el tratamiento con CHX a 20 µM durante 0 h, 3 h, 6 h y 9 h en células con control negativo y células con reducción de METTL3. GAPDH sirvió como control de carga. Los experimentos en células se realizaron en tres réplicas biológicas independientes. Las transferencias mostradas son representativas; los números debajo de las bandas indican las relaciones normalizadas de FOXD1/GAPDH. *P < 0,05; **P < 0,01; ***P < 0,001. Haga clic aquí para ver una versión ampliada de esta figura.

Figura 6: La inhibición de FOXD1 suprime el crecimiento del xenoinjerto SK-MES-1. (A) Tumores extirpados de los grupos control negativo y con inhibición de FOXD1; barra de escala = 1 cm. (B) Crecimiento tumoral medido cada 4 días durante 20 días, comenzando una semana después del implante, con recolección al día 27 tras la inoculación. (C) Peso tumoral en el momento de la recolección. (D) Inmunohistoquímica representativa de FOXD1 y Ki-67 en tejido de xenoinjerto; barra de escala = 50 µm. n = 8 ratones por grupo; **P < 0,01. Haga clic aquí para ver una versión ampliada de esta figura.

Figura 7: Diseño general del estudio y límites de la evidencia. Se utilizaron tejidos de CUSC pareados (n = 12 pares) y líneas celulares para el perfilado de FOXD1, METTL3 y m6A global; después de la perturbación genética y el rescate se realizaron ensayos de viabilidad, apoptosis, migración y estabilidad. Se monitoreó durante 20 días un modelo de xenoinjerto con silenciamiento de FOXD1 (n = 8 ratones por grupo). El diagrama distingue la evidencia basada en células para METTL3 de la evidencia específica de FOXD1 in vivo e identifica la unión directa, el m6A sitio-específico, la identidad del lector y la actividad de METTL3 in vivo como aspectos no probados. Haga clic aquí para ver una versión más grande de esta figura.
| Gen | Sentido | Secuencia (5′–3′) |
| FOXD1 | Directo | TCTGAGTAACGCGTTATGCTT |
| FOXD1 | Reverso | CGGAATCGGACGCAAGTCTT |
| GAPDH | Directo | GGAGCGAGATCCCTCCAAAAT |
| GAPDH | Reverso | GGCTGTTGTCATACTTCTCATGG |
| METTL14 | Directo | CTGAAAGTGCCGACAGCATTGG |
| METTL14 | Reverso | CTCTCCTTCATCCAGATACTTACG |
| METTL3 | Directo | CCAGCACAGCTTCAGCAGTTCC |
| METTL3 | Reverso | GCGTGGAGATGGCAAGACAGATG |
| RBM15 | Directo | GTTGTGGCTTATGTGGAGTTTAC |
| RBM15 | Reverso | CACTTAAAACACCGGCATTGG |
| WTAP | Directo | GCAACAACAGCAGGAGTCTGCA |
| WTAP | Reverso | CTGCTGGACTTGCTTGAGGTAC |
Tabla 1: Secuencias de cebadores para la reacción en cadena de la polimerasa con transcripción inversa cuantitativa. Las secuencias se indican en la dirección 5′–3′.