Next Article in Journal
Treatment of Recurrent Giant Cell Tumor of Bones: A Systematic Review
Next Article in Special Issue
Proteomic and Metabolomic Analysis of Bone Marrow and Plasma from Patients with Extramedullary Multiple Myeloma Identifies Distinct Protein and Metabolite Signatures
Previous Article in Journal
Atomic Force Microscopy Methods to Measure Tumor Mechanical Properties
Previous Article in Special Issue
Long-Chain Acyl Coenzyme A Dehydrogenase, a Key Player in Metabolic Rewiring/Invasiveness in Experimental Tumors and Human Mesothelioma Cell Lines
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Cancers
Volume 15
Issue 13
10.3390/cancers15133284
cancers-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Gene Expression Analysis Links Autocrine Vasoactive Intestinal Peptide and ZEB1 in Gastrointestinal Cancers

by
Ishani H. Rao
1,
Edmund K. Waller
1,
Rohan K. Dhamsania
2 and
Sanjay Chandrasekaran
3,*
1
Department of Hematology and Medical Oncology, Emory University School of Medicine, Atlanta, GA 30322, USA
2
Philadelphia College of Osteopathic Medicine (PCOM)-Georgia Campus, Suwanee, GA 30024, USA
3
Harold C. Simmons Cancer Center, University of Texas Southwestern Medical Center, Dallas, TX 75390, USA
*
Author to whom correspondence should be addressed.
Cancers 2023, 15(13), 3284; https://doi.org/10.3390/cancers15133284
Submission received: 30 May 2023 / Revised: 12 June 2023 / Accepted: 19 June 2023 / Published: 22 June 2023
(This article belongs to the Special Issue Biomarkers of Tumor Metastasis and Invasiveness)
Browse Figures
Versions Notes

Abstract

:

Simple Summary

The downstream signaling mechanisms and importance of the autocrine secretion of the vasoactive intestinal peptide (VIP) in cancer remains poorly understood. We hypothesized that VIP expression may promote cancer-associated signaling pathways. We analyzed gene sequencing data from cancer and healthy tissues based on the co-expression data of the VIP with 760 cancer-related genes. We identified a meaningful and novel association between the VIP and transcription factor ZEB1 in healthy and malignant human gastrointestinal tissues. ZEB1 is a known regulator of cancer EMT (epithelial–mesenchymal transition). Gene set analysis further supports the overlap in the EMT and cell cycle pathways. Our results identify a potentially novel function of the autocrine VIP as an important signaling peptide and biomarker of ZEB1-mediated EMT.

Abstract

VIP (vasoactive intestinal peptide) is a 28-amino acid peptide hormone expressed by cancer and the healthy nervous system, digestive tract, cardiovascular, and immune cell tissues. Many cancers express VIP and its surface receptors VPAC1 and VPAC2, but the role of autocrine VIP signaling in cancer as a targetable prognostic and predictive biomarker remains poorly understood. Therefore, we conducted an in silico gene expression analysis to study the mechanisms of autocrine VIP signaling in cancer. VIP expression from TCGA PANCAN tissue samples was analyzed against the expression levels of 760 cancer-associated genes. Of the 760 genes, 10 (MAPK3, ZEB1, TEK, NOS2, PTCH1 EIF4G1, GMPS, CDK2, RUVBL1, and TIMELESS) showed statistically meaningful associations with the VIP (Pearson’s R-coefficient > |0.3|; p < 0.05) across all cancer histologies. The strongest association with the VIP was for the epithelial–mesenchymal transition regulator ZEB1 in gastrointestinal malignancies. Similar positive correlations between the VIP and ZEB1 expression were also observed in healthy gastrointestinal tissues. Gene set analysis indicates the VIP is involved in the EMT and cell cycle pathways, and a high VIP and ZEB1 expression is associated with higher median estimate and stromal scores These findings uncover novel mechanisms for VIP- signaling in cancer and specifically suggest a role for VIP as a biomarker of ZEB1-mediated EMT. Further studies are warranted to characterize the specific mechanism of this interaction.

1. Introduction

Vasoactive Intestinal Peptide (VIP) is the final 28-amino acid (AA) product following multiple processing steps starting from prepro-VIP, a 170 AA precursor peptide [1,2] to which VIP primarily binds and signals through using two GPCRs (G-protein coupled receptors), VPAC1 and VPAC2 [3]. VIP, VPAC1, and VPAC2 are all broadly expressed in multiple healthy tissue types, in particular, tissues of the nervous system, digestive tract, and the respiratory, cardiovascular, and immune systems [4,5], and physiologically, VPAC/VIP signaling regulates embryonal development and growth, neuronal and epithelial cell signaling, gut absorption and motility, endocrine-mediated glycemia and circadian rhythms, immune cell responses, and carcinogenesis [6,7]. In the immune system, the VIP functions as a Type 2 cytokine [8]. VPAC1 is constitutively expressed in various immune cell populations, including T lymphocytes, macrophages, and dendritic cells [1,9]. Studies using transgenic mice have demonstrated that Th2 CD4+ T cells up-regulate VIP expression upon antigen presentation, while Th1 CD4+ T cells do not [10].
The importance of autocrine VIP as a biomarker or therapeutic target in cancer has remained elusive for many years. An in vitro study of >400 human primary tumors, tumor metastases, and normal tissues using subtype selective VPAC receptor autoradiography has established that most cancers express VPAC, with VPAC1 being the predominant receptor subtype in cancers and healthy tissues [11]. In clinical studies, radiolabeled-VIP binds VPAC and detects primary and metastatic tumors in a variety of cancer histologies [12,13], and serum VIP levels are > two-fold higher in patients with colon adenocarcinoma [14,15].
To date, the primary function of autocrine VIP signaling was thought to drive cancer cell proliferation. Mechanistically, VIP binds to surface VPAC1 and VPAC2 and generates downstream signals through adenylate cyclase to catalyze cyclic AMP (cAMP) synthesis and protein kinase A (PKA) activation [16]. The PKA phosphorylation of the CREB (cAMP-response element binding) transcription factor subsequently activates oncogenes such as c-Myc to drive cancer proliferation. In practice, in vitro studies with exogenous VIP or the pharmacologic inhibition of VIP signaling using VPAC1 and VPAC2 peptide antagonists have shown mixed results. No consistent increase in the proliferation rates of colon and pancreatic cancer cell lines was seen when adding exogenous VIP to cancer cells in vitro, and, in certain cell lines, exogenous VIP even inhibited growth [14]. Similarly, the pharmacologic inhibition of VIP signaling using VPAC1 and VPAC2 peptide antagonists in vitro demonstrated variable anti-tumor responses depending on cancer histology and the cell line tested, with positive responses primarily seen in pancreatic, colorectal, gastric, and breast cancer models [5,17,18,19,20,21]. Recent studies showed exogenous VIP feeds back to induce EGFR (epidermal growth factor receptor) and HER2 (human epidermal growth factor receptor 2) phosphorylation in breast cancer cells and inhibits the apoptotic effects of the RAS/RAF inhibitor sorafenib on cancer stem cells, suggesting there remain unknown downstream mechanisms of VIP [22,23,24,25].
The consistent co-expression of VIP and VPAC in cancer and the availability of radiotracers to follow the VIP/VPAC as predictive or prognostic biomarkers, coupled with the ability to inhibit VIP signaling via VPAC1 and VPAC2 peptide antagonists, supports the need to further characterize the VIP signaling landscape and determine its importance in the properties of cancer that have been described as the “hallmarks of cancer” [24,25]. In this exploratory analysis, we utilized an in silico gene expression model to uncover new downstream signaling pathways of VIP in cancer and corroborate these pathways in analyses of healthy tissues.

2. Materials and Methods

Gene association and expression: Cancer mRNA expression data for VIP and Pearson’s correlation (R) for VIP versus 760 other genes (gene panel based on the Nanostring nCounter® Tumor Signaling 360™ profiling panel) were abstracted from the TCGA (The Cancer Genome Atlas) PANCAN dataset using the University of California Santa Cruz (UCSC) Xena platform and recorded into Microsoft Excel by cancer hallmark (per Nanostring designation) (Supplemental Table S1) [26,27,28]. The date of cutoff for data abstraction was 1 December 2020 and of the 12,839 samples available for analysis, 2036 null values were excluded (Supplemental Figure S1.1). Individual tissue sample mRNA expression levels and R coefficients for VIP and 10 identified lead genes were further abstracted for all TCGA tumor histologies. Histologies were also organized into sub-groups by the tissue germ layer of origin (ectoderm, endoderm, and mesoderm) (Supplemental Table S2). Healthy tissue gene expression data was obtained from the Genotype-Tissue Expression (GTEx) Project on 19 April 2023 (https://www.gtexportal.org/home/ (19 April 2023)) [29]. Cancer cell line-specific gene expression data was downloaded on 4/19/23 from the CCLE (cancer cell line encyclopedia) using the DepMap resource (Public Expression Dataset 22Q4) (https://depmap.org/portal (accessed on 19 April 2023)).
Statistical analysis: The strength of association between VIP and comparator genes’ expression data from the TCGA, GTEX, and CCLE datasets was determined using a Pearson’s correlation test (R). For initial screening, R-coefficients > |0.3| and corresponding p-value < 0.05 were considered statistically meaningful [30]. Gene expression data (RNALog2(norm_count+1)) for VIP and ZEB1 in healthy (GTEX) and cancer (TCGA) tissues were analyzed per a Mann–Whitney U test for non-parametric data with p < 0.05 defined as statistically significant using GraphPad Prism (v9.5.1).
VIP and ZEB1 Expression and Cancer Pathway Activity: Gene set analysis was performed to assess the pathway activity of VIP and ZEB1 in gastrointestinal malignancies, including COAD (colon adenocarcinoma), ESCA (esophageal carcinoma), PAAD (pancreatic adenocarcinoma), and STAD (stomach adenocarcinoma). The analysis was conducted using the GSCA (Gene Set Cancer Analysis) tool, as previously described [31], with pathway activity scores assigned in ten critical cancer-related pathways (TSC/mTOR, RTK, RAS/MAPK, PI3K/AKT, hormone ER, hormone AR, EMT, DNA damage response, cell cycle, and apoptosis) based on a total of 7876 TCGA samples across 32 different cancer types. Pathway activity was determined by stratifying samples into high and low-expression groups based on the median gene expression, and a score was then calculated using a student t-test, with p-values adjusted by the false discovery rate (FDR). Reported results were considered significant when the FDR was <0.05, indicating a low likelihood of false positives.
Estimate, Stromal, and Immune Scoring: Estimate, stromal, and immune scores were used to compare tumor purity, stromal cell presence, and immune cell infiltration in stomach adenocarcinoma samples based on VIP and ZEB1 high vs. low tumor expression using the ESTIMATE scoring system, as previously described (https://bioinformatics.mdanderson.org/estimate/index.html (accessed on 8 June 2023)) [32]. A total of 377 of the 415 STAD samples with scoring available were evaluable (nulls removed and primary tumor only). The median estimate, stromal, and immune scores were analyzed per the Mann–Whitney U test for non-parametric data with p < 0.05 defined as statistically significant using GraphPad Prism (v9.5.1)

3. Results

3.1. VIP Expression in Cancer

The protein expression of VIP, VPAC1, and VPAC2 was investigated using immunohistochemistry data from the Human Protein Atlas (HPA) (Figure 1) [33,34]. Consistent with previous research [11], the HPA protein analysis revealed that VPAC1 is expressed in nearly all cancer histologies. In fact, VPAC1 expression was observed in over 70% of tissue samples from breast, carcinoid, head and neck, ovarian, and pancreatic cancers. On the other hand, the expression of VIP varied among different tumor histologies, with several histologies showing no detectable protein expression. However, it is noteworthy that many of the VIP non-expressing histologies still exhibited the expression of VPAC1 or VPAC2 receptors. For example, cancers of the breast, endometrium, prostate, testis, urothelium/bladder, and melanoma showed an expression of VPAC1 or VPAC2 despite lacking VIP expression. The differential expression patterns of VIP and its receptors across different histologies suggest distinct roles and potential interactions between VIP-receptor signaling and other pathways in cancer. Therefore, we undertook further investigations to elucidate the functional significance of these observations in the context of specific cancers.
To further characterize the signaling landscape, VIP mRNA expression was compared against a panel of 760 cancer-pertinent genes across all evaluable tissue samples in the TCGA PANCAN database. The selected gene panel was based on the Nanostring nCounter® Tumor Signaling 360TM testing panel [28]. A total of 285 genes showed positive R associations with VIP expression, while 475 were negatively associated (Figure 2A). Only ~10% of each cohort of positively or negatively associated genes demonstrated an R > |0.2| (25/285 positive correlation, 52/475 negative correlation). Based on our pre-specified cutoff of R > |0.3|, five positively correlated and five negatively correlated lead genes were identified (positive correlation: MAPK3, ZEB1, NOS2, TEK, and PTCH1; negative correlation: EIF4G1, GMPS, CDK2, RUVBL1, and TIMELESS) (all p < 0.05, Supplemental Figure S1.2).

3.2. mRNA Expression of VIP Associated Genes by Cancer Hallmark and Tissue Histology

To provide functional context to our findings, we also examined the association of the genes of interest with the pathways involved in the cancer hallmarks. Eight genes (CPA3, FAM30A, FOXP3, HDC, HSD11B1, PNIC, SH2D1A, and TCL1A) were excluded due to their primary involvement in the functions of non-malignant immune cells. Amongst the 752 genes remaining, the majority are involved in cancer hallmarks for activating sustaining proliferative signaling (n = 223), tumor-promoting inflammation (n = 178), invasion and metastasis (n = 146), and avoiding immune destruction (n = 122), with some overlap in certain genes across multiple cancer hallmarks (Supplemental Figure S1.3). Further focusing on lead gene functions, mapping by cancer hallmark, shows that VIP may be involved in the up- or down-regulation of genes relevant to nearly all cancer hallmarks (other than resisting cell death), with the higher number of lead gene associations seen in the hallmark pathway of “sustaining proliferative signaling” (RUVBL1, TEK, MAPK3, and PTCH1) (Figure 2B).
To explore the cancer histologies where the associations between VIP and lead genes are most significant, we employed increasingly stringent R-coefficient cutoffs (>|0.4| and >|0.5|). The R-coefficients between VIP and each lead gene were determined for all assessable samples within each TCGA cancer histology (Supplemental Figure S2). When using R > |0.4|, nine out of the ten lead genes exhibited histology-specific associations with VIP. Among them, TEK showed the highest number of histology-specific associations (n = 13), followed by ZEB1 (n = 8) (Figure 3). When further increasing the cutoff to R > |0.5|, only four out of the ten lead genes retained significant associations with VIP. These genes were ZEB1, TEK, GMPS, and TIMELESS. These associations showed histologic overlap in cancers of the gastrointestinal tract (STAD and COAD), lung (LUSC), and kidney (KICH). Notably, the strongest association observed was between VIP and ZEB1 in stomach adenocarcinoma (R = 0.76, Figure 4A). These findings highlight the cancer histologies where the association between the expression of VIP and lead genes is particularly meaningful, based on progressively stricter R-coefficient criteria. The strong association between VIP and ZEB1 in stomach adenocarcinoma warranted further investigation.
Given the known importance of VIP in embryogenesis, we also sub-grouped individual cancer histologies by germ layer of origin and looked for germ layer-specific expression associations (Supplemental Figure S3A, Supplemental Table S2). In this context, we also evaluated whether cross-reactive peptide signaling influenced our findings. The PACAP (pituitary adenylate cyclase-activating polypeptide) peptide shares 67% homology with VIP and is also able to bind to VPAC1 and VPAC2, in addition to its specific PAC1 receptor (Supplemental Figure S3B) [5,35]. Associations between lead genes and VIP or PACAP stratified by germinal layer did not indicate confounding functional overlap between VIP and PACAP. Germinal layer findings for VIP revealed tumors of ectodermal and mesodermal origin demonstrated correlations with R > |0.4| with non-meaningful R < |0.2| signals seen in PTCH1 (mesoderm) and GMPS, CDK2, RUVBL1, and TIMELESS (endoderm).

3.3. VIP and ZEB1 Expression in Healthy Tissue and Cancer Cell Lines

Given the compelling association between VIP and ZEB1 identified in gastrointestinal malignancies, especially of the stomach and colon, we tested VIP and ZEB1 expression in their corresponding healthy tissues. The healthy gastrointestinal tissue expression of VIP and ZEB1 was analyzed using the GTEX dataset [29]. Healthy gastric, colon, esophageal, and pancreatic tissues had R coefficients between VIP and ZEB1, an expression similar to those seen in TCGA cancer tissues (p < 0.0001 for all analyses) (Figure 4A). The biphasic distribution of VIP/ZEB1 (low) vs. VIP/ZEB1 (high) seen in the colon and esophageal tissues is thought to be related to the differences in the presence of mucosal and submucosal tissue in the samples. VIP expression in healthy pancreatic tissue is likely more susceptible to hormonal regulation of the endocrine pancreatic tissue, supporting the findings of variable VIP expression with somewhat consistent ZEB1 expression.
To determine if VIP and ZEB1 expression is up-regulated or down-regulated in cancer vs. healthy tissue, the total mRNA levels of VIP and ZEB1 between healthy and cancer tissue were compared. VIP and ZEB1 expression were both lower in stomach, colon, and esophageal cancer tissue compared to healthy tissue (p < 0.0001, Figure 4B). There were minimally significant differences in VIP expression between normal and pancreatic cancer tissue (p = 0.012) and increased ZEB1 expression in cancer tissue. We next analyzed associations between VIP and ZEB1 expression in the human stomach, colon, esophageal, and pancreatic adenocarcinoma cell lines from the CCLE (cancer cell line encyclopedia) using the DepMap resource. Interestingly, the cell line expression data did not recapitulate the VIP and ZEB1 associations seen in healthy and cancer tissues and showed only a statistically significant association between VIP and ZEB1 in a limited number of esophageal adenocarcinoma cases (p = 0.05) (Figure 4C).

3.4. Gene Set Analysis of VIP and ZEB1

The consistent association between VIP and ZEB1 expression in both healthy and cancerous tissues suggests a potential linkage of the signaling pathways of these two genes. To investigate this further, gene set analysis (GSA) was conducted in colon, gastric, esophageal, and pancreatic cancers using the GSCA tool for the following cancer-related pathways: TSC/mTOR, RTK, RAS/MAPK, PI3K/AKT, hormone ER, hormone AR, EMT, DNA damage response, cell cycle, and apoptosis. GSA revealed that VIP and ZEB1 had activating effects on EMT and inhibitory effects on the cell cycle pathway in both colon and gastric cancers. Additionally, in gastric cancer, VIP and ZEB1 showed inhibitory effects on the apoptosis pathway (FDR < 0.05, Figure 5A, Supplemental Table S3). Comparisons between the high and low tissue expression of VIP and ZEB1 were performed based on the GSA results, focusing on the EMT and cell cycle pathways in gastric and colon cancer (Figure 5B). These findings indicate that VIP and ZEB1 may jointly influence the EMT and cell cycle pathways in gastric and colon cancers. Moreover, the inhibitory effects on the apoptosis pathway observed in gastric cancer suggest potential roles for VIP and ZEB1 in modulating cell survival mechanisms in this specific cancer type.

3.5. Estimate, Stromal, and Immune Scoring

VIP is also known to exert paracrine effects in cancer by supporting the immunosuppressive activity of CD4+CD25+ regulatory T cells and tolerogenic dendritic cells [8,36]. In vivo studies have demonstrated that treatment with VPAC peptide antagonists can reverse these immune effects and down-regulate T cell expression of the inhibitory marker PD-1, promote cytotoxic T cell differentiation and expansion, and enhance intratumoral T cell infiltration, ultimately leading to tumor elimination [15,37,38,39], and we have recently shown that inhibiting VIP in murine pancreatic cancer models improves the response to T cell checkpoint inhibitors [15]. To determine if the VIP and ZEB1 associations seen in our analysis were confounded by immune cell-mediated paracrine VIP, we compared whether high versus low VIP and ZEB1 expression correlated with tumor purity, stromal tissue presence, or immune cell infiltration using the estimate, stromal, and immune scores, respectively. ZEB1 and VIP high-expressing tumors (n = 189) were associated with higher median estimate and stromal scores, with high expression most associated with positive scores (Figure 6A,B). Median immune scores were positive for low and high ZEB1 and VIP expression, with higher scores associated with increased VIP and ZEB1 expression (Figure 6C).

4. Discussion

To comprehensively characterize the signaling landscape of autocrine VIP in cancer, we conducted an in silico gene expression analysis. Determining the downstream effects of VIP in cancer is crucial to understanding its role as a circulating biomarker and VPAC as a diagnostic, prognostic, or predictive target. The utility of radiolabeled-VIP and the development of VIP-targeted theranostics is severely limited by this lack of understanding.
Comprehensive gene expression studies of cancer are valuable when studying peptides such as VIP. The VIP sequence is highly conserved across mammals and rarely mutated in cancers, in support of a critical role for VIP-signaling in normal physiology [2]. Furthermore, protein expression studies, including techniques such as Western blot and immunohistochemistry (IHC), may not fully capture peptide expression data due to limiting factors such as the multiple cleavage steps of the prepro peptide, variability in the antibody binding, and the short half-life of the peptide.
Among the initial screening panel of 760 genes, our in silico model identified 10 lead genes that showed strong associations with VIP expression. Some of these genes, such as MAPK, NOS2, CDK2, and TIMELESS, have been previously associated with VIP in non-cancer models [40,41,42,43]. However, our analysis also revealed several novel associations, the most compelling being the association with the ZEB1 transcription factor.
The association between VIP and ZEB1 in gastrointestinal tissues appears to be present in both healthy and malignant tissues. ZEB1 is a well-characterized transcription factor known to play a role in the regulation of epithelial–mesenchymal transition (EMT) in cancer [44]. EMT is a process in which primary tumor cells down-regulate the expression of structural and cell adhesion molecules, such as cytokeratin and E-cadherin, and lose their epithelial markers to acquire a more aggressive, spindle cell-like mesenchymal phenotype to drive tumor metastasis and distant invasion [45].
Detecting the early signs of EMT and targeting the process pharmacologically has been a long-standing focus of cancer research, but the heterogeneity across multiple cancer histologies and redundant downstream signaling pathways involved in E-cadherin loss have limited the success of these efforts. Transcription factors critical to EMT, including ZEB1, STAT, SNAIL, and TWIST, are not ideal targets for pharmacological interventions or readily assessed through imaging techniques. Therefore, identifying upstream cell surface targets such as VPAC1 or VPAC2 could aid in the development of EMT-targeted cancer therapies [45]. Recently, Colangelo et al. uncovered an axis between TIMELESS and ZEB1 in colorectal cancers, demonstrating that the loss of TIMELESS promotes tumor progression and poor prognosis by inducing ZEB1 expression and EMT. These findings align with the negative correlation we observed between VIP and TIMELESS in our analysis [46].
However, the current study has some design limitations. The 760-gene panel used is a limited dataset, however, this selection was deliberate to focus on well-elucidated genes with known functions in cancer to facilitate the further translation of our findings. The under-representation of certain histologies in the TCGA dataset may have led to histology-specific errors and potential missed associations due to a smaller sample size. The initial screening of lead genes by R > {0.3} may appear to be a less stringent threshold, however, this was done only when evaluating the VIP-gene associations across all TCGA histologies and samples. Of note, all R in the initial screen were < {0.4}.
Findings from in silico analyses can be challenging to translate into in vitro studies or in vivo therapies. The divergent findings in VIP and ZEB1 expression associations between tissues (healthy and malignant) and cancer cell lines highlight the need for future studies in both settings. Factors such as culture conditions may affect VIP, ZEB1, and EMT-related phenotypes in vitro. Heterogeneity in the tumor microenvironment and between sequenced samples may also affect data interpretation, however, the estimate, stromal, and immune scoring data support that our findings are not confounded by the paracrine effects of immune cell infiltrates and are cancer cell and cancer stroma dependent. The differences seen by the cancer-mediated stromal cell presence may also indicate why the patterns vary between healthy tissues and cancer tissues.
The identified association between ZEB1 and VIP expression may not be causal. However, the consistent association of the expression levels for these genes across multiple cancer histologies, healthy tissues, and gene set analyses indicates potential significance and warrants further investigation. Furthermore, while the transition from a proliferative to invasive phenotype is considered fundamental to EMT, our findings that VIP is involved in both processes support more recent arguments that EMT and, conversely, MET (mesenchymal–epithelial transition) represent a continuum of states and not discrete phenotypes [47].
Taken together, the VIP–VPAC axis offers a widely expressed and targetable biomarker in cancer. The presence of VPAC in both VIP-expressing and non-expressing tissues suggests that cancers may be simultaneously susceptible to both autocrine and paracrine VIP effects. Whether targeting VPAC can be diagnostic, prognostic, and/or predictive remains to be fully explored. Our findings indicate that in gastrointestinal malignancies, further in vitro and in vivo validation is necessary to elucidate the mechanism of interaction between VIP and ZEB1 and to determine how VIP supports ZEB1-mediated EMT involved in tumor invasion and metastasis.

5. Conclusions

We performed hypothesis-generating in silico analyses of associations of autocrine VIP gene expression with gene targets and pathways in cancer. We identified a potential novel association of the VIP–VPAC axis with ZEB1 expression in gastrointestinal malignancies. Further studies are necessary to understand how inhibiting or activating the VIP–VPAC axis may serve as a prognostic biomarker of cancer metastasis or a predictive biomarker for response to VIP-targeted therapies on EMT in gastrointestinal malignancies.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/cancers15133284/s1, Supplemental Figure S1: Analyzed sample numbers, p-values, and pertinent cancer hallmarks for the 760 target gene panel; Supplemental Figure S2: Heat map of R-values for lead gene associations by cancer histology; Supplemental Figure S3: Heat map of lead gene associations by tissue germ layer; Supplemental Table S1: R-Values in target gene set by cancer hallmark with heat map; Supplemental Table S2. TCGA tumor histologies and abbreviations by germ layer; and Supplemental Table S3: Gene set analysis data for VIP and ZEB1.

Author Contributions

Conceptualization, S.C. and I.H.R.; methodology, S.C., I.H.R. and R.K.D.; software, S.C.; formal analysis, S.C., I.H.R. and E.K.W.; investigation, I.H.R.; resources, S.C. and E.K.W.; data curation, S.C.; writing—original draft preparation, I.H.R. and R.K.D.; writing—review and editing, S.C., I.H.R. and E.K.W.; visualization, S.C.; supervision, S.C.; project administration, I.H.R.; funding acquisition, S.C. and E.K.W. All authors have read and agreed to the published version of the manuscript.

Funding

These studies were supported by funding from the Winship Cancer Institute Elkin Fellowship (SC), NIH T32 CA160040 (SC), and Abraham J. and Phyllis Katz Foundation (EKW).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

TCGA data can be accessed via UCSC Xena at (http://xena.ucsc.edu/ (accessed on 1 December 2020)). Healthy tissue data are available at the GTEX portal (https://www.gtexportal.org/home/ (accessed on 19 April 2023)). Cancer cell line data are available through the cancer cell line encyclopedia (CCCLE) available via the DepMap resource (https://depmap.org/portal (accessed on 19 April 2023)).

Acknowledgments

The results are based upon data generated from the TCGA Research Network: https://www.cancer.gov/tcga (accessed on 1 December 2020)). The TCGA data used for the analyses were obtained from the UCSC Xena Portal (https://xena.ucsc.edu (accessed on 1 December 2020)). The Genotype-Tissue Expression (GTEx) Project was supported by the Common Fund of the Office of the Director of the National Institutes of Health, and by NCI, NHGRI, NHLBI, NIDA, NIMH, and NINDS. The healthy tissue data used for the analyses described in this manuscript were obtained from the GTEx Portal on 19 April 2023.

Conflicts of Interest

E.K.W. holds intellectual property covering VIP antagonists and is a founder and equity partner in Cambium Oncology. The other authors have no financial conflict of interest.

References

  1. Delgado, M.; Pozo, D.; Ganea, D. The Significance of Vasoactive Intestinal Peptide in Immunomodulation. Pharmacol. Rev. 2004, 56, 249–290. [Google Scholar] [CrossRef]
  2. Ng, S.Y.L.; Chow, B.K.C.; Kasamatsu, J.; Kasahara, M.; Lee, L.T.O. Agnathan VIP, PACAP and Their Receptors: Ancestral Origins of Today’s Highly Diversified Forms. PLoS ONE 2012, 7, e44691. [Google Scholar] [CrossRef] [Green Version]
  3. Harmar, A.J.; Arimura, A.; Gozes, I.; Journot, L.; Laburthe, M.; Pisegna, J.R.; Rawlings, S.R.; Robberecht, P.; Said, S.I.; Sreedharan, S.P.; et al. International Union of Pharmacology. XVIII. Nomenclature of receptors for vasoactive intestinal peptide and pituitary adenylate cyclase-activating polypeptide. Pharmacol. Rev. 1998, 50, 265–270. [Google Scholar]
  4. Dickson, L.; Finlayson, K. VPAC and PAC receptors: From ligands to function. Pharmacol. Ther. 2009, 121, 294–316. [Google Scholar] [CrossRef]
  5. Moody, T.W.; Nuche-Berenguer, B.; Jensen, R.T. VIP/PACAP, and their receptors and cancer. Curr. Opin. Endocrinol. Diabetes Obes. 2016, 23, 38–47. [Google Scholar] [CrossRef] [Green Version]
  6. Iwasaki, M.; Akiba, Y.; Kaunitz, J.D. Recent advances in vasoactive intestinal peptide physiology and pathophysiology: Focus on the gastrointestinal system. F1000Research 2019, 8, 1629. [Google Scholar] [CrossRef] [Green Version]
  7. Fabricius, D.; Karacay, B.; Shutt, D.; Leverich, W.; Schafer, B.; Takle, E.; Thedens, D.; Khanna, G.; Raikwar, S.; Yang, B.; et al. Characterization of Intestinal and Pancreatic Dysfunction in VPAC1-Null Mutant Mouse. Pancreas 2011, 40, 861–871. [Google Scholar] [CrossRef]
  8. Delgado, M.; Ganea, D. Vasoactive intestinal peptide: A neuropeptide with pleiotropic immune functions. Amino Acids 2013, 45, 25–39. [Google Scholar] [CrossRef]
  9. Smalley, S.G.R.; Barrow, P.A.; Foster, N. Immunomodulation of innate immune responses by vasoactive intestinal peptide (VIP): Its therapeutic potential in inflammatory disease. Clin. Exp. Immunol. 2009, 157, 225–234. [Google Scholar] [CrossRef]
  10. Vassiliou, E.; Jiang, X.; Delgado, M.; Ganea, D. TH2 Lymphocytes Secrete Functional VIP upon Antigen Stimulation. Arch. Physiol. Biochem. 2001, 109, 365–368. [Google Scholar] [CrossRef]
  11. Reubi, J.C.; Läderach, U.; Waser, B.; Gebbers, J.O.; Robberecht, P.; Laissue, J.A. Vasoactive intestinal peptide/pituitary adenylate cyclase-activating peptide receptor subtypes in human tumors and their tissues of origin. Cancer Res 2000, 60, 3105–3112. [Google Scholar] [PubMed]
  12. Virgolini, I.; Raderer, M.; Kurtaran, A.; Angelberger, P.; Yang, Q.; Radosavljevic, M.; Leimer, M.; Kaserer, K.; Li, S.; Kornek, G.; et al. 123I-vasoactive intestinal peptide (VIP) receptor scanning: Update of imaging results in patients with adenocarcinomas and endocrine tumors of the gastrointestinal tract. Nucl. Med. Biol. 1996, 23, 685–692. [Google Scholar] [CrossRef] [PubMed]
  13. Virgolini, I. Mack Forster Award Lecture Receptor nuclear medicine: Vasointestinal peptide and somatostatin receptor scintigraphy for diagnosis and treatment of tumour patients. Eur. J. Clin. Investig. 1997, 27, 793–800. [Google Scholar] [CrossRef] [PubMed]
  14. Hejna, M.; Hamilton, G.; Brodowicz, T.; Haberl, I.; Fiebiger, W.C.; Scheithauer, W.; Virgolini, I.; Köstler, W.J.; Oberhuber, G.; Raderer, M. Serum levels of vasoactive intestinal peptide (VIP) in patients with adenocarcinomas of the gastrointestinal tract. Anticancer. Res. 2001, 21, 1183–1187. [Google Scholar]
  15. Ravindranathan, S.; Passang, T.; Li, J.-M.; Wang, S.; Dhamsania, R.; Ware, M.B.; Zaidi, M.Y.; Zhu, J.; Cardenas, M.; Liu, Y.; et al. Targeting vasoactive intestinal peptide-mediated signaling enhances response to immune checkpoint therapy in pancreatic ductal adenocarcinoma. Nat. Commun. 2022, 13, 6418. [Google Scholar] [CrossRef] [PubMed]
  16. Harmar, A.J.; Fahrenkrug, J.; Gozes, I.; Laburthe, M.; May, V.; Pisegna, J.R.; Vaudry, D.; Vaudry, H.; Waschek, J.A.; Said, S.I. Pharmacology and functions of receptors for vasoactive intestinal peptide and pituitary adenylate cyclase-activating polypeptide: IUPHAR Review 1. Br. J. Pharmacol. 2012, 166, 4–17. [Google Scholar] [CrossRef] [Green Version]
  17. Laburthe, M.; Rousset, M.; Chevalier, G.; Boissard, C.; Dupont, C.; Zweibaum, A.; Rosselin, G. Vasoactive intestinal peptide control of cyclic adenosine 3′:5′-monophosphate levels in seven human colorectal adenocarcinoma cell lines in culture. Cancer Res. 1980, 40, 2529–2533. [Google Scholar]
  18. Laburthe, M.; Rousset, M.; Boissard, C.; Chevalier, G.; Zweibaum, A.; Rosselin, G. Vasoactive intestinal peptide: A potent stimulator of adenosine 3′:5′-cyclic monophosphate accumulation in gut carcinoma cell lines in culture. Proc. Natl. Acad. Sci. USA 1978, 75, 2772–2775. [Google Scholar] [CrossRef] [Green Version]
  19. Moody, T.W.; Hill, J.M.; Jensen, R.T. VIP as a trophic factor in the CNS and cancer cells. Peptides 2003, 24, 163–177. [Google Scholar] [CrossRef]
  20. Moody, T.W.; Zia, F.; Draoui, M.; Brenneman, D.E.; Fridkin, M.; Davidson, A.; Gozes, I. A vasoactive intestinal peptide antagonist inhibits non-small cell lung cancer growth. Proc. Natl. Acad. Sci. USA 1993, 90, 4345–4349. [Google Scholar] [CrossRef] [Green Version]
  21. Moody, T.W.; Leyton, J.; Chan, D.; Brenneman, D.C.; Fridkin, M.; Gelber, E.; Levy, A.; Gozes, I. VIP receptor antagonists and chemotherapeutic drugs inhibit the growth of breast cancer cells. Breast Cancer Res. Treat. 2001, 68, 55–64. [Google Scholar] [CrossRef] [PubMed]
  22. Valdehita, A.; Bajo, A.M.; Schally, A.V.; Varga, J.L.; Carmena, M.J.; Prieto, J.C. Vasoactive intestinal peptide (VIP) induces transactivation of EGFR and HER2 in human breast cancer cells. Mol. Cell. Endocrinol. 2009, 302, 41–48. [Google Scholar] [CrossRef] [PubMed]
  23. Sastry, K.S.; Chouchane, A.I.; Wang, E.; Kulik, G.; Marincola, F.M.; Chouchane, L. Cytoprotective effect of neuropeptides on cancer stem cells: Vasoactive intestinal peptide-induced antiapoptotic signaling. Cell Death Dis. 2017, 8, e2844. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  24. Hanahan, D.; Weinberg, R.A. Hallmarks of cancer: The next generation. Cell 2011, 144, 646–674. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  25. Hanahan, D. Hallmarks of Cancer: New Dimensions. Cancer Discov. 2022, 12, 31–46. [Google Scholar] [CrossRef]
  26. The Cancer Genome Atlas n.d. Available online: https://www.cancer.gov/about-nci/organization/ccg/research/structural-genomics/tcga (accessed on 1 January 2020).
  27. Goldman, M.J.; Craft, B.; Hastie, M.; Repečka, K.; McDade, F.; Kamath, A.; Banerjee, A.; Luo, Y.; Rogers, D.; Brooks, A.N.; et al. Visualizing and interpreting cancer genomics data via the Xena platform. Nat. Biotechnol. 2020, 38, 675–678. [Google Scholar] [CrossRef]
  28. Nanostring Tumor Signaling 360 n.d. Available online: https://nanostring.com/support-documents/tumor-signaling-360-panel/ (accessed on 1 January 2020).
  29. Carithers, L.J.; Moore, H.M. The Genotype-Tissue Expression (GTEx) Project. Biopreservation Biobanking 2015, 13, 307–308. [Google Scholar] [CrossRef]
  30. Schober, P.; Boer, C.; Schwarte, L.A. Correlation Coefficients: Appropriate Use and Interpretation. Anesth. Analg. 2018, 126, 1763–1768. [Google Scholar] [CrossRef]
  31. Liu, C.-J.; Hu, F.-F.; Xie, G.-Y.; Miao, Y.-R.; Li, X.-W.; Zeng, Y.; Guo, A.-Y. GSCA: An integrated platform for gene set cancer analysis at genomic, pharmacogenomic and immunogenomic levels. Brief. Bioinform. 2023, 24, bbac558. [Google Scholar] [CrossRef]
  32. Yoshihara, K.; Shahmoradgoli, M.; Martínez, E.; Vegesna, R.; Kim, H.; Torres-Garcia, W.; Treviño, V.; Shen, H.; Laird, P.W.; Levine, D.A.; et al. Inferring tumour purity and stromal and immune cell admixture from expression data. Nat. Commun. 2013, 4, 2612. [Google Scholar] [CrossRef] [Green Version]
  33. The Human Protein Atlas, n.d. Available online: https://www.proteinatlas.org/ (accessed on 1 January 2020).
  34. Uhlén, M.; Fagerberg, L.; Hallström, B.M.; Lindskog, C.; Oksvold, P.; Mardinoglu, A.; Sivertsson, Å.; Kampf, C.; Sjöstedt, E.; Asplund, A.; et al. Proteomics. Tissue-Based Map of the Human Proteome. Science 2015, 347, 1260419. [Google Scholar] [CrossRef] [PubMed]
  35. Vaudry, D.; Falluel-Morel, A.; Bourgault, S.; Basille, M.; Burel, D.; Wurtz, O.; Fournier, A.; Chow, B.K.C.; Hashimoto, H.; Galas, L.; et al. Pituitary Adenylate Cyclase-Activating Polypeptide and Its Receptors: 20 Years after the Discovery. Pharmacol. Rev. 2009, 61, 283–357. [Google Scholar] [CrossRef] [PubMed]
  36. Gonzalez-Rey, E.; Chorny, A.; Fernandez-Martin, A.; Ganea, D.; Delgado, M. Vasoactive intestinal peptide generates human tolerogenic dendritic cells that induce CD4 and CD8 regulatory T cells. Blood 2006, 107, 3632–3638. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  37. Li, J.-M.; Petersen, C.T.; Li, J.-X.; Panjwani, R.; Chandra, D.J.; Giver, C.R.; Blazar, B.R.; Waller, E.K. Modulation of Immune Checkpoints and Graft-versus-Leukemia in Allogeneic Transplants by Antagonizing Vasoactive Intestinal Peptide Signaling. Cancer Res. 2016, 76, 6802–6815. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  38. Petersen, C.T.; Hassan, M.; Morris, A.B.; Jeffery, J.; Lee, K.; Jagirdar, N.; Staton, A.D.; Raikar, S.S.; Spencer, H.T.; Sulchek, T.; et al. Improving ex vivo T cell expansion from DLBCL patients for T cell therapies via antagonism of PI3K δ and VIP. Blood 2017, 130, 3195. [Google Scholar]
  39. Petersen, C.T.; Hassan, M.; Morris, A.B.; Jeffery, J.; Lee, K.; Jagirdar, N.; Staton, A.D.; Raikar, S.S.; Spencer, H.T.; Sulchek, T.; et al. Improving T-cell expansion and function for adoptive T-cell therapy using ex vivo treatment with PI3Kδ inhibitors and VIP antagonists. Blood Adv. 2018, 2, 210–223. [Google Scholar] [CrossRef] [Green Version]
  40. Fernández, M.; Sánchez-Franco, F.; Palacios, N.; Sánchez, I.; Cacicedo, L. IGF-I and vasoactive intestinal peptide (VIP) regulate cAMP-response element-binding protein (CREB)-dependent transcription via the mitogen-activated protein kinase (MAPK) pathway in pituitary cells: Requirement of Rap1. J. Mol. Endocrinol. 2005, 34, 699–712. [Google Scholar] [CrossRef]
  41. MacEachern, S.J.; Patel, B.A.; Keenan, C.M.; Dicay, M.; Chapman, K.; McCafferty, D.-M.; Savidge, T.C.; Beck, P.L.; MacNaughton, W.K.; Sharkey, K.A. Inhibiting Inducible Nitric Oxide Synthase in Enteric Glia Restores Electrogenic Ion Transport in Mice With Colitis. Gastroenterology 2015, 149, 445–455.e3. [Google Scholar] [CrossRef] [Green Version]
  42. Anderson, P.; Gonzalez-Rey, E. Vasoactive Intestinal Peptide Induces Cell Cycle Arrest and Regulatory Functions in Human T Cells at Multiple Levels. Mol. Cell. Biol. 2010, 30, 2537–2551. [Google Scholar] [CrossRef] [Green Version]
  43. Muraro, N.; Pírez, N.; Ceriani, M. The circadian system: Plasticity at many levels. Neuroscience 2013, 247, 280–293. [Google Scholar] [CrossRef]
  44. Caramel, J.; Ligier, M.; Puisieux, A. Pleiotropic Roles for ZEB1 in Cancer. Cancer Res. 2018, 78, 30–35. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  45. Davis, F.M.; Stewart, T.A.; Thompson, E.W.; Monteith, G.R. Targeting EMT in cancer: Opportunities for pharmacological intervention. Trends Pharmacol. Sci. 2014, 35, 479–488. [Google Scholar] [CrossRef] [PubMed]
  46. Colangelo, T.; Carbone, A.; Mazzarelli, F.; Cuttano, R.; Dama, E.; Nittoli, T.; Albanesi, J.; Barisciano, G.; Forte, N.; Palumbo, O.; et al. Loss of circadian gene Timeless induces EMT and tumor progression in colorectal cancer via Zeb1-dependent mechanism. Cell Death Differ. 2022, 29, 1552–1568. [Google Scholar] [CrossRef] [PubMed]
  47. Yang, J.; Antin, P.; Berx, G.; Blanpain, C.; Brabletz, T.; Bronner, M.; Campbell, K.; Cano, A.; Casanova, J.; Christofori, G.; et al. Guidelines and definitions for research on epithelial–mesenchymal transition. Nat. Rev. Mol. Cell Biol. 2020, 21, 341–352, Correction in Nat. Rev. Mol. Cell Biol. 2021, 22, 834. [Google Scholar] [CrossRef] [PubMed] [Green Version]
Cancers 15 03284 g001 550
Figure 1. Protein expression of VIP and VIP Receptors In Cancer: VIP expression is highest in gliomas and colorectal, pancreas, thyroid cancers. VPAC1 is generally expressed on all cancer types and VPAC2 expression is limited to certain cancer types including those of the head and neck, lung, skin, kidney, thyroid, bladder, and melanoma.
Figure 1. Protein expression of VIP and VIP Receptors In Cancer: VIP expression is highest in gliomas and colorectal, pancreas, thyroid cancers. VPAC1 is generally expressed on all cancer types and VPAC2 expression is limited to certain cancer types including those of the head and neck, lung, skin, kidney, thyroid, bladder, and melanoma.
Cancers 15 03284 g001
Cancers 15 03284 g002 550
Figure 2. 10 Lead Genes Identified and Associated Cancer Hallmarks: (A) 10 lead-genes identified with R values > |0.3| for the entire analyzed dataset. (B) Further analysis by cancer hallmark shows potential associations of VIP with hallmarks other than proliferative signaling.
Figure 2. 10 Lead Genes Identified and Associated Cancer Hallmarks: (A) 10 lead-genes identified with R values > |0.3| for the entire analyzed dataset. (B) Further analysis by cancer hallmark shows potential associations of VIP with hallmarks other than proliferative signaling.
Cancers 15 03284 g002
Cancers 15 03284 g003 550
Figure 3. Lead Gene associations by tumor histology and increasing R-value Strength: Lead genes with corresponding histologies by increasing R-value strength. At R > |0.5|, only 4/10 lead genes retain meaningful associations with VIP.
Figure 3. Lead Gene associations by tumor histology and increasing R-value Strength: Lead genes with corresponding histologies by increasing R-value strength. At R > |0.5|, only 4/10 lead genes retain meaningful associations with VIP.
Cancers 15 03284 g003
Cancers 15 03284 g004 550
Figure 4. VIP and ZEB1 mRNA expression in gastrointestinal cancer and health tissues: (A) Similar VIP vs ZEB1 associations are seen in gastric, colon, esophageal, and pancreatic cancer (TCGA) and healthy (GTEX) tissues (p < 0.0001 for all analyses). (B) Total mRNA expression of VIP and ZEB1 is reduced in cancers of the stomach, colon, and esophagus (p < 0.0001 (****)), increased ZEB1 in cancer tissue (p < 0.0001), and similar in the pancreas for VIP p = 0.012 (*). (C) VIP and ZEB1 mRNA expression (log2(TPM+1)) in human gastric (stomach), colon, and pancreatic adenocarcinoma cancer cell lines (CCLE) does not show meaningful expression correlations. Positive correlation is seen in esophageal cell lines (R = 0.8, p = 0.05).
Figure 4. VIP and ZEB1 mRNA expression in gastrointestinal cancer and health tissues: (A) Similar VIP vs ZEB1 associations are seen in gastric, colon, esophageal, and pancreatic cancer (TCGA) and healthy (GTEX) tissues (p < 0.0001 for all analyses). (B) Total mRNA expression of VIP and ZEB1 is reduced in cancers of the stomach, colon, and esophagus (p < 0.0001 (****)), increased ZEB1 in cancer tissue (p < 0.0001), and similar in the pancreas for VIP p = 0.012 (*). (C) VIP and ZEB1 mRNA expression (log2(TPM+1)) in human gastric (stomach), colon, and pancreatic adenocarcinoma cancer cell lines (CCLE) does not show meaningful expression correlations. Positive correlation is seen in esophageal cell lines (R = 0.8, p = 0.05).
Cancers 15 03284 g004
Cancers 15 03284 g005 550
Figure 5. Gene Set Analysis of VIP and ZEB1 in colon, esophageal, gastric, and pancreatic cancers: (A) VIP and ZEB1 share activating effects on EMT and inhibitory effects on cell cycle in gastric and colon cancer. No pathway effects seen in pancreatic cancer (B) Representative box plots showing impact of high vs low VIP and ZEB1 expression in gastric and colon cancer in EMT activation and cell cycle inhibition (FDR < 0.25).
Figure 5. Gene Set Analysis of VIP and ZEB1 in colon, esophageal, gastric, and pancreatic cancers: (A) VIP and ZEB1 share activating effects on EMT and inhibitory effects on cell cycle in gastric and colon cancer. No pathway effects seen in pancreatic cancer (B) Representative box plots showing impact of high vs low VIP and ZEB1 expression in gastric and colon cancer in EMT activation and cell cycle inhibition (FDR < 0.25).
Cancers 15 03284 g005
Cancers 15 03284 g006 550
Figure 6. Estimate, Stromal, and Immune scores grouped by low versus high VIP and ZEB1 Expression: (A) Median Estimate score for ZEB1 high (1569) and ZEB1 low (−348.7) and VIP high (1428) and VIP low (−66.05) (B). Median Stromal score for ZEB1 high (568.2) and ZEB1 low (−621.5) and VIP high (387.2) and VIP low (−397.6) (C). Median Immune score for ZEB1 high (864) and ZEB1 low (215) and VIP high (814.2) and VIP low (264.1). p < 0.0001 for all comparisons (****).
Figure 6. Estimate, Stromal, and Immune scores grouped by low versus high VIP and ZEB1 Expression: (A) Median Estimate score for ZEB1 high (1569) and ZEB1 low (−348.7) and VIP high (1428) and VIP low (−66.05) (B). Median Stromal score for ZEB1 high (568.2) and ZEB1 low (−621.5) and VIP high (387.2) and VIP low (−397.6) (C). Median Immune score for ZEB1 high (864) and ZEB1 low (215) and VIP high (814.2) and VIP low (264.1). p < 0.0001 for all comparisons (****).
Cancers 15 03284 g006
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Rao, I.H.; Waller, E.K.; Dhamsania, R.K.; Chandrasekaran, S. Gene Expression Analysis Links Autocrine Vasoactive Intestinal Peptide and ZEB1 in Gastrointestinal Cancers. Cancers 2023, 15, 3284. https://doi.org/10.3390/cancers15133284

AMA Style

Rao IH, Waller EK, Dhamsania RK, Chandrasekaran S. Gene Expression Analysis Links Autocrine Vasoactive Intestinal Peptide and ZEB1 in Gastrointestinal Cancers. Cancers. 2023; 15(13):3284. https://doi.org/10.3390/cancers15133284

Chicago/Turabian Style

Rao, Ishani H., Edmund K. Waller, Rohan K. Dhamsania, and Sanjay Chandrasekaran. 2023. "Gene Expression Analysis Links Autocrine Vasoactive Intestinal Peptide and ZEB1 in Gastrointestinal Cancers" Cancers 15, no. 13: 3284. https://doi.org/10.3390/cancers15133284

APA Style

Rao, I. H., Waller, E. K., Dhamsania, R. K., & Chandrasekaran, S. (2023). Gene Expression Analysis Links Autocrine Vasoactive Intestinal Peptide and ZEB1 in Gastrointestinal Cancers. Cancers, 15(13), 3284. https://doi.org/10.3390/cancers15133284

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop