Table of Contents Author Guidelines Submit a Manuscript
International Journal of Breast Cancer
Volume 2015 (2015), Article ID 539721, 10 pages
http://dx.doi.org/10.1155/2015/539721
Research Article

SEMA6D Expression and Patient Survival in Breast Invasive Carcinoma

1Division of Preventive Medicine, Department of Medicine, University of Alabama at Birmingham, Birmingham, AL 35294, USA
2Comprehensive Cancer Center, Department of Medicine, University of Alabama at Birmingham, Birmingham, AL 35294, USA
3Research Division, Department of Genetics, University of Alabama at Birmingham, Birmingham, AL 35294, USA

Received 12 February 2015; Accepted 26 March 2015

Academic Editor: Ian S. Fentiman

Copyright © 2015 Dongquan Chen et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract

Breast cancer (BC) is the second most common cancer diagnosed in American women and is also the second leading cause of cancer death in women. Research has focused heavily on BC metastasis. Multiple signaling pathways have been implicated in regulating BC metastasis. Our knowledge of regulation of BC metastasis is, however, far from complete. Identification of new factors during metastasis is an essential step towards future therapy. Our labs have focused on Semaphorin 6D (SEMA6D), which was implicated in immune responses, heart development, and neurogenesis. It will be interesting to know SEMA6D-related genomic expression profile and its implications in clinical outcome. In this study, we examined the public datasets of breast invasive carcinoma from The Cancer Genome Atlas (TCGA). We analyzed the expression of SEMA6D along with its related genes, their functions, pathways, and potential as copredictors for BC patients’ survival. We found 6-gene expression profile that can be used as such predictors. Our study provides evidences for the first time that breast invasive carcinoma may contain a subtype based on SEMA6D expression. The expression of SEMA6D gene may play an important role in promoting patient survival, especially among triple negative breast cancer patients.

1. Introduction

Breast cancer (BC) is the second most common cancer diagnosed in American women and is also the second leading cause of cancer death in women [1, 2]. It is estimated that, in the developed world, one in eight women will develop breast cancer in her lifetime [3, 4]. BC lethality is mainly caused by metastasis, which accounts for approximately 90% of BC deaths [510]. Metastatic BC can be treated, sometimes for many years, but cannot be cured. BC primarily metastasizes to the bone, lungs, regional lymph nodes, liver, and brain, with the most common site being the bone. Research has focused heavily on BC metastasis for many years. Multiple signaling pathways, such as TGFβ, Wnt, Notch, and EGF, have been implicated in regulating metastasis of BC cells [510]. However, our knowledge of regulation of BC metastasis is far from complete. Identification of new factors that play critical roles in driving/inhibiting metastatic progression is an essential step toward fully understanding BC metastasis and will also provide novel therapeutic targets/reagents against BC. Our labs have focused on Semaphorins (SEMAs), especially SEMA6D. SEMAs were implicated in immune responses, heart development, and neurogenesis [1113] and recently in BC metastases [1417].

Semaphorins were initially recognized as phylogenetically conserved neuronal guidance cues, and their critical regulatory roles in BC metastasis have rapidly emerged in recent years. Based on their sequence similarity, Semaphorins are classified into eight classes: classes 1-2 are found in invertebrates, classes 3–7 comprise the vertebrate Semaphorins, and class V is encoded by viruses. Class 2, 3, and V Semaphorins are secreted, while all other members are membrane tethered through a single transmembrane domain [1725]. The signature structure of Semaphorins is the ~500 amino acid (aa) Sema domain, which is a variant of the β-propeller fold revealed from structural studies [22]. All Semaphorins, except class V, contain a PSI (Plexin-semaphorin-integrin) domain immediately to the C-terminal side of the Sema domain. Different classes of Semaphorins may also contain additional functional motifs.

Published studies regarding functions of Semaphorins in BC have mainly focused on class 3 secreted Semaphorins and SEMA4D [1417]. Our recent bioinformatic analysis by using public datasets reveals for the first time the potential role of SEMA6D in regulating BC pathogenesis. The extracellular region of SEMA6D, which contains the Sema and PSI domains, can be released from the cell surface to act as a secreted cytokine with unknown molecular mechanisms [26, 27]. Thus, SEMA6D may act both locally through cell-cell contacts and more distantly through diffusion of its cleaved ectodomain. The primary receptors of SEMA6D are PLXNA1 and PLXNA4, which belong to the Plexin receptor family. The extracellular domains of Plexins interact with Semaphorins, while the C-terminal tails of Plexins mediate intracellular signal transduction. Unlike many signaling pathways (such as the TGFβ/SMAD pathway), there is no “canonical” intracellular transduction cascade mediating activities of Semaphorin-Plexin signaling. Many intracellular signaling molecules, such as GTPase activating proteins, GTP/GDP exchange factors, and various tyrosine kinases, can be activated and/or inactivated by Semaphorin-Plexin signaling in a context-dependent manner [1825].

In this study, we examine the public datasets from The Cancer Genome Atlas (TCGA), National Cancer Institute (NCI) for expression of SEMA6D along with genes that interact with SEMA6D. Other genes coregulated with SEMA6D were analyzed for their function, pathway, and potential as copredictors for BC patients’ survival. We found 6-gene expression profile that can be used as such predictors. We also found that SEMA6D expression correlated with the cancer status of triple negative (TNBC) markers (ER, PR, and Her2 genes). The study shows the role of SEMA6D as potential survival predictor especially in TNBC patients.

2. Methods

2.1. Datasets

The Cancer Genome Atlas (TCGA) Data Portal was used to download breast invasive carcinoma (BRCA) samples (). The RNAseqV2 level 3 data, which includes fragments per kilobase of exon per million fragments mapped- (FPKM-) normalized gene level data, were used before statistics. In addition, idf and sdrf files were also downloaded for sample mapping and annotation. Clinical outcomes data were downloaded for correlation and survival analysis.

2.2. Gene Expression Data Analysis and Annotation

Gene-level normalized expression data were used in Partek Genomic Suite (PGS, St. Louis, MO) for additional normalization, statistics, and annotation. The analysis of variance (ANOVA) methods were used for group comparisons. False discovery rate (FDR) correction (Benjamini-Hochberg methods) was applied for multiple hypothesis testing purpose. Other statistical tools such as SAS (Cary, NC) and Ingenuity Pathway Analysis (IPA, Redwood City, CA) were used for pathway analysis and building gene-gene interaction network. Heatmap was generated by using hierarchical clustering methods after z-normalization.

2.3. Survival Analysis

A total of 140 patients with clinical outcomes data available (survival status, months of survival, demographics, and ER, PR, and HER2 status, etc.) were included in the analysis. Among significant genes after SEMA6D-high versus SEMA6D-low expression comparisons, we selected top 20 genes with the highest or lowest expression levels to correlate with clinical outcomes. Logarithm 2 based transformation of each gene was performed prior to any analysis. The correlation among these 20 genes was evaluated using Pearson correlation coefficient, and summary statistics were presented including mean with standard deviation, median, and range. Associations between level of genes and overall survival (OS) were assessed with Kaplan-Meier (K-M) curves and log-rank tests. Each gene was dichotomized as above or below median level of expression in the survival analysis. Significant association was determined at 5% type I error level. Multiple comparisons were not explicitly controlled for due to the small sample size and exploratory nature of the analysis.

3. Results and Discussion

Semaphorins, including members in subclass 3 and SEMA4D, have emerged as critical signaling molecules in regulating BC pathology [1417]. The potential roles of other members of the Semaphorin family in BC have not been well addressed in the literature. Our ongoing studies suggest that SEMA6D plays a critical role in mediating Bone Morphogenetic Protein (BMP) signaling to regulate epithelial-mesenchymal-transition (EMT) by endocardial cells in developing hearts (Kai Jiao et al.’s manuscript in preparation). EMT is an essential step for initiating metastasis in BC and other cancers. We thus decided to apply a bioinformatic approach to examine the potential role of SEMA6D in BC tumorigenesis and progression using publically available datasets. After testing a few smaller datasets form Gene Ontology Omnibus (GEO), we chose to use a breast invasive carcinoma (BRCA) dataset from TCGA as it represents a major public data source with clinical outcomes information, such as overall survival after diagnosis. We divided all samples () into three roughly equal sized groups based on SEMA6D expression (high, medium, and low). The comparison of SEMA6D-high versus SEMA6D-low expression group will reveal genes that are coregulated with SEMA6D.

Gene Expression Profile by Principle Components Analysis (PCA) and Hierarchical Clustering. To examine overall gene expression profile and sample similarities, we perform the PCA analysis of all samples. As showed in Figure 1, the PCA showed a clear separation among SEMA6D-high (H), SEMA6D-medium (M), and SEMA6D-low (L) groups. This indicates different gene expression profiles among the three groups.

Figure 1: Principle component analysis (PCA) of all samples.

Based on the genes that are differentially expressed in SEMA6D-high versus SEMA6D-low expression groups, we then performed a hierarchical clustering (Figure 2) using log2-transformed normalized signal intensities. The observed separation of samples in each group (SEMA6D-high, SEMA6D-medium, and SEMA6D-low) indicates different expression profiles among these groups.

Figure 2: Hierarchical clustering of significant genes of SEMA6D-H versus -L expression. Genes (vertical: high expression in red and low expression in green) and samples (horizontal: SEMA6D-high in green, SEMA6D-medium in blue, and SEMA6D-low in brown) were clustered based on Euclidean dissimilarity matrix.

Consistent with the PCA analysis, the high SEME6D expression samples showed a congregation in the lower part of the figure, which indicates a clear separation of samples based on SEMA6D expression. In other words, BC samples may contain a subtype with high SEMA6D expression.

We further examined SEMA6D levels by including SEMA6D-medium versus SEMA6D-low expression group comparison. As shown in the Venn diagram in Figure 3, 58 unique genes are significant between SEMA6D-medium and SEMA6D-low patients, while 2,357 genes are significant in SEMA6D-high versus SEMA6D-low expression comparison. This suggests that higher level of SEMA6D expression may lead to more significant changes, as evidenced by an increased number of genes with significant changes. This also suggests that the function of SEMA6D may be dependent on the expression level or that SEMA6D-induced functional effects are dose-dependent.

Figure 3: Number of significant genes between the two comparisons: H versus L and M versus L. FC: fold change.

Among significant genes of SEMA6D-high versus SEMA6D-low comparison, Gene Ontology (GO) analysis for biological processes revealed that multicellular organismal development and G-protein coupled receptor protein signaling pathway are among the top changed GO biological processes (Table 1). It has been shown that SEMA6D may play a role not only during heart development [26, 27] but also during development of retina [28] and axon in zebra fish [29]. Our results support these findings and showed that SEMA6D played a major role in organismal development in addition to the G-protein coupled receptor signaling (Table 1).

Table 1: Biological process: SEMA6D high versus low comparison.

The semaphorins and their receptors, the neuropilins and the Plexins, are constituents of a complex regulatory system that controls axonal guidance [30]. It was suggested that SEMA6D may bind to different receptor components and thus exert distinct functions during cardiac morphogenesis [26]. Our results suggested a broad function of SEMA6D to initiate signaling events that link to G-protein coupled receptor (GPCR) signaling (Table 1) and overall receptor activities (Table 2). GPCRs represent a super family of cell surface signaling proteins and play essential roles in cancer metastasis [2, 31, 32] and are one of the most promising targets of metastatic breast cancer therapy [3337].

Table 2: Molecular function: SEMA6D high versus low comparison.

In cell adhesion, cell-cell interactions between cancer cells with endothelium determine the metastatic spread. There are two major cell adhesions, including selectin and integrin, and accumulating evidence confirms that tumor cell interactions through them actively contribute to the metastatic spread of tumor cells [38]. Our results suggest a pivotal role for SEMA6D in tumor metastasis especially receptor activities (Table 2) and GPCR signaling (Table 1).

The GO-molecular functions also reveal that receptor activity, sequence-specific DNA binding, and voltage-gated sodium channel activities are among top affected molecular functions when SEMA6D level is high (Table 2). These results suggest that SEMA6D may initiate member receptors activation as a ligand. The gene-gene interactions among those genes that directly or indirectly interact with SEMA6D partially confirmed this hypothesis. As shown in Figure 4, elevated PLXNA4 may lead to an increase of SEMA6D expression and trigger transcriptions by the general transcription factors FOS and FOXO1. This may lead to a cascade of activations of membrane receptors including G-protein coupled receptors (Table 2).

Figure 4: Activation of SEMA6D and transcription. The gene-gene interaction network was built based on direct interactions by using Ingenuity Pathway Analysis (IPA) suite. Color indicates increased (in red) expression when SEMA6D-high samples were compared with SEMA6D-low samples. The number indicated the fold changes of this comparison.

As reported, Plexin-B1 is a receptor for the transmembrane semaphorin SEMA4D (CD100) [39], and PLXNA4 negatively regulates T lymphocyte responses [40]. It has been shown that SEMA6D induces NF-κB transcriptional activity in nonmalignant mesothelial cells [30]. Two potential targets of SEMA6D, the general transcription factors FOS and FOXO1, were both increased in SEME6D-high patients. FOXO1 has been widely reported in tumor oncogenesis and metastasis [4143]. This suggests an important role for SEMA6D in promoting general transcription through FOS coupled with FOXO1 as previously reported [44]. The balance of transcriptions of both tumor suppressors and oncogenes may be the key to understand the underlining mechanism.

SEMA6D and Tumor Metastasis. SEMA6D plays an important role in tissue development and differentiation, a process involving epithelial-mesenchymal-transition (EMT); it will be interesting to know if EMT-related genes are coregulated in SEMA6D-high patients. As shown in Table 3, major tumor metastatic promoter- (MMP-) 9 was dramatically reduced among SEMA6D-high samples, while several tumor metastatic promoters such as TGF-β-related factors, ZEB2, ZEB1, and GNG11, however, were elevated corresponding to a high level of SEMA6D. In addition, the expressions of all these genes (except for DSC2) are highly correlated with the expression of SEMA6D (Table 4). As SEMA6D was implicated in VEGF-dependent and anchorage-independent cell growth [30], we also included VEGF genes in the correlation analysis. High levels of correlations between SEMA6D level and VEGFs were found as well (Table 4) suggesting a role of VEGF family genes in mediating SEMA6D signaling.

Table 3: Expression of top EMT-related genes in SEMA6D-H versus L comparison.
Table 4: Correlation of SEMA6D with EMT gene expressions.

MMP family proteins, especially MMP9, were suggested to be involved in the process of metastasis of breast cancer to the brain [45], CD147-mediated metastasis in MCF7 cells [19], TGFβ-mediated signaling at the tumor-bone interface [46], and L2-mediated matrix remodeling in metastasis and mammary gland involution [47]. Decreased autocrine EGFR signaling in metastatic breast cancer cells inhibits tumor growth in bone and mammary fat pad through MMP9-dependent pathways [48]. By using an RNA interference approach, the reduced levels of MMP-9 mRNA and protein correlated with inhibited phenotype of tumor invasion and metastasis [14]. Our results are in line with these findings and suggest a tumor suppressor function for SEMA6D.

On the other hand, our results also showed an increased expression among SEMA6D-high samples of some important tumor promoters such as ZEB1/2, which had been reported to promote EMT by modulating Zeb1/2 and TGFβ expression [49]. Our results thus strongly suggest that the balance between tumor suppressors and promoters is the key to understand the role of SEMA6D during EMT. Another explanation is that the increased expression of SEMA6D may be the results, not the cause ZEB1/2 changes and vice versa.

SEMA6D Expression and Affected Signaling Pathways. Although roles of SEMAs have been suggested in breast cancer [1417], prostate cancer [50], and malignant mesothelioma [30], the underlying functional mechanisms including pathways are largely unknown. Nevertheless, SEMA6D has been reported to play a role in immune responses [12], NF-κB signaling [30], and stromal expression of SEMA6D [50]. Our results implicated top canonical pathways (Table 5), which partially confirmed previous reports such as cAMP-mediated signaling in cervical cancer cell migration [51] and in lung cancer [52], G-protein coupled receptor signaling in breast cancer [34, 53], and adhesion and diapedesis in a breast cancer cell line [54]. Therefore SEMA6D may play multiple roles during these processes although additional studies may be needed to further delineate SEMA6D functions in these pathways.

Table 5: Canonical signaling pathway by SEME6D high expression.

SEMA6D Expression Correlates with Patients’ Survival. In order to determine if SEMA6D and SEMA6D-related genes are correlated with overall patient survival, we filtered the significant gene list after SEMA6D-high versus SEMA6D-low expression comparison by choosing the top 10 most upregulated genes and top 10 most downregulated genes and conducted a survival analysis by using K-M methods. We found that 6 candidate genes were significantly associated with overall survival. These genes are SEMA6D, CLEC9A, C10orf107, DONSON, CHAC1, and TUBA1C. Two more genes COL4A6 and CBX2 genes are in borderline to be significant. A summary of these 8 genes is listed in Table 6 and Figure 5.

Table 6: Correlation of gene expression with patients’ survival.
Figure 5: SEMA6D correlates with patient survival.

In addition, increased expressions of SEMA6D, CLEC9A, COL4A6, and C10orf107 are associated with better survival while decreased expressions of DONSON, CHAC1, TUBA1C, and CBX2 also correlate to better survival. Figure 5 showed survival probability based on SEMA6D expression (≥median or ≤median expressions). Similar significant separation trends were also observed for both CLEC9A and C10orf107 as positive predictors and DONSON, CHAC1, and TUBA1C as negative predictors (data not shown).

Correlation of Expressions of SEMA6D and Other Genes with TNBC Status. As only 30% of women with metastases survive five years and virtually all TNBC women will ultimately die of their disease despite systemic therapy [55], we further explore the role of SEMA6D in promoting survival in TNBC patients. We found not only that these genes associated with survival, but also that they are interacting with TNBC status (Yes or No) significantly for SEMA6D, for example, with a log-rank , as shown in Figure 6. It is clearly shown that TNBC patients (Figure 6, SEMA6D-high in brown relative to SEMA6D-low in green) show larger survival differences as compared with non-TNBC patients (SEMA6D-high in red relative to SEMA6D-low in blue). Other genes such as CLEC9A () and C10orf107 () are similarly associated with TNBC status (data not shown).

Figure 6: Interaction of SEMA6D with TNBC in patients’ survival.

These results strongly suggest that SEMA6D expression levels correlate with overall survival (Figure 5), especially in TNBC patients (Figure 6).

4. Conclusions

Our study provides evidences that breast invasive carcinoma (BRCA) may contain a subtype based on SEMA6D expression. The expression of SEMA6D gene may play an important role in promoting patient survival, especially among triple negative breast cancer (TNBC) patients.

Abbreviations

TNBC:Triple negative breast cancer
BRCA:Breast invasive carcinoma
BC:Breast cancer
SEMA6D:Semaphorin 6D
BMP:Bone Morphogenetic Protein
EMT:Epithelial-mesenchymal-transition
PCA:Principle component analysis
GO:Gene Ontology
TCGA:The Cancer Genome Atlas.

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

Authors’ Contribution

Dongquan Chen and Yufeng Li contributed equally to the study.

Acknowledgments

The study was partially supported by institutional funding by University of Alabama at Birmingham (UAB) to Dongquan Chen and a Faculty Development Grant from UAB Comprehensive Cancer Center to Kai Jiao, an R01 (HL095783) to Kai Jiao, and an R21 (CA179282) to Lizhong Wang.

References

  1. C. Desantis, R. Siegel, P. Bandi, and A. Jemal, “Breast cancer statistics, 2011,” CA: Cancer Journal for Clinicians, vol. 61, no. 6, pp. 409–418, 2011. View at Publisher · View at Google Scholar · View at Scopus
  2. C. DeSantis, J. Ma, L. Bryan, and A. Jemal, “Breast cancer statistics, 2013,” CA Cancer Journal for Clinicians, vol. 64, no. 1, pp. 52–62, 2014. View at Publisher · View at Google Scholar · View at Scopus
  3. J. Ferlay, D. M. Parkin, and E. Steliarova-Foucher, “Estimates of cancer incidence and mortality in Europe in 2008,” European Journal of Cancer, vol. 46, no. 4, pp. 765–781, 2010. View at Publisher · View at Google Scholar · View at Scopus
  4. F. Lalloo and D. G. Evans, “Familial breast cancer,” Clinical Genetics, vol. 82, no. 2, pp. 105–114, 2012. View at Publisher · View at Google Scholar · View at Scopus
  5. M. Bendre, D. Gaddy, R. W. Nicholas, and L. J. Suva, “Breast cancer metastasis to bone: it is not all about PTHrP,” Clinical Orthopaedics and Related Research, no. 415, pp. S39–S45, 2003. View at Google Scholar · View at Scopus
  6. A. J. Redig and S. S. McAllister, “Breast cancer as a systemic disease: a view of metastasis,” Journal of Internal Medicine, vol. 274, pp. 113–126, 2013. View at Google Scholar
  7. O. J. Scully, B.-H. Bay, G. Yip, and Y. Yu, “Breast cancer metastasis,” Cancer Genomics & Proteomics, vol. 9, no. 5, pp. 311–320, 2012. View at Google Scholar · View at Scopus
  8. S. Byler, S. Goldgar, S. Heerboth et al., “Genetic and epigenetic aspects of breast cancer progression and therapy,” Anticancer Research, vol. 34, no. 3, pp. 1071–1077, 2014. View at Google Scholar · View at Scopus
  9. M. A. Blanco and Y. Kang, “Signaling pathways in breast cancer metastasis—novel insights from functional genomics,” Breast Cancer Research, vol. 13, no. 2, p. 206, 2011. View at Publisher · View at Google Scholar · View at Scopus
  10. J. Izrailit and M. Reedijk, “Developmental pathways in breast cancer and breast tumor-initiating cells: therapeutic implications,” Cancer Letters, vol. 317, no. 2, pp. 115–126, 2012. View at Publisher · View at Google Scholar · View at Scopus
  11. R. J. Pasterkamp, “Getting neural circuits into shape with semaphorins,” Nature Reviews Neuroscience, vol. 13, no. 9, pp. 605–618, 2012. View at Publisher · View at Google Scholar · View at Scopus
  12. B. P. O'Connor, S.-Y. Eun, Z. Ye et al., “Semaphorin 6D regulates the late phase of CD4+ T cell primary immune responses,” Proceedings of the National Academy of Sciences of the United States of America, vol. 105, no. 35, pp. 13015–13020, 2008. View at Publisher · View at Google Scholar · View at Scopus
  13. J. R. Leslie, F. Imai, K. Fukuhara et al., “Ectopic myelinating oligodendrocytes in the dorsal spinal cord as a consequence of altered semaphorin 6D signaling inhibit synapse formation,” Development, vol. 138, no. 18, pp. 4085–4095, 2011. View at Publisher · View at Google Scholar · View at Scopus
  14. L. Tamagnone, “Emerging role of semaphorins as major regulatory signals and potential therapeutic targets in cancer,” Cancer Cell, vol. 22, no. 2, pp. 145–152, 2012. View at Publisher · View at Google Scholar · View at Scopus
  15. M. Rehman and L. Tamagnone, “Semaphorins in cancer: biological mechanisms and therapeutic approaches,” Seminars in Cell & Developmental Biology, vol. 24, no. 3, pp. 179–189, 2013. View at Publisher · View at Google Scholar · View at Scopus
  16. C. Gu and E. Giraudo, “The role of semaphorins and their receptors in vascular development and cancer,” Experimental Cell Research, vol. 319, no. 9, pp. 1306–1316, 2013. View at Publisher · View at Google Scholar · View at Scopus
  17. G. C. Harburg and L. Hinck, “Navigating breast cancer: axon guidance molecules as breast cancer tumor suppressors and oncogenes,” Journal of Mammary Gland Biology and Neoplasia, vol. 16, no. 3, pp. 257–270, 2011. View at Publisher · View at Google Scholar · View at Scopus
  18. U. Yazdani and J. R. Terman, “The semaphorins,” Genome Biology, vol. 7, no. 3, article 211, 2006. View at Publisher · View at Google Scholar · View at Scopus
  19. T. S. Tran, A. L. Kolodkin, and R. Bharadwaj, “Semaphorin regulation of cellular morphology,” Annual Review of Cell and Developmental Biology, vol. 23, pp. 263–292, 2007. View at Publisher · View at Google Scholar · View at Scopus
  20. A. Casazza, P. Fazzari, and L. Tamagnone, “Semaphorin signals in cell adhesion and cell migration: functional role and molecular mechanisms,” Advances in Experimental Medicine and Biology, vol. 600, pp. 90–108, 2007. View at Publisher · View at Google Scholar · View at Scopus
  21. R. P. Kruger, J. Aurandt, and K.-L. Guan, “Semaphorins command cells to move,” Nature Reviews Molecular Cell Biology, vol. 6, no. 10, pp. 789–800, 2005. View at Publisher · View at Google Scholar · View at Scopus
  22. E. Gherardi, C. A. Love, R. M. Esnouf, and E. Y. Jones, “The sema domain,” Current Opinion in Structural Biology, vol. 14, no. 6, pp. 669–678, 2004. View at Publisher · View at Google Scholar · View at Scopus
  23. L. Tamagnone and P. M. Comoglio, “To move or not to move? Semaphorin signalling in cell migration,” EMBO Reports, vol. 5, no. 4, pp. 356–361, 2004. View at Publisher · View at Google Scholar · View at Scopus
  24. A. Ahmed and B. J. Eickholt, “Intracellular kinases in semaphorin signaling,” Advances in Experimental Medicine and Biology, vol. 600, pp. 24–37, 2007. View at Publisher · View at Google Scholar · View at Scopus
  25. T. Toyofuku and H. Kikutani, “Semaphorin signaling during cardiac development,” Advances in Experimental Medicine and Biology, vol. 600, pp. 109–117, 2007. View at Publisher · View at Google Scholar · View at Scopus
  26. T. Toyofuku, H. Zhang, A. Kumanogoh et al., “Dual roles of Sema6D in cardiac morphogenesis through region-specific association of its receptor, Plexin-A1, with off-track and vascular endothelial growth factor receptor type 2,” Genes and Development, vol. 18, no. 4, pp. 435–447, 2004. View at Publisher · View at Google Scholar · View at Scopus
  27. T. Toyofuku, H. Zhang, A. Kumanogoh et al., “Guidance of myocardial patterning in cardiac development by Sema6D reverse signalling,” Nature Cell Biology, vol. 6, no. 12, pp. 1204–1211, 2004. View at Publisher · View at Google Scholar · View at Scopus
  28. R. L. Matsuoka, L. O. Sun, K.-I. Katayama, Y. Yoshida, and A. L. Kolodkin, “Sema6B, Sema6C, and Sema6D expression and function during mammalian retinal development,” PLoS ONE, vol. 8, no. 4, Article ID e63207, 2013. View at Publisher · View at Google Scholar · View at Scopus
  29. M. Kimura, M. Taniguchi, Y. Mikami et al., “Identification and characterization of zebrafish semaphorin 6D,” Biochemical and Biophysical Research Communications, vol. 363, no. 3, pp. 762–768, 2007. View at Publisher · View at Google Scholar · View at Scopus
  30. A. Catalano, R. Lazzarini, S. D. Nuzzo, S. Orciari, and A. Procopio, “The plexin-A1 receptor activates vascular endothelial growth factor-receptor 2 and nuclear factor-kappaB to mediate survival and anchorage-independent growth of malignant mesothelioma cells,” Cancer Research, vol. 69, no. 4, pp. 1485–1493, 2009. View at Publisher · View at Google Scholar · View at Scopus
  31. S. Majumder, M. P. Sowden, S. A. Gerber et al., “G-protein-coupled receptor-2-interacting protein-1 is required for endothelial cell directional migration and tumor angiogenesis via cortactin-dependent lamellipodia formation,” Arteriosclerosis, Thrombosis, and Vascular Biology, vol. 34, no. 2, pp. 419–426, 2014. View at Publisher · View at Google Scholar · View at Scopus
  32. H.-Y. Li, X.-Y. Cui, W. Wu et al., “Pyk2 and Src mediate signaling to CCL18-induced breast cancer metastasis,” Journal of Cellular Biochemistry, vol. 115, no. 3, pp. 596–603, 2014. View at Publisher · View at Google Scholar · View at Scopus
  33. S.-H. Chan, W.-C. Huang, J.-W. Chang et al., “MicroRNA-149 targets GIT1 to suppress integrin signaling and breast cancer metastasis,” Oncogene, vol. 33, pp. 4496–4507, 2014. View at Publisher · View at Google Scholar · View at Scopus
  34. W. Wei, Z.-J. Chen, K.-S. Zhang et al., “The activation of G protein-coupled receptor 30 (GPR30) inhibits proliferation of estrogen receptor-negative breast cancer cells in vitro and in vivo,” Cell Death & Disease, vol. 5, Article ID e1428, 2014. View at Publisher · View at Google Scholar
  35. J.-M. Renoir, V. Marsaud, and G. Lazennec, “Estrogen receptor signaling as a target for novel breast cancer therapeutics,” Biochemical Pharmacology, vol. 85, no. 4, pp. 449–465, 2013. View at Publisher · View at Google Scholar · View at Scopus
  36. N. Dey, B. R. Smith, and B. Leyland-Jones, “Targeting basal-like breast cancers,” Current Drug Targets, vol. 13, no. 12, pp. 1510–1524, 2012. View at Publisher · View at Google Scholar · View at Scopus
  37. J.-M. Renoir, “Estradiol receptors in breast cancer cells: associated co-factors as targets for new therapeutic approaches,” Steroids, vol. 77, no. 12, pp. 1249–1261, 2012. View at Publisher · View at Google Scholar · View at Scopus
  38. H. Läubli and L. Borsig, “Selectins promote tumor metastasis,” Seminars in Cancer Biology, vol. 20, no. 3, pp. 169–177, 2010. View at Publisher · View at Google Scholar · View at Scopus
  39. L. Tamagnone, S. Artigiani, H. Chen et al., “Plexins are a large family of receptors for transmembrane, secreted, and GPI-anchored semaphorins in vertebrates,” Cell, vol. 99, no. 1, pp. 71–80, 1999. View at Google Scholar
  40. M. Yamamoto, K. Suzuki, T. Okuno et al., “Plexin-A4 negatively regulates T lymphocyte responses,” International Immunology, vol. 20, no. 3, pp. 413–420, 2008. View at Publisher · View at Google Scholar · View at Scopus
  41. J. Yang, T. Li, C. Gao et al., “FOXO1 3′UTR functions as a ceRNA in repressing the metastases of breast cancer cells via regulating miRNA activity,” FEBS Letters, vol. 588, no. 17, pp. 3218–3224, 2014. View at Publisher · View at Google Scholar
  42. X. Wang, C. Lin, X. Zhao et al., “Acylglycerol kinase promotes cell proliferation and tumorigenicity in breast cancer via suppression of the FOXO1 transcription factor,” Molecular Cancer, vol. 13, article 106, 2014. View at Publisher · View at Google Scholar · View at Scopus
  43. X. Feng, Z. Wu, Y. Wu et al., “Cdc25A regulates matrix metalloprotease 1 through Foxo1 and mediates metastasis of breast cancer cells,” Molecular and Cellular Biology, vol. 31, no. 16, pp. 3457–3471, 2011. View at Publisher · View at Google Scholar · View at Scopus
  44. D. V. Skarra, D. J. Arriola, C. A. Benson, and V. G. Thackray, “Forkhead box O1 is a repressor of basal and GnRH-induced Fshb transcription in gonadotropes,” Molecular Endocrinology, vol. 27, no. 11, pp. 1825–1839, 2013. View at Publisher · View at Google Scholar · View at Scopus
  45. O. Mendes, H.-T. Kim, and G. Stoica, “Expression of MMP2, MMP9 and MMP3 in breast cancer brain metastasis in a rat model,” Clinical and Experimental Metastasis, vol. 22, no. 3, pp. 237–246, 2005. View at Publisher · View at Google Scholar · View at Scopus
  46. K. C. Nannuru, M. Futakuchi, M. L. Varney, T. M. Vincent, E. G. Marcusson, and R. K. Singh, “Matrix metalloproteinase (MMP)-13 regulates mammary tumor-induced osteolysis by activating MMP9 and transforming growth factor-β signaling at the tumor-bone interface,” Cancer Research, vol. 70, no. 9, pp. 3494–3504, 2010. View at Publisher · View at Google Scholar · View at Scopus
  47. H. E. Barker, J. Chang, T. R. Cox et al., “LOXL2-mediated matrix remodeling in metastasis and mammary gland involution,” Cancer Research, vol. 71, no. 5, pp. 1561–1572, 2011. View at Publisher · View at Google Scholar · View at Scopus
  48. N. K. Nickerson, K. S. Mohammad, J. L. Gilmore et al., “Decreased autocrine EGFR signaling in metastatic breast cancer cells inhibits tumor growth in bone and mammary fat pad,” PLoS ONE, vol. 7, no. 1, Article ID e30255, 2012. View at Publisher · View at Google Scholar · View at Scopus
  49. L. Bakiri, S. Macho-Maschler, I. Custic et al., “Fra-1/AP-1 induces EMT in mammary epithelial cells by modulating Zeb1/2 and TGFβ expression,” Cell Death & Differentiation, vol. 22, pp. 336–350, 2015. View at Publisher · View at Google Scholar
  50. A. Henke, O. C. Grace, G. R. Ashley et al., “Stromal expression of decorin, Semaphorin6D, SPARC, Sprouty1 and Tsukushi in developing prostate and decreased levels of decorin in prostate cancer,” PLoS ONE, vol. 7, Article ID e42516, 2012. View at Publisher · View at Google Scholar · View at Scopus
  51. J. W. Lee, J. Lee, and E. Y. Moon, “HeLa human cervical cancer cell migration is inhibited by treatment with dibutyryl-cAMP,” Anticancer Research, vol. 34, no. 7, pp. 3447–3455, 2014. View at Google Scholar
  52. S. D. Kim, Y. J. Lee, J. S. Baik et al., “Baicalein inhibits agonist- and tumor cell-induced platelet aggregation while suppressing pulmonary tumor metastasis via cAMP-mediated VASP phosphorylation along with impaired MAPKs and PI3K-Akt activation,” Biochemical Pharmacology, vol. 92, no. 2, pp. 251–265, 2014. View at Publisher · View at Google Scholar
  53. S. Tao, H. He, Q. Chen, and W. Yue, “GPER mediated estradiol reduces miR-148a to promote HLA-G expression in breast cancer,” Biochemical and Biophysical Research Communications, vol. 451, no. 1, pp. 74–78, 2014. View at Publisher · View at Google Scholar
  54. C. Strell, K. Lang, B. Niggemann, K. S. Zaenker, and F. Entschladen, “Surface molecules regulating rolling and adhesion to endothelium of neutrophil granulocytes and MDA-MB-468 breast carcinoma cells and their interaction,” Cellular and Molecular Life Sciences, vol. 64, no. 24, pp. 3306–3316, 2007. View at Publisher · View at Google Scholar · View at Scopus
  55. I. A. Mayer, V. G. Abramson, B. D. Lehmann, and J. A. Pietenpol, “New strategies for triple-negative breast cancer—deciphering the heterogeneity,” Clinical Cancer Research, vol. 20, no. 4, pp. 782–790, 2014. View at Publisher · View at Google Scholar · View at Scopus