VEGF and SEMA4D have synergistic effects on the promotion of angiogenesis in epithelial ovarian cancer
Cellular & Molecular Biology Lettersvolume 23, Article number: 2 (2018)
Anti-angiogenesis therapy that targets VEGF is one of the important treatment strategies in advanced ovarian cancer. However, depending on the pharmaceutical agent, treatment can have undesirable side effects. SEMA4D has recently gained interest for its role in promoting angiogenesis. Here, we try to further understand the mechanism by which SEMA4D promotes angiogenesis in ovarian cancer.
Correlation and western blot assaya were used to detect the relationship between VEGF and SEMA4D in clinical tissues and cells. Vasculogenic mimicry and transwell migration analyses were used to detect the roles of VEGF, SEMA4D and plexin-B1 on vasculogenic mimicry and migration. Vascular density and SEMA4D expression was determined using immunofluorescence staining in clinical tissues of EOC. Western blot was used to detect the expressions of CD31, MMP2 and VE-cadherin. We also analyzed the relationship between VEGF-SEMA4D and malignant tumor prognosis.
We found that knockdown of VEGF could suppress SEMA4D expression and that the expressions of VEGF and SEMA4D have a positive correlation in EOC cancer tissues. Vasculogenic mimicry and transwell migration analyses showed that SEMA4D and VEGF have a synergistic effect on the promotion of angiogenesis in A2780 and HUVEC cells. Soluble SEMA4D (sSEMA4D) could promote VM and migration in A2780 and HUVEC cells via the SEMA4D/plexin-B1 pathway, but the effect was not noted in stably transfected shR-plexin-B1 cells. In clinical tissues of EOC, the vascular density and SEMA4D/plexin-B1 expression were higher. When VEGF, SEMA4D and plexin-B1 was knocked down, the expression of CD31, MMP2 and VE-cadherin, which are the markers and initiators of angiogenesis and the epithelial–mesenchymal transition (EMT) process were reduced. VEGF and SEMA4D had a positive correlation with the malignant degree of ovarian cancer, and SEMA4D can serve as an independent prognostic factor.
VEGF and SEMA4D have synergistic effects on the promotion of angiogenesis in epithelial ovarian cancer. Targeting VEGF and the SEMA4D signaling pathway could be important for the therapy for EOC.
Epithelial ovarian cancer (EOC) is the most common malignant gynecological tumor. The EOC mortality rate ranks first among the gynecological malignancies and the overall 5-year survival rate is only ~30% . As with most malignant solid tumors, the rich microvascular environment encourages the growth of EOC and the prognosis is generally poor. The formation of large numbers of microvessels is the basis for the growth and metastasis of ovarian cancer in general and EOC in particular .
Vascular endothelial growth factor (VEGF) plays a critical role in the formation of vessels under both physiological and pathological conditions. VEGF increases the formation of nascent blood vessels by promoting the proliferation and migration of endothelial cells, increasing vascular permeability, and inhibiting endothelial cell apoptosis [3, 4].
Angiogenesis is the process of generating new blood vessels from existing blood vessels. In 1971, Folkman proposed that the growth and metastasis of tumors are dependent on tumor angiogenesis . Vascular growth has two different phases: slow and rapid. Without angiogenesis, the primary tumor does not grow. Tumor blood vessels provide the necessary oxygen and nutrients to the tumor tissue, allowing it to grow rapidly, and also facilitate distant metastases [6,7,8]. Inhibiting angiogenesis is a relatively new anti-tumor strategy . The role of tumor angiogenesis in the development of cancer has also attracted the attention of scholars both at home and abroad .
Anti-angiogenesis therapy, mainly targeting anti-VEGF, is an important treatment strategy in numerous malignancies . The targeted therapy bevacizumab (Avastin) inhibits tumor progression by binding to VEGF and blocking its biological activity. However, its clinical effect is not satisfactory. For example, Phase III clinical trials GOG218  and ICON7  showed that chemotherapy plus bevacizumab did not significantly improve overall survival for patients with ovarian cancer after first-line treatment despite the slight increase in progression-free survival. Fischer et al.  suggest that once patients with advanced cancer are administered anti-angiogenic therapy, the expression of other pro-angiogenic factors (e.g., SDF-1, FGF-1, FGF-2, etc.) is significantly increased, eventually leading to the failure of anti-VEGF treatment. Therefore, it is important to explore the molecular mechanism of anti-VEGF treatment failure to improve the efficacy of treatments for epithelial ovarian cancer.
The semaphorin superfamily is a class of proteins sharing a common Sema region. One semaphorin, axon guidance protein, which is part of the SEMA family, participates in axon guidance, spindle formation and nerve signal transduction pathways. It was initially discovered in the embryonic nervous system. Recent studies indicate that members of the SEMA family are also involved in the process of tumor angiogenesis [15,16,17].
The SEMA IV subfamily, a subcategory of the SEMA family, promotes angiogenesis. One important member is semaphorin4D (SEMA4D), which mediates endothelial cell chemotaxis and promotes the process of microvascular formation and distant metastasis in head and neck squamous cell carcinoma . Basile et al.  reported that membrane-type matrix metalloproteinase 1 (MT1-MMP) is secreted into the tumor microenvironment and acts in concert with SEMA4D to promote tumor angiogenesis, thus also promoting distant metastasis. Plexins are the only membrane-bound receptors of SEMA4D. Plexin-B1 combines with SEMA4D secreted into the microenvironment from tumor cells to promote tumor angiogenesis and endothelial cell migration and vessel formation [18,19,20].
Here, we mainly explore the role of SEMA4D in the process of angiogenesis in epithelial ovarian cancer. VEGF and SEMA4D expression exhibit a positive correlation in clinical tissues. Vasculogenic mimicry (VM) and transwell migration analyses showed the synergistic promotional effect of SEMA4D and VEGF on angiogenesis in A2780 and HUVEC cells. Soluble SEMA4D (sSEMA4D) promotes VM and migration in A2780 and HUVEC cells via the SEMA4D/plexin-B1 pathway. This effect was not noted in stably transfected shR-plexin-B1 cells. In EOC tissues, where the vascular network is denser, increased SEMA4D/plexin-B1 expression can be observed. When VEGF, SEMA4D and plexin-B1 were knocked down, the expressions of CD31, MMP2 and VE-Cadherin, which are the markers and initiators of angiogenesis and epithelial–mesenchymal transition (EMT), were reduced. VEGF and SEMA4D were positively correlated with the malignant degree of ovarian cancer. SEMA4D can serve as an independent prognostic factor.
In conclusion, SEMA4D may be an important molecule for anti-angiogenesis therapy, and further development of anti-SEMA4D therapy may give rise to a significant breakthrough in the field of EOC therapy.
Cell lines and cell culture
Human umbilical vein endothelial cells (HUVECs, ScienCell, Cat. #8000) were cultured in endothelial cell medium (ECM, Cat. #1001). The human ovarian cancer cell line A2780 was cultured in DMEM as the basal medium with 10% fetal bovine serum (FBS) and double antibiotics added for cell culture.
Transwell migration analysis
The transwell migration assay was performed in a Boyden chamber. Serum-free medium containing sSEMA4D and VEGF separately or in combination was placed in the bottom layer as a chemotactic. A solution of 0.1% BSA and 10% FBS was used as a control. HUVECs transfected with shRNA were placed in the upper layer. The control group was used to determine changes in HUVEC migration ability in each group.
Vasculogenic mimicry analysis
Vasculogenic tube formation was tested using a commercial Matrigel assay kit (BD Biosciences). We coated 24-well culture plates with Matrigel matrix (0.1 ml/well; BD Biosciences) and allowed it to solidify at 37 °C for 30 min. HUVECs and A2780 cells were added to the plates 24 h after transfection. The cells were dissociated into single cells and resuspended at 6 × 104 cells/ml in endothelial basal medium containing 2% fetal calf serum. The cells in each group were then plated at 0.5 ml/well onto the Matrigel surface and incubated at 37 °C. Tube networks were quantified as the total number of pixels (×200) in thresholded images using Image Pro Plus 6.0 software. A total of 3 wells were randomly selected.
Total RNA was reverse transcribed into cDNA using a reverse transcription kit (Bao Biological Engineering (Dalian) Co., Ltd.) and quantitative real-time PCR was performed on a Bio-Rad CFX96 instrument. The following solutions were successively added into the PCR mixture (20 μl total): 1 μl of cDNA, 10 μl of QIAGEN SYBR Green, 5 μl of RNase-free H2O, and 2 μl of upstream and downstream primers. The PCR conditions and procedures were: pre-denaturation at 95 °C for 5 min, 40 cycles of denaturation at 95 °C for 10 s and annealing at 60 °C for 30 s. β-actin levels were used to normalize expression.
Cells were lysed in lysis buffer using the RIPA method. The lysate was subject to SDS-PAGE, electrically transferred onto nitrocellulose membrane and incubated at 4 °C overnight in primary antibody. Membranes were washed with TBST and incubated in secondary antibody for 1 h. The membrane was then washed with TBST and exposed.
Human EOC tissues and adjacent noncancerous tissues were collected from the Tianjin Medical University Cancer Institute and Hospital. Informed consent was obtained from each patient. The Ethics Committee of Tianjin Medical University granted ethical approval. Immunohistochemistry was performed according to previously described methods. The sections were pre-treated using microwave irradiation, blocked, and incubated with polyclonal antibodies. The staining intensity was then assessed.
Statistical analysis was done using SPSS 18.0 statistical software. Data were represented as means and standard deviation. p < 0.05 was considered statistically significant.
All experiments were performed in accordance with relevant guidelines and regulations.
Results and discussion
Positive correlation between VEGF and SEMA4D expression
To explore the change in SEMA4D expression when VEGF expression is inhibited in ovarian cancer cells, we constructed a VEGF-knockdown plasmid and screened the stably transfected cells. Quantitative RT-PCR and western blot were used to verify the effectiveness of transfection and stability of the cells (Fig. 1a). Figure 1b shows that the expressions of VEGF and SEMA4D were significantly reduced in these cells compared with the expressions in A2780 cells or HUVECs (by approximately 50% and 40%, respectively).
To explore the relationship between VEGF and SEMA4D in A2780 cells and HUVECs, 124 samples of epithelial ovarian cancer tissues and 20 samples of normal ovarian tissues were subjected to immunohistochemical examination. As shown in Fig. 1c and Table 1, the expression levels of SEMA4D and VEGF in epithelial ovarian cancer tissues were significantly higher than those in normal ovarian tissue. We also applied correlation analysis to the expressions of SEMA4D and VEGF in epithelial ovarian cancer tissues. The results revealed a positive correlation between VEGF and SEMA4D (Table 2).
Vasculogenic mimicry promoted by the synergistic action of SEMA4D and VEGF
To study the role of SEMA4D in vasculogenic mimicry (VM), cells stably transfected with siR-SEMA4D were constructed and subjected to VM analysis. Figure 2a shows that vessel formation decreased by approximately 40% or 60% when SEMA4D was knocked down in A2780 cells or HUVECs, respectively. We prepared sSEMA4D significantly promoted VM of HUVECs (Fig. 2b).
When HUVECs and A2780 cells were co-transfected with vectors to overexpress SEMA4D and VEGF, this largely suppressed vessel formation (Fig. 2c). When sSEMA4D and VEGF were added to cells stably transfected with We also constructed VEGFR2-shRNA and plexin-B1-shRNA, no suppression recovery was seen in VM (Fig. 2d). These results indicate that SEMA4D and VEGF synergistically promote angiogenesis in HUVEC cells and in cells of the epithelial ovarian cancer line A2780.
SEMA4D and VEGF synergistically promote migration
We further explored the effects of SEMA4D and VEGF on migration, using the transwell migration assay to assess migration ability. As shown in Fig. 3a, migration decreased by approximately 30% or 45% when SEMA4D was knocked down in A2780 cells or HUVECs, respectively.
SSEMA4D were found to significantly promote migration of HUVECs (Fig. 3b). In addition, co-transfection of SEMA4D and VEGF into HUVECs and A2780 cells largely suppressed migration (Fig. 3c). Subsequently, we added sSEMA4D and VEGF to cells stably transfected with VEGFR2-shRNA and plexin-B1-shRNA but found this did not rescue the suppression of VEGFR and plexin-B1 on migration (Fig. 3d). These results indicate that SEMA4D and VEGF synergistically promote the migration ability of epithelial ovarian cancer cells.
To confirm the role of SEMA4D in angiogenesis, we applied immunofluorescence staining to examine the location of SEMA4D in ovarian cancer tissues and normal tissues. We also examined the endothelial cell marker CD31, which targets vascular endothelial cells. Immunofluorescence staining showed that the vascular density in ovarian cancer tissues increased significantly, and that the expression of SEMA4D considerably increased in the endothelial cells in blood vessels. SEMA4D was found to co-localizes with CD31 (Fig. 4a). We also detected the expression and location of plexin-B1 in clinical ovarian cancer tissues and normal tissues. The results indicate that plexin-B1 localizes in the area of endothelial cells (Fig. 4b). SEMA4D and plexin-B1 expression levels were considerably higher in ovarian cancer tissues than in normal tissues.
The molecular mechanism of VEGF–SEMA4D promotion in angiogenesis
We examined the expressions of the markers of angiogenesis. When VEGF or SEMA4D were inhibited, CD31 expression reduced by 61.1% and 52.8%, respectively. We also found that CD31 expression was inhibited significantly when plexin-B1 expression was inhibited (Fig. 4c). When VEGF, SEMA4D and plexin-B1 were knocked down, the expressions of MMP2 and VE-cadherin, which are involved in the initiation of angiogenesis, were reduced. Finally, when VEGF, SEMA4D and plexin-B1 were knocked down, the expression of E-cadherin reduced, while that of vimentin increased (Fig. 4c). These results indicate that VEGF, SEMA4D and plexin-B1 could inhibit EMT, which could also explain their inhibition of migration.
Relationship of VEGF–SEMA4D to malignancy and prognosis
We further analyzed the correlation between SEMA4D–VEGF and the malignancy of clinical samples. We used the chi-square test to demonstrate that SEMA4D and VEGF expression levels increase in epithelial ovarian carcinomas with a higher degree of malignancy or clinical aggressiveness (Table 3). Univariate analysis indicated that the overall survival and progression-free survival of epithelial ovarian cancer patients with positive expression of SEMA4D were lower than those for epithelial ovarian cancer patients with negative SEMA4D expression (Table 4 and Fig. 5). A multivariate proportional hazards model indicated that SEMA4D is an independent prognostic factor based on the analysis of epithelial ovarian cancer patient prognosis.
VEGF is an important molecule in ovarian cancer and its expression is negatively correlated with prognosis. Orre et al. found that the vascular density in ovarian malignant tumors is significantly higher than that in borderline and benign tumors. In ovarian clear cell carcinoma, the expression of VEGF is significantly correlated with vascular density . Through immunohistochemical analysis of 110 cases of epithelial ovarian cancer, Yamamo et al. found that the VEGF-positive rate in poorly differentiated ovarian cancer specimens was significantly higher than that of well-differentiated ovarian cancer and benign ovarian tumors. ELISA assessment revealed that VEGF levels were significantly higher in the serum of ovarian cancer patients with malignant tumors than in those with borderline and benign tumors. After tumor resection, VEGF levels decreased significantly, whereas tumor recurrence increased again, suggesting that VEGF levels in the blood could be used as prognostic indicators of ovarian cancer .
Here, we found that knockdown of VEGF expression led to the inhibition of VM and migration. We also used VEGF to detect the formation of VM and migration in ovarian cancer cells and HUVECs.
Despite the importance of VEGF inhibition in clinical trials, the outcome was not perfect. For example, GOG218  and ICON7  phase III clinical trials showed that chemotherapy plus bevacizumab gave no significant improvement in overall survival for ovarian cancer patients after first-line treatment despite a slight reduction in progression-free survival. Thus, some pro-angiogenic activity may be increased during anti-VEGF therapy. In this study, we found that SEMA4D may participate in the failure of anti-VEGF therapy.
SEMA4D, a member of the SEMA family within the semaphorin superfamily, plays an important role in the nervous system, the immune system, angiogenesis, and tumor invasion and metastasis. We found that SEMA4D is highly expressed in colon, breast and lung cancer tissues and plays an important role in the development of squamous cell carcinoma of the head and neck and the induction of angiogenesis . In the tumor microenvironment, SEMA4D expression plays a critical role in tumor angiogenesis and vascular maturation and increases the tumorigenicity of tumor cells. In cells not expressing SEMA4D, tumorigenicity and tumor metastasis are seriously weakened .
We confirmed that SEMA4D promotes the angiogenesis of ovarian cancer cells. SEMA4D is also known to bind to its receptor plexin-B1 to induce and promote tumor angiogenesis via the Met signaling pathway . Here, we found that SEMA4D promotes tumor VM and migration, partly explaining the promotional effects of SEMA4D on angiogenesis. Moreover, we found that the addition of sSEMA4D to the medium of ovarian cancer cells and HUVECs promotes the formation of VM and migration. This ability was largely inhibited when the receptors of VEGF and SEMA4D were knocked down. Furthermore, we found that VEGF and SEMA4D were significantly correlated with the malignant degree of ovarian cancer through data analysis of clinical specimens. SEMA4D can be used as an independent prognostic factor.
A variety of genes and bio-processes are involved in angiogenesis, including CD31, VE-cadherin, MMP2 and EMT. CD31 is a bio-marker for endothelial cells. VE-cadherin is a transmembrane cadherin protein located in the endothelial cell adherence region. It can maintain the vascular permeability and stability of the intravascular environment, which is a starting element in angiogenesis. Hendrix et al.  found that VE-cadherin promotes the formation of VM in highly invasive melanoma cells.
The MMP family of proteins plays an important role in remodeling the pericellular environment and degrading the extracellular matrix (ECM) . Cell behaviors such as cell proliferation, adhesion, migration, invasion, apoptosis and host defense are all related to the MMPs . MMP2 is necessary for classical angiogenesis . The EMT process could mimic the onset of angiogenesis, which is an important molecular mechanism.
Here, we found that knockdown of VEGF, SEMA4D and plexin-B1 could suppress the expression of CD31, VE-cadherin and MMP2, and the process of EMT which indicated that VEGF, SEMA4D and plexin-B1 could promote the onset of angiogenesis.
Our research aims to define the role of SEMA4D in anti-VEGF failure. We have shown that SEMA4D coordinates with VEGF during angiogenesis via plexin-B1 in EOC. This finding is significant for improving our understanding of chemotherapy resistance and the prognosis for patients with advanced epithelial ovarian cancer, for whom traditional means of treatment have proven ineffective.
VEGF and SEMA4D have synergistic effects on the promotion of angiogenesis in epithelial ovarian cancer. This provides a reference for therapy that should be investigated further.
Bovine serum albumin
Platelet endothelial cell adhesion molecule-1
Epithelial ovarian cancer
Fetal bovine serum
human umbilical vein endothelial cells
Matrix metalloprotein 2
Matrix metalloprotein 14
Vascular endothelial cadherin
Vascular endothelial growth factor
Bookman MA. Should studies of maintenance therapy be maintained in women with ovarian cancer? J Gynecol Oncol. 2013;24:105–7.
Folkman J. Incipient angiogenesis. Natl Cancer Inst. 2000;92:94–5.
Gianni-Barrera R, Trani M, Reginato S, Banfi A. To sprout or to split? VEGF, notch and vascular morphogenesis. Biochem Soc Trans. 2011;39:1644–8. https://doi.org/10.1042/BST20110650.
Scartozzi M, et al. VEGF and VEGFR genotyping in the prediction of clinical outcome for HCC patients receiving sorafenib: the ALICE-1 study. Int J Cancer. 2014;135:1247–56.
Folkman J, et al. Isolation of a tumor factor responsible for angiogenesis. J ExpMed. 1971;133(2):275–88.
Folkman J. Tumor angiogenesis: therapeutic implications. N Engl J Med. 1971;285(21):1182–6.
Carmeliet P, Jain RK. Angiogenesis in cancer and other diseases. Nature. 2000;407(6801):249–57.
Carmeliet P, Baes M. Metabolism and therapeutic angiogenesis. N Engl J Med. 2008;358(23):2511–2.
Burstein HJ, Schwartz RS. Molecular origins of cancer. N Engl J Med. 2008;358(5):527.
Avail-Lundqvist EH, et al. The impact of tumour angiogenesis, p53 overexpression and proliferative activity (M1B-1) on survival in squamous cervical carcinoma. Eur J Cancer. 1997;33(11):1799–804.
Sitohy B, Nagy JA, Dvorak HF. Anti-VEGF/VEGFR therapy for cancer: reassessing the target. Cancer Res. 2012;72:1909–14.
Burger RA, Brady MF, Bookman MA. Phase III trial of bevacizumab (BEV) in the primary treatment of advanced epithelial ovarian cancer (EOC), primary peritoneal cancer (PPC), or fallopian tube cancer (FTC): a gynecologic oncology study. J Clin Oncol. 2010;28(Suppl):abstr LBA1.
Kristensen G, Perren T, Qian W. Result of interim analysis of overall survival in the GCIG ICON7 phase III randomized trial of bevacizumab in women with newly diagnosed ovarian cancer. J Clin Oncol. 2011;29(Suppl):abstr LBA5006.
Fischer C, Jonckx B, Mazzone M, Zacchigna S, Loges S, Pattarini L, Chorianopoulos E, Liesenborghs L, Koch M, De Mol M, Autiero M, Wyns S, Plaisance S, Moons L, van Rooijen N, Giacca M, Stassen JM, Dewerchin M, Collen D, Carmeliet P. Anti-PlGF inhibits growth of VEGF(R)-inhibitor-resistant tumors without affecting healthy vessels. Cell. 2007;131:463–75.
Kong-Beltran M, Stamos J, Wickramasinghe D. The Sema domain of met is necessary for receptor dimerization and activation. Cancer Cell. 2004;6:75–84.
Serini G, Maione F, Giraudo E, Bussolino F. Semaphorins and tumor angiogenesis. Angiogenesis. 2009;12:187–93.
Gu C, Giraudo E. The role of semaphorins and their receptors in vascular development and cancer. Exp Cell Res. 2013. [Epub ahead of print].
Basile JR, Castilho RM, Williams VP, Gutkind JS. Semaphorin 4D provides a link between axon guidance processes and tumor-induced angiogenesis. Proc Natl Acad Sci U S A. 2006;103:9017–22.
Basile JR, Barac A, Zhu T, Guan KL, Gutkind JS. Class IV semaphorins promote angiogenesis by stimulating rho-initiated pathways through plexin-B. Cancer Res. 2004;64:5212–24.
Conrotto P, Valdembri D, Corso S, Serini G, Tamagnone L, Comoglio PM, Bussolino F, Giordano S. Sema4D induces angiogenesis through met recruitment by Plexin B1. Blood. 2005;105:4321–9.
Orre M, Rogers PA. VEGF, VEGFR-1, VEGFR-2, microvessel density and endothelial cell proliferation in tumours of the ovary. Int J Cancer. 1999;84(2):101–8.
Takeshita N, Hasegawa M, Sasaki K. In vivo expression and regulation of genes associated with vascularization during early response of sutures to tensile force. J Bone Miner Metab. 2016;29.
Maniotis AJ, Folberg R, Hess A, Seftor EA, Gardner LM, Pe’er J, Trent JM, Meltzer PS, Hendrix MJ. Vascular channel formation by human melanoma cells in vivo and in vitro: vasculogenic mimicry. Am J Pathol. 1999;155:739–52.
Sierra JR, Corso S, Caione L, et al. Tumor angiogenesis and progression are enhanced by Sema4D produced by tumor-associated macrophages. J Exp Med. 2008;205:1673–85.
Sun T, Zhao N, Zhao XL, Gu Q, Zhang SW, Che N, Wang XH, Du J, Liu YX, Sun BC. Expression and functional significance of Twist1 in hepatocellular carcinoma: its role in vasculogenic mimicry. Hepatology. 2010;51(2):545–56.
This study was funded by Tianjin Health Bureau of Science and Technology Funds (2012KZ073) and the National Natural Science Foundation (81302250).
Availability of data and materials
The data and materials involved in this report are available to readers.
Ying Chen and Lei Zhang are from Tianjin Medical University Cancer Institute and Hospital, Tianjin Key Laboratory of Cancer Prevention and Therapy and Tianjin National Clinical Research Centre of Cancer,
Wen-xin Liu and Ke Wang are from Tianjin Medical University Cancer Institute and Hospital.
Ethics approval and consent to participate
Human EOC tissues and adjacent non-cancerous tissues were collected from patients at the Tianjin Medical University Cancer Institute and Hospital. Informed consent was obtained from each patient. The Ethics Committee of Tianjin Medical University granted ethics approval (Code of Ethics: TMUhMEC2015021). All patients were cleared and consented to participate according to the parameters defined in ‘Ethics, consent and permissions’ and rights were given to publish individual patient data by each participant .
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.