Cell phenotypic switching alters expression of TRP channels in healthy cells
Decidualization of ESCs entails drastic remodeling of the fibroblast phenotype into epithelioid-like secretory cells. Decidual cells secrete several factors important for embryonic growth, such as prolactin (PRL) and insulin-like growth factor-binding protein 1 (IGFBP-1). To induce the MET-like decidualization process, healthy primary ESCs were supplemented with 8-Br-cAMP and MPA for an incubation period of five days. Afterwards, the decidualization of ESC was verified by the significant upregulation of PRL mRNA expression compared to vehicle stimulated ESCs (supplementary Fig. 2A). Next, the expression pattern of several TRP channels was investigated in control and decidualized ESC. We only assessed expression of TRP channels that were shown in previous work to be expressed in the endometrium, namely TRPV2, TRPV4, TRPC1, TRPC4, TRPC6, TRPM4 and TRPM7 (supplementary Fig. 1). These results showed a significant downregulation of the mRNA expression levels of TRPV2, TRPC1 and TRPC6 after decidualization, whereas expression levels of TRPV4, TRPC4, and TRPM7 were unaltered (Fig. 1A). TRPV2 was most abundantly expressed in ESC and underwent the most pronounced reduction in expression levels during decidualization (±10-fold). Interestingly, the reduced levels of TRPV2 mRNA were both time- and 8-Br-cAMP dose-dependent (supplementary Fig. 2B and C), underlining a strong association with the fibroblast phenotype. Next, the functional expression of TRPV2 was investigated using calcium microfluorimetry. Stimulation of primary ESCs by the TRPV2 agonist, THC (50µM) with co-application of CB receptor blockers, induced a robust calcium influx (D[Ca2+]i = 239.5 ± 56.9 nM (mean ± SEM)) in 34.4 ± 6.0 % of the cells. In contrast, application of THC in decidualized ESC resulted in significantly diminished calcium influx ((D[Ca2+]i = 161 ± 37.61 nM) and a lower number of responding cells (~11.8 ± 4.2 %) (Fig 1B-D). Finally, decreased TRPV2 protein expression after decidualization was validated on western blot (supplementary Fig. 2D). These findings suggest that during MET-related processes like decidualization, the TRP channel signature undergoes alterations, both at the mRNA and functional level.
Next, the TRP channel expression was investigated in EEC undergoing EMT. To induce EMT, EEC were incubated EMT-inducing medium supplement for 5 days. Successful EMT was verified by increased expression of mesenchymal markers MMP2, MMP9, ZEB1, ACTA, CDH2 and VIM, in combination with decreased expression of epithelial markers CDH1, KRT18, EPCAM and MMP7 (supplementary Fig. 3). Interestingly, inducing EMT in EEC was associated with significant upregulation of TRPV2 mRNA levels (4-fold), whereas expression levels of TRPM4 and TRPV6 were significantly downregulated. The expression levels of TRPV4 and TRPM7 were unaltered during the cell phenotype switch (Fig. 2A). The increased mRNA expression of TRPV2 was further validated at the functional level via calcium imaging experiments. These experiments showed the significantly enhanced TRPV2 functionality, demonstrated by increased percentage of THC-responding cells compared to vehicle treated cells (17.2 ± 3.6% vs. 1.5 ± 0.4%, respectively) and an increased calcium amplitude (Fig. 2B-D). Taken together, these results suggest a cell phenotypic TRP channel expression signature, where the mesenchymal phenotype is characterized by increased expression of TRPV2, TRPC1, TRPC4, while the epithelial phenotype is characterized by a pronounced expression of TRPM4 and TRPV6.
TRP channel expression in primary and metastatic endometrial cancer positively correlates with mesenchymal marker gene expression
Phenotypic switching is a malignant hallmark of endometrial cancer pathophysiology. EEC acquire mesenchymal characteristics, resulting in a more invasive cell phenotype, which, together with tumor microenvironment, contributes to disease progression and dissemination. Given the strong association between cell phenotype and the typical TRP channel expression, the representation of TRP channels in both primary and metastatic EC biopsies was investigated.
First, the mRNA expression of TRP channels, as well as of mesenchymal and epithelial markers was assessed in 54 patient biopsies collected from the primary tumor site via RT-qPCR. High expression of TRPM4 and TRPM7 was detected in primary tumor biopsies, while TRPV4, TRPC4, TRPC6 and especially TRPV6 expression was limited, while the expression of TRPC1 and TRPV2 are moderately expressed (Fig. 3A). Next, correlations between the expression levels of the different TRP channels and mesenchymal and epithelial markers was investigated. Expression of TRPV2, TRPC1, TRPC4, TRPC6 and TRPM7 significantly correlated with mesenchymal markers (Fig. 3B). Interestingly, multiple linear regression analysis identified that the gene expression of TRPV2, TRPC1, TRPC4, TRPC6 and TRPM7 was significantly related to the mesenchymal marker genes present in the model. Notably, expression levels of TRPV2 and TRPC1 were most strongly correlated with mesenchymal marker genes (R2 = 0.703 and R2 = 0.713, respectively) and were further considered to identify significant contributors to predict their expression (supplementary table 1). Interestingly, MMP2 and MMP9 expression were significant contributors in the model for predicting TRPV2 expression, whereas MMP7, CTSB, MMP9 and ZEB1 were significantly contributing to predict TRPC1 expression (supplementary table 2).
Next, we aimed to evaluate the expression of TRP channels and marker genes in tumor biopsies obtained from metastatic sites compared to primary tumor biopsies. Metastatic biopsies were acquired from 11 of the 54 patients analyzed above, allowing direct comparison of gene expression between primary and metastatic tumors from the same patient. As expected, metastatic biopsies exhibited significantly decreased expression of epithelial marker genes and showed an increased expression of typical mesenchymal genes, suggesting an increased EMT status (Fig. 4A). Univariate linear regression revealed a significant negative association of KRT18 and EPCAM expression, and a significant positive association of ZEB1 and ACTA2 expression with metastatic biopsies (Fig 4B). Interestingly, expression levels of TRPV2, TRPC1 and TRPM7 were significantly increased in the metastatic biopsy compared to the primary biopsy (Fig. 4C and D).
Similar, but more pronounced, correlation patterns were observed in metastatic biopsies (Fig. 4E). As such, TRPV2 and TRPC1 expression was positively correlated with mesenchymal marker gene expression and negatively correlated with epithelial marker gene expression. Contrary, TRPM4 expression correlated positively with epithelial marker genes and negatively with mesenchymal marker genes. In accordance with the results obtained in primary biopsies, a high correlation for TRPV2 and TRPC1 with mesenchymal gene expression was observed (R2 = 0.867 and R2 = 0.728, respectively) (supplementary table 3). In metastatic biopsies, ZEB1 was the significant predictor for both TRPV2 and TRPC1 in the regression model (supplementary table 4).
Relationship between TRP channel expression and clinical parameters
Next, the potential correlation between the variations in TRP channel expression levels and typical marker genes was assessed to explain the pathophysiological characteristics of the tumor biopsy. First, expression of epithelial and mesenchymal marker genes was assessed between several clinical parameters, including FIGO stage, grade, tumor histology, molecular classification, steroid dependency, presence of myometrial and lymph vascular invasion (LVI), and disease recurrence (Fig. 5A and B). Patient statistics are presented in supplementary Fig 4. Median expression values were used to construct a mesenchymal/epithelial (M/E) ratio for each biopsy, a value indicative for high-risk EC and (Fig. 5C). Despite the big overlap in the 95% CI intervals between the groups, a trend towards and increased M/E ratio was observed in late stage (from IB onwards) and high-grade tumors, as well as in tumors with serous and carcinosarcoma histology and p53 abnormal molecular classification (abnormal p53 expression, high risk). Furthermore, an increased M/E ratio was observed in tumors that expressed steroid receptors, tumors that invade up to 50% of the myometrium and displayed LVI. Finally, patients who experienced recurrence showed an elevated M/E ratio in the primary tumor biopsy.
Subsequently, univariate linear regression was used to identify significant associations between epithelial and mesenchymal marker genes and clinical parameters (Fig. 5D-H). Late-stage tumors showed increased expression of CDH2, ACTA2, MMP2 and ZEB1 compared to stage IA tumors (Fig. 5D). Reduced expression of CDH1 and KRT18 was significantly associated with high grade tumors and non-endometrioid tumor histology, while increased expression CTSB, CDH2 and ZEB1 could be observed in serous tumors (Fig. 5E and F). No significant associations were observed between marker genes and molecular subclasses, nor when comparing steroid receptor positive and negative tumors and presence of LVI. Increased expression of CDH2, CTSB, ACTA2, MMP2 and ZEB1 was associated with increased myometrial invasion, and expression of CTSB and MMP9 was significantly associated with tumor recurrence (Fig. 5G and H). These results suggest that high-risk EC tumors display an increased EMT status.
Interestingly, evaluating TRP channel expression signatures across clinical parameters reveals that TRPV2 and TRPC1 expression is associated with high-risk EC biopsies characterized by a high EMT status, while TRPM4 expression is associated with low-risk EC biopsies characterized by a low EMT status (Fig. 6A). Indeed, increased TRPV2 expression was significantly associated with stage IB, stage III and stage IV tumors compared to stage IA tumors (Fig. 6B), deep myometrial invasion (Fig. 6G) and recurrence (Fig. 6H). Expression of TRPC1 was significantly associated with high grade tumors (Fig. 6C), tumors with serous and carcinosarcoma histology (Fig. 6D) and p53 abnormal (abnormal p53 expression, high risk) classified tumors (Fig. 6E). In contrast, high TRPM4 expression was significantly associated with low grade tumors (Fig. 6C), endometrioid histology (Fig. 6D) and lower risk p53 WT (wild type p53 expression) and MMRd tumors (mismatch repair deficient) (Fig. 6E).
TRPV2 expression in primary biopsies as predictor for recurrence
Interestingly, previous results showed that recurrence of the disease was significant associated with increased levels of CTSB and MMP9, and TRPV2 expression (Fig. 5H and 6H). Thus, univariate logistic regression with recurrence as dependent variable and gene expression as independent variable was used to assess whether these genes in primary tumors would be able to predict recurrence. Remarkably, CTSB, MMP9 and TRPV2 were able to significantly predict disease recurrence. Moreover, tumor grade was also able to significantly predict disease recurrence (supplementary table 6). As such, ROC curves showed that TRPV2 and tumor grade were the best predictors for recurrence in this data set, with AUC values of 0.707 and 0.708 respectively (Fig. 7). These parameters were further used in a multivariate logistic regression model to predict disease recurrence (supplementary table 5). In this model, both TRPV2 gene expression and tumor grade remained independent predictors for disease recurrence. Altogether, these data suggest TRPV2 gene expression in the primary tumor biopsy as a valuable predictor for disease recurrence.
TRP channel expression in primary endometrial cancer cells
The gene expression data suggest that the expression of specific TRP channels correlates with the EMT status of the tumor in both primary biopsies and metastatic biopsies. However, tumor biopsies comprise two distinct but interdependent compartments: the neoplastic cancer cells (of epithelial origin) and the surrounding stroma, both contributing to the observed gene expression in our samples.
Hence, TRP channel expression was evaluated in primary endometrial epithelial cancer cells, isolated from tumor biopsies with different characteristics regarding tumor stage, histology and EMT status (Supplementary Fig. 5A-D). Validation and characterization of the cells is described elsewhere [17]. M/E ratios were determined for all four cell cultures (supplementary table 7 and Fig. 5I). The EM018a cells were derived from a stage IA tumor with mixed histopathology, comprising both endometrioid and serous components. Only the EM018a cells displayed obvious CDH1 expression, both at the protein (supplementary Fig. 5E) and RNA level (Fig. 8A) and displayed a low EMT status (M/E ratio < 1). The EM033 cells, derived from a stage IA tumor with dedifferentiated endometrioid histology, displayed a high EMT status but no expression of MMPs (Fig. 8A). In contrast, EM046 cells, isolated from a stage IB undifferentiated tumor, displayed high MMP2 expression, which suggests an active invasion phenotype in these cells. Lastly, the EM012 cells were isolated from a stage IV high grade endometrioid tumor. These cells display a moderate EMT status, with retained expression of KRT18 and MMP7, but high expression of mesenchymal markers, including MMP9 (Fig. 8A).
Interestingly, in accordance with the results in the tumor biopsies, in which TRPV2 expression was associated with increased tumor stage, TRPV2 expression was higher (10-fold increase) in stage IB (EM046 cells) and stage IV (EM012) cells, compared to the stage IA cell lines EM018a and EM033 (Fig. 8B). Similarly, TRPC1 expression was the highest in stage IB and stage IV cells, moderate in EM033 cells that underwent EMT and low in EM018 cells with a low EMT status. In contrast, TRPM4 expression was high in the EM018 cell line (stage IA with CDH1 expression intact), but undetectable in other cell lines. The TRPC4, TRPC6 and TRPV6 expression levels were low in all cells, while TRPM7 was highly expressed in all different cell cultures, with a 2-fold higher expression in EM046 cell types compared to the other cell lines (Fig. 8B). Overall, these results suggest that, like in tissue biopsies, the expression of TRPV2 and TRPC1 are associated with increased EMT status (or invasiveness), while TRPM4 expression is correlated with a low EMT status in endometrial cancer cells.