Comprehensive copy number profiles of breast cancer cell model genomes

Caskin2 is a novel talin- and Abi1-binding protein that promotes cell motility

Talin (herein referring collectively to talin 1 and 2) couples the actomyosin cytoskeleton to integrins and transmits tension to the extracellular matrix. Talin also interacts with numerous additional proteins capable of modulating the actin-integrin linkage and thus downstream mechanosignaling cascades. Here, we demonstrate that the scaffold protein Caskin2 interacts directly with the R8 domain of talin through its C-terminal LD motif. Caskin2 also associates with the WAVE regulatory complex to promote cell migration in an Abi1-dependent manner. Furthermore, we demonstrate that the Caskin2-Abi1 interaction is regulated by growth factor-induced phosphorylation of Caskin2 on serine 878. In MCF7 and UACC893 cells, which contain an amplification of CASKIN2, Caskin2 localizes in plasma membrane-associated plaques and around focal adhesions in cortical microtubule stabilization complexes. Taken together, our results identify Caskin2 as a novel talin-binding protein that might not only connect integrin-mediated adhesion to actin polymerization but could also play a role in crosstalk between integrins and microtubules.

Functional analysis of a putative HER2-associated expressed enhancer, Her2-Enhancer1, in breast cancer cells

HER-2/neu (HER2) is a member of the epidermal growth factor receptors family, encoding a protein with tyrosine kinase activity. Following the gene amplification or increased HER2 transcription, carcinogenesis has been observed in some cancers. Genetic and epigenetic changes occurring in enhancer sequences can deeply affect the expression and transcriptional regulation of downstream genes, which can cause some physiological and pathological changes, including tumor progression. A therapeutic approach that directly targets the genomic sequence alterations is of high importance, with low side effects on healthy cells. Here, we employed the CRISPR/Cas9 method to genetically knockout an expressed putative enhancer (GH17J039694; we coined it as Her2-Enhancer1) located within the HER2 gene, 17q12: 39,694,339-39,697,219 (UCSC-hg38). We then investigated the potential regulatory effect of Her2-Enhancer1 on HER2 and HER2-interacting genes. To evaluate the cis and trans effects of Her2-Enhancer1, genetic manipulation of this region was performed in HER2-positive and -negative breast cancer cells. Our bioinformatics and real-time PCR data revealed that this putative enhancer region is indeed expressed, and acts as an expressed enhancer. Further functional analysis on edited and unedited cells revealed a significant alteration in the expression of HER2 variants, as well as some other target genes of HER2. Moreover, the apoptosis rate was considerably elevated within the edited cells. As we expected, Western blot analysis confirmed a reduction in protein levels of HER2, GRB7, the gene interacting with HER2, and P-AKT in the PI3K/AKT pathway. Altogether, our findings revealed an enhancer regulatory role for Her2-Enhancer1 on HER2 and HER2-interacting genes; and that this region has a potential for targeted therapy of HER2-positive cancers.

Integrative multi-omic sequencing reveals the MMTV-Myc mouse model mimics human breast cancer heterogeneity

BACKGROUND

Breast cancer is a complex and heterogeneous disease with distinct subtypes and molecular profiles corresponding to different clinical outcomes. Mouse models of breast cancer are widely used, but their relevance in capturing the heterogeneity of human disease is unclear. Previous studies have shown the heterogeneity at the gene expression level for the MMTV-Myc model, but have only speculated on the underlying genetics.

METHODS

Tumors from the microacinar, squamous, and EMT histological subtypes of the MMTV-Myc mouse model of breast cancer underwent whole genome sequencing. The genomic data obtained were then integrated with previously obtained matched sample gene expression data and extended to additional samples of each histological subtype, totaling 42 gene expression samples. High correlation was observed between genetic copy number events and resulting gene expression by both Spearman's rank correlation coefficient and the Kendall rank correlation coefficient. These same genetic events are conserved in humans and are indicative of poor overall survival by Kaplan-Meier analysis. A supervised machine learning algorithm trained on METABRIC gene expression data was used to predict the analogous human breast cancer intrinsic subtype from mouse gene expression data.

RESULTS

Herein, we examine three common histological subtypes of the MMTV-Myc model through whole genome sequencing and have integrated these results with gene expression data. Significantly, key genomic alterations driving cell signaling pathways were well conserved within histological subtypes. Genomic changes included frequent, co-occurring mutations in KIT and RARA in the microacinar histological subtype as well as SCRIB mutations in the EMT subtype. EMT tumors additionally displayed strong KRAS activation signatures downstream of genetic activating events primarily ascribed to KRAS activating mutations, but also FGFR2 amplification. Analogous genetic events in human breast cancer showed stark decreases in overall survival. In further analyzing transcriptional heterogeneity of the MMTV-Myc model, we report a supervised machine learning model that classifies MMTV-Myc histological subtypes and other mouse models as being representative of different human intrinsic breast cancer subtypes.

CONCLUSIONS

We conclude the well-established MMTV-Myc mouse model presents further opportunities for investigation of human breast cancer heterogeneity.

Characteristics of MYC Amplification and Their Association with Clinicopathological and Molecular Factors in Patients with Breast Cancer

MYC amplification is detected in ∼15% of breast tumors and is associated with poor prognosis by mediating acquired resistance to anticancer therapies. This study aimed to determine the prevalence of MYC amplifications in Chinese women with breast cancer (BRCA) and investigate the correlation between MYC amplification and clinicopathological and molecular characteristics and its clinical implications. We analyzed MYC alterations in tissue specimens from 410 women diagnosed with BRCA in our hospital from June 1, 2017 to September 27, 2018. We compared our results with publicly available data from The Cancer Genome Atlas (TCGA) BRCA cohort (n = 1079). MYC amplification was identified in 12.4% (51/410) of our cohort, with mean copy number (CN) of 4.42 (range: 2.84-11.27). In TCGA cohort, MYC amplification was identified in 21.2% (229/1079) and was associated with age, estrogen receptor status, progesterone receptor status, human epidermal growth factor receptor 2 (HER2) status, and molecular subtype, whereas in our cohort, MYC amplification was associated with smaller tumor size (T1-2, p = 0.023) and higher Ki-67 levels (≥20%; p = 0.031). Analysis of molecular profiles revealed that MYC-amplified breast tumors had significantly more concurrent CN variations compared with MYC nonamplified BRCA in both Guangdong Provincial People's Hospital (GDPH) and TCGA cohorts (p < 0.001). Pathway mapping analysis demonstrated that MYC-amplified tumors had more mutations involved in 15 different but interrelated pathways critical in DNA repair, cell cycle, and cell proliferation. Patients in TCGA cohort with MYC-amplified hormone receptor (HR)-positive/HER2-positive BRCA (p = 0.038) and MYC nonamplified triple-negative BRCA (p = 0.027) had significantly shorter overall survival. In conclusion, this study contributes to a better understanding that MYC-amplified breast tumors had distinct clinicopathological and molecular features compared with MYC nonamplified breast tumors. Further research with a larger sample size is necessary to further elucidate the clinical and survival implications of MYC amplifications.

A pipeline for copy number profiling of single circulating tumor cells to assess intra-patient tumor heterogeneity

Intrapatient tumour heterogeneity is likely a major determinant of clinical outcome in cancer patients. To assess heterogeneity in a minimally invasive manner, methods to perform single circulating tumour cell (CTC) genomics at high resolution are necessary. However, due to the rarity of CTCs, development of such methods is challenging. Here, we developed a modular single CTC analysis pipeline to assess intrapatient heterogeneity by copy number (CN) profiling. To optimize this pipeline, spike‐in experiments using MCF‐7 breast cancer cells were performed. The VyCAP puncher system was used to isolate single cells. The quality of whole genome amplification (WGA) products generated by REPLI‐g and Ampli1™ methods, as well as the results from the Illumina Truseq and the Ampli1™ LowPass library preparation techniques, was compared. Moreover, a bioinformatic pipeline was designed to generate CN profiles from single CTCs. The optimal combination of Ampli1™ WGA and Illumina Truseq library preparation was successfully validated on patient‐derived CTCs. In conclusion, we developed a novel modular pipeline to isolate single CTCs and subsequently generate detailed patient‐derived CN profiles that allow assessment of intrapatient heterogeneity in future studies.

The mechanisms involved in the resistance of estrogen receptor-positive breast cancer cells to palbociclib are multiple and change over time

PURPOSE

Cyclin-dependent kinase 4 and 6 (CDK4/6) inhibitors are widely used for the treatment of advanced estrogen receptor (ER)-positive breast cancer. To develop a treatment strategy for cancers resistant to CDK4/6 inhibitors, here, we established palbociclib-resistant sublines and analyzed their resistance mechanisms.

METHODS

Palbociclib-resistant sublines were established from T47D and MCF7 cells. Sensitivity to other drugs was assessed via the WST assay. Altered expression/phosphorylation of proteins related to signal transduction and cell cycle regulation was examined using western blotting. Copy number alterations and mutations in the retinoblastoma (RB1) gene were also analyzed.

RESULTS

Although an increase in CDK6 and decrease in retinoblastoma protein (Rb) expression/phosphorylation were commonly observed in the resistant sublines, changes in other cell cycle-related proteins were heterogeneous. Upon extended exposure to palbociclib, the expression/phosphorylation of these proteins became altered, and the long-term removal of palbociclib did not restore the Rb expression/phosphorylation patterns. Consistently a copy number decrease, as well as RB1 mutations were detected. Moreover, although the resistant sublines exhibited cross-resistance to abemaciclib, their response to dinaciclib was the same as that of wild-type cells. Of note, the cell line exhibiting increased mTOR phosphorylation also showed a higher sensitivity to everolimus. However, the sensitivity to chemotherapeutic agents was unchanged in palbociclib-resistant sublines.

CONCLUSION

ER-positive breast cancer cells use multiple molecular mechanisms to survive in the presence of palbociclib, suggesting that targeting activated proteins may be an effective strategy to overcome resistance. Additionally, palbociclib monotherapy induces mutations and copy number alterations in the RB1 gene.

Gli activation by the estrogen receptor in breast cancer cells: Regulation of cancer cell growth by Gli3

BACKGROUND

Gli is an oncogenic transcription factor family thought to be involved in breast cancer (BrCa) cell growth. Gli activity is regulated by a post-translational proteolytic process that is suppressed by Hedgehog signaling. In prostate cancer cells, however, Gli activation is mediated by an interaction of active androgen receptor proteins with Gli3 that stabilizes Gli3 in its un-proteolyzed form. Here we show that the estrogen receptor (ER), ERα, also binds Gli3 and activates Gli in BrCa cells. Moreover, we show that ER + BrCa cells are dependent on Gli3 for cancer cell growth.

METHODS

Transfection with Gli-luciferase reporter was used to report Gli activity in 293FT or BrCa cells (MCF7, T47D, MDA-MB-453) with or without steroid ligands. Co-immunoprecipitation and proximity ligation were used to show association of Gli3 with ERα. Gli3 stability was determined by western blots of BrCa cell extracts. ERα knockdown or destabilization (by fulvestrant) was used to assess how loss of ERα affects estradiol-induced Gli reporter activity, formation of intranuclear ERα-Gli3 complexes and Gli3 stability. Expression of Gli1 and/or other endogenous Gli-target genes in BrCa cells were measured by qPCR in the presence or absence of estradiol. Gli3 knockdown was assessed for effects on BrCa cell growth using the Cyquant assay.

RESULTS

ERα co-transfection increased Gli reporter activity in 293FT cells that was further increased by estradiol. Gli3 co-precipitated in ERα immunoprecipitates. Acute (2 h) estradiol increased Gli reporter activity and the formation of intranuclear ERα-Gli3 complexes in ER + BrCa cells but more chronic estradiol (48 h) reduced ERα-Gli complexes commensurate with reduced ERα levels. Gli3 stability and endogenous activity was only increased by more chronic estradiol treatment. Fulvestrant or ERα knockdown suppressed E2-induction of Gli activity, intranuclear ERα-Gli3 complexes and stabilization of Gli3. Gli3 knockdown significantly reduced the growth of BrCa cells.

CONCLUSIONS

ERα interacts with Gli3 in BrCa cells and estradiol treatment leads to Gli3 stabilization and increased expression of Gli-target genes. Furthermore, we found tthat Gli3 is necessary for BrCa cell growth. These results support the idea that the ERα-Gli interaction and Gli3 may be novel targets for effective control of BrCa growth.

The transcriptomic response of cells to a drug combination is more than the sum of the responses to the monotherapies

Our ability to discover effective drug combinations is limited, in part by insufficient understanding of how the transcriptional response of two monotherapies results in that of their combination. We analyzed matched time course RNAseq profiling of cells treated with single drugs and their combinations and found that the transcriptional signature of the synergistic combination was unique relative to that of either constituent monotherapy. The sequential activation of transcription factors in time in the gene regulatory network was implicated. The nature of this transcriptional cascade suggests that drug synergy may ensue when the transcriptional responses elicited by two unrelated individual drugs are correlated. We used these results as the basis of a simple prediction algorithm attaining an AUROC of 0.77 in the prediction of synergistic drug combinations in an independent dataset.

A Pilot Study for the Feasibility of Exome-Sequencing in Circulating Tumor Cells Versus Single Metastatic Biopsies in Breast Cancer

The comparison of the landscape of somatic alterations in circulating tumor cells (CTCs) versus metastases is challenging. Here, we comprehensively characterized the somatic landscape in bulk (amplified and non-amplified), spike-in breast cancer cells, CTCs, and metastases from breast cancer patients using whole-exome sequencing (WES). We determined the level of genomic concordance for somatic nucleotide variants (SNVs), copy number alterations (CNAs), and structural variants (SVs). The variant allele fractions (VAFs) of somatic variants were remarkably similar between amplified and non-amplified cell line samples as technical replicates. In clinical samples, a significant fraction of somatic variants had low VAFs in CTCs compared to metastases. The most frequently recurrent gene mutations in clinical samples were associated with an elevated C > T mutational signature. We found complex rearrangement patterns including intra- and inter-chromosomal rearrangements, singleton, and recurrent gene fusions, and tandem duplications. We observed high molecular discordance for somatic alterations between paired samples consistent with marked heterogeneity of the somatic landscape. The most prevalent copy number calls were focal deletion events in CTCs and metastases. Our results demonstrate the feasibility of an integrated workflow for the identification of a complete repertoire of somatic alterations and highlight the intrapatient genomic differences that occur between CTCs and metastases.

Evolution of a Landscape Phage Library in a Mouse Xenograft Model of Human Breast Cancer

Peptide-displayed phage libraries are billion-clone collections of diverse chimeric bacteriophage particles, decorated by genetically fused peptides built from a random combination of natural amino acids. Studying the molecular evolution of peptide-displayed libraries in mammalian model systems, using in vivo phage display techniques, can provide invaluable knowledge about the underlying physiology of the vasculature system, allow recognition of organ- and tissue-specific networks of protein-protein interactions, and provide ligands for targeted diagnostics and therapeutics. Recently, we discovered that landscape phage libraries, a specific type of multivalent peptide phage display library, expose on their surface comprehensive collections of elementary binding units (EBUs), which can form short linear motifs (SLiMs) that interact with functional domains of physiologically relevant proteins. Because of their unique structural and functional features, landscape phages can use an alternative mechanism of directed molecular evolution, i.e., combinatorial avidity selection. These discoveries fueled our interest in revisiting the in vivo evolution of phage displayed libraries using another format of display, i.e., landscape phages. In this study, we monitored the evolution of a landscape phage library in a mouse model with and without an implanted human breast cancer tumor xenograft. As expected, the multivalent architecture of landscape phage displayed proteins provided strong tissue selectivity and resulted in a huge diversity of tissue penetrating, chimeric phage particles. We identified several types of EBU interactions that evolved during the course of tissue distribution, which included interactions of EBUs with all tissue types, those EBUs that interacted selectively with specific organs or tissues with shared gene expression profiles or functionalities, and other EBUs that interacted in a tissue-selective manner. We demonstrated that landscape phage libraries are a rich collection of unique nanobioparticles that can be used to identify functional organ and tissue-binding elements after the evolution of a phage display library in vivo.

Activation of HIF-1α by δ-Opioid Receptors Induces COX-2 Expression in Breast Cancer Cells and Leads to Paracrine Activation of Vascular Endothelial Cells

Opioids promote tumor angiogenesis in mammary malignancies, but the underlying signaling mechanism is largely unknown. The current study investigated the hypothesis that stimulation of δ-opioid receptors (DOR) in breast cancer (BCa) cells activates the hypoxia-inducible factor 1α (HIF-1α), which triggers synthesis and release of diverse angiogenic factors. Immunoblotting revealed that incubation of human MCF-7 and T47D breast cancer cells with the DOR agonist d-Ala²,d-Leu⁵-enkephalin (DADLE) resulted in a transient accumulation and thus activation of HIF-1α DADLE-induced HIF-1α activation preceded PI3K/Akt stimulation and was blocked by the DOR antagonist naltrindole and naloxone, pertussis toxin, different phosphoinositide 3-kinase (PI3K) inhibitors, and the Akt inhibitor Akti-1/2. Whereas DADLE exposure had no effect on the expression and secretion of vascular endothelial growth factor (VEGF) in BCa cells, an increased abundance of cyclooxygenase-2 (COX-2) and release of prostaglandin E2 (PGE₂) was detected. DADLE-induced COX-2 expression was also observed in three-dimensional cultured MCF-7 cells and impaired by PI3K/Akt inhibitors and the HIF-1α inhibitor echinomycin. Supernatant from DADLE-treated MCF-7 cells triggered sprouting of endothelial (END) cells, which was blocked when MCF-7 cells were pretreated with echinomycin or the COX-2 inhibitor celecoxib. Also no sprouting was observed when END cells were exposed to the PGE₂ receptor antagonist PF-04418948. The findings together indicate that DOR stimulation in BCa cells leads to PI3K/Akt-dependent HIF-1α activation and COX-2 expression, which trigger END cell sprouting by paracrine activation of PGE₂ receptors. These findings provide a potential mechanism of opioid-driven tumor angiogenesis and thus therapeutic targets to combat the tumor-angiogenic opioid effect. SIGNIFICANCE STATEMENT: Opioids are indispensable analgesics for treating cancer-related pain. However, opioids were found to promote tumor growth and metastasis, which questions the use of these potent pain-relieving drugs in cancer patients. Enhanced tumor vascularization after opioid treatment implies that tumor progression results from angiogenic opioid effects. Thus, understanding the signaling mechanism of opioid-driven tumor angiogenesis helps to identify therapeutic targets to combat these undesired tumor effects. The present study reveals that stimulation of δ-opioid receptors in breast cancer cells leads to an activation of HIF-1α and expression of COX-2 via PI3K/Akt stimulation, which results in a paracrine activation of vascular endothelial cells by prostaglandin E₂ receptors.

SALMON: Survival Analysis Learning With Multi-Omics Neural Networks on Breast Cancer

Improved cancer prognosis is a central goal for precision health medicine. Though many models can predict differential survival from data, there is a strong need for sophisticated algorithms that can aggregate and filter relevant predictors from increasingly complex data inputs. In turn, these models should provide deeper insight into which types of data are most relevant to improve prognosis. Deep Learning-based neural networks offer a potential solution for both problems because they are highly flexible and account for data complexity in a non-linear fashion. In this study, we implement Deep Learning-based networks to determine how gene expression data predicts Cox regression survival in breast cancer. We accomplish this through an algorithm called SALMON (Survival Analysis Learning with Multi-Omics Neural Networks), which aggregates and simplifies gene expression data and cancer biomarkers to enable prognosis prediction. The results revealed improved performance when more omics data were used in model construction. Rather than use raw gene expression values as model inputs, we innovatively use eigengene modules from the result of gene co-expression network analysis. The corresponding high impact co-expression modules and other omics data are identified by feature selection technique, then examined by conducting enrichment analysis and exploiting biological functions, escalated the interpretation of input feature from gene level to co-expression modules level. Our study shows the feasibility of discovering breast cancer related co-expression modules, sketch a blueprint of future endeavors on Deep Learning-based survival analysis. SALMON source code is available at https://github.com/huangzhii/SALMON/.

Clinical Importance of Myc Family Oncogene Aberrations in Epithelial Ovarian Cancer

Background

The Myc oncogene family has been implicated in many human malignancies and is often associated with particularly aggressive disease, suggesting Myc as an attractive prognostic marker and therapeutic target. However, for epithelial ovarian cancer (EOC), there is little consensus on the incidence and clinical relevance of Myc aberrations. Here we comprehensively investigated alterations in gene copy number, expression, and activity for Myc and evaluated their clinical significance in EOC.

Methods

To address inconsistencies in the literature regarding the definition of copy number variations, we developed a novel approach using quantitative polymerase chain reaction (qPCR) coupled with a statistical algorithm to estimate objective thresholds for detecting Myc gain/amplification in large cohorts of serous (n = 150) and endometrioid (n = 80) EOC. MYC, MYCN, and MYCL1 mRNA expression and Myc activity score for each case were examined by qPCR. Kaplan–Meier and Cox-regression analyses were conducted to assess clinical significance of Myc aberrations.

Results

Using a large panel of cancer cell lines (n = 34), we validated the statistical algorithm for determining clear thresholds for Myc gain/amplification. MYC was the most predominantly amplified of the Myc oncogene family members, and high MYC mRNA expression levels were associated with amplification in EOC. However, there was no association between prognosis and increased copy number or gene expression of MYC/MYCN/MYCL1 or with a pan-Myc transcriptional activity score, in EOC, although MYC amplification was associated with late stage and high grade in endometrioid EOC.

Conclusion

A systematic and comprehensive analysis of Myc genes, transcripts, and activity levels using qPCR revealed that although such aberrations commonly occur in EOC, overall they have limited impact on outcome, suggesting that the biological relevance of Myc oncogene family members is limited to certain subsets of this disease.

Nanopipettes as Monitoring Probes for the Single Living Cell: State of the Art and Future Directions in Molecular Biology

Examining the behavior of a single cell within its natural environment is valuable for understanding both the biological processes that control the function of cells and how injury or disease lead to pathological change of their function. Single-cell analysis can reveal information regarding the causes of genetic changes, and it can contribute to studies on the molecular basis of cell transformation and proliferation. By contrast, whole tissue biopsies can only yield information on a statistical average of several processes occurring in a population of different cells. Electrowetting within a nanopipette provides a nanobiopsy platform for the extraction of cellular material from single living cells. Additionally, functionalized nanopipette sensing probes can differentiate analytes based on their size, shape or charge density, making the technology uniquely suited to sensing changes in single-cell dynamics. In this review, we highlight the potential of nanopipette technology as a non-destructive analytical tool to monitor single living cells, with particular attention to integration into applications in molecular biology.

Stress-inducible Protein-1 promotes metastasis of gastric cancer via Wnt/β-catenin signaling pathway

BACKGROUND

Stress-Inducible Protein-1 (STIP1) is a co-chaperone that associates directly with heat shock proteins, and regulates motility of various types of cancer. In the present study, we investigated the role of STIP1 on metastasis of gastric cancer (GC).

METHODS

In vivo metastatic experimental model was employed to investigate the effect of STIP1 on metastasis of GC cells. Loss-of-function and gain-of-function experiments were performed to examine the role of STIP1 on metastasis of GC cells. Western blot, immunofluorescence staining, migration and invasion assays, microarray and KEGG pathway analysis were applied to explore the underlying mechanism.

RESULTS

In current study, we demonstrated that STIP1 promoted lung metastasis of GC cells in vivo. Furthermore, STIP1 significantly enhanced migration and invasion abilities of GC cells. In contrast, knock-down of STIP1 yielded the opposite effects on these phenotypes in vitro. STIP1 promoted tumor metastasis through inducing epithelial-to-mesenchymal transition in GC cells. Mechanistically, STIP1 promoted GC metastasis via up-regulation of targeted genes in Wnt/β-catenin signaling pathway, including c-Myc and Cyclin D1, and accompanied with nuclear translocation of β-catenin.

CONCLUSIONS

Our findings indicate that elevated expression of STIP1 exhibited a metastasis-promoting effect in GC cells through activation of Wnt/β-catenin signaling pathway. STIP1 may be served as a potential therapeutic target for preventing GC metastasis.

Transcriptional signature of lymphoblastoid cell lines of BRCA1, BRCA2 and non-BRCA1/2 high risk breast cancer families

Approximately 25% of hereditary breast cancer cases are associated with a strong familial history which can be explained by mutations in BRCA1 or BRCA2 and other lower penetrance genes. The remaining high-risk families could be classified as BRCAX (non-BRCA1/2) families.

Gene expression involving alternative splicing represents a well-known mechanism regulating the expression of multiple transcripts, which could be involved in cancer development. Thus using RNA-seq methodology, the analysis of transcriptome was undertaken to potentially reveal transcripts implicated in breast cancer susceptibility and development.

RNA was extracted from immortalized lymphoblastoid cell lines of 117 women (affected and unaffected) coming from BRCA1, BRCA2 and BRCAX families. Anova analysis revealed a total of 95 transcripts corresponding to 85 different genes differentially expressed (Bonferroni corrected p-value <0.01) between those groups. Hierarchical clustering allowed distinctive subgrouping of BRCA1/2 subgroups from BRCAX individuals. We found 67 transcripts, which could discriminate BRCAX from BRCA1/BRCA2 individuals while 28 transcripts discriminate affected from unaffected BRCAX individuals.

To our knowledge, this represents the first study identifying transcripts differentially expressed in lymphoblastoid cell lines from major classes of mutation-related breast cancer subgroups, namely BRCA1, BRCA2 and BRCAX. Moreover, some transcripts could discriminate affected from unaffected BRCAX individuals, which could represent potential therapeutic targets for breast cancer treatment.

Autocrine and paracrine STIP1 signaling promote osteolytic bone metastasis in renal cell carcinoma

Bone metastases are responsible for some of the most devastating complications of renal cell carcinoma (RCC). However, pro-metastatic factors leading to the highly osteolytic characteristics of RCC bone metastasis have barely been explored. We previously developed novel bone-seeking RCC cell lines by the in vivo selection strategy and performed a comparative proteome analysis on their total cell lysate. Here, we focused on STIP1 (stress-induced phosphoprotein 1), the high up-regulated protein in the bone-seeking cells, and explored its clinical relevance and functions in RCC bone metastasis. We observed high levels of both intracellular and extracellular STIP1 protein in bone metastatic tissue samples. Elevated STIP1 mRNA in the primary RCC tumors remarkably correlated with worse clinical outcomes. Furthermore, both human recombinant STIP1 protein and anti-STIP1 neutralizing antibody were used in the functional studies. We found that 1) STIP1 protein on the extracellular surface of tumor cells promoted the proliferation and migration/invasion of RCC tumor cells through the autocrine STIP1-ALK2-SMAD1/5 pathway; and 2) STIP1 protein secreted into the extracellular tumor stromal area, promoted the differentiation of osteoclasts through the paracrine STIP1-PrPc-ERK1/2 pathway. Increased cathepsin K (CTSK), the key enzyme secreted by osteoclasts to degrade collagen and other matrix proteins during bone resorption was further detected in the differentiated osteoclasts. These results provide evidence of the great potential of STIP1 as a novel biomarker and therapeutic target in RCC bone metastasis.

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.

Anti-breast cancer effects of live, heat-killed and cytoplasmic fractions of Enterococcus faecalis and Staphylococcus hominis isolated from human breast milk

Development of tumour that is resistant to chemotherapeutics and synthetic drugs, coupled with their life-threatening side effects and the adverse effects of surgery and hormone therapies, led to increased research on probiotics' anticancer potentials. The current study investigated the potential of live, heat-killed cells (HKC) and the cytoplasmic fractions (CF) of Enterococcus faecalis and Staphylococcus hominis as anti-breast cancer agents. MCF-7 cell line was treated with 25, 50, 100 and 200 μg/mL each of live, HKC and CF of the bacteria; and cytotoxicity was evaluated for 24, 48 and 72 h using MTT assay. The morphological features of the treated cells were examined by fluorescence microscopy. The stage of cell cycle arrest and apoptosis were quantified by flow cytometry. The bacterial effect on non-malignant breast epithelial cell line, MCF-10A, was assessed using MTT assay for 24, 48 and 72 h. All the three forms of the bacteria caused a significant decrease in MCF-7 (up to 33.29%) cell proliferation in concentration- and time-dependent manner. Morphological features of apoptosis like cell death, cell shrinkage and membrane blebbing were observed. Flow cytometry analyses suggested that about 34.60% of treated MCF-7 was undergoing apoptosis. A strong anti-proliferative activity was efficiently induced through sub-G1 accumulation (up to 83.17%) in treated MCF-7 and decreased number in the G0/G1 phase (74.39%). MCF-10A cells treated with both bacteria showed no significant difference with the untreated (>90% viability). These bacteria can be used as good alternative nutraceutical with promising therapeutic indexes for breast cancer because of their non-cytotoxic effects to normal cells.

A novel gammaretroviral shuttle vector insertional mutagenesis screen identifies SHARPIN as a breast cancer metastasis gene and prognostic biomarker

Breast cancer (BC) is the second leading cause of malignancy among U.S. women. Metastasis results in a poor prognosis and increased mortality, but the molecular mechanisms by which metastatic tumors occur are not well understood. Identifying the genes that drive the metastatic process could provide targets for improved therapy and biomarkers to improve BC patient outcomes. Using a forward mutagenesis screen, BC cells mutagenized with a replication-incompetent gammaretroviral vector (γRV) were xenotransplanted into the mammary fat pad of immunodeficient mice. In this approach the vector provirus dysregulates nearby genes, providing a selective advantage to transduced cells to form metastases. Metastatic tumors were analyzed for proviral integration sites to identify nearby candidate metastasis genes. The γRV has a transgene cassette that allows for rescue in bacteria and rapid identification of vector integration sites. Using this approach, we identified the previously described metastasis gene WWTR1 (TAZ), and three other novel candidate metastasis genes including SHARPIN. SHARPIN was independently validated in vivo as a BC metastasis gene. Analysis of patient data showed that SHARPIN expression predicts metastasis-free survival after adjuvant therapy. Our approach has broad potential to identify genes involved in oncogenic processes for BC and other cancers. We show here it can identify both known (WWTR1) and novel (SHARPIN) BC metastasis genes.

Akt Phosphorylates Wnt Coactivator and Chromatin Effector Pygo2 at Serine 48 to Antagonize Its Ubiquitin/Proteasome-mediated Degradation

Pygopus 2 (Pygo2/PYGO2) is an evolutionarily conserved coactivator and chromatin effector in the Wnt/β-catenin signaling pathway that regulates cell growth and differentiation in various normal and malignant tissues. Although PYGO2 is highly overexpressed in a number of human cancers, the molecular mechanism underlying its deregulation is largely unknown. Here we report that Pygo2 protein is degraded through the ubiquitin/proteasome pathway and is posttranslationally stabilized through phosphorylation by activated phosphatidylinositol 3-kinase/Akt signaling. Specifically, Pygo2 is stabilized upon inhibition of the proteasome, and its intracellular level is regulated by Cullin 4 (Cul4) and DNA damage-binding protein 1 (DDB1), components of the Cul4-DDB1 E3 ubiquitin ligase complex. Furthermore, Pygo2 is phosphorylated at multiple residues, and Akt-mediated phosphorylation at serine 48 leads to its decreased ubiquitylation and increased stability. Finally, we provide evidence that Akt and its upstream growth factors act in parallel with Wnt to stabilize Pygo2. Taken together, our findings highlight chromatin regulator Pygo2 as a common node downstream of oncogenic Wnt and Akt signaling pathways and underscore posttranslational modification, particularly phosphorylation and ubiquitylation, as a significant mode of regulation of Pygo2 protein expression.

Cross-species DNA copy number analyses identifies multiple 1q21-q23 subtype-specific driver genes for breast cancer

A large number of DNA copy number alterations (CNAs) exist in human breast cancers, and thus characterizing the most frequent CNAs is key to advancing therapeutics because it is likely that these regions contain breast tumor ‘drivers’ (i.e., cancer causal genes). This study aims to characterize the genomic landscape of breast cancer CNAs and identify potential subtype-specific drivers using a large set of human breast tumors and genetically engineered mouse (GEM) mammary tumors. Using a novel method called SWITCHplus, we identified subtype-specific DNA CNAs occurring at a 15 % or greater frequency, which excluded many well-known breast cancer-related drivers such as amplification of ERBB2, and deletions of TP53 and RB1. A comparison of CNAs between mouse and human breast tumors identified regions with shared subtype-specific CNAs. Additional criteria that included gene expression-to-copy number correlation, a DawnRank network analysis, and RNA interference functional studies highlighted candidate driver genes that fulfilled these multiple criteria. Numerous regions of shared CNAs were observed between human breast tumors and GEM mammary tumor models that shared similar gene expression features. Specifically, we identified chromosome 1q21-23 as a Basal-like subtype-enriched region with multiple potential driver genes including PI4KB, SHC1, and NCSTN. This step-wise computational approach based on a cross-species comparison is applicable to any tumor type for which sufficient human and model system DNA copy number data exist, and in this instance, highlights that a single region of amplification may in fact harbor multiple driver genes.

Tumor protein D52-like 2 contributes to proliferation of breast cancer cells

Breast cancer is a major cause of cancer-related death among women. Tumor protein D52-like 2 (TPD52L2) is one member of the TPD52 family, which has been shown to function in mediating cell proliferation, apoptosis, and vehicle trafficking. TPD52 was originally identified in human breast carcinoma. In this study, the authors found that TPD52L2 is extensively expressed in multiple human breast cancer cell lines. To elucidate the functional role of TPD52L2 in breast cancer, the authors employed lentivirus-mediated short hairpin RNA (shRNA) to knock down TPD52L2 in one breast cancer cell line, ZR-75-30, which showed high TPD52L2 expression. The shRNA-mediated TPD52L2 knockdown inhibited the proliferation and colony formation in ZR-75-30 cells, as determined by MTT and colony formation assays. Knockdown of TPD52L2 led to an accumulation of cells in the G0/G1 phase of the cell cycle. Furthermore, knockdown of TPD52L2 promoted GSK3β phosphorylation in ZR-75-30 cells. This investigation indicates that TPD52L2 plays an essential role in the growth of breast cancer cells, which may contribute to provide gene therapy for breast cancer treatment.

Identification of genomic signatures in circulating tumor cells from breast cancer

Levels of circulating tumor cells (CTCs) in blood have prognostic value in early and metastatic breast cancer. CTCs also show varying degrees of concordance with molecular markers of primary tumors they originate from. It is expected that individual cells reflect the heterogeneity and evolution of tumor cells as they acquire new functions and differential responses to chemotherapy. However, a degree of commonality is also plausible, highlighting alterations that allow tumor cells to perform CTC-defining activities such as invasion and intravasation. Using a matched tumor-normal approach, we performed high-resolution copy number profiling of CTCs from breast cancer to identify occult changes occurring during progression to metastasis. We identified a signature of recurrent gain in CTCs, consisting of 90 minimal common regions (MCRs) of copy number gain. These were predominantly found across chromosome 19 and were identified at low frequencies (3-4%) in 787 primary breast carcinomas examined. CTC genomic signatures clustered into two groups independent of subtype: a dormancy-related signature with 16 MCRs (AKT2, PTEN, CADM2); and a tumor-aggressiveness related signature with 358 MCRs (ANGPTL4, BSG, MIR-373). There were two MCRs in common between the groups on 19q13 and 21q21, containing genes involved in resistance to anoikis, TGFβ-signaling and metastasis (TFF3, LTBP4, NUMBL). Furthermore, a region harboring the ERBB2 gene was gained in a majority of patients. Regions 20q13 and 15q24 were associated with distant metastasis. The distinctiveness of CTC signatures highlights cell populations with different functional or metastatic potential. Such novel targets could help to specifically identify and block dissemination.

p21-activated kinases in Erbb2-positive breast cancer: A new therapeutic target?

The activation of receptor tyrosine kinases, particularly ErbB2, has been linked to the genesis and progression of breast cancer. Two of the central signaling pathways activated by ErbB2 are the Ras/Raf-1/Mek/Erk pathway, which plays an important role in tumor cell growth and migration, and the PI3K/Akt pathway, which plays an important role in cell survival. Recently, we and others have shown that signaling through the Ras-Erk pathway can be influenced by p21-activated kinase 1 (Pak1), an effector of the Rho family GTP ases Rac and Cdc42. Expression of activated forms of Rac promotes activation of Erk through mechanisms involving Pak1 phosphorylation of Raf-1 and Mek1. In addition, Pak1 has also been implicated in the activation of Akt. However, our understanding regarding the degree to which Rho GTPases, and their effectors such as Pak1, contribute to ErbB2-mediated signaling is very limited.Recent results from our laboratory indicate that ErbB2 expression correlates with Pak activation in estrogen receptor negative human breast tumor samples. Using a three-dimensional (3D) culture of human MCF-10A mammary epithelial cells, we found that activation of Rac-Pak pathway by ErbB2 induces growth factor independent proliferation and promotes disruption of acini-like structures through the activation of the Erk and Akt pathways. We also observed that blocking Pak1 activity by small molecule inhibitors impeded the ability of activated ErbB2 to transform these cells and to activate its associated downstream signaling targets. In addition, we found that suppressing Pak activity in ErbB2-amplified breast cancer cells delayed tumor formation and downregulated Erk and Akt signaling in vivo. These results support a model in which Pak, by activating Erk and Akt, cooperates with ErbB2 in transforming mammary epithelial cells.

GAP-Seq: a method for identification of DNA palindromes

BACKGROUND

Closely spaced long inverted repeats, also known as DNA palindromes, can undergo intrastrand annealing to form DNA hairpins. The ability to form these hairpins results in genome instability, difficulties in maintaining clones in Escherichia coli and major problems for most DNA sequencing approaches. Because of their role in genomic instability and gene amplification in some human cancers, it is important to develop systematic approaches to detect and characterize DNA palindromes.

RESULTS

We developed a new protocol to identify palindromes that couples the S1 nuclease treated Cot0 DNA (GAPF) with high-throughput sequencing (GAP-Seq). Unlike earlier protocols, it does not involve restriction enzymatic digestion prior to DNA snap-back thereby preserving longer DNA sequences. It also indicates the location of the novel junction, which can then be recovered. Using MCF-7 breast cancer cell line as the proof-of-principle analysis, we have identified 35 palindrome candidates and physically characterized the top 5 candidates and their junctions. Because this protocol eliminates many of the false positives that plague earlier techniques, we have improved palindrome identification.

CONCLUSIONS

The GAP-Seq approach underscores the importance of developing new tools for identifying and characterizing palindromes, and provides a new strategy to systematically assess palindromes in genomes. It will be useful for studying human cancers and other diseases associated with palindromes.

Overexpression of the HER2/neu Gene: A New Therapeutic Possibility for Patients With Advanced Gallbladder Cancer

BACKGROUND

The HER2/neu gene is a proto-oncogene that can predict the response to treatment with trastuzumab, pertuzumab, and lapatinib. This study was conducted to determine the frequency of HER2/neu overexpression and to identify a subgroup of patients with gallbladder cancer who would benefit from targeted therapy.

METHODS

Patients with gallbladder cancer (n = 187; 165 women and 22 men) with a recorded follow-up of at least 5 years were included, along with control subjects (n = 75). An automated immunohistochemical technique was used with an anti-ErbB2 antibody. Scoring was conducted according to the CAP/ASCO (College of American Pathologists/American Society of Clinical Oncology) criteria for breast cancer.

RESULTS

Overexpression of HER2/neu was observed in 12.8% of the cases. Of those, 0% were mucosal, 14.3% muscular, 12.8% subserosal, and 10.6% serosal. In 20% of the cases, equivocal staining was observed. Overexpression was more frequent in the advanced cancers and in the better differentiated tumors (13.8% and 17.4%, respectively), but the difference was nonsignificant. The patients with overexpression of HER2/neu had a worse overall survival, when compared with those who had no expression at 5 years (34% vs. 41%).

CONCLUSION

This is the single largest study of HER2/neu expression in gallbladder cancer to use commonly accepted scoring criteria. The results indicate that HER2/neu overexpression occurred in 14% of the advanced gallbladder cancer cases. This subgroup may benefit from inhibitors of the HER2/neu pathway.

Breast cancer biomarker measurements and standards

Cancer is a heterogeneous disease characterized by changes in the levels and activities of important cellular proteins, including oncogenes and tumor suppressors. Genetic mutations cause changes in protein activity and protein expression levels that result in the altered metabolism, proliferation, and metastasis seen in cancer cells. The identification of the critical biochemical changes in cancer has led to advances in its detection and treatment. An important example of this is the measurement of human epidermal growth factor receptor 2 (HER2), where increased expression occurs in approximately 20-30% of breast cancer tumors. HER2 is a member of the epidermal growth factor receptor family and is an important biomarker expressed on the cell surface. Measurement of the HER2 levels in tumor cells provides diagnostic, prognostic, and treatment information, because a targeted therapeutic is available. The most common methods to measure HER2 levels are immunohistochemistry and in situ hybridization assays. The accurate and reliable measurements of the specific changes in protein biomarkers for detection and treatment of cancer are important challenges. This review is focused on efforts to improve the quantitation and reliability of cancer biomarkers by using standards and reference materials.

Differences and homologies of chromosomal alterations within and between breast cancer cell lines: a clustering analysis

Background

The MCF7 (ER+/HER2-), T47D (ER+/HER2-), BT474 (ER+/HER2+) and SKBR3 (ER-/HER2+) breast cancer cell lines are widely used in breast cancer research as paradigms of the luminal and HER2 phenotypes. Although they have been subjected to cytogenetic analysis, their chromosomal abnormalities have not been carefully characterized, and their differential cytogenetic profiles have not yet been established. In addition, techniques such as comparative genomic hybridization (CGH), microarray-based CGH and multiplex ligation-dependent probe amplification (MLPA) have described specific regions of gains, losses and amplifications of these cell lines; however, these techniques cannot detect balanced chromosomal rearrangements (e.g., translocations or inversions) or low frequency mosaicism.

Results

A range of 19 to 26 metaphases of the MCF7, T47D, BT474 and SKBR3 cell lines was studied using conventional (G-banding) and molecular cytogenetic techniques (multi-color fluorescence in situ hybridization, M-FISH). We detected previously unreported chromosomal changes and determined the content and frequency of chromosomal markers. MCF7 and T47D (ER+/HER2-) cells showed a less complex chromosomal make up, with more numerical than structural alterations, compared to BT474 and SKBR3 (HER2+) cells, which harbored the highest frequency of numerical and structural aberrations. Karyotype heterogeneity and clonality were determined by comparing all metaphases within and between the four cell lines by hierarchical clustering. The latter analysis identified five main clusters. One of these clusters was characterized by numerical chromosomal abnormalities common to all cell lines, and the other four clusters encompassed cell-specific chromosomal abnormalities. T47D and BT474 cells shared the most chromosomal abnormalities, some of which were shared with SKBR3 cells. MCF7 cells showed a chromosomal pattern that was markedly different from those of the other cell lines.

Conclusions

Our study provides a comprehensive and specific characterization of complex chromosomal aberrations of MCF7, T47D, BT474 and SKBR3 cell lines.

The chromosomal pattern of ER+/HER2- cells is less complex than that of ER+/HER2+ and ER-/HER2+ cells. These chromosomal abnormalities could influence the biologic and pharmacologic response of cells. Finally, although gene expression profiling and aCGH studies have classified these four cell lines as luminal, our results suggest that they are heterogeneous at the cytogenetic level.

A Single-Array-Based Method for Detecting Copy Number Variants Using Affymetrix High Density SNP Arrays and its Application to Breast Cancer

Cumulative evidence has shown that structural variations, due to insertions, deletions, and inversions of DNA, may contribute considerably to the development of complex human diseases, such as breast cancer. High-throughput genotyping technologies, such as Affymetrix high density single-nucleotide polymorphism (SNP) arrays, have produced large amounts of genetic data for genome-wide SNP genotype calling and copy number estimation. Meanwhile, there is a great need for accurate and efficient statistical methods to detect copy number variants. In this article, we introduce a hidden-Markov-model (HMM)-based method, referred to as the PICR-CNV, for copy number inference. The proposed method first estimates copy number abundance for each single SNP on a single array based on the raw fluorescence values, and then standardizes the estimated copy number abundance to achieve equal footing among multiple arrays. This method requires no between-array normalization, and thus, maintains data integrity and independence of samples among individual subjects. In addition to our efforts to apply new statistical technology to raw fluorescence values, the HMM has been applied to the standardized copy number abundance in order to reduce experimental noise. Through simulations, we show our refined method is able to infer copy number variants accurately. Application of the proposed method to a breast cancer dataset helps to identify genomic regions significantly associated with the disease.

Amplification of MPZL1/PZR promotes tumor cell migration through Src-mediated phosphorylation of cortactin in hepatocellular carcinoma

We have previously identified 1 241 regions of somatic copy number alterations (CNAs) in hepatocellular carcinoma (HCC). In the present study, we found that a novel recurrent focal amplicon, 1q24.1-24.2, targets the MPZL1 gene in HCC. Notably, there is a positive correlation between the expression levels of MPZL1 and intrahepatic metastasis of the HCC specimens. MPZL1 can significantly enhance the migratory and metastatic potential of the HCC cells. Moreover, we found that one of the mechanisms by which MPZL1 promotes HCC cell migration is by inducing the phosphorylation and activation of the pro-metastatic protein, cortactin. Additionally, we found that Src kinase mediates the phosphorylation and activation of cortactin induced by MPZL1 overexpression. Taken together, these findings suggest that MPZL1 is a novel pro-metastatic gene targeted by a recurrent region of copy number amplification at 1q24.1-24.2 in HCC.

Genetic and epigenetic analysis of putative breast cancer stem cell models

Background

Cancer stem cell model hypothesizes existence of a small proportion of tumor cells capable of sustaining tumor formation, self-renewal and differentiation. In breast cancer, these cells were found to be associated with CD44⁺CD24^-low and ALDH⁺ phenotype. Our study was performed to evaluate the suitability of current approaches for breast cancer stem cell analyses to evaluate heterogeneity of breast cancer cells through their extensive genetic and epigenetic characterization.

Methods

Breast cancer cell lines MCF7 and SUM159 were cultured in adherent conditions and as mammospheres. Flow cytometry sorting for CD44, CD24 and ALDH was performed. Sorted and unsorted populations, mammospheres and adherent cell cultures were subjected to DNA profiling by array CGH and methylation profiling by Epitect Methyl qPCR array. Methylation status of selected genes was further evaluated by pyrosequencing. Functional impact of methylation was evaluated by mRNA analysis for selected genes.

Results

Array CGH did not reveal any genomic differences. In contrast, putative breast cancer stem cells showed altered methylation levels of several genes compared to parental tumor cells.

Conclusions

Our results underpin the hypothesis that epigenetic mechanisms seem to play a major role in the regulation of CSCs. However, it is also clear that more efficient methods for CSC enrichment are needed. This work underscores requirement of additional approaches to reveal heterogeneity within breast cancer.

Alu-Alu fusion sequences identified at junction sites of copy number amplified regions in cancer cell lines

Alu elements are short, ∼300-bp stretches of DNA and are the most abundant repetitive elements in the human genome. A large number of chromosomal rearrangements mediated by Alu-Alu recombination have been reported in germline cells, but only a few in somatic cells. Cancer development is frequently accompanied by various chromosomal rearrangements including gene amplification. To explore an involvement of Alu-Alu fusion in gene amplification events, we determined 20 junction site sequences of 5 highly amplified regions in 4 cancer cell lines. The amplified regions exhibited a common copy number profile: a stair-like increase with multiple segments, which is implicated in the breakage-fusion-bridge (BFB) cycle-mediated amplification. All of the sequences determined were characterized as head-to-head or tail-to-tail fusion of sequences separated by 1-5 kb in the genome sequence. Of these, 4 junction site sequences were identified as Alu-Alu fusions between inverted, paired Alu elements with relatively long overlapping sequences of 17, 21, 22, and 24 bp. Together with genome mapping data of Alu elements, these findings suggest that when breakages occur at or near inverted, paired Alu elements in the process of BFB cycle-mediated amplification, sequence homology of Alu elements is frequently used to repair the broken ends.

Dysfunctional telomeres promote genomic instability and metastasis in the absence of telomerase activity in oncogene induced mammary cancer

Telomerase is a ribonucleoprotein that maintains the ends of chromosomes (telomeres). In normal cells lacking telomerase activity, telomeres shorten with each cell division because of the inability to completely synthesize the lagging strand. Critically shortened telomeres elicit DNA damage responses and limit cellular division and lifespan, providing an important tumor suppressor function. Most human cancer cells express telomerase which contributes significantly to the tumor phenotype. In human breast cancer, telomerase expression is predictive of clinical outcomes such as lymph node metastasis and survival. In mouse models of mammary cancer, telomerase expression is also upregulated. Telomerase overexpression resulted in spontaneous mammary tumor development in aged female mice. Increased mammary cancer also was observed when telomerase deficient mice were crossed with p53 null mutant animals. However, the effects of telomerase and telomere length on oncogene driven mammary cancer have not been completely characterized. To address these issues we characterized neu proto-oncogene driven mammary tumor formation in G1 Terc-/- (telomerase deficient with long telomeres), G3 Terc-/- (telomerase deficient with short telomeres), and Terc+/+ mice. Telomerase deficiency reduced the number of mammary tumors and increased tumor latency regardless of telomere length. Decreased tumor formation correlated with increased apoptosis in Terc deficient tumors. Short telomeres dramatically increased lung metastasis which correlated with increased genomic instability, and specific alterations in DNA copy number and gene expression. We concluded that short telomeres promote metastasis in the absence of telomerase activity in neu oncogene driven mammary tumors.

Integrative model of genomic factors for determining binding site selection by estrogen receptor-α

Besides the estrogen receptor response elements (EREs), estrogen receptor-α binding is augmented by FOXA1 co-binding, the presence of the histone mark, histone 3 monomethylated at the lysine 4 position and the presence of open chromatin.

The major determinant of ER binding is the strength of the ERE.

The differences in estrogen receptor-binding profiles between breast cancer cell lines appear to be at sites with less ‘attractive' EREs but modulated by the non-sequence factors.

A major question in transcription factor (TF) biology is why a TF binds to only a small fraction of motif eligible binding sites in the genome. Estrogen receptor-α (ERα) binds only to 2% of optimal sequence eligible sites. Using ERα as a model ligand-inducible TF, we sought to explicitly define parameters that determine TF-binding site selection on a genomic scale in an inducible system that minimizes confounding chromatin effects by the TF itself. By examining 12 separate genetic and epigenetic parameters, we find that an energetically favorable estrogen response element (ERE) motif sequence, evidence of occupancy of a ‘pioneering' TF FOXA1, the presence of the enhancer mark, H3K4me1 and an open chromatin configuration at the pre-ligand state provide specificity for ER binding. These factors can model estrogen-induced ER binding with high accuracy (ROC-AUC=0.95 using random genomic regions as background). Moreover, when ERα and FOXA1 binding as well as H3K4me1 and FAIRE localization (measure of open chromatin) was assessed in another estrogen-responsive cell line, T47D, this model was shown to be highly predictive for ER binding (ROC-AUC=0.86). Variance in binding site selection between MCF-7 and T47D appeared to reside in sites with suboptimal ERE-binding motifs but modulated by the chromatin configuration. Taken together, these results suggest definable interplay between sequence motifs and local chromatin characteristics in determining TF-binding site selection.

A major question in transcription factor (TF) biology is why a TF binds to only a small fraction of motif eligible binding sites in the genome. Using the estrogen receptor-α as a model system, we sought to explicitly define parameters that determine TF-binding site selection. By examining 12 genetic and epigenetic parameters, we find that an energetically favorable estrogen response element (ERE) motif sequence, co-occupancy by the TF FOXA1, the presence of the H3K4me1 mark and an open chromatin configuration in the pre-ligand state provide specificity for ER binding. These factors can model estrogen-induced ER binding with high accuracy (ROC-AUC=0.95 and 0.88 using different genomic backgrounds). Moreover, when assessed in another estrogen-responsive cell line, this model was highly predictive for ERα binding (ROC-AUC=0.86). Variance in binding site selection between MCF-7 and T47D resides in sites with suboptimal ERE motifs, but modulated by the chromatin configuration. These results suggest a definable interplay between sequence motifs and local chromatin in selecting TF binding.

Quantification and clinical relevance of gene amplification at chromosome 17q12-q21 in human epidermal growth factor receptor 2-amplified breast cancers

Introduction

Human epidermal growth factor receptor 2 (HER2)-amplified breast cancers represent a tumor subtype with chromosome 17q rearrangements that lead to frequent gene amplifications. The aim of this study was to quantify the amplification of genes located on chromosome 17q and to analyze the relations between the pattern of gene amplifications and the patients' characteristics and survival.

Methods

Patients with HER2-positive breast tumors (HER2 score of 3+ by immunohistochemistry or positive for HER2 amplification by fluorescence in situ hybridization (FISH)) (n = 86) and with HER2-negative breast tumors (n = 40) (negative controls) were included in this study. Using a quantitative polymerase chain reaction method and DNA extracted from frozen tumor specimens, 11 genes (MED1, STARD3, HER2, GRB7, THRA, RARA, TOP2A, IGFBP4, CCR7, KRT20, KRT19 and GAS), which are localized within Chr17q12-q21 and have a putative role in breast cancer development, were quantified. Relapse-free and overall survival rates were estimated from the date of surgery to the date of the event of interest (recurrence or death) using the Kaplan-Meier method.

Results

Gene amplification was observed only in HER2-positive tumors, and the frequency of amplification decreased with the distance of the gene from HER2. HER2 presented the highest level of amplification. TOP2A was not included in the smallest region of amplification involving HER2. Amplification of RARA, KRT20 and KRT19 was significantly associated with node-positive breast cancer (P = 0.030, P = 0.002 and P = 0.033, respectively). During a median follow-up period of 55 months (range, 6 to 81 months), the subgroup of patients with hormone receptor-negative cancer and without TOP2A amplification showed the worst survival (relapse-free survival: hazard ratio (HR) = 0.29, 95% confidence interval (95% CI), 0.13 to 0.65, P = 0.001; and overall survival: HR = 0.28, 95% CI, 0.10 to 0.76, P = 0.008).

Conclusions

HER2 amplification seems to drive genomic instability along chromosome 17q, leading to different patterns of gene amplification. This study confirms the clinical importance of identifying, among patients with HER2-positive breast tumors, the subgroup of patients with hormone receptor-negative and nonamplified TOP2A cancers as they have the worst prognosis.

ZNF217 confers resistance to the pro-apoptotic signals of paclitaxel and aberrant expression of Aurora-A in breast cancer cells

BACKGROUND

ZNF217 is a candidate oncogene located at 20q13, a chromosomal region frequently amplified in breast cancers. The precise mechanisms involved in ZNF217 pro-survival function are currently unknown, and utmost importance is given to deciphering the role of ZNF217 in cancer therapy response.

RESULTS

We provide evidence that stable overexpression of ZNF217 in MDA-MB-231 breast cancer cells conferred resistance to paclitaxel, stimulated cell proliferation in vitro associated with aberrant expression of several cyclins, and increased tumor growth in mouse xenograft models. Conversely, siRNA-mediated silencing of ZNF217 expression in MCF7 breast cancer cells, which possess high endogenous levels of ZNF217, led to decreased cell proliferation and increased sensitivity to paclitaxel. The paclitaxel resistance developed by ZNF217-overexpressing MDA-MB-231 cells was not mediated by the ABCB1/PgP transporter. However, ZNF217 was able to counteract the apoptotic signals mediated by paclitaxel as a consequence of alterations in the intrinsic apoptotic pathway through constitutive deregulation of the balance of Bcl-2 family proteins. Interestingly, ZNF217 expression levels were correlated with the oncogenic kinase Aurora-A expression levels, as ZNF217 overexpression led to increased expression of the Aurora-A protein, whereas ZNF217 silencing was associated with low Aurora-A expression levels. We showed that a potent Aurora-A kinase inhibitor was able to reverse paclitaxel resistance in the ZNF217-overexpressing cells.

CONCLUSION

Altogether, these data suggest that ZNF217 might play an important role in breast neoplastic progression and chemoresistance, and that Aurora-A might be involved in ZNF217-mediated effects.

Reinventing diagnostics for personalized therapy in oncology

Human cancers are still diagnosed and classified using the light microscope. The criteria are based upon morphologic observations by pathologists and tend to be subject to interobserver variation. In preoperative biopsies of non-small cell lung cancers, the diagnostic concordance, even amongst experienced pulmonary pathologists, is no better than a coin-toss. Only 25% of cancer patients, on average, benefit from therapy as most therapies do not account for individual factors that influence response or outcome. Unsuccessful first line therapy costs Canada CAN$1.2 billion for the top 14 cancer types, and this extrapolates to $90 billion globally. The availability of accurate drug selection for personalized therapy could better allocate these precious resources to the right therapies. This wasteful situation is beginning to change with the completion of the human genome sequencing project and with the increasing availability of targeted therapies. Both factors are giving rise to attempts to correlate tumor characteristics and response to specific adjuvant and neoadjuvant therapies. Static cancer classification and grading systems need to be replaced by functional classification systems that not only account for intra- and inter- tumor heterogeneity, but which also allow for the selection of the correct chemotherapeutic compounds for the individual patient. In this review, the examples of lung and breast cancer are used to illustrate the issues to be addressed in the coming years, as well as the emerging technologies that have great promise in enabling personalized therapy.

An integrative multi-dimensional genetic and epigenetic strategy to identify aberrant genes and pathways in cancer

Background

Genomics has substantially changed our approach to cancer research. Gene expression profiling, for example, has been utilized to delineate subtypes of cancer, and facilitated derivation of predictive and prognostic signatures. The emergence of technologies for the high resolution and genome-wide description of genetic and epigenetic features has enabled the identification of a multitude of causal DNA events in tumors. This has afforded the potential for large scale integration of genome and transcriptome data generated from a variety of technology platforms to acquire a better understanding of cancer.

Results

Here we show how multi-dimensional genomics data analysis would enable the deciphering of mechanisms that disrupt regulatory/signaling cascades and downstream effects. Since not all gene expression changes observed in a tumor are causal to cancer development, we demonstrate an approach based on multiple concerted disruption (MCD) analysis of genes that facilitates the rational deduction of aberrant genes and pathways, which otherwise would be overlooked in single genomic dimension investigations.

Conclusions

Notably, this is the first comprehensive study of breast cancer cells by parallel integrative genome wide analyses of DNA copy number, LOH, and DNA methylation status to interpret changes in gene expression pattern. Our findings demonstrate the power of a multi-dimensional approach to elucidate events which would escape conventional single dimensional analysis and as such, reduce the cohort sample size for cancer gene discovery.

Activation of estrogen-responsive genes does not require their nuclear co-localization

The spatial organization of the genome in the nucleus plays a role in the regulation of gene expression. Whether co-regulated genes are subject to coordinated repositioning to a shared nuclear space is a matter of considerable interest and debate. We investigated the nuclear organization of estrogen receptor alpha (ERalpha) target genes in human breast epithelial and cancer cell lines, before and after transcriptional activation induced with estradiol. We find that, contrary to another report, the ERalpha target genes TFF1 and GREB1 are distributed in the nucleoplasm with no particular relationship to each other. The nuclear separation between these genes, as well as between the ERalpha target genes PGR and CTSD, was unchanged by hormone addition and transcriptional activation with no evidence for co-localization between alleles. Similarly, while the volume occupied by the chromosomes increased, the relative nuclear position of the respective chromosome territories was unaffected by hormone addition. Our results demonstrate that estradiol-induced ERalpha target genes are not required to co-localize in the nucleus.

Breast cancer cell lines carry cell line-specific genomic alterations that are distinct from aberrations in breast cancer tissues: comparison of the CGH profiles between cancer cell lines and primary cancer tissues

Background

Cell lines are commonly used in various kinds of biomedical research in the world. However, it remains uncertain whether genomic alterations existing in primary tumor tissues are represented in cell lines and whether cell lines carry cell line-specific genomic alterations. This study was performed to answer these questions.

Methods

Array-based comparative genomic hybridization (CGH) was employed with 4030 bacterial artificial chromosomes (BACs) that cover the genome at 1.0 megabase resolution to analyze DNA copy number aberrations (DCNAs) in 35 primary breast tumors and 24 breast cancer cell lines. DCNAs were compared between these two groups. A tissue microdissection technique was applied to primary tumor tissues to reduce the contamination of samples by normal tissue components.

Results

The average number of BAC clones with DCNAs was 1832 (45.3% of spotted clones) and 971 (24.9%) for cell lines and primary tumor tissues, respectively. Gains of 1q and 8q and losses of 8p, 11q, 16q and 17p were detected in >50% of primary cancer tissues. These aberrations were also frequently detected in cell lines. In addition to these alterations, the cell lines showed recurrent genomic alterations including gains of 5p14-15, 20q11 and 20q13 and losses of 4p13-p16, 18q12, 18q21, Xq21.1 and Xq26-q28 that were barely detected in tumor tissue specimens. These are considered to be cell line-specific DCNAs. The frequency of the HER2 amplification was high in both cell lines and tumor tissues, but it was statistically different between cell lines and primary tumors (P = 0.012); 41.3 ± 29.9% for the cell lines and 15.9 ± 18.6% for the tissue specimens.

Conclusions

Established cell lines carry cell lines-specific DCNAs together with recurrent aberrations detected in primary tumor tissues. It must therefore be emphasized that cell lines do not always represent the genotypes of parental tumor tissues.

High frequency of common DNA copy number abnormalities detected by bacterial artificial chromosome array comparative genomic hybridization in 24 breast cancer cell lines

Breast cancer is a widespread disease in Japan and across the world. Breast cancer cells, as well as most other types of cancer cells, have diverse chromosomal aberrations. Clarifying the character of these chromosomal aberrations should contribute to the development of more suitable therapies, along with the predictions of metastasis and prognosis. Twenty-four breast cancer cell lines were analyzed by bacterial artificial chromosome (BAC) array comparative genomic hybridization (CGH). The array slide contained duplicate spots of 4030 BAC clone DNAs covering the entire human genome with 1 Mbp resolution. In all 24 breast cancer cell lines, frequent and significant amplifications as well as deletions were detected by BAC array CGH. Common DNA copy number gains, detected in 60% (above 15 cell lines) of the 24 breast cancer cell lines were found in 76 BAC clones, located at 1q, 5p, 8q, 9p, 16p, 17q, and 20q. Moreover, common DNA copy number loss was detected in 136 BAC clones, located at 1q, 2q, 3p, 4p, 6q, 8p, 9p, 11p, 13q, 17p, 18q, 19p, Xp, and Xq. The DNA copy number abnormalities found included abnormality of the well-known oncogene cMYC (8q24.21); however, most of them were not reported to relate to breast cancer. BAC array CGH has great potential to detect DNA copy number abnormalities, and has revealed that breast cancer cell lines have substantial heterogeneity.

Chimeric transcript discovery by paired-end transcriptome sequencing

Recurrent gene fusions are a prevalent class of mutations arising from the juxtaposition of 2 distinct regions, which can generate novel functional transcripts that could serve as valuable therapeutic targets in cancer. Therefore, we aim to establish a sensitive, high-throughput methodology to comprehensively catalog functional gene fusions in cancer by evaluating a paired-end transcriptome sequencing strategy. Not only did a paired-end approach provide a greater dynamic range in comparison with single read based approaches, but it clearly distinguished the high-level "driving" gene fusions, such as BCR-ABL1 and TMPRSS2-ERG, from potential lower level "passenger" gene fusions. Also, the comprehensiveness of a paired-end approach enabled the discovery of 12 previously undescribed gene fusions in 4 commonly used cell lines that eluded previous approaches. Using the paired-end transcriptome sequencing approach, we observed read-through mRNA chimeras, tissue-type restricted chimeras, converging transcripts, diverging transcripts, and overlapping mRNA transcripts. Last, we successfully used paired-end transcriptome sequencing to detect previously undescribed ETS gene fusions in prostate tumors. Together, this study establishes a highly specific and sensitive approach for accurately and comprehensively cataloguing chimeras within a sample using paired-end transcriptome sequencing.

The transcriptional induction of PIK3CA in tumor cells is dependent on the oncoprotein Y-box binding protein-1

PIK3CA, which codes for the p110alpha catalytic subunit of phosphatidylinositol-3-kinase (PI3K), is implicated as an oncogene. Despite importance of PIK3CA in cancer, little is known about what drives up its expression in tumor cells. We recently characterized the PIK3CA promoter and reported that it is transcriptionally silenced by the tumor suppressor protein p53. In the present study, we demonstrate that PIK3CA can be induced by the oncogenic transcription factor Y-box binding protein-1 (YB-1). Three YB-1-responsive elements were identified on the PIK3CA promoter using chromatin immunoprecipitation and electrophoretic mobility shift assays. Interestingly, silencing YB-1 with siRNA in models of basal-like breast cancer decreased p110alpha protein levels regardless of whether PIK3CA was wild type, amplified or mutated. This decrease in p110alpha led to a reduction in PI3K activity and the downstream signaling primarily through p90 ribosomal S6 kinase and S6 ribosomal protein. Disruption in PIK3CA-dependent signaling suppressed cellular invasion correlative with loss of urokinase plasminogen activator (uPA). Similarly, silencing YB-1 suppressed invasion and uPA production however this was reversible through the introduction of constitutively active PIK3CA. In conclusion, YB-1 is the first reported oncogene to induce the expression of PIK3CA through transcriptional control of its promoter.

Familial occurrence of dysembryoplastic neuroepithelial tumor-like neoplasm of the septum pellucidum: case report

OBJECTIVE

Dysembryoplastic neuroepithelial tumor (DNT)-like neoplasms of the septum pellucidum are extremely rare. In this article, we report the familial occurrence of these neoplasms.

CLINICAL PRESENTATION

We report two cases of such neoplasms: Patient 1, a 42-year-old woman, and Patient 2, the 20-year-old nephew of Patient 1. Patient 1 experienced headache and worsening dizziness; Patient 2 experienced headache and worsening dizziness and also had partial seizures. In both cases, magnetic resonance imaging (MRI) revealed an intraventricular tumor adjacent to the septum pellucidum. Both tumors appeared as a hypointense region on T1-weighted MRI, and both appeared as a hyperintense region on T2-weighted MRI without gadolinium enhancement. Interestingly, both tumors had a high apparent diffusion coefficient.

INTERVENTION

Both tumors were subtotally removed and had common histological findings, such as alveolar structures with oligodendroglia-like cells and "specific glioneuronal element." These findings are consistent with a dysembryoplastic neuroepithelial tumor-like neoplasm. After tumor removal, the symptoms disappeared. The postoperative course was uneventful, and the patients did not require adjuvant therapy. MRI showed no regrowth of residual tumors at 4 years (Patient 1) and 2 years (Patient 2) postoperatively.

CONCLUSION

The familial occurrence of this rare tumor suggests that both of these cases arose from a common germline mutation. Identification of this rare tumor in this rare location is important to avoid unnecessary adjuvant therapy. A markedly high apparent diffusion coefficient and histological findings of specific glioneuronal element can facilitate the differential diagnosis of dysembryoplastic neuroepithelial tumor-like neoplasms. Genetic study of affected patients in this family may provide clues to its molecular pathogenesis.

ChIP-Seq of ERalpha and RNA polymerase II defines genes differentially responding to ligands

We used ChIP-Seq to map ERalpha-binding sites and to profile changes in RNA polymerase II (RNAPII) occupancy in MCF-7 cells in response to estradiol (E2), tamoxifen or fulvestrant. We identify 10 205 high confidence ERalpha-binding sites in response to E2 of which 68% contain an estrogen response element (ERE) and only 7% contain a FOXA1 motif. Remarkably, 596 genes change significantly in RNAPII occupancy (59% up and 41% down) already after 1 h of E2 exposure. Although promoter proximal enrichment of RNAPII (PPEP) occurs frequently in MCF-7 cells (17%), it is only observed on a minority of E2-regulated genes (4%). Tamoxifen and fulvestrant partially reduce ERalpha DNA binding and prevent RNAPII loading on the promoter and coding body on E2-upregulated genes. Both ligands act differently on E2-downregulated genes: tamoxifen acts as an agonist thus downregulating these genes, whereas fulvestrant antagonizes E2-induced repression and often increases RNAPII occupancy. Furthermore, our data identify genes preferentially regulated by tamoxifen but not by E2 or fulvestrant. Thus (partial) antagonist loaded ERalpha acts mechanistically different on E2-activated and E2-repressed genes.

Integrative genomic and gene expression analysis of chromosome 7 identified novel oncogene loci in non-small cell lung cancer

Lung cancer accounts for over a quarter of cancer deaths, with non-small cell lung cancer (NSCLC) accounting for approximately 80% of cases. Several genome studies have been undertaken in both cell models of NSCLC and clinical samples to identify alterations underlying disease behaviour, and many have identified recurring aberrations of chromosome 7. The presence of recurring chromosome 7 alterations that do not span the well-studied oncogenes EGFR (at 7p11.2) and MET (at 7q31.2) has raised the hypothesis of additional genes on this chromosome that contribute to tumourigenesis. In this study, we demonstrated that multiple loci on chromosome 7 are indeed amplified in NSCLC, and through integrative analysis of gene dosage alterations and parallel gene expression changes, we identified new lung cancer oncogene candidates, including FTSJ2, NUDT1, TAF6, and POLR2J. Activation of these key genes was confirmed in panels of clinical lung tumour tissue as compared with matched normal lung tissue. The detection of gene activation in multiple cohorts of samples strongly supports the presence of key genes involved in lung cancer that are distinct from the EGFR and MET loci on chromosome 7.