Die Studie wurde von der Ethikkommission des Provinzkrankenhauses Shaanxi (Genehmigungsnummer 2024-R055) genehmigt und gemäß den Grundsätzen der Deklaration von Helsinki durchgeführt. Von jedem Teilnehmer wurde vor der Blutentnahme eine schriftliche, informierte Einwilligung eingeholt. Alle menschlichen Blut- und Serumproben wurden gemäß den institutionellen Richtlinien der Abteilung für Labormedizin gesammelt, bearbeitet und gelagert.
Alle venösen Blut- und Serumproben wurden als potenziell infektiöses biologisches Material behandelt. Die Probenaufbereitung erfolgte in einem Labor der Biosicherheitsstufe 2 durch geschultes Personal, das die entsprechende persönliche Schutzausrüstung trug. Einweg-Plastikwaren, Pipettenspitzen und verwendete Mikroplatten, die mit biologischem Material in Kontakt gekommen waren, wurden entweder autoklaviert oder als biohazardöser Abfall gesammelt und gemäß lokalen Vorschriften über den lizenzierten medizinischen Abfallstrom der Einrichtung entsorgt.
Überblick zum Studiendesign:
Der Arbeitsablauf dieser Studie ist in einem schematischen Studiendesign dargestellt (Supplementary Figure S1) und in zwei komplementäre Arme unterteilt. Im in silico-Arm wurden öffentliche transkriptomische, Einzelzell-, Protein-, posttranslationale Modifikations-, Immun- und Arzneimittelsensitivitäts-Datensätze analysiert, um den biologischen Kontext von hWAPL zu charakterisieren und mechanistische Hypothesen zu generieren. Im experimentellen Arm wurden hWAPL und SCC-Ag im Serum in einem unabhängigen Einzelzentrum-Kohortenkollektiv gemessen, und ihre diagnostische Leistung – einzeln und in Kombination – wurde mittels ROC-Analyse (Receiver Operating Characteristic) und logistischer Regression bewertet. Die beiden Arme wurden ausschließlich auf der Ebene der Interpretation zusammengeführt, wobei der Serum-Arm die einzige experimentelle Validierung lieferte.
Genexpressionsprofilierung bei Pan-Krebs und Plattenepithelkarzinom des Zervix
RNA-Sequenzierungsdaten aus dem Cancer Genome Atlas (TCGA) wurden über den UCSC Xena Browser bezogen. Alle in dieser Studie verwendeten öffentlichen Datenbanken und Web-Tools wurden zwischen September und Oktober 2025 abgerufen, und die jeweiligen Versionen sowie Zugriffsdaten sind in der Tabelle der Materialien aufgeführt. Die Expressionsniveaus von WAPL, quantifiziert als Fragmente pro Kilobase Transkriptlänge pro Million kartierte Reads (FPKM), wurden einer log₂-Transformation unterzogen, wobei ein Pseudocount von 1 addiert wurde. Diese transformierten Werte wurden anschließend verwendet, um Vergleichsanalysen zwischen malignen Geweben und jeweils passenden angrenzenden normalen Proben über mehrere Krebsarten hinweg durchzuführen. Für die gezielte Untersuchung von Zervixmalignitäten wurden WAPL-Expressionsdaten spezifisch für das Plattenepithelkarzinom und das Endozervixkarzinom (CESC) extrahiert. Das Expressionsprofil in diesen Tumorproben wurde anschließend visualisiert und mit den Expressionsniveaus in normalen zervikalen epithelialen Kontrollgeweben verglichen.
Analyse der WAPL-Transkript-Isoformen
Die Expressionsprofile der WAPL-Transkript-Isoformen wurden mithilfe der GEPIA2-Datenbank untersucht. Spezifische Isoformen wurden anhand ihrer entsprechenden Ensembl-Transkript-Identifikatoren (ENST) abgegrenzt. Die Expressionsniveaus wurden als Logarithmus zur Basis 2 der Transkripte pro Million plus eins [log₂(TPM+1)] berechnet und grafisch als Violin-Plots dargestellt.
Analyse der Einzelzell-RNA-Sequenzierung
Der scRNA-seq-Datensatz GSE168652 wurde aus dem Gene Expression Omnibus (GEO)-Repository bezogen. Die Datenvorverarbeitung, Qualitätsbewertung, Normalisierung und Clusterung wurden mit der Plattform scCancerExplorer durchgeführt. Zellen mit weniger als 200 exprimierten Genen sowie Gene, die in weniger als drei Zellen nachgewiesen wurden, wurden ausgeschlossen. Die Normalisierung erfolgte mithilfe des Log-Normalisierungsansatzes mit einem Skalierungsfaktor von 10.000.
Die 2.000 variabelsten Gene wurden mithilfe einer varianzstabilisierenden Transformation identifiziert und anschließend für die Hauptkomponentenanalyse (PCA) verwendet. Zur Dimensionsreduktion wurde die uniforme Mannigfaltigkeitsapproximation und -projektion (UMAP) basierend auf den ersten 30 Hauptkomponenten angewandt. Die graphenbasierte Clusterung erfolgte über den Louvain-Algorithmus mit einer Auflösung von 0,8. Zellcluster wurden anhand etablierter Markergene annotiert. Die Expressionsniveaus von WAPL wurden in verschiedenen Zellpopulationen gemessen, und die Zellzyklusphasen wurden mittels der CellCycleScoring-Methode bestimmt.
Proteinausdruck und Analyse posttranslationaler Modifikationen
Der Ausdruck und die subzelluläre Verteilung von WAPL in CSCC wurden mithilfe der Human Protein Atlas-Datenbank bewertet. Posttranslationale Modifikationsstellen (PTM) wurden mit PhosphoSitePlus identifiziert. Darüber hinaus wurden potenzielle Assoziationen zwischen den WAPL-Expressionsniveaus und Genen, die an Proteinmodifikationsprozessen beteiligt sind – insbesondere AKT1, AKT2, KAT2B, EP300, USP7 und KAT5 – im Kontext von CSCC mit dem Xiantao Academic-Tool13,14,15,16,17 untersucht.
Ko-Expressionsnetzwerk- und funktionelle Anreicherungsanalyse
Die mittels des STAR-Aligners verarbeiteten und auf Transkripte pro Million (TPM) normalisierten RNA-Sequenzierungsdaten für TCGA-CSCC wurden vom Genomic Data Commons (GDC)-Portal bezogen. Um Gene zu identifizieren, die ko-exprimiert sind mit WAPL, wurde ein Screening-Schwellenwert basierend auf einem absoluten Pearson-Korrelationskoeffizienten von mehr als 0,6 und einem statistischen Signifikanzniveau von P < 0,05 angewandt.
Um die mit WAPL assoziierten biologischen Funktionen zu untersuchen, wurde zunächst eine Analyse der differentiellen Expression durchgeführt, bei der die Patienten anhand des medianen Expressionswerts als Grenzwert in Gruppen mit hoher und niedriger WAPL-Expression eingeteilt wurden. Differenziell exprimierte Gene wurden mit der Xiantao-Academic-Plattform identifiziert und anschließend einer Anreicherungsanalyse der Genontologie unter Einbeziehung der Kategorien biologischer Prozess, zelluläre Komponente und molekulare Funktion sowie einer Pfadanalyse der Kyoto-Enzyklopädie für Gene und Genome (KEGG) unterzogen. Unabhängig davon wurden Protein-Interaktionsbeziehungen, an denen WAPL beteiligt ist, mittels STRING untersucht, und das resultierende Interaktionsnetzwerk wurde in Cytoscape rekonstruiert und analysiert.
Die Gen-Set-Enrichment-Analyse (GSEA) wurde unter Verwendung der Datenbanken Reactome, WikiPathways und KEGG durchgeführt, wobei die Screening-Kriterien bei einem absoluten normalisierten Anreicherungsscore (|NES|) > 1 und einer falschen Entdeckungsrate (FDR) < 0,25 festgelegt wurden.
Analysen der Immuninfiltration
Die rechnergestützten Werkzeuge TIMER 2.0 und xCell wurden verwendet, um transkriptomische Daten aus dem TCGA-CESC-Kollektiv zu analysieren, wodurch die quantitative Beurteilung des Infiltrationsstatus für 24 unterschiedliche Immunzellpopulationen ermöglicht wurde.
Vorhersage der Arzneimittelsensitivität
Vorhersagen zur Arzneimittelsensitivität wurden anhand der GDSC2-Datenbank mithilfe des oncoPredict-R-Pakets abgeleitet. Die halbmaximale Hemmkonzentration (IC50) wurde für ein Panel von 86 chemotherapeutischen Verbindungen berechnet.
Überlebens- und Diagnoseanalyse
Kurven für das Gesamtüberleben (OS) und das progressionsfreie Intervall (PFI) wurden mithilfe der Xiantao-Analyseplattform erstellt, wobei die Patientenkohorten anhand des medianen WAPL-Expressionsniveaus in zwei Gruppen eingeteilt wurden. Die prognostische Bedeutung der Variablen wurde anhand univariater und multivariater Cox-Regressionmodelle mit proportionalen Hazards bewertet. Um die diagnostische Leistungsfähigkeit von WAPL und SCC (SERPINB3) bei der Unterscheidung von CSCC und angrenzendem Normalgewebe zu prüfen, wurde eine ROC-Kurvenanalyse (Receiver Operating Characteristic) auf Transkriptomdaten angewandt, die als Transkripte pro Million (TPM) aus dem Cancer Genome Atlas stammten.
Sammlung von Serumproben und ELISA
Diese Studie erhielt die Genehmigung der Institutional Review Board des Provinzialkrankenhauses Shaanxi (Genehmigungsnummer 2024-R055) und wurde in Übereinstimmung mit der Deklaration von Helsinki durchgeführt. Insgesamt wurden 89 Serumproben von Patientinnen mit histopathologisch bestätigtem CSCC sowie 89 altersangepasste gesunde weibliche Kontrollen gewonnen. Das mittlere Alter betrug in der CSCC-Gruppe 54 Jahre (Spannweite 31–77 Jahre) und in der Kontrollgruppe 41 Jahre (Spannweite 30–52 Jahre). Alle Gebärmutterhalskrebsfälle waren zum Zeitpunkt der Blutentnahme neu diagnostiziert und noch nicht therapiert; die Proben wurden vor Beginn jeglicher Operation, Chemotherapie oder Strahlentherapie entnommen. Wir weisen darauf hin, dass, obwohl die beiden Gruppen aus einem vergleichbaren klinischen Umfeld stammten, der Altersunterschied zwischen den Gruppen 13 Jahre betrug. Dieser Unterschied wird direkt im Ergebnisteil und in einer speziellen altersadjustierten Sensitivitätsanalyse untersucht und im Diskussionsteil als Limitation anerkannt.
Venöses Blut wurde in röhrchen ohne Antikoagulanzien gesammelt, bei Raumtemperatur 10–20 min gerinnen gelassen und anschließend 20 min bei 2.000–3.000 × g zentrifugiert. Das Serum wurde innerhalb von 1 h nach der Entnahme abgetrennt, unverzüglich aliquotiert und bei −20 °C gelagert; das Zeitintervall zwischen der Phlebotomie und dem Einfrieren überschritt bei keiner Probe 2 h. Die Aliquots wurden vor der Analyse nur einem einzigen Gefrier-Tau-Zyklus unterzogen, und alle Messungen wurden innerhalb von drei Monaten nach der Lagerung abgeschlossen. Das erhaltene Serum wurde in Aliquots aufgeteilt und bei −20 °C aufbewahrt. Hämolytische oder lipämische Proben wurden von der Analyse ausgeschlossen.
Die Quantifizierung der Serumspiegel des hWAPL-Proteins erfolgte mithilfe eines kommerziell erhältlichen enzymgekoppelten Immunosorbens-Assays (ELISA) nach dem Prinzip des Doppelantikörpersandwichs. Kurz gesagt, wurden sowohl Kalibrierstandards als auch Testproben verdoppelt in Mikrotiterplattenmulden pipettiert, die zuvor mit einem für hWAPL spezifischen Erfassungsantikörper beschichtet worden waren, gefolgt von einer 30-minütigen Inkubation bei 37 °C. Nach einer Reihe von Waschschritten zur Entfernung nicht gebundenen Materials wurde ein gegen hWAPL gerichteter Detektionsantikörper, markiert mit Meerrettichperoxidase (HRP), hinzugegeben und unter identischen Bedingungen inkubiert.
Nach weiterem Waschen wurden jeweils 50 µL der chromogenen Substratlösungen A und B in die Vertiefungen gegeben und 10 min im Dunkeln bei 37 °C zur Entwicklung inkubiert. Anschließend wurden 50 µL Stopplösung zugegeben, um die enzymatische Reaktion zu beenden. Die optische Dichte wurde innerhalb von 15 min nach Reaktionsstopp bei einer primären Wellenlänge von 450 nm mit einer Referenzablesung bei 630 nm gemessen. Unter Verwendung seriell verdünnter Kalibratoren, die einen Konzentrationsbereich von 10–450 pg/mL abdeckten, wurde eine Standardkurve erstellt. Die hWAPL-Konzentration jeder unbekannten Probe wurde durch Interpolation aus dieser Standardkurve ermittelt und entsprechend eventueller vorgängiger Verdünnungen angepasst. Der Test wurde streng gemäß den Herstelleranweisungen durchgeführt. Gemäß den Spezifikationen des Herstellers liegt der funktionale Nachweisbereich des Kits bei 10–450 pg/mL. Die Proben wurden abschließend fünffach verdünnt gemessen; Proben, deren optische Dichte über der des höchsten Kalibrators lag, wurden zusätzlich verdünnt und erneut analysiert, wobei der gesamte Verdünnungsfaktor bei der Rückrechnung berücksichtigt wurde. Die analytische Leistung wurde im Serum-Matrix vor der Analyse der Studienproben überprüft: Ein seriell verdünnter Serum-Pool mit hoher Konzentration zeigte eine akzeptable Linearität über den gesamten Kalibrierbereich, und die Wiederfindungsrate des rekombinanten Kalibrators in drei Serum-Pools lag innerhalb eines akzeptablen Bereichs, was mit dem Double-Antibody-Sandwich-Design des Kits übereinstimmt. Die untere Quantifizierungsgrenze entsprach dem niedrigsten Kalibrator (10 pg/mL). Wir weisen darauf hin, dass dieses Kit laut Herstellerangabe ausschließlich für Forschungszwecke bestimmt ist und kein klinisch validiertes In-vitro-Diagnostikum darstellt; die absoluten Konzentrationen und die abgeleiteten Schwellenwerte sollten daher, wie unten diskutiert, als gerätespezifisch und vorläufig betrachtet werden.
Parallel dazu wurden die Serum-SCC-Ag-Spiegel mit der Abbott Alinity i-Plattform mittels eines chemilumineszenten mikropartikulären Immunoassays (CMIA) quantifiziert, wobei die obere Referenzgrenze bei 1,5 ng/mL lag. Alle analytischen Verfahren wurden im Duplikat von Laborpersonal durchgeführt, das bezüglich des Krankheitsstatus und aller klinischen Variablen – einschließlich klinischem Stadium, histologischem Grad und HPV-Status – verblindet war. Sowohl die interassay- als auch die intraassay-Variationskoeffizienten der eingesetzten Methoden lagen unter 10 %.
Statistische Analyse
Da die Konzentrationen von hWAPL im Serum gemäß dem Shapiro-Wilk-Test nicht der Normalverteilungsannahme genügten (alle P < 0,05), wurden die nachfolgenden Analysen mit nichtparametrischen Methoden durchgeführt. Stetige Variablen werden daher als Mediane mit Interquartilsbereichen (IQRs) angegeben. Unterschiede zwischen zwei unabhängigen Gruppen wurden mittels des Mann-Whitney-U-Tests bewertet, während Vergleiche zwischen drei oder mehr Gruppen mit dem Kruskal-Wallis-H-Test durchgeführt wurden. Zusammenhänge zwischen Variablen wurden mithilfe der Rangkorrelationsanalyse nach Spearman für stetige oder ordinale Daten sowie mit dem punktbiserialen Korrelationskoeffizienten für binäre Variablen, wie den HPV-Infektionsstatus, untersucht. Überlebensendpunkte wurden mit dem Log-Rang-Test verglichen, und Hazard-Ratios mit 95 % Konfidenzintervallen wurden anhand von Cox-Regressionen mit proportionalen Hazards geschätzt.
Die diagnostische Leistung wurde mittels ROC-Kurvenanalyse bewertet, wobei der optimale Cut-off-Wert für jeden Biomarker durch Maximierung des Youden-Index bestimmt wurde. Unterschiede in der Fläche unter der Kurve (AUC) wurden mit der DeLong-Methode geprüft. Um die kombinierte diagnostische Aussagekraft zu untersuchen, wurde ein binäres logistisches Regressionsmodell mit dem Krankheitsstatus (CSCC = 1; Kontrolle = 0) als abhängige Variable und dem serumgebundenen hWAPL sowie dem SCC-Antigen als unabhängige Prädiktoren erstellt. Die resultierenden vorhergesagten Wahrscheinlichkeiten wurden als kompositiver Score für die anschließende ROC-Analyse gespeichert.
Die unterschiedliche Empfindlichkeit gegenüber Medikamenten zwischen den Gruppen mit hohem und niedrigem WAPL-Spiegel wurde mit dem limma-Paket bewertet, wobei Schwellenwerte von einem adjustierten P < 0,05 und einem |log₂-Fold-Change| > 0,3 angewandt wurden. Statistische Signifikanz wurde als zweiseitiger P-Wert < 0,05 definiert. Alle Analysen wurden mit SPSS, R oder den angegebenen Online-Tools durchgeführt. Zu den verwendeten R-Paketen gehörten survival, survminer, pROC, ggplot2, oncoPredict und limma. Bei allen hochdimensionalen Analysen, bei denen viele Hypothesen gleichzeitig geprüft wurden – nämlich den pan-kutanen Vergleichen, dem Ko-Expressions-Screening, den Analysen zur differentiellen Expression und funktionalen Anreicherung sowie der Analyse der Medikamentenempfindlichkeit – wurden die P-Werte mittels der False-Discovery-Rate nach Benjamini–Hochberg für multiple Vergleiche korrigiert, und es werden, wo zutreffend, die adjustierten P-Werte angegeben; die paarweisen Korrelationen zwischen hWAPL und einzelnen Genen oder Anteilen von Immunzellen werden mit nominalen P-Werten berichtet, und diese Korrelationsergebnisse gelten als hypothesebildend, nicht als bestätigend. Da einige klinisch-pathologische Felder im öffentlichen Kollektiv unvollständig besetzt waren, wurden Fälle mit fehlenden Daten paarweise aus der jeweiligen Analyse ausgeschlossen, anstatt sie zu imputieren, und die Anzahl der auswertbaren Fälle wird für jede Variable in Tabelle 1 angegeben. Um den Altersunterschied zwischen den Serumkollektiven zu berücksichtigen, wurde eine Sensitivitätsanalyse mittels altersadjustierter binärer logistischer Regression durchgeführt, und die Spearman-Korrelation zwischen serumgebundenem hWAPL und Alter wurde innerhalb jeder Gruppe berechnet.