A gene expression signature that defines breast cancer metastases

Copper metabolism-related signature for prognosis prediction and MMP13 served as malignant factor for breast cancer

Objectives

To comprehensively analyze the copper metabolism in Breast cancer, we established a prognostic signature for breast cancer (BC) related to copper metabolism.

Methods

Copper metabolism-related genes were sourced from previous literatures and were selected by the Univariate Cox regression. Cu-enrichment scores were calculated via ssGSEA. Differentially expressed genes were identified with limma between high and low Cu-enrichment scores group, then we used the Random Survival Forest and LASSO to build the CuScore for BC. Kaplan-Meier analysis, ROC curves, and Cox regression were used to evaluate CuScore. Genomic mutations were analyzed with GISTIC. Immune cells were examined using ESTIMATE, ssGSEA and TIMER. Enrichment analysis used clusterProfiler and GSVA. The GDSC database and oncoPredict package analyzed chemotherapeutic sensitivity. MMP13 was selected for in vitro assays.

Results

Four copper metabolism-related genes (UBE2D2, SLC31A1, ATP7A, and MAPK1) with prognostic value were identified. Higher expression levels of these genes were associated with higher Cu-enrichment scores, a factor of malignancy in breast cancer. Among 115 differentially expressed genes, 19 prognostic genes were identified, with three (CEACAM5, MMP13, and CRISP3) highlighted by Random Survival Forest and LASSO. Higher CuScores correlated with worse prognoses and were effective in predicting breast cancer outcomes. CuScore and metastasis were independent prognostic factors. Tumor-infiltrating immune cells were associated with lower CuScores. GO-GSEA analysis indicated six immune-related pathways might be regulated by CuScore. Patients with higher CuScores had lower TMB and were more sensitive to Sapitinib and LCL161, while those with lower CuScores might respond better to anti-PD1 therapy. High MMP13 expression in breast cancer was linked to malignancy, affecting cell proliferation and migration.

Conclusion

The identified copper metabolism-related gene signature has the potential to predict prognosis and guide clinical treatment for BC. Among these genes, MMP13 may act as a malignant factor in BC.

Cytokines-activated nuclear IKKα-FAT10 pathway induces breast cancer tamoxifen-resistance

Endocrine therapy that blocks estrogen signaling is the most effective treatment for patients with estrogen receptor positive (ER+) breast cancer. However, the efficacy of agents such as tamoxifen (Tam) is often compromised by the development of resistance. Here we report that cytokines-activated nuclear IKKα confers Tam resistance to ER+ breast cancer by inducing the expression of FAT10, and that the expression of FAT10 and nuclear IKKα in primary ER+ human breast cancer was correlated with lymphotoxin β (LTB) expression and significantly associated with relapse and metastasis in patients treated with adjuvant mono-Tam. IKKα activation or enforced FAT10 expression promotes Tam-resistance while loss of IKKα or FAT10 augments Tam sensitivity. The induction of FAT10 by IKKα is mediated by the transcription factor Pax5, and coordinated via an IKKα-p53-miR-23a circuit in which activation of IKKα attenuates p53-directed repression of FAT10. Thus, our findings establish IKKα-to-FAT10 pathway as a new therapeutic target for the treatment of Tam-resistant ER+ breast cancer.

Biosensors for metastatic cancer cell detection

Early detection and effective cancer treatment are critical to improving metastatic cancer cell diagnosis and management today. In particular, accurate qualitative diagnosis of metastatic cancer cell represents an important step in the diagnosis of cancer. Today, biosensors have been widely developed due to the daily need to measure different chemical and biological species. Biosensors are utilized to quantify chemical and biological phenomena by generating signals that are directly proportional to the quantity of the analyte present in the reaction. Biosensors are widely used in disease control, drug delivery, infection detection, detection of pathogenic microorganisms, and markers that indicate a specific disease in the body. These devices have been especially popular in the field of metastatic cancer cell diagnosis and treatment due to their portability, high sensitivity, high specificity, ease of use and short response time. This article examines biosensors for metastatic cancer cells. It also studies metastatic cancer cells and the mechanism of metastasis. Finally, the function of biosensors and biomarkers in metastatic cancer cells is investigated.

Collagen 1A1 (COL1A1) and Collagen11A1(COL11A1) as diagnostic biomarkers in Breast, colorectal and gastric cancers

PURPOSE

Collagen family genes (CFGs) play a significant role in the pathogenesis of cancers. This study aimed to evaluate changes in the expression levels (Els) of CFGs related to epithelial-mesenchymal transition (EMT) and metastasis in gastric (GC), breast (BC), and colorectal (CRC) cancers to introduce these genes as potential diagnostic biomarkers for these three types of cancer.

METHODS

The Cancer Genome Atlas (TCGA) examined ELS changes in CFGs associated with EMT and metastasis to determine their diagnostic value for GC, BC, and CRC. InteractiVenn was used to find genes shared by these three cancers. The biomarker role of CFGs was determined using the receiver operating characteristic (ROC) analysis. GC, BC, and CRC samples were analyzed using the RT-qPCR method to verify the bioinformatics results and evaluate the EL of the selected genes as biomarkers for these cancers.

RESULTS

The in-silico results showed a significant increase in the EL of several CFGs involved in EMT and metastasis in GC, BC, and CRC samples compared to healthy samples. Six common genes (COL11A1, COL12A1, COL1A1, COL1A2, COL5A1, and COL5A2) showed significantly increased in these three cancers, therebysupporting their oncogenic role. Furthermore, the biomarker-related analyses indicated that COL11A1 and COL1A1 were common diagnostic biomarkers for the three cancers. The RT-qPCR method confirmed that the ELs of COL11A1 and COL1A1 in the GC, BC, and CRC samples increased significantly compared to the adjacent normal samples.

CONCLUSION

CFGs in EMT and metastasis of GC, BC, and CRC are strong common diagnostic biomarkers for these cancers.

In silico analysis of differentially expressed-aberrantly methylated genes in breast cancer for prognostic and therapeutic targets

Breast cancer (BC) is the leading cause of death among women across the globe. Abnormal gene expression plays a crucial role in tumour progression, carcinogenesis and metastasis of BC. The alteration of gene expression may be through aberrant gene methylation. In the present study, differentially expressed genes which may be regulated by DNA methylation and their pathways associated with BC have been identified. Expression microarray datasets GSE10780, GSE10797, GSE21422, GSE42568, GSE61304, GSE61724 and one DNA methylation profile dataset GSE20713 were downloaded from Gene Expression Omnibus database (GEO). Differentially expressed-aberrantly methylated genes were identified using online Venn diagram tool. Based on fold change expression of differentially expressed-aberrantly methylated genes were chosen through heat map. Protein-protein interaction (PPI) network of the hub genes was constructed by Search Tool for the Retrieval of Interacting Genes (STRING). Gene expression and DNA methylation level of the hub genes were validated through UALCAN. Overall survival analysis of the hub genes was analysed through Kaplan-Meier plotter database for BC. A total of 72 upregulated-hypomethylated genes and 92 downregulated-hypermethylated genes were obtained from GSE10780, GSE10797, GSE21422, GSE42568, GSE61304, GSE61724, and GSE20713 datasets by GEO2R and Venn diagram tool. PPI network of the upregulated-hypomethylated hub genes (MRGBP, MANF, ARF3, HIST1H3D, GSK3B, HJURP, GPSM2, MATN3, KDELR2, CEP55, GSPT1, COL11A1, and COL1A1) and downregulated-hypermethylated hub genes were constructed (APOD, DMD, RBPMS, NR3C2, HOXA9, AMKY2, KCTD9, and EDN1). All the differentially expressed hub genes expression was validated in UALCAN database. 4 in 13 upregulated-hypomethylated and 5 in 8 downregulated-hypermethylated hub genes to be significantly hypomethylated or hypermethylated in BC were confirmed using UALCAN database (p < 0.05). MANF, HIST1H3D, HJURP, GSK3B, GPSM2, MATN3, KDELR2, CEP55, COL1A1, APOD, RBPMS, NR3C2, HOXA9, ANKMY2, and EDN1 were significantly (p < 0.05) associated with poor overall survival (OS). The identified aberrantly methylated-differentially expressed genes and their related pathways and function in BC can serve as novel diagnostic and prognostic biomarkers and therapeutic targets.Please confirm if the author names are presented accurately and in the correct sequence (given name, middle name/initial, family name). Author 4 Given name: [Jeewan Ram] Last name [Vishnoi]. Also, kindly confirm the details in the metadata are correct.It is correct.

Comprehensive analysis of the oncogenic and immunological role of SPON2 in human tumors

BACKGROUND

Sapiens spondin-2 (SPON2) is a protein found in the extracellular matrix that plays a role in a number of processes, including immune reactions and cell adhesion, and is closely linked to the emergence of a number of tumor types. However, we know very little about Sapiens spondin-2. Therefore, we performed a systematic pan-carcinogenic analysis to explore the relationship between Sapiens spondin-2 and cancers.

MATERIALS AND METHODS

By comprehensive use of datasets from TCGA, GEO, GTEx, HPA, CPTAC, GEPIA2, TIMER2, cBioPortal, STRING, we adopted bioinformatics methods to dig up the potential carcinogenesis of SPON2, including dissecting the correlation between SPON2 and gene expression, prognosis, gene mutation, Immunohistochemistry staining, immune cell infiltration, and constructed the interaction network of a total of 54 SPON2-binding proteins as well as explored the enrichment analysis of SPON2-related partners.

RESULTS

The expression of Sapiens spondin-2 in most tumor tissues was higher than that of normal tissues. In addition, SPON2 showed the early diagnostic value in 33 kinds of tumors and was positively or negatively associated with the prognosis of different tumors. It also validates that SPON2 is the gene associated with the majority of immune-infiltrating cells in pan-cancer. High SPON2 expression is associated with tumor progression related pathways.

CONCLUSION

We found and validated the potential use of SPON2 in cancer detection for the first time through pan-cancer analysis. The expression levels of SPON2 in various tumors were quite different from those in normal tissues. Furthermore, the performance of SPON2 in tumorigenesis and tumor immunity verified our hypothesis. At the same time, it has high specificity and sensitivity in cancer detection. Therefore, SPON2 can be employed as an auxiliary index for the initial diagnosis of tumors and a prognostic marker for various types of tumors.

Usefulness of COL11A1 as a Prognostic Marker of Tumor Infiltration

BACKGROUND

Determining the infiltration of carcinomas is essential for the proper follow-up and treatment of cancer patients. However, it continues to be a diagnostic challenge for pathologists in multiple types of tumors. In previous studies (carried out in surgical specimens), the protein COL11A1 has been postulated as an infiltration marker mainly expressed in the extracellular matrix (ECM). We hypothesized that a differential expression of COL11A1 may exist in the peritumoral stroma of tumors that have acquired infiltrating properties and that it may be detected in the small biopsies usually available in normal clinical practice.

MATERIAL AND METHODS

In our study, we performed immunohistochemical staining in more than 350 invasive and noninvasive small samples obtained via core needle biopsy (CNB), colonoscopy, or transurethral resection of bladder tumor (TURBT) of breast, colorectal, bladder, and ovarian cancer.

RESULTS

Our results revealed that COL11A1 immunostaining had a sensitivity to classify the samples into infiltrative vs. noninfiltrative tumors of 94% (breast), 97% (colorectal), >90% (bladder), and 74% (ovarian); and a specificity of 97% (breast), 100% (colorectal), and >90% (bladder). In ovarian cancer, the negative predictive value (0.59) did not present improvement over the usual histopathological markers. In all samples tested, the cumulative sensitivity was 86% and the specificity 96% (p < 0.0001).

CONCLUSIONS

COL11A1-positive immunostaining in small biopsies of breast, colon, bladder and ovarian cancer is an accurate predictive marker of tumor infiltration that can be easily implemented in daily clinical practice.

The Pleiotropy of PAX5 Gene Products and Function

PAX5, a member of the Paired Box (PAX) transcription factor family, is an essential factor for B-lineage identity during lymphoid differentiation. Mechanistically, PAX5 controls gene expression profiles, which are pivotal to cellular processes such as viability, proliferation, and differentiation. Given its crucial function in B-cell development, PAX5 aberrant expression also correlates with hallmark cancer processes leading to hematological and other types of cancer lesions. Despite the well-established association of PAX5 in the development, maintenance, and progression of cancer disease, the use of PAX5 as a cancer biomarker or therapeutic target has yet to be implemented. This may be partly due to the assortment of PAX5 expressed products, which layers the complexity of their function and role in various regulatory networks and biological processes. In this review, we provide an overview of the reported data describing PAX5 products, their regulation, and function in cellular processes, cellular biology, and neoplasm.

Comprehensive analysis of gene expression profiles to identify differential prognostic factors of primary and metastatic breast cancer

Breast cancer accounts for nearly half of all cancer-related deaths in women worldwide. However, the molecular mechanisms that lead to tumour development and progression remain poorly understood and there is a need to identify candidate genes associated with primary and metastatic breast cancer progression and prognosis. In this study, candidate genes associated with prognosis of primary and metastatic breast cancer were explored through a novel bioinformatics approach. Primary and metastatic breast cancer tissues and adjacent normal breast tissues were evaluated to identify biomarkers characteristic of primary and metastatic breast cancer. The Cancer Genome Atlas-breast invasive carcinoma (TCGA-BRCA) dataset (ID: HS-01619) was downloaded using the mRNASeq platform. Genevestigator 8.3.2 was used to analyse TCGA-BRCA gene expression profiles between the sample groups and identify the differentially-expressed genes (DEGs) in each group. For each group, Gene Ontology and Kyoto Encyclopedia of Genes and Genomes pathway enrichment analyses were used to determine the function of DEGs. Networks of protein-protein interactions were constructed to identify the top hub genes with the highest degree of interaction. Additionally, the top hub genes were validated based on overall survival and immunohistochemistry using The Human Protein Atlas. Of the top 20 hub genes identified, four (KRT14, KIT, RAD51, and TTK) were considered as prognostic risk factors based on overall survival. KRT14 and KIT expression levels were upregulated while those of RAD51 and TTK were downregulated in patients with breast cancer. The four proposed candidate hub genes might aid in further understanding the molecular changes that distinguish primary breast tumours from metastatic tumours as well as help in developing novel therapeutics. Furthermore, they may serve as effective prognostic risk markers based on the strong correlation between their expression and patient overall survival.

Distinct Gene Expression Profiles of Matched Primary and Metastatic Triple-Negative Breast Cancers

Background: Although triple-negative breast cancer (TNBC) is associated with an increased risk of recurrence and metastasis, the molecular mechanisms underlying metastasis in TNBC remain unknown. To identify transcriptional changes and genes regulating metastatic progression in TNBC, we compared the transcriptomic profiles of primary and matched metastatic tumors using massively parallel RNA sequencing. Methods: We performed gene expression profiling using formalin-fixed paraffin-embedded (FFPE) TNBC tissues of patients from two cohorts: the Zurich cohort (n = 31) and the Stavanger cohort (n = 5). Among the 31 patients in the Zurich cohort, 18 had primary TNBC tumors that did not metastasize, and 13 had primary tumors that metastasized (11 paired primary and locoregional recurrences). The Stavanger cohort included five matched primary and metastatic TNBC tumors. Significantly differentially expressed genes (DEGs; absolute fold change ≥2, p < 0.05) were identified and subjected to functional analyses. We investigated if there was any overlap between DEGs from both the cohorts with epithelial-to-mesenchymal-to-amoeboid transition (EMAT) gene signature. xCell was used to estimate relative fractions of 64 immune and stromal cell types in each RNA-seq sample. Results: In the Zurich cohort, we identified 1624 DEGs between primary TNBC tumors and matched metastatic lesions. xCell analysis revealed a significantly higher immune scores for metastatic lesions compared to paired primary tumors in the Zurich cohort. We also found significant upregulation of three MammaPrint signature genes (HRASLS, TGFB3 and RASSF7) in primary tumors that metastasized compared to primary tumors that remained metastasis-free. In the Stavanger cohort, we identified 818 DEGs between primary tumors and matched metastatic lesions. No significant differences in xCell immune scores were observed. We found that 21 and 14 DEGs from Zurich and Stavanger cohort, respectively, overlapped with the EMAT gene signature. In both cohorts, genes belonging to the MMP, FGF, and PDGFR families were upregulated in primary tumors compared to matched metastatic lesions. Conclusions: Our results suggest that distinct gene expression patterns exist between primary TNBCs and matched metastatic tumors. Further studies are warranted to explore whether these discrete expression profiles underlie or result from disease status.

Multiple roles for basement membrane proteins in cancer progression and EMT

Metastasis or the progression of malignancy poses a major challenge in cancer therapy and is the principal reason for increased mortality. The epithelial-Mesenchymal transition (EMT) of the Basement Membrane (BM) allows cells of epithelial phenotype to transform into a mesenchymal-like (quasi-mesenchymal) phenotype and metastasize via the lymphovascular system through a metastatic cascade by intravasation and extravasation. This helps in the progression of carcinoma from the primary site to distant organs. Collagen, laminin, and integrin are the prime components of BM and help in tumor cell metastasis, which makes them ideal cancer drug targets. Further, recent studies have shown that collagen, laminin, and integrin can be used as a biomarker for metastatic cells. In this review, we have summarized the current knowledge of such therapeutics, which are either currently in preclinical or clinical stages and could be promising cancer therapeutics. DATA AVAILABILITY: Not applicable.

Collagen XI alpha 1 chain, a potential therapeutic target for cancer

Extracellular matrix (ECM) plays an important role in the progression of cancer. Collagen is the most abundant component in ECM, and it is involved in the biological formation of cancer. Although type XI collagen is a minor fibrillar collagen, collagen XI alpha 1 chain (COL11A1) has been found to be upregulated in a variety of cancers including ovarian cancer, breast cancer, thyroid cancer, pancreatic cancer, non-small-cell lung cancer, and transitional cell carcinoma of the bladder. High levels of COL11A1 usually predict poor prognosis, while COL11A1 is related to angiogenesis, invasion, and drug resistance of cancer. However, little is known about the specific mechanism by which COL11A1 regulates tumor progression. Here, we have organized and summarized the recent developments regarding elucidation of the relationship between COL11A1 and various cancers, as well as the interaction between COL11A1 and intracellular signaling pathways. In addition, we have selected therapeutic agents targeting COL11A1. All these indicate the possibility of using COL11A1 as a target for cancer treatment.

MiR-144-3p inhibits the proliferation, migration and invasion of lung adenocargen cancer cells by targeting COL11A1

This study aimed to investigate the regulatory relationship between miR-144-3p and COL11A1, and to explore its effect on the proliferation, migration and invasion of lung adenocarcinoma (LUAD) cells. A series of methods and experiments were applied. miR-144-3p was downregulated in LUAD tissue and cells, whereas COL11A1 was highly expressed. Overexpressing miR-144-3p inhibited the proliferation, migration and invasion of LUAD cells, which could be reversed by overexpression of COL11A1. Overexpressing miR-144-3p inhibits the proliferation, migration and invasion of LUAD cells by silencing COL11A1.

Targeting PVT1 Exon 9 Re-Expresses Claudin 4 Protein and Inhibits Migration by Claudin-Low Triple Negative Breast Cancer Cells

PVT1 is a long non-coding RNA transcribed from a gene located at the 8q24 chromosomal region that has been implicated in multiple cancers including breast cancer (BC). Amplification of the 8q24 chromosomal region is a common event in BC and is associated with poor clinical outcomes. Claudin-low (CL) triple negative breast cancer (TNBC) is a subtype of BC with a particularly dismal outcome. We assessed PVT1 exon 9 expression in the T47D estrogen receptor positive BC cell line, and in the MDA MB 468 and MDA MB 231 TNBC cell lines, followed by the assessment of the expression of claudins 1, 3, 4 and 7, in MDA MB 468 and MDA MB 231 (TNBC) cells. We found that MDA MB 231 TNBC cells significantly express less claudin 1, 3, 4, and 7 than MDA MB 468 TNBC cells. PVT1 exon 9 is significantly upregulated in MDA MB 231 CL TNBC cells, and significantly downregulated in MDA MB 468 claudin high (CH) TNBC cells, in comparison to T47D estrogen receptor positive BC cells. We then analyzed the functional consequences of siRNA targeting of PVT1 exon 9 expression in the MDA MB 231 CL TNBC cells. Notably, siRNA targeting of PVT1 exon 9 expression in the MDA MB 231 CL TNBC cells led to a significant reduction in migration and the re-expression of claudin 4. Taken together, our data indicate that PVT1 exon 9 regulates claudin 4 expression and migration in CL TNBC cells, and may have clinical implications in CL TNBC.

Novel somatic alterations underlie Chinese papillary thyroid carcinoma

To characterize the somatic alterations of papillary thyroid carcinomas (PTC) in Chinese patients, we performed the next-generation-sequencing (NGS) study of the tumor-normal pairs of DNA and RNA samples extracted from 16 Chinese PTC patients. The whole genome sequencing (WGS) and transcriptome sequencing (RNA-seq) were conducted for 6 patients who were either current or former smokers and the whole exome sequencing (WES) and RNA-seq were conducted for another 10 patients who were never smokers. The NGS data were analyzed to identify somatic alteration events that may underlie PTC in Chinese patients. We identified a number of PTC driver genes harboring somatic driver mutations with significant functional impact such as COL11A1, TP53, PLXNA4, UBA1, AHNAK, CSMD2 and TTLL5 etc. Significant driver pathways underlying PTC were found, namely, the metabolic pathway, the pathway in cancer, the olfactory transduction pathway and the calcium signaling pathway. In addition, this study revealed genes with significant somatic copy number aberrations and corresponding somatic gene expression changes in PTC tumors, the most promising ones being BRD9, TRIP13, FZD3, and TFDP1 etc. We also identified several structural variants of PTCs, especially the novel in-frame fusion proteins such as TRNAU1AP-RCC1, RAB3GAP1-R3HDM1, and ENAH-ZSWIM5. Our study provided a list of novel PTC candidate genes with somatic alterations that may function as biomarkers for PTC in Chinese patients. The follow-up mechanism studies may be conducted based on the findings from this study.

Wnt signaling in breast cancer: biological mechanisms, challenges and opportunities

Wnt signaling is a highly conserved signaling pathway that plays a critical role in controlling embryonic and organ development, as well as cancer progression. Genome-wide sequencing and gene expression profile analyses have demonstrated that Wnt signaling is involved mainly in the processes of breast cancer proliferation and metastasis. The most recent studies have indicated that Wnt signaling is also crucial in breast cancer immune microenvironment regulation, stemness maintenance, therapeutic resistance, phenotype shaping, etc. Wnt/β-Catenin, Wnt-planar cell polarity (PCP), and Wnt-Ca2+ signaling are three well-established Wnt signaling pathways that share overlapping components and play different roles in breast cancer progression. In this review, we summarize the main findings concerning the relationship between Wnt signaling and breast cancer and provide an overview of existing mechanisms, challenges, and potential opportunities for advancing the therapy and diagnosis of breast cancer.

Studying Lymphatic Metastasis in Breast Cancer: Current Models, Strategies, and Clinical Perspectives

Breast cancer is the most commonly diagnosed cancer in women and the second most common cause of cancer-related deaths in the United States. Although early detection has significantly decreased breast cancer mortality, patients diagnosed with distant metastasis still have a very poor prognosis. The most common site that breast cancer spreads to are local lymph nodes. Therefore, the presence of lymph node metastasis remains one of most important prognostic factors in breast cancer patients. Given its significant clinical implications, increased efforts have been dedicated to better understand the molecular mechanism governing lymph node metastasis in breast cancer. The identification of lymphatic-specific biomarkers, including podoplanin and LYVE-1, has propelled the field of lymphatic metastasis forward. In addition, several animal models such as cell line-derived xenografts, patient-derived xenografts, and spontaneous tumor models have been developed to recreate the process of lymphatic metastasis. Moreover, the incorporation of various -omic platforms have provided further insight into the genetic drivers facilitating lymphatic metastasis, as well as potential biomarkers and therapeutic targets. Here, we highlight various models of lymphatic metastasis, their potential pitfalls, and other tools available to study lymphatic metastasis including imaging modalities and -omic studies.

Prediction of tumor metastasis from sequencing data in the era of genome sequencing

Tumor metastasis is the key reason for the high mortality rate of tumor. Growing number of scholars have begun to pay attention to the research on tumor metastasis and have achieved satisfactory results in this field. The advent of the era of sequencing has enabled us to study cancer metastasis at the molecular level, which is essential for understanding the molecular mechanism of metastasis, identifying diagnostic markers and therapeutic targets and guiding clinical decision-making. We reviewed the metastasis-related studies using sequencing data, covering detection of metastasis origin sites, determination of metastasis potential and identification of distal metastasis sites. These findings include the discovery of relevant markers and the presentation of prediction tools. Finally, we discussed the challenge of studying metastasis considering the difficulty of obtaining metastatic cancer data, the complexity of tumor heterogeneity and the uncertainty of sample labels.

Differentially Expressed Genes in Matched Normal, Cancer, and Lymph Node Metastases Predict Clinical Outcomes in Patients With Breast Cancer

Genome-wide screening of transcriptional changes among normal, cancer, and nodal metastases provides insights into the molecular basis of breast cancer (BC) progression and metastasis. To identify transcriptional changes and differentially expressed genes (DEGs) in the metastatic progression of BC and to determine the prognostic role of these DEGs in clinical outcome, we compared transcriptome profiling in matched normal, cancer, and lymph node metastatic tissues of 7 patients with estrogen receptor-positive, HER2-negative BC by using massive parallel RNA sequencing. The global profiles of gene expression in cancer and nodal metastases were highly correlated (r=0.962, P<0.001). In 6 (85.8%) patients, cancer and corresponding nodal metastases from the same patient clustered together. We identified 1522 and 664 DEGs between normal and cancer and between cancer and nodal metastases, respectively. The DEGs in normal versus cancer and cancer versus nodal metastases were significantly clustered in 1 and 8 Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways, respectively. The chemokine signaling pathway was the most significant pathway in the cancer-to-nodal metastasis transition (false discovery rate=2.15E-13). The expression of 2 dysregulated RAC2 and PTGDS genes was confirmed by quantitative real-time polymerase chain reaction and immunohistochemistry. Interestingly, the lower RAC2 and PTGDS expression were associated with significantly worse disease-free survival in patients with BC. Our results show a high concordance of gene expression in BC and their nodal metastases, and identify DEGs associated with the metastatic progression of BC. The DEGs identified in this study represent novel biomarkers for predicting the prognosis of patients with BC.

Comprehensive Profiling of Poor-Risk Paired Primary and Recurrent Triple-Negative Breast Cancers Reveals Immune Phenotype Shifts

PURPOSE

Emerging data suggest immune checkpoint inhibitors have reduced efficacy in heavily pretreated triple-negative breast cancers (TNBC), but underlying mechanisms are poorly understood. To better understand the phenotypic evolution of TNBCs, we studied the genomic and transcriptomic profiles of paired tumors from patients with TNBC.

EXPERIMENTAL DESIGN

We collected paired primary and metastatic TNBC specimens from 43 patients and performed targeted exome sequencing and whole-transcriptome sequencing. From these efforts, we ascertained somatic mutation profiles, tumor mutational burden (TMB), TNBC molecular subtypes, and immune-related gene expression patterns. Stromal tumor-infiltrating lymphocytes (stromal TIL), recurrence-free survival, and overall survival were also analyzed.

RESULTS

We observed a typical TNBC mutational landscape with minimal shifts in copy number or TMB over time. However, there were notable TNBC molecular subtype shifts, including increases in the Lehmann/Pietenpol-defined basal-like 1 (BL1, 11.4%-22.6%) and mesenchymal (M, 11.4%-22.6%) phenotypes, and a decrease in the immunomodulatory phenotype (IM, 31.4%-3.2%). The Burstein-defined basal-like immune-activated phenotype was also decreased (BLIA, 42.2%-17.2%). Among downregulated genes from metastases, we saw enrichment of immune-related Kyoto Encyclopedia of Genes and Genomes pathways and gene ontology (GO) terms, and decreased expression of immunomodulatory gene signatures (P < 0.03) and percent stromal TILs (P = 0.03). There was no clear association between stromal TILs and survival.

CONCLUSIONS

We observed few mutational shifts, but largely consistent transcriptomic shifts in longitudinally paired TNBCs. Transcriptomic and IHC analyses revealed significantly reduced immune-activating gene expression signatures and TILs in recurrent TNBCs. These data may explain the observed lack of efficacy of immunotherapeutic agents in heavily pretreated TNBCs. Further studies are ongoing to better understand these initial observations.See related commentary by Savas and Loi, p. 526.

Aberrant WNT/CTNNB1 Signaling as a Therapeutic Target in Human Breast Cancer: Weighing the Evidence

WNT signaling is crucial for tissue morphogenesis during development in all multicellular animals. After birth, WNT/CTNNB1 responsive stem cells are responsible for tissue homeostasis in various organs and hyperactive WNT/CTNNB1 signaling is observed in many different human cancers. The first link between WNT signaling and breast cancer was established almost 40 years ago, when Wnt1 was identified as a proto-oncogene capable of driving mammary tumor formation in mice. Since that discovery, there has been a dedicated search for aberrant WNT signaling in human breast cancer. However, much debate and controversy persist regarding the importance of WNT signaling for the initiation, progression or maintenance of different breast cancer subtypes. As the first drugs designed to block functional WNT signaling have entered clinical trials, many questions about the role of aberrant WNT signaling in human breast cancer remain. Here, we discuss three major research gaps in this area. First, we still lack a basic understanding of the function of WNT signaling in normal human breast development and physiology. Second, the overall extent and precise effect of (epi)genetic changes affecting the WNT pathway in different breast cancer subtypes are still unknown. Which underlying molecular and cell biological mechanisms are disrupted as a result also awaits further scrutiny. Third, we survey the current status of targeted therapeutics that are aimed at interfering with the WNT pathway in breast cancer patients and highlight the importance and complexity of selecting the subset of patients that may benefit from treatment.

Computational Detection of Breast Cancer Invasiveness with DNA Methylation Biomarkers

Breast cancer is the most common female malignancy. It has high mortality, primarily due to metastasis and recurrence. Patients with invasive and noninvasive breast cancer require different treatments, so there is an urgent need for predictive tools to guide clinical decision making and avoid overtreatment of noninvasive breast cancer and undertreatment of invasive cases. Here, we divided the sample set based on the genome-wide methylation distance to make full use of metastatic cancer data. Specifically, we implemented two differential methylation analysis methods to identify specific CpG sites. After effective dimensionality reduction, we constructed a methylation-based classifier using the Random Forest algorithm to categorize the primary breast cancer. We took advantage of breast cancer (BRCA) HM450 DNA methylation data and accompanying clinical data from The Cancer Genome Atlas (TCGA) database to validate the performance of the classifier. Overall, this study demonstrates DNA methylation as a potential biomarker to predict breast tumor invasiveness and as a possible parameter that could be included in the studies aiming to predict breast cancer aggressiveness. However, more comparative studies are needed to assess its usability in the clinic. Towards this, we developed a website based on these algorithms to facilitate its use in studies and predictions of breast cancer invasiveness.

Collagen (XI) alpha-1 chain is an independent prognostic factor in breast ductal carcinoma in situ

Collagen11A1 (COL11A1) is a fibrillary type collagen constituting a minor component of the extracellular matrix and plays role in tissue tensile strength. Overexpression of COL11A1 expression is associated with aggressive behavior and poor outcome in several human malignancies. In this study, we evaluated the association between COL11A1 expression and clinicopathological parameters of the breast ductal carcinoma in situ (DCIS) and its prognostic value. COL11A1 protein expression was assessed immunohistochemically in a large well-characterized cohort of DCIS including pure (n = 776) and DCIS associated with invasive carcinoma (DCIS-mixed, n = 239). COL11A1 expression was assessed in tumor cells and surrounding stromal cells, and correlated with clinicopathological parameters, immunoprofile and disease outcome. In pure DCIS, high COL11A1 expression was observed in tumor cells and surrounding stromal cells in 25 and 13% of cases, respectively. Higher COL11A1 expression within the stromal cells was associated with hormone receptor negative, HER2 enriched and triple negative molecular subtypes and showed a positive linear correlation with proliferation index, dense tumor infiltrating lymphocytes and hypoxia-inducible factor 1 alpha. COL11A1 expression in tumor and stromal cells was significantly higher in DCIS associated with invasive carcinoma than in pure DCIS, and within the DCIS-mixed cohort, the invasive component showed higher COL11A1 expression than the DCIS component (all, p < 0.0001). Overexpression of stromal COL11A1 was an independent predictor of shorter local recurrence-free interval for all recurrences (HR = 13.2, 95% CI = 6.9-25.4, p < 0.0001) and for invasive recurrences (HR = 11.2, 95% CI = 4.9-25.8, p < 0.0001). When incorporated with other risk factors, stromal COL11A1 provided better patient risk stratification. DCIS with higher stromal COL11A1 expression showed poor outcome even with adjuvant radiotherapy management. In conclusion, overexpression of stromal COL11A1 is associated with invasive recurrence in DCIS and is a potential marker to predict the response to radiotherapy.

Novel role of the SRY-related high-mobility-group box D gene in cancer

The SRY-related high-mobility-group box (Sox) gene family encodes a set of transcription factors and is defined by the presence of highly conserved domains. The Sox gene can be divided into 10 groups (A-J). The SoxD subpopulation consists of Sox5, Sox6, Sox13 and Sox23, which are involved in the transcriptional regulation of developmental processes, including embryonic development, nerve growth and cartilage formation. Recently, the SoxD gene family was recognized as important transcriptional regulators associated with many types of cancer. In addition, Sox5 and Sox6 are representatives of the D subfamily, and there are many related studies; however, there are few reports on Sox13 and Sox23. In this review, we first introduce the structures of the SoxD genes. Next, we summarize the latest research progress on SoxD in various types of cancer. Finally, we discuss the potential direction of future SoxD research. In general, the information reviewed here may contribute to future experimental design and increase the potential of SoxD as a cancer treatment target.

Genomic comparison of paired primary breast carcinomas and lymph node macrometastases using the Oncotype DX Breast Recurrence Score® test

PURPOSE

Adjuvant therapy decisions may in part be based on results of Oncotype DX Breast Recurrence Score^® (RS) testing of primary tumors. When necessary, lymph node metastases may be considered as a surrogate. Here we evaluate the concordance in gene expression between primary breast cancers and synchronous lymph node metastases, based on results from quantitative RT-PCR-based RS testing between matched primary tumors and synchronous nodal metastases.

METHODS

This retrospective, exploratory study included patients (≥ 18 years old) treated at our center (2005-2009) who had ER+ , HER2-negative invasive breast cancer and synchronous nodal metastases with available tumor blocks from both sites. Paired tissue blocks underwent RS testing, and RS and single-gene results for ER, PR, and HER2 were explored between paired samples.

RESULTS

A wide distribution of RS results in tumors and in synchronous nodal metastases were modestly correlated between 84 paired samples analyzed (Pearson correlation 0.69 [95% CI 0.55-0.78]). Overall concordance in RS group classification between samples was 63%. ER, PR, and HER2 by RT-PCR between the primary tumor and lymph node were also modestly correlated (Pearson correlation [95% CI] 0.64 [0.50-0.75], 0.64 [0.49-0.75], and 0.51 [0.33-0.65], respectively). Categorical concordance (positive or negative) was 100% for ER, 77% for PR, and 100% for HER2.

CONCLUSIONS

There is modest correlation in continuous gene expression, as measured by the RS and single-gene results for ER, PR, and HER2 between paired primary tumors and synchronous nodal metastases. RS testing for ER+ breast cancer should continue to be based on analysis of primary tumors.

Differential gene expression and AKT targeting in triple negative breast cancer

Background: Metastatic triple negative breast cancer (mTNBC) is a heterogeneous disease with poor prognosis. Molecular evolution of TNBC through chemotherapy selection pressure is well recognized but poorly understood. PI3K/AKT/mTOR is one of the most commonly identified oncogenic-driver pathways in breast cancer. The current study is designed to understand the genomic and transcriptomic changes, focusing on the PI3K/AKT/mTOR pathway alterations in paired primary and metastatic TNBCs. Results: Genomic analysis of 7 paired specimens identified 67 known mutations including those from the following signaling pathways: cell cycle, p53, PI3K/AKT/mTOR, RAS/MAPK, and RTK/GF. Principle coordinate analysis (PCoA) identified 4 distinctive molecular groups based on the gene expression patterns of PI3K/AKT/mTOR pathway. Key differentially-expressed genes included AKT3, GSK3B, GNA11, PI3KR1, and GNAQ. Importantly, AKT-targeted therapy showed efficacy in a patient-derived xenograft (PDX) model of TNBC in vivo. Conclusion: Genomic discordance of paired primary and metastatic TNBCs was identified, with significant increase in tumor proliferation pathways seen in metastases. Among the differentially expressed genes, AKT3 can potentially serve as a target for novel combination therapy for treatment of metastatic TNBC. Methods: Paired specimens from 10 patients with TNBCs were identified through an IRB-approved protocol (2002-2015). FoundationOneTM sequencing was performed for genomic profiling, and Affymetrix Human Genechip 2.0st was used for mRNA expression profiling. The similarity among samples was calculated based on Pearson correlation coefficients, which were used to construct hierarchical clustering and heat maps.

Collagens and Cancer associated fibroblasts in the reactive stroma and its relation to Cancer biology

The extracellular matrix (ECM) plays an important role in cancer progression. It can be divided into the basement membrane (BM) that supports epithelial/endothelial cell behavior and the interstitial matrix (IM) that supports the underlying stromal compartment. The major components of the ECM are the collagens. While breaching of the BM and turnover of e.g. type IV collagen, is a well described part of tumorigenesis, less is known regarding the impact on tumorigenesis from the collagens residing in the stroma. Here we give an introduction and overview to the link between tumorigenesis and stromal collagens, with focus on the fibrillar collagens type I, II, III, V, XI, XXIV and XXVII as well as type VI collagen. Moreover, we discuss the impact of the cells responsible for this altered stromal collagen remodeling, the cancer associated fibroblasts (CAFs), and how these cells are key players in orchestrating the tumor microenvironment composition and tissue microarchitecture, hence also driving tumorigenesis and affecting response to treatment. Lastly, we discuss how specific collagen-derived biomarkers reflecting the turnover of stromal collagens and CAF activity may be used as tools to non-invasively interrogate stromal reactivity in the tumor microenvironment and predict response to treatment.

FUS/circ_002136/miR-138-5p/SOX13 feedback loop regulates angiogenesis in Glioma

Background

Angiogenesis plays a critical role in the progression of glioma. Previous studies have indicated that RNA-binding proteins (RBPs) interact with RNAs and participate in the regulation of the malignant behaviors of tumors. As a type of endogenous non-coding RNAs, circular RNAs (circRNAs) are abnormally expressed in various cancers and are involved in diverse tumorigeneses including angiogenesis.

Methods

The expression levels of FUS, circ_002136, miR-138-5p, SOX13, and SPON2 were determined using quantitative real-time PCR (qRT-PCR) and western blot. Transient cell transfection was performed using the Lipofectamine 3000 reagent. The RNA-binding protein immunoprecipitation (RNA-IP) and the RNA pull-down assays were used to detect the interaction between FUS and circ_002136. The dual-luciferase reporter assay system was performed to detect the binding sites of circ_002136 and miR-138-5p, miR-138-5p and SOX13. The chromatin immunoprecipitation (ChIP) assays were used to examine the interactions between transcription factor SOX13 and its target proteins .

Results

We demonstrated that down-regulation of FUS or circ_002136 dramatically inhibited the viability, migration and tube formation of U87 glioma-exposed endothelial cells (GECs). MiR-138-5p was down-regulated in GECs and circ_002136 functionally targeted miR-138-5p in an RNA-induced silencing complex (RISC). Inhibition of circ_002136, combined with the restoration of miR-138-5p, robustly reduced the angiogenesis of GECs. As a target gene of miR-138-5p, SOX13 was overexpressed in GECs and was proved to be involved in circ_002136 and miR-138-5p-mediated angiogenesis in gliomas. In addition, we found that SOX13 was directly associated with and activated the SPON2 promoter, thereby up-regulating the expression of SPON2 at the transcriptional level. Knockdown of SPON2 suppressed the angiogenesis in GECs. More important, SOX13 activated the FUS promoter and increased its expression, forming a feedback loop.

Conclusion

Our data suggests that the feedback loop of FUS/circ_002136/miR-138-5p/SOX13 played a crucial role in the regulation of angiogenesis in glioma. This also provides a potential target and an alternative strategy for combined glioma therapy.

Electronic supplementary material

The online version of this article (10.1186/s13046-019-1065-7) contains supplementary material, which is available to authorized users.

CCNE1 amplification is associated with poor prognosis in patients with triple negative breast cancer

BACKGROUND

Triple negative breast cancer (TNBC) is aggressive with limited treatment options upon recurrence. Molecular discordance between primary and metastatic TNBC has been observed, but the degree of biological heterogeneity has not been fully explored. Furthermore, genomic evolution through treatment is poorly understood. In this study, we aim to characterize the genomic changes between paired primary and metastatic TNBCs through transcriptomic and genomic profiling, and to identify genomic alterations which may contribute to chemotherapy resistance.

METHODS

Genomic alterations and mRNA expression of 10 paired primary and metastatic TNBCs were determined through targeted sequencing, microarray analysis, and RNA sequencing. Commonly mutated genes, as well as differentially expressed and co-expressed genes were identified. We further explored the clinical relevance of differentially expressed genes between primary and metastatic tumors to patient survival using large public datasets.

RESULTS

Through gene expression profiling, we observed a shift in TNBC subtype classifications between primary and metastatic TNBCs. A panel of eight cancer driver genes (CCNE1, TPX2, ELF3, FANCL, JAK2, GSK3B, CEP76, and SYK) were differentially expressed in recurrent TNBCs, and were also overexpressed in TCGA and METABRIC. CCNE1 and TPX2 were co-overexpressed in TNBCs. DNA mutation profiling showed that multiple mutations occurred in genes comprising a number of potentially targetable pathways including PI3K/AKT/mTOR, RAS/MAPK, cell cycle, and growth factor receptor signaling, reaffirming the wide heterogeneity of mechanisms driving TNBC. CCNE1 amplification was associated with poor overall survival in patients with metastatic TNBC.

CONCLUSIONS

CCNE1 amplification may confer resistance to chemotherapy and is associated with poor overall survival in TNBC.

Spondin-2 is a novel diagnostic biomarker for laryngeal squamous cell carcinoma

Spondin-2, belongs to the SOX (SRY-related HMG box) gene family, plays a vital role in the development of malignancy, however, the role of Spondin-2 in laryngeal squamous cell carcinoma (LSCC) remains unknown. The aim of this study is to investigate the prognostic significance of and probable mechanism of Spondin-2 in LSCC. qRT-PCR, western blotting assays and IHC analysis demonstrated that Spondin-2 was significantly increased in LSCC tissues compared with adjacent non-tumorous tissues. In addition, high levels of Spondin-2 was associated with clinical stage, lymph node metastasis and pathology grade of LSCC patients (P ＜0.05). Kaplan-Meier analysis showed that patients with high expression of Spondin-2 had a lower overall survival rate (P＜0.05) than that with low expression of Spondin-2. Moreover, spondin-2 silencing inhibited the proliferation of LSCC cells through inhibiting the activation of PI3K/AKT signaling. In conclusion, spondin-2 might be a novel therapeutic target and prognostic biomarker for LSCC patients.

Molecular patterns of cancer colonisation in lymph nodes of breast cancer patients

Lymph node (LN) metastasis is an important prognostic parameter in breast carcinoma, a crucial site for tumour–immune cell interaction and a gateway for further dissemination of tumour cells to other metastatic sites. To gain insight into the underlying molecular changes from the pre-metastatic, via initial colonisation to the fully involved LN, we reviewed transcriptional research along the evolving microenvironment of LNs in human breast cancers patients. Gene expression studies were compiled and subjected to pathway-based analyses, with an emphasis on immune cell-related genes. Of 366 studies, 14 performed genome-wide gene expression comparisons and were divided into six clinical-biological scenarios capturing different stages of the metastatic pathway in the LN, as follows: metastatically involved LNs are compared to their patient-matched primary breast carcinomas (scenario 1) or the normal breast tissue (scenario 2). In scenario 3, uninvolved LNs were compared between LN-positive patients and LN-negative patients. Scenario 4 homed in on the residual uninvolved portion of involved LNs and compared it to the patient-matched uninvolved LNs. Scenario 5 contrasted uninvolved and involved LNs, whilst in scenario 6 involved (sentinel) LNs were assessed between patients with other either positive or negative LNs (non-sentinel).

Gene lists from these chronological steps of LN metastasis indicated that gene patterns reflecting deficiencies in dendritic cells and hyper-proliferation of B cells parallel to tumour promoting pathways, including cell adhesion, extracellular matrix remodelling, cell motility and DNA repair, play key roles in the changing microenvironment of a pro-metastatic to a metastatically involved LN. Similarities between uninvolved LNs and the residual uninvolved portion of involved LNs hinted that LN alterations expose systemic tumour-related immune responses in breast cancer patients. Despite the diverse settings, gene expression patterns at different stages of metastatic colonisation in LNs were recognised and may provide potential avenues for clinical interventions to counteract disease progression for breast cancer patients.

Insulin-like growth factor receptor signaling in breast tumor epithelium protects cells from endoplasmic reticulum stress and regulates the tumor microenvironment

BACKGROUND

Early analyses of human breast cancer identified high expression of the insulin-like growth factor type 1 receptor (IGF-1R) correlated with hormone receptor positive breast cancer and associated with a favorable prognosis, whereas low expression of IGF-1R correlated with triple negative breast cancer (TNBC). We previously demonstrated that the IGF-1R acts as a tumor and metastasis suppressor in the Wnt1 mouse model of TNBC. The mechanisms for how reduced IGF-1R contributes to TNBC phenotypes is unknown.

METHODS

We analyzed the METABRIC dataset to further stratify IGF-1R expression with patient survival and specific parameters of TNBC. To investigate molecular events associated with the loss of IGF-1R function in breast tumor cells, we inhibited IGF-1R in human cell lines using an IGF-1R blocking antibody and analyzed MMTV-Wnt1-mediated mouse tumors with reduced IGF-1R function through expression of a dominant-negative transgene.

RESULTS

Our analysis of the Molecular Taxonomy of Breast Cancer International Consortium (METABRIC) dataset revealed association between low IGF-1R and reduced overall patient survival. IGF-1R expression was inversely correlated with patient survival even within hormone receptor-positive breast cancers, indicating reduced overall patient survival with low IGF-1R was not due simply to low IGF-1R expression within TNBCs. Inhibiting IGF-1R in either mouse or human tumor epithelial cells increased reactive oxygen species (ROS) production and activation of the endoplasmic reticulum stress response. IGF-1R inhibition in tumor epithelial cells elevated interleukin (IL)-6 and C-C motif chemokine ligand 2 (CCL2) expression, which was reversed by ROS scavenging. Moreover, the Wnt1/dnIGF-1R primary tumors displayed a tumor-promoting immune phenotype. The increased CCL2 promoted an influx of CD11b⁺ monocytes into the primary tumor that also had increased matrix metalloproteinase (MMP)-2, MMP-3, and MMP-9 expression. Increased MMP activity in the tumor stroma was associated with enhanced matrix remodeling and collagen deposition. Further analysis of the METABRIC dataset revealed an increase in IL-6, CCL2, and MMP-9 expression in patients with low IGF-1R, consistent with our mouse tumor model and data in human breast cancer cell lines.

CONCLUSIONS

Our data support the hypothesis that reduction of IGF-1R function increases cellular stress and cytokine production to promote an aggressive tumor microenvironment through infiltration of immune cells and matrix remodeling.

Upregulation of Spondin-2 protein expression correlates with poor prognosis in hepatocellular carcinoma

OBJECTIVE

The aim of this study was to measure the extracellular matrix protein Spondin-2 (SPON2) in hepatocellular carcinoma (HCC) tissues and to determine its potential value as a prognostic indicator by assessing its correlation with clinicopathological variables and survival.

METHODS

SPON2 mRNA expression was assessed in 20 matched pairs of HCC and non-cancerous liver tissues by quantitative reverse transcription-polymerase chain reaction analysis. SPON2 protein expression was determined in 107 matched pairs of HCC and normal liver tissue by immunohistochemical staining of tissue microarrays.

RESULTS

Analysis of patient tissues and Oncomine datasets showed that SPON2 mRNA and SPON2 protein expression were both significantly upregulated in HCC tissues, compared with non-cancerous liver tissue; moreover, both correlated significantly with tumor size. Kaplan-Meier analysis revealed that HCC patients who showed high levels of cytoplasmic SPON2 protein had poorer survival following curative resection, compared with HCC patients who exhibited low protein expression levels. Multivariate Cox regression analysis showed that tumor thrombus and SPON2 protein expression both independently correlated with reduced survival in HCC patients.

CONCLUSION

Upregulated expression of SPON2 protein in tumor tissue could be an effective prognostic indicator for patients with HCC.

Quantifying critical states of complex diseases using single-sample dynamic network biomarkers

Dynamic network biomarkers (DNB) can identify the critical state or tipping point of a disease, thereby predicting rather than diagnosing the disease. However, it is difficult to apply the DNB theory to clinical practice because evaluating DNB at the critical state required the data of multiple samples on each individual, which are generally not available, and thus limit the applicability of DNB. In this study, we developed a novel method, i.e., single-sample DNB (sDNB), to detect early-warning signals or critical states of diseases in individual patients with only a single sample for each patient, thus opening a new way to predict diseases in a personalized way. In contrast to the information of differential expressions used in traditional biomarkers to “diagnose disease”, sDNB is based on the information of differential associations, thereby having the ability to “predict disease” or “diagnose near-future disease”. Applying this method to datasets for influenza virus infection and cancer metastasis led to accurate identification of the critical states or correct prediction of the immediate diseases based on individual samples. We successfully identified the critical states or tipping points just before the appearance of disease symptoms for influenza virus infection and the onset of distant metastasis for individual patients with cancer, thereby demonstrating the effectiveness and efficiency of our method for quantifying critical states at the single-sample level.

The concept of dynamic network biomarkers (DNB) was proposed for detecting the critical state or tipping point of a complex disease (a pre-disease state immediately preceding the disease state), and has been applied to study the mechanism of cell fate decision and immune checkpoint blockade. But DNB cannot be used to identify the critical state or tipping point for a single patient because evaluating DNB for critical state required the data of multiple samples. The proposed method can identify the critical state of a complex disease for a single patient by implementing the concept of DNB. This method not only can be applied to detect the critical state or tipping point of a single sample, but also can be used to study the mechanism of complex disease at a single sample level. The ability of accurately and efficiently identifying the critical state for a single sample can benefit the development of personalized medicine.

P-cadherin: a useful biomarker for axillary-based breast cancer decisions in the clinical practice

Axillary lymph node metastases represent the most powerful breast cancer prognostic factor, dictating disease staging and clinical therapeutic decisions. Nonetheless, breast cancer patients with positive lymph nodes still exhibit a heterogeneous behavior regarding disease progression. Stem-like subpopulations of cancer cells show high migratory and metastatic capacity, thus we hypothesize that breast cancer stem cell markers evaluation in metastasized lymph nodes could provide a more accurate prediction of patient's prognosis. Therefore, the expression profile of P-cadherin, CD44, and CD49f, which have been already associated to stem cell properties in breast cancer, has been evaluated by immunohistochemistry in a series of 135 primary tumors and matched axillary lymph node metastases from 135 breast cancer patients. Taking in consideration the expression of the stem cell markers only in axillary nodes, P-cadherin was the only biomarker significantly associated with poor disease-free and overall patient's survival. Moreover, although a concordant expression between primary tumors and matched lymph nodes has been found in the majority of the cases, a small but significant percentage displayed divergent expression (18.2-26.2%). Remarkably, although CD44 and CD49f changes between primary tumors and lymph node metastasis did not impact survival, the cases that were positive for P-cadherin in lymph node metastases being negative in the primary tumor, presented the worst disease-free and overall survival of the whole series. Accordingly, negative cases for this marker in the lymph nodes with positive expression in the matched breast carcinoma demonstrated a better prognosis, which overlapped with tumors that were negative in both sites. P-cadherin and CD49f gain of expression was mainly found in triple-negative carcinomas. Our results indicate for the first time that the evaluation of P-cadherin expression in lymph node metastases is an important predictor of disease outcome, being a putative valuable marker for axillary-based breast cancer decisions in the clinical practice.

Pax-5 is a potent regulator of E-cadherin and breast cancer malignant processes

Pax-5, an essential transcription factor for B lymphocyte development, has been linked with the development and progression of lymphoid cancers and carcinoma. In contrast to B-cell cancer lesions, the specific expression signatures and roles of Pax-5 in breast cancer progression are relatively unknown. In the present study, we set out to profile Pax-5 expression in mammary tissues and elucidate the cellular and molecular roles of Pax-5 in breast cancer processes. Using immunohistology on mammary tissue arrays, Pax-5 was detected in a total of 298/306 (97.6%) samples tested. Interestingly, our studies reveal that Pax-5 inhibits aggressive features and confers anti-proliferative effects in breast carcinoma cells in contrast to its oncogenic properties in B cell cancers. More precisely, Pax-5 suppressed breast cancer cell migration, invasion and tumor spheroid formation while concomitantly promoting cell adhesion properties. We also observed that Pax-5 inhibited and reversed breast cancer epithelial to mesenchymal phenotypic transitioning. Mechanistically, we found that the Pax-5 transcription factor binds and induces gene expression of E-cadherin, a pivotal regulator of epithelialisation. Globally, we demonstrate that Pax-5 is predominant expressed factor in mammary epithelial cells. We also present an important role for Pax-5 in the phenotypic transitioning processes and aggressive features associated with breast cancer malignancy and disease progression.

Breast Cancer Metastasis to the Axillary Lymph Nodes: Are Changes to the Lymph Node "Soil" Localized or Systemic?

Metastasis is a multistep process that is not well understood. Colonization of a secondary organ requires specific molecular alterations of the host microenvironment. To determine the temporal and spatial changes associated with metastatic dissemination to the axillary lymph nodes, gene expression profiles were compared between histologically normal lymph nodes from node-positive patients and tumor-free nodes from node-negative patients. Using a stringent false discovery rate correction (<0.05) for multiple hypothesis testing, we did not detect any differentially expressed genes between the lymph node groups. Thus, the presence of metastatic cells within the lymphatic system does not elicit widespread changes in gene expression through the axillary basin; rather, lymph nodes independently respond to disseminated tumor cells.

Analysis of molecular markers as predictive factors of lymph node involvement in breast carcinoma

Nodal status is the most significant independent prognostic factor in breast cancer. Identification of molecular markers would allow stratification of patients who require surgical assessment of lymph nodes from the large numbers of patients for whom this surgical procedure is unnecessary, thus leading to a more accurate prognosis. However, up to now, the reported studies are preliminary and controversial, and although hundreds of markers have been assessed, few of them have been used in clinical practice for treatment or prognosis in breast cancer. The purpose of the present study was to determine whether protein phosphatase Mg2+/Mn2+ dependent 1D, β-1,3-N-acetylglucosaminyltransferase, neural precursor cell expressed, developmentally down-regulated 9, prohibitin, phosphoinositide-3-kinase regulatory subunit 5 (PIK3R5), phosphatidylinositol-5-phosphate 4-kinase type IIα, TRF1-interacting ankyrin-related ADP-ribose polymerase 2, BCL2 associated agonist of cell death, G2 and S-phase expressed 1 and PAX interacting protein 1 genes, described as prognostic markers in breast cancer in a previous microarray study, are also predictors of lymph node involvement in breast carcinoma Reverse transcription-quantitative polymerase chain reaction analysis was performed on primary breast tumor tissues from women with negative lymph node involvement (n=27) compared with primary tumor tissues from women with positive lymph node involvement (n=23), and was also performed on primary tumors and paired lymph node metastases (n=11). For all genes analyzed, only the PIK3R5 gene exhibited differential expression in samples of primary tumors with positive lymph node involvement compared with primary tumors with negative lymph node involvement (P=0.0347). These results demonstrate that the PIK3R5 gene may be considered predictive of lymph node involvement in breast carcinoma. Although the other genes evaluated in the present study have been previously characterized to be involved with the development of distant metastases, they did not have predictive potential.

Intrinsic disorder in spondins and some of their interacting partners

Spondins, which are proteins that inhibit and promote adherence of embryonic cells so as to aid axonal growth are part of the thrombospondin-1 family. Spondins function in several important biological processes, such as apoptosis, angiogenesis, etc. Spondins constitute a thrombospondin subfamily that includes F-spondin, a protein that interacts with Aβ precursor protein and inhibits its proteolytic processing; R-spondin, a 4-membered group of proteins that regulates Wnt pathway and have other functions, such as regulation of kidney proliferation, induction of epithelial proliferation, the tumor suppressant action; M-spondin that mediates mechanical linkage between the muscles and apodemes; and the SCO-spondin, a protein important for neuronal development. In this study, we investigated intrinsic disorder status of human spondins and their interacting partners, such as members of the LRP family, LGR family, Frizzled family, and several other binding partners in order to establish the existence and importance of disordered regions in spondins and their interacting partners by conducting a detailed analysis of their sequences, finding disordered regions, and establishing a correlation between their structure and biological functions.

Breast Cancer Malignant Processes are Regulated by Pax-5 Through the Disruption of FAK Signaling Pathways

The study of genetic factors regulating breast cancer malignancy is a top priority to mitigate the morbidity and mortality associated with this disease. One of these factors, Pax-5, modulates cancer aggressiveness through the regulation of various components of the epithelial to mesenchymal transitioning (EMT) process. We have previously reported that Pax-5 expression profiles in cancer tissues inversely correlate with those of the Focal Adhesion Kinase (FAK), a potent activator of breast cancer malignancy. In this study, we set out to elucidate the molecular and regulatory relationship between Pax-5 and FAK in breast cancer processes. Interestingly, we found that Pax-5 mediated suppression of breast cancer cell migration is dependent of FAK activity. Our mechanistic examination revealed that Pax-5 inhibits FAK expression and activation. We also demonstrate that Pax-5 is a potent modulator of FAK repressors (p53 and miR-135b) and activator (NFκB) which results in the overall suppression of FAK-mediated signaling cascades. Altogether, our findings bring more insight to the molecular triggers regulating phenotypic transitioning process and signaling cascades leading to breast cancer progression.

Integrative analysis of mutational and transcriptional profiles reveals driver mutations of metastatic breast cancers

Despite the explosion in the numbers of cancer genomic studies, metastasis is still the major cause of cancer mortality. In breast cancer, approximately one-fifth of metastatic patients survive 5 years. Therefore, detecting the patients at a high risk of developing distant metastasis at first diagnosis is critical for effective treatment strategy. We hereby present a novel systems biology approach to identify driver mutations escalating the risk of metastasis based on both exome and RNA sequencing of our collected 78 normal-paired breast cancers. Unlike driver mutations occurring commonly in cancers as reported in the literature, the mutations detected here are relatively rare mutations occurring in less than half metastatic samples. By supposing that the driver mutations should affect the metastasis gene signatures, we develop a novel computational pipeline to identify the driver mutations that affect transcription factors regulating metastasis gene signatures. We identify driver mutations in ADPGK, NUP93, PCGF6, PKP2 and SLC22A5, which are verified to enhance cancer cell migration and prompt metastasis with in vitro experiments. The discovered somatic mutations may be helpful for identifying patients who are likely to develop distant metastasis.

A contrasting function for miR-137 in embryonic mammogenesis and adult breast carcinogenesis

MicroRNAs are differentially expressed in breast cancer cells and have been implicated in cancer formation, tumour invasion and metastasis. We investigated the miRNA expression profiles in the developing mammary gland. MiR-137 was expressed prominently in the developing mammary gland. When the miR-137 was over-expressed in the embryo, the mammary epithelium became thickened. Moreover, genes associated with mammary gland formation such as Tbx3 and Lef1 were not expressed. This suggests that miR-137 induces gland formation and invasion. When miR-137 was over-expressed in MDA-MB-231 cells, their ability to form tumours in adult mice was significantly reduced. These data support miR-137 decides epithelial cell behavior in the human breast cancer. It also suggests that miR-137 is a potential therapeutic target for amelioration of breast cancer progression.

COL11A1 is overexpressed in recurrent non-small cell lung cancer and promotes cell proliferation, migration, invasion and drug resistance

Collagen type XI α1 (COL11A1), a minor fibrillar collagen, has been demonstrated to be involved in cell proliferation, migration and the tumorigenesis of many human malignancies. Previous studies have shown that COL11A1 may be a valuable diagnostic marker for non-small cell lung carcinoma (NSCLC). However, its biological function in NSCLC progression remains largely unclear. In the present study, we investigated the expression levels of COL11A1 in different human NSCLC samples, and found that COL11A1 was overexpressed in NSCLC with lymph node metastasis and in recurrent NSCLC tissues. We also revealed that COL11A1 promoted the cell proliferation, migration and invasion of NSCLC cell lines in vitro. Furthermore, our results highlighted the importance of COL11A1 in chemoresistance to cisplatin. Mechanistically, we found that the effects of the overexpression of COL11A1 in NSCLC cells were mediated by Smad signaling. Collectively, our findings suggest that COL11A1 may sever as a biomarker for metastatic NSCLC, and can be used to predict recurrence after surgical resection. Therapeutic approaches targeting COL11A1 may facilitate the optimization of cisplatin treatment of NSCLC by overcoming chemoresistance.

An exact algorithm for finding cancer driver somatic genome alterations: the weighted mutually exclusive maximum set cover problem

Background

The mutual exclusivity of somatic genome alterations (SGAs), such as somatic mutations and copy number alterations, is an important observation of tumors and is widely used to search for cancer signaling pathways or SGAs related to tumor development. However, one problem with current methods that use mutual exclusivity is that they are not signal-based; another problem is that they use heuristic algorithms to handle the NP-hard problems, which cannot guarantee to find the optimal solutions of their models.

Method

In this study, we propose a novel signal-based method that utilizes the intrinsic relationship between SGAs on signaling pathways and expression changes of downstream genes regulated by pathways to identify cancer signaling pathways using the mutually exclusive property. We also present a relatively efficient exact algorithm that can guarantee to obtain the optimal solution of the new computational model.

Results

We have applied our new model and exact algorithm to the breast cancer data. The results reveal that our new approach increases the capability of finding better solutions in the application of cancer research. Our new exact algorithm has a time complexity of \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$O^{*}(1.325^{m})$$\end{document}O∗(1.325m)(Note: Following the recent convention, we use a star * to represent that the polynomial part of the time complexity is neglected), which has solved the NP-hard problem of our model efficiently.

Conclusion

Our new method and algorithm can discover the true causes behind the phenotypes, such as what SGA events lead to abnormality of the cell cycle or make the cell metastasis lose control in tumors; thus, it identifies the target candidates for precision (or target) therapeutics.

Electronic supplementary material

The online version of this article (doi:10.1186/s13015-016-0073-9) contains supplementary material, which is available to authorized users.

Upregulation of spondin-2 predicts poor survival of colorectal carcinoma patients

Colorectal cancer (CRC) is the third and second most common cancer in males and females worldwide, respectively. Spondin-2 is a conserved secreted extracellular matrix protein and a candidate cancer biomarker. Here we found that Spondin-2 mRNA was upregulated in CRC tissues using quantitative RT-PCR and data-mining of public Oncomine microarray datasets. Spondin-2 protein was increased in CRC tissues, as revealed by immunohistochemistry analyses of two tissue microarrays containing 180 cases. Spondin-2 gene expression was significantly associated with CRC stage, T stage, M stage and Dukes stage, while its protein was associated with age and M stage. Kaplan-Meier analysis revealed that the upregulated Spondin-2 mRNA and protein predicted poor prognosis of CRC patients. Univariate and multivariate Cox regression analyses indicated that grade, recurrence, N stage and high Spondin-2 were independent predictors of overall survival of CRC patients. ELISA revealed that plasma Spondin-2 was upregulated in CRC and dropped after surgery. Receiver operating characteristic curve analysis demonstrated that plasma Spondin-2 has superior predictive performance for CRC with an area under the curve of 0.959 and the best sensitivity/specificity of 100%/90%. Furthermore, ectopic expression of Spondin-2 enhanced colon cancer cell proliferation. Spondin-2 could be an independent diagnostic and prognostic biomarker of colon cancer.

Molecular Concordance Between Primary Breast Cancer and Matched Metastases

Clinical management of breast cancer is increasingly personalized and based on molecular profiling. Often, primary tumors are used as proxies for systemic disease at the time of recurrence. However, recent studies have revealed substantial discordances between primary tumors and metastases, both with respect to traditional clinical treatment targets and on the genomic and transcriptomic level. With the increasing use of molecularly targeted therapy, discordance of actionable molecular targets between primary tumors and recurrences can result in nonoptimal treatment or unnecessary side effects. The purpose of this review is to illuminate the extent of cancer genome evolution through disease progression and the degree of molecular concordance between primary breast cancers and matched metastases. We present an overview of the most prominent studies investigating the expression of endocrine receptors, transcriptomics, and genome aberrations in primary tumors and metastases. In conclusion, biopsy of metastatic lesions at recurrence of breast cancer is encouraged to provide optimal treatment of the disease. Furthermore, molecular profiling of metastatic tissue provides invaluable mechanistic insight into the biology underlying metastatic progression and has the potential to identify novel, potentially druggable, drivers of progression.

Survival of breast cancer patients according to changes in expression of selected markers between primary tumor and lymph node metastases

BACKGROUND

The differences between primary and metastatic tumor cells might be important for treatment selection and prognostication.

MATERIALS & METHODS

Expression of ER, PR, HER2, CK5/6, EGFR, Ki-67, Ep-CAM, P-cadherin, CD24, CD44, ALDH was assessed immunohistochemically in primary tumor (T) and corresponding synchronous nodal metastases (LNM) in 156 invasive ductal breast cancer patients (T ≥1, N ≥1, M0).

RESULTS

Independent negative prognostic factors for disease-free survival were pN3, ALDH immunopositivity in LNM, nonluminal A subtype in LNM, reduction of Ep-CAM expression in LNM, lack of changes or enhancement of CK5/6 and ALDH expression in LNM.

DISCUSSION

Our results suggest that in some cases expression of markers in lymph node metastases might bring additional prognostic information to that obtained from primary tumor.

Development of novel real-time PCR methodology for quantification of COL11A1 mRNA variants and evaluation in breast cancer tissue specimens

Background

Collagen XI is a key structural component of the extracellular matrix and consists of three alpha chains. One of these chains, the α1 (XI), is encoded by the COL11A1 gene and is transcribed to four different variants at least (A, B, C and E) that differ in the propensity to N-terminal domain proteolysis and potentially in the way the extracellular matrix is arranged. This could affect the ability of tumor cells to invade the remodeled stroma and metastasize. No study in the literature has so far investigated the expression of these four variants in breast cancer nor does a method for their accurate quantitative detection exist.

Methods

We developed a conventional PCR for the general detection of the general COL11A1 transcript and real-time qPCR methodologies with dual hybridization probes in the LightCycler platform for the quantitative determination of the variants. Data from 90 breast cancer tissues with known histopathological features were collected.

Results

The general COL11A1 transcript was detected in all samples. The developed methodologies for each variant were rapid as well as reproducible, sensitive and specific. Variant A was detected in 30 samples (33 %) and variant E in 62 samples (69 %). Variants B and C were not detected at all. A statistically significant correlation was observed between the presence of variant E and lymph nodes involvement (p = 0.037) and metastasis (p = 0.041).

Conclusions

With the newly developed tools, the possibility of inclusion of COL11A1 variants as prognostic biomarkers in emerging multiparameter technologies examining tissue RNA expression should be further explored.

Electronic supplementary material

The online version of this article (doi:10.1186/s12885-015-1725-8) contains supplementary material, which is available to authorized users.