מחקר זה אושר על ידי ועדת ההגנה על נסייני אדם (IRB) של בית החולים המחוזי שאןשי (מספר אישור 2024-R055) ונערך בהתאם לעקרונות הצהרת הלסינקי. הסכמה מדעת בכתב התקבלה מכל משתתף לפני דגימת דם. כל דגימות הדם והסרום האנושיות נאספו, טופלו ואוחסנו בהתאם להנחיות המוסדיות של המחלקה לרפואה מעבדתית.
כל דגימות דם ורידי וסרום טופלו כחומר ביולוגי שעלול להיות זיהומי. עיבוד הדגימות בוצע במעבדה ברמת בטיחות ביולוגית 2 (biosafety level 2) על ידי כוח אדם מיומן שהשתמש בציוד מיגון אישי מתאים. כלי פלסטיק חד-פעמיים, טיפים של פיפטות ומיקרו-פלטות משומשות שבאו במגע עם חומר ביולוגי עברו עיקור באוטוקלאב או נאספו כפסולת ביו-סיכונית והושלכו דרך ערוץ פסולת הרפואית המורשה של המוסד, בהתאם לתקנות המקומיות.
סקירת תכנון המחקר:
זרימת העבודה של מחקר זה מסוכמת בסכימה של תכנון המחקר (Supplementary Figure S1) ומחולקת לשני זרועות משלימות. בזרוע ה-in silico, נערכה אנליזה של מאגרי נתונים ציבוריים של תמלול (transcriptomic), תאים בודדים, חלבונים, מודיפיקציות לאחר תרגום, מערכת החיסון ורגישות לתרופות, במטרה לאפיין את ההקשר הביולוגי של hWAPL ולהעלות היפותזות מכניסטיות. בזרוע הניסויית, נמדדו serum hWAPL ו-SCC-Ag בקוהורט עצמאי ממרכז רפואי יחיד, וביצועיהם האבחנתיים, כביגוד נפרד ובשילוב, הוערכו באמצעות אנליזת עקומי ROC (receiver operating characteristic) ורגרסיה לוגיסטית. שתי הזרועות שולבו רק בשלב הפרשנות, כאשר זרוע הסרום סיפקה את התיקוף הניסויי היחיד.
פרופיל ביטוי גנים של סרטן צוואר הרחם מסוג קרצינומה של תאי קשקש וסרטנים אחרים
נתוני ריצוף RNA מתוך The Cancer Genome Atlas (TCGA) הושגו באמצעות ה-UCSC Xena Browser. לכל מאגרי המידע הציבוריים וכלי הרשת ששימשו במחקר זה ניגשנו בין ספטמבר לאוקטובר 2025, והגרסאות הספציפיות ותאריכי הגישה מפורטות ב-טבלת החומרים. רמות הביטוי של WAPL, שכומי הן כfragments per kilobase of transcript per million mapped reads (FPKM), עברו טרנספורמציית log₂ עם תוספת של pseudocount של 1. ערכים שונו אלו שימשו לביצוע ניתוחים השוואתיים בין רקמות ממאירות לדגימות נורמליות סמוכות תואמות במספר סוגי סרטנים. לצורך החקירה הממוקדת של ממאירות צוואר הרחם, בודדו נתוני ביטוי של WAPL ספציפיים לקרצינומה של תאי קשקש ואדנוקרצינומה אנדו-צווארית של צוואר הרחם (CESC). פרופיל הביטוי בדגימות גידולים אלו הוצג ויזואלית לאחר מכן והושווה לרמות הביטוי ברקמות בקרה של אפיתל נורמלי בצוואר הרחם.
ניתוח של איזופורמים של תמלילי WAPL
פרופילי הביטוי של איזופורמים של תמלילי WAPL נבחנו באמצעות מאגר הנתונים GEPIA2. איזופורמים ספציפיים הוגדרו על ידי מזהי התמלילים המתאימים שלהם ב-Ensembl (ENST). רמות הביטוי חושבו כערך הלוגריתם בבסיס 2 של תמלילים למיליון פלוס אחת [log₂(TPM+1)] והוצגו גרפית כתרשימי כינור (violin plots).
ניתוח רצף RNA של תא בודד
מערכת הנתונים של scRNA-seq שמספרה GSE168652 הופקה ממאגר Gene Expression Omnibus (GEO). עיבוד מקדמי של נתונים, הערכת איכות, נרמול ואשכול (clustering) בוצעו באמצעות פלטפורמת scCancerExplorer. תאים שהביעו פחות מ-200 גנים וגנים שזוהו בפחות משלושה תאים הוסרו בסינון. הנרמול בוצע בשיטת Log-Normalize עם מקדם קנה מידה (scaling factor) של 10,000.
2,000 הגנים בעלי השונות הגבוהה ביותר זוהו באמצעות טרנספורמציה למייצבי שונות (variance-stabilizing transformation) ושימשו לאחר מכן לניתוח רכיבים ראשיים (PCA). להפחתת ממדיות הוחל שימוש ב-Uniform Manifold Approximation and Projection (UMAP) על בסיס 30 הרכיבים הראשיים הראשונים. אשכולה (Clustering) מבוססת גרף בוצע באמצעות אלגוריזם Louvain ברזולוציה של 0.8. אשכולות התאים סווגו על פי גנים מסמנים (marker genes) מובהקים. רמות הביטוי של WAPL נמדדו באוכלוסיות תאים שונות, ושלבי מחזור התא נקבעו באמצעות שיטת CellCycleScoring.
ניתוח ביטוי חלבונים ומודיפיקציות לאחר תרגום
הביטוי וההתפלגות התת-תאית של WAPL ב-CSCC הוערכו באמצעות מאגר הנתונים Human Protein Atlas. אתרי מודיפיקציה לאחר תרגום (PTM) זוהו באמצעות PhosphoSitePlus. יתרה מכך, קשרים פוטנציאליים בין רמות הביטוי של WAPL לבין גנים המעורבים בתהליכי מודיפיקציה של חלבונים, ובסגול AKT1, AKT2, KAT2B, EP300, USP7 ו-KAT5, בתוך CSCC, נחקרו באמצעות כלי Xiantao Academic13,14,15,16,17.
רשת ביטוי משותף ואנליזת העשרה תפקודית
נתוני רצף ה-RNA עבור TCGA-CSCC, שעובדו באמצעות מתאם STAR ונורמלו למספר תמלילים למיליון (TPM), הושגו מפורטל Genomic Data Commons (GDC). כדי לזהות גנים המובעים במשותף עם WAPL, הוחל סף סינון המבוסס על מקדם מתאם פירסון מוחלט הגבוה מ-0.6 ורמת מובהקות סטטיסטית של P < 0.05.
על מנת לחקור את התפקודים הביולוגיים הקשורים ל-WAPL, בוצע תחילה ניתוח ביטוי דיפרנציאלי לאחר חלוקת החולים לקבוצות של ביטוי WAPL גבוה ונמוך, תוך שימוש בערך הביטוי החציוני כנקודת חיתוך. גנים עם ביטוי דיפרנציאלי זוהו באמצעות פלטפורמת Xiantao Academic והועמדו לאחר מכן לניתוח העשרה של Gene Ontology, כולל קטגוריות של תהליך ביולוגי, רכיב תאי ותפקוד מולקולרי, יחד עם ניתוח מסלולים של Kyoto Encyclopedia of Genes and Genomes. באופן בלתי תלוי, נבחנו קשרי אינטראקציה של חלבונים המערבים את WAPL באמצעות STRING, ורשת האינטראקציות שהתקבלה שוחזרה ונבחנה ב-Cytoscape.
ניתוח העשרה של קבוצות גנים (GSEA) בוצע באמצעות מאגרי המידע Reactome, WikiPathways ו-KEGG, כאשר קריטריוני הסינון נקבעו על מדד העשרה מנורמל מוחלט (|NES|) > 1 ושיעור גילוי שגוי (FDR) < 0.25.
ניתוח חדירה חיסונית
הכליםB החישובים TIMER 2.0 ו-xCell שימשו לניתוח נתונים טרנסקריפטומיים מקבוצת ה-TCGA-CESC, ובכך אפשרו הערכה כמותית של מצב החדירה עבור 24 אוכלוסיות נבדלות של תאי מערכת החיסון.
חיזוי רגישות לתרופות
חיזויים של רגישות לתרופות הופקו ממאגר ה-GDSC2 באמצעות חבילת ה-R שנקראת oncoPredict. ריכוז העכבה החצאי-מקסימלי (IC50) חושב עבור פאנל של 86 תרכובות כימותרפיוטיות.
ניתוח הישרדות ואבחנה
עקומי הישרדות כוללת (OS) ומרווח ללא התקדמות (PFI) נבנו באמצעות פלטפורמת הניתוח של Xiantao, כאשר קבוצות המטופלים חולקו לשתי קבוצות על בסיס רמת הביטוי החציונית של WAPL. המשמעות הפרוגנוסטית של המשתנים הוערכה באמצעות מודלים של רגרסיית סיכונים פרופורציונלית של Cox, הן בניתוח חד-משתני והן בניתוח רב-משתני. כדי להעריך את היכולת האבחנתית של WAPL ו-SCC (SERPINB3) בהבחנה בין CSCC לבין רקמה נורמלית סמוכה, הוחל ניתוח עקומות מאפיין תגובה (ROC) על נתונים טרנסקריפטומיים, המבוטאים כמספר טרנסקריפטים למיליון (TPM), שהופקו מ-The Cancer Genome Atlas.
איסוף דגימות סרום ו-ELISA
מחקר זה קיבל אישור מוועדה לאתיקה במחקר אנושי (IRB) של בית החולים המחוזי שאנשי (מספר אישור 2024-R055) ובוצע בהתאם להצהרת הלסינקי. סך כולל של 89 דגימות סרום התקבלו ממטופלות עם CSCC שהוכח היסטופתולוגית, לצד 89 נשויות בריאות כביקורת, שהותאמו מבחינת גיל. הגיל החציוני היה 54 שנים, בטווח של 31–77 שנים, בקבוצת ה-CSCC, ו-41 שנים, בטווח של 30–52 שנים, בקבוצת הביקורת. כל מקרי סרטן צוואר הרחם אובחנו לאחרונה והיו נאיביים לטיפול במועד איסוף הדם; הדגימות נלקחו לפני תחילת כל ניתוח, כימותרפיה או רדיותרפיה. נציין כי למרות ששתי הקבוצות נלקחו מהסביבה הקלינית דומה, הגיל החציוני השתנה ב-13 שנים. הבדל זה נבחן ישירות בפרק התוצאות ובניתוח רגישות ייעודי המתואם לגיל, והוא מוכר כמגבלה בפרק הדיון.
דם ורידי נאסף למבחנות ללא נוגדי קרישה, הושאר להיקרש בטמפרטורת החדר למשך 10–20 min, ולאחר מכן הופרד בצנטריפוגה ב-2,000–3,000 × g למשך 20 min. הסרום הופרד בתוך שעה מהאיסוף, חולק למכסות (aliquots) באופן מיידי ואושרש ב-−20 °C; פרק הזמן בין הדקירה הוורידית להקפאה לא עלה על 2 h עבור אף דגימה. המכסות עברו מחזור הקפאה-הפשרה בודד לפני האנליזה, וכל המדידות הושלמו בתוך שלושה חודשים מהאחסון. הסרום שהתקבל חולק למכסות ונשמר ב-−20 °C. דגימות המולטיות או ליפמיות הוצאו מהאנליזה.
כימות רמות חלבון hWAPL בסרום בוצע באמצעות ערכה מסחרית של ספיחה חיסונית מקישור אנזימטי (ELISA), המבוססת על עיקרון "סנדוויץ'" של נוגדויות כפולות. בקצרה, הן סטנדרטי הכיול והן דגימות הבדיקה חולקו בשכפולים לבארי microplate מצופים מראש בנוגדן לכידה ספציפי ל-hWAPL, ולאחריהם דגירה של 30 min בטמפרטורה של 37 °C. לאחר סדרת שלבי שטיפה להסרת חומרים שאינם קשורים, הוסף נוגדן זיהוי המסומן ב-horseradish peroxidase (HRP) המכוון ל-hWAPL, והוא דגור תחת תנאים זהים.
לאחר שטיפה נוספת, תמיסו תת-מצע כרומוגניות A ו-B בנפח של 50 µL כל אחת בכל באר, והן הושארו להתפתח בחשכה בטמפרטורה של 37 °C למשך 10 min. תמיסת עצירה (50 µL) הוספה כדי להפסיק את התגובה האנזימטית. הצפיפות האופטית נמדדה באורך גל ראשי של 450 nm, עם קריאת ייחוס ב-630 nm, בתוך 15 min מהפסקת התגובה. עקומת סטנדרט נבנתה באמצעות כיליורים בדילול סדרתי בטווח ריכוזים של 10–450 pg/mL. ריכוז ה-hWAPL עבור כל דגימה לא ידוע הופק באמצעות אינטרפולציה מעקומת סטנדרט זו והותאם לכל דילול פרה-אנליטי. הבדיקה בוצעה בהחלט לפי הוראות היצרן. על פי מפרטי היצרן, טווח הגילוי העובד של הערכה הוא 10–450 pg/mL. הדגימות נבדקו בדילול סופי של 5-fold, ודגימות שהצפיפות האופטית שלהן עלתה מזו של הכיליור הגבוה ביותר דוללו באופן נוסף ונבדקו שוב, תוך החלת מקדם הדילול הכולל במהלך חישוב הערכים חזרה. הביצועים האנליטיים אומתו במטריקס סרום לפני הרצת דגימות המחקר: מאגר סרום בריכוז גבוה בדילול סדרתי הראה לינאריות מקובלת לאורך טווח הכיול, ואחוזי התאוששות (spike recovery) של הכיליור הרקומביננטי בשלושה מאגרי סרום היו בטווח מקובל, בהתאם למבנה ה-double-antibody sandwich של הערכה. הגבול התחתון של הכמות (LLQ) תאם לכיליור הנמוך ביותר (10 pg/mL). אנו מדגישים כי הערכה זו מסומנת על ידי היצרן לשימוש למחקר בלבד ואינה מכשיר אבחון in vitro המקבל תיקוף קליני; לפיכך, יש לפרש את הריכוזים המוחלטים ואת ערך הסף (cutoff) הנגזר מהם כספציפיים לפלטפורמה וכזמניים, כפי שנדון להלן.
במקביל, רמות SCC-Ag בסרום נקבעו באמצעות פלטפורמת Abbott Alinity i באמצעות בדיקת חיסון של חלקיקים מיקרוסקופיים עם פליטה כימית (CMIA), עם גבול ייחוס עליון של 1.5 ng/mL. כל ההליכים האנליטיים בוצעו בכפילות על ידי צוות מעבדה שהיה עיוור למצב המחלה ולכל המשתנים הקליניים, כולל השלב הקליני, הדרגה ההיסטולוגית ומצב ה-HPV. מקדמי השונות, הן בין-בדיקתיים (inter-assay) והן תוך-בדיקתיים (intra-assay), עבור המתודולוגיות שבהן נעשה שימוש נשמרו מתחת ל-10%.
ניתוח סטטיסטי
מאחר שריכוזי hWAPL בסרום לא עמדו בהנחת הנורמליות לפי מבחן Shapiro–Wilk (כל ה-P < 0.05), ניתוחים המשכיים בוצעו באמצעות שיטות לא-פרמטריות. לפיכך, משתנים רציפים מוצגים כחציונים עם טווחי בין-רבעוניים (IQRs). הבדלים בין שתי קבוצות בלתי תלויות הוערכו באמצעות מבחן Mann–Whitney U, בעוד שהשוואות בין שלוש קבוצות או יותר בוצעו באמצעות מבחן Kruskal–Wallis H. מתאמים בין משתנים הוערכו באמצעות ניתוח מתאם דרגות של Spearman עבור נתונים רציפים או סודרים, ומקדם מתאם point-biserial עבור משתנים בינאריים, כגון מצב זיהום ב-HPV. תוצאות הישרדות הושוו באמצעות מבחן log-rank, ויחסי סיכון (hazard ratios) עם רווחי סמך של 95% נאמדו מתוך מודלים של רגרסיית סיכונים פרופורציונליים של Cox.
הביצועים האבחנתיים הוערכו באמצעות ניתוח עקומות ROC, כאשר ערך הסף האופטימלי עבור כל סמן ביולוגי נקבע על ידי מקסום מדד Youden. הבדלים בשטח תחת העקומה (AUC) נבחנו באמצעות שיטת DeLong. כדי לבחון את התועלת האבחנתית המשולבת, נבנה מודל רגרסיה לוגיסטית בינארית שבו מצב המחלה (CSCC = 1; ביקורת = 0) שימש כמשתנה התלוי, ו-hWAPL בנסיוב ואנטיגן SCC שימשו כמנבאים בלתי תלויים. ההסתברויות החזויות שהתקבלו נשמרו כציון מורכב לצורך ניתוח ROC המשכי.
רגישות דיפרנציאלית לתרופות בין קבוצות WAPL-high ו-WAPL-low הוערכה באמצעות חבילת limma, תוך יישום ספי שער של adjusted P < 0.05 ו- |log₂ fold change| > 0.3. מובהקות סטטיסטית הוגדרה כ-P-value דו-זנבי < 0.05. כל הניתוחים בוצעו באמצעות SPSS, R, או הכלים המקוונים שצוינו. חבילות ה-R ששימשו למחקר כללו את survival, survminer, pROC, ggplot2, oncoPredict, ו-limma. עבור כל ניתוחי הממדים הגבוהים שבהם נבחנו היפותזות רבות בו-זמנית, דהיינו השוואות הפן-סרטן (pan-cancer), סריקת הביטוי המשותף (co-expression), ניתוחי הביטוי הדיפרנציאלי והעשרה תפקודית, וניתוח הרגישות לתרופות, ערכי ה-P-values הותאמו להשוואות מרובות באמצעות שיטת Benjamini–Hochberg false discovery rate, וערכי P-values מתוקננים מדווחים במקומות המתאימים; המתאמים הזוגיים בין hWAPL לבין גנים בודדים או פракציות של תאי מערכת החיסון מדווחות עם P-values נומינליים, ותוצאות מתאם אלו נחשבות כמייצרות היפותזות ולא כמאששות. מכיוון שמספר קטן של שדות קליניקו-פתולוגיים בקוהורט הציבורי לא הושלמו במלואם, מקרים עם נתונים חסרים הוצאו מהניתוח המתאים על בסיס זוגי במקום שיושלמו (imputation), ומספר המקרים שניתן להערכה מופיע עבור כל משתנה ב- Table 1. כדי להתייחס להבדלי הגיל בין קוהורט הסרום, בוצע ניתוח רגישות באמצעות רגרסיה לוגיסטית בינארית מותאמת לגיל, וחושב מתאם ספירמן (Spearman correlation) בין hWAPL בסרום לבין הגיל בתוך כל קבוצה.