Die Studie wurde von der Ethikkommission des Zweiten Angeschlossenen Krankenhauses der Nanchang-Universität (CDEFYYLK 3-05) genehmigt. Von allen Patienten wurde vor der Probenentnahme schriftlich informierte Einwilligung gemäß der Deklaration von Helsinki eingeholt. Die Rohdaten wurden mit Hilfe der in der Tabelle der Materialien aufgelisteten Software und Online-Ressourcen exportiert und analysiert.
Klinische Gewebesammlung
Gepaarte Tumor- und entsprechende entfernt liegende gesunde Lungengewebe wurden von Patienten gesammelt, die sich zwischen Januar und Juni 2024 einer chirurgischen Resektion unterzogen hatten und bei denen pathologisch ein invasives Lungenadenokarzinom diagnostiziert wurde. Die Proben wurden vom Operateur entfernt, innerhalb von 5 Minuten nach der Resektion aufgearbeitet, sofort in flüssigem Stickstoff schockgefroren und bis zur Analyse gelagert. Anonymisierte klinisch-pathologische Merkmale sind in Zusatz-Tabelle 1 angegeben.
Innere extraktive Elektrospray-Ionisation-Massenspektrometrie (iEESI-MS) zur Unterscheidung von LUAD und normalem Lungengewebe
Gewebeproben (1 mm3) wurden mittels iEESI-MS unter Verwendung eines hybriden linearen Ionenfallen-elektrostatischen Orbitrap-Massenspektrometers im positiven Ionenmodus analysiert. Als Extraktionslösung diente Methanol (100 %), das mit 3,0 µL/min und einer Vorspannung von +5 kV zugeführt wurde. Der iEESI-MS1-Scanbereich betrug m/z 50–2.000. Bei iEESI-MS2 kam eine kollisionsinduzierte Dissoziation mit einer normalisierten Kollisionsenergie von 25–35 %, einem Isolationsfenster von m/z 2,0 und einer dynamischen Ausschlusszeit von 30 s zum Einsatz. Eine tägliche externe Kalibrierung gemäß der instrumentenüblichen Vorgehensweise gewährleistete einen Massenfehler von <5 ppm. Nach medianer Normalisierung, log₂-Transformation und Autoskalierung wurde eine OPLS-DA mit 200 Permutationstests durchgeführt. Die Metaboliten wurden mithilfe der HMDB annotiert, und die Pfadverdünnungsanalyse erfolgte anhand von KEGG-Anmerkungen.
Verarbeitung von Massenspektrometriedaten und Identifizierung von Metaboliten
Die Roh-Massenspektrometriedaten wurden durch Median-Normalisierung normalisiert, gefolgt von einer log₂-Transformation und Autoskalierung (Zentrierung um den Mittelwert und Division durch die Standardabweichung jeder Variablen). Die orthogonale partielle kleinste-Quadrate-Diskriminanzanalyse (OPLS-DA) wurde verwendet, um diskriminierende m/z-Merkmale zu identifizieren. Merkmale wurden als differentiell betrachtet, wenn VIP > 1,0, |log₂(FC)| > 0,58 und p < 0,05 betrug. Die differentiellen Merkmale wurden einer MS2-Analyse unterzogen und mittels der Human Metabolome Database annotiert. Anschließend wurde eine auf KEGG basierende Pfad-Enrichment-Analyse verwendet, um veränderte Stoffwechselwege zu identifizieren und wegbezogene Enzyme zu priorisieren. Die annotierten differentiellen Metaboliten sind in Zusatzdatei 2 enthalten.
Analysen der TCGA-basierten SHMT2-Expression, Überlebens und klinisch-pathologischen Assoziationen
RNA-seq-Expressionsprofile und entsprechende klinische Annotationen für TCGA-LUAD wurden am 18. Dezember 2025 aus dem Genomic Data Commons heruntergeladen; Tumor (n = 542) und normal (n = 59). Fragmente pro Kilobase Transkript pro Million kartierte Reads (FPKM)-Werte wurden als log2(FPKM + 1) transformiert. Für prognostische Analysen wurden nur Proben mit vollständigen Überlebensinformationen einbezogen. Tumor- und normale SHMT2-Expressionswerte wurden mittels Wilcoxon-Rangsummentest verglichen. Die Patienten wurden anhand des medianen Expressionswerts in die Gruppen SHMT2-hoch und SHMT2-niedrig unterteilt. Die Gesamtüberlebenszeit wurde mittels Kaplan-Meier-Analyse und Log-Rank-Test ausgewertet. Assoziationen mit dem Tumor-Node-Metastase (TNM)-Stadium, dem Alter und dem Geschlecht wurden je nach Eignung mittels Chi-Quadrat- oder Wilcoxon-Tests untersucht. Zeitabhängige Receiver-Operating-Characteristic-(ROC)-Kurven wurden für das 1-, 3- und 5-Jahres-Überleben erstellt. Alle statistischen Tests waren zweiseitig.
Differenzielle Expressionsanalyse und funktionelle Anreicherung (GO/KEGG)
Differenziell exprimierte Gene (DEGs) zwischen den SHMT2-hohen und SHMT2-niedrigen Gruppen wurden mittels einer False-Discovery-Rate < 0,05 und |log₂(FC)| ≥ 1 identifiziert. Die Ergebnisse wurden als Volcano-Plot und als Heatmap der 40 stärksten DEGs mit zeilenweiser z-Score-Normalisierung dargestellt. Bei den Gene-Ontology-(GO)- und KEGG-Anreicherungsanalysen wurde ein adjustierter p < 0,05 verwendet. Die Annotationen menschlicher Gene stammten aus Ensembl vom Dezember 2025.
Immuninfiltration, Vorhersage der Immuntherapieantwort und Ableitung der Arzneimittelsensitivität
Tumor-Mikroumgebung-Scores, einschließlich ImmuneScore, StromalScore und ESTIMATEScore, wurden mithilfe von ESTIMATE berechnet. Die Zusammensetzung der Immunzellen wurde mittels CIBERSORT mit dem Referenzsatz LM22 (heruntergeladen am 5. Dezember 2025), 1.000 Permutationen und unter Beibehaltung von Proben mit einer Entmischungs-p < 0,05 abgeleitet. Gruppenweise Vergleiche wurden zwischen SHMT2-reichen und SHMT2-armen Tumoren durchgeführt. Spearman-Korrelationen zwischen SHMT2 und immunologischen Checkpoint-assoziierten Genen wurden als Korrelationsmatrix dargestellt. Immunophenoscores wurden von TCIA bezogen, und TIDE wurde zur ergänzenden Vorhersage der Ansprechbarkeit auf eine Immuntherapie verwendet. Die Arzneimittelsensitivität wurde mittels oncoPredict abgeschätzt, das auf dem GDSC-Release vom Dezember 2023 trainiert wurde, um IC50-Werte zu ermitteln. Für Gruppenunterschiede kamen Wilcoxon-Tests zum Einsatz, und für kontinuierliche Zusammenhänge wurden Spearman-Korrelationen verwendet.
Einzelzell-RNA-seq-Analyse von LUAD
Einzelzelltranskriptomische Daten19 wurden am 20. Dezember 2025 aus dem Gene Expression Omnibus heruntergeladen und analysiert. Zellen mit 200–6.000 detektierten Merkmalen und < 10 % mitochondrialem Transkript wurden beibehalten. Die Zählungen wurden mittels LogNormalize mit einem Skalierungsfaktor von 10.000 normalisiert; die 2.000 variabelsten Gene wurden ausgewählt, und alle Gene wurden skaliert. Die Hauptkomponentenanalyse verwendete die ersten 20 Hauptkomponenten, UMAP nutzte die Dimensionen 1–20, und die Clusterung erfolgte mit einer Auflösung von 0,5. Die Zelltypen wurden mithilfe des Referenzdatensatzes Human Primary Cell Atlas annotiert. Die Anteile der Zelltypen wurden zwischen Tumor- und Normalproben verglichen, und Spearman-Korrelationsmatrizen wurden erstellt. Die Verteilung von SHMT2 wurde über die annotierten Zelltypen hinweg untersucht. Die Zellzykluszustände wurden mittels cyclone abgeleitet und mit Ridge-Plots visualisiert.
Zellkultur
BEAS-2B-, H1299- und A549-Zellen wurden bei 37 °C in einem befeuchteten Inkubator mit 5 % CO2 und einer relativen Luftfeuchtigkeit von etwa 95 % unter Verwendung der in der Tabelle der Materialien angegebenen Medien und Zusätze kultiviert.
Quantitative Echtzeit-Reverse-Transkriptions-PCR (RT-qPCR)
Gesamt-RNA wurde unter Verwendung eines auf Phenol basierenden RNA-Extraktionsreagenz extrahiert. Die extrahierte RNA wurde unter Verwendung eines Mastermixes für die Reverse Transkription in cDNA umgeschrieben. Die Reverse Transkription wurde in einem 20 µL Ansatz durchgeführt, wobei eine genomische DNA-Entfernung bei 42 °C für 2 min, die Reverse Transkription bei 50 °C für 15 min und die Beendigung bei 85 °C für 5 s erfolgte. Die quantitative PCR wurde mit einem auf DNA-Bindungs-Farbstoff basierenden qPCR-Mastermix in einem Echtzeit-PCR-System durchgeführt, wobei eine Initialdenaturierung bei 95 °C für 3 min, gefolgt von 40 Zyklen mit 95 °C für 10 s und 60 °C für 30 s sowie eine Schmelzkurvenanalyse mit 95 °C für 15 s, 60 °C für 60 s und erneut 95 °C für 15 s durchgeführt wurde. Alle Reaktionen wurden im Dreifachansatz durchgeführt. Die relative Expression wurde nach der 2−ΔΔCt-Methode unter Verwendung von GAPDH als Housekeeping-Gen berechnet. Stammlösungen der Primer wurden in doppelt destilliertem Wasser auf 10 µM angesetzt und bei -20 °C im Dunkeln gelagert. Die Primersequenzen sind in Zusatz-Tabelle 3 angegeben.
Immunoblotting
Gesamtzelluläres Protein wurde mit Radioimmunpräzipitationspuffer auf Eis für 30 min extrahiert und bei 4 °C für 15 min bei 12.000 × g zentrifugiert. Die Proteine wurden mittels 12%iger Natriumdodecylsulfat-Polyacrylamid-Gelelektrophorese getrennt und auf 0,22 µm Polyvinylidenfluorid-Membranen übertragen. Die Membranen wurden mit 5% (w/v) fettfreier Milch in Tris-gepuffertem Salzlösung mit 0,1% Tween 20 (TBST) für 2 h bei Raumtemperatur blockiert und anschließend über Nacht bei 4 °C mit den primären Antikörpern gegen SHMT2 (1:1.000) und gegen GAPDH (1:2.000) inkubiert, die in 5% fettfreier Milch-TBST verdünnt waren. Nach dreimaligem Waschen jeweils 8 min mit TBST (20 mM Tris-HCl, 150 mM NaCl, 0,1% Tween-20) wurden die Membranen für 1 h bei Raumtemperatur mit mit Peroxidase konjugierten sekundären Antikörpern (1:5.000) inkubiert. Die Signale wurden mit einem durch Mischen gleicher Volumina der Lösungen A und B hergestellten verbesserten Chemilumineszenz-Substrat entwickelt, für 2 min bei Raumtemperatur im Dunkeln inkubiert und mit einem Chemilumineszenz-Bildgebungssystem erfasst.
Validierung der HPA-Immunhistochemie
Die Expressionsmuster des SHMT2-Proteins in normalem Lungengewebe und bei LUAD wurden mithilfe der HPA analysiert. „SHMT2“ wurde in der HPA abgefragt, und repräsentative Immunhistochemie-(IHC-)Bilder sowie zugehörige pathologische Annotationen wurden zum Vergleich heruntergeladen. Um Verzerrungen aufgrund der Antikörperspezifität zu minimieren, wurden Antikörper mit höheren Validierungsstufen in der HPA bevorzugt. Wo verfügbar, wurden IHC-Ergebnisse mehrerer Antikörper gegenseitig überprüft (HPA020543 und HPA020549).
Verringerung der SHMT2-Expression in LUAD-Zellen
Um die tumorspezifische Rolle von SHMT2 bei LUAD zu untersuchen, wurde eine SHMT2-Knockdown-Behandlung in H1299- und A549-Zellen mittels shRNA durchgeführt. Die shRNA wurde von einem über einen H1-Promotor gesteuerten, mit Puromycin selektierbaren lentiviralen Vektor getragen; Verpackungsplasmide wurden während der Virusproduktion cotransfiziert. sh-SHMT2 zielte auf den kodierenden Bereich des menschlichen SHMT2 ab, während sh-NC eine zufällig zusammengestellte negative Kontrolle darstellte. Die hergestellten Viruspräparationen wiesen eine Titerkonzentration von ≥ 1 × 108 transduzierenden Einheiten/mL auf und wurden bei -80 °C gelagert. Die Zellen wurden mit den in der Materialtabelle beschriebenen sh-SHMT2- oder sh-NC-lentiviralen Partikeln infiziert. Nach 48 Stunden wurde die Effizienz des Knockdowns mittels RT-qPCR und Western Blot analysiert. Nur Zellen mit nachgewiesener SHMT2-Silenzierung wurden für Proliferations-, Migrations- und Koloniebildungs-Assays verwendet.
Proliferationsassay für Zellen
Die Zellproliferation wurde mittels CCK-8-Assay bewertet. Transfizierte H1299- und A549-Zellen wurden in 96-Well-Platten in Sextuplikaten mit je 5 x 103 Zellen pro Well ausgesät. An den Tagen 1, 2, 3, 4 und 5 wurden jeweils 10 µL CCK-8-Reagenz (10 % des Kulturmediumvolumens) zu jedem Well hinzugefügt und 2 h bei 37 °C inkubiert. Die Extinktion wurde bei 450 nm gemessen und als Parameter für die Zellviabilität herangezogen.
Wundheilungs-Assay
Die Zellmigration wurde mithilfe eines Wundheilungsassays bewertet. Transfizierte H1299- und A549-Zellen wurden in 6-Well-Platten ausgesät und bis zur Konfluenz von 90–95 % kultiviert. Eine lineare Wunde wurde mit einer sterilen 200-µL-Pipettenspitze erzeugt, und abgelöste Zellen wurden durch vorsichtiges Waschen mit phosphatgepufferter Salzlösung entfernt. Die Zellen wurden anschließend in serumfreiem oder niedrigserumhaltigem Medium kultiviert, und der Wundbereich wurde nach 0 h und 24 h unter einem inversen Mikroskop dokumentiert. Die Wundschlussrate wurde mit der im Materialverzeichnis aufgeführten Bildanalyse-Software quantifiziert und als prozentuale Abnahme im Vergleich zur anfänglichen Wundfläche berechnet.
Test zur Koloniebildung
Für die Koloniebildungsanalyse wurden transfizierte H1299- und A549-Zellen in 6-Loch-Platten mit jeweils 800 Zellen pro Loch ausgesät und etwa 10–14 Tage lang kultiviert, bis sichtbare Kolonien entstanden, wobei die Kultur mikroskopisch überwacht wurde. Die Kolonien wurden 30 Minuten bei Raumtemperatur mit 4 % Paraformaldehyd fixiert und anschließend 15 Minuten bei Raumtemperatur mit 0,1 % Kristallviolett angefärbt. Kolonien mit mehr als 50 Zellen wurden unter dem Lichtmikroskop gezählt.
Statistische Analyse für in vitro-Experimente
Alle in vitro-Experimente wurden unabhängig voneinander mindestens dreimal durchgeführt. Die Daten werden als Mittelwert ± Standardabweichung angegeben. Für den Vergleich von zwei Gruppen wurde der t-Test nach Student verwendet, und ein p < 0,05 wurde als statistisch signifikant angesehen.