Genetic profiling of chromosome 1 in breast cancer: mapping of regions of gains and losses and identification of candidate genes on 1q

HyGAnno: hybrid graph neural network-based cell type annotation for single-cell ATAC sequencing data

Reliable cell type annotations are crucial for investigating cellular heterogeneity in single-cell omics data. Although various computational approaches have been proposed for single-cell RNA sequencing (scRNA-seq) annotation, high-quality cell labels are still lacking in single-cell sequencing assay for transposase-accessible chromatin (scATAC-seq) data, because of extreme sparsity and inconsistent chromatin accessibility between datasets. Here, we present a novel automated cell annotation method that transfers cell type information from a well-labeled scRNA-seq reference to an unlabeled scATAC-seq target, via a parallel graph neural network, in a semi-supervised manner. Unlike existing methods that utilize only gene expression or gene activity features, HyGAnno leverages genome-wide accessibility peak features to facilitate the training process. In addition, HyGAnno reconstructs a reference-target cell graph to detect cells with low prediction reliability, according to their specific graph connectivity patterns. HyGAnno was assessed across various datasets, showcasing its strengths in precise cell annotation, generating interpretable cell embeddings, robustness to noisy reference data and adaptability to tumor tissues.

Low-Pass Genomic Sequencing Reveals Novel Subtypes of Pancreatic Cystic Neoplasms

BACKGROUND

Over the years, the detection rate of pancreatic cystic neoplasms (PCNs) has significantly increased; however, the differential diagnosis and identification of high-risk PCNs remain challenging. We sought to investigate whether chromosomal instability (CIN) features in cell-free DNA in the cystic fluid of PCNs could help to identify high-risk PCNs.

METHODS

Pancreatic cystic fluid samples from 102 patients with PCNs were intraoperatively collected for detection of CIN using an ultrasensitive chromosomal aneuploidy detector. Clinical and imaging data were retrospectively collected, and statistical analysis was performed to assess the potential role of CIN in clinical practice.

RESULTS

CIN was investigated in a total of 100 patients. Sixteen of 26 serous cystic cystadenomas (SCAs) harbored deletions of chr3p and/or chr6p, whereas low rates of CIN were detected in mucinous cystic neoplasms. Most malignant PCNs presented with more than one type of CIN; amplification of chr1q and chr8q found in nine and seven of 11 malignant PCNs (81.8% and 63.6%), respectively, could aid in distinguishing high-risk IPMNs from low-risk ones, with a higher sensitivity than imaging. A combination of the mural nodule imaging feature and amplification of chr1q and chr8q achieved a sensitivity of 70.0% and a specificity of 82.4% in identifying high-risk IPMNs.

CONCLUSIONS

Our work revealed the distinct CIN signature of different types of PCNs. Deletions of chr3p and chr6p defined a subtype of SCAs. Gains of chr1q and chr8q were associated with insidious malignant PCNs and helped identify high-risk IPMNs.

KDEL Receptor Trafficking to the Plasma Membrane Is Regulated by ACBD3 and Rab4A-GTP

KDEL receptor-1 maintains homeostasis in the early secretory pathway by capturing and retrieving ER chaperones to the ER during heavy secretory activity. Unexpectedly, a fraction of the receptor is also known to reside in the plasma membrane (PM), although it is largely unknown exactly how the KDEL receptor gets exported from the Golgi and travels to the PM. We have previously shown that a Golgi scaffolding protein (ACBD3) facilitates KDEL receptor localization at the Golgi via the regulating cargo wave-induced cAMP/PKA-dependent signaling pathway. Upon endocytosis, surface-expressed KDEL receptor undergoes highly complex itineraries through the Golgi and the endo-lysosomal compartments, where the endocytosed receptor utilizes Rab14A- and Rab11A-positive recycling endosomes and clathrin-decorated tubulovesicular carriers. In this study, we sought to investigate the mechanism through which the KDEL receptor gets exported from the Golgi en route to the PM. We report here that ACBD3 depletion results in greatly increased trafficking of KDEL receptor to the PM via Rab4A-positive tubular carriers emanating from the Golgi. Expression of constitutively activated Rab4A mutant (Q72L) increases the surface expression of KDEL receptor up to 2~3-fold, whereas Rab4A knockdown or the expression of GDP-locked Rab4A mutant (S27N) inhibits KDEL receptor targeting of the PM. Importantly, KDELR trafficking from the Golgi to the PM is independent of PKA- and Src kinase-mediated mechanisms. Taken together, these results reveal that ACBD3 and Rab4A play a key role in regulating KDEL receptor trafficking to the cell surface.

Effects of Human Papilloma Virus E6/E7 Oncoproteins on Genomic Structure in Head and Neck Squamous Cell Carcinoma

Human Papilloma Virus (HPV) is highly prevalent within the U.S., with studies estimating that over 80% of individuals will contract the virus in their lifetime. HPV is considered a primary risk factor for the development and progression of oropharyngeal cancers. The impact of the HPV virus's E6 and E7 oncoproteins on cellular signaling pathways and genomic integration has been extensively characterized. Indirect genomic effects; however, remain relatively unidentified. In this study, we analyzed 83 HPV+ Head and Neck Squamous Cell Carcinoma (HNSCC) patients of varying HPV types. Expression counts of the HPV E6 and E7 oncogenes were estimated across samples and correlated with genomic mutational classes. High expression of E6 and E7 oncoproteins was associated with a greater number of total point mutations, especially on chromosomes 1, 11, and 17, which have been implicated in HPV-mediated cancers in previous studies. Samples with high E6 and E7 expression also exhibited more frequent non-clustered structural variation and a lack of clustered variation altogether. Copy number segments were present with fewer number of repeats in high E6 and E7 expression samples, which is known to correlate with decreased expression of affected genes. E6 and E7 expression was associated with increased activity of several cellular pathways associated in oncogenesis and telomere maintenance. In comprehensively characterizing the effects of the HPV oncoproteins on the human genome, potential mechanisms of HNSCC pathogenesis may be further elucidated.

ACBD3 Bioinformatic Analysis and Protein Expression in Breast Cancer Cells

ACBD3 overexpression has previously been found to correlate with worse prognosis for breast cancer patients and, as an incredibly diverse protein in both function and cellular localisation, ACBD3 may have a larger role in breast cancer than previously thought. This study further investigated ACBD3′s role in breast cancer. Bioinformatic databases were queried to characterise ACBD3 expression and mutation in breast cancer and to investigate how overexpression affects breast cancer patient outcomes. Immunohistochemistry was carried out to examine ACBD3 location within cells and tissue structures. ACBD3 was more highly expressed in breast cancer than in any other cancer or matched normal tissue, and expression over the median level resulted in reduced relapse-free, overall, and distant metastasis-free survival for breast cancer patients as a whole, with some differences observed between subtypes. IHC analysis found that ACBD3 levels varied based on hormone receptor status, indicating that ACBD3 could be a candidate biomarker for poor patient prognosis in breast cancer and may possibly be a biomarker for ER signal reprogramming of precancerous breast tissue.

PPP2R2B downregulation is associated with immune evasion and predicts poor clinical outcomes in triple-negative breast cancer

Background

Although immune checkpoint blockade has emerged as a novel promising strategy for triple-negative breast cancer (TNBC), many patients fail response or acquire resistance to current agents. Consequently, our focus need to shift toward alternative inhibitory targets, predictor for responsiveness, and immune suppressive mechanisms.

Methods

In this study, we performed systematic bioinformatics analyses to identify PPP2R2B as a robust tumor suppressor in TNBC. Meanwhile, breast cancer progression cell line model was applied in our research. Quantitative real-time PCR assay (Q-PCR) was carried out to assess the role of PPP2R2B in the onset and progression of breast cancer. Furthermore, we validated the effect of PPP2R2B on immune activity via in vitro experiments based on macrophages. To further decipher the roles of PPP2R2B in TNBC, we investigated the transcriptome level, genomic profiles, and its clinical prognostic value.

Results

In TNBC tissues, PPP2R2B expression was significantly downregulated compared to normal breast tissues. Kaplan‐Meier survival analysis revealed that patients with low PPP2R2B expression had shorter survival time than those with high PPP2R2B expression. Q-PCR analysis suggested that PPP2R2B downregulation could play a key role in breast-cancer initiation and progression. Additionally, our findings showed that PPP2R2B was positively related with CD8 T cells, CD4 Th1 helper cells, and M1 macrophages, but negatively related with M2 macrophages. Subsequent results identified that PPP2R2B was strongly related with immune inhibitor genes (GZMA, PRF1, and IFNG), which could improve T lymphocytes antitumor function and restrict immune evasion. Meanwhile, T cell receptor signaling pathway and antigen processing and presentation signaling pathway were significantly suppressed in low PPP2R2B expression group. Afterwards, distinct subgroups based on PPP2R2B expression exhibited several unique features in somatic mutations, copy numbers alterations, extent of copy number burden, and promoter methylation level.

Conclusion

Our results indicated that PPP2R2B could serve as a promising biomarker for TNBC, and help predict immunotherapeutic response and guide personalized strategies in TNBC treatment.

Pan-cancer driver copy number alterations identified by joint expression/CNA data analysis

Analysis of large gene expression datasets from biopsies of cancer patients can identify co-expression signatures representing particular biomolecular events in cancer. Some of these signatures involve genomically co-localized genes resulting from the presence of copy number alterations (CNAs), for which analysis of the expression of the underlying genes provides valuable information about their combined role as oncogenes or tumor suppressor genes. Here we focus on the discovery and interpretation of such signatures that are present in multiple cancer types due to driver amplifications and deletions in particular regions of the genome after doing a comprehensive analysis combining both gene expression and CNA data from The Cancer Genome Atlas.

Novel Insights Into Epigenetic Reprogramming and Destabilization of Pericentromeric Heterochromatin in Cancer

Pericentromeric heterochromatin is maintained in a condensed structure by repressive epigenetic control mechanisms and perturbation of these may cause diseases. The chromosome 1q12 region harbors the largest pericentromeric heterochromatin domain in the genome and is among the most common breakpoints in both solid and hematopoietic cancers. Furthermore, the 1q arm is frequently amplified in cancer and this may support tumorigenesis by increasing the dosage of the many oncogenes of this genomic region. Recent studies have provided insight into the mechanisms leading to loss of 1q12 stability and 1q amplification and DNA hypomethylation seems to play a prominent role. This may be the result of decreased activity of DNA methyltransferases and instrumental for 1q12 destabilization or arise secondary to perturbation of other important epigenetic mechanisms that control repression of pericentromeric heterochromatin. Polycomb proteins were recently demonstrated to epigenetically reprogram demethylated 1q12 pericentromeric heterochromatin in premalignant and malignant cells to form large subnuclear structures known as polycomb bodies. This may influence the regulation and stability of 1q12 pericentromeric heterochromatin and/or the distribution of polycomb factors to support tumorigenesis. This review will discuss recent insight into the epigenetic perturbations causing the destabilization of 1q12 pericentromeric heterochromatin and its possible implications for tumor biology.

Single-Cell NGS-Based Analysis of Copy Number Alterations Reveals New Insights in Circulating Tumor Cells Persistence in Early-Stage Breast Cancer

Simple Summary

Circulating tumor cells (CTCs) are crucial for the identification of patients with a higher risk of relapse, including those diagnosed with breast cancer (BC). The aim of this study was to explore their molecular aspects in 11 early-stage BC patients during patient management, focusing on copy number alterations (CNAs) and exploiting a single-CTC next-generation sequencing approach. CTCs showed different degrees of aberration based on access time. Moreover, CTCs, in particular those persisting even months after tumor resection, shared CNAs with matched tumor tissue. Enrichment analyses of CNAs on CTCs highlighted peculiar aberrations, especially associated with interferon (IFN)-associated terms. The study of CTCs CNAs can provide information about the molecular mechanisms involving CTC-related processes and their survival ability in occult niches, supporting the goal of exploiting their application in patients’ surveillance and follow-up.

Abstract

Circulating tumor cells (CTCs) are a rare population of cells representing a key player in the metastatic cascade. They are recognized as a validated tool for the identification of patients with a higher risk of relapse, including those diagnosed with breast cancer (BC). However, CTCs are characterized by high levels of heterogeneity that also involve copy number alterations (CNAs), structural variations associated with gene dosage changes. In this study, single CTCs were isolated from the peripheral blood of 11 early-stage BC patients at different time points. A label-free enrichment of CTCs was performed using OncoQuick, and single CTCs were isolated using DEPArray. Libraries were prepared from single CTCs and DNA extracted from matched tumor tissues for a whole-genome low-coverage next-generation sequencing (NGS) analysis using the Ion Torrent S5 System. The analysis of the CNA burden highlighted that CTCs had different degrees of aberration based on the time point and subtype. CTCs were found even six months after surgery and shared CNAs with matched tumor tissue. Tumor-associated CNAs that were recurrent in CTCs were patient-specific, and some alterations involved regions associated with BC and survival (i.e., gains at 1q21-23 and 5p15.33). The enrichment analysis emphasized the involvement of aberrations of terms, associated in particular with interferon (IFN) signaling. Collectively, our findings reveal that these aberrations may contribute to understanding the molecular mechanisms involving CTC-related processes and their survival ability in occult niches, supporting the goal of exploiting their application in patients’ surveillance and follow-up.

A Pygopus 2-Histone Interaction Is Critical for Cancer Cell Dedifferentiation and Progression in Malignant Breast Cancer

Pygopus 2 (Pygo2) is a coactivator of Wnt/β-catenin signaling that can bind bi- or trimethylated lysine 4 of histone-3 (H3K4me^2/3) and participate in chromatin reading and writing. It remains unknown whether the Pygo2-H3K4me^2/3 association has a functional relevance in breast cancer progression in vivo. To investigate the functional relevance of histone-binding activity of Pygo2 in malignant progression of breast cancer, we generated a knock-in mouse model where binding of Pygo2 to H3K4me^2/3 was rendered ineffective. Loss of Pygo2-histone interaction resulted in smaller, differentiated, and less metastatic tumors, due, in part, to decreased canonical Wnt/β-catenin signaling. RNA- and ATAC-sequencing analyses of tumor-derived cell lines revealed downregulation of TGFβ signaling and upregulation of differentiation pathways such as PDGFR signaling. Increased differentiation correlated with a luminal cell fate that could be reversed by inhibition of PDGFR activity. Mechanistically, the Pygo2-histone interaction potentiated Wnt/β-catenin signaling, in part, by repressing the expression of Wnt signaling antagonists. Furthermore, Pygo2 and β-catenin regulated the expression of miR-29 family members, which, in turn, repressed PDGFR expression to promote dedifferentiation of wild-type Pygo2 mammary epithelial tumor cells. Collectively, these results demonstrate that the histone binding function of Pygo2 is important for driving dedifferentiation and malignancy of breast tumors, and loss of this binding activates various differentiation pathways that attenuate primary tumor growth and metastasis formation. Interfering with the Pygo2-H3K4me^2/3 interaction may therefore serve as an attractive therapeutic target for metastatic breast cancer. SIGNIFICANCE: Pygo2 represents a potential therapeutic target in metastatic breast cancer, as its histone-binding capability promotes β-catenin-mediated Wnt signaling and transcriptional control in breast cancer cell dedifferentiation, EMT, and metastasis.

BCL9/STAT3 regulation of transcriptional enhancer networks promote DCIS progression

The molecular processes by which some human ductal carcinoma in situ (DCIS) lesions advance to the more aggressive form, while others remain indolent, are largely unknown. Experiments utilizing a patient-derived (PDX) DCIS Mouse INtraDuctal (MIND) animal model combined with ChIP-exo and RNA sequencing revealed that the formation of protein complexes between B Cell Lymphoma-9 (BCL9), phosphoserine 727 STAT3 (PS-727-STAT3) and non-STAT3 transcription factors on chromatin enhancers lead to subsequent transcription of key drivers of DCIS malignancy. Downregulation of two such targets, integrin β3 and its associated metalloproteinase, MMP16, resulted in a significant inhibition of DCIS invasive progression. Finally, in vivo targeting of BCL9, using rosemary extract, resulted in significant inhibition of DCIS malignancy in both cell line and PDX DCIS MIND animal models. As such, our studies provide compelling evidence for future testing of rosemary extract as a chemopreventive agent in breast cancer.

Banding Together: A Systematic Comparison of The Cancer Genome Atlas and the Mitelman Databases

Cytogenetic aberrations at the single-cell level represent an important characteristic of cancer cells relevant to tumor evolution and prognosis. However, with the advent of The Cancer Genome Atlas (TCGA), there has been a major shift in cancer research to the use of data from aggregate cell populations. Given that tumor cells harbor hundreds to thousands of biologically relevant genetic alterations that manifest as intratumor heterogeneity, these aggregate analyses may miss alterations readily observable at single-cell resolution. Using the Mitelman Database of Chromosome Aberrations and Gene Fusions in Cancer, we developed an algorithm to parse International System for Cytogenetic Nomenclature notation for quantitative abnormalities. Comparison of the Mitelman database and TCGA demonstrated that the Mitelman database is a powerful resource, and that cytogenetic aberrations captured by traditional approaches used in Mitelman database are on par with population-based genomic analyses used in TCGA. This algorithm will help nonspecialists to overcome the challenges associated with the format and syntax of the Mitelman database. SIGNIFICANCE: A novel in silico approach compares cytogenetic data between the Mitelman database and TCGA, highlighting the advantages and limitations of both datasets.

The Great Escape: how phosphatidylinositol 4-kinases and PI4P promote vesicle exit from the Golgi (and drive cancer)

Phosphatidylinositol 4-phosphate (PI4P) is a membrane glycerophospholipid and a major regulator of the characteristic appearance of the Golgi complex as well as its vesicular trafficking, signalling and metabolic functions. Phosphatidylinositol 4-kinases, and in particular the PI4KIIIβ isoform, act in concert with PI4P to recruit macromolecular complexes to initiate the biogenesis of trafficking vesicles for several Golgi exit routes. Dysregulation of Golgi PI4P metabolism and the PI4P protein interactome features in many cancers and is often associated with tumour progression and a poor prognosis. Increased expression of PI4P-binding proteins, such as GOLPH3 or PITPNC1, induces a malignant secretory phenotype and the release of proteins that can remodel the extracellular matrix, promote angiogenesis and enhance cell motility. Aberrant Golgi PI4P metabolism can also result in the impaired post-translational modification of proteins required for focal adhesion formation and cell–matrix interactions, thereby potentiating the development of aggressive metastatic and invasive tumours. Altered expression of the Golgi-targeted PI 4-kinases, PI4KIIIβ, PI4KIIα and PI4KIIβ, or the PI4P phosphate Sac1, can also modulate oncogenic signalling through effects on TGN-endosomal trafficking. A Golgi trafficking role for a PIP 5-kinase has been recently described, which indicates that PI4P is not the only functionally important phosphoinositide at this subcellular location. This review charts new developments in our understanding of phosphatidylinositol 4-kinase function at the Golgi and how PI4P-dependent trafficking can be deregulated in malignant disease.

Copy number alterations associated with clinical features in an underrepresented population with breast cancer

Background

As the most incident tumor among women worldwide, breast cancer is a heterogeneous disease. Tremendous efforts have been made to understand how tumor characteristics as histological type, molecular subtype, and tumor microenvironment collectively influence disease diagnosis to treatment, which impact outcomes. Differences between populations and environmental and cultural factors have impacts on the origin and evolution of the disease, as well as the therapeutic challenges that arise due to these factors. We, then, compared copy number variations (CNVs) in mucinous and nonmucinous luminal breast tumors from a Brazilian cohort to investigate major CNV imbalances in mucinous tumors versus non‐mucinous luminal tumors, taking into account their clinical and pathological features.

Methods

48 breast tumor samples and 48 matched control blood samples from Brazilian women were assessed for CNVs by chromosome microarray. Logistic regression and random forest models were used in order to assess CNVs in chromosomal regions from tumors.

Results

CNVs that were identified in chromosomes 1, 5, 8, 17, 19, and 21 classify tumors according to their histological type, ethnicity, disease stage, and familial history.

Conclusion

Copy number alterations described in this study provide a better understanding of the landscape of genomic aberrations in mucinous breast cancers that are associated with clinical features.

Prognostic and clinicopathological significance of Cks1 in cancer: Evidence from a meta-analysis

BACKGROUND

Cyclin-dependent kinase subunit 1 (Cks1), as a highly conserved regulatory protein, has pleiotropic roles in cell cycle progression. As research progresses, increasingly more statistics show that Cks1 may be involved in the occurrence, development, and prognosis of a variety of tumors but the conclusions remain controversial. In addition, there has been no meta-analysis demonstrating the correlation between Cks1 and cancer. Therefore, this meta-analysis was performed to determine the prognostic and clinicopathological significance of Cks1 in various cancers.

METHODS

Systematic computer literature retrieval was conducted on the Web of Science, Embase, PubMed, CNKI, and Wanfang databases. Stata SE12.0 software was used in the quantitative meta-analysis. The hazard ratio (HR) and relative risk (RR) were pooled to assess the relationship between Cks1 expression and overall survival (OS), disease-free survival (DFS), and clinicopathological parameters.

RESULTS

Nineteen studies were included, totaling 2,224 participants. High expression of Cks1 was significantly correlated with worse OS (HR, 2.62; 95% confidence interval [CI], 2.18-3.14; p < 0.001) and poorer DFS (HR, 2.73; 95% CI, 1.83-4.08; p < 0.001). In addition, high expression of Cks1 was related to lymph node metastasis (RR, 1.59; 95% CI, 1.22-2.07; p = 0.001) and advanced T stage (RR, 1.14; 95% CI, 1.04-1.25; p = 0.005).

CONCLUSIONS

High Cks1 expression predicted poorer prognosis and worse clinicopathological parameters in various cancers. Increased Cks1 could be a significant prognostic biomarker for poor survival in patients with various cancers.

Tumor mutational burden is a determinant of immune-mediated survival in breast cancer

Mounting evidence supports a role for the immune system in breast cancer outcomes. The ability to distinguish highly immunogenic tumors susceptible to anti-tumor immunity from weakly immunogenic or inherently immune-resistant tumors would guide development of therapeutic strategies in breast cancer. Genomic, transcriptomic and clinical data from The Cancer Genome Atlas (TCGA) and Molecular Taxonomy of Breast Cancer International Consortium (METABRIC) breast cancer cohorts were used to examine statistical associations between tumor mutational burden (TMB) and the survival of patients whose tumors were assigned to previously-described prognostic immune subclasses reflecting favorable, weak or poor immune-infiltrate dispositions (FID, WID or PID, respectively). Tumor immune subclasses were associated with survival in patients with high TMB (TMB-Hi, P < 0.001) but not in those with low TMB (TMB-Lo, P = 0.44). This statistical relationship was confirmed in the METABRIC cohort (TMB-Hi, P = 0.047; TMB-Lo, P = 0.39), and also found to hold true in the more-indolent Luminal A tumor subtype (TMB-Hi, P = 0.011; TMB-Lo, P = 0.91). In TMB-Hi tumors, the FID subclass was associated with prolonged survival independent of tumor stage, molecular subtype, age and treatment. Copy number analysis revealed the reproducible, preferential amplification of chromosome 1q immune-regulatory genes in the PID immune subclass. These findings demonstrate a previously unappreciated role for TMB as a determinant of immune-mediated survival of breast cancer patients and identify candidate immune-regulatory mechanisms associated with immunologically cold tumors. Immune subtyping of breast cancers may offer opportunities for therapeutic stratification.

Circulating tumor cells: potential markers of minimal residual disease in ovarian cancer? a study of the OVCAD consortium

Purpose

In 75% of ovarian cancer patients the tumor mass is completely eradicated by established surgical and cytotoxic treatment; however, the majority of the tumors recur within 24 months. Here we investigated the role of circulating tumor cells (CTCs) indicating occult tumor load, which remains inaccessible by established diagnostics.

Experimental design

Blood was taken at diagnosis (baseline samples, n = 102) and six months after completion of adjuvant first-line chemotherapy (follow-up samples; n = 78). CTCs were enriched by density gradient centrifugation. A multi-marker immunostaining was established and further complemented by FISH on CTCs and tumor/metastasis tissues using probes for stem-cell like fusion genes MECOM and HHLA1.

Results

CTCs were observed in 26.5% baseline and 7.7% follow-up blood samples at a mean number of 12.4 and 2.8 CTCs per ml blood, respectively. Baseline CTCs indicated a higher risk of death in R0 patients with complete gross resection (univariate: HR 2.158, 95% CI 1.111–4.191, p = 0.023; multivariate: HR 2.720, 95% CI 1.340–5.522, p = 0.006). At follow-up, the presence of CTCs was associated with response to primary treatment as assessed using RECIST criteria. Chromosomal gains at MECOM and HHLA1 loci suggest that the observed cells were cancer cells and reflect pathophysiological decisive chromosomal aberrations of the primary and metastatic tumors.

Conclusions

Our data suggest that CTCs detected by the multi-marker protein panel and/or MECOM/HHLA1 FISH represent minimal residual disease in optimally debulked ovarian cancer patients. The role of CTCs cells especially for clinical therapy stratification of the patients has to be validated in consecutive larger studies applying standardized treatment schemes.

Predicting pathway cross-talks in ankylosing spondylitis through investigating the interactions among pathways

Given that the pathogenesis of ankylosing spondylitis (AS) remains unclear, the aim of this study was to detect the potentially functional pathway cross-talk in AS to further reveal the pathogenesis of this disease. Using microarray profile of AS and biological pathways as study objects, Monte Carlo cross-validation method was used to identify the significant pathway cross-talks. In the process of Monte Carlo cross-validation, all steps were iterated 50 times. For each run, detection of differentially expressed genes (DEGs) between two groups was conducted. The extraction of the potential disrupted pathways enriched by DEGs was then implemented. Subsequently, we established a discriminating score (DS) for each pathway pair according to the distribution of gene expression levels. After that, we utilized random forest (RF) classification model to screen out the top 10 paired pathways with the highest area under the curve (AUCs), which was computed using 10-fold cross-validation approach. After 50 bootstrap, the best pairs of pathways were identified. According to their AUC values, the pair of pathways, antigen presentation pathway and fMLP signaling in neutrophils, achieved the best AUC value of 1.000, which indicated that this pathway cross-talk could distinguish AS patients from normal subjects. Moreover, the paired pathways of SAPK/JNK signaling and mitochondrial dysfunction were involved in 5 bootstraps. Two paired pathways (antigen presentation pathway and fMLP signaling in neutrophil, as well as SAPK/JNK signaling and mitochondrial dysfunction) can accurately distinguish AS and control samples. These paired pathways may be helpful to identify patients with AS for early intervention.

Characterization of immortalized human mammary epithelial cell line HMEC 2.6

Primary human mammary epithelial cells have a limited life span which makes it difficult to study them in vitro for most purposes. To overcome this problem, we have developed a cell line that was immortalized using defined genetic elements, and we have characterized this immortalized non-tumorigenic human mammary epithelial cell line to establish it as a potential model system. human mammary epithelial cells were obtained from a healthy individual undergoing reduction mammoplasty at SIU School of Medicine. The cells were transduced with CDK4R24C followed by transduction with human telomerase reverse transcriptase. Post all manipulation, the cells displayed a normal cell cycle phase distribution and were near diploid in nature, which was confirmed by flow cytometry and karyotyping. In vitro studies showed that the cells were anchorage dependent and were non-invasive in nature. The cell line expressed basal epithelial markers such as cytokeratin 7, CD10, and p63 and was negative for the expression of estrogen receptor and progesterone receptor. Upon G-band karyotyping, the cell line displayed the presence of a few cytogenic abnormalities, including trisomy 20 and trisomy 7, which are also commonly present in other immortalized mammary cell lines. Furthermore, the benign nature of these cells was confirmed by multiple in vitro and in vivo experiments. Therefore, we think that this cell line could serve as a good model to understand the molecular mechanisms involved in the development and progression of breast cancer and to also assess the effect of novel therapeutics on human mammary epithelial cells.

Amplification of 1q32.1 Refines the Molecular Classification of Endometrial Carcinoma

Purpose: Molecular classification of endometrial cancer identified distinct molecular subgroups. However, the largest subset of endometrial cancers remains poorly characterized and is referred to as the "nonspecific molecular profile" (NSMP) subgroup. Here, we aimed at refining the classification of this subgroup by profiling somatic copy-number aberrations (SCNAs).Experimental Design: SCNAs were analyzed in 141 endometrial cancers using whole-genome SNP arrays and pooled with 361 endometrial cancers from The Cancer Genome Atlas. Genomic Identification of Significant Targets in Cancer (GISTIC) identified statistically enriched SCNAs and penalized Cox regression assessed survival effects. The prognostic significance of relevant SCNAs was validated using multiplex ligation-dependent probe amplification in 840 endometrial cancers from the PORTEC-1/2 trials. Copy-number status of genes was correlated with gene expression to identify potential cancer drivers. One plausible oncogene was validated in vitro using antisense oligonucleotide-based strategy.Results: SCNAs affecting chromosome 1q32.1 significantly correlated with worse relapse-free survival (RFS) in the NSMP subgroup (HR, 2.12; 95% CI, 1.26-3.59; P = 0.005). This effect was replicated in NSMP endometrial cancers from PORTEC-1/2 (HR, 2.34; 95% CI, 1.17-4.70; P = 0.017). A new molecular classification including the 1q32.1 amplification improved risk prediction of recurrence. MDM4 gene expression strongly correlated with 1q32.1 amplification. Silencing MDM4 inhibited cell growth in cell lines carrying 1q32.1 amplification, but not in those without MDM4 amplification. Vice versa, increasing MDM4 expression in nonamplified cell lines stimulated cell proliferation.Conclusions: 1q32.1 amplification was identified as a prognostic marker for poorly characterized NSMP endometrial cancers, refining the molecular classification of this subgroup. We functionally validated MDM4 as a potential oncogenic driver in the 1q32.1 region. Clin Cancer Res; 23(23); 7232-41. ©2017 AACR.

Deconvolution of DNA methylation identifies differentially methylated gene regions on 1p36 across breast cancer subtypes

Breast cancer is a complex disease consisting of four distinct molecular subtypes. DNA methylation-based (DNAm) studies in tumors are complicated further by disease heterogeneity. In the present study, we compared DNAm in breast tumors with normal-adjacent breast samples from The Cancer Genome Atlas (TCGA). We constructed models stratified by tumor stage and PAM50 molecular subtype and performed cell-type reference-free deconvolution to control for cellular heterogeneity. We identified nineteen differentially methylated gene regions (DMGRs) in early stage tumors across eleven genes (AGRN, C1orf170, FAM41C, FLJ39609, HES4, ISG15, KLHL17, NOC2L, PLEKHN1, SAMD11, WASH5P). These regions were consistently differentially methylated in every subtype and all implicated genes are localized to the chromosomal cytoband 1p36.3. Seventeen of these DMGRs were independently validated in a similar analysis of an external data set. The identification and validation of shared DNAm alterations across tumor subtypes in early stage tumors advances our understanding of common biology underlying breast carcinogenesis and may contribute to biomarker development. We also discuss evidence of the specific importance and potential function of 1p36 in cancer.

Identification of novel targets for multiple myeloma through integrative approach with Monte Carlo cross-validation analysis

More than one pathway is involved in disease development and progression, and two or more pathways may be interconnected to further affect the disease onset, as functional proteins participate in multiple pathways. Thus, identifying cross-talk among pathways is necessary to understand the molecular mechanisms of multiple myeloma (MM). Based on this, this paper looked at extracting potential pathway cross-talk in MM through an integrative approach using Monte Carlo cross-validation analysis. The gene expression library of MM (accession number: GSE6477) was downloaded from the Gene Expression Omnibus (GEO) database. The integrative approach was then used to identify potential pathway cross-talk, and included four steps: Firstly, differential expression analysis was conducted to identify differentially expressed genes (DEGs). Secondly, the DEGs obtained were mapped to the pathways downloaded from an ingenuity pathways analysis (IPA), to reveal the underlying relationship between the DEGs and pathways enriched by these DEGs. A subset of pathways enriched by the DEGs was then obtained. Thirdly, a discriminating score (DS) value for each paired pathway was computed. Lastly, random forest (RF) classification was used to identify the paired pathways based on area under the curve (AUC) and Monte Carlo cross-validation, which was repeated 50 times to explore the best paired pathways. These paired pathways were tested with another independently published MM microarray data (GSE85837), using in silico validation. Overall, 60 DEGs and 19 differential pathways enriched by DEGs were extracted. Each pathway was sorted based on their AUC values. The paired pathways, inhibition of matrix metalloproteases and EIF2 signaling pathway, indicated the best AUC value of 1.000. Paired pathways consisting of IL-8 and EIF2 signaling pathways with higher AUC of 0.975, were involved in 7 runs. Furthermore, it was validated consistently in separate microarray data sets (GSE85837). Paired pathways (inhibition of matrix metalloproteases and EIF2 signaling, IL-8 signaling and EIF2 signaling) exhibited the best AUC values and higher frequency of validation. Two paired pathways (inhibition of matrix metalloproteases and EIF2 signaling, IL-8 signaling and EIF2 signaling) were used to accurately classify MM and control samples. These paired pathways may be potential bio-signatures for diagnosis and management of MM.

Network analysis of HBV‑ and HCV‑induced hepatocellular carcinoma based on Random Forest and Monte Carlo cross‑validation

Hepatocellular carcinoma (HCC) is one of the leading causes of cancer‑associated mortality worldwide. Hepatitis B virus (HBV) and hepatitis C virus (HCV) are two common risk factors for HCC. The majority of patients with HCC present at an advanced stage and are refractory to therapy. It is important to identify a method for efficient diagnosis at early stage. In the present study gene expression profile data, generated from microarray data, were pretreated according to the annotation files. The genes were mapped to pathways of Ingenuity Pathways Analysis. Dysregulated pathways and dysregulated pathway pairs were identified and constructed into individual networks, and a main network was constructed from individual networks with several edges. Random Forest (RF) classification was introduced to calculate the area under the curve (AUC) value of this network. Subsequently, 50 runs of Monte Carlo cross‑validation were used to screen the optimal main network. The results indicated that a total of 4,929 genes were identified in the pathways and gene expression profile. By combining dysregulated pathways with Z<0.05 and dysregulated pathway pairs with Z<0.2, individual networks were constructed. The optimal main network with the highest AUC value was identified. In the HCV group, the network was identified with an AUC value of 0.98, including 41 pairs of pathways, and in the HBV group, the network was identified with an AUC value of 0.94, including eight pairs of pathways. In addition, four pairs were identified in both groups. In conclusion, the optimal networks of HCV and HBV groups were identified with the highest AUC values. The use of these networks is expected to assist in diagnosing patients effectively at an early stage.

Pathway analysis based on Monte Carlo Cross-Validation in polyarticular juvenile idiopathic arthritis

INTRODUCTION

Juvenile idiopathic arthritis (JIA) is a common chronic disease with onset before the 16 years old in a child. Polyarticular JIA has been reported as the main form of JIA in several locations. Until now, understanding of the genetic basis of JIA is incomplete. The purpose of this study was to identify pathway pairs of great potential functional relevance in the progression of polyarticular JIA.

MATERIALS AND METHODS

Microarray data of 59 peripheral blood samples from healthy children and 61 samples from polyarticular JIA were transformed to gene expression data. Differential expressed genes (DEG) between patients and normal controls were identified using Linear Models for Microarray Analysis. After performed enrichment of DEG, differential pathways were identified with Fisher's test and false discovery rate. Differential pathway pairs were constructed with random two differential pathways, and were evaluated by Random Forest classification. Monte Carlo Cross-Validation was introduced to screen the best pathway pair.

RESULTS

42 DEG with P-values<0.01 were identified. 19 differential pathways with P-values<0.01 were identified. Area under the curve (AUC) of pathway pairs was generated with RF classification. After 50 bootstraps of Monte Carlo Cross-Validation, the best pathway pair with the highest AUC value was identified, and it was the pair of tumoricidal function of hepatic natural killer cells pathway and erythropoietin signaling pathway.

CONCLUSION

These identified pathway pairs may play pivotal roles in the progress of polyarticular JIA and can be applied for diagnosis. Particular attention can be focused on them for further research.

Characterizing Dysregulated Networks in Individual Patients with Ischemic Stroke Based on Monte Carlo Cross-Validation

The purpose of this study was to introduce a new method to elucidating the molecular mechanisms in ischemic stroke. Genes from microarray data were performed enrichment to biological pathways. Dysregulated pathways and dysregulated pathway pairs were identified and constructed into networks. After Random Forest classification was performed, area under the curve (AUC) value of main network was calculated. After 50 bootstraps of Monte Carlo Cross-Validation, six pairs of pathways were found for >40 times. The best main network with AUC value = 0.735 was identified, including 14 pairs of pathways. Compared with the traditional method (gene set enrichment analysis), although a small part of pathways were shared, most of the pathways were closely related with ischemic stroke. The best network may give new insights into the underlying molecular mechanisms in ischemic stroke. It may play pivotal roles in the progression of ischemic stroke and particular attention should be focused on them for further research.

Effect of low doses of estradiol and tamoxifen on breast cancer cell karyotypes

Evidence supports a role of 17&-estradiol (E2) in carcinogenesis and the large majority of breast carcinomas are dependent on estrogen. The anti-estrogen tamoxifen (TAM) is widely used for both treatment and prevention of breast cancer; however, it is also carcinogenic in human uterus and rat liver, highlighting the profound complexity of its actions. The nature of E2- or TAM-induced chromosomal damage has been explored using relatively high concentrations of these agents, and only some numerical aberrations and chromosomal breaks have been analyzed. This study aimed to determine the effects of low doses of E2 and TAM (10(&8 )mol L(&1) and 10(&6 )mol L(&1) respectively) on karyotypes of MCF7, T47D, BT474, and SKBR3 breast cancer cells by comparing the results of conventional karyotyping and multi-FISH painting with cell proliferation. Estrogen receptor (ER)-positive (+) cells showed an increase in cell proliferation after E2 treatment (MCF7, T47D, and BT474) and a decrease after TAM treatment (MCF7 and T47D), whereas in ER& cells (SKBR3), no alterations in cell proliferation were observed, except for a small increase at 96 h. Karyotypes of both ER+ and ER& breast cancer cells increased in complexity after treatments with E2 and TAM leading to specific chromosomal abnormalities, some of which were consistent throughout the treatment duration. This genotoxic effect was higher in HER2+ cells. The ER&/HER2+ SKBR3 cells were found to be sensitive to TAM, exhibiting an increase in chromosomal aberrations. These in vitro results provide insights into the potential role of low doses of E2 and TAM in inducing chromosomal rearrangements in breast cancer cells.

Quantification of read species behavior within whole genome sequencing of cancer genomes for the stratification and visualization of genomic variation

The cancer genome is abnormal genome, and the ability to monitor its sequence had undergone a technological revolution. Yet prognosis and diagnosis remain an expert-based decision, with only limited abilities to provide machine-based decisions. We introduce a heterogeneity-based method for stratifying and visualizing whole-genome sequencing (WGS) reads. This method uses the heterogeneity within WGS reads to markedly reduce the dimensionality of next-generation sequencing data; it is available through the tool HiBS (Heterogeneity-Based Subclassification) that allows cancer sample classification. We validated HiBS using >200 WGS samples from nine different cancer types from The Cancer Genome Atlas (TCGA). With HiBS, we show progress with two WGS related issues: (i) differentiation between normal (NB) and tumor (TP) samples based solely on the information structure of their WGS data, and (ii) identification of specific regions of chromosomal amplification/deletion and their association with tumor stage. By comparing results to those obtained through available WGS analyses tools, we demonstrate some of the novelties obtained by the approach implemented in HiBS and also show nearly perfect normal/tumor classification, used to identify known and unknown chromosomal aberrations. Finally, the HiBS index has been associated with breast cancer tumor stage.

Genetic alterations of histone lysine methyltransferases and their significance in breast cancer

Histone lysine methyltransferases (HMTs), a large class of enzymes that catalyze site-specific methylation of lysine residues on histones and other proteins, play critical roles in controlling transcription, chromatin architecture, and cellular differentiation. However, the genomic landscape and clinical significance of HMTs in breast cancer remain poorly characterized. Here, we conducted a meta-analysis of approximately 50 HMTs in breast cancer and identified associations among recurrent copy number alterations, mutations, gene expression, and clinical outcome. We identified 12 HMTs with the highest frequency of genetic alterations, including 8 with high-level amplification, 2 with putative homozygous deletion, and 2 with somatic mutation. Different subtypes of breast cancer have different patterns of copy number and expression for each HMT gene. In addition, chromosome 1q contains four HMTs that are concurrently or independently amplified or overexpressed in breast cancer. Copy number or mRNA expression of several HMTs was significantly associated with basal-like breast cancer and shorter patient survival. Integrative analysis identified 8 HMTs (SETDB1, SMYD3, ASH1L, SMYD2, WHSC1L1, SUV420H1, SETDB2, and KMT2C) that are dysregulated by genetic alterations, classifying them as candidate therapeutic targets. Together, our findings provide a strong foundation for further mechanistic research and therapeutic options using HMTs to treat breast cancer.

Integrating genetics and epigenetics in breast cancer: biological insights, experimental, computational methods and therapeutic potential

BACKGROUND

Development of human cancer can proceed through the accumulation of different genetic changes affecting the structure and function of the genome. Combined analyses of molecular data at multiple levels, such as DNA copy-number alteration, mRNA and miRNA expression, can clarify biological functions and pathways deregulated in cancer. The integrative methods that are used to investigate these data involve different fields, including biology, bioinformatics, and statistics.

RESULTS

These methodologies are presented in this review, and their implementation in breast cancer is discussed with a focus on integration strategies. We report current applications, recent studies and interesting results leading to the identification of candidate biomarkers for diagnosis, prognosis, and therapy in breast cancer by using both individual and combined analyses.

CONCLUSION

This review presents a state of art of the role of different technologies in breast cancer based on the integration of genetics and epigenetics, and shares some issues related to the new opportunities and challenges offered by the application of such integrative approaches.

Expression profiling of in vivo ductal carcinoma in situ progression models identified B cell lymphoma-9 as a molecular driver of breast cancer invasion

Introduction

There are an estimated 60,000 new cases of ductal carcinoma in situ (DCIS) each year. A lack of understanding in DCIS pathobiology has led to overtreatment of more than half of patients. We profiled the temporal molecular changes during DCIS transition to invasive ductal carcinoma (IDC) using in vivo DCIS progression models. These studies identified B cell lymphoma-9 (BCL9) as a potential molecular driver of early invasion. BCL9 is a newly found co-activator of Wnt-stimulated β-catenin-mediated transcription. BCL9 has been shown to promote progression of multiple myeloma and colon carcinoma. However BCL9 role in breast cancer had not been previously recognized.

Methods

Microarray and RNA sequencing were utilized to characterize the sequential changes in mRNA expression during DCIS invasive transition. BCL9-shRNA knockdown was performed to assess the role of BCL9 in in vivo invasion, epithelial-mesenchymal transition (EMT) and canonical Wnt-signaling. Immunofluorescence of 28 patient samples was used to assess a correlation between the expression of BCL9 and biomarkers of high risk DCIS. The cancer genome atlas data were analyzed to assess the status of BCL9 gene alterations in breast cancers.

Results

Analysis of BCL9, by RNA and protein showed BCL9 up-regulation to be associated with DCIS transition to IDC. Analysis of patient DCIS revealed a significant correlation between high nuclear BCL9 and pathologic characteristics associated with DCIS recurrence: Estrogen receptor (ER) and progesterone receptor (PR) negative, high nuclear grade, and high human epidermal growth factor receptor2 (HER2). In vivo silencing of BCL9 resulted in the inhibition of DCIS invasion and reversal of EMT. Analysis of the TCGA data showed BCL9 to be altered in 26 % of breast cancers. This is a significant alteration when compared to HER2 (ERBB2) gene (19 %) and estrogen receptor (ESR1) gene (8 %). A significantly higher proportion of basal like invasive breast cancers compared to luminal breast cancers showed BCL9 amplification.

Conclusion

BCL9 is a molecular driver of DCIS invasive progression and may predispose to the development of basal like invasive breast cancers. As such, BCL9 has the potential to serve as a biomarker of high risk DCIS and as a therapeutic target for prevention of IDC.

Electronic supplementary material

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

Integrative Analysis with Monte Carlo Cross-Validation Reveals miRNAs Regulating Pathways Cross-Talk in Aggressive Breast Cancer

In this work an integrated approach was used to identify functional miRNAs regulating gene pathway cross-talk in breast cancer (BC). We first integrated gene expression profiles and biological pathway information to explore the underlying associations between genes differently expressed among normal and BC samples and pathways enriched from these genes. For each pair of pathways, a score was derived from the distribution of gene expression levels by quantifying their pathway cross-talk. Random forest classification allowed the identification of pairs of pathways with high cross-talk. We assessed miRNAs regulating the identified gene pathways by a mutual information analysis. A Fisher test was applied to demonstrate their significance in the regulated pathways. Our results suggest interesting networks of pathways that could be key regulatory of target genes in BC, including stem cell pluripotency, coagulation, and hypoxia pathways and miRNAs that control these networks could be potential biomarkers for diagnostic, prognostic, and therapeutic development in BC. This work shows that standard methods of predicting normal and tumor classes such as differentially expressed miRNAs or transcription factors could lose intrinsic features; instead our approach revealed the responsible molecules of the disease.

Amplification of Chromosome 1q Genes Encoding the Phosphoinositide Signalling Enzymes PI4KB, AKT3, PIP5K1A and PI3KC2B in Breast Cancer

Little is known about the possible oncogenic roles of genes encoding for the phosphatidylinositol 4-kinases, a family of enzymes that regulate an early step in phosphoinositide signalling. To address this issue, the mutational status of all four human phosphatidylinositol 4-kinases genes was analyzed across 852 breast cancer samples using the COSMIC data resource. Point mutations in the phosphatidylinositol 4-kinase genes were uncommon and appeared in less than 1% of the patient samples however, 62% of the tumours had increases in gene copy number for PI4KB which encodes the phosphatidylinositol 4-kinase IIIbeta isozyme. Extending this analysis to subsequent enzymes in the phosphoinositide signalling cascades revealed that the only PIP5K1A, PI3KC2B and AKT3 genes exhibited similar patterns of gene copy number variation. By comparison, gene copy number increases for established oncogenes such as EGFR and HER2/Neu were only evident in 20% of the samples. The PI4KB, PIP5K1A, PI3KC2B and AKT3 genes are related in that they all localize to chromosome 1q which is often structurally and numerically abnormal in breast cancer. These results demonstrate that a gene quartet encoding a potential phosphoinositide signalling pathway is amplified in a subset of breast cancers.

Clinical value of isoform-specific detection and targeting of AKT1, AKT2 and AKT3 in breast cancer

SUMMARY 

Overactivation of the PI3K/AKT signaling pathway is frequently reported in breast cancer, consequently inhibitors targeting this pathway are clinically useful. AKT constitutes a hub in the regulation of several cancer hallmarks, such as proliferation, survival and migration. Three AKT isoforms, named AKT1, AKT2 and AKT3, are identified in humans. AKT alterations, mainly upregulation of phosphorylated AKT in tumors may have prognostic and predictive value. Moreover, the AKT isoforms may possess partly divergent cellular functions and be upregulated in certain breast cancer subtypes, suggesting the importance of isoform-specific analyses. In conclusion, AKT isoform-specific detection and targeting in different tumor subtypes will hopefully result into a further developed personalized medicine.

Integration of mRNA expression profile, copy number alterations, and microRNA expression levels in breast cancer to improve grade definition

Defining the aggressiveness and growth rate of a malignant cell population is a key step in the clinical approach to treating tumor disease. The correct grading of breast cancer (BC) is a fundamental part in determining the appropriate treatment. Biological variables can make it difficult to elucidate the mechanisms underlying BC development. To identify potential markers that can be used for BC classification, we analyzed mRNAs expression profiles, gene copy numbers, microRNAs expression and their association with tumor grade in BC microarray-derived datasets. From mRNA expression results, we found that grade 2 BC is most likely a mixture of grade 1 and grade 3 that have been misclassified, being described by the gene signature of either grade 1 or grade 3. We assessed the potential of the new approach of integrating mRNA expression profile, copy number alterations, and microRNA expression levels to select a limited number of genomic BC biomarkers. The combination of mRNA profile analysis and copy number data with microRNA expression levels led to the identification of two gene signatures of 42 and 4 altered genes (FOXM1, KPNA4, H2AFV and DDX19A) respectively, the latter obtained through a meta-analytical procedure. The 42-based gene signature identifies 4 classes of up- or down-regulated microRNAs (17 microRNAs) and of their 17 target mRNA, and the 4-based genes signature identified 4 microRNAs (Hsa-miR-320d, Hsa-miR-139-5p, Hsa-miR-567 and Hsa-let-7c). These results are discussed from a biological point of view with respect to pathological features of BC. Our identified mRNAs and microRNAs were validated as prognostic factors of BC disease progression, and could potentially facilitate the implementation of assays for laboratory validation, due to their reduced number.

Novel genetic aberrations in breast phyllodes tumours: comparison between prognostically distinct groups

Phyllodes tumours of the breast are uncommon fibroepithelial neoplasms which pose management challenges due to difficulties in accurate prediction of clinical behaviour, as histological assessment has its limitations. Molecular studies have improved the understanding of these rare tumours but such findings are scant. We aimed to investigate genetic aberrations in phyllodes tumours stratified according to clinical behaviour, to identify potential genes contributing to disease progression. Twenty phyllodes tumours were separated into prognostically distinct categories depending on whether they had recurred/metastasized within the follow-up period. DNA extracted from FFPE materials was subjected to Affymetrix OncoScan™ FFPE Express molecular inversion probe microarray platform for analysis of copy number changes and mutational status. Results were cross validated with Sanger sequencing, FISH and immunohistochemistry. A higher number of chromosomal aberrations were observed in cases which recurred/metastasized, with median events of 19 compared to 3.5 in cases which did not recur/metastasize. High-level amplification and homozygous deletions were detected exclusively in the former group. Regions of high-level amplification included MDM4 (1q32.1), RAF1 (3p25), EGFR (7p12) and PDZD2 (5p13.3). EGFR amplification was confirmed on FISH and accompanied by intense EGFR immunostaining. Regions of homozygous deletion included CDKN2A (9p21) and MACROD2 (20p12.1). Homozygous deletion of 9p21 which involved CDKN2A was accompanied by loss of protein expression. No mutations were identified in all samples. These findings provide insights into identifying target genes and pathways exploited by phyllodes tumours, which would aid future development of individualised therapy.

Impact of chromosomal instability on colorectal cancer progression and outcome

Background

It remains presently unclear whether disease progression in colorectal carcinoma (CRC), from early, to invasive and metastatic forms, is associated to a gradual increase in genetic instability and to a scheme of sequentially occurring Copy Number Alterations (CNAs).

Methods

In this work we set to determine the existence of such links between CRC progression and genetic instability and searched for associations with patient outcome. To this aim we analyzed a set of 162 Chromosomal Instable (CIN) CRCs comprising 131 primary carcinomas evenly distributed through stage 1 to 4, 31 metastases and 14 adenomas by array-CGH. CNA profiles were established according to disease stage and compared. We, also, asked whether the level of genomic instability was correlated to disease outcome in stage 2 and 3 CRCs. Two metrics of chromosomal instability were used; (i) Global Genomic Index (GGI), corresponding to the fraction of the genome involved in CNA, (ii) number of breakpoints (nbBP).

Results

Stage 1, 2, 3 and 4 tumors did not differ significantly at the level of their CNA profiles precluding the conventional definition of a progression scheme based on increasing levels of genetic instability. Combining GGI and nbBP,we classified genomic profiles into 5 groups presenting distinct patterns of chromosomal instability and defined two risk classes of tumors, showing strong differences in outcome and hazard risk (RFS: p = 0.012, HR = 3; OS: p < 0.001, HR = 9.7). While tumors of the high risk group were characterized by frequent fractional CNAs, low risk tumors presented predominantly whole chromosomal arm CNAs. Searching for CNAs correlating with negative outcome we found that losses at 16p13.3 and 19q13.3 observed in 10% (7/72) of stage 2–3 tumors showed strong association with early relapse (p < 0.001) and death (p < 0.007, p < 0.016). Both events showed frequent co-occurrence (p < 1x10-8) and could, therefore, mark for stage 2–3 CRC susceptible to negative outcome.

Conclusions

Our data show that CRC disease progression from stage 1 to stage 4 is not paralleled by increased levels of genetic instability. However, they suggest that stage 2–3 CRC with elevated genetic instability and particularly profiles with fractional CNA represent a subset of aggressive tumors.

Breast tumors with elevated expression of 1q candidate genes confer poor clinical outcome and sensitivity to Ras/PI3K inhibition

Genomic aberrations are common in cancers and the long arm of chromosome 1 is known for its frequent amplifications in breast cancer. However, the key candidate genes of 1q, and their contribution in breast cancer pathogenesis remain unexplored. We have analyzed the gene expression profiles of 1635 breast tumor samples using meta-analysis based approach and identified clinically significant candidates from chromosome 1q. Seven candidate genes including exonuclease 1 (EXO1) are consistently over expressed in breast tumors, specifically in high grade and aggressive breast tumors with poor clinical outcome. We derived a EXO1 co-expression module from the mRNA profiles of breast tumors which comprises 1q candidate genes and their co-expressed genes. By integrative functional genomics investigation, we identified the involvement of EGFR, RAS, PI3K / AKT, MYC, E2F signaling in the regulation of these selected 1q genes in breast tumors and breast cancer cell lines. Expression of EXO1 module was found as indicative of elevated cell proliferation, genomic instability, activated RAS/AKT/MYC/E2F1 signaling pathways and loss of p53 activity in breast tumors. mRNA-drug connectivity analysis indicates inhibition of RAS/PI3K as a possible targeted therapeutic approach for the patients with activated EXO1 module in breast tumors. Thus, we identified seven 1q candidate genes strongly associated with the poor survival of breast cancer patients and identified the possibility of targeting them with EGFR/RAS/PI3K inhibitors.

Elevated transcription of the gene QSOX1 encoding quiescin Q6 sulfhydryl oxidase 1 in breast cancer

The q arm of chromosome 1 is frequently amplified at the gene level in breast cancer. Since the significance of this is unclear we investigated whether 1q genes are overexpressed in this disease. The cDNA levels of 1q-located genes were analysed in a search for overexpressed genes. 26 genes mapping to the 1q arm show highly significant (P≤0.01) overexpression of transcripts in breast cancer compared to normal breast tissue. Amongst those showing the highest levels of overexpression in both expressed sequence tag (EST) and serial analysis of gene expression (SAGE) databases was enzyme quiescin Q6 sulfhydryl oxidase 1 (QSOX1). We investigated QSOX1 cDNA derived from T47D breast carcinoma cells by RT-PCR and 3′-RACE PCR and identified a novel extended form of QSOX1 transcript, containing a long 3′UTR, nearly double the size of the previously reported QSOX1 cDNA, and confirmed its 3′ end nucleotide sequence using RACE-PCR. We also used quantitative real-time PCR to analyse a panel of cDNAs derived from 50 clinically-graded normal and malignant breast tissue samples for the expression of QSOX1 mRNAs. QSOX1 transcription was elevated in an increasing proportion in the grade 2 and grade 3 tumours (graded according to the Nottingham prognostic index), with 10 of the 15 grade 3 tumours (67%) examined exceeding the normal range. There was a significant correlation between relative transcript level and clinical grade (P≤0.01) for all qPCR primer sets tested. QSOX1 mRNA levels, based on SAGE expression data, did not correlate with either Estrogen Receptor (ER) or Epidermal Growth Factor Receptor 2 (ErbB-2 or HER2/neu) expression. Our data indicate that QSOX1 is a potential new prognostic marker which may prove of use in the staging of breast tumours and the stratification of breast cancer patients.

Frequent copy number gains at 1q21 and 1q32 are associated with overexpression of the ETS transcription factors ETV3 and ELF3 in breast cancer irrespective of molecular subtypes

Several ETS transcription factors are involved in the pathogenesis of human cancers by different mechanisms. As gene copy number gain/amplification is an alternative mechanism of oncogenic activation and 1q gain is the most common copy number change in breast carcinoma, we investigated how that genomic change impacts in the expression of the three 1q ETS family members ETV3, ELK4, and ELF3. We have first evaluated 141 breast carcinomas for genome-wide copy number changes by chromosomal CGH and showed that 1q21 and 1q32 were the two chromosome bands with most frequent genomic copy number gains. Second, we confirmed by FISH with locus-specific BAC clones that cases showing 1q gain/amplification by CGH showed copy number increase of the ETS genes ETV3 (located in 1q21~23), ELF3, and ELK4 (both in 1q32). Third, gene expression levels of the three 1q ETS genes, as well as their potential targets MYC and CRISP3, were evaluated by quantitative real-time PCR. We here show for the first time that the most common genomic copy number gains in breast cancer, 1q21 and 1q32, are associated with overexpression of the ETS transcription factors ETV3 and ELF3 (but not ELK4) at these loci irrespective of molecular subtypes. Among the three 1q ETS genes, ELF3 has a relevant role in breast carcinogenesis and is also the most likely target of the 1q copy number increase. The basal-like molecular subtype presented the worst prognosis regarding disease-specific survival, but no additional prognostic value was found for 1q copy number status or ELF3 expression. In addition, we show that there is a correlation between the expression of the oncogene MYC, irrespectively of copy number gain at its loci in 8q24, and the expression of both the transcriptional repressor ETV3 and the androgen respondent ELK4.

Deubiquitination of EGFR by Cezanne-1 contributes to cancer progression

Once stimulated, the epidermal growth factor receptor (EGFR) undergoes self-phosphorylation, which, on the one hand, instigates signaling cascades, and on the other hand, recruits CBL ubiquitin ligases, which mark EGFRs for degradation. Using RNA interference screens, we identified a deubiquitinating enzyme, Cezanne-1, that opposes receptor degradation and enhances EGFR signaling. These functions require the catalytic and ubiquitin-binding domains of Cezanne-1, and they involve physical interactions and trans-phosphorylaton of Cezanne-1 by EGFR. In line with the ability of Cezanne-1 to augment EGF-induced growth and migration signals, the enzyme is overexpressed in breast cancer. Congruently, the corresponding gene is amplified in approximately one third of mammary tumors, and high transcript levels predict an aggressive disease course. In conclusion, deubiquitination by Cezanne-1 curtails degradation of growth factor receptors, thereby promotes oncogenic growth signals.

Haploinsufficiency in the prometastasis Kiss1 receptor Gpr54 delays breast tumor initiation, progression, and lung metastasis

Activation of KISS1 receptor (KISS1R or GPR54) by its ligands (Kisspeptins) regulates a diverse function both in normal physiology and pathophysiology. In cancer, KISS1R has been implicated in tumor angiogenesis and metastasis, but a broader evaluation of KISS1R in tumorigenesis and tumor progression is yet to be conducted. In this study, we used mouse models of Kiss1r gene knockout and mouse mammary tumor virus-polyoma virus middle T antigen (MMTV-PyMT)-induced breast cancer to conduct such an evaluation. Kiss1r heterozygosity in MMTV-PyMT mice was sufficient to attenuate breast cancer initiation, growth, latency, multiplicity, and lung metastasis. To confirm these effects and assess possible contributions of endogenous ligands, we isolated primary tumor cells from PyMT/Kiss1r(+/+) and PyMT/Kiss1r(+/-) mice and compared their phenotypes by in vitro and in vivo assays. Kiss1r loss attenuated in vitro tumorigenic properties as well as tumor growth in vivo in immunocompromised NOD.SCID/NCr mice. Kiss1r activation in these cells, resulting from the addition of its ligand Kisspeptin-10, resulted in RhoA activation and RhoA-dependent gene expression through the Gαq-p63RhoGEF signaling pathway. Anchorage-independent growth was tightly linked to dose-dependent regulation of RhoA by Kiss1r. In support of these results, siRNA-mediated knockdown of KISS1R or inactivation of RhoA in human MCF10A breast epithelial cells overexpressing H-RasV12 was sufficient to reduce Ras-induced anchorage-independent growth. In summary, we concluded that Kiss1r attenuation was sufficient to delay breast tumor initiation, progression, and metastasis through inhibitory effects on the downstream Gαq-p63RhoGEF-RhoA signaling pathway.

Prognostic significance of TRIM24/TIF-1α gene expression in breast cancer

In this study, we have analyzed the expression of TRIM24/TIF-1α, a negative regulator of various transcription factors (including nuclear receptors and p53) at the genomic, mRNA, and protein levels in human breast tumors. In breast cancer biopsy specimens, TRIM24/TIF-1α mRNA levels (assessed by Real-Time Quantitative PCR or microarray expression profiling) were increased as compared to normal breast tissues. At the genomic level, array comparative genomic hybridization analysis showed that the TRIM24/TIF-1α locus (7q34) exhibited both gains and losses that correlated with mRNA levels. By re-analyzing a series of 238 tumors, high levels of TRIM24/TIF-1α mRNA significantly correlated with various markers of poor prognosis (such as the molecular subtype) and were associated with worse overall survival. By using a rabbit polyclonal antibody for immunochemistry, the TRIM24/TIF-1α protein was detected in nuclei of normal luminal epithelial breast cells, but not in myoepithelial cells. Tissue microarray analysis confirmed that its expression was increased in epithelial cells from normal to breast infiltrating duct carcinoma and correlated with worse overall survival. Altogether, this work is the first study that shows that overexpression of the TRIM24/TIF-1α gene in breast cancer is associated with poor prognosis and worse survival, and it suggests that this transcription coregulator may play a role in mammary carcinogenesis and represent a novel prognostic marker.