CAMK1D amplification implicated in epithelial-mesenchymal transition in basal-like breast cancer

Transformation of an olfactory placode-derived cell into one with stem cell characteristics by disrupting epigenetic barriers

The mammalian olfactory neuronal lineage is regenerative, and accordingly, maintains a population of pluripotent cells that replenish olfactory sensory neurons and other olfactory cell types during the life of the animal. Moreover, in response to acute injury, the early transit amplifying cells along the olfactory sensory neuronal lineage are able to de-differentiate to shift resources in support of tissue restoration. In order to further explore plasticity of various cellular stages along the olfactory sensory neuronal lineage, we challenged the epigenetic stability of two olfactory placode-derived cell lines that model immature olfactory sensory neuronal stages. We found that perturbation of the Ehmt2 chromatin modifier transformed the growth properties, morphology, and gene expression profiles towards states with several stem cell characteristics. This transformation was dependent on continued expression of the large T-antigen, and was enhanced by Sox2 over-expression. These findings may provide momentum for exploring inherent cellular plasticity within early cell types of the olfactory lineage, as well as potentially add to our knowledge of cellular reprogramming.

SUMMARY STATEMENT

Discovering how epigenetic modifications influence olfactory neuronal lineage plasticity offers insights into regenerative potential and cellular reprogramming.

Repositioning VU-0365114 as a novel microtubule-destabilizing agent for treating cancer and overcoming drug resistance

Microtubule-targeting agents represent one of the most successful classes of anticancer agents. However, the development of drug resistance and the appearance of adverse effects hamper their clinical implementation. Novel microtubule-targeting agents without such limitations are urgently needed. By employing a gene expression-based drug repositioning strategy, this study identifies VU-0365114, originally synthesized as a positive allosteric modulator of human muscarinic acetylcholine receptor M5 (M5 mAChR), as a novel type of tubulin inhibitor by destabilizing microtubules. VU-0365114 exhibits a broad-spectrum in vitro anticancer activity, especially in colorectal cancer cells. A tumor xenograft study in nude mice shows that VU-0365114 slowed the in vivo colorectal tumor growth. The anticancer activity of VU-0365114 is not related to its original target, M5 mAChR. In addition, VU-0365114 does not serve as a substrate of multidrug resistance (MDR) proteins, and thus, it can overcome MDR. Furthermore, a kinome analysis shows that VU-0365114 did not exhibit other significant off-target effects. Taken together, our study suggests that VU-0365114 primarily targets microtubules, offering potential for repurposing in cancer treatment, although more studies are needed before further drug development.

Bone-matrix mineralization dampens integrin-mediated mechanosignalling and metastatic progression in breast cancer

In patients with breast cancer, lower bone mineral density increases the risk of bone metastasis. Although the relationship between bone-matrix mineralization and tumour-cell phenotype in breast cancer is not well understood, mineralization-induced rigidity is thought to drive metastatic progression via increased cell-adhesion forces. Here, by using collagen-based matrices with adjustable intrafibrillar mineralization, we show that, unexpectedly, matrix mineralization dampens integrin-mediated mechanosignalling and induces a less proliferative stem-cell-like phenotype in breast cancer cells. In mice with xenografted decellularized physiological bone matrices seeded with human breast tumour cells, the presence of bone mineral reduced tumour growth and upregulated a gene-expression signature that is associated with longer metastasis-free survival in patients with breast cancer. Our findings suggest that bone-matrix changes in osteogenic niches regulate metastatic progression in breast cancer and that in vitro models of bone metastasis should integrate organic and inorganic matrix components to mimic physiological and pathologic mineralization.

Adipocyte-derived CCDC3 promotes tumorigenesis in epithelial ovarian cancer through the Wnt/ß-catenin signalling pathway

Introduction

Epithelial ovarian cancer (EOC) is a highly aggressive disease whose unique metastatic site is the omentum. Coiled-coil domain containing 3 (CCDC3) is an adipocyte-derived secreted protein that is specifically elevated in omental adipose tissue. However, its function is still unknown.

Material and methods

Initially, a Kaplan-Meier plot was applied to evaluate the prognostic value of CCDC3 expression in patients with EOC. A bioinformatics analysis was next used to explore the biological function of CCDC3 in EOC. Western blot, quantitative real-time polymerase chain reaction, and in vitro invasion and migration assays were performed using SKOV3 cells and CCDC3 secreted by rat adipocytes to analyzes the impact of CCDC3 on EOC and the underlying mechanism.

Results

Overexpression of CCDC3 was associated with poor prognosis of EOC. CCDC3 interacted with multiple key signalling pathways, including the Wnt/β-catenin pathway. EOC cellular proliferation, migration, and invasion were promoted in vitro when co-cultured with CCDC3 enriched conditioned medium, and this tumour-promoting effect was induced by activating the Wnt/β-catenin pathway. Furthermore, the epithelial-mesenchymal transition of EOC cells was reversed after CCDC3 silencing.

Conclusions

Our results support that CCDC3 promotes EOC tumorigenesis through the Wnt/β-catenin pathway and that CCDC3 may serve as a novel prognostic biomarker and therapeutic target for metastatic EOC.

Early detection of colorectal cancer based on circular DNA and common clinical detection indicators

BACKGROUND

Colorectal cancer (CRC) is the third most common cancer worldwide, and it is the second leading cause of death from cancer in the world, accounting for approximately 9% of all cancer deaths. Early detection of CRC is urgently needed in clinical practice.

AIM

To build a multi-parameter diagnostic model for early detection of CRC.

METHODS

Total 59 colorectal polyps (CRP) groups, and 101 CRC patients (38 early-stage CRC and 63 advanced CRC) for model establishment. In addition, 30 CRP groups, and 62 CRC patients (30 early-stage CRC and 32 advanced CRC) were separately included to validate the model. 51 commonly used clinical detection indicators and the 4 extrachromosomal circular DNA markers NDUFB7, CAMK1D, PIK3CD and PSEN2 that we screened earlier. Four multi-parameter joint analysis methods: binary logistic regression analysis, discriminant analysis, classification tree and neural network to establish a multi-parameter joint diagnosis model.

RESULTS

Neural network included carcinoembryonic antigen (CEA), ischemia-modified albumin (IMA), sialic acid (SA), PIK3CD and lipoprotein a (LPa) was chosen as the optimal multi-parameter combined auxiliary diagnosis model to distinguish CRP and CRC group, when it differentiated 59 CRP and 101 CRC, its overall accuracy was 90.8%, its area under the curve (AUC) was 0.959 (0.934, 0.985), and the sensitivity and specificity were 91.5% and 82.2%, respectively. After validation, when distinguishing based on 30 CRP and 62 CRC patients, the AUC was 0.965 (0.930-1.000), and its sensitivity and specificity were 66.1% and 70.0%. When distinguishing based on 30 CRP and 32 early-stage CRC patients, the AUC was 0.960 (0.916-1.000), with a sensitivity and specificity of 87.5% and 90.0%, distinguishing based on 30 CRP and 30 advanced CRC patients, the AUC was 0.970 (0.936-1.000), with a sensitivity and specificity of 96.7% and 86.7%.

CONCLUSION

We built a multi-parameter neural network diagnostic model included CEA, IMA, SA, PIK3CD and LPa for early detection of CRC, compared to the conventional CEA, it showed significant improvement.

Identification of epithelial and mesenchymal circulating tumor cells in clonal lineage of an aggressive prostate cancer case

Little is known about the complexity and plasticity of circulating tumor cell (CTC) biology in different compartments of the fluid microenvironment during tumor metastasis. Here we integrated phenomics, genomics, and targeted proteomics to characterize CTC phenotypic and genotypic heterogeneity in paired peripheral blood (PB) and bone marrow aspirate (BMA) from a metastatic prostate cancer patient following the rapid disease progression, using the High-Definition Single Cell Assay 3.0 (HDSCA3.0). Uniquely, we identified a subgroup of genetically clonal CTCs that acquired a mesenchymal-like state and its presence was significantly associated with one subclone that emerged along the clonal lineage. Higher CTC abundance and phenotypic diversity were observed in the BMA than PB and differences in genomic alterations were also identified between the two compartments demonstrating spatial heterogeneity. Single cell copy number profiling further detected clonal heterogeneity within clusters of CTCs (also known as microemboli or aggregates) as well as phenotypic variations by targeted proteomics. Overall, these results identify epithelial and mesenchymal CTCs in the clonal lineage of an aggressive prostate cancer case and also demonstrate a single cell multi-omic approach to deconvolute the heterogeneity and association of CTC phenotype and genotype in multi-medium liquid biopsies of metastatic prostate cancer.

CAMK1D Inhibits Glioma Through the PI3K/AKT/mTOR Signaling Pathway

Calcium/calmodulin-dependent protein ID (CAMK1D) is widely expressed in many tissues and involved in tumor cell growth. However, its role in gliomas has not yet been elucidated. This study aimed to investigate the roles of CAMK1D in the proliferation, migration, and invasion of glioma. Through online datasets, Western blot, and immunohistochemical analysis, glioma tissue has significantly lower CAMK1D expression levels than normal brain (NB) tissues, and CAMK1D expression was positively correlated with the WHO classification. Kaplan-Meier survival analysis shows that CAMK1D can be used as a potential prognostic indicator to predict the overall survival of glioma patients. In addition, colony formation assay, cell counting Kit-8, and xenograft experiment identified that knockdown of CAMK1D promotes the proliferation of glioma cells. Transwell and wound healing assays identified that knockdown of CAMK1D promoted the invasion and migration of glioma cells. In the above experiments, the results of overexpression of CAMK1D were all contrary to those of knockdown. In terms of mechanism, this study found that CAMK1D regulates the function of glioma cells by the PI3K/AKT/mTOR pathway. In conclusion, these findings suggest that CAMK1D serves as a prognostic predictor and a new target for developing therapeutics to treat glioma.

PDSS1-Mediated Activation of CAMK2A-STAT3 Signaling Promotes Metastasis in Triple-Negative Breast Cancer

Genomic alterations are crucial for the development and progression of human cancers. Copy-number gains found in genes encoding metabolic enzymes may induce triple-negative breast cancer (TNBC) adaptation. However, little is known about how metabolic enzymes regulate TNBC metastasis. Using our previously constructed multiomic profiling of a TNBC cohort, we identified decaprenyl diphosphate synthase subunit 1 (PDSS1) as an essential gene for TNBC metastasis. PDSS1 expression was significantly upregulated in TNBC tissues compared with adjacent normal tissues and was positively associated with poor survival among patients with TNBC. PDSS1 knockdown inhibited TNBC cell migration, invasion, and distant metastasis. Mechanistically, PDSS1, but not a catalytically inactive mutant, positively regulated the cellular level of coenzyme Q10 (CoQ10) and intracellular calcium levels, thereby inducing CAMK2A phosphorylation, which is essential for STAT3 phosphorylation in the cytoplasm. Phosphorylated STAT3 entered the nucleus, promoting oncogenic STAT3 signaling and TNBC metastasis. STAT3 phosphorylation inhibitors (e.g., Stattic) effectively blocked PDSS1-induced cell migration and invasion in vitro and tumor metastasis in vivo. Taken together, our study highlights the importance of targeting the previously uncharacterized PDSS1/CAMK2A/STAT3 oncogenic signaling axis, expanding the repertoire of precision medicine in TNBC. SIGNIFICANCE: A novel metabolic gene PDSS1 is highly expressed in triple-negative breast cancer tissues and contributes to metastasis, serving as a potential therapeutic target for combating metastatic disease.

Biochemical characterization of four splice variants of mouse Ca2+/calmodulin-dependent protein kinase Iδ

Ca2+/calmodulin (CaM)-dependent protein kinase Iδ (CaMKIδ) is a Ser/Thr kinase that plays pivotal roles in Ca2+ signalling. CaMKIδ is activated by Ca2+/CaM-binding and phosphorylation at Thr180 by CaMK kinase (CaMKK). In this study, we characterized four splice variants of mouse CaMKIδ (mCaMKIδs: a, b, c and d) found by in silico analysis. Recombinant mCaMKIδs expressed in Escherichia coli were phosphorylated by CaMKK; however, only mCaMKIδ-a and c showed protein kinase activities towards myelin basic protein in vitro, with mCaMKIδ-b and mCaMKIδ-d being inactive. Although mCaMKIδ-a and mCaMKIδ-c underwent autophosphorylation in vitro, only mCaMKIδ-c underwent autophosphorylation in 293T cells. Site-directed mutagenesis showed that the autophosphorylation site is Ser349, which is found in the C-terminal region of only variants c and b (Ser324). Furthermore, phosphorylation of these sites (Ser324 and Ser349) in mCaMKIδ-b and c was more efficiently catalyzed by cAMP-dependent protein kinase in vitro and in cellulo as compared to the autophosphorylation of mCaMKIδ-c. Thus, variants of mCaMKIδ possess distinct properties in terms of kinase activities, autophosphorylation and phosphorylation by another kinase, suggesting that they play physiologically different roles in murine cells.

LncRNA ILF3-AS1 Promotes the Progression of Colon Adenocarcinoma Cells Through the miR-619-5p/CAMK1D Axis

Introduction

Colon adenocarcinoma (COAD) is the third most common tumor of the digestive tract. Recent studies reported that lncRNA’s abnormal expression might play a vital role in the occurrence and development of COAD.

Methods

In the present study, we investigated the expression of ILF3-AS1 in COAD cell lines, human normal colon epithelial cell line, patient tumor tissues and adjacent normal tissues by real-time quantitative PCR (RT-qPCR). Small interfering RNAs (siRNAs) were transfected into COAD cells to inhibit the expression of ILF3-AS1. The effects of ILF3-AS1 on cell proliferation, migration, invasion and apoptosis were measured by CCK-8 assay, transwell migration and invasion assay, and flow cytometry apoptosis assay, respectively. The direct binding of ILF3-AS1 and miR-619-5p/CAMK1D was validated by the luciferase reporter assay. The expression of CAMK1D and epithelial-mesenchymal transformation (EMT)-related proteins was detected by Western Blot analysis. Besides, in vivo experiments were conducted to demonstrate the oncogenic role of ILF3-AS1 in COAD.

Results

The results showed that the expression of ILF3-AS1 was significantly higher in COAD tissue than in normal adjacent samples, and this conclusion was confirmed in the available public datasets. After ILF3-AS1 knockdown, the proliferation of COAD cell lines SW480 and HT29 was significantly inhibited. At the same time, the apoptosis was increased, and the invasion and migration abilities were decreased. After further exploring the mechanisms, we found that ILF3-AS1 serves as a competitive endogenous RNA of mir-619-5p. It can bind to mir-619-5p and reduce its expression, thus regulating the target gene CAMK1D. In addition, we found that high expression of ILF3-AS1 was significantly associated with tumor grade, tumor size, and distant metastasis in COAD samples. In vivo experiments confirmed that ILF3-AS1 promotes tumor growth in COAD models.

Conclusion

The present study demonstrated that ILF3-AS1 plays an oncogenic role in COAD through regulating the miR-619-5p/CAMK1D axis, and inhibition of ILF3-AS1 may pave the way for COAD treatment.

Backbone resonance assignments of the catalytic and regulatory domains of Ca2+/calmodulin-dependent protein kinase 1D

The CaMK subfamily of Ser/Thr kinases are regulated by calmodulin interactions with their C-terminal regions. They are exemplified by Ca²⁺/calmodulin dependent protein kinase 1δ which is known as CaMK1D, CaMKIδ or CKLiK. CaMK1D mediates intracellular signalling downstream of Ca²⁺ influx and thereby exhibits amplifications of Ca²⁺signals and polymorphisms that have been implicated in breast cancer and diabetes. Here we report the backbone ¹H, ¹³C, ¹⁵N assignments of the 38 kDa human CaMK1D protein in its free state, including both the canonical bi-lobed kinase fold as well as the autoinhibitory and calmodulin binding domains.

CAMK1D Triggers Immune Resistance of Human Tumor Cells Refractory to Anti-PD-L1 Treatment

The success of cancer immunotherapy is limited by resistance to immune checkpoint blockade. We therefore conducted a genetic screen to identify genes that mediated resistance against CTLs in anti-PD-L1 treatment-refractory human tumors. Using PD-L1-positive multiple myeloma cells cocultured with tumor-reactive bone marrow-infiltrating CTL as a model, we identified calcium/calmodulin-dependent protein kinase 1D (CAMK1D) as a key modulator of tumor-intrinsic immune resistance. CAMK1D was coexpressed with PD-L1 in anti-PD-L1/PD-1 treatment-refractory cancer types and correlated with poor prognosis in these tumors. CAMK1D was activated by CTL through Fas-receptor stimulation, which led to CAMK1D binding to and phosphorylating caspase-3, -6, and -7, inhibiting their activation and function. Consistently, CAMK1D mediated immune resistance of murine colorectal cancer cells in vivo The pharmacologic inhibition of CAMK1D, on the other hand, restored the sensitivity toward Fas-ligand treatment in multiple myeloma and uveal melanoma cells in vitro Thus, rapid inhibition of the terminal apoptotic cascade by CAMK1D expressed in anti-PD-L1-refractory tumors via T-cell recognition may have contributed to tumor immune resistance.

Discovery of Highly Selective Inhibitors of Calmodulin-Dependent Kinases That Restore Insulin Sensitivity in the Diet-Induced Obesity in Vivo Mouse Model

Polymorphisms in the region of the calmodulin-dependent kinase isoform D (CaMK1D) gene are associated with increased incidence of diabetes, with the most common polymorphism resulting in increased recognition by transcription factors and increased protein expression. While reducing CaMK1D expression has a potentially beneficial effect on glucose processing in human hepatocytes, there are no known selective inhibitors of CaMK1 kinases that can be used to validate or translate these findings. Here we describe the development of a series of potent, selective, and drug-like CaMK1 inhibitors that are able to provide significant free target cover in mouse models and are therefore useful as in vivo tool compounds. Our results show that a lead compound from this series improves insulin sensitivity and glucose control in the diet-induced obesity mouse model after both acute and chronic administration, providing the first in vivo validation of CaMK1D as a target for diabetes therapeutics.

Comprehensive Proteomic Characterization Reveals Subclass-Specific Molecular Aberrations within Triple-negative Breast Cancer

Triple-negative breast cancer (TNBC) is the most aggressive subtype of breast cancer lacking targeted therapies. This is attributed to its high heterogeneity that complicates elucidation of its molecular aberrations. Here, we report identification of specific proteome expression profiles pertaining to two TNBC subclasses, basal A and basal B, through in-depth proteomics analysis of breast cancer cells. We observed that kinases and proteases displayed unique expression patterns within the subclasses. Systematic analyses of protein-protein interaction and co-regulation networks of these kinases and proteases unraveled dysregulated pathways and plausible targets for each TNBC subclass. Among these, we identified kinases AXL, PEAK1, and TGFBR2 and proteases FAP, UCHL1, and MMP2/14 as specific targets for basal B subclass, which represents the more aggressive TNBC cell lines. Our study highlights intricate mechanisms and distinct targets within TNBC and emphasizes that these have to be exploited in a subclass-specific manner rather than a one-for-all TNBC therapy.

•

Proteome profiling reveals functionally distinct subclasses within TNBC

•

Kinases and proteases underlie unique functional signatures among the subclasses

•

Kinase-protease-centric networks highlight subclass-specific molecular rewiring

•

Protein association dysregulations reveal TNBC subclass-specific protein targets

The Role of Calmodulin in Tumor Cell Migration, Invasiveness, and Metastasis

Calmodulin (CaM) is the principal Ca²⁺ sensor protein in all eukaryotic cells, that upon binding to target proteins transduces signals encoded by global or subcellular-specific changes of Ca²⁺ concentration within the cell. The Ca²⁺/CaM complex as well as Ca²⁺-free CaM modulate the activity of a vast number of enzymes, channels, signaling, adaptor and structural proteins, and hence the functionality of implicated signaling pathways, which control multiple cellular functions. A basic and important cellular function controlled by CaM in various ways is cell motility. Here we discuss the role of CaM-dependent systems involved in cell migration, tumor cell invasiveness, and metastasis development. Emphasis is given to phosphorylation/dephosphorylation events catalyzed by myosin light-chain kinase, CaM-dependent kinase-II, as well as other CaM-dependent kinases, and the CaM-dependent phosphatase calcineurin. In addition, the role of the CaM-regulated small GTPases Rac1 and Cdc42 (cell division cycle protein 42) as well as CaM-binding adaptor/scaffold proteins such as Grb7 (growth factor receptor bound protein 7), IQGAP (IQ motif containing GTPase activating protein) and AKAP12 (A kinase anchoring protein 12) will be reviewed. CaM-regulated mechanisms in cancer cells responsible for their greater migratory capacity compared to non-malignant cells, invasion of adjacent normal tissues and their systemic dissemination will be discussed, including closely linked processes such as the epithelial–mesenchymal transition and the activation of metalloproteases. This review covers as well the role of CaM in establishing metastatic foci in distant organs. Finally, the use of CaM antagonists and other blocking techniques to downregulate CaM-dependent systems aimed at preventing cancer cell invasiveness and metastasis development will be outlined.

The novel long noncoding RNA LOC283070 is involved in the transition of LNCaP cells into androgen-independent cells via its interaction with PHB2

We sought to investigate the underlying mechanism of action of the long noncoding RNA (lncRNA) LOC283070 in the development of androgen independence in prostate cancer. The interactions between LOC283070 and target proteins were investigated by RNA pull-down and RNA-binding protein immunoprecipitation (RIP) assays. Subcellular fractionation and quantitative reverse transcription-polymerase chain reaction (qRT-PCR) were used to detect the subcellular localization of LOC283070. Western blotting was performed to detect the expression of prohibitin 2 (PHB2). Luciferase activity assays were performed to evaluate the effects of LOC283070 and PHB2 on the androgen receptor (AR) signaling pathway. A methyl thiazolyl tetrazolium (MTT) assay and a growth curve assay were used to test cell viability. Flow cytometry was performed to analyze cell cycles. A transwell assay was employed to test cell migration. We identified PHB2 as an interaction partner of LOC283070 in the pull-down and RIP experiments. Furthermore, we confirmed that the enrichment of LOC283070 with PHB2 in androgen-independent LNCaP (LNCaP-AI) cells was much greater than that in LNCaP cells. Moreover, the expression of PHB2 was not significantly different between the two cell lines, and the expression of LOC283070 in the nuclei of the LNCaP-AI cells was significantly greater than that in the LNCaP cells. In vitro data revealed that PHB2 overexpression significantly inhibited AR activity and cell proliferation and migration and induced accumulation of prostate cancer cells in G0/G1 phase. Moreover, the overexpression of LOC283070 fully abrogated the effects of PHB2 overexpression. In conclusion, we found that LOC283070 can bind to PHB2 located in the nucleus and inhibit its effect, and this is one of the mechanisms by which LOC283070 is involved in the transition of LNCaP cells into androgen-independent cells.

Role of the calcium toolkit in cancer stem cells

Cancer stem cells are a subpopulation of tumor cells that proliferate, self-renew and produce more differentiated tumoral cells building-up the tumor. Responsible for the sustained growth of malignant tumors, cancer stem cells are proposed to play significant roles in cancer resistance to standard treatment and in tumor recurrence. Among the mechanisms dysregulated in neoplasms, those related to Ca²⁺ play significant roles in various aspects of cancers. Ca²⁺ is a ubiquitous second messenger whose fluctuations of its intracellular concentrations are tightly controlled by channels, pumps, exchangers and Ca²⁺ binding proteins. These components support the genesis of Ca²⁺ signals with specific spatio-temporal characteristics that define the cell response. Being involved in the coupling of extracellular events with intracellular responses, the Ca²⁺ toolkit is often hijacked by cancer cells to promote notably their proliferation and invasion. Growing evidence obtained during the last decade pointed to a role of Ca²⁺ handling and mishandling in cancer stem cells. In this review, after a general overview of the concept of cancer stem cells we analyse and discuss the studies and current knowledge regarding the complex roles of Ca²⁺ toolkit and signaling in these cells. We highlight that numbers of Ca²⁺ signaling actors promote cancer stem cell state and are associated with cell resistance to current cancer treatments and thus may represent promising targets for potential clinical applications.

Analysis of microRNA expression profiles reveals a 5‑microRNA prognostic signature for predicting overall survival time in patients with gastric adenocarcinoma

There is growing evidence supporting dysregulated microRNAs (miRNAs) as potential prognostic biomarkers in cancer. The present study aimed to identify an miRNA model set with prognostic power for patients with gastric adenocarcinoma. miRNA‑seq data from 155 patients and 37 controls were downloaded from The Cancer Genome Atlas (TCGA) database for a comprehensive analysis of miRNA expression profiles and were used as training data. A total of 5 prognostic miRNAs, which have not been previously reported, were identified using univariate and multivariate Cox regression analyses. A separate 155‑patient TCGA cohort was used as a validation set for evaluation of the risk model. Patients in the training set were assigned into high‑ and low‑risk groups according to the 5‑miRNA signature risk scores. Kaplan‑Meier survival analyses demonstrated that patients with high risk scores had significantly shorter survival times than those with low risk scores. The risk model validation confirmed the prognostic ability of this 5‑miRNA signature in predicting the risk status of patients. Stratification analysis for clinical prognostic variables demonstrated recurrence and age were significant prognostic factors in the low‑ and high‑risk groups, respectively. In conclusion, the present 5‑miRNA signature is a potential independent risk factor for patient outcomes. The risk model based on the 5‑miRNA signature performed well in predicting overall survival time in patients with gastric adenocarcinoma.

Gene-Specific Intron Retention Serves as Molecular Signature that Distinguishes Melanoma from Non-Melanoma Cancer Cells in Greek Patients

Background: Skin cancer represents the most common human malignancy, and it includes BCC, SCC, and melanoma. Since melanoma is one of the most aggressive types of cancer, we have herein attempted to develop a gene-specific intron retention signature that can distinguish BCC and SCC from melanoma biopsy tumors. Methods: Intron retention events were examined through RT-sqPCR protocols, using total RNA preparations derived from BCC, SCC, and melanoma Greek biopsy specimens. Intron-hosted miRNA species and their target transcripts were predicted via the miRbase and miRDB bioinformatics platforms, respectively. Ιntronic ORFs were recognized through the ORF Finder application. Generation and visualization of protein interactomes were achieved by the IntAct and Cytoscape softwares, while tertiary protein structures were produced by using the I-TASSER online server. Results: c-MYC and Sestrin-1 genes proved to undergo intron retention specifically in melanoma. Interaction maps of proteins encoded by genes being potentially targeted by retained intron-accommodated miRNAs were generated and SRPX2 was additionally delivered to our melanoma-specific signature. Novel ORFs were identified in MCT4 and Sestrin-1 introns, with potentially critical roles in melanoma development. Conclusions: The property of c-MYC, Sestrin-1, and SRPX2 genes to retain specific introns could be clinically used to molecularly differentiate non-melanoma from melanoma tumors.

An integrated investigation of oocyte developmental competence: expression of key genes in human cumulus cells, morphokinetics of early divisions, blastulation, and euploidy

PURPOSE

To investigate the association of cumulus cell (CC)-related expression of a selected cluster of key genes (PTGS2, CAMK1D, HAS2, STC1, and EFNB2) with embryo development to blastocyst.

METHODS

Exploratory study at a private clinic. Eighteen advanced maternal age patients were enrolled (37.3 ± 4.0 years). Seventy-five cumuli were collected, whose oocytes resulted in either developmental arrest (N = 33) or blastocyst formation (N = 42). The noninvasive CC gene expression was combined with time-lapse morphokinetic parameters and, for blastocysts, with qPCR-based aneuploidy testing on trophectoderm biopsies.

RESULTS

The detection rate was 100% for all transcripts, but STC1 (96%) and CAMK1D (89%). Among amplified assays, CC mean expression levels of CAMK1D, PTGS2, and HAS2 were higher from oocytes that developed to blastocyst. No difference in CC key gene expression was reported between euploid (N = 21) and aneuploid (N = 21) blastocysts. Some timings of early embryo development were faster in embryos developing to blastocyst (time of pronuclei appearance and fading, division to two- and four-cells, first and second cell cycles). However, the generalized linear models outlined increasing CAMK1D expression levels as the strongest parameter associated with oocytes' developmental potential from both a general (AUC = 0.78 among amplified samples) and an intrapatient perspectives (AUC = 0.9 among patients obtaining ≥ 2 zygotes from the cohort with different developmental outcomes).

CONCLUSIONS

CAMK1D level of expression in CCs associated with blastocyst development. If confirmed from larger studies in wider populations of patients, the investigation of CC key gene expression might suit IVF clinics not adopting blastocyst culture. Future investigations should clarify the role of CAMK1D in ovarian physiology and could provide novel insights on how oocytes gain competence during folliculogenesis.

Deduction of novel genes potentially involved in hypoxic AC16 human cardiomyocytes using next-generation sequencing and bioinformatics approaches

Atherosclerotic cardiovascular disease and acute myocardial infarction are the leading causes of mortality worldwide, and apoptosis is the major pathway of cardiomyocyte death under hypoxic conditions. Although studies have reported changes in the expression of certain pro‑apoptotic and anti‑apoptotic genes in hypoxic cardiomyocytes, genetic regulations are complex in human cardiomyocytes and there is much that remains to be fully elucidated. The present study aimed to identify differentially expressed genes in hypoxic human AC16 cardiomyocytes using next‑generation sequencing and bioinformatics. A total of 24 genes (15 upregulated and 9 downregulated) with potential micro (mi)RNA‑mRNA interactions were identified in the miRmap database. Utilising the Gene Expression Omnibus database of cardiac microvascular endothelial cells, tensin 1, B‑cell lymphoma 2‑interacting protein 3 like, and stanniocalcin 1 were found to be upregulated, and transferrin receptor and methyltransferase like 7A were found to be downregulated in response to hypoxia. Considering the results from miRmap, TargetScan and miRDB together, two potential miRNA‑mRNA interactions were identified: hsa‑miRNA (miR)‑129‑5p/CDC42EP3 and hsa‑miR‑330‑3p/HELZ. These findings contribute important insights into possible novel diagnostic or therapeutic strategies for targeting cardiomyocytes under acute hypoxic stress in conditions, including acute myocardial infarction. The results of the present study also introduce an important novel approach in investigating acute hypoxic pathophysiology.

Imaging Cobalamin Uptake within Malignant Breast Tumors In Vivo

PURPOSE

To image the uptake of cobalamin (Cbl) within malignant breast tumors in vivo.

PROCEDURES

Prior to surgery 20 female patients with clinically suspected breast tumors were intravenously administered 0.25 μg of an In-111 labeled 5-deoxyadenosylcobalamin (AC) analog ([¹¹¹In]AC) and sequentially imaged with whole-body planar (WBP) and single-photon emission computed tomography (SPECT) between 2-5 h and 20-24 h post-injection (P.I.). The tumor to background (T/B) ratio for [¹¹¹In]AC in breast tumors at 2-5 h was correlated to its expression of estrogen (ER), progesterone (PR), and human epidermal growth factor 2 (HER2) receptors. Subsequent pulse chase (PC) experiments in nude mice burdened with the MDA-MB-231 triple-negative (TN) breast tumor xenograft measured the effect that pulses of AC or dexamethasone (DEX) had on [¹¹¹In]AC uptake in both normal murine tissue and the TN breast tumor.

RESULTS

The mean [¹¹¹In]AC T/B ratio of the patients' 18 resected tumors was 5.8. Comparing ER- and PR-positive tumors (n = 11) to TN and HER2-positive tumors (n = 7), the mean [¹¹¹In]AC T/B ratios at 2-5 h P.I. were 3.2 (range 1.8-5.6) and 10.4 (range 3.3-22.5), respectively. Pulses of 2.0 μg of AC at 2, 8, or 24 h; or 40.0 μg of DEX at 24 h prior to injecting 0.5 μg of [¹¹¹In]AC, increased mean tracer uptake in the MDA-MB-231 tumors by 26.4, 71.5, 92.6, and 49.1 %, respectively. Only the 2- and 24-h PC intervals concomitantly suppressed [¹¹¹In]AC uptake in normal murine tissue while enhancing [¹¹¹In]AC uptake in MDA-MB-231 tumors.

CONCLUSION

The uptake of Cbl within malignant breast tumors can be imaged clinically. Cbl uptake is greatest in TN and HER2-positive breast tumors. A solitary bolus of AC or DEX increases the [¹¹¹In]AC uptake within a breast tumor in vivo. Investigating the cytogenetic mechanisms controlling the endocytosis of Cbl in malignant breast tumors is warranted.

FGF4 induces epithelial-mesenchymal transition by inducing store-operated calcium entry in lung adenocarcinoma

Several fibroblast growth factor (FGF) isoforms act to stimulate epithelial-mesenchymal transition (EMT) during cancer progression. FGF4 and FGF7 are two ligands of FGF receptor 2 (FGFR2). Using two lung adenocarcinoma (ADC) cell lines, A549 and H1299, we showed that FGF4, but not FGF7, altered cell morphology, promoted EMT-associated protein expression, and enhanced cell proliferation, migration/invasion and colony initiation. In addition, FGF4 increased store-operated calcium entry (SOCE) and expression of the calcium signal-associated protein Orai1. The SOCE inhibitor 2,5-di-tert-butylhydroquinone (BHQ) or Orai1 knockdown reversed all of the EMT-promoting effects of FGF4. BHQ also inhibited FGF4-induced EMT in a mouse xenograft model. Finally, 60 human lung ADC samples and 21 sets of matched specimens (primary and metastatic foci in lymph nodes from one patient) were used to confirm the clinicopathologic significance of FGF4 and its correlation with E-cadherin, Vimentin and Orai1 expression. Our study thus shows that FGF4 induces EMT by elevating SOCE in lung ADC.

Differential expression and clinical significance of epithelial-mesenchymal transition markers among different histological types of triple-negative breast cancer

Background: Triple-negative breast cancer (TNBC) is a heterogeneous disease closely associated with epithelial-to-mesenchymal transition (EMT). This study aimed to investigate the role of EMT in metaplastic carcinoma.

Methods: E-cadherin, Slug, Twist and Vimentin levels were detected by immunohistochemistry in 167 TNBC tumors, including 145 invasive carcinomas of no special type (ICONSTs), 14 spindle cell carcinomas (SpCCs) and 8 matrix-producing carcinomas (MPCs).

Results: Nuclear Slug and Twist were more frequently detected in SpCC and MPC tumors than that in ICONST tumors (p<0.001). The rate of E-cadherin loss was much lower in the ICONST tumors than that in the SpCC and MPC tumors (p<0.001). Vimentin was expressed in all SpCC and MPC tumors. Furthermore, nuclear expression of Slug and Twist was positively associated with the cytoplasmic localization of Vimentin (p<0.001) and was inversely associated with membranous staining of E-cadherin (p<0.001). These trends were more apparent in the SpCC and MPC tumors than in the ICONST tumors. Follow-up data were available for 151 patients. The follow-up times ranged from 1 month to 11 years (mean: 74 m; median: 21 m). The median progression-free survival and overall survival times were 24 months (mean: 32 months) and 22 months (mean: 35 months), respectively. Tumor size, TNM stage and E-cadherin were found to be independent prognostic factors of TNBC.

Conclusions: EMT may play an important role in TNBC, especially in MPC and SpCC. Further researches are needed to confirm this finding. The results of this study may facilitate the future development of targeted therapies based on alterations in the EMT and stem cell markers.

Identifying global expression patterns and key regulators in epithelial to mesenchymal transition through multi-study integration

Background

Epithelial to mesenchymal transition (EMT) is the process by which stationary epithelial cells transdifferentiate to mesenchymal cells with increased motility. EMT is integral in early stages of development and wound healing. Studies have shown that EMT could be a critical early event in tumor metastasis that is involved in acquisition of migratory and invasive properties in multiple carcinomas.

Methods

In this study, we used 15 published gene expression microarray datasets from Gene Expression Omnibus (GEO) that represent 12 cell lines from 6 cancer types across 95 observations (45 unique samples and 50 replicates) with different modes of induction of EMT or the reverse transition, mesenchymal to epithelial transition (MET). We integrated multiple gene expression datasets while considering study differences, batch effects, and noise in gene expression measurements. A universal differential EMT gene list was obtained by normalizing and correcting the data using four approaches, computing differential expression from each, and identifying a consensus ranking. We confirmed our discovery of novel EMT genes at mRNA and protein levels in an in vitro EMT model of prostate cancer – PC3 epi, EMT and Taxol resistant cell lines. We validate our discovery of C1orf116 as a novel EMT regulator by siRNA knockdown of C1orf116 in PC3 epithelial cells.

Results

Among differentially expressed genes, we found known epithelial and mesenchymal marker genes such as CDH1 and ZEB1. Additionally, we discovered genes known in a subset of carcinomas that were unknown in prostate cancer. This included epithelial specific LSR and S100A14 and mesenchymal specific DPYSL3. Furthermore, we also discovered novel EMT genes including a poorly-characterized gene C1orf116. We show that decreased expression of C1orf116 is associated with poor prognosis in lung and prostate cancer patients. We demonstrate that knockdown of C1orf116 expression induced expression of mesenchymal genes in epithelial prostate cancer cell line PC3-epi cells, suggesting it as a candidate driver of the epithelial phenotype.

Conclusions

This comprehensive approach of statistical analysis and functional validation identified global expression patterns in EMT and candidate regulatory genes, thereby both extending current knowledge and identifying novel drivers of EMT.

Electronic supplementary material

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

Long non-coding RNA LOC283070 mediates the transition of LNCaP cells into androgen-independent cells possibly via CAMK1D

AIMS

The present study is to investigate the role of long non-coding RNAs (lncRNAs) in the development of androgen independence in prostate cancer and its underlying mechanism.

METHODS

We established an androgen-independent prostate carcinoma (AIPC) cell line LNCaP-AI from androgen-dependent prostate carcinoma (ADPC) cell line LNCaP. Different expression profiles of lncRNAs and mRNAs between LNCaP and LNCaP-AI cells were investigated using microarray analysis. The expression of RNAs was determined using quantitative real-time polymerase chain reaction. Protein levels were measured using Western blotting. MTT assay was used to test cell viability. Tumor formation assay was performed in nude mice to detect tumor growth in vivo. Flow cytometry was performed to detect cell cycles. Transwell assay was employed to test cell migration and invasion.

RESULTS

According to bioinformatics prediction, lncRNA LOC283070 could possibly play an important role in the transition of LNCaP cells into LNCaP-AI cells. LOC283070 was up-regulated in LNCaP-AI cells and frequently up-regulated in AIPC cell lines. Overexpression of LOC283070 in LNCaP cells accelerated cell proliferation and migration, even under androgen-independent circumstances. Knockdown of LOC283070 inhibited LNCaP-AI cell proliferation and migration. Moreover, overexpression of LOC283070 promoted tumor growth in vivo in both normal mice and castrated mice. CAMK1D overexpression had similar effect with LOC283070, and CAMK1D knockdown fully abrogated the effect of LOC283070 overexpression on the transition of LNCaP cells into androgen-independent cells.

CONCLUSIONS

The present study shows that overexpression of LOC283070 mediates the transition of LNCaP cells into androgen-independent LNCaP-AI cells possibly via CAMK1D.

Control of cell mechanics by RhoA and calcium fluxes during epithelial scattering

Epithelial tissues use adherens junctions to maintain tight interactions and coordinate cellular activities. Adherens junctions are remodeled during epithelial morphogenesis, including instances of epithelial-mesenchymal transition, or EMT, wherein individual cells detach from the tissue and migrate as individual cells. EMT has been recapitulated by growth factor induction of epithelial scattering in cell culture. In culture systems, cells undergo a highly reproducible series of cell morphology changes, most notably cell spreading followed by cellular compaction and cell migration. These morphology changes are accompanied by striking actin rearrangements. The current evidence suggests that global changes in actomyosin-based cellular contractility, first a loss of contractility during spreading and its activation during cell compaction, are the main drivers of epithelial scattering. In this review, we focus on how spreading and contractility might be controlled during epithelial scattering. While we propose a central role for RhoA, which is well known to control cellular contractility in multiple systems and whose role in epithelial scattering is well accepted, we suggest potential roles for additional cellular systems whose role in epithelial cell biology has been less well documented. In particular, we propose critical roles for vesicle recycling, calcium channels, and calcium-dependent kinases.

Intra-tumor genetic heterogeneity and alternative driver genetic alterations in breast cancers with heterogeneous HER2 gene amplification

Background

HER2 is overexpressed and amplified in approximately 15% of invasive breast cancers, and is the molecular target and predictive marker of response to anti-HER2 agents. In a subset of these cases, heterogeneous distribution of HER2 gene amplification can be found, which creates clinically challenging scenarios. Currently, breast cancers with HER2 amplification/overexpression in just over 10% of cancer cells are considered HER2-positive for clinical purposes; however, it is unclear as to whether the HER2-negative components of such tumors would be driven by distinct genetic alterations. Here we sought to characterize the pathologic and genetic features of the HER2-positive and HER2-negative components of breast cancers with heterogeneous HER2 gene amplification and to define the repertoire of potential driver genetic alterations in the HER2-negative components of these cases.

Results

We separately analyzed the HER2-negative and HER2-positive components of 12 HER2 heterogeneous breast cancers using gene copy number profiling and massively parallel sequencing, and identified potential driver genetic alterations restricted to the HER2-negative cells in each case. In vitro experiments provided functional evidence to suggest that BRF2 and DSN1 overexpression/amplification, and the HER2 I767M mutation may be alterations that compensate for the lack of HER2 amplification in the HER2-negative components of HER2 heterogeneous breast cancers.

Conclusions

Our results indicate that even driver genetic alterations, such as HER2 gene amplification, can be heterogeneously distributed within a cancer, and that the HER2-negative components are likely driven by genetic alterations not present in the HER2-positive components, including BRF2 and DSN1 amplification and HER2 somatic mutations.

Electronic supplementary material

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

The many faces of calmodulin in cell proliferation, programmed cell death, autophagy, and cancer

Calmodulin (CaM) is a ubiquitous Ca(2+) receptor protein mediating a large number of signaling processes in all eukaryotic cells. CaM plays a central role in regulating a myriad of cellular functions via interaction with multiple target proteins. This review focuses on the action of CaM and CaM-dependent signaling systems in the control of vertebrate cell proliferation, programmed cell death and autophagy. The significance of CaM and interconnected CaM-regulated systems for the physiology of cancer cells including tumor stem cells, and processes required for tumor progression such as growth, tumor-associated angiogenesis and metastasis are highlighted. Furthermore, the potential targeting of CaM-dependent signaling processes for therapeutic use is discussed.

Complex temporal changes in TGFβ oncogenic signaling drive thyroid carcinogenesis in a mouse model

Despite recent advances, understanding of molecular genetic alterations underlying thyroid carcinogenesis remains unclear. One key question is how dynamic temporal changes in global genomic expression affect carcinogenesis as the disease progresses. To address this question, we used a mouse model that spontaneously develops follicular thyroid cancer similar to human cancer (Thrb (PV/PV) mice). Using complementary DNA microarrays, we compared global gene expression profiles of thyroid tumors of Thrb (PV/PV) mice with the age- and gender-matched thyroids of wild-type mice at 3 weeks and at 2, 4, 6 and 14 months. These time points covered the pathological progression from early hyperplasia to capsular invasion, vascular invasion and eventual metastasis. Microarray data indicated that 462 genes were upregulated (Up-cluster genes) and 110 genes were downregulated (Down-cluster genes). Three major expression patterns (trending up, cyclical and spiking up and then down) and two (trending down and cyclical) were apparent in the Up-cluster and Down-cluster genes, respectively. Functional clustering of tumor-related genes followed by Ingenuity Pathways Analysis identified the transforming growth factor β (TGF β)-mediated network as key signaling pathways. Further functional analyses showed sustained activation of TGFβ receptor-pSMAD2/3 signaling, leading to decreased expression of E-cadherin and increased expression of fibronectin, vimentin, collagens and laminins. These TGFβ-induced changes facilitated epithelial-to-mesenchymal transition, which promotes cancer invasion and migration. Thus, complex temporal changes in gene expression patterns drive thyroid cancer progression, and persistent activation of TGFβ-TGFRβII-pSMAD2/3 signaling leads to EMT, thus promoting metastasis. This study provides new understanding of progression and metastatic spread of human thyroid cancer.

E47 and Id1 interplay in epithelial-mesenchymal transition

E12/E47 proteins (encoded by E2A gene) are members of the class I basic helix-loop-helix (bHLH) transcription factors (also known as E proteins). E47 has been described as repressor of E-cadherin and inducer of epithelial-mesenchymal transition (EMT). We reported previously that EMT mediated by E47 in MDCK cells occurs with a concomitant overexpression of Id1 and Id3 proteins. Id proteins belong to class V of HLH factors that lack the basic domain; they dimerise with E proteins and prevent their DNA interaction, thus, acting as dominant negative of E proteins. Here, we show that E47 interacts with Id1 in E47 overexpressing MDCK cells that underwent a full EMT as well as in mesenchymal breast carcinoma and melanoma cell lines. By conducting chromatin immunoprecipitation assays we demonstrate that E47 binds directly to the endogenous E-cadherin promoter of mesenchymal MDCK-E47 cells in a complex devoid of Id1. Importantly, our data suggest that both E47 and Id1 are required to maintain the mesenchymal phenotype of MDCK-E47 cells. These data support the collaboration between E47 and Id1 in the maintenance of EMT by mechanisms independent of the dominant negative action of Id1 on E47 binding to E-cadherin promoter. Finally, the analysis of several N0 breast tumour series indicates that the expression of E47 and ID1 is significantly associated with the basal-like phenotype supporting the biological significance of the present findings.

Digital transcriptome profiling of normal and glioblastoma-derived neural stem cells identifies genes associated with patient survival

Background

Glioblastoma multiforme, the most common type of primary brain tumor in adults, is driven by cells with neural stem (NS) cell characteristics. Using derivation methods developed for NS cells, it is possible to expand tumorigenic stem cells continuously in vitro. Although these glioblastoma-derived neural stem (GNS) cells are highly similar to normal NS cells, they harbor mutations typical of gliomas and initiate authentic tumors following orthotopic xenotransplantation. Here, we analyzed GNS and NS cell transcriptomes to identify gene expression alterations underlying the disease phenotype.

Methods

Sensitive measurements of gene expression were obtained by high-throughput sequencing of transcript tags (Tag-seq) on adherent GNS cell lines from three glioblastoma cases and two normal NS cell lines. Validation by quantitative real-time PCR was performed on 82 differentially expressed genes across a panel of 16 GNS and 6 NS cell lines. The molecular basis and prognostic relevance of expression differences were investigated by genetic characterization of GNS cells and comparison with public data for 867 glioma biopsies.

Results

Transcriptome analysis revealed major differences correlated with glioma histological grade, and identified misregulated genes of known significance in glioblastoma as well as novel candidates, including genes associated with other malignancies or glioma-related pathways. This analysis further detected several long non-coding RNAs with expression profiles similar to neighboring genes implicated in cancer. Quantitative PCR validation showed excellent agreement with Tag-seq data (median Pearson r = 0.91) and discerned a gene set robustly distinguishing GNS from NS cells across the 22 lines. These expression alterations include oncogene and tumor suppressor changes not detected by microarray profiling of tumor tissue samples, and facilitated the identification of a GNS expression signature strongly associated with patient survival (P = 1e-6, Cox model).

Conclusions

These results support the utility of GNS cell cultures as a model system for studying the molecular processes driving glioblastoma and the use of NS cells as reference controls. The association between a GNS expression signature and survival is consistent with the hypothesis that a cancer stem cell component drives tumor growth. We anticipate that analysis of normal and malignant stem cells will be an important complement to large-scale profiling of primary tumors.

Established breast cancer stem cell markers do not correlate with in vivo tumorigenicity of tumor-initiating cells

The tumor-initiating capacity of primary human breast cancer cells is maintained in vitro by culturing these cells as spheres/aggregates. Inoculation of small cell numbers derived from these non-adherent cultures leads to rapid xenograft tumor formation in mice. Accordingly, injection of more differentiated monolayer cells derived from spheres results in significantly decelerated tumor growth. For our study, two breast cancer cell lines were generated from primary tumors and cultured as mammospheres or as their adherent counterparts. We examined the in vivo tumorigenicity of these cells by injecting serial dilutions into immunodeficient mice. Inoculation of 10⁶ cells per mouse led to rapid tumor formation, irrespective of cell line or culture conditions. However, after injection of only 10³ cells, solely sphere cells were highly tumorigenic. In vitro, we investigated differentiation markers, established breast CSC markers and conducted mRNA profiling. Cytokeratin 5 and 18 were increased in both monolayer cell types, indicating a more differentiated phenotype. All cell lines were CD24⁻/CD44⁺ and did not express CD133, CD326 or E-cadherin. ALDH1 activity was not detectable in any cell line. A verapamil-sensitive Hoechst side population was present in sphere cells, but there was no correlation with tumorigenicity in vivo. mRNA profiling did not reveal upregulation of relevant transcription factors. In vitro cell cycle kinetics and in vivo tumor doubling times displayed no difference between sphere and monolayer cultures. Our data indicate that intrinsic genetic and functional markers investigated are not indicative of the in vivo tumori-genicity of putative breast tumor-initiating cells.

Basal-like Breast cancer DNA copy number losses identify genes involved in genomic instability, response to therapy, and patient survival

Breast cancer is a heterogeneous disease with known expression-defined tumor subtypes. DNA copy number studies have suggested that tumors within gene expression subtypes share similar DNA Copy number aberrations (CNA) and that CNA can be used to further sub-divide expression classes. To gain further insights into the etiologies of the intrinsic subtypes, we classified tumors according to gene expression subtype and next identified subtype-associated CNA using a novel method called SWITCHdna, using a training set of 180 tumors and a validation set of 359 tumors. Fisher's exact tests, Chi-square approximations, and Wilcoxon rank-sum tests were performed to evaluate differences in CNA by subtype. To assess the functional significance of loss of a specific chromosomal region, individual genes were knocked down by shRNA and drug sensitivity, and DNA repair foci assays performed. Most tumor subtypes exhibited specific CNA. The Basal-like subtype was the most distinct with common losses of the regions containing RB1, BRCA1, INPP4B, and the greatest overall genomic instability. One Basal-like subtype-associated CNA was loss of 5q11-35, which contains at least three genes important for BRCA1-dependent DNA repair (RAD17, RAD50, and RAP80); these genes were predominantly lost as a pair, or all three simultaneously. Loss of two or three of these genes was associated with significantly increased genomic instability and poor patient survival. RNAi knockdown of RAD17, or RAD17/RAD50, in immortalized human mammary epithelial cell lines caused increased sensitivity to a PARP inhibitor and carboplatin, and inhibited BRCA1 foci formation in response to DNA damage. These data suggest a possible genetic cause for genomic instability in Basal-like breast cancers and a biological rationale for the use of DNA repair inhibitor related therapeutics in this breast cancer subtype.

Epithelial-mesenchymal transition in breast cancer correlates with high histological grade and triple-negative phenotype

AIMS

  Epithelial-mesenchymal transition (EMT) is characterized by a loss of epithelial nature and the acquisition of a mesenchymal form. The aim of this study was to assess the role of EMT in human mammary carcinogenesis, by performing immunohistochemical studies of EMT markers with tissue microarrays.

METHODS AND RESULTS

  A total of 492 cases were evaluated and classified as hormone receptor (HR)-positive type, HER2 type and triple-negative (TN) type by the use of immunohistochemistry and in-situ hybridization. We compared these groups in terms of epithelial and mesenchymal marker expression patterns. Of the 102 cases of TN-type breast cancer, 24.5% expressed vimentin, 13.7% expressed N-cadherin, and 9.8% expressed smooth muscle actin (SMA). Of the 221 cases of HR-type breast cancer, 4.1% expressed vimentin, 5.9% expressed N-cadherin, and 0.4% expressed SMA. Regarding epithelial markers, decreased expression was seen in 16.7% of cases for E-cadherin, in 45.1% for cytokeratin (CK)19 and in 60.8% for CK8 and CK18 (CAM5.2) in TN-type breast cancer cases. Decreased expression was seen in 11.8% of cases for E-cadherin, in 6.8% for CK19 and in 3.2% for CAM5.2 in HR-type cases.

CONCLUSIONS

  EMT features were particularly seen in TN-type breast cancer (P < 0.001). EMT was also significantly associated with high histological grade (P < 0.001).

Delineation of commonly deleted chromosomal regions in meningiomas by high-density single nucleotide polymorphism genotyping arrays

Despite recent advances in the identification of the cytogenetic profiles of meningiomas, a significant group of tumors still show normal karyotypes or few chromosomal changes. The authors analyzed the cytogenetic profile of 50 meningiomas using fluorescence in situ hybridization and high-density (500 K) single nucleotide polymorphism (SNP) arrays. Our results confirm that del(22q) (52%) and del(1p) (16%) (common deleted regions: 22q11.21-22q13.3. and 1p31.2-p36.33) are the most frequent alterations. Additionally, recurrent monosomy 14 (8%), del(6q) (10%), del(7p) (10%), and del(19q) (4%) were observed, while copy number patterns consistent with recurrent chromosomal gains, gene amplification, and copy number neutral loss of heterozygosity (cnLOH) were either absent or rare. Based on their overall SNP profiles, meningiomas could be classified into: (i) diploid cases, (ii) meningiomas with a single chromosomal change [e.g., monosomy 22/del(22q)] and (iii) tumors with ≥2 altered chromosomes. In summary, our results confirm and extend on previous observations showing that the most recurrent chromosomal abnormalities in meningiomas correspond to chromosome losses localized in chromosomes 1, 22 and less frequently in chromosomes 6, 7, 14, and 19, while chromosomal gains and cnLOH are restricted to a small proportion of cases. Finally, a set of cancer-associated candidate genes associated with the TP53, MYC, CASP3, HDAC1, and TERT signaling pathways was identified, in cases with coexisting monosomy 14 and del(1p).

Cdk9 T-loop phosphorylation is regulated by the calcium signaling pathway

Eukaryotic RNA polymerase II transcriptional elongation is a tightly regulated process and is dependent upon positive transcription elongation factor-b (P-TEFb). The core P-TEFb complex is composed of Cdk9 and Cyclin T and is essential for the expression of most protein coding genes. Cdk9 kinase function is dependent upon phosphorylation of Thr186 in its T-loop. In this study, we examined kinases and signaling pathways that influence Cdk9 T-loop phosphorylation. Using an RNAi screen in HeLa cells, we found that Cdk9 T-loop phosphorylation is regulated by Ca(2+)/calmodulin-dependent kinase 1D (CaMK1D). Using small molecules inhibitors in HeLa cells and primary CD4(+) T lymphocytes, we found that the Ca(2+) signaling pathway is required for Cdk9 T-loop phosphorylation. Inhibition of Ca(2+) signaling led to dephosphorylation of Thr186 on Cdk9. In reporter plasmid assays, inhibition of the Ca(2+) signaling pathway repressed the PCNA promoter and HIV-1 Tat transactivation of the HIV-1 LTR, but not HTLV-1 Tax transactivation of the HTLV-1 LTR, suggesting that perturbation of the Ca(2+) pathway and reduction of Cdk9 T-loop phosphorylation inhibits transcription units that have a rigorous requirement for P-TEFb function.

Lessons from a decade of integrating cancer copy number alterations with gene expression profiles

Over the last decade, multiple functional genomic datasets studying chromosomal aberrations and their downstream effects on gene expression have accumulated for several cancer types. A vast majority of them are in the form of paired gene expression profiles and somatic copy number alterations (CNA) information on the same patients identified using microarray platforms. In response, many algorithms and software packages are available for integrating these paired data. Surprisingly, there has been no serious attempt to review the currently available methodologies or the novel insights brought using them. In this work, we discuss the quantitative relationships observed between CNA and gene expression in multiple cancer types and biological milestones achieved using the available methodologies. We discuss the conceptual evolution of both, the step-wise and the joint data integration methodologies over the last decade. We conclude by providing suggestions for building efficient data integration methodologies and asking further biological questions.

Genome-wide copy number alterations in subtypes of invasive breast cancers in young white and African American women

Genomic copy number alterations (CNA) are common in breast cancer. Identifying characteristic CNAs associated with specific breast cancer subtypes is a critical step in defining potential mechanisms of disease initiation and progression. We used genome-wide array comparative genomic hybridization to identify distinctive CNAs in breast cancer subtypes from 259 young (diagnosed with breast cancer at <55 years) African American (AA) and Caucasian American (CA) women originally enrolled in a larger population-based study. We compared the average frequency of CNAs across the whole genome for each breast tumor subtype and found that estrogen receptor (ER)-negative tumors had a higher average frequency of genome-wide gain (P < 0.0001) and loss (P = 0.02) compared to ER-positive tumors. Triple-negative (TN) tumors had a higher average frequency of genome-wide gain (P < 0.0001) and loss (P = 0.003) than non-TN tumors. No significant difference in CNA frequency was observed between HER2-positive and -negative tumors. We also identified previously unreported recurrent CNAs (frequency >40%) for TN breast tumors at 10q, 11p, 11q, 16q, 20p, and 20q. In addition, we report CNAs that differ in frequency between TN breast tumors of AA and CA women. This is of particular relevance because TN breast cancer is associated with higher mortality and young AA women have higher rates of TN breast tumors compared to CA women. These data support the possibility that higher overall frequency of genomic alteration events as well as specific focal CNAs in TN breast tumors might contribute in part to the poor breast cancer prognosis for young AA women.

Common variants associated with breast cancer in genome-wide association studies are modifiers of breast cancer risk in BRCA1 and BRCA2 mutation carriers

Recent studies have identified single nucleotide polymorphisms (SNPs) that significantly modify breast cancer risk in BRCA1 and BRCA2 mutation carriers. Since these risk modifiers were originally identified as genetic risk factors for breast cancer in genome-wide association studies (GWASs), additional risk modifiers for BRCA1 and BRCA2 may be identified from promising signals discovered in breast cancer GWAS. A total of 350 SNPs identified as candidate breast cancer risk factors (P < 1 x 10(-3)) in two breast cancer GWAS studies were genotyped in 3451 BRCA1 and 2006 BRCA2 mutation carriers from nine centers. Associations with breast cancer risk were assessed using Cox models weighted for penetrance. Eight SNPs in BRCA1 carriers and 12 SNPs in BRCA2 carriers, representing an enrichment over the number expected, were significantly associated with breast cancer risk (P(trend) < 0.01). The minor alleles of rs6138178 in SNRPB and rs6602595 in CAMK1D displayed the strongest associations in BRCA1 carriers (HR = 0.78, 95% CI: 0.69-0.90, P(trend) = 3.6 x 10(-4) and HR = 1.25, 95% CI: 1.10-1.41, P(trend) = 4.2 x 10(-4)), whereas rs9393597 in LOC134997 and rs12652447 in FBXL7 showed the strongest associations in BRCA2 carriers (HR = 1.55, 95% CI: 1.25-1.92, P(trend) = 6 x 10(-5) and HR = 1.37, 95% CI: 1.16-1.62, P(trend) = 1.7 x 10(-4)). The magnitude and direction of the associations were consistent with the original GWAS. In subsequent risk assessment studies, the loci appeared to interact multiplicatively for breast cancer risk in BRCA1 and BRCA2 carriers. Promising candidate SNPs from GWAS were identified as modifiers of breast cancer risk in BRCA1 and BRCA2 carriers. Upon further validation, these SNPs together with other genetic and environmental factors may improve breast cancer risk assessment in these populations.

Genomic instability in breast cancer: pathogenesis and clinical implications

Breast cancer is a heterogeneous disease, appreciable by molecular markers, gene-expression profiles, and most recently, patterns of genomic alteration. In particular, genomic profiling has revealed three distinct patterns of DNA copy-number alteration: a "simple" type with few gains or losses of whole chromosome arms, an "amplifier" type with focal high-level DNA amplifications, and a "complex" type marked by numerous low-amplitude changes and copy-number transitions. The three patterns are associated with distinct gene-expression subtypes, and preferentially target different loci in the genome (implicating distinct cancer genes). Moreover, the different patterns of alteration imply distinct underlying mechanisms of genomic instability. The amplifier pattern may arise from transient telomere dysfunction, although new data suggest ongoing "amplifier" instability. The complex pattern shows similarity to breast cancers with germline BRCA1 mutation, which also exhibit "basal-like" expression profiles and complex-pattern genomes, implicating a possible defect in BRCA1-associated repair of DNA double-strand breaks. As such, targeting presumptive DNA repair defects represents a promising area of clinical investigation. Future studies should clarify the pathogenesis of breast cancers with amplifier and complex-pattern genomes, and will likely identify new therapeutic opportunities.

LYN is a mediator of epithelial-mesenchymal transition and a target of dasatinib in breast cancer

Epithelial-mesenchymal transition (EMT), a switch of polarized epithelial cells to a migratory, fibroblastoid phenotype, is considered a key process driving tumor cell invasiveness and metastasis. Using breast cancer cell lines as a model system, we sought to discover gene expression signatures of EMT with clinical and mechanistic relevance. A supervised comparison of epithelial and mesenchymal breast cancer lines defined a 200-gene EMT signature that was prognostic across multiple breast cancer cohorts. The immunostaining of LYN, a top-ranked EMT signature gene and Src-family tyrosine kinase, was associated with significantly shorter overall survival (P = 0.02) and correlated with the basal-like ("triple-negative") phenotype. In mesenchymal breast cancer lines, RNAi-mediated knockdown of LYN inhibited cell migration and invasion, but not proliferation. Dasatinib, a dual-specificity tyrosine kinase inhibitor, also blocked invasion (but not proliferation) at nanomolar concentrations that inhibit LYN kinase activity, suggesting that LYN is a likely target and that invasion is a relevant end point for dasatinib therapy. Our findings define a prognostically relevant EMT signature in breast cancer and identify LYN as a mediator of invasion and a possible new therapeutic target (and theranostic marker for dasatinib response), with particular relevance to clinically aggressive basal-like breast cancer.

Regularized Multivariate Regression for Identifying Master Predictors with Application to Integrative Genomics Study of Breast Cancer

In this paper, we propose a new method remMap - REgularized Multivariate regression for identifying MAster Predictors - for fitting multivariate response regression models under the high-dimension-low-sample-size setting. remMap is motivated by investigating the regulatory relationships among different biological molecules based on multiple types of high dimensional genomic data. Particularly, we are interested in studying the influence of DNA copy number alterations on RNA transcript levels. For this purpose, we model the dependence of the RNA expression levels on DNA copy numbers through multivariate linear regressions and utilize proper regularization to deal with the high dimensionality as well as to incorporate desired network structures. Criteria for selecting the tuning parameters are also discussed. The performance of the proposed method is illustrated through extensive simulation studies. Finally, remMap is applied to a breast cancer study, in which genome wide RNA transcript levels and DNA copy numbers were measured for 172 tumor samples. We identify a trans-hub region in cytoband 17q12-q21, whose amplification influences the RNA expression levels of more than 30 unlinked genes. These findings may lead to a better understanding of breast cancer pathology.

Focal amplification and oncogene dependency of GAB2 in breast cancer

DNA amplifications in breast cancer are frequent on chromosome 11q, in which multiple driver oncogenes likely reside in addition to cyclin D1 (CCND1). One such candidate, the scaffolding adapter protein, GRB2-associated binding protein 2 (GAB2), functions in ErbB signaling and was recently shown to enhance mammary epithelial cell proliferation, and metastasis of ERBB2 (HER2/neu)-driven murine breast cancer. However, the amplification status and function of GAB2 in the context of amplification remain undefined. In this study, by genomic profiling of 172 breast tumors, and fluorescence in situ hybridization validation in an independent set of 210 scorable cases, we observed focal amplification spanning GAB2 (11q14.1) independent of CCND1 (11q13.2) amplification, consistent with a driver role. Further, small interfering RNA (siRNA)-mediated knockdown of GAB2 in breast cancer lines (SUM52, SUM44PE and MDA468) with GAB2 amplification revealed a dependency on GAB2 for cell proliferation, cell-cycle progression, survival and invasion, likely mediated through altered phosphatidylinositol 3-kinase (PI3K) and mitogen-activated protein kinase (MAPK) signaling. GAB2 knockdown also reduced proliferation and survival in a cell line (BT474) with ERBB2 amplification, consistent with the possibility that GAB2 can function downstream of ERBB2. Our studies implicate focal amplification of GAB2 in breast carcinogenesis, and underscore an oncogenic role of scaffolding adapter proteins, and a potential new point of therapeutic intervention.

Molecular profiling of breast cancer cell lines defines relevant tumor models and provides a resource for cancer gene discovery

Background

Breast cancer cell lines have been used widely to investigate breast cancer pathobiology and new therapies. Breast cancer is a molecularly heterogeneous disease, and it is important to understand how well and which cell lines best model that diversity. In particular, microarray studies have identified molecular subtypes–luminal A, luminal B, ERBB2-associated, basal-like and normal-like–with characteristic gene-expression patterns and underlying DNA copy number alterations (CNAs). Here, we studied a collection of breast cancer cell lines to catalog molecular profiles and to assess their relation to breast cancer subtypes.

Methods

Whole-genome DNA microarrays were used to profile gene expression and CNAs in a collection of 52 widely-used breast cancer cell lines, and comparisons were made to existing profiles of primary breast tumors. Hierarchical clustering was used to identify gene-expression subtypes, and Gene Set Enrichment Analysis (GSEA) to discover biological features of those subtypes. Genomic and transcriptional profiles were integrated to discover within high-amplitude CNAs candidate cancer genes with coordinately altered gene copy number and expression.

Findings

Transcriptional profiling of breast cancer cell lines identified one luminal and two basal-like (A and B) subtypes. Luminal lines displayed an estrogen receptor (ER) signature and resembled luminal-A/B tumors, basal-A lines were associated with ETS-pathway and BRCA1 signatures and resembled basal-like tumors, and basal-B lines displayed mesenchymal and stem/progenitor-cell characteristics. Compared to tumors, cell lines exhibited similar patterns of CNA, but an overall higher complexity of CNA (genetically simple luminal-A tumors were not represented), and only partial conservation of subtype-specific CNAs. We identified 80 high-level DNA amplifications and 13 multi-copy deletions, and the resident genes with concomitantly altered gene-expression, highlighting known and novel candidate breast cancer genes.

Conclusions

Overall, breast cancer cell lines were genetically more complex than tumors, but retained expression patterns with relevance to the luminal-basal subtype distinction. The compendium of molecular profiles defines cell lines suitable for investigations of subtype-specific pathobiology, cancer stem cell biology, biomarkers and therapies, and provides a resource for discovery of new breast cancer genes.

CAMK1D amplification implicated in epithelial-mesenchymal transition in basal-like breast cancer

Breast cancer exhibits clinical and molecular heterogeneity, where expression profiling studies have identified five major molecular subtypes. The basal-like subtype, expressing basal epithelial markers and negative for estrogen receptor (ER), progesterone receptor (PR) and HER2, is associated with higher overall levels of DNA copy number alteration (CNA), specific CNAs (like gain on chromosome 10p), and poor prognosis. Discovering the molecular genetic basis of tumor subtypes may provide new opportunities for therapy. To identify the driver oncogene on 10p associated with basal-like tumors, we analyzed genomic profiles of 172 breast carcinomas. The smallest shared region of gain spanned just seven genes at 10p13, including calcium/calmodulin-dependent protein kinase ID (CAMK1D), functioning in intracellular signaling but not previously linked to cancer. By microarray, CAMK1D was overexpressed when amplified, and by immunohistochemistry exhibited elevated expression in invasive carcinomas compared to carcinoma in situ. Engineered overexpression of CAMK1D in non-tumorigenic breast epithelial cells led to increased cell proliferation, and molecular and phenotypic alterations indicative of epithelial-mesenchymal transition (EMT), including loss of cell-cell adhesions and increased cell migration and invasion. Our findings identify CAMK1D as a novel amplified oncogene linked to EMT in breast cancer, and as a potential therapeutic target with particular relevance to clinically unfavorable basal-like tumors.