Navigating breast cancer: axon guidance molecules as breast cancer tumor suppressors and oncogenes

Multiple anti-tumor programs are activated by blocking BAD phosphorylation

The Bcl-2 family member BAD is a candidate disease modulator because it stimulates apoptosis in a cell context basis and inhibits cell migration during normal mammary gland morphogenesis. This activity depends on 3 key regulatory serines (S75, 99, 118) in the unphosphorylated state. Given that developmental programs are often hijacked in cancer, we hypothesized that BAD would impede breast cancer progression. We generated breast cancer mouse models representing loss-of-function or phosphorylation deficient mutations (PyMT-Bad-/- and PyMT-Bad3SA/3SA, respectively). Preventing BAD phosphorylation significantly decreased breast cancer progression and metastasis. The knock-out phenocopied the control PyMT-Bad+/+ suggesting that phosphorylated BAD protein was inert. Thus, the BAD3SA mutation unmasked latent anti-tumor activity. Indeed, transcriptomics showed PyMT-Bad3SA/3SA activated multiple anti-tumor programs including apoptosis, inflammation, cellular differentiation, and diminished cell migration. This anti-tumor effect associated with clinical survival of breast cancer patients whose tumors had high levels of unphosphorylated BAD. Kinase screens identified ERK as the major BAD kinase in breast cells, and ERK inhibition impeded tumoroid invasion. Our data suggest that unphosphorylated BAD modulates anti-tumor pathways that contribute to excellent patient prognosis. Thus, targeting ERK to dephosphorylate BAD may be an exciting therapeutic opportunity in the future.

The neuroscience in breast cancer: Current insights and clinical opportunities

The involvement of nerves in the development of breast cancer has emerged as a significant factor. Interaction between the nervous system and breast cancer can influence tumor initiation, growth, invasion, metastasis, reverse resistance to drugs, promote inflammation in tumors, and impair the immune system's ability to combat cancer. This review examined the intricate relationship linking the nervous system with breast cancer, emphasizing both central and peripheral aspects of the nervous system. Moreover, we reviewed neural cell factors and their impact on breast cancer progression, alongside the interactions between nerves and immunology, microbiota in breast cancer. Furthermore, the study discussed the potential of nerves as biomarkers for diagnosing and prognosticating breast cancer, and evaluated prospects for improving chemotherapy and immunotherapy therapeutic outcomes in breast cancer treatment. We hope to provide a deeper understanding of the neurobiological underpinnings of breast cancer and pave the way for the discovery of innovative therapeutic targets and prognostic markers.

Single-nuclei sequencing of uterine serous carcinoma reveals racial differences in immune signaling

Significant racial disparities exist between Black and White patients with uterine serous carcinoma (USC). While the reasons for these disparities are unclear, several studies have demonstrated significantly different rates of driver mutations between racial groups, including TP53. However, limited research has investigated the transcriptional differences of tumors or the composition of the tumor microenvironment (TME) between these groups. Here, we report the single-nuclei RNA-sequencing profiles of primary USC tumors from diverse racial backgrounds. We find that there are significant differences between the tumors of Black and White patients. Tumors from Black patients exhibited higher expression of specific genes associated with aggressiveness, such as PAX8, and axon guidance and synaptic signaling pathways. We also demonstrated that T cell populations are reduced in the tumor tissue compared to matched benign, while anti-inflammatory macrophage populations are retained within the TME. Furthermore, we investigated the connection between PAX8 overexpression and immunosuppression in USC through regulation of several cytokines and chemokines. Notably, we show that PAX8 activity can influence macrophage gene expression and protein secretion. These studies provide a detailed understanding of the USC transcriptome and TME, and identify differences in tumor biology from patients of different racial backgrounds.

The SEMA3F-NRP1/NRP2 axis is a key factor in the acquisition of invasive traits in in situ breast ductal carcinoma

BACKGROUND

A better understanding of ductal carcinoma in situ (DCIS) is urgently needed to identify these preinvasive lesions as distinct clinical entities. Semaphorin 3F (SEMA3F) is a soluble axonal guidance molecule, and its coreceptors Neuropilin 1 (NRP1) and NRP2 are strongly expressed in invasive epithelial BC cells.

METHODS

We utilized two cell line models to represent the progression from a healthy state to the mild-aggressive or ductal carcinoma in situ (DCIS) stage and, ultimately, to invasive cell lines. Additionally, we employed in vivo models and conducted analyses on patient databases to ensure the translational relevance of our results.

RESULTS

We revealed SEMA3F as a promoter of invasion during the DCIS-to-invasive ductal carcinoma transition in breast cancer (BC) through the action of NRP1 and NRP2. In epithelial cells, SEMA3F activates epithelialmesenchymal transition, whereas it promotes extracellular matrix degradation and basal membrane and myoepithelial cell layer breakdown.

CONCLUSIONS

Together with our patient database data, these proof-of-concept results reveal new SEMA3F-mediated mechanisms occurring in the most common preinvasive BC lesion, DCIS, and represent potent and direct activation of its transition to invasion. Moreover, and of clinical and therapeutic relevance, the effects of SEMA3F can be blocked directly through its coreceptors, thus preventing invasion and keeping DCIS lesions in the preinvasive state.

Dysregulated gene subnetworks in breast invasive carcinoma reveal novel tumor suppressor genes

Breast invasive carcinoma (BRCA) is the most malignant and leading cause of death in women. Global efforts are ongoing for improvement in early detection, prevention, and treatment. In this milieu, a comprehensive analysis of RNA-sequencing data of 1097 BRCA samples and 114 normal adjacent tissues is done to identify dysregulated genes in major molecular classes of BRCA in various clinical stages. Significantly enriched pathways in distinct molecular classes of BRCA have been identified. Pathways such as interferon signaling, tryptophan degradation, granulocyte adhesion & diapedesis, and catecholamine biosynthesis were found to be significantly enriched in Estrogen/Progesterone Receptor positive/Human Epidermal Growth Factor Receptor 2 negative, pathways such as RAR activation, adipogenesis, the role of JAK1/2 in interferon signaling, TGF-β and STAT3 signaling intricated in Estrogen/Progesterone Receptor negative/Human Epidermal Growth Factor Receptor 2 positive and pathways as IL-1/IL-8, TNFR1/TNFR2, TWEAK, and relaxin signaling were found in triple-negative breast cancer. The dysregulated genes were clustered based on their mutation frequency which revealed nine mutated clusters, some of which were well characterized in cancer while others were less characterized. Each cluster was analyzed in detail which led to the identification of NLGN3, MAML2, TTN, SYNE1, ANK2 as candidate genes in BRCA. They are central hubs in the protein-protein-interaction network, indicating their important regulatory roles. Experimentally, the Real-Time Quantitative Reverse Transcription PCR and western blot confirmed our computational predictions in cell lines. Further, immunohistochemistry corroborated the results in ~ 100 tissue samples. We could experimentally show that the NLGN3 & ANK2 have tumor-suppressor roles in BRCA as shown by cell viability assay, transwell migration, colony forming and wound healing assay. The cell viability and migration was found to be significantly reduced in MCF7 and MDA-MB-231 cell lines in which the selected genes were over-expressed as compared to control cell lines. The wound healing assay also demonstrated a significant decrease in wound closure at 12 h and 24 h time intervals in MCF7 & MDA-MB-231 cells. These findings established the tumor suppressor roles of NLGN3 & ANK2 in BRCA. This will have important ramifications for the therapeutics discovery against BRCA.

Effects of Eribulin on the RNA Content of Extracellular Vesicles Released by Metastatic Breast Cancer Cells

Extracellular vesicles (EVs) are small lipid particles secreted by almost all human cells into the extracellular space. They perform the essential function of cell-to-cell communication, and their role in promoting breast cancer progression has been well demonstrated. It is known that EVs released by triple-negative and highly aggressive MDA-MB-231 breast cancer cells treated with paclitaxel, a microtubule-targeting agent (MTA), promoted chemoresistance in EV-recipient cells. Here, we studied the RNA content of EVs produced by the same MDA-MB-231 breast cancer cells treated with another MTA, eribulin mesylate. In particular, we analyzed the expression of different RNA species, including mRNAs, lncRNAs, miRNAs, snoRNAs, piRNAs and tRNA fragments by RNA-seq. Then, we performed differential expression analysis, weighted gene co-expression network analysis (WGCNA), functional enrichment analysis, and miRNA-target identification. Our findings demonstrate the possible involvement of EVs from eribulin-treated cells in the spread of chemoresistance, prompting the design of strategies that selectively target tumor EVs.

The FLRT3-UNC5B checkpoint pathway inhibits T cell-based cancer immunotherapies

Cancers exploit coinhibitory receptors on T cells to escape tumor immunity, and targeting such mechanisms has shown remarkable clinical benefit, but in a limited subset of patients. We hypothesized that cancer cells mimic noncanonical mechanisms of early development such as axon guidance pathways to evade T cell immunity. Using gain-of-function genetic screens, we profiled axon guidance proteins on human T cells and their cognate ligands and identified fibronectin leucine-rich transmembrane protein 3 (FLRT3) as a ligand that inhibits T cell activity. We demonstrated that FLRT3 inhibits T cells through UNC5B, an axon guidance receptor that is up-regulated on activated human T cells. FLRT3 expressed in human cancers favored tumor growth and inhibited CAR-T and BiTE + T cell killing and infiltration in humanized cancer models. An FLRT3 monoclonal antibody that blocked FLRT3-UNC5B interactions reversed these effects in an immune-dependent manner. This study supports the concept that axon guidance proteins mimic T cell checkpoints and can be targeted for cancer immunotherapy.

Gene Signature Associated with Nervous System in an Experimental Radiation- and Estrogen-Induced Breast Cancer Model

Breast cancer is frequently the most diagnosed female cancer in the world. The experimental studies on cancer seldom focus on the relationship between the central nervous system and cancer. Despite extensive research into the treatment of breast cancer, chemotherapy resistance is an important issue limiting the efficacy of treatment. Novel biomarkers to predict prognosis or sensitivity to chemotherapy are urgently needed. This study examined nervous-system-related genes. The profiling of differentially expressed genes indicated that high-LET radiation, such as that emitted by radon progeny, in the presence of estrogen, induced a cascade of events indicative of tumorigenicity in human breast epithelial cells. Bioinformatic tools allowed us to analyze the genes involved in breast cancer and associated with the nervous system. The results indicated that the gene expression of the Ephrin A1 gene (EFNA1), the roundabout guidance receptor 1 (ROBO1), and the kallikrein-related peptidase 6 (KLK6) was greater in T2 and A5 than in the A3 cell line; the LIM domain kinase 2 gene (LIMK2) was greater in T2 than A3 and A5; the kallikrein-related peptidase 7 (KLK7), the neuroligin 4 X-linked gene (NLGN4X), and myelin basic protein (MBP) were greater than A3 only in T2; and the neural precursor cell expressed, developmentally down-regulated 9 gene (NEDD9) was greater in A5 than in the A3 and E cell lines. Concerning the correlation, it was found a positive correlation between ESR1 and EFNA1 in BRCA-LumA patients; with ROBO1 in BRCA-Basal patients, but this correlation was negative with the kallikrein-related peptidase 6 (KLK6) in BRCA-LumA and -LumB, as well as with LIMK2 and ROBO1 in all BRCA. It was also positive with neuroligin 4 X-linked (NLGN4X) in BRCA-Her2 and BRCA-LumB, and with MBP in BRCA-LumA and -LumB, but negative with KLK7 in all BRCA and BRCA-LumA and NEDD9 in BRCA-Her2. The differential gene expression levels between the tumor and adjacent tissue indicated that the ROBO1, KLK6, LIMK2, KLK7, NLGN4X, MBP, and NEDD9 gene expression levels were higher in normal tissues than in tumors; however, EFNA1 was higher in the tumor than the normal ones. EFNA1, LIMK2, ROBO1, KLK6, KLK7, and MBP gene expression had a negative ER status, whereas NEDD9 and NLGN4X were not significant concerning ER status. In conclusion, important markers have been analyzed concerning genes related to the nervous system, opening up a new avenue of studies in breast cancer therapy.

A parasite DNA binding protein with potential to influence disease susceptibility acts as an analogue of mammalian HMGA transcription factors

Intracellular pathogens construct their environmental niche, and influence disease susceptibility, by deploying factors that manipulate infected host cell gene expression. Theileria annulata is an important tick-borne parasite of cattle that causes tropical theileriosis. Excellent candidates for modulating host cell gene expression are DNA binding proteins bearing AT-hook motifs encoded within the TashAT gene cluster of the parasite genome. In this study, TashAT2 was transfected into bovine BoMac cells to generate three expressing and three non-expressing (opposite orientation) cell lines. RNA-Seq was conducted and differentially expressed (DE) genes identified. The resulting dataset was compared with genes differentially expressed between infected cells and non-infected cells, and DE genes between infected cell lines from susceptible Holstein vs tolerant Sahiwal cattle. Over 800 bovine genes displayed differential expression associated with TashAT2, 209 of which were also modulated by parasite infection. Network analysis showed enrichment of DE genes in pathways associated with cellular adhesion, oncogenesis and developmental regulation by mammalian AT-hook bearing high mobility group A (HMGA) proteins. Overlap of TashAT2 DE genes with Sahiwal vs Holstein DE genes revealed that a significant number of shared genes were associated with disease susceptibility. Altered protein levels encoded by one of these genes (GULP1) was strongly linked to expression of TashAT2 in BoMac cells and was demonstrated to be higher in infected Holstein leucocytes compared to Sahiwal. We conclude that TashAT2 operates as an HMGA analogue to differentially mould the epigenome of the infected cell and influence disease susceptibility.

DriverMP enables improved identification of cancer driver genes

BACKGROUND

Cancer is widely regarded as a complex disease primarily driven by genetic mutations. A critical concern and significant obstacle lies in discerning driver genes amid an extensive array of passenger genes.

FINDINGS

We present a new method termed DriverMP for effectively prioritizing altered genes on a cancer-type level by considering mutated gene pairs. It is designed to first apply nonsilent somatic mutation data, protein‒protein interaction network data, and differential gene expression data to prioritize mutated gene pairs, and then individual mutated genes are prioritized based on prioritized mutated gene pairs. Application of this method in 10 cancer datasets from The Cancer Genome Atlas demonstrated its great improvements over all the compared state-of-the-art methods in identifying known driver genes. Then, a comprehensive analysis demonstrated the reliability of the novel driver genes that are strongly supported by clinical experiments, disease enrichment, or biological pathway analysis.

CONCLUSIONS

The new method, DriverMP, which is able to identify driver genes by effectively integrating the advantages of multiple kinds of cancer data, is available at https://github.com/LiuYangyangSDU/DriverMP. In addition, we have developed a novel driver gene database for 10 cancer types and an online service that can be freely accessed without registration for users. The DriverMP method, the database of novel drivers, and the user-friendly online server are expected to contribute to new diagnostic and therapeutic opportunities for cancers.

Crosstalk between the peripheral nervous system and breast cancer influences tumor progression

Recent studies have shown that peripheral nerves play an important role in the progression of breast cancer. Breast cancer cells (BCCs) promote local peripheral nerve growth and branching by secreting neuroactive molecules, including neurotrophins and axon guidance molecules (AGMs). Sympathetic nerves promote breast cancer progression, while parasympathetic and sensory nerves mainly have anti-tumor effects in the progression of breast cancer. Specifically, peripheral nerves can influence the progression of breast cancer by secreting neurotransmitters not only directly binding to the corresponding receptors of BCCs, but also indirectly acting on immune cells to modulate anti-tumor immunity. In this review, we summarize the crosstalk between breast cancer and peripheral nerves and the roles of important neuroactive molecules in the progression of breast cancer. In addition, we summarize indicators, including nerve fiber density and perineural invasion (PNI), that may help determine the prognosis of breast cancer based on current research results, as well as potential therapeutic approaches, such as β-blockers and retroviral-mediated genetic neuroengineering techniques, that may enhance the prognosis of breast cancer. In addition, we propose suggestions for future research priorities based on a current lack of knowledge in this area.

The Value of EphB2 Receptor and Cognate Ephrin Ligands in Prognostic and Predictive Assessments of Human Breast Cancer

Cell–cell communication proteins Eph and ephrin constitute the largest family of receptor tyrosine kinases (RTKs). They are distinguished by the fact that both receptors and ligands are membrane-bound, and both can drive intracellular signaling in their respective cells. Ever since these RTKs have been found to be involved in cancer development, strategies to target them therapeutically have been actively pursued. However, before this goal can be rationally achieved, the contributions of either Eph receptors or their ephrin ligands to cancer development and progression should be scrutinized in depth. To assess the clinical pertinence of this concern, we performed a systematic review and meta-analysis of the prognostic/predictive value of EphB2 and its multiple cognate ephrin ligands in breast cancer. We found that EphB2 has prognostic value, as indicated by the association of higher EphB2 expression levels with lower distant metastasis-free survival (DMFS), and the association of lower EphB2 expression levels with poorer relapse-free survival (RFS). We also found that higher EphB2 expression could be a prognostic factor for distant metastasis, specifically in the luminal subtypes of breast cancer. EFNB2 showed a marked correlation between higher expression levels and shorter DMFS. EFNA5 or EFNB1 overexpression is correlated with longer RFS. Increased EFNB1 expression is correlated with longer OS in lymph node (LN)-negative patients and the luminal B subtype. Higher levels of EFNB2 or EFNA5 are significantly correlated with shorter RFS, regardless of LN status. However, while this correlation with shorter RFS is true for EFNB2 in all subtypes except basal, it is also true for EFNA5 in all subtypes except HER2+. The analysis also points to possible predictive value for EphB2. In systemically treated patients who have undergone either endocrine therapy or chemotherapy, we found that higher expression of EphB2 is correlated with better rates of RFS. Bearing in mind the limitations inherent to any mRNA-based profiling method, we complemented our analysis with an immunohistochemical assessment of expression levels of both the EphB2 receptor and cognate ephrin ligands. We found that the latter are significantly more expressed in cancers than in normal tissues, and even more so in invasive and metastatic samples than in ductal carcinoma in situ (DCIS). Finally, in an in vitro cellular model of breast cancer progression, based on H-Ras-transformation of the MCF10A benign mammary cell line, we observed dramatic increases in the mRNA expression of EphB2 receptor and EFNB1 and EFNB2 ligands in transformed and invasive cells in comparison with their benign counterparts. Taken together, these data show the clinical validity of a model whereby EphB2, along with its cognate ephrin ligands, have dual anti- and pro-tumor progression effects. In so doing, they reinforce the necessity of further biological investigations into Ephs and ephrins, prior to using them in targeted therapies.

Metabolomic and Transcriptomic Analysis of MCF-7 Cells Exposed to 23 Chemicals at Human-Relevant Levels: Estimation of Individual Chemical Contribution to Effects

BACKGROUND

Humans are constantly being exposed to various xenobiotics at relatively low concentrations. To date, limited evidence is available to ascertain whether a complex xenobiotic mixture at human-relevant levels causes any health effect. Moreover, there is no effective method to pinpoint the contribution of each chemical toward such an effect.

OBJECTIVES

This study aims to understand the responses of cells to a mixture containing 23 xenobiotics at human-relevant levels and develop a feasible method to decipher the chemical(s) that contribute significantly to the observed effect.

METHODS

We characterized the metabolome and transcriptome of breast cancer cells (MCF-7) before and after exposure to the mixture at human-relevant levels; preexposure levels were derived from existing large-scale biomonitoring data. A high-throughput metabolomics-based "leave-one-out" method was proposed to understand the relative contribution of each component by comparing the metabolome with and without the particular chemical in the mixture.

RESULTS

The metabolomic analysis suggested that the mixture altered metabolites associated with cell proliferation and oxidative stress. For the transcriptomes, gene ontology terms and pathways including "cell cycle," "cell proliferation," and "cell division" were significantly altered after mixture exposure. The mixture altered genes associated with pathways such as "genotoxicity" and "nuclear factor erythroid 2-related factor 2 (Nrf2)." Through joint pathways analysis, metabolites and genes were observed to be well-aligned in pyrimidine and purine metabolisms. The leave-one-out results showed that many chemicals made their contributions to specific metabolic pathways. The overall metabolome pattern of the absence of 2,4-dihyroxybenzophenone (DHB) or bisphenol A (BPA) showed great resemblance to controls, suggesting their higher relative contribution to the observed effect.

DISCUSSION

The omics results showed that exposure to the mixture at human-relevant levels can induce significant in vitro cellular changes. Also, the leave one out method offers an effective approach for deconvoluting the effects of the mixture. https://doi.org/10.1289/EHP6641.

Identification of Differentially Methylated Regions Associated with a Knockout of SUV39H1 in Prostate Cancer Cells

Epigenetic alterations, such as histone methylations, affect the pathogenesis of tumors including prostate cancer (PCa). Previously, we reported that metformin reduced SUV39H1, a histone methyltransferase of H3 Lys9, to inhibit the migration of PCa cells. Since histone methylation is functionally linked to DNA methylation, we speculate that the knockout of the SUV39H1 gene will affect the genomic DNA methylation profile to regulate PCa cell migration and invasion. The genome-wide DNA methylation level is lower in SUV39H1 knockout (KO) cells than wild-type (WT) ones. However, the methylation levels in functional regions of CpG Islands (CGI), 5' untranslated region (UTR5), and exon regions are higher in KO cells than WT cells. Analysis of differentially methylated regions (DMRs) identified 1241 DMR genes that have differential methylation on CG sites when comparing the KO and WT samples. Gene ontology enrichment and Kyoto Encyclopedia of Genes and Genomes Pathways analysis showed that knockout of SUV39H1 affects gene sets and pathways that are heavily involved in cell shapes, cell recognition, adhesion, motility, and migration. Our study suggests that SUV39H1 plays an important role in PCa migration via the epigenetic regulation of methylation on CG sites, and is a novel and legitimate target to inhibit PCa cell migration.

Small RNA deep sequencing reveals the expressions of microRNAs in ovine mammary gland development at peak-lactation and during the non-lactating period

MicroRNAs (miRNAs) are small non-coding RNAs that are involved in mammary gland development and lactation in livestock. Little is known about the roles of miRNAs in ovine mammary gland development, hence in this study the expression profiles of miRNAs of the mammary gland tissues of ewes at peak-lactation and during the non-lactating period were investigated using RNA sequencing. A total of 147 mature miRNAs were expressed in the two periods. Compared with peak-lactation, eight miRNAs in the non-lactating ewe mammary gland were significantly up-regulated, whereas fifteen miRNAs were down-regulated. A KEGG analysis revealed that the target genes of the up-regulated miRNAs were significantly enriched in lysosome, Wnt and MAPK signaling pathways, while the target genes of down-regulated miRNAs were significantly enriched in the PI3K-Akt signaling pathway, protein processing in endoplasmic reticulum and axon guidance. These results suggest that further study of the differentially expressed miRNAs could provide a better understanding of the molecular mechanisms of mammary development and lactation in sheep.

Construction of a 14-lncRNA risk score system predicting survival of children with acute myelocytic leukemia

Acute myelocytic leukemia (AML) is a frequent type of acute leukemia. The present study was performed to build a risk score system for the prognostic prediction of AML. AML RNA-sequencing data from samples from 111 children were downloaded from The Cancer Genome Atlas database. Using the DEseq and edgeR packages, the differentially expressed long non-coding RNAs (DE-lncRNAs) between bad and good prognosis groups were identified. A survival package was used to screen prognosis-associated lncRNAs and clinical factors. The optimal lncRNA combination was selected using the penalized package, and the risk-score system was built and evaluated. After the lncRNA-mRNA expression correlation network was constructed, the potential pathways involving the key lncRNAs were enriched using Gene Set Enrichment Analysis. Among the 61 DE-lncRNAs, 48 lncRNAs were significantly associated with prognosis. Relapse was an independent prognostic factor. The optimal 14-lncRNA risk score system was constructed. After 730 differentially expressed mRNAs were identified between the good and bad prognosis groups divided using a prognostic index, the lncRNA-mRNA expression correlation network was constructed. Enrichment analysis showed that semaphorin-3C [SEMA3C; regulated by probable leucine-tRNA ligase, mitochondrial (LARS2-AS1)] and secreted frizzled-related protein 5 [SFRP5; mediated by WASH complex subunit 5 (WASHC5)-antisense RNA 1 (AS1)] were involved in axon guidance and the Wnt signaling pathway, respectively. A 14-lncRNA (including paired box protein Pax8-AS1 and MYB AS1) risk-score system might be effective in predicting the prognosis of AML. Axon guidance (involving SEMA3C and LARS2-AS1) and the Wnt signaling pathway (involving SFRP5 and WASHC5-AS1) might be two important pathways affecting the prognosis of AML.

Pathway Mutations in Breast Cancer Using Whole-Exome Sequencing

The genomic landscape of breast cancer (BC) is complex. The purpose of this study was to decipher the mutational profiles of Taiwanese patients with BC using next-generation sequencing. We performed whole-exome sequencing on DNA from 24 tumor tissue specimens from BC patients. Sanger sequencing was used to validate the identified variants. Sanger sequencing was also performed on paired adjacent nontumor tissues. After genotype calling and algorithmic annotations, we identified 49 deleterious variants in canonical cancer-related genes in our BC cohort. The most frequently mutated genes were PIK3CA (16.67%), FKBP9 (12.5%), TP53 (12.5%), ATM (8.33%), CHEK2 (8.33%), FOXO3 (8.33%), NTRK1 (8.33%), and NUTM2B (8.33%). Seven mutated variants (ATR p.V1581fs, CSF1R p.R579Q, GATA3 p.T356delinsTMKS, LRP5 p.W389*, MAP3K1 p.T918fs, MET p.K1161fs, and MTR p.P1178S) were novel variants that are not present in any gene mutation database. After grouping the samples according to molecular subtype, we found that the cell cycle, MAPK, and chemokine signaling pathways in the luminal A subtype of BC; the focal adhesion, axon guidance, and endocytosis pathways in the luminal B subtype; and amyotrophic lateral sclerosis in the basal-like subtype were exclusively altered. Survival curve analysis showed that the presence of the MAPK signaling pathway and endocytosis mutations were correlated with a poor prognosis. These survival data were consistent with cBioPortal analyses of 2,051 BC cases. We discovered novel mutations in patients with BC. These results have implications for developing strategic, adjuvant, and gene-targeted therapies.

In Silico Approach for the Definition of radiomiRNomic Signatures for Breast Cancer Differential Diagnosis

Personalized medicine relies on the integration and consideration of specific characteristics of the patient, such as tumor phenotypic and genotypic profiling.

BACKGROUND

Radiogenomics aim to integrate phenotypes from tumor imaging data with genomic data to discover genetic mechanisms underlying tumor development and phenotype.

METHODS

We describe a computational approach that correlates phenotype from magnetic resonance imaging (MRI) of breast cancer (BC) lesions with microRNAs (miRNAs), mRNAs, and regulatory networks, developing a radiomiRNomic map. We validated our approach to the relationships between MRI and miRNA expression data derived from BC patients. We obtained 16 radiomic features quantifying the tumor phenotype. We integrated the features with miRNAs regulating a network of pathways specific for a distinct BC subtype.

RESULTS

We found six miRNAs correlated with imaging features in Luminal A (miR-1537, -205, -335, -337, -452, and -99a), seven miRNAs (miR-142, -155, -190, -190b, -1910, -3617, and -429) in HER2+, and two miRNAs (miR-135b and -365-2) in Basal subtype. We demonstrate that the combination of correlated miRNAs and imaging features have better classification power of Luminal A versus the different BC subtypes than using miRNAs or imaging alone.

CONCLUSION

Our computational approach could be used to identify new radiomiRNomic profiles of multi-omics biomarkers for BC differential diagnosis and prognosis.

A network-based pathway-extending approach using DNA methylation and gene expression data to identify altered pathways

Pathway analysis allows us to gain insights into a comprehensive understanding of the molecular mechanisms underlying cancers. Currently, high-throughput multi-omics data and various types of large-scale biological networks enable us to identify cancer-related pathways by comprehensively analyzing these data. Combining information from multidimensional data, pathway databases and interaction networks is a promising strategy to identify cancer-related pathways. Here we present a novel network-based approach for integrative analysis of DNA methylation and gene expression data to extend original pathways. The results show that the extension of original pathways can provide a basis for discovering new components of the original pathway and understanding the crosstalk between pathways in a large-scale biological network. By inputting the gene lists of the extended pathways into the classical gene set analysis (ORA and FCS), we effectively identified the altered pathways which are correlated well with the corresponding cancer. The method is evaluated on three datasets retrieved from TCGA (BRCA, LUAD and COAD). The results show that the integration of DNA methylation and gene expression data through a network of known gene interactions is effective in identifying altered pathways.

A temporal and spatial map of axons in developing mouse prostate

Prostate autonomic and sensory axons control glandular growth, fluid secretion, and smooth muscle contraction and are remodeled during cancer and inflammation. Morphogenetic signaling pathways reawakened during disease progression may drive this axon remodeling. These pathways are linked to proliferative activities in prostate cancer and benign prostate hyperplasia. However, little is known about which developmental signaling pathways guide axon investment into prostate. The first step in defining these pathways is pinpointing when axon subtypes first appear in prostate. We accomplished this by immunohistochemically mapping three axon subtypes (noradrenergic, cholinergic, and peptidergic) during fetal, neonatal, and adult stages of mouse prostate development. We devised a method for peri-prostatic axon density quantification and tested whether innervation is uniform across the proximo-distal axis of dorsal and ventral adult mouse prostate. Many axons directly interact with or innervate neuroendocrine cells in other organs, so we examined whether sensory or autonomic axons innervate neuroendocrine cells in prostate. We first detected noradrenergic, cholinergic, and peptidergic axons in prostate at embryonic day (E) 14.5. Noradrenergic and cholinergic axon densities are uniform across the proximal-distal axis of adult mouse prostate while peptidergic axons are denser in the periurethral and proximal regions. Peptidergic and cholinergic axons are closely associated with prostate neuroendocrine cells whereas noradrenergic axons are not. These results provide a foundation for understanding mouse prostatic axon development and organization and, provide strategies for quantifying axons during progression of prostate disease.

Integrating proteomic and phosphoproteomic data for pathway analysis in breast cancer

Background

As protein is the basic unit of cell function and biological pathway, shotgun proteomics, the large-scale analysis of proteins, is contributing greatly to our understanding of disease mechanisms. Proteomics study could detect the changes of both protein expression and modification. With the releases of large-scale cancer proteome studies, how to integrate acquired proteomic and phosphoproteomic data in more comprehensive pathway analysis becomes implemented, but remains challenging. Integrative pathway analysis at proteome level provides a systematic insight into the signaling network adaptations in the development of cancer.

Results

Here we integrated proteomic and phosphoproteomic data to perform pathway prioritization in breast cancer. We manually collected and curated breast cancer well-known related pathways from the literature as target pathways (TPs) or positive control in method evaluation. Three different strategies including Hypergeometric test based over-representation analysis, Kolmogorov-Smirnov (K-S) test based gene set analysis and topology-based pathway analysis, were applied and evaluated in integrating protein expression and phosphorylation. In comparison, we also assessed the ranking performance of the strategy using information of protein expression or protein phosphorylation individually. Target pathways were ranked more top with the data integration than using the information from proteomic or phosphoproteomic data individually. In the comparisons of pathway analysis strategies, topology-based method outperformed than the others. The subtypes of breast cancer, which consist of Luminal A, Luminal B, Basal and HER2-enriched, vary greatly in prognosis and require distinct treatment. Therefore we applied topology-based pathway analysis with integrating protein expression and phosphorylation profiles on four subtypes of breast cancer. The results showed that TPs were enriched in all subtypes but their ranks were significantly different among the subtypes. For instance, p53 pathway ranked top in the Basal-like breast cancer subtype, but not in HER2-enriched type. The rank of Focal adhesion pathway was more top in HER2- subtypes than in HER2+ subtypes. The results were consistent with some previous researches.

Conclusions

The results demonstrate that the network topology-based method is more powerful by integrating proteomic and phosphoproteomic in pathway analysis of proteomics study. This integrative strategy can also be used to rank the specific pathways for the disease subtypes.

Electronic supplementary material

The online version of this article (10.1186/s12918-018-0646-y) contains supplementary material, which is available to authorized users.

Analysis of miRNA signature differentially expressed in exosomes from adriamycin-resistant and parental human breast cancer cells

A major cause of failure in chemotherapy is drug resistance of cancer cells. Exosomes have been introduced to spread chemoresistance through delivering miRNAs. However, a systematic evaluation of the exosomal miRNA expression profiles responsible for chemoresistance is still lacking. In the present study, miRNA signature differentially expressed in exosomes derived from adriamycin-resistant (A/exo) and parental breast cancer cells (S/exo) were analyzed by microarray and the results were confirmed by PCR. A total of 309 miRNAs were increased and 66 miRNAs were decreased significantly in A/exo compared with S/exo. Specifically, 52 novel miRNAs with increased expression levels >16.0-fold in A/exo were identified. After prediction of target genes for 13 of 52 selected novel miRNAs, pathway analysis, gene ontology (GO) terms, and protein–protein interactions (PPIs) were constructed. The results implied that these selected exosomal miRNAs inhibited target genes involved in transcriptional misregulation in cancer, MAPK, and Wnt signaling pathways. Functional enrichment analysis demonstrated that the target genes were mainly responsible for protein phosphorylation, transcription regulation, molecular binding, and kinase activity. In summary, the current bioinformatics study of exosomal miRNAs may offer a new understanding into mechanisms of chemoresistance, which is helpful to find potential exosomal miRNAs to overcome drug insensitivity in future breast cancer treatment.

RNA Sequencing of Carboplatin- and Paclitaxel-Resistant Endometrial Cancer Cells Reveals New Stratification Markers and Molecular Targets for Cancer Treatment

Despite advances in surgical technique and adjuvant treatment, endometrial cancer has recently seen an increase in incidence and mortality in the USA. The majority of endometrial cancers can be cured by surgery alone or in combination with adjuvant chemo- or radiotherapy; however, a subset of patients experience recurrence for reasons that remain unclear. Recurrence is associated with chemoresistance to carboplatin and paclitaxel and consequentially, high mortality. Understanding the pathways involved in endometrial cancer chemoresistance is paramount for the identification of biomarkers and novel molecular targets for this disease. Here, we generated the first matched pairs of carboplatin-sensitive/carboplatin-resistant and paclitaxel-sensitive/paclitaxel-resistant endometrial cancer cells and subjected them to bulk RNA sequencing analysis. We found that 45 genes are commonly upregulated in carboplatin- and paclitaxel-resistant cells as compared to controls. Of these, the leukemia inhibitory factor, (LIF), the protein tyrosine phosphatase type IVA, member 3 (PTP4A3), and the transforming growth factor beta 1 (TGFB1) showed a highly significant correlation between expression level and endometrial cancer overall survival (OS) and can stratify the 545 endometrial cancer patients in the TCGA cohort into a high-risk and low-risk-cohorts. Additionally, four genes within the 45 upregulated chemoresistance-associated genes are ADAMTS5, MICAL2, STAT5A, and PTP4A3 codes for proteins for which small-molecule inhibitors already exist. We identified these proteins as molecular targets for chemoresistant endometrial cancer and showed that treatment with their correspondent inhibitors effectively killed otherwise chemoresistant cells. Collectively, these findings underline the utility of matched pair of chemosensitive and chemoresistant cancer cells to identify markers for endometrial cancer risk stratification and to serve as a pharmacogenomics model for identification of alternative chemotherapy approaches for treatment of patients with recurrent disease.

The Epigenetics of Triple-Negative and Basal-Like Breast Cancer: Current Knowledge

Breast cancer has the highest incidence among all malignancies diagnosed in women. Therapies have significantly improved over the years due to extensive molecular and clinical research; in a large number of cases, targeted therapies have provided better prognosis. However, one specific subtype remains elusive to targeted therapies–the triple-negative breast cancer. This immunohistochemically defined subtype is resistant to both endocrine and targeted therapies, leading to its poor prognosis. A field that is of great promise in current cancer research is epigenetics. By studying the epigenetic mechanisms underlying tumorigenesis–DNA methylation, histone modifications, and noncoding RNAs–advances in cancer treatment, diagnosis, and prevention are possible. This review aims to synthesize the epigenetic discoveries that have been made related to the triple-negative breast cancer.

Identification of pathway-based recurrence-associated signatures in optimally debulked patients with serous ovarian cancer

Serous ovarian cancer (SOC) is the most common form of the histological subtype of epithelial ovarian cancer, with the worst clinical outcome. Despite improvements in surgery and chemotherapy, most patients with SOC experience recurrence within 12-18 months of first-line treatment. Current studies are unable to robustly predict the recurrence of SOC, and more accurate predictive models are urgently required. We have, therefore, developed a novel pathway-structured model to predict the recurrence of SOC. We trained the model on a set of 333 patients and validated it in 3 diversified validation datasets of 403 patients. Genes significantly associated with recurrence within each pathway were identified using a Cox proportional hazards model based on LASSO estimation in the training dataset. Next, a pathway-structured scoring matrix was obtained after computation of the prognostic score for each pathway by fitting to the Cox proportional hazards model. With the pathway-structure scoring matrix as an input, the pathway-based recurrent signatures were identified using the Cox proportional hazards model based on LASSO estimation and the significant pathway-based signatures were externally validated in 3 independent datasets. Meanwhile, our pathway-structured model was compared with a commonly used gene-based model. Our results revealed that our 12 pathway-based signatures successfully predicted the recurrence of SOC with high accuracy in the training dataset and in the 3 validation datasets. Moreover, our pathway-structured model was superior to the gene-based model in 4 datasets. The pathways selected in our study will provide new insights into the pathogenesis and clinical treatments of SOC.

Identification of six key miRNAs associated with breast cancer through screening large-scale microarray data

Breast Cancer (BC) is one of the most common primary malignant tumors, which is life threatening. Previous studies have demonstrated that microRNAs (miRNA) may regulate or affect the incidence of BC. However, results of these studies are inconsistent, due to factors including the different sequencing platforms and sample selection methods used. To explore the key miRNAs involved in the pathogenesis of BC, and to use these miRNAs to monitor the tumor progression of BC, a systematic review was performed on the previous studies examining BC miRNA; the function of the target genes that were modulated by these key miRNAs were also analyzed. A total of 8 representative miRNA datasets examining the pathogenesis of BC were selected. Key miRNAs were identified by comparing the overlap between these datasets. Then, the target genes of these key miRNAs were predicted through TargetScan. Furthermore, functional enrichment analysis of target genes and transcription factor (TF) binding analysis was also performed using the Database for Annotation, Visualization and Integrated Discovery and Tfacts database, respectively. A total of 6 key miRNAs were identified by comparing the differentially expressed miRNAs datasets in the pathogenesis of BC. Compared with normal tissues, 3 miRNAs were upregulated: Hsa-miR-21b; hsa-miR-29b; and hsa-miR-155; and 3 miRNAs were downregulated: Hsa-miR-10b; hsa-miR-125; and hsa-miR-145. The target genes regulated by the up- and downregulated miRNAs were significantly enriched in the biological processes of ‘transcriptional regulation’, and these target genes depended on RNA polymerase II promoter and DNA template, respective to the up- and downregulated genes. The downregulated key miRNAs were specifically enriched in the biological processes of ‘ephrin receptor signaling pathway’ (GO: 0048013) and ‘axon guidance’ (GO: 0007411). TF analysis of the key miRNA target genes revealed that 104 TFs interacted with the 319 target genes of the upregulated miRNAs, while the 92 TFs interacted with the 254 target genes of the downregulated miRNAs. In total, there were 133 TFs and 63 (47.3%) TFs shared by the 2 types (up- and downregulated) of target genes. In summary, 6 key miRNAs in BC were identified by systematic review; the corresponding target genes and TFs that bind to these target genes were also identified, and the potential functions of target genes were revealed. These data may be beneficial to increasing the accuracy of BC treatment through monitoring miRNA.

Progesterone down-regulates SLIT/ROBO expression in mouse corpus luteum

BACKGROUND

Progesterone produced by the corpus luteum (CL) is essential for preparation, implantation and maintenance of gestation. Furthermore, progesterone plays a protective role against luteolysis in rodents. It has been reported that Slit/Robo family members expressed in the CL and involved in prostaglandin F_2α (PGF_2α) induced luteolysis. However, the interactions between progesterone and Slits/Robos in CL are not clear. This study was designed to examine whether or not luteolysis is regulated by the interaction of progesterone and Slits/Robos in mouse CL.

METHODS

In the current study, we used Real-time PCR to identify the effect of progesterone on Slit2/Robo1 expression in cultured luteal cells in vitro, and the exogenous progesterone injection on mouse luteolysis and Slit/Robo expression in vivo was studied via Real-time PCR and Western bolt.

RESULTS

Our in vitro experiment revealed that 1μM progesterone significantly decreased Slit2/Robo1 mRNA levels at 6h, 12h and 24h. Our in vivo experiment showed that the mRNA and protein levels of Slit2 and Robo1 decreased significantly 7days after progesterone supplement.

CONCLUSION

These findings indicate that progesterone maintains CL function and resists luteolysis possibly through down-regulating Slit/Robo signaling pathway in the CL.

UTX promotes hormonally responsive breast carcinogenesis through feed-forward transcription regulation with estrogen receptor

UTX is implicated in embryonic development and lineage specification. However, how this X-linked histone demethylase contributes to the occurrence and progression of breast cancer remains to be clarified. Here we report that UTX is physically associated with estrogen receptor (ER) and functions in ER-regulated transcription. We showed that UTX coordinates with JHDM1D and CBP to direct H3K27 methylation-acetylation transition and to create a permissive chromatin state on ER targets. Genome-wide analysis of the transcriptional targets of UTX by ChIP-seq identified a set of genes such as chemokine receptor CXCR4 that are intimately involved in breast cancer tumorigenesis and metastasis. We demonstrated that UTX promotes the proliferation and migration of ER⁺ breast cancer cells. Interestingly, UTX itself is transactivated by ER, forming a feed-forward loop in the regulation of hormone response. Indeed, UTX is upregulated during ER⁺ breast cancer progression, and the expression level of UTX is positively correlated with that of CXCR4 and negatively correlated with the overall survival of ER⁺ breast cancer patients. Our study identified a feed-forward loop between UTX and ER in the regulation of hormonally responsive breast carcinogenesis, supporting the pursuit of UTX as an emerging therapeutic target for the intervention of certain ER⁺ breast cancer with specific epigenetic vulnerability.

Genetic association of deleted in colorectal carcinoma variants with breast cancer risk: A case-control study

Deleted in colorectal carcinoma (DCC), a netrin-1 dependence receptor, is correlated with cell progression, migration, and adhesion. Evidence indicated that DCC was frequently down-regulated in many cancers. However, the association of DCC with breast cancer remains uncertain. We conducted a case-control study to investigate the impact of three DCC gene variants (rs2229080, rs7504990, and rs4078288) on breast cancer susceptibility in Chinese women. This study included 560 breast cancer patients and 583 age-matched healthy controls from Northwest China. The three gene variants were genotyped via Sequenom MassARRAY. Odds ratios (ORs) and 95% confidence intervals (CIs) were utilized to evaluate the associations. We found that individuals with the rs2229080 C/G, C/C, and C/G-CC genotypes had a higher breast cancer risk, and the minor allele C was associated with increased breast cancer risk in an allele model. We observed a significantly decreased breast cancer risk with the rs7504990 C/T, T/T, and C/T-T/T genotypes, and the minor allele T was protective against breast cancer in an allele model. In addition, rs2229080 was associated with the axillary lymph node (LN) metastasis status. An age-stratified analysis revealed an association between rs2229080 and reduced breast cancer risk among older patients (≥ 49 years). Furthermore, the haplotype analysis showed that the Crs2229080Crs7504990Ars4078288 haplotype was associated with a decreased breast cancer risk. However, the results indicated a lack of association between rs4078288 and breast cancer risk. These findings affirmed that rs2229080 and rs7504990 polymorphisms in DCC might be related with breast cancer susceptibility in Chinese women.

Loss of miR-203 regulates cell shape and matrix adhesion through ROBO1/Rac/FAK in response to stiffness

Breast tumor progression is accompanied by changes in the surrounding extracellular matrix (ECM) that increase stiffness of the microenvironment. Mammary epithelial cells engage regulatory pathways that permit dynamic responses to mechanical cues from the ECM. Here, we identify a SLIT2/ROBO1 signaling circuit as a key regulatory mechanism by which cells sense and respond to ECM stiffness to preserve tensional homeostasis. We observed that Robo1 ablation in the developing mammary gland compromised actin stress fiber assembly and inhibited cell contractility to perturb tissue morphogenesis, whereas SLIT2 treatment stimulated Rac and increased focal adhesion kinase activity to enhance cell tension by maintaining cell shape and matrix adhesion. Further investigation revealed that a stiff ECM increased Robo1 levels by down-regulating miR-203. Consistently, patients whose tumor expressed a low miR-203/high Robo1 expression pattern exhibited a better overall survival prognosis. These studies show that cells subjected to stiffened environments up-regulate Robo1 as a protective mechanism that maintains cell shape and facilitates ECM adherence.

MiRNA Transcriptome Profiling of Spheroid-Enriched Cells with Cancer Stem Cell Properties in Human Breast MCF-7 Cell Line

Breast cancer is the second leading cause of cancer-related mortality worldwide as most patients often suffer cancer relapse. The reason is often attributed to the presence of cancer stem cells (CSCs). Recent studies revealed that dysregulation of microRNA (miRNA) are closely linked to breast cancer recurrence and metastasis. However, no specific study has comprehensively characterised the CSC characteristic and miRNA transcriptome in spheroid-enriched breast cells. This study described the generation of spheroid MCF-7 cell in serum-free condition and the comprehensive characterisation for their CSC properties. Subsequently, miRNA expression differences between the spheroid-enriched CSC cells and their parental cells were evaluated using next generation sequencing (NGS). Our results showed that the MCF-7 spheroid cells were enriched with CSCs properties, indicated by the ability to self-renew, increased expression of CSCs markers, and increased resistance to chemotherapeutic drugs. Additionally, spheroid-enriched CSCs possessed greater cell proliferation, migration, invasion, and wound healing ability. A total of 134 significantly (p<0.05) differentially expressed miRNAs were identified between spheroids and parental cells using miRNA-NGS. MiRNA-NGS analysis revealed 25 up-regulated and 109 down-regulated miRNAs which includes some miRNAs previously reported in the regulation of breast CSCs. A number of miRNAs (miR-4492, miR-4532, miR-381, miR-4508, miR-4448, miR-1296, and miR-365a) which have not been previously reported in breast cancer were found to show potential association with breast cancer chemoresistance and self-renewal capability. The gene ontology (GO) analysis showed that the predicted genes were enriched in the regulation of metabolic processes, gene expression, DNA binding, and hormone receptor binding. The corresponding pathway analyses inferred from the GO results were closely related to the function of signalling pathway, self-renewability, chemoresistance, tumorigenesis, cytoskeletal proteins, and metastasis in breast cancer. Based on these results, we proposed that certain miRNAs identified in this study could be used as new potential biomarkers for breast cancer stem cell diagnosis and targeted therapy.

Old friends, new story: The role of Slit2C signaling through PlexinA1

Growth cone guidance is driven by attractive and repulsive signaling cues. Until recently, repulsive signaling by semaphorins was thought to be mediated through Plexin receptors, whereas Slits-induced repulsion was solely mediated through Robo receptors. In a recent report published in Nature Neuroscience, Celine Delloye-Bourgeois and colleagues (2015) combined phenotypic analyses of transgenic mouse lines and in vitro biochemical experiments to identify PlexinA1 as a novel receptor for Slits. Strikingly, they uncovered for the very first time that the Slit2C-terminal fragment possesses some unique biological activity as binding partner for PlexinA1. Even more excitingly, the signaling cascade triggered by SlitC binding to PlexinA1 mediates growth cone collapse of commissural axons both in vivo and ex vivo and nicely complements Robo-Slit signaling in the developing spinal cord midline to prevent midline recrossing.

An essential role for Gα(i2) in Smoothened-stimulated epithelial cell proliferation in the mammary gland

Hedgehog (Hh) signaling is critical for organogenesis, tissue homeostasis, and stem cell maintenance. The gene encoding Smoothened (SMO), the primary effector of Hh signaling, is expressed aberrantly in human breast cancer, as well as in other cancers. In mice that express a constitutively active form of SMO that does not require Hh stimulation in mammary glands, the cells near the transgenic cells proliferate and participate in hyperplasia formation. Although SMO is a seven-transmembrane receptor like G protein-coupled receptors (GPCRs), SMO-mediated activation of the Gli family of transcription factors is not known to involve G proteins. However, data from Drosophila and mammalian cell lines indicate that SMO functions as a GPCR that couples to heterotrimeric G proteins of the pertussis toxin (PTX)-sensitive Gαi class. Using genetically modified mice, we demonstrated that SMO signaling through G proteins occurred in the mammary gland in vivo. SMO-induced stimulation of proliferation was PTX-sensitive and required Gαi2, but not Gαi1, Gαi3, or activation of Gli1 or Gli2. Our findings show that activated SMO functions as a GPCR to stimulate proliferation in vivo, a finding that may have clinical importance because most SMO-targeted agents have been selected based largely on their ability to block Gli-mediated transcription.

Nerve fibers infiltrate the tumor microenvironment and are associated with nerve growth factor production and lymph node invasion in breast cancer

Infiltration of the tumor microenvironment by nerve fibers is an understudied aspect of breast carcinogenesis. In this study, the presence of nerve fibers was investigated in a cohort of 369 primary breast cancers (ductal carcinomas in situ, invasive ductal and lobular carcinomas) by immunohistochemistry for the neuronal marker PGP9.5. Isolated nerve fibers (axons) were detected in 28% of invasive ductal carcinomas as compared to only 12% of invasive lobular carcinomas and 8% of ductal carcinomas in situ (p = 0.0003). In invasive breast cancers, the presence of nerve fibers was observed in 15% of lymph node negative tumors and 28% of lymph node positive tumors (p = 0.0031), indicating a relationship with the metastatic potential. In addition, there was an association between the presence of nerve fibers and the expression of nerve growth factor (NGF) in cancer cells (p = 0.0001). In vitro, breast cancer cells were able to induce neurite outgrowth in PC12 cells, and this neurotrophic activity was partially inhibited by anti-NGF blocking antibodies. In conclusion, infiltration by nerve fibers is a feature of the tumor microenvironment that is associated with aggressiveness and involves NGF production by cancer cells. The potential participation of nerve fibers in breast cancer progression needs to be further considered.

SEMA6D Expression and Patient Survival in Breast Invasive Carcinoma

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.

The Neural Stem Cell Microenvironment: Focusing on Axon Guidance Molecules and Myelin-Associated Factors

Neural stem cells (NSCs) could produce various cell phenotypes in the subventricular zone (SVZ) and dentate gyrus of the hippocampus in the central nervous system (CNS), where neurogenesis has been determined to occur. The extracellular microenvironment also influences the behaviors of NSCs during development and at CNS injury sites. Our previous study indicates that myelin, a component of the CNS, could regulate the differentiation of NSCs in vitro. Recent reports have implicated three myelin-derived inhibitors, NogoA, myelin-associated glycoprotein (MAG), and oligodendrocyte-myelin glycoprotein (OMgp), as well as several axon guidance molecules as regulators of NSC survival, proliferation, migration, and differentiation. However, the molecular mechanisms underlying the behavior of NSCs are not fully understood. In this study, we summarize the current literature on the effects of different extrinsic factors on NSCs and discuss possible mechanisms, as well as future possible clinical applications.

A co-culture genome-wide RNAi screen with mammary epithelial cells reveals transmembrane signals required for growth and differentiation

INTRODUCTION

The extracellular signals regulating mammary epithelial cell growth are of relevance to understanding the pathophysiology of mammary epithelia, yet they remain poorly characterized. In this study, we applied an unbiased approach to understanding the functional role of signalling molecules in several models of normal physiological growth and translated these results to the biological understanding of breast cancer subtypes.

METHODS

We developed and utilized a cytogenetically normal clonal line of hTERT immortalized human mammary epithelial cells in a fibroblast-enhanced co-culture assay to conduct a genome-wide small interfering RNA (siRNA) screen for evaluation of the functional effect of silencing each gene. Our selected endpoint was inhibition of growth. In rigorous postscreen validation processes, including quantitative RT-PCR, to ensure on-target silencing, deconvolution of pooled siRNAs and independent confirmation of effects with lentiviral short-hairpin RNA constructs, we identified a subset of genes required for mammary epithelial cell growth. Using three-dimensional Matrigel growth and differentiation assays and primary human mammary epithelial cell colony assays, we confirmed that these growth effects were not limited to the 184-hTERT cell line. We utilized the METABRIC dataset of 1,998 breast cancer patients to evaluate both the differential expression of these genes across breast cancer subtypes and their prognostic significance.

RESULTS

We identified 47 genes that are critically important for fibroblast-enhanced mammary epithelial cell growth. This group was enriched for several axonal guidance molecules and G protein-coupled receptors, as well as for the endothelin receptor PROCR. The majority of genes (43 of 47) identified in two dimensions were also required for three-dimensional growth, with HSD17B2, SNN and PROCR showing greater than tenfold reductions in acinar formation. Several genes, including PROCR and the neuronal pathfinding molecules EFNA4 and NTN1, were also required for proper differentiation and polarization in three-dimensional cultures. The 47 genes identified showed a significant nonrandom enrichment for differential expression among 10 molecular subtypes of breast cancer sampled from 1,998 patients. CD79A, SERPINH1, KCNJ5 and TMEM14C exhibited breast cancer subtype-independent overall survival differences.

CONCLUSION

Diverse transmembrane signals are required for mammary epithelial cell growth in two-dimensional and three-dimensional conditions. Strikingly, we define novel roles for axonal pathfinding receptors and ligands and the endothelin receptor in both growth and differentiation.

Genome-wide DNA methylation patterns in pancreatic ductal adenocarcinoma reveal epigenetic deregulation of SLIT-ROBO, ITGA2 and MET signaling

The importance of epigenetic modifications such as DNA methylation in tumorigenesis is increasingly being appreciated. To define the genome-wide pattern of DNA methylation in pancreatic ductal adenocarcinomas (PDAC), we captured the methylation profiles of 167 untreated resected PDACs and compared them to a panel of 29 adjacent nontransformed pancreata using high-density arrays. A total of 11,634 CpG sites associated with 3,522 genes were significantly differentially methylated (DM) in PDAC and were capable of segregating PDAC from non-malignant pancreas, regardless of tumor cellularity. As expected, PDAC hypermethylation was most prevalent in the 5' region of genes (including the proximal promoter, 5'UTR and CpG islands). Approximately 33% DM genes showed significant inverse correlation with mRNA expression levels. Pathway analysis revealed an enrichment of aberrantly methylated genes involved in key molecular mechanisms important to PDAC: TGF-β, WNT, integrin signaling, cell adhesion, stellate cell activation and axon guidance. Given the recent discovery that SLIT-ROBO mutations play a clinically important role in PDAC, the role of epigenetic perturbation of axon guidance was pursued in more detail. Bisulfite amplicon deep sequencing and qRT-PCR expression analyses confirmed recurrent perturbation of axon guidance pathway genes SLIT2, SLIT3, ROBO1, ROBO3, ITGA2 and MET and suggests epigenetic suppression of SLIT-ROBO signaling and up-regulation of MET and ITGA2 expression. Hypomethylation of MET and ITGA2 correlated with high gene expression, which was associated with poor survival. These data suggest that aberrant methylation plays an important role in pancreatic carcinogenesis affecting core signaling pathways with potential implications for the disease pathophysiology and therapy.

Exosomes from docetaxel-resistant breast cancer cells alter chemosensitivity by delivering microRNAs

Breast cancer (BCa) remains chemo-unresponsive by inevitable progression of resistance to first-line treatment with docetaxel (doc). Emerging studies indicate that exosomes act as mediators of intercellular communication between heterogeneous populations of tumor cells, engendering a transmitted drug resistance for cancer development. Such modulatory effects have been related to the constant shuttle of biologically active molecules including microRNAs (miRNAs). Here, we aimed to investigate the relevance of exosome-mediated miRNA delivery in resistance transmission of BCa subpopulations. Using microarray and polymerase chain reaction, we found that exosomes from doc-resistant BCa cells (D/exo) loaded cellular miRNAs. Following D/exo transfer to the fluorescent sensitive cells (GFP-S), some miRNAs were significantly increased in recipient GFP-S. Target gene prediction and pathway analysis revealed the involvement of the top 20 most abundant miRNAs of D/exo in pathways implicated in therapy failure. Coculture assays showed that miRNA-containing D/exo increased the overall resistance of GFP-S to doc exposure. Moreover, D/exo was able to alter gene expression in GFP-S. Our results open up an intriguing possibility that drug-resistant BCa cells may spread chemoresistance to sensitive ones by releasing exosomes and that the effects could be partly attributed to the intercellular transfer of specific miRNAs.

Genome-wide pathway analysis of breast cancer

The aim of this study was to identify candidate single-nucleotide polymorphisms (SNPs) that might affect susceptibility to breast cancer and then elucidate their potential mechanisms and generate SNP-to-gene-to-pathway hypotheses. A genome-wide association study (GWAS) dataset of breast cancer that included 453,852 SNPs from 1,145 breast cancer patients and 1,142 control subjects of European descent was used in this study. The identify candidate causal SNPs and pathways (ICSNPathway) method was applied to the GWAS dataset. ICSNPathway analysis identified 16 candidate SNPs, 13 genes, and 7 pathways, which together revealed 7 hypothetical biological mechanisms. The strongest hypothetical biological mechanism was that rs3168891 and rs2899849 alter the role of MBIP in the inactivation of mitogen-activated protein kinase (MAPK) (p < 0.001; false discovery rate (FDR) = 0.038). The second strongest mechanism was that rs2229714 modulates RPS6KA1 to affect its role in growth hormone signaling (p = 0.001; FDR = 0.039). The third strongest mechanism was that rs2230394 modulates ITGB1 to regulate the PTEN pathway and hsa04360 (axon guidance pathway) (p < 0.001; FDR = 0.039, 0.041). Use of the ICSNPathway to analyze breast cancer GWAS data identified 16 candidate SNPs, 13 genes (including MBIP, RPS6KA1, and ITGB1), and 7 pathways that might contribute to the susceptibility of patients to breast cancer.

Exosomes from drug-resistant breast cancer cells transmit chemoresistance by a horizontal transfer of microRNAs

Adriamycin and docetaxel are two agents commonly used in treatment of breast cancer, but their efficacy is often limited by the emergence of chemoresistance. Recent studies indicate that exosomes act as vehicles for exchange of genetic cargo between heterogeneous populations of tumor cells, engendering a transmitted drug resistance for cancer development and progression. However, the specific contribution of breast cancer-derived exosomes is poorly understood. Here we reinforced other's report that human breast cancer cell line MCF-7/S could acquire increased survival potential from its resistant variants MCF-7/Adr and MCF-7/Doc. Additionally, exosomes of the latter, A/exo and D/exo, significantly modulated the cell cycle distribution and drug-induced apoptosis with respect to S/exo. Exosomes pre-treated with RNase were unable to regulate cell cycle and apoptosis resistance, suggesting an RNA-dependent manner. Microarray and polymerase chain reaction for the miRNA expression profiles of A/exo, D/exo, and S/exo demonstrated that they loaded selective miRNA patterns. Following A/exo and D/exo transfer to recipient MCF-7/S, the same miRNAs were significantly increased in acquired cells. Target gene prediction and pathway analysis showed the involvement of miR-100, miR-222, and miR-30a in pathways implicated in cancer pathogenesis, membrane vesiculation and therapy failure. Furthermore, D/exo co-culture assays and miRNA mimics transfection experiments indicated that miR-222-rich D/exo could alter target gene expression in MCF-7/S. Our results suggest that drug-resistant breast cancer cells may spread resistance capacity to sensitive ones by releasing exosomes and that such effects could be partly attributed to the intercellular transfer of specific miRNAs.

Semaphorin7A promotes tumor growth and exerts a pro-angiogenic effect in macrophages of mammary tumor-bearing mice

Semaphorins are a large family of molecules involved in axonal guidance during the development of the nervous system and have been recently shown to have both angiogenic and anti-angiogenic properties. Specifically, semaphorin 7A (SEMA7A) has been reported to have a chemotactic activity in neurogenesis and to be an immune modulator through α1β1integrins. SEMA7A has been shown to promote monocyte chemotaxis and induce them to produce proinflammatory mediators. In this study we explored the role of SEMA7A in a murine model of breast cancer. We show that SEMA7A is highly expressed by DA-3 murine mammary tumor cells in comparison to normal mammary cells (EpH4), and that peritoneal elicited macrophages from mammary tumor-bearing mice also express SEMA7A at higher levels compared to those derived from normal mice. We also show that murine macrophages treated with recombinant murine SEMA7A significantly increased their expression of proangiogenic molecule CXCL2/MIP-2. Gene silencing of SEMA7A in peritoneal elicited macrophages from DA-3 tumor-bearing mice resulted in decreased CXCL2/MIP-2 expression. Mice implanted with SEMA7A silenced tumor cells showed decreased angiogenesis in the tumors compared to the wild type tumors. Furthermore, peritoneal elicited macrophages from mice bearing SEMA7A-silenced tumors produce significantly (p < 0.01) lower levels of angiogenic proteins, such as CXCL2/MIP-2, CXCL1, and MMP-9, compared to those from control DA-3 mammary tumors. We postulate that SEMA7A in mammary carcinomas may skew monocytes into a pro-tumorigenic phenotype to support tumor growth. SEMA7A could prove to be valuable in establishing new research avenues toward unraveling important tumor-host immune interactions in breast cancer patients.

Microarray expression profile analysis of long non-coding RNAs of advanced stage human gastric cardia adenocarcinoma

Gastric cardia adenocarcinoma (GCA) is a unique malignant tumor for its characteristics different from gastric and esophageal cancer epidemiologically and pathologically. The incidence of GCA has steadily increased for the last three decades and many patients are diagnosed with advanced stage because of the lack of typical and obvious symptoms at an early stage. To gain insight into the molecular mechanisms of GCA of advanced stage, we investigated the microarray expression profile of long non-coding RNAs of 12 advanced stage GCA patients. Long non-coding RNAs (lncRNAs) lack protein-coding potential and are over 200 bp in length. LncRNAs are known to be involved in the multifactor and multistep processes of tumor development and metastasis. In this study, we performed lncRNA transcriptome profiling of GCA biopsy tissue from 12 GCA patients who were confirmed by pathology to have developed lymph node metastasis and 12 paired non-cancerous gastric cardia tissues to determine if a gene expression profile unique to the lymph node metastasis group could be detected. Comparison of differentially expressed transcripts between the groups identified eight pathways that corresponded to down-regulated transcripts and 18 pathways that corresponded to up-regulated transcripts (p value cut-off 0.05). Gene ontology analysis showed that the up-regulated transcripts were most highly enriched in SRP-dependent cotranslational protein targeting to membrane, cytosolic ribosome, and structural constituent of ribosome, and the down-regulated transcripts were highly enriched in carboxylic acid transport, focal adhesion, and cation binding. This study shows that lncRNAs dysregulation exerts important roles in human GCA lymph node metastasis, indicating that lncRNAs are novel candidate biomarkers for the clinical diagnosis of advanced stage GCA and that could be targets for further therapy.

Role of the EphB2 receptor in autophagy, apoptosis and invasion in human breast cancer cells

The Eph and Ephrin proteins, which constitute the largest family of receptor tyrosine kinases, are involved in normal tissue development and cancer progression. Here, we examined the expression and role of the B-type Eph receptor EphB2 in breast cancers. By immunohistochemistry using a progression tissue microarray of human clinical samples, we found EphB2 to be expressed in benign tissues, but strongly increased in cancers particularly in invasive and metastatic carcinomas. Subsequently, we found evidence that EphB2, whose expression varies in established cell breast lines, possesses multiple functions. First, the use of a DOX-inducible system to restore EphB2 function to low expressers resulted in decreased tumor growth in vitro and in vivo, while its siRNA-mediated silencing in high expressers increased growth. This function involves the onset of apoptotic death paralleled by caspases 3 and 9 activation. Second, EphB2 was also found to induce autophagy, as assessed by immunofluorescence and/or immunoblotting examination of the LC3, ATG5 and ATG12 markers. Third, EphB2 also has a pro-invasive function in breast cancer cells that involves the regulation of MMP2 and MMP9 metalloproteases and can be blocked by treatment with respective neutralizing antibodies. Furthermore, EphB2-induced invasion is kinase-dependent and is impeded in cells expressing a kinase-dead mutant EphB2. In summary, we identified a mechanism involving a triple role for EphB2 in breast cancer progression, whereby it regulates apoptosis, autophagy, and invasion.

Stage-specific methylome screen identifies that NEFL is downregulated by promoter hypermethylation in breast cancer

Aberrant hypermethylation of promoter regions in specific genes is a key event in the formation and progression of breast cancers, and an increasing number of marker genes have been identified. However, few genes which show methylation change in accordance with the progression of breast cancer have been identified. To identify genes which consistently undergo promoter methylation alterations as the tumor develops from a benign to a malignant form, genome-wide methylation databases of breast cancer cell lines from stage I to stage IV were analyzed. Heatmap and cluster analysis revealed that the genome-wide methylation changes showed a good accordance with tumor progression. Seven out of 14,495 genes were found to be consistently increased alongside the promoter methylation level through the normal cell line to the cancer stage IV cell lines. NEFL, one of the in silico hypermethylated genes in cancer, showed hypermethylation and lower expression in the cancer cell line MDA-MB-231, as well as in cancer tissues (methylation, p<0.05; expression, p<0.01). The expression was restored by inducing demethylation of the promoter in MDA-MB-231 cells. Our findings may lend credence to the possibility of using tumor stage-specific alterations in methylation patterns as biomarkers for estimating prognosis and assessing treatment options for breast cancer.

Support Vector Machine classifier for estrogen receptor positive and negative early-onset breast cancer

Two major breast cancer sub-types are defined by the expression of estrogen receptors on tumour cells. Cancers with large numbers of receptors are termed estrogen receptor positive and those with few are estrogen receptor negative. Using genome-wide single nucleotide polymorphism genotype data for a sample of early-onset breast cancer patients we developed a Support Vector Machine (SVM) classifier from 200 germline variants associated with estrogen receptor status (p<0.0005). Using a linear kernel Support Vector Machine, we achieved classification accuracy exceeding 93%. The model indicates that polygenic variation in more than 100 genes is likely to underlie the estrogen receptor phenotype in early-onset breast cancer. Functional classification of the genes involved identifies enrichment of functions linked to the immune system, which is consistent with the current understanding of the biological role of estrogen receptors in breast cancer.

Riding the crest of the wave: parallels between the neural crest and cancer in epithelial-to-mesenchymal transition and migration

The neural crest (NC) is first induced as an epithelial population of cells at the neural plate border requiring complex signaling between bone morphogenetic protein, Wnt, and fibroblast growth factors to differentiate the neural and NC fate from the epidermis. Remarkably, following induction, these cells undergo an epithelial-to-mesenchymal transition (EMT), delaminate from the neural tube, and migrate through various tissue types and microenvironments before reaching their final destination where they undergo terminal differentiation. This process is mirrored in cancer metastasis, where a primary tumor will undergo an EMT before migrating and invading other cell populations to create a secondary tumor site. In recent years, as our understanding of NC EMT and migration has deepened, important new insights into tumorigenesis and metastasis have also been achieved. These discoveries have been driven by the observation that many cancers misregulate developmental genes to reacquire proliferative and migratory states. In this review, we examine how the NC provides an excellent model for studying EMT and migration. These data are discussed from the perspective of the gene regulatory networks that control both NC and cancer cell EMT and migration. Deciphering these processes in a comparative manner will expand our knowledge of the underlying etiology and pathogenesis of cancer and promote the development of novel targeted therapeutic strategies for cancer patients. © 2013 Wiley Periodicals, Inc.

Role of the nervous system in cancer metastasis

The notion that tumors lack innervation was proposed several years ago. However, nerve fibers are irregulatedly found in some tumor tissues. Their terminals interaction with cancer cells are considered to be neuro-neoplastic synapses. Moreover, neural-related factors, which are important players in the development and activity of the nervous system, have been found in cancer cells. Thus, they establish a direct connection between the nervous system and tumor cells. They modulate the process of metastasis, including degradation of base membranes, cancer cell invasion, migration, extravasation and colonization. Peripheral nerve invasion provides another pathway for the spread of cancer cells when blood and lymphatic metastases are absent, which is based on the interactions between the microenvironments of nerve fibers and tumor cells. The nervous system also modulates angiogenesis, the tumor microenvironment, bone marrow, immune functions and inflammatory pathways to influence metastases. Denervation of the tumor has been reported to enhance cancer metastasis. Stress, social isolation and other emotional factors may increase distant metastasis through releasing hormones from the brain, the hypothalamic-pituitary-adrenal axis and autonomic nervous system. Disruption of circadian rhythms will also promote cancer metastasis through direct and indirect actions of the nervous system. Therefore, the nervous system plays an important role in cancer metastasis.