L'espressione di FOXD1 è aumentata nei tessuti umani di LUSC
Un'analisi interattiva del profilo di espressione genica (GEPIA) ha mostrato un'abbondanza maggiore di trascritto FOXD1 nel LUSC rispetto ai campioni di polmone normale, mentre il confronto con l'adenocarcinoma polmonare non ha raggiunto la significatività statistica (Figura 1A). L'analisi di 12 coppie di tessuti di LUSC e tessuti adiacenti non tumorali ha prodotto lo stesso andamento di variazione mediante qRT-PCR e immunoblotting (Figura 1B,C). La colorazione immunoistochimica ha inoltre evidenziato un segnale FOXD1 più intenso nelle sezioni tumorali rispetto al tessuto adiacente abbinato (Figura 1D).
La riduzione di FOXD1 diminuisce la vitalità e la migrazione e aumenta l'apoptosi nelle cellule di LUSC
L'abbondanza del trascritto FOXD1 era maggiore nelle cellule SK-MES-1 e NCI-H226 rispetto alle cellule BEAS-2B (Figura 2A). Sono stati analizzati tre piccoli RNA interferenti mirati a FOXD1; l'siRNA-1 ha prodotto la riduzione più ampia, rilevata mediante qRT-PCR e immunoblotting, ed è stato selezionato per i saggi successivi (Figura 2B,C). Rispetto al controllo negativo, la riduzione di FOXD1 ha diminuito il segnale del Cell Counting Kit-8 (Figura 2D), aumentato le cellule positive per Annexina V (Figura 2E) e ridotto la migrazione in Transwell (Figura 2F) in entrambe le linee cellulari di LUSC. È stato utilizzato un singolo reagente selezionato, pertanto i risultati funzionali non costituiscono un'esclusione indipendente dalla sequenza degli effetti fuori bersaglio.
L'espressione di METTL3 è correlata a FOXD1 nei campioni di LUSC
La metilazione m6A è una modifica chimica dinamica e reversibile ampiamente presente negli eucarioti, regolata da metiltransferasi e demetilasi, e riconosciuta ed elaborata da proteine leganti la m6A18. Il segnale globale di m6A era più elevato nei 12 tessuti di LUSC rispetto ai tessuti adiacenti abbinati ed era inoltre più alto in entrambe le linee cellulari di LUSC rispetto alle cellule BEAS-2B (Figura 3A). I trascritti che codificano per METTL3, METTL14, WTAP e RBM15 erano aumentati nei campioni tumorali abbinati (Figura 3B). Nei dati di GEPIA e nel set di 12 tumori, FOXD1 mostrava la sua associazione positiva più chiara con METTL3; l'associazione con RBM15 era più debole (Figura 3C,D). Il saggio di immunoblotting ha ulteriormente mostrato un'abbondanza maggiore di METTL3 nei tessuti tumorali (Figura 3E). Queste analisi indicano una variazione concomitante, non una regolazione diretta.
L'abbattimento di METTL3 modifica l'abbondanza di FOXD1 e i fenotipi delle cellule di LUSC
Sono stati analizzati tre RNA interferenti specifici per METTL3 in cellule SK-MES-1 e NCI-H226 mediante qRT-PCR e immunoblotting (Figura 4A,B). L'abbattimento di METTL3 ha ridotto il segnale globale di m6A (Figura 4C). L'overespressione di FOXD1 ha aumentato i livelli di FOXD1 senza incrementare quelli di METTL3 (Figura 4D), e l'esperimento di recupero con quattro gruppi ha confermato una minore abbondanza di FOXD1 dopo l'abbattimento di METTL3 e il suo ripristino mediante il costrutto per FOXD1 (Figura 4E,F). L'abbattimento di METTL3 ha ridotto la vitalità e la migrazione cellulare e ha aumentato l'apoptosi; la ri-espressione di FOXD1 ha in parte contrastato questi effetti (Figura 4G–I). Pertanto, i dati di recupero sono coerenti con un ruolo di FOXD1 a valle di METTL3 nelle cellule in coltura.
La perturbazione di METTL3 è associata alla stabilità dell'RNA e della proteina FOXD1
Crescenti evidenze indicano che la modificazione dell'mRNA mediante m6A regola vari aspetti del ciclo dell'mRNA, inclusi lo splicing dell'pre-mRNA, la traduzione, l'esportazione nucleocitoplasmatica, il degrado e lo splicing19,20. SRAMP ha identificato motivi candidati di m6A nell'RNA di FOXD1, mentre m6A2Target ha indicato METTL3 come scrittore predetto associato a FOXD1 nelle cellule H1299 (Figura 5A,B). Questi risultati ottenuti dai database hanno valore ipotetico-generativo e non dimostrano il legame o la metilazione sito-specifica nelle cellule di LUSC. Dopo il blocco trascrizionale con attinomicina D, l'RNA di FOXD1 è diminuito più rapidamente dopo il silenziamento di METTL3 rispetto alle cellule di controllo (Figura 5C). L'abbondanza di FOXD1 era maggiore nelle cellule con sovraespressione di METTL3 trattate con MG132 rispetto ai controlli con vettore trattati con MG132, mentre questo andamento non è stato osservato dopo il trattamento con CHX (Figura 5D,E). Il corso temporale con CHX ha mostrato una perdita più rapida della proteina FOXD1 dopo il silenziamento di METTL3 (Figura 5F). Nel complesso, questi saggi sono compatibili con variazioni della stabilità dell'RNA e della proteina FOXD1 associate a METTL3, ma non identificano l'adenosina modificata, non dimostrano un legame diretto né definiscono un lettore di m6A.
La riduzione di FOXD1 inibisce la crescita dell'impianto xenotrapiantato di LUSC
Per verificare la funzione di FOXD1 in vivo, cellule SK-MES-1 di controllo e con silenziamento di FOXD1 sono state impiantate per via sottocutanea in topi nudi (n = 8 per gruppo). A partire da una settimana dopo l'impianto, i tumori sono stati misurati ogni 4 giorni per 20 giorni e sono stati prelevati al giorno 27 dopo l'inoculo. I tumori asportati erano più piccoli nel gruppo con silenziamento di FOXD1 (Figura 6A), le loro curve di crescita aumentavano più lentamente (Figura 6B) e i loro pesi finali erano inferiori (Figura 6C). La colorazione per FOXD1 e Ki-67 era inoltre ridotta nei tumori con silenziamento (Figura 6D). Questi esperimenti valutano esclusivamente il silenziamento di FOXD1; non verificano l'inibizione di METTL3 in vivo.
Nel complesso, FOXD1 era sovraespresso nei tessuti e nelle cellule del carcinoma polmonare a cellule squamose e favoriva la proliferazione, la migrazione e la crescita dell'xenotrapianto. L'abbattimento di METTL3 ha ridotto i livelli globali di m6A, diminuito l'espressione di FOXD1 e compromesso i fenotipi cellulari maligni, mentre la ri-espressione di FOXD1 ha in parte invertito questi effetti. I saggi di stabilità supportano una regolazione post-trascrizionale di FOXD1 associata a METTL3, ma non dimostrano una metilazione diretta sito-specifica o il ripristino in vivo dell'intero percorso. Il disegno dello studio è riassunto nella Figura 7.
DISPONIBILITÀ DEI DATI:
Tutti i dati analizzati in questo studio sono inclusi nell'articolo e nelle relative figure. I dati numerici anonimi alla base dei riassunti riportati sono disponibili presso l'autore corrispondente su richiesta ragionevole, previa approvazione secondo le normative istituzionali e i requisiti di protezione dei partecipanti. Non sono stati utilizzati file di simulazione né codice personalizzato per l'analisi.

Figura 1: Espressione di FOXD1 in campioni di cancro polmonare umano. (A) Confronto tramite GEPIA dell'abbondanza dei trascritti di FOXD1 nel carcinoma polmonare a cellule adenocarcinoma (LUAD) e nel carcinoma polmonare a cellule squamose (LUSC) rispetto ai campioni normali di polmone. (B) Abbondanza dei trascritti di FOXD1 in 12 coppie di tessuti di LUSC e tessuti adiacenti non tumorali. (C) Abbondanza della proteina FOXD1 negli stessi tessuti accoppiati. (D) Esempio rappresentativo di immunohistochimica per FOXD1. Numero di coppie di tessuti umani = 12; le ripetizioni tecniche di tipo in vitro non modificano il numero a livello paziente. Barra della scala nel pannello D = 50 µm. ***P < 0,001. Cliccare qui per visualizzare una versione ingrandita di questa figura.

Figura 2: Effetti dell'abbattimento di FOXD1 nelle cellule SK-MES-1 e NCI-H226. (A) Abbondanza del trascritto di FOXD1 nelle cellule BEAS-2B, SK-MES-1 e NCI-H226. (B) Screening mediante qRT-PCR di tre piccoli RNA interferenti mirati a FOXD1, eseguito 48 h dopo la trasfezione. (C) Screening mediante Western blot degli stessi tre piccoli RNA interferenti mirati a FOXD1. (D) Saggio con il kit per il conteggio cellulare Cell Counting Kit-8 dopo l'abbattimento di FOXD1. (E) Analisi dell'apoptosi mediante Annexin V–fluorescein isothiocyanate/propidio ioduro a 48 h. (F) Migrazione in Transwell misurata 24 h dopo la semina delle cellule trasfettate; le immagini archiviate del campo visivo sono presentate senza barra di scala retrospettiva poiché i metadati di calibrazione non erano disponibili. Sono stati eseguiti tre esperimenti biologici indipendenti. *P < 0.05; **P < 0.01; ***P < 0.001. Cliccare qui per visualizzare una versione ingrandita di questa figura.

Figura 3: Metilazione globale N6-adenosina e associazione tra METTL3 e FOXD1. (A) Segnale globale di m6A in tessuti accoppiati e nelle cellule BEAS-2B, SK-MES-1 e NCI-H226. (B) Trascritti di METTL3, METTL14, WTAP e RBM15 in 12 coppie di tessuti di LUSC e tessuti adiacenti. (C) Correlazioni GEPIA tra FOXD1 e ciascun enzima metiltransferasi. (D) Correlazioni tra i trascritti di FOXD1 e quelli degli enzimi metiltransferasi in 12 tumori di LUSC. (E) Proteina METTL3 in tessuti di LUSC e tessuti adiacenti accoppiati. Numero di campioni di tessuto umano = 12 coppie; i saggi cellulari rappresentano tre esperimenti biologici indipendenti. *P < 0.05; **P < 0.01; ***P < 0.001. Cliccare qui per visualizzare una versione ingrandita di questa figura.

Figura 4: Silenziamento di METTL3, recupero di FOXD1 e fenotipi delle cellule LUSC. (A) Screening mediante qRT-PCR di tre piccoli RNA interferenti mirati a METTL3, eseguito 48 h dopo la trasfezione. (B) Screening mediante western blot degli stessi tre piccoli RNA interferenti mirati a METTL3. (C) Segnale globale di m6A dopo il silenziamento di METTL3. (D) Abbondanza del trascritto e della proteina FOXD1 dopo sovraespressione di FOXD1. (E) Abbondanza dei trascritti di METTL3 e FOXD1 nelle condizioni di controllo, si-METTL3, ov-FOXD1 e si-METTL3 + ov-FOXD1 combinate, misurata mediante qRT-PCR. (F) Abbondanza proteica di METTL3 e FOXD1 negli stessi gruppi di recupero, misurata mediante western blot. (G) Saggio con Cell Counting Kit-8. (H) Apoptosi mediante citometria a flusso. (I) Migrazione in Transwell. Sono stati eseguiti tre esperimenti biologici indipendenti. I numeri al di sotto delle bande del western blot indicano i rapporti normalizzati tra proteina bersaglio e GAPDH. ns, non significativo; *P < 0,05; **P < 0,01; ***P < 0,001. Cliccare qui per visualizzare una versione ingrandita di questa figura.

Figura 5: METTL3 è associato alla regolazione del trascritto e della proteina FOXD1. (A) Siti di N6-metiladenosina previsti da SRAMP all'interno del trascritto FOXD1. (B) Predizione tramite m6A2Target che identifica METTL3 come scrittore candidato associato a FOXD1. (C) Abbondanza relativa di mRNA di FOXD1 dopo trattamento con actinomicina D a 5 µg/mL in cellule con controllo negativo e in cellule con silenziamento di METTL3 a 0 h, 2 h, 4 h e 6 h. (D) Abbondanza della proteina FOXD1 in cellule con vettore di controllo non trattate, in cellule con vettore di controllo trattate con MG132 e in cellule che sovraesprimono METTL3 trattate con MG132 (10 µM MG132). (E) Abbondanza della proteina FOXD1 in cellule con vettore di controllo non trattate, in cellule con vettore di controllo trattate con CHX e in cellule che sovraesprimono METTL3 trattate con CHX (20 µM CHX). (F) Espressione della proteina FOXD1 dopo trattamento con CHX a 20 µM per 0 h, 3 h, 6 h e 9 h in cellule con controllo negativo e in cellule con silenziamento di METTL3. GAPDH è stata utilizzata come controllo di caricamento. Gli esperimenti cellulari sono stati eseguiti in tre replicati biologici indipendenti. I blot mostrati sono rappresentativi; i numeri sotto le bande indicano i rapporti normalizzati FOXD1/GAPDH. *P < 0,05; **P < 0,01; ***P < 0,001. Cliccare qui per visualizzare una versione ingrandita di questa figura.

Figura 6: L'abbattimento di FOXD1 inibisce la crescita dell'impianto eterotrasplantato di SK-MES-1. (A) Tumori escissi dai gruppi di controllo negativo e con abbattimento di FOXD1; barra della scala = 1 cm. (B) Crescita tumorale misurata ogni 4 giorni per 20 giorni, a partire da una settimana dopo l'innesto, con raccolta al giorno 27 dopo l'inoculazione. (C) Peso del tumore al momento del prelievo. (D) Immunohistochimica rappresentativa per FOXD1 e Ki-67 nei tessuti dell'impianto eterotrasplantato; barra della scala = 50 µm. n = 8 topi per gruppo; **P < 0,01. Cliccare qui per visualizzare una versione ingrandita di questa figura.

Figura 7: Progettazione complessiva dello studio e limiti delle evidenze. Tessuti di LUSC accoppiati (n = 12 coppie) e linee cellulari sono stati utilizzati per l'analisi di FOXD1, METTL3 e dell'm6A globale; dopo perturbazione genetica e recupero sono stati eseguiti saggi di vitalità, apoptosi, migrazione e stabilità. Un modello di xenotrapianto con silenziamento di FOXD1 (n = 8 topi per gruppo) è stato monitorato per 20 giorni. Il diagramma distingue le evidenze relative a METTL3 ottenute in cellule da quelle specifiche di FOXD1 ottenute in vivo e identifica il legame diretto, l'm6A sito-specifico, l'identità del lettore e l'attività di METTL3 in vivo come aspetti non testati. Cliccare qui per visualizzare una versione ingrandita di questa figura.
| Gene | Direzione | Sequenza (5′–3′) |
| FOXD1 | Forward | TCTGAGTAACGCGTTATGCTT |
| FOXD1 | Reverse | CGGAATCGGACGCAAGTCTT |
| GAPDH | Forward | GGAGCGAGATCCCTCCAAAAT |
| GAPDH | Reverse | GGCTGTTGTCATACTTCTCATGG |
| METTL14 | Forward | CTGAAAGTGCCGACAGCATTGG |
| METTL14 | Reverse | CTCTCCTTCATCCAGATACTTACG |
| METTL3 | Forward | CCAGCACAGCTTCAGCAGTTCC |
| METTL3 | Reverse | GCGTGGAGATGGCAAGACAGATG |
| RBM15 | Forward | GTTGTGGCTTATGTGGAGTTTAC |
| RBM15 | Reverse | CACTTAAAACACCGGCATTGG |
| WTAP | Forward | GCAACAACAGCAGGAGTCTGCA |
| WTAP | Reverse | CTGCTGGACTTGCTTGAGGTAC |
Tabella 1: Sequenze dei primer per la reazione a catena della polimerasi con trascrizione inversa quantitativa. Le sequenze sono riportate nella direzione 5′–3′.