High expression of NNMT in stromal CAFs is closely associated with tumor angiogenesis in OSCC tissues
In our previous study, we proved the high stromal expression pattern of NNMT in CAFs and its correlation with poor prognosis in OSCC patients [15]. To further determine the correlation between high stromal NNMT expression and tumor angiogenesis, we first analyzed the existing TCGA database using the Gene Expression Profiling Interactive Analysis (GEPIA) web tool. The results showed that, in patients with head and neck squamous cell carcinoma (HNSCC), the expression of NNMT was positively correlated with EC markers, PECAM1, and CD34 (Fig. 1a-b,). Additionally, we performed gene set enrichment analysis (GSEA) from TCGA-HNSC datasets, and the results showed that the upregulation of angiogenesis signatures (the gene set: HALLMARK) was strongly correlated with NNMT expression (Fig. 1c). Next, to calculate the microvessel density, we performed immunohistochemistry (IHC) to evaluate CD31 expression level in 54 samples from patients with stage I–IV oral cancer, Angio tool software was used to determine the MVD in the tissues (Fig. 1d-e, FigS1a). Next, IHC and the multiplex immunofluorescence staining assay were performed to analyze the NNMT protein levels and MVD in the OSCC tissue microarray. The results showed that NNMT expression was highly restricted to the stromal component of OSCC samples, and higher NNMT protein levels were significantly correlated with higher MVD (Fig. 1f-g, Fig.S1a). Overall, these results indicated that stromal NNMT expression was significantly correlated with tumor angiogenesis in OSCC.
NNMT regulates the pro-angiogenic phenotypes of fibroblasts in OSCC
To verify the contribution of NNMT to regulating the expression of pro-angiogenetic markers, immunoblotting assays were performed in CAFs stably infected with sh-Ctrl and sh-NNMT (Fig.S2a). The result demonstrated that the expression of VEGFA was significantly decreased in CAFs when NNMT was knocked down (Fig. 2a). Next, ELISA assays showed that the secretion of VEGFA was decreased in the culture supernatant of NNMT knockdown CAFs (Fig. 2b). To explore whether NNMT in CAFs could contribute to promoting angiogenesis in vitro, HUVEC was used. And matrigel tube formation assays indicated that the number of meshes formed by HUVEC were significantly decreased when incubated with the conditioned medium (CM) from NNMT knockdown CAFs (Fig. 2c-d). As in our previous study, we have verified that PFs, the potential precursor cells of CAFs, could be activated by NNMT overexpression [15]. To identify whether NNMT overexpression in PFs could contribute to angiogenesis, we conducted several experiments. The results showed that NNMT overexpression in PFs could help to upregulate VEGFA protein expression and VEGFA secretion in PFs (Fig. 2e-f). Matrigel tube formation assays indicated that the number of meshes formed by HUVEC were significantly increased when incubated with the CM from NNMT overexpression PFs (Fig. 2g-h).
In addition, to further clarify the function of stromal NNMT in promoting angiogenesis, we established a three-dimensional (3D) co-culture system of fibroblasts with HUVECs (Fig. 2i). (see methods for more details). In this 3D co-culture system, fibroblasts co-cultured with endothelial cells (ECs) could help ECs form lumen-like structure (Fig. 2i), which was further verified in the confocal scanning 3D model (Fig. 2j). Interestingly, the lumen formation ability of HUVEC was significantly decreased when co-cultured with NNMT knockdown CAFs (Fig. 2k, Fig.S2b). In addition, IF results showed knockdown of NNMT could decrease the number of ECs in the co-culture system (Fig. 2l-m). Taken together, these results indicated that stomal NNMT could promote the pro-angiogenic phenotypes of fibroblasts in OSCC in vitro.
Stromal NNMT accelerates angiogenesis by increasing VEGFA expression in vivo and in vitro
To clarify whether NNMT in fibroblasts contributes to tumor angiogenesis, an assembled organoid model and a co-implanted xenograft model, both integrating the patient-derived tumor organoids (PDOs) and corresponding fibroblasts, as well as HUVEC, were established for further experiments (Fig. 3a). To establish the assembled organoid model, we mixed the endothelial cells, fibroblasts, and organoid at a ratio of 1:5:5 to generate clusters, 24 hours later, clusters were mixed with matrigel and seeded into 24-well plates (Fig. 3b). As a result, the volume of assembled organoids integrated with NNMT overexpression PFs was significantly increased (Fig. 3c-d). In addition, the number of endothelial cells in assembled organoids integrated with NNMT overexpression PFs was significantly increased (Fig. 3e). More importantly, IF staining for CD31 expression in assembled organoids further confirmed these results (Fig. 3f, Fig.S3a). Next, a co-implanted xenograft model was established (Fig. 3a). We found that xenografts integrated with NNMT knockdown CAFs possessed worse tumor formation ability, while xenografts integrated with NNMT overexpression PFs exhibited enhanced tumor formation ability (Fig. 3g). Both the volume and weight of tumors integrated with NNMT knockdown CAFs were lower than those in the negative control group (Fig. 3h-i). As expected, the volume and weight of tumors of integrated with PFs with NNMT overexpression were higher than those in the negative control group (Fig. 3j-k). In addition, tumors with NNMT overexpression PFs displayed increased expression of CD31 and VEGFA. In contrast, expression levels of CD31 and VEGFA in tumors with NNMT knockdown CAFs were greatly suppressed (Fig.S3b-d). Overall, these results demonstrated that stromal NNMT accelerated angiogenesis by increasing VEGFA expression in vivo and in vitro, thus promoting tumor progression.
NNMT increases VEGFA expression through epigenetically regulating ETS2 expression
The above results indicated that NNMT could increase VEGFA expression, and VEGFA plays a crucial role in NNMT-induced tumor angiogenesis. To investigate the exact mechanism of NNMT increasing VEGFA expression, we bioinformatically analyzed the potential transcription factors that could bind to the promoter region of VEGFA (Fig. 4a-b). Then we evaluated the mRNA expression of these transcription factors in NNMT knockdown CAFs and their scramble nonsense transfectants. The results showed that only ETS2 showed low mRNA expression in NNMT knockdown cells, while others showed no significant differences (Fig. 4c). Consistently, ETS2 protein expression was significantly downregulated in NNMT knockdown CAFs, while significantly upregulated in NNMT overexpression CAFs (Fig. 4e). These results suggested that ETS2 might be involved in the regulation of VEGFA in fibroblasts by NNMT. To further confirm whether ETS2 directly regulates VEGFA, we analyzed the existing TCGA database using the GEPIA web tool. The results showed that the expression of ETS2 was positively correlated with VEGFA expression in HNSCC tissues (Fig.S4a). CAFs and PFs were then infected with ETS2 knockdown and overexpression lentiviruses and their scramble nonsense transfectants, respectively. And we found that VEGFA protein expression was significantly decreased in ETS2 knockdown CAFs, while increased in ETS2 overexpression PFs (Fig. 4f-g).
To further confirm that ETS2 serves as a transcription factor of VEGFA, we searched for the putative binding motif of the transcription factor ETS2 (Fig. 4h) and found two putative binding sites of ETS2 in the promoter region (between − 2000 and 0) of VEGFA from JASPAR website (http://jaspar.genereg.net) (Fig. 4i) [22]. To identify exact regulatory elements in the VEGFA promoter for ETS2, a series of deletions of the potential VEGFA promoter sequences were cloned into pGL3-basic vectors. The reduced promoter activity implies that the truncated region might encompass crucial binding sites for transcription factors. Transient transfections in 293T cells revealed the significance of the − 1159 to 0 bp 5' flanking sequence in conferring VEGFA promoter activity (Fig. 4j).
To further explore the effect of ETS2 on VEGFA promoter activation, we conducted transfections in 293T cells using a luciferase reporter system driven by the VEGFA promoter, with a mutation in the E1 binding region. The VEGFA promoter activity exhibited a significant increase when the wild-type construct was co-transfected with ETS2. In contrast, co-transfection of the mutated E1 construct did not result in a substantial increase in promoter activity. Therefore, overexpression of ETS2 upregulates VEGFA‐luciferase activity (Fig. 4k). Moreover, the ChIP assay showed that ETS2 could directly bind to the VEGFA promoter and the E1 site was significantly enriched in comparison with the E2 site. (Fig. 4l). These results demonstrate that ETS2 does act as a transcription factor of the VEGFA promoter. Moreover, ELISA assays showed that secretion of VEGFA was decreased in the culture supernatant of ETS2 knockdown CAFs (Fig. 4m). Matrigel tube formation assays indicated the number of meshes formed by HUVEC was significantly decreased when incubated with the CM from ETS2 knockdown CAFs (Fig. 4n, Fig.S4b). Interestingly, in the 3D co-culture model, the lumen formation ability of HUVEC co-cultured with ETS2 knockdown CAFs was significantly decreased compared to control (Fig.S4c-d). Collectively, These results demonstrated that ETS2 serves as a transcription factor of VEGFA and participates in NNMT-mediated tumor angiogenesis.
NNMT modulates the expression of ETS2 by suppressing the tri-methylation of H3K27
As NNMT drives gene expression changes in CAFs through attenuation of the S-adenosyl methionine (SAM): S-adenosyl homocysteine(SAH) ratio, which regulates the methylation of histones [16]. We performed targeted immunoblotting to detect relative levels of histone lysine methylation, and the results confirmed that H3K27 tri-methylation was upregulated upon knockdown of NNMT in CAFs of OSCC. In addition, H3K27 tri-methylation was downregulated upon NNMT overexpression in PFs. Treatment of NNMT knockdown CAFs with the EZH2 histone methyltransferase inhibitor DZNep [16] restored H3K27me3 expression (Fig. 5a). In addition, treatment of NNMT overexpression PFs with GSK-J1, an inhibitor of H3K27me3 demethylases [23], restored H3K27me3 expression (Fig. 5b). To Further explore whether NNMT regulates ETS2 expression through H3K27 trimethylation, we performed CHIP assays and found that NNMT activity reduced H3K27me3 occupancy at the promoter site of ETS2 (Fig. 5e-f). Therefore, these results demonstrated that NNMT modulates ETS2 expression by suppressing the tri-methylation of H3K27me3.
NNMT-ETS2-VEGFA signal axis orchestrates the pro-angiogenic properties of CAFs to promote tumor angiogenesis
As ETS2, which is under epigenetic regulation by NNMT, serves as a transcription factor of VEGFA, we speculated that ETS2 is involved in high NNMT expression-mediated tumor angiogenesis in CAFs of OSCC. To address this question, we used CAFs stably expressing sh-NNMT plus ETS2 overexpression. Immunoblotting results showed that knockdown of NNMT significantly decreased VEGFA expression, while overexpression of ETS2 dramatically increased VEGFA expression (Fig. 6a). ELISA assays also showed that the VEGFA secretion was decreased in the culture supernatant of CAFs with NNMT knockdown, While overexpression of ETS2 rescued VEGFA secretion (Fig. 6b). Moreover, the angiogenic effect of CMs from CAFs with sh-NNMT was attenuated, while the angiogenic function of sh-NNMT was reversed by overexpression of ETS2 (Fig. 6c-d). Consistently, HUVECs co-cultured with CAFs with NNMT knockdown exhibited a worse lumen formation ability, and this effect was recovered by ETS2 overexpression (Fig. 6e-f). In addition, IF staining showed knockdown of NNMT markedly decreased the number of ECs, while these results were reversed by ETS2 overexpression (Fig. 6g-h,). Interestingly, the volume of assembled organoids and the number of endothelial cells were significantly decreased when integrated with NNMT knockdown CAFs, whereas this effect was significantly rescued by overexpression of ETS2 (Fig. 6i, Fig.S5a-b). To assess the effect of ETS2 involved in high NNMT expression of CAFs mediated tumor angiogenesis in vivo, a co-implanted xenograft model integrated with CAFs, HUVEC, and organoid was applied, and tumor masses were harvested on day 28. The results showed knockdown of NNMT in CAFs strongly suppressed tumor volume and weight of xenograft, while overexpression of ETS2 abolished the inhibition of tumor formation ability of sh-NNMT (Fig. 6j-l). Furthermore, IF staining of tumor tissues showed that VEGFA and CD31 expression was decreased with sh-NNMT, but these results were reversed by ETS2 overexpression (Fig. 6m-o). Collectively, these data demonstrated that stromal NNMT could regulate VEGFA through epigenetically regulating the transcription of ETS2, thus promoting tumor angiogenesis of OSCC (Fig. 7).