FOXD1-expressie is verhoogd in menselijke LUSC-weefsels
Een Gene Expression Profiling Interactive Analysis (GEPIA) query toonde een hogere abundantie van FOXD1-transcripten in LUSC dan in normale longmonsters, terwijl de vergelijking met longadenocarcinoom geen significantie bereikte (Figuur 1A). Analyse van 12 gepaarde LUSC- en aangrenzende niet-tumorweefsels leverde via qRT-PCR en western blotting dezelfde trend in verandering op (Figuur 1B,C). Immunohistochemische kleuring toonde eveneens een sterker FOXD1-signaal in tumorsecties dan in het gepaarde aangrenzende weefsel (Figuur 1D).
FOXD1-knockdown vermindert de levensvatbaarheid en migratie en verhoogt apoptose in LUSC-cellen
De abundantie van FOXD1-transcripten was groter in SK-MES-1- en NCI-H226-cellen dan in BEAS-2B-cellen (Figuur 2A). Er werden drie small interfering RNA's gericht tegen FOXD1 gescreend; siRNA-1 veroorzaakte de grootste reductie gedetecteerd via qRT-PCR en western blotting en werd geselecteerd voor daaropvolgende assays (Figuur 2B,C). Ten opzichte van de negatieve controle verminderde FOXD1-knockdown het Cell Counting Kit-8-signaal (Figuur 2D), verhoogde het aantal Annexine V-positieve cellen (Figuur 2E) en verlaagde het de Transwell-migratie (Figuur 2F) in beide LUSC-cellijnen. Er is één geselecteerde reagent gebruikt, waardoor de functionele resultaten geen sequentie-onafhankelijke uitsluiting van off-target effecten vormen.
METTL3-expressie correleert met FOXD1 in LUSC-monsters
m6A-methylering is een dynamische en reversibele chemische modificatie die veel voorkomt in eukaryoten, wordt gereguleerd door methyltransferases en demethylases, en wordt herkend en verwerkt door m6A-bindende eiwitten18. Het globale m6A-signaal was hoger in de 12 LUSC-weefsels dan in de gepaarde aangrenzende weefsels en was ook hoger in beide LUSC-cellijnen dan in BEAS-2B-cellen (Figuur 3A). Transcripten die coderen voor METTL3, METTL14, WTAP en RBM15 waren verhoogd in de gepaarde tumorstalen (Figuur 3B). In zowel de GEPIA-data als de set van 12 tumoren vertoonde FOXD1 de duidelijkste positieve associatie met METTL3; de associatie met RBM15 was zwakker (Figuur 3C,D). Western blotting toonde verder een hogere METTL3-abundantie aan in tumorweefsels (Figuur 3E). Deze analyses stellen co-variatie vast, geen directe regulatie.
METTL3-knockdown verandert de abundantie van FOXD1 en de fenotypen van LUSC-cellen
Drie small interfering RNA's die gericht zijn op METTL3 werden gescreend in SK-MES-1 en NCI-H226 cellen via qRT-PCR en western blotting (Figuur 4A,B). METTL3 knockdown verlaagde het globale m6A-signaal (Figuur 4C). Overexpressie van FOXD1 verhoogde FOXD1 zonder METTL3 te verhogen (Figuur 4D), en het rescue-experiment met vier groepen bevestigde een lagere abundantie van FOXD1 na METTL3 knockdown en de herstel hiervan door het FOXD1-construct (Figuur 4E,F). METTL3 knockdown verminderde de viabiliteit en migratie en verhoogde de apoptose; re-expressie van FOXD1 maakte deze veranderingen gedeeltelijk ongedaan (Figuur 4G–I). De rescue-data zijn hiermee consistent met de stelling dat FOXD1 stroomafwaarts van METTL3 werkt in gekweekte cellen.
Perturbatie van METTL3 is geassocieerd met de stabiliteit van FOXD1 RNA en eiwit
Toenemend bewijs suggereert dat m6A mRNA-modificatie verschillende aspecten van de mRNA-cyclus moduleert, waaronder pre-mRNA splicing, translatie, nucleocytoplasmatische export, afbraak en splicing19,20. SRAMP identificeerde kandidaat m6A-motieven in FOXD1 RNA, en m6A2Target vermeldde METTL3 als een voorspelde FOXD1-geassocieerde writer in H1299-cellen (Figuur 5A,B). Deze database-outputs zijn hypothese-genererend en tonen geen binding of site-specifieke methylering aan in LUSC-cellen. Na transcriptionele blokkade met actinomycine D nam FOXD1 RNA sneller af na METTL3-knockdown dan in controlecellen (Figuur 5C). De abundantie van FOXD1 was hoger in met MG132 behandelde METTL3-overexpressiecellen dan in met MG132 behandelde vectorcontroles, terwijl dit patroon niet werd waargenomen na CHX-behandeling (Figuur 5D,E). De CHX-tijdreeks liet een sneller verlies van FOXD1-eiwit zien na METTL3-knockdown (Figuur 5F). Collectief zijn deze assays consistent met METTL3-geassocieerde veranderingen in de stabiliteit van FOXD1 RNA en -eiwit, maar ze identificeren de gemodificeerde adenosine niet, tonen geen directe binding aan en definiëren geen m6A-reader.
FOXD1 knockdown onderdrukt de groei van LUSC-xenotransplantaten
Om de functie van FOXD1 in vivo te testen, werden controle- en FOXD1-knockdown SK-MES-1 cellen subcutaan geïmplanteerd in naaktmuizen (n = 8 per groep). Vanaf 1 week na implantatie werden de tumoren elke 4 dagen gedurende 20 dagen gemeten en werden ze verzameld op dag 27 na inoculatie. De uitgesneden tumoren waren kleiner in de FOXD1-knockdown groep (Figuur 6A), hun groeicurves stegen langzamer (Figuur 6B) en hun uiteindelijke gewichten waren lager (Figuur 6C). De FOXD1- en Ki-67-kleuring was ook verminderd in de knockdown-tumoren (Figuur 6D). Deze experimenten evalueren uitsluitend de knockdown van FOXD1; ze testen de METTL3-inhibitie in vivo niet.
Samengevat was FOXD1 overexpressie aanwezig in weefsels en cellen van plaveiselcelcarcinoom van de long en ondersteunde het proliferatie, migratie en xenograft-groei. METTL3-knockdown verminderde de globale m6A-niveaus, verlaagde de FOXD1-expressie en verslechterde maligne cellulaire fenotypes, terwijl re-expressie van FOXD1 deze effecten gedeeltelijk terugdraaide. De stabiliteitsassays ondersteunen METTL3-geassocieerde post-transcriptionele regulatie van FOXD1, maar stellen geen directe site-specifieke methylering of in vivo rescue van het volledige pad vast. Het onderzoeksontwerp is samengevat in Figuur 7.
DATABESCHIKBAARHEID:
Alle in deze studie geanalyseerde gegevens zijn opgenomen in het artikel en de bijbehorende figuren. Geanonimiseerde numerieke gegevens die ten grondslag liggen aan de gerapporteerde samenvattingen zijn op redelijk verzoek verkrijgbaar bij de corresponderende auteur, onderhevig aan institutionele regelgeving en vereisten voor de bescherming van deelnemers. Er zijn geen simulatiebestanden of aangepaste analysecodes gebruikt.

Figuur 1: FOXD1-expressie in menselijke longkankermonsters. (A) GEPIA-vergelijking van de abundantie van FOXD1-transcripten in longadenocarcinoom (LUAD) en longplaveiselcelcarcinoom (LUSC) versus normale longmonsters. (B) Abundantie van FOXD1-transcripten in 12 gepaarde LUSC- en aangrenzende niet-tumorweefsels. (C) Abundantie van FOXD1-eiwit in dezelfde gepaarde weefsels. (D) Representatieve FOXD1-immunohistochemie. Menselijk weefsel n = 12 paren; in vitro-stijl technische herhaling verandert de n op patiëntniveau niet. Schaalbalk in paneel D = 50 µm. ***P < 0,001. Klik hier om een grotere versie van deze figuur te bekijken.

Figuur 2: Effecten van FOXD1 knockdown in SK-MES-1 en NCI-H226 cellen. (A) FOXD1 transcript- abundantie in BEAS-2B, SK-MES-1 en NCI-H226 cellen. (B) qRT-PCR screening van drie FOXD1-gerichte small interfering RNA's 48 u na transfectie. (C) Western-blot screening van dezelfde drie FOXD1-gerichte small interfering RNA's. (D) Cell counting kit-8 assay na FOXD1 knockdown. (E) Annexine V–fluoresceïne-isothiocyanaat/propidiumjodide apoptose-analyse na 48 u. (F) Transwell-migratie gemeten 24 u na het uitplaten van getransfecteerde cellen; de gearchiveerde beelden van het gezichtsveld worden gepresenteerd zonder retrospectieve schaalbalk omdat kalibratiemetadata niet beschikbaar waren. Drie onafhankelijke biologische experimenten zijn uitgevoerd. *P < 0.05; **P < 0.01; ***P < 0.001. Klik hier om een grotere versie van deze figuur te bekijken.

Figuur 3: Globale N6-methyladenosine en de associatie tussen METTL3 en FOXD1. (A) Globaal m6A-signaal in gepaarde weefsels en in BEAS-2B, SK-MES-1 en NCI-H226-cellen. (B) METTL3-, METTL14-, WTAP- en RBM15-transcripten in 12 gepaarde LUSC- en aangrenzende weefsels. (C) GEPIA-correlaties tussen FOXD1 en elke writer. (D) Correlaties tussen FOXD1 en writer-transcripten in 12 LUSC-tumoren. (E) METTL3-eiwit in gepaarde LUSC- en aangrenzende weefsels. Menselijk weefsel n = 12 paren; celassays vertegenwoordigen drie onafhankelijke biologische experimenten. *P < 0.05; **P < 0.01; ***P < 0.001. Klik hier om een grotere versie van deze figuur te bekijken.

Figuur 4: METTL3-knockdown, FOXD1-rescue en LUSC-cel-fenotypes. (A) qRT-PCR screening van drie kleine interfererende RNA's gericht tegen METTL3, 48 u na transfectie. (B) Western-blot screening van dezelfde drie kleine interfererende RNA's gericht tegen METTL3. (C) Globaal m6A-signaal na METTL3-knockdown. (D) FOXD1-transcript- en eiwitovervloed na FOXD1-overexpressie. (E) METTL3- en FOXD1-transcriptovervloed in controle-, si-METTL3-, ov-FOXD1- en gecombineerde si-METTL3 + ov-FOXD1-condities, gemeten via qRT-PCR. (F) METTL3- en FOXD1-eiwitovervloed in dezelfde rescue-groepen, gemeten via western blotting. (G) Cell Counting Kit-8 assay. (H) Apoptose via flowcytometrie. (I) Transwell-migratie. Er zijn drie onafhankelijke biologische experimenten uitgevoerd. Getallen onder de western-blot banden geven genormaliseerde target-eiwit/GAPDH-ratio's aan. ns, niet significant; *P < 0.05; **P < 0.01; ***P < 0.001. Klik hier om een grotere versie van deze figuur te bekijken.

Figuur 5: METTL3 is geassocieerd met de regulatie van FOXD1-transcript en -eiwit. (A) Door SRAMP voorspelde N6-methyladenosine-sites binnen het FOXD1-transcript. (B) m6A2Target-voorspelling waarbij METTL3 is geïdentificeerd als een kandidaat-writer geassocieerd met FOXD1. (C) Relatieve FOXD1-mRNA-overvloed na behandeling met actinomycine D (5 µg/mL) in negatieve-controlecellen en METTL3-knockdowncellen op 0 h, 2 h, 4 h en 6 h. (D) FOXD1-eiwitovervloed in onbehandelde vector-controlecellen, met MG132 behandelde vector-controlecellen en met MG132 behandelde METTL3-overexpressiecellen (10 µM MG132). (E) FOXD1-eiwitovervloed in onbehandelde vector-controlecellen, met CHX behandelde vector-controlecellen en met CHX behandelde METTL3-overexpressiecellen (20 µM CHX). (F) FOXD1-eiwitexpressie na behandeling met CHX (20 µM) gedurende 0 h, 3 h, 6 h en 9 h in negatieve-controlecellen en METTL3-knockdowncellen. GAPDH diende als ladingscontrole. Celexperimenten werden uitgevoerd in drie onafhankelijke biologische replica's. De getoonde blots zijn representatief; getallen onder de banden geven genormaliseerde FOXD1/GAPDH-ratio's aan. *P < 0,05; **P < 0,01; ***P < 0,001. Klik hier om een grotere versie van deze figuur te bekijken.

Figuur 6: FOXD1-knockdown onderdrukt de groei van SK-MES-1 xenografts. (A) Uitgesneden tumoren van de negatieve-controlegroep en de FOXD1-knockdowngroep; schaalbalk = 1 cm. (B) Tumorgroei gemeten elke 4 dagen gedurende 20 dagen, beginnend 1 week na implantatie, met verzameling op dag 27 na inoculatie. (C) Tumorgewicht bij verzameling. (D) Representatieve immunohistochemie voor FOXD1 en Ki-67 in xenograft-weefsel; schaalbalk = 50 µm. n = 8 muizen per groep; **P < 0.01. Klik hier om een grotere versie van deze figuur te bekijken.

Figuur 7: Algemeen onderzoeksontwerp en bewijsgrenzen. Gepaarde LUSC-weefsels (n = 12 paren) en cellijnen werden gebruikt voor de profilering van FOXD1, METTL3 en globale m6A; genetische perturbatie en rescue werden gevolgd door assays voor levensvatbaarheid, apoptose, migratie en stabiliteit. Een FOXD1-knockdown xenograft-model (n = 8 muizen per groep) werd gedurende 20 dagen gemonitord. Het diagram maakt onderscheid tussen celgebaseerd METTL3-bewijs en FOXD1-specifiek in vivo bewijs en identificeert directe binding, site-specifieke m6A, reader-identiteit en in vivo METTL3-activiteit als niet getest. Klik hier om een grotere versie van deze figuur te bekijken.
| Gen | Richting | Sequence (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 |
Tabel 1: Primersequenties voor kwantitatieve reverse-transcriptie polymerasekettingreactie. Sequenties worden vermeld in de 5′–3′ richting.