
==== Front
Sci Rep
Sci Rep
Scientific Reports
2045-2322
Nature Publishing Group UK London

39232098
71616
10.1038/s41598-024-71616-8
Article
Semaphorin 3F (SEMA3F) influences patient survival in esophageal adenocarcinoma
Knipper Karl karl.knipper@uk-koeln.de

1
Lyu Su Ir 2
Jung Jin-On 1
Alich Niklas 3
Popp Felix C. 1
Schröder Wolfgang 1
Fuchs Hans F. 1
Bruns Christiane J. 1
Quaas Alexander 2
Nienhueser Henrik 3
Schmidt Thomas thomas.schmidt1@uk-koeln.de

1
1 https://ror.org/00rcxh774 grid.6190.e 0000 0000 8580 3777 Faculty of Medicine and University Hospital of Cologne, Department of General, Visceral and Cancer Surgery, University of Cologne, Cologne, Germany
2 https://ror.org/00rcxh774 grid.6190.e 0000 0000 8580 3777 Faculty of Medicine and University Hospital of Cologne, Institute of Pathology, University of Cologne, Cologne, Germany
3 https://ror.org/038t36y30 grid.7700.0 0000 0001 2190 4373 Department of General, Visceral and Transplant Surgery, University of Heidelberg, Heidelberg, Germany
4 9 2024
4 9 2024
2024
14 2058917 12 2023
29 8 2024
© The Author(s) 2024
2024
https://creativecommons.org/licenses/by/4.0/ Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.
In esophageal adenocarcinoma, the presence of lymph node metastases predicts patients' survival even after curative resection. Currently, there is no highly accurate marker for detecting the presence of lymph node metastasis. The SEMA3F/NRP2 axis was initially characterized in axon guidance and recent evidence has revealed its significant involvement in lymphangiogenesis, angiogenesis, and carcinogenesis. Hence, the objective of this study was to elucidate the roles of SEMA3F and its receptor NRP2 in esophageal adenocarcinoma. We conducted an immunohistochemical evaluation of SEMA3F and NRP2 protein expression in 776 patients with esophageal adenocarcinoma who underwent Ivor-Lewis esophagectomy at the University Hospital of Cologne. Total and positive cancer cell counts were digitally analyzed using QuPath and verified by experienced pathologists to ensure accuracy. Positive expression was determined as a cell percentage exceeding the 50th percentile threshold. In our cohort, patients exhibiting SEMA3F positive expression experience significantly lower pT- and pN-stages. In contrast, positive NRP2 expression is associated with the presence of lymph node metastases. Survival analyses showed that the expression status of NRP2 had no impact on patient survival. However, SEMA3F positivity was associated with a favorable patient survival outcome (median OS: 38.9 vs. 26.5 months). Furthermore, SEMA3F could be confirmed as an independent factor for better patient survival in patients with early tumor stage (pT1N0-3: HR = 0.505, p = 0.014, pT1-4N0: HR = 0.664, p = 0.024, pT1N0: HR = 0.483, p = 0.040). In summary, SEMA3F emerges as an independent predictor for a favorable prognosis in patients with early-stage esophageal adenocarcinoma. Additionally, NRP2 expression is linked to a higher risk of lymph node metastases occurrence. We hypothesize that low SEMA3F expression could identify patients with early-stage tumors who might benefit from more aggressive treatment options or intensified follow-up. Furthermore, SEMA3F and its associated pathways should be explored as potential tumor-suppressing agents.

Keywords

Esophageal adenocarcinoma
SEMA3F
Semaphorin
NRP2
Neuropilin
Lymph node metastasis
Subject terms

Tumour biomarkers
Oesophageal cancer
http://dx.doi.org/10.13039/501100017515 Heidelberger Stiftung Chirurgie Universitätsklinikum Köln (8977)Open Access funding enabled and organized by Projekt DEAL.

issue-copyright-statement© Springer Nature Limited 2024
==== Body
pmcIntroduction

The prognosis for patients with esophageal adenocarcinoma remains bleak despite extensive research efforts and multimodal therapy regimens, often accompanied by high morbidity rates1. Lymph node status is recognized as a critical prognostic indicator for patients with esophageal adenocarcinoma2. In this context, not only the number of metastasized lymph nodes is pivotal, but also the total number of lymph nodes resected and the lymph node ratio (the proportion of metastasized to resected lymph nodes) are significant predictors of patient survival3,4. Furthermore, the lymph node status at the initial staging serves as a parameter for clinical decision-making. Patients with low cT-stage and negative cN-stage qualify for sole endoscopic resection or primary tumor resection5. In endoscopic resection patients do not undergo node resection. Therefore, an accurate assessment of the lymph node status during staging is of a crucial importance. It must also be noted that patients with a low T-stage can still have metastasized lymph nodes6. Consequently, there is a pressing need for a biomarker capable of predicting lymph node status and additionally identifying high-risk individuals within patient groups typically considered low-risk.

Semaphorins are a family of proteins that were first described as important molecules for axonal guidance during neural development7. Since then, numerous subtypes have been identified across various organisms, including vertebrates, invertebrates, and viruses7. Semaphorin 3 refers to secreted molecules that bind to their receptors, neuropilin7,8.

Semaphorins have been detected in the development of various organs such as the parathyroid, thymus, and many others9,10. Furthermore, semaphorins have been implicated not only in physiological mechanisms but also in influencing cell proliferation, angiogenesis, and the development of metastases in cancer11. Semaphorin 3F (SEMA3F) plays a complex role in cancer progression that is not yet fully understood. Transcriptome profiling of hepatocellular carcinomas suggests that SEMA3F may promote metastasis development by activating focal adhesion pathways12. SEMA3F acts mainly through its receptor Neuropilin 2 (NRP2)13. Conversely, SEMA3F is characterized as a tumor suppressor gene in esophageal squamous cell carcinoma, where low mRNA expression levels of SEMA3F are associated with poor survival outcomes and a positive lymph node status. Furthermore, the SEMA3F expression correlated negatively with the NRP2 expression14. However, lymph node metastatic patterns seem to differ between esophageal adenocarcinoma and esophageal squamous cell carcinoma15. The role of SEMA3F and NRP2 in esophageal adenocarcinoma remains unclear. Consequently, the objective of this study was to assess NRP2 and its ligand SEMA3F as potential biomarkers for patient survival and their ability to predict lymph node metastases in esophageal adenocarcinoma.

Materials and methods

Patients and tumor samples

In this study, 776 patients were included. The inclusion criteria were esophageal adenocarcinoma, and Ivor-Lewis esophagectomy conducted with curative intention. Patients with survival of less than 30 days after surgery or insufficient tumor sample quality were excluded. All patients underwent surgery from 1998 until 2019 at the University Hospital of Cologne. Written informed consent was obtained from every patient. The study was approved by the Ethics Committee of the University Hospital of Cologne (ethics committee number: 21-1146) and was conducted in accordance with the declaration of Helsinki. The specimens’ pathological tumor ((y)pT) as well as lymph node status (following neoadjuvant therapy) ((y)pN), lymphatic invasion (L), vascular invasion (V), perineural invasion (Pn), and grading (G) were analyzed according to the 7th edition of the Union for International Cancer Control at our department of pathology. Tissue cylinders (1.2 mm) were punched out of each tumor sample with a semi-automated precision instrument and transferred to a paraffin-embedded tissue microarray (TMA). Stainings were conducted on 4 µm thick slices.

Immunohistochemistry

Immunohistochemical stainings and control stainings were performed according to the manufacturers' advice. Stainings were done for SEMA3F (HPA035008, Sigma-Aldrich, St. Louis, USA) and NRP2 (HPA039980, Sigma-Aldrich, St. Louis, USA). The automatic staining system Leica BOND-MAX (Leica Biosystems, Wetzlar, Germany) was used for all stainings. The stained slides were digitalized using Aperio GT 450 DX (Leica Biosystems, Wetzlar, Germany). The stainings were analyzed digitally with QuPath v0.3.2 as described before (Suppl. Fig. 1)16. The percentage of positive cells in tumor tissue was calculated using a positive cell detector and tissue classifier. Quality control was performed by two experienced pathologists (S.L. and A.Q.). The total study cohort was split up into two cohorts: negative and positive marker expression. Positivity was defined as values above the 50th percentile.

Statistical analysis

Statistical analyses were performed using IBM SPSS Statistics (Version 28.0.1.1). p values below 0.05 were considered statistically significant. Evaluation of qualitative value was performed with the Chi-Square test. Survival analyses were made with Kaplan-Meier curves and the log-rank test. Univariate and multivariate Cox regression analyses were used to detect interdependencies of clinicopathologic and survival data. Only clinicopathologic values, which received p values below 0.20 in univariate Cox regression analyses, were included in multivariate Cox regression analyses. All clinical data was collected prospectively and analyzed retrospectively. Overall survival was measured from date of surgery until death or loss of follow-up.

Results

All 776 included patients underwent Ivor-Lewis esophagectomy with curative intention. Clinicopathological values are depicted in Table 1. Of these patients, 685 (88.3%) were male. 70.6% (n = 548) of the included patients received perioperative/neoadjuvant therapy prior to resection. Lymph node metastases were detected in 458 patients (59.0%). SEMA3F and NRP2 immunohistochemical stainings were performed and analyzed digitally (Suppl. Fig. 1). The cohort was divided into two groups according to the negative or positive marker expression. SEMA3F-positive patients showed significantly lower pT-stages, lower pN-stages, fewer lymphovascular and perineural invasions (p(T) = 0.016, p(N) = 0.041, p(L) = 0.026, p(Pn) = 0.033, Table 1). On the contrary, patients with positive NRP2-stainings showed significantly more lymph node metastases, and more often lymphovascular invasion (p(N) = 0.041, p(L) = 0.002, Table 1). Table 1 General clinicopathological values of the total study population and patients with negative or positive SEMA3F as well as NRP2 marker expression.

Characteristic	Total	SEMA3F	p value	NRP2	p value	
	Negative	Positive	Negative	Positive	
n (%)	n (%)	n (%)	n (%)	n (%)	
No. of patients	776 (100)	388 (100)	388 (100)		388 (100)	388 (100)		
Sex				0.147			0.220	
 Male	685 (88.3)	336 (86.6)	349 (89.9)		348 (89.7)	337 (86.9)		
 Female	91 (11.7)	52 (13.4)	39 (10.1)		40 (10.3)	51 (13.1)		
Age				0.219			0.718	
 < 65	433 (55.8)	208 (53.6)	225 (58.0)		214 (55.2)	219 (56.4)		
 ≥ 65	343 (44.2)	180 (46.4)	163 (42.0)		174 (44.8)	169 (43.6)		
Median overall survival (months)	31.6	26.5	38.9		30.6	32.1		
(95% confidence interval)	(26.9–36.6)	(22.3–30.7)	(24.5–53.4)		(24.7–36.5)	(24.8–39.4)		
Perioperative/neoadjuvant therapy				0.753			0.636	
 No	228 (29.4)	112 (28.9)	116 (29.9)		111 (28.6)	117 (30.2)		
 Yes	548 (70.6)	276 (71.1)	272 (70.1)		277 (71.4)	271 (69.8)		
Perioperative/neoadjuvant				0.050			0.004	
treatment regime						
 CROSS	288 (29.6)	140 (36.1)	148 (38.1)		164 (42.3)	124 (32.0)		
 FLOT	92 (11.9)	38 (9.8)	54 (13.9)		46 (11.9)	46 (11.9)		
 Other	166 (21.4)	97 (25.0)	69 (17.8)		65 (16.8)	101 (26.0)		
 None	230 (29.6)	113 (29.1)	117 (30.2)		113 (29.1)	117 (30.2		
(y)pT				0.016			0.264	
 1	143 (18.4)	56 (14.4)	87 (22.4)		68 (17.5)	75 (19.3)		
 2	128 (16.5)	70 (18.0)	58 (14.9)		60 (15.5)	68 (17.5)		
 3	482 (62.1)	247 (63.7)	235 (60.6)		252 (64.9)	230 (59.3)		
 4	23 (3.0)	15 (3.9)	8 (2.1)		8 (2.1)	15 (3.9)		
(y)pN				0.041			0.041	
 0	318 (41.0)	145 (37.4)	173 (44.6)		173 (44.6)	145 (37.4)		
 1–3	458 (28.7)	243 (62.6)	215 (55.4)		215 (55.4)	243 (62.6)		
L				0.026			0.002	
 0	363 (46.8)	166 (42.8)	197 (50.8)		203 (52.3)	160 (41.2)		
 1	413 (53.2)	222 (57.2)	191 (49.2)		185 (47.7)	228 (58.8)		
V				0.110			0.161	
 0	583 (75.1)	290 (74.7)	293 (75.5)		303 (78.1)	280 (72.2)		
 1	77 (9.9)	32 (8.3)	45 (11.6)		34 (8.8)	43 (11.1)		
 2	116 (14.9)	66 (17.0)	50 (12.9)		51 (13.1)	65 (16.8)		
Pn				0.033			0.258	
 0	508 (65.5)	259 (66.8)	249 (64.2)		259 (66.8)	249 (64.2)		
 1	156 (20.1)	65 (16.8)	91 (23.5)		81 (20.9)	75 (19.3)		
 2	112 (14.4)	64 (16.5)	48 (12.4)		48 (12.4)	64 (16.5)		
G				0.995			0.762	
 1	2 (0.3)	1 (0.3)	1 (0.3)		1 (0.3)	1 (0.3)		
 2	118 (15.2)	58 (14.9)	60 (15.5)		56 (14.4)	62 (16.0)		
 3	103 (13.3)	50 (12.9)	53 (13.7)		54 (13.9)	49 (12.6)		
 4	0 (0.0)	0 (0.0)	0 (0.0)		0 (0.0)	0 (0.0)		
 Not applicable/unknown	553 (71.2)	279 (71.9)	274 (70.6)		277 (71.3)	276 (71.1)		
Bold print marks p values below 0.05.

CROSS Chemoradiotherapy for Oesophageal Cancer followed by Surgery Study, FLOT fluorouracil, leucovorin, oxaliplatin and docetaxel, NRP2 neuropilin 2, SEMA3F semaphorine 3F.

To evaluate the impact of SEMA3F and NRP2 expression on patients' survival, we carried out survival analyses. Here, we could not find an impact of NRP2 expression on overall survival (p = 0.820, Fig. 1A). Positive SEMA3F expression correlated with favorable patient survival (38.9 months vs. 26.5 months, p = 0.014, Fig. 1B). Since perioperative/neoadjuvant therapy could influence the proteome of tumors, we conducted our survival analyses depending on perioperative/neoadjuvant therapy or primary surgery. NRP2 expression proved to have no impact on patients' survival in both subgroups (p(perioperative/neoadjuvant) = 0.741, p(primary surgery) = 0.457, Fig. 1C and E). SEMA3F showed also no impact in patients after perioperative/neoadjuvant therapy (p = 0.703, Fig. 1D). However, patients, who received primary surgery and showed a positive SEMA3F expression, showed a significantly better overall survival compared to a negative expression (140.9 vs. 26.5 months, p < 0.001, Fig. 1F).Fig. 1 Kaplan-Meier curves for overall survival of the total cohort depending on (A) NRP2 expression (n(negative) = 388, n(positive) = 388, p = 0.820) and (B) SEMA3F expression (n(negative) = 388, n(positive) = 388, p = 0.014). Furthermore, Kaplan-Meier curves for overall survival for the perioperative/neoadjuvant treated subgroup depending on (C) NRP2 expression (n(negative) = 277, n(positive) = 271, p = 0.741) and (D) SEMA3F expression (n(negative) = 276, n(positive) = 272, p = 0.703) as well as for the primary resected subgroup depending on (E) NRP2 expression (n(negative) = 111, n(positive) = 117, p = 0.457) and (F) SEMA3F expression (n(negative) = 112, n(positive) = 116, p < 0.001).

Additionally, we conducted survival analyses in a patient subgroup, that received adjuvant therapy due to the perioperative FLOT treatment regime. Here, no significant survival difference could be detected depending on the SEMA3F or NRP2 expression (p(NRP2) = 0.737, p(SEMA3F) = 0.367). In line, patients, who received no adjuvant therapy as a primary therapy regime, also did not show a significant survival difference depending on the SEMA3F or NRP2 expression (p(NRP2) = 0.333, p(SEMA3F) = 0.065).

Following these analyses, we performed multivariate Cox regression analyses for the total cohort and the primary resected subgroup. Here, SEMA3F expression did not prove to be an independent factor for patients' survival (total cohort: HR = 0.785, 95% CI 0.555–1.111, p = 0.172, Table 2; primary surgery: HR = 0.747, 95% CI 0.498–1.121, p = 0.159, Table 3). Higher age, (y) pT-stage, (y) pN-stage, and lymphovascular invasion were shown to be independent risk factors for worse survival in our cohort (Tables 2 and 3). Female sex and vascular invasion are significant factors for favorable survival (Table 2). Table 2 Multivariate Cox regression analyses of the total study cohort.

Characteristic	Borders	Hazard ratio	95% confidence interval	p value	
Sex	Female vs male	0.469	0.261–0.841	0.011	
Age	≥ 65 vs < 65	1.473	1.011–2.146	0.044	
Perioperative/neoadjuvant therapy	Yes vs no	1.199	0.768–1.874	0.425	
(y)pT	≥ 2 vs 1	1.465	1.129–1.903	0.004	
(y)pN	≥ 1 vs 0	1.701	1.413–2.047	 < 0.001	
L	1 vs 0	1.951	1.293–2.943	0.001	
V	≥ 1 vs 0	0.563	0.372–0.773	 < 0.001	
G	≥ 2 vs 1	1.176	0.822–1.683	0.376	
SEMA3F	Positive vs negative	0.785	0.555–1.111	0.172	
Bold print marks p values below 0.05.

SEMA3F semaphorine 3F.

Table 3 Multivariate Cox regression analyses of the primarily resected subcohort.

Characteristic	Borders	Hazard ratio	95% confidence interval	p value	
Age	≥ 65 vs < 65	1.411	0.904–2.202	0.129	
(y)pT	≥ 2 vs 1	1.554	1.169–2.065	0.002	
(y)pN	≥ 1 vs 0	1.663	1.348–2.051	 < 0.001	
V	≥ 1 vs 0	0.877	0.677–1.137	0.322	
G	≥ 2 vs 1	1.237	0.827–1.849	0.300	
SEMA3F	Positive vs negative	0.747	0.498–1.121	0.159	
Bold print marks p values below 0.05.

SEMA3F semaphorine 3F.

Since patients with positive SEMA3F expression showed a favorable survival but SEMA3F expression was not an independent risk factor in the total cohort, we carried out further subgroup analyses depending on the (y)pT- and (y)pN-stage (pT1N0-3: Suppl. Table 1, pT1-4N0: Suppl. Table 2, pT1N0: Suppl. Table 3). Here, patients with positive NRP2 expression showed significantly more often lymph node metastases as well as lymphovascular invasion in the subgroup of patients with early (y)pT-stage (p((y)pN) = 0.024, p(L) = 0.019, Suppl. Table 1). No further differences in clinicopathological values depending on the NRP2 or SEMA3F expression could be detected. Additional survival analyses showed no correlation between NRP2 expression and patient survival (p(pT1N0-3) = 0.259, p(pT1-4N0) = 0.575, p(pT1N0) = 0.055, Fig. 2A–C). However, patients with positive SEMA3F expression showed significantly better overall survival compared to patients with negative SEMA3F expression in all three subgroups (p(pT1N0-3) = 0.007, p(pT1-4N0) = 0.012, p(pT1N0) = 0.010, Fig. 2D–F).Fig. 2 Kaplan-Meier curves for overall survival depending on NRP2 expression in the (A) pT1N0-3 (n(negative) = 68, n(positive) = 75, p = 0.259), (B) pT1-4N0 (n(negative) = 173, n(positive) = 145, p = 0.575), and (C) pT1N0 subgroup (n(negative) = 58, n(positive) = 52, p = 0.055). Furthermore, Kaplan-Meier curves for overall survival depending on SEMA3F expression in the (D) pT1N0-3 (n(negative) = 56, n(positive) = 87, p = 0.007), (E) pT1-4N0 (n(negative) = 145, n(positive) = 173, p = 0.012), and (F) pT1N0 subgroup (n(negative) = 45, n(positive) = 65, p = 0.010).

Moreover, multivariate Cox regression analyses were carried out. Here, SEMA3F expression could be proven as a significant, independent factor for favorable overall survival in these three subgroups (pT1N0-3: HR = 0.505, 95% CI 0.293–0.870, p = 0.014, Suppl. Table 4; pT1-4N0: HR = 0.664, 95% CI 0.465–0.947, p = 0.024, Table 4; pT1N0: HR = 0.483, 95% CI 0.242–0.967, p = 0.040, Suppl. Table 5). Furthermore, high patient age, perioperative/neoadjuvant therapy, and the occurrence of lymph node metastases were shown to be independent risk factors for worse patient survival (Table 4, Suppl. Tables 4–5). Table 4 Multivariate Cox regression analyses of the pT1-4N0 study cohort.

Characteristic	Borders	Hazard ratio	95% confidence interval	p value	
Sex	Female vs male	0.611	0.343–1.090	0.095	
Age	≥ 65 vs < 65	1.441	1.007–2.063	0.046	
Perioperative/neoadjuvant therapy	yes vs no	1.494	0.932–2.397	0.096	
(y)pT	≥ 2 vs 1	1.207	0.958–1.521	0.111	
V	≥ 1 vs 0	0.853	0.680–1.071	0.171	
SEMA3F	Positive vs negative	0.664	0.465–0.947	0.024	
Bold print marks p values below 0.05.

SEMA3F semaphorine 3F.

In summary, SEMA3F expression could be identified as an independent factor for favorable patient survival in patients with early-stage esophageal adenocarcinoma. Furthermore, positive NRP2 expression was associated with the occurrence of lymph node metastases as well as lymphovascular invasion.

Discussion

In this study, we investigated the role of SEMA3F and its receptor, NRP2, in patients diagnosed with esophageal adenocarcinoma. We observed a significant correlation between SEMA3F expression and patients' survival. Conversely, there was no observed correlation between NRP2 expression and patient survival in the same cohort. SEMA3F positivity was linked to improved patient survival (median OS: 38.9 months vs. 26.5 months), aligning with similar findings observed in esophageal and oral squamous cell carcinoma cases14,17. Furthermore, in-vitro experiments revealed an inhibitory effect of SEMA3F on cell migration and invasion in, for example, oral squamous cell cancer cells or breast cancer cells18,19. SEMA3F influences not only the tumor cells directly but also impacts angiogenesis. Here, SEMA3F expression resulted in decreased intratumoral microvessel density in lung cancer models in mice20. These effects of SEMA3F are mainly established through the inhibition of mTOR signaling13.

Interestingly, divergent roles of SEMA3F are described in various cancer entities. In hepatocellular carcinoma (HCC), a high SEMA3F expression could be correlated with an unfavorable survival21. This can be elucidated by the distinct macromolecules measured in the mentioned publications—mRNA versus protein—for SEMA3F. In our study, we measured protein expression using immunohistochemistry.

Additionally, SEMA3F is known to modify the tumor microenvironment22. High SEMA3F expression levels are positively correlated with infiltrating B cells, T cells, macrophages, and neutrophils in HCC21. The tumor microenvironment consistency is recognized for its variability across different cancer types23. Consequently, the influence of SEMA3F may change based on the specific cancer entity and the characteristics of the surrounding tumor microenvironment.

The impact of SEMA3F on patients' survival did not prove to be an independent factor in Cox regression analyses in our total patient cohort or in patients, treated with primary surgery. In subgroup analyses of patients with early tumor stage (pT1N0-3, pT1-4N0, pT1N0), positive SEMA3F expression correlated with better patient survival. Additionally, SEMA3F proves to be an independent factor for a favorable prognosis in these patient subgroups (p(pT1N0-3) = 0.014, p(pT1-4N0) = 0.024, p(pT1N0) = 0.040). Interestingly, SEMA3F does not have an impact on patients' survival in more advanced tumor stages. This could imply an impact on the carcinogenesis of SEMA3F in early tumor stages. Confirming this thesis, SEMA3F expression gets lost during cancer progression in head and neck squamous cell carcinoma (HNSCC)24. Inheriting this, SEMA3F negativity could serve as a marker for a high-risk patient cohort in—according to classically used risk definitions—low-risk patient cohorts.

SEMA3F positivity showed a significant correlation with improved survival outcomes in patients who underwent primary surgery, in contrast to those who received multimodal therapy including perioperative/neoadjuvant treatment. This trend was evident both in the overall cohort and among patients diagnosed with early-stage tumors. Additional fundamental studies are required to elucidate the underlying pathway behind this phenomenon.

In our patient cohort, high SEMA3F expression was correlated with a less frequent occurrence of lymph node metastases in esophageal adenocarcinoma (p = 0.041). In line with our data, high SEMA3F expression is correlated with a lower risk of occult lymph node metastases in HNSCC25. Furthermore, patients with high SEMA3F expressions did not only show fewer lymph node metastases but also less distant metastases in colorectal cancer26. Interestingly, we could show that patients with positive NRP2 expression have significantly more often lymph node metastases (p = 0.041). Similar effects could be observed in HNSCC25. In-vitro analyses showed, that SEMA3F predominantly interacts with lymphatic endothelial cells through its receptor NRP224. Lymphadenectomy is recommended for patients undergoing surgical resection5. However, the extensiveness of mediastinal lymphadenectomy correlates with the occurrence of resection-depending morbidity27. Clinical trials to minimize morbidity while maximizing diagnostic certainty of complete resection of all tumor cells including lymph nodes are ongoing and recruiting. In the realm of breast cancer and melanoma treatment, the practice of sentinel lymph node resection has demonstrated promise when applied to patients with esophageal cancer28. Despite its promising potential, this workflow has not yet been integrated into routine clinical practice and is not recommended by current guidelines5. Therefore, biomarkers for the occurrence of lymph node metastases and thus markers for the determination of the individually required extent of lymphadenectomy are needed.

SEMA3F could not only be used as a biomarker for patients' survival but also possibly in therapy planning. For example, the chemosensitivity to fluorouracil of human colorectal cancer cells was investigated in vitro and in vivo. Cells with higher SEMA3F expression showed higher chemosensitivity and a higher rate of cancer cell apoptosis29. Moreover, rats injected with lung cancer cells expressing SEMA3F exhibited significantly prolonged survival compared to rats injected with control lung cancer cells30. Vice versa the knockdown of SEMA3F reduced the apoptotic activity of rituximab in diffuse large B-cell lymphoma (DLBCLC)31. Based on this SEMA3F expression could help to modify the treatment selection.

In summary, our study has identified SEMA3F as a factor associated with favorable patient survival, particularly in those with early-stage esophageal adenocarcinomas. Furthermore, NRP2 expression could be correlated with a higher rate of lymph node metastases. Although retrospective in design, our findings have the potential to greatly inform future clinical decision-making, subject to validation through additional fundamental and prospective clinical studies. Low SEMA3F expression may identify a higher-risk subgroup within the traditionally defined low-risk patient category. This subgroup could benefit from more aggressive treatment options or intensified post-treatment care.

Conclusions

We demonstrated that positive SEMA3F protein expression is correlated with favorable patient survival in esophageal adenocarcinoma. Moreover, SEMA3F could be identified as an independent factor for better patient survival in patients with early tumor stage (pT1N0-3, pT1-4N0, pT1N0). The SEMA3F receptor NRP2 is associated with a higher likelihood of lymph node metastases. This highlights the necessity for further prospective clinical research to investigate whether negative SEMA3F expression and positive NRP2 expression could identify a subset of patients who are at higher risk. Such markers hold the potential to recalibrate risk assessment in patients who are currently categorized as low-risk by conventional pathological classifications, thereby necessitating more tailored monitoring and management strategies. These patients could potentially benefit from intensified follow-up, more aggressive systematic therapy options, or even extended lymphadenectomy. Additionally, the SEMA3F/NRP2 pathway may present novel therapeutic avenues if subsequent foundational studies corroborate its purported tumor-suppressive capabilities.

Supplementary Information

Supplementary Information.

Supplementary Information

The online version contains supplementary material available at 10.1038/s41598-024-71616-8.

Author contributions

K.K., H.N., A.Q. and T.S. contributed to the study's conception and design. Material preparation and data collection were performed by S.L. Analysis was performed by K.K. The first draft of the manuscript was written by K.K. All authors commented on previous versions of the manuscript. All authors read and approved the final manuscript.

Funding

Open Access funding enabled and organized by Projekt DEAL. This project received funding from the “Heidelberger Stiftung Chirurgie”.

Data availability

The datasets generated and analyzed during the current study are available from the corresponding author on reasonable request.

Competing interests

The authors declare no competing interests.

Ethical approval

This study was performed in line with the principles of the Declaration of Helsinki. Approval was granted by the Ethics Committee of the University of Cologne.

Consent to participate

Informed consent was obtained from all individual participants included in the study.

Publisher's note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

These authors contributed equally: Karl Knipper and Su Ir Lyu.
==== Refs
References

1. Eyck BM van Lanschot JJB Hulshof M van der Wilk BJ Shapiro J van Hagen P van Berge Henegouwen MI Wijnhoven BPL van Laarhoven HWM Nieuwenhuijzen GAP Hospers GAP Bonenkamp JJ Cuesta MA Blaisse RJB Busch OR Creemers GM Punt CJA Plukker JTM Verheul HMW Spillenaar Bilgen EJ van der Sangen MJC Rozema T Ten Kate FJW Beukema JC Piet AHM van Rij CM Reinders JG Tilanus HW Steyerberg EW van der Gaast A CROSS Study Group Ten-year outcome of neoadjuvant chemoradiotherapy plus surgery for esophageal cancer: The randomized controlled CROSS trial J. Clin. Oncol. 2021 39 18 1995 2004 10.1200/JCO.20.03614 33891478
Eyck, B. M. et al. Ten-year outcome of neoadjuvant chemoradiotherapy plus surgery for esophageal cancer: The randomized controlled CROSS trial. J. Clin. Oncol. 39(18), 1995–2004. 10.1200/JCO.20.03614 (2021).33891478 10.1200/JCO.20.03614
2. Lyu SI Popp FC Simon AG Schultheis AM Zander T Fretter C Schroder W Bruns CJ Schmidt T Quaas A Knipper K Copy-number-gain of telomerase reverse transcriptase (hTERT) is associated with an unfavorable prognosis in esophageal adenocarcinoma Sci. Rep. 2023 13 1 17699 10.1038/s41598-023-44844-7 37848472
Lyu, S. I. et al. Copy-number-gain of telomerase reverse transcriptase (hTERT) is associated with an unfavorable prognosis in esophageal adenocarcinoma. Sci. Rep. 13(1), 17699. 10.1038/s41598-023-44844-7 (2023).37848472 10.1038/s41598-023-44844-7
3. Akutsu Y Matsubara H The significance of lymph node status as a prognostic factor for esophageal cancer Surg. Today 2011 41 9 1190 1195 10.1007/s00595-011-4542-y 21874413
Akutsu, Y. & Matsubara, H. The significance of lymph node status as a prognostic factor for esophageal cancer. Surg. Today 41(9), 1190–1195. 10.1007/s00595-011-4542-y (2011).21874413 10.1007/s00595-011-4542-y
4. Henckens SPG Hagens ERC van Berge Henegouwen MI Meijer SL Eshuis WJ Gisbertz SS Impact of increasing lymph node yield on staging, morbidity and survival after esophagectomy for esophageal adenocarcinoma Eur. J. Surg. Oncol. 2023 49 1 89 96 10.1016/j.ejso.2022.07.012 35933270
Henckens, S. P. G. et al. Impact of increasing lymph node yield on staging, morbidity and survival after esophagectomy for esophageal adenocarcinoma. Eur. J. Surg. Oncol. 49(1), 89–96. 10.1016/j.ejso.2022.07.012 (2023).35933270 10.1016/j.ejso.2022.07.012
5. Obermannova R Alsina M Cervantes A Leong T Lordick F Nilsson M van Grieken NCT Vogel A Smyth EC ESMO Guidelines Committee Oesophageal cancer: ESMO Clinical Practice Guideline for diagnosis, treatment and follow-up Ann. Oncol. 2022 10.1016/j.annonc.2022.07.003 35914638
Obermannova, R. et al. Oesophageal cancer: ESMO Clinical Practice Guideline for diagnosis, treatment and follow-up. Ann. Oncol.10.1016/j.annonc.2022.07.003 (2022).35914638 10.1016/j.annonc.2022.07.003
6. Ancona E Rampado S Cassaro M Battaglia G Ruol A Castoro C Portale G Cavallin F Rugge M Prediction of lymph node status in superficial esophageal carcinoma Ann. Surg. Oncol. 2008 15 11 3278 3288 10.1245/s10434-008-0065-1 18726651
Ancona, E. et al. Prediction of lymph node status in superficial esophageal carcinoma. Ann. Surg. Oncol. 15(11), 3278–3288. 10.1245/s10434-008-0065-1 (2008).18726651 10.1245/s10434-008-0065-1
7. Alto LT Terman JR Semaphorins and their signaling mechanisms Methods Mol. Biol. 2017 1493 1 25 10.1007/978-1-4939-6448-2_1 27787839
Alto, L. T. & Terman, J. R. Semaphorins and their signaling mechanisms. Methods Mol. Biol. 1493, 1–25. 10.1007/978-1-4939-6448-2_1 (2017).27787839 10.1007/978-1-4939-6448-2_1
8. He Z Tessier-Lavigne M Neuropilin is a receptor for the axonal chemorepellent Semaphorin III Cell 1997 90 4 739 751 10.1016/s0092-8674(00)80534-6 9288753
He, Z. & Tessier-Lavigne, M. Neuropilin is a receptor for the axonal chemorepellent Semaphorin III. Cell 90(4), 739–751. 10.1016/s0092-8674(00)80534-6 (1997).9288753 10.1016/s0092-8674(00)80534-6
9. Takahashi K Ishida M Hirokawa K Takahashi H Expression of the semaphorins Sema 3D and Sema 3F in the developing parathyroid and thymus Dev. Dyn. 2008 237 6 1699 1708 10.1002/dvdy.21556 18489001
Takahashi, K., Ishida, M., Hirokawa, K. & Takahashi, H. Expression of the semaphorins Sema 3D and Sema 3F in the developing parathyroid and thymus. Dev. Dyn. 237(6), 1699–1708. 10.1002/dvdy.21556 (2008).18489001 10.1002/dvdy.21556
10. Hinck L The versatile roles of “axon guidance” cues in tissue morphogenesis Dev. Cell 2004 7 6 783 793 10.1016/j.devcel.2004.11.002 15572123
Hinck, L. The versatile roles of “axon guidance” cues in tissue morphogenesis. Dev. Cell 7(6), 783–793. 10.1016/j.devcel.2004.11.002 (2004).15572123 10.1016/j.devcel.2004.11.002
11. Bica C Tirpe A Nutu A Ciocan C Chira S Gurzau ES Braicu C Berindan-Neagoe I Emerging roles and mechanisms of semaphorins activity in cancer Life Sci. 2023 318 121499 10.1016/j.lfs.2023.121499 36775114
Bica, C. et al. Emerging roles and mechanisms of semaphorins activity in cancer. Life Sci. 318, 121499. 10.1016/j.lfs.2023.121499 (2023).36775114 10.1016/j.lfs.2023.121499
12. Ye K Ouyang X Wang Z Yao L Zhang G SEMA3F promotes liver hepatocellular carcinoma metastasis by activating focal adhesion pathway DNA Cell Biol. 2020 39 3 474 483 10.1089/dna.2019.4904 31968181
Ye, K., Ouyang, X., Wang, Z., Yao, L. & Zhang, G. SEMA3F promotes liver hepatocellular carcinoma metastasis by activating focal adhesion pathway. DNA Cell Biol. 39(3), 474–483. 10.1089/dna.2019.4904 (2020).31968181 10.1089/dna.2019.4904
13. Nakayama H Bruneau S Kochupurakkal N Coma S Briscoe DM Klagsbrun M Regulation of mTOR signaling by semaphorin 3F-neuropilin 2 interactions in vitro and in vivo Sci. Rep. 2015 5 11789 10.1038/srep11789 26156437
Nakayama, H. et al. Regulation of mTOR signaling by semaphorin 3F-neuropilin 2 interactions in vitro and in vivo. Sci. Rep. 5, 11789. 10.1038/srep11789 (2015).26156437 10.1038/srep11789
14. Xie Z Li T Huang B Liu S Zhang L Zhang Q Semaphorin 3F serves as a tumor suppressor in esophageal squamous cell carcinoma and is associated with lymph node metastasis in disease progression Technol. Cancer Res. Treat. 2020 19 1533033820928117 10.1177/1533033820928117 32441221
Xie, Z. et al. Semaphorin 3F serves as a tumor suppressor in esophageal squamous cell carcinoma and is associated with lymph node metastasis in disease progression. Technol. Cancer Res. Treat. 19, 1533033820928117. 10.1177/1533033820928117 (2020).32441221 10.1177/1533033820928117
15. Ma Z Chen C Shang X Yue J Jiang H Comparison of lymph node metastasis pattern from esophagogastric junction adenocarcinoma versus very low thoracic esophageal squamous cancer: A propensity-matched analysis J. Thorac. Dis. 2023 15 2 442 451 10.21037/jtd-22-1028 36910061
Ma, Z., Chen, C., Shang, X., Yue, J. & Jiang, H. Comparison of lymph node metastasis pattern from esophagogastric junction adenocarcinoma versus very low thoracic esophageal squamous cancer: A propensity-matched analysis. J. Thorac. Dis. 15(2), 442–451. 10.21037/jtd-22-1028 (2023).36910061 10.21037/jtd-22-1028
16. Bankhead P Loughrey MB Fernandez JA Dombrowski Y McArt DG Dunne PD McQuaid S Gray RT Murray LJ Coleman HG James JA Salto-Tellez M Hamilton PW QuPath: Open source software for digital pathology image analysis Sci. Rep. 2017 7 1 16878 10.1038/s41598-017-17204-5 29203879
Bankhead, P. et al. QuPath: Open source software for digital pathology image analysis. Sci. Rep. 7(1), 16878. 10.1038/s41598-017-17204-5 (2017).29203879 10.1038/s41598-017-17204-5
17. Zhang B Gao Z Sun M Li H Fan H Chen D Zheng J Prognostic significance of VEGF-C, semaphorin 3F, and neuropilin-2 expression in oral squamous cell carcinomas and their relationship with lymphangiogenesis J. Surg. Oncol. 2015 111 4 382 388 10.1002/jso.23842 25475162
Zhang, B. et al. Prognostic significance of VEGF-C, semaphorin 3F, and neuropilin-2 expression in oral squamous cell carcinomas and their relationship with lymphangiogenesis. J. Surg. Oncol. 111(4), 382–388. 10.1002/jso.23842 (2015).25475162 10.1002/jso.23842
18. Liu Y Li R Yin K Ren G Zhang Y The crucial role of SEMA3F in suppressing the progression of oral squamous cell carcinoma Cell Mol. Biol. Lett. 2017 22 32 10.1186/s11658-017-0064-y 29299034
Liu, Y., Li, R., Yin, K., Ren, G. & Zhang, Y. The crucial role of SEMA3F in suppressing the progression of oral squamous cell carcinoma. Cell Mol. Biol. Lett. 22, 32. 10.1186/s11658-017-0064-y (2017).29299034 10.1186/s11658-017-0064-y
19. Nasarre P Kusy S Constantin B Castellani V Drabkin HA Bagnard D Roche J Semaphorin SEMA3F has a repulsing activity on breast cancer cells and inhibits E-cadherin-mediated cell adhesion Neoplasia 2005 7 2 180 189 10.1593/neo.04481 15802023
Nasarre, P. et al. Semaphorin SEMA3F has a repulsing activity on breast cancer cells and inhibits E-cadherin-mediated cell adhesion. Neoplasia 7(2), 180–189. 10.1593/neo.04481 (2005).15802023 10.1593/neo.04481
20. Potiron VA Sharma G Nasarre P Clarhaut JA Augustin HG Gemmill RM Roche J Drabkin HA Semaphorin SEMA3F affects multiple signaling pathways in lung cancer cells Cancer Res. 2007 67 18 8708 8715 10.1158/0008-5472.CAN-06-3612 17875711
Potiron, V. A. et al. Semaphorin SEMA3F affects multiple signaling pathways in lung cancer cells. Cancer Res. 67(18), 8708–8715. 10.1158/0008-5472.CAN-06-3612 (2007).17875711 10.1158/0008-5472.CAN-06-3612
21. Lou W Wang W Chen J Wang S Huang Y ncRNAs-mediated high expression of SEMA3F correlates with poor prognosis and tumor immune infiltration of hepatocellular carcinoma Mol. Ther. Nucleic Acids 2021 24 845 855 10.1016/j.omtn.2021.03.014 34026328
Lou, W., Wang, W., Chen, J., Wang, S. & Huang, Y. ncRNAs-mediated high expression of SEMA3F correlates with poor prognosis and tumor immune infiltration of hepatocellular carcinoma. Mol. Ther. Nucleic Acids 24, 845–855. 10.1016/j.omtn.2021.03.014 (2021).34026328 10.1016/j.omtn.2021.03.014
22. Nakayama H Kusumoto C Nakahara M Fujiwara A Higashiyama S Semaphorin 3F and netrin-1: The novel function as a regulator of tumor microenvironment Front. Physiol. 2018 9 1662 10.3389/fphys.2018.01662 30532711
Nakayama, H., Kusumoto, C., Nakahara, M., Fujiwara, A. & Higashiyama, S. Semaphorin 3F and netrin-1: The novel function as a regulator of tumor microenvironment. Front. Physiol. 9, 1662. 10.3389/fphys.2018.01662 (2018).30532711 10.3389/fphys.2018.01662
23. de Visser KE Joyce JA The evolving tumor microenvironment: From cancer initiation to metastatic outgrowth Cancer Cell 2023 41 3 374 403 10.1016/j.ccell.2023.02.016 36917948
de Visser, K. E. & Joyce, J. A. The evolving tumor microenvironment: From cancer initiation to metastatic outgrowth. Cancer Cell 41(3), 374–403. 10.1016/j.ccell.2023.02.016 (2023).36917948 10.1016/j.ccell.2023.02.016
24. Doci CL Mikelis CM Lionakis MS Molinolo AA Gutkind JS Genetic identification of SEMA3F as an antilymphangiogenic metastasis suppressor gene in head and neck squamous carcinoma Cancer Res. 2015 75 14 2937 2948 10.1158/0008-5472.CAN-14-3121 25952650
Doci, C. L., Mikelis, C. M., Lionakis, M. S., Molinolo, A. A. & Gutkind, J. S. Genetic identification of SEMA3F as an antilymphangiogenic metastasis suppressor gene in head and neck squamous carcinoma. Cancer Res. 75(14), 2937–2948. 10.1158/0008-5472.CAN-14-3121 (2015).25952650 10.1158/0008-5472.CAN-14-3121
25. Meler-Claramonte C Aviles-Jurado FX Vilaseca I Terra X Bragado P Fuster G Leon Vintro X Camacho M Semaphorin-3F/neuropilin-2 transcriptional expression as a predictive biomarker of occult lymph node metastases in HNSCC Cancers 2022 14 9 2259 10.3390/cancers14092259 35565388
Meler-Claramonte, C. et al. Semaphorin-3F/neuropilin-2 transcriptional expression as a predictive biomarker of occult lymph node metastases in HNSCC. Cancers 14(9), 2259. 10.3390/cancers14092259 (2022).35565388 10.3390/cancers14092259
26. Zhou ZH Rao J Yang J Wu F Tan J Xu SL Ding Y Zhan N Hu XG Cui YH Zhang X Dong W Liu XD Bian XW SEMA3F prevents metastasis of colorectal cancer by PI3K-AKT-dependent down-regulation of the ASCL2-CXCR4 axis J. Pathol. 2015 236 4 467 478 10.1002/path.4541 25866254
Zhou, Z. H. et al. SEMA3F prevents metastasis of colorectal cancer by PI3K-AKT-dependent down-regulation of the ASCL2-CXCR4 axis. J. Pathol. 236(4), 467–478. 10.1002/path.4541 (2015).25866254 10.1002/path.4541
27. D'Journo XB Doddoli C Michelet P Loundou A Trousse D Giudicelli R Fuentes PA Thomas PA Transthoracic esophagectomy for adenocarcinoma of the oesophagus: Standard versus extended two-field mediastinal lymphadenectomy? Eur. J. Cardiothorac. Surg. 2005 27 4 697 704 10.1016/j.ejcts.2004.12.022 15784377
D’Journo, X. B. et al. Transthoracic esophagectomy for adenocarcinoma of the oesophagus: Standard versus extended two-field mediastinal lymphadenectomy?. Eur. J. Cardiothorac. Surg. 27(4), 697–704. 10.1016/j.ejcts.2004.12.022 (2005).15784377 10.1016/j.ejcts.2004.12.022
28. Nagaraja V Eslick GD Cox MR Sentinel lymph node in oesophageal cancer-a systematic review and meta-analysis J. Gastrointest. Oncol. 2014 5 2 127 141 10.3978/j.issn.2078-6891.2014.008 24772341
Nagaraja, V., Eslick, G. D. & Cox, M. R. Sentinel lymph node in oesophageal cancer-a systematic review and meta-analysis. J. Gastrointest. Oncol. 5(2), 127–141. 10.3978/j.issn.2078-6891.2014.008 (2014).24772341 10.3978/j.issn.2078-6891.2014.008
29. Tao M Ma H Fu X Wang C Li Y Hu X Lv R Zhou G Wang J Liu R Zhou M Xu G Wang Z Qin X Long Y Huang Q Chen M Zhou Q Semaphorin 3F induces colorectal cancer cell chemosensitivity by promoting P27 nuclear export Front. Oncol. 2022 12 899927 10.3389/fonc.2022.899927 36119535
Tao, M. et al. Semaphorin 3F induces colorectal cancer cell chemosensitivity by promoting P27 nuclear export. Front. Oncol. 12, 899927. 10.3389/fonc.2022.899927 (2022).36119535 10.3389/fonc.2022.899927
30. Kusy S Nasarre P Chan D Potiron V Meyronet D Gemmill RM Constantin B Drabkin HA Roche J Selective suppression of in vivo tumorigenicity by semaphorin SEMA3F in lung cancer cells Neoplasia 2005 7 5 457 465 10.1593/neo.04721 15967098
Kusy, S. et al. Selective suppression of in vivo tumorigenicity by semaphorin SEMA3F in lung cancer cells. Neoplasia 7(5), 457–465. 10.1593/neo.04721 (2005).15967098 10.1593/neo.04721
31. Li Q Ma N Li X Yang C Zhang W Xiong J Zhu L Li J Wen Q Gao L Yang C Rao L Gao L Zhang X Rao J Reverse effect of Semaphorin-3F on rituximab resistance in diffuse large B-cell lymphoma via the Hippo pathway Chin. Med. J. (Engl.) 2023 136 12 1448 1458 10.1097/CM9.0000000000002686 37114652
Li, Q. et al. Reverse effect of Semaphorin-3F on rituximab resistance in diffuse large B-cell lymphoma via the Hippo pathway. Chin. Med. J. (Engl.) 136(12), 1448–1458. 10.1097/CM9.0000000000002686 (2023).37114652 10.1097/CM9.0000000000002686
