BIRC5 Genomic Copy Number Variation in Early-Onset Breast Cancer

Determining expression changes of ANO7 and SLC38A4 membrane transporters in colorectal cancer

Membrane transporters are proteins responsible for facilitating the movement of molecules within biological membranes. They play a vital role in maintaining cellular homeostasis by regulating the transport of nutrients, ions, and other molecules into and out of cells. Our aim is to identify biomarkers in colorectal cancer using membrane transporter proteins. We utilized COAD TCGA data for this purpose. Subsequently, we conducted differential gene analysis and feature selection using membrane transporter proteins. Furthermore, we identified two potential genes, including ANO7 and SLC38A4. To validate the expression profiles of ANO7 and SLC38A4, key genes in this context, RT-qPCR was employed on colorectal cancer samples and adjacent normal tissues. Additionally, utilizing GEPIA2, Kaplan-Meier survival analysis, and cBioPortal, we assessed the status of these genes in various cancers, examining their methylation and mutation patterns. In conclusion, we suggest that ANO7 and SLC38A4 serve as prognostic biomarkers in colorectal cancer.

Identification and verification of three autophagy-related genes as potential biomarkers for the diagnosis of psoriasis

Psoriasis vulgaris is the most common form of the four clinical types. However, early diagnosis of psoriasis vulgaris is difficult due to the lack of effective biomarkers. The aim of this study was to screen potential biomarkers for the diagnosis of psoriasis. In our study, we downloaded the original data from GSE30999 and GSE41664, and the autophagy-related genes list from human autophagy database to identify differentially expressed autophagy-related genes (DERAGs) by R software. Then Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analysis were performed for DERAGs. DERAGs were validated by the other four databases (GSE13355, GSE14905, GSE6710, and GSE55201) to screen biomarkers with high diagnostic value for the early diagnosis of psoriasis vulgaris. Finally, DERAGs were verified in our clinical blood samples by ELISA. A total of 12 DERAGs were identified between 123 paired non-lesional and lesional skin samples from patients with psoriasis vulgaris. GO and KEGG enrichment analysis indicated the TORC2 complex was more enriched and the NOD-like receptor signaling pathway was mostly enriched. Three autophagy-related genes (BIRC5, NAMPT and BCL2) were identified through bioinformatics analysis and verified by ELISA in clinical blood samples. And these genes showed high diagnostic value for the early diagnosis of psoriasis vulgaris. We identified three autophagy-related genes (BIRC5, NAMPT and BCL2) with high diagnostic value for the early diagnosis of psoriasis vulgaris through bioinformatics analysis and clinical samples. Therefore, we proposed that BIRC5, NAMPT and BCL2 may be as potential biomarkers for the early diagnosis of psoriasis vulgaris. In addition, BIRC5, NAMPT and BCL2 may affect the development of psoriasis by regulating autophagy.

Multi-omics characteristics of tumor-associated macrophages in the tumor microenvironment of gastric cancer and their exploration of immunotherapy potential

The incidence and mortality rate of gastric cancer (GC) have remained high worldwide. Although some progress has been made in immunotargeted therapy, the treatment effect remains limited. With more attention has been paid to the immune potential of tumor-associated macrophages (TAMs), but the specific mechanisms of tumor immunity are still unclear. Thus, we screened marker genes in TAMs differentiation (MDMs) through single-cell RNA sequencing, and combined with GC transcriptome data from TCGA and GEO databases, the clinical and TME characteristics, prognostic differences, immune infiltration, and drug sensitivity among different subtypes of patients with GC in different data sets were analyzed. A prognostic model of GC was constructed to evaluate the prognosis and immunotherapy response of patients with GC. In this study, we extensively studied the mutations in MDMs such as CGN, S100A6, and C1QA, and found differences in the infiltration of immune cells and immune checkpoints including M2 TAMs, T cells, CD274, and CTLA4 in different GC subtypes. In the model, we constructed a predictive scoring system with high accuracy and screened out key MDMs-related genes associated with prognosis and M2 TAMs, among which VKORC1 may be involved in GC progression and iron death in tumor cells. Therefore, this study explores the therapeutic strategy of TAMs reprogramming in-depth, providing new ideas for the clinical diagnosis, treatment, and prognosis assessment of GC.

Single-cell and genetic multi-omics analysis combined with experiments confirmed the signature and potential targets of cuproptosis in hepatocellular carcinoma

Background: Cuproptosis, as a recently discovered type of programmed cell death, occupies a very important role in hepatocellular carcinoma (HCC) and provides new methods for immunotherapy; however, the functions of cuproptosis in HCC are still unclear. Methods: We first analyzed the transcriptome data and clinical information of 526 HCC patients using multiple algorithms in R language and extensively described the copy number variation, prognostic and immune infiltration characteristics of cuproptosis related genes (CRGs). Then, the hub CRG related genes associated with prognosis through LASSO and Cox regression analyses and constructed a prognostic prediction model including multiple molecular markers and clinicopathological parameters through training cohorts, then this model was verified by test cohorts. On the basis of the model, the clinicopathological indicators, immune infiltration and tumor microenvironment characteristics of HCC patients were further explored via bioinformation analysis. Then, We further explored the key gene biological function by single-cell analysis, cell viability and transwell experiments. Meantime, we also explored the molecular docking of the hub genes. Results: We have screened 5 hub genes associated with HCC prognosis and constructed a prognosis prediction scoring model. And the model results showed that patients in the high-risk group had poor prognosis and the expression levels of multiple immune markers, including PD-L1, CD276 and CTLA4, were higher than those patients in the low-risk group. We found a significant correlation between risk score and M0 macrophages and memory CD4+ T cells. And the single-cell analysis and molecular experiments showed that BEX1 were higher expressed in HCC tissues and deletion inhibited the proliferation, invasion and migration and EMT pathway of HCC cells. Finally, it was observed that BEX1 could bind to sorafenib to form a stable conformation. Conclusion: The study not only revealed the multiomics characteristics of CRGs in HCC but also constructed a new high-accuracy prognostic prediction model. Meanwhile, BEX1 were also identified as hub genes that can mediate the cuproptosis of hepatocytes as potential therapeutic targets for HCC.

Bioinformatics and network-based screening and discovery of potential molecular targets and small molecular drugs for breast cancer

Accurate identification of molecular targets of disease plays an important role in diagnosis, prognosis, and therapies. Breast cancer (BC) is one of the most common malignant cancers in women worldwide. Thus, the objective of this study was to accurately identify a set of molecular targets and small molecular drugs that might be effective for BC diagnosis, prognosis, and therapies, by using existing bioinformatics and network-based approaches. Nine gene expression profiles (GSE54002, GSE29431, GSE124646, GSE42568, GSE45827, GSE10810, GSE65216, GSE36295, and GSE109169) collected from the Gene Expression Omnibus (GEO) database were used for bioinformatics analysis in this study. Two packages, LIMMA and clusterProfiler, in R were used to identify overlapping differential expressed genes (oDEGs) and significant GO and KEGG enrichment terms. We constructed a PPI (protein–protein interaction) network through the STRING database and identified eight key genes (KGs) EGFR, FN1, EZH2, MET, CDK1, AURKA, TOP2A, and BIRC5 by using six topological measures, betweenness, closeness, eccentricity, degree, MCC, and MNC, in the Analyze Network tool in Cytoscape. Three online databases GSCALite, Network Analyst, and GEPIA were used to analyze drug enrichment, regulatory interaction networks, and gene expression levels of KGs. We checked the prognostic power of KGs through the prediction model using the popular machine learning algorithm support vector machine (SVM). We suggested four TFs (TP63, MYC, SOX2, and KDM5B) and four miRNAs (hsa-mir-16-5p, hsa-mir-34a-5p, hsa-mir-1-3p, and hsa-mir-23b-3p) as key transcriptional and posttranscriptional regulators of KGs. Finally, we proposed 16 candidate repurposing drugs YM201636, masitinib, SB590885, GSK1070916, GSK2126458, ZSTK474, dasatinib, fedratinib, dabrafenib, methotrexate, trametinib, tubastatin A, BIX02189, CP466722, afatinib, and belinostat for BC through molecular docking analysis. Using BC cell lines, we validated that masitinib inhibits the mTOR signaling pathway and induces apoptotic cell death. Therefore, the proposed results might play an effective role in the treatment of BC patients.

The Unique Biology behind the Early Onset of Breast Cancer

Breast cancer commonly affects women of older age; however, in developing countries, up to 20% of breast cancer cases present in young women (younger than 40 years as defined by oncology literature). Breast cancer in young women is often defined to be aggressive in nature, usually of high histological grade at the time of diagnosis and negative for endocrine receptors with poor overall survival rate. Several researchers have attributed this aggressive nature to a hidden unique biology. However, findings in this aspect remain controversial. Thus, in this article, we aimed to review published work addressing somatic mutations, chromosome copy number variants, single nucleotide polymorphisms, differential gene expression, microRNAs and gene methylation profile of early-onset breast cancer, as well as its altered pathways resulting from those aberrations. Distinct biology behind early-onset of breast cancer was clear among estrogen receptor-positive and sporadic cases. However, further research is needed to determine and validate specific novel markers, which may help in customizing therapy for this group of patients.

Molecular Cytogenetic Characterization of the Murine Melanoma Cell Lines S91 Clone M3 and B16-F1 with Variant B16-4A5

Melanoma is considered to be one of the most aggressive human tumors. Thus, early molecular diagnosis with risk factor stratification could be an efficacious strategy to increase the survival rates in affected patients. Murine cell lines B16-F1, B16-4A5, and S91 clone M3 are the ones most commonly applied in melanoma research. However, genetic peculiarities of these 3 cell lines have not been studied in detail before. Here, we closed this gap by molecular cytogenetic and array-comparative genomic hybridization studies and the translation of the characterized imbalances into the human genome. This study revealed severely rearranged karyotypes with in parts similar imbalances for all 3 cell lines. Interestingly, they involve genes known to play major roles in human melanoma. These are specifically the oncogenes and tumor suppressor genes, being associated with aggressive forms of melanoma. B16-F1, B16-4A5, and S91 clone M3 revealed aberrations which were similarly observed in human eye and skin but not in human uveal melanoma. Thus, they can be considered as model systems for advanced eye and skin melanoma.

Identification of key genes and pathways in endometriosis by integrated expression profiles analysis

The purpose of this study was to integrate the existing expression profile data on endometriosis (EM)-related tissues in order to identify the differentially expressed genes. In this study, three series of raw expression data were downloaded from GEO database. Differentially expressed genes (DEGs) in three tissue types were screened. GO, KEGG pathway enrichment analysis, core differential genes (CDGs) protein–protein interaction (PPI) network and weighted gene co-expression network analysis (WGCNA) were performed, finally, the dysregulation of Hippo pathway in ectopic endometrium (EC) was detected by Western blotting. A total of 1,811 DEGs between eutopic (EU) and normal endometrium (NE), 5,947 DEGs between EC and EU, and 3,192 DEGs between EC and NE datasets were identified. After screening, 394 CDGs were obtained, and 5 hub genes identified in the PPI network. CDGs enrichment and WGCNA network analysis revealed cell proliferation, differentiation, migration and other biological processes, Hippo and Wnt signaling pathways, and a variety of tumor-related pathways. Western blotting results showed that YAP/TAZ was upregulated, and MOB1, pMOB1, SAV1, LATS1 and LATS2 were downregulated in EC. Moreover, CDGs, especially the hub genes, are potential biomarkers and therapeutic targets. Finally, the Hippo pathway might play a key role in the development of endometriosis.

Gene Expression Alterations Associated with Oleuropein-Induced Antiproliferative Effects and S-Phase Cell Cycle Arrest in Triple-Negative Breast Cancer Cells

It is known that the Mediterranean diet is effective in reducing the risk of several chronic diseases, including cancer. A critical component of the Mediterranean diet is olive oil, and the relationship between olive oil consumption and the reduced risk of cancer has been established. Oleuropein (OL) is the most prominent polyphenol component of olive fruits and leaves. This compound has been shown to have potent properties in various types of cancers, including breast cancer. In the present study, the molecular mechanism of OL was examined in two racially different triple-negative breast cancer (TNBC) cell lines-African American (AA, MDA-MB-468) and Caucasian American (CA, MDA-MB-231). The data obtained showed that OL effectively inhibits cell growth in both cell lines, concomitant with S-phase cell cycle arrest-mediated apoptosis. The results also showed that OL-treated MDA-MB-468 cells were two-fold more sensitive to OL antiproliferative effect than MDA-MB-231 cells were. At lower concentrations, OL modified the expression of many apoptosis-involved genes. OL was more effective in MDA-MB-468, compared to MDA-MB-231 cells, in terms of the number and the fold-change of the altered genes. In MDA-MB-468 cells, OL induced a noticeable transcription activation in fourteen genes, including two members of the caspase family: caspase 1 (CASP1) and caspase 14 (CASP14); two members of the TNF receptor superfamily: Fas-associated via death domain (FADD) and TNF receptor superfamily 21 (TNFRSF21); six other proapoptotic genes: growth arrest and DNA damage-inducible 45 alpha (GADD45A), cytochrome c somatic (CYCS), BCL-2 interacting protein 2 (BNIP2), BCL-2 interacting protein 3 (BNIP3), BH3 interacting domain death agonist (BID), and B-cell lymphoma/leukemia 10 (BCL10); and the CASP8 and FADD-like apoptosis regulator (CFLAR) gene. Moreover, in MDA-MB-468 cells, OL induced a significant upregulation in two antiapoptotic genes: bifunctional apoptosis regulator (BFAR) and B-Raf proto-oncogene (BRAF) and a baculoviral inhibitor of apoptosis (IAP) repeat-containing 3 (BIRC3). On the contrary, in MDA-MB-231 cells, OL showed mixed impacts on gene expression. OL significantly upregulated the mRNA expression of four genes: BIRC3, receptor-interacting serine/threonine kinase 2 (RIPK2), TNF receptor superfamily 10A (TNFRSF10A), and caspase 4 (CASP4). Additionally, another four genes were repressed, including caspase 6 (CASP6), pyrin domain (PYD), and caspase recruitment domain (CARD)-containing (PAYCARD), baculoviral IAP repeat-containing 5 (BIRC5), and the most downregulated TNF receptor superfamily member 11B (TNFRSF11B, 16.34-fold). In conclusion, the data obtained indicate that the two cell lines were markedly different in the anticancer effect and mechanisms of oleuropein's ability to alter apoptosis-related gene expressions. The results obtained from this study should also guide the potential utilization of oleuropein as an adjunct therapy for TNBC to increase chemotherapy effectiveness and prevent cancer progression.

Integrated Genomic Characterization of the Human Immunome in Cancer

Alterations in immune-related pathways are common hallmarks of cancer. A comprehensive understanding of how cancer mutations rewire immune signaling networks and functional output across cancer types is instrumental to realize the full potential of immunotherapy. Here, we systematically interrogated somatic mutations involved in immune signaling that alter immune responses in patients with cancer. To do so, we developed a Network-based Integrative model to Prioritize Potential immune respondER genes (NIPPER). Identified mutations were enriched in essential protein domains and genes identified by NIPPER were associated with responsiveness to multiple immunotherapy modalities. These genes were used to devise an interactome network propagation framework integrated with drug-associated gene signatures to identify potential immunomodulatory drug candidates. Together, our systems-level analysis results help interpret the heterogeneous immune responses among patients and serve as a resource for future functional studies and targeted therapeutics. SIGNIFICANCE: This study demonstrates that integration of multi-omics data can help identify critical molecular determinants for effective targeted therapeutics.

Hybrid Stem Cell States: Insights Into the Relationship Between Mammary Development and Breast Cancer Using Single-Cell Transcriptomics

Similarities between stem cells and cancer cells have implicated mammary stem cells in breast carcinogenesis. Recent evidence suggests that normal breast stem cells exist in multiple phenotypic states: epithelial, mesenchymal, and hybrid epithelial/mesenchymal (E/M). Hybrid E/M cells in particular have been implicated in breast cancer metastasis and poor prognosis. Mounting evidence also suggests that stem cell phenotypes change throughout the life course, for example, through embryonic development and pregnancy. The goal of this study was to use single cell RNA-sequencing to quantify cell state distributions of the normal mammary (NM) gland throughout developmental stages and when perturbed into a stem-like state in vitro using conditional reprogramming (CR). Using machine learning based dataset alignment, we integrate multiple mammary gland single cell RNA-seq datasets from human and mouse, along with bulk RNA-seq data from breast tumors in the Cancer Genome Atlas (TCGA), to interrogate hybrid stem cell states in the normal mammary gland and cancer. CR of human mammary cells induces an expanded stem cell state, characterized by increased expression of embryonic stem cell associated genes. Alignment to a mouse single-cell transcriptome atlas spanning mammary gland development from in utero to adulthood revealed that NM cells align to adult mouse cells and CR cells align across the pseudotime trajectory with a stem-like population aligning to the embryonic mouse cells. Three hybrid populations emerge after CR that are rare in NM: KRT18+/KRT14+ (hybrid luminal/basal), EPCAM+/VIM+ (hybrid E/M), and a quadruple positive population, expressing all four markers. Pseudotime analysis and alignment to the mouse developmental trajectory revealed that E/M hybrids are the most developmentally immature. Analyses of single cell mouse mammary RNA-seq throughout pregnancy show that during gestation, there is an enrichment of hybrid E/M cells, suggesting that these cells play an important role in mammary morphogenesis during lactation. Finally, pseudotime analysis and alignment of TCGA breast cancer expression data revealed that breast cancer subtypes express distinct developmental signatures, with basal tumors representing the most “developmentally immature” phenotype. These results highlight phenotypic plasticity of normal mammary stem cells and provide insight into the relationship between hybrid cell populations, stemness, and cancer.

Runt-related transcription factor 1 promotes apoptosis and inhibits neuroblastoma progression in vitro and in vivo

Background

Runt-related transcription factor 1 (RUNX1) is a heterodimeric transcription factor that binds to the core element of many enhancers and promoters and can accelerate apoptosis in various tumors. However, the regulatory mechanisms underlying RUNX1 expression in neuroblastoma (NB), a highly malignant tumor in childhood, remain largely unclear. In this study, we aimed to assess the role of RUNX1 in NB and to reveal the underlying mechanisms that may contribute to finding a potential therapeutics strategy against NB.

Methods

Growth, invasion, metastasis and angiogenesis were assessed using Cell Counting Kit-8 (CCK-8) immunocytochemistry, and studies involving soft agar, cell invasion, tube formation and whole animals. The levels of expression were measured using real-time quantitative PCR for RNA, Western blot and immunostaining analyses for proteins. Luciferase reporter and chromatin immunoprecipitation assays indicated that RUNX1 directly binds within the BIRC5, CSF2RB and NFKBIA promoter regions to facilitate transcription. The level of apoptosis was assessed by determining mitochondrial membrane potential and flow cytometry.

Results

RUNX1 was highly expressed in ganglioneuroma (GN) and well-differentiated (WD) tissues relative to the poorly differentiated (PD) and undifferentiated (UD) ones. Moreover, RUNX1 effectively reduced cell viability, invasion, metastasis, angiogenesis, and promoted apoptosis in vitro and in vivo. RUNX1 reduced BIRC5 transcription and increased CSF2RB and NFKBIA transcription by directly binding BIRC5, CSF2RB and NFKBIA promoters. In addition, cytotoxic drugs, especially cisplatin, significantly increased RUNX1 expression in NB cells and promoted apoptosis.

Conclusions

These data show that RUNX1 is an independent surrogate marker for the progression of NB and it can be used for monitoring NB prognosis during therapy.

Integrated bioinformatics analysis for the screening of hub genes and therapeutic drugs in ovarian cancer

Background

Ovarian cancer (OC) ranks fifth as a cause of gynecological cancer-associated death globally. Until now, the molecular mechanisms underlying the tumorigenesis and prognosis of OC have not been fully understood. This study aims to identify hub genes and therapeutic drugs involved in OC.

Methods

Four gene expression profiles (GSE54388, GSE69428, GSE36668, and GSE40595) were downloaded from the Gene Expression Omnibus (GEO), and the differentially expressed genes (DEGs) in OC tissues and normal tissues with an adjusted P-value < 0.05 and a |log fold change (FC)| > 1.0 were first identified by GEO2R and FunRich software. Next, Gene Ontology (GO) and Kyoto Encyclopaedia of Genes and Genomes (KEGG) analyses were performed for functional enrichment analysis of these DEGs. Then, the hub genes were identified by the cytoHubba plugin and the other bioinformatics approaches including protein-protein interaction (PPI) network analysis, module analysis, survival analysis, and miRNA-hub gene network construction was also performed. Finally, the GEPIA2 and DGIdb databases were utilized to verify the expression levels of hub genes and to select the candidate drugs for OC, respectively.

Results

A total of 171 DEGs were identified, including 114 upregulated and 57 downregulated DEGs. The results of the GO analysis indicated that the upregulated DEGs were mainly involved in cell division, nucleus, and protein binding, whereas the biological functions showing enrichment in the downregulated DEGs were mainly negative regulation of transcription from RNA polymerase II promoter, protein complex and apicolateral plasma membrane, and glycosaminoglycan binding. As for the KEGG-pathway, the upregulated DEGs were mainly associated with metabolic pathways, biosynthesis of antibiotics, biosynthesis of amino acids, cell cycle, and HTLV-I infection. Additionally, 10 hub genes (KIF4A, CDC20, CCNB2, TOP2A, RRM2, TYMS, KIF11, BIRC5, BUB1B, and FOXM1) were identified and survival analysis of these hub genes showed that OC patients with the high-expression of CCNB2, TYMS, KIF11, KIF4A, BIRC5, BUB1B, FOXM1, and CDC20 were statistically more likely to have poorer progression free survival. Meanwhile, the expression levels of the hub genes based on GEPIA2 were in accordance with those based on GEO. Finally, DGIdb database was used to identify 62 small molecules as the potentially targeted drugs for OC treatment.

Conclusions

In summary, the data may produce new insights regarding OC pathogenesis and treatment. Hub genes and candidate drugs may improve individualized diagnosis and therapy for OC in future.

Genetic and molecular biology of breast cancer among Iranian patients

Abstract

Background, Breast cancer (BC) is one of the leading causes of cancer related deaths in Iran. This high ratio of mortality had a rising trend during the recent years which is probably associated with late diagnosis.

Main body

Therefore it is critical to define a unique panel of genetic markers for the early detection among our population. In present review we summarized all of the reported significant genetic markers among Iranian BC patients for the first time, which are categorized based on their cellular functions.

Conclusions

This review paves the way of introducing a unique ethnic specific panel of diagnostic markers among Iranian BC patients. Indeed, this review can also clarify the genetic and molecular bases of BC progression among Iranians.

Effects of gossypol on apoptosis‑related gene expression in racially distinct triple‑negative breast cancer cells

Apoptosis is a gene‑directed mechanism that regulates cell proliferation and maintains homeostasis. Moreover, an aberrant apoptotic process can lead to several pathological conditions, such as tumorigenesis and cancer metastasis. In the present study, the apoptotic effect of the natural polyphenol compound gossypol GOSS) was investigated in triple‑negative breast cancer TNBC) cells. The effect of GOSS was evaluated in two cell lines representative of a Caucasian‑American and African‑American origin, MDA‑MB‑231 MM‑231) and MDA‑MB‑468 MM‑468), respectively. A similar response to both cytotoxicity and proliferation was observed in the two cell lines. However, MM‑468 cells were 2‑fold more sensitive to the apoptotic effect of the compound, which was accompanied by a longer delay in colony formation. Furthermore, GOSS was found to alter the mRNA expression of many apoptosis‑related genes. The compound significantly upregulated growth arrest and DNA damage‑inducible 45 alpha protein (GADD45A), tumor necrosis factor receptor superfamily 9 (TNFRSF9) and BCL2 interacting protein 3 BNIP3) in MM‑231 cells. Similarly, GADD45A and BNIP3 were upregulated in MM‑468 cells. A significant finding in this study is the profound 159‑fold increase in TNF gene expression that was observed in MM‑468 cells. Moreover, the apoptosis‑suppressor gene baculoviral IAP repeat containing 5 BIRC5) was significantly repressed (by more than 90%) in both cell lines, as well as death‑associated protein kinase 1 (DAPK1) in MM‑231 cells and tumor protein 73 (TP73) in MM‑468 cells. In conclusion, the data obtained in this study provide a molecular understanding of the GOSS‑induced apoptosis effect and suggest the importance of this polyphenol compound targeted towards TNBC treatment, particularly in African‑American women.

Integrated bioinformatics analysis to identify 15 hub genes in breast cancer

The aim of the present study was to identify the hub genes and provide insight into the tumorigenesis and development of breast cancer. To examine the hub genes in breast cancer, integrated bioinformatics analysis was performed. Gene expression profiles were obtained from the Gene Expression Omnibus (GEO) database and the differentially expressed genes (DEGs) were identified using the ‘limma’ package in R. Gene Ontology enrichment analysis and Kyoto Encyclopedia of Genes and Genomes pathway analysis was used to determine the functional annotations and potential pathways of the DEGs. Subsequently, a protein-protein interaction network analysis and weighted correlation network analysis (WGCNA) were conducted to identify hub genes. To confirm the reliability of the identified hub genes, RNA gene expression profiles were obtained from The Cancer Genome Atlas (TCGA)-breast cancer database, and WGCNA was used to screen for genes that were markedly correlated with breast cancer. By combining the results from the GEO and TCGA datasets, 15 hub genes were identified to be associated with breast cancer pathophysiology. Overall survival analysis was performed to examine the association between the expression of hub genes and the overall survival time of patients with breast cancer. Higher expression of all hub genes was associated with significantly shorter overall survival in patients with breast cancer compared with patients with lower levels of expression of the respective gene.

Polymorphisms in Survivin (BIRC5 Gene) Are Associated with Age of Onset in Breast Cancer Patients

Survivin, encoded by BIRC5 gene (baculoviral IAP repeat containing 5), belongs to the family of inhibitors of apoptosis proteins (IAPs). In mammalian cells it participates in the control of mitosis, apoptosis regulation, and cellular stress response. Its expression is increased in almost all types of cancers. The aim of this study was to investigate the role of BIRC5 polymorphisms in breast cancer (BC) and to connect survivin expression with various clinicopathological characteristics of BC patients. Blood and archival tumour tissue samples were collected from 26 BC patients from Croatia. Survivin expression was determined immunohistochemically. BIRC5 promoter, coding region, and 3'UTR were genotyped. DNA from 74 healthy women was used as control. BIRC5 polymorphisms and survivin expression were tested against age of onset, histological grade, tumour type and size, lymph node status, oestrogen, progesterone, Her2, and Ki67 status. Numbers of samples with weak, moderate, and strong survivin expression were 9 (33.3%), 11 (40.7%), and 7 (25.9%), respectively. Most patients had nuclear survivin staining (92.6%). High survivin expression was significantly associated with negative oestrogen receptor status (p=0.007) and positive Ki67 expression (p=0.032). Ki67 expression was also positively correlated with histological grade (p=0.0009). Fourteen polymorphisms were found in BC samples, located mostly in promoter and 3'UTR of BIRC5. There was no significant difference in the distribution of polymorphisms between BC and control samples. Among clinicopathological characteristics of BC patients, alleles of five BIRC5 polymorphisms were associated with younger age of onset: c.-644T>C (55.8 years [y] vs. 48.1 y; p=0.006), c.-241C>T (54.2 y vs. 45.0; p=0.029), c.9809T>C (55.8 y vs. 48.1 y; p=0.006), c.-1547C>T (58.3 y vs. 50.9 y; p=0.011), and c.9386T>C (50.8 y vs. 59.5 y; p=0.004). To assess the significance of BIRC5 polymorphisms and survivin expression as predictive and prognostic biomarkers for BC further research with a larger sample size is needed.

Mitochondrial DNA copy number instability in ERBB2-amplified breast cancer tumors

Increase in the copy number of ERBB2, a Tyrosine Kinase Receptor (TKR) leads to the overexpression of oncogene product and consequently uncontrolled cell proliferation which has been reported in different aggressive cancers with mitochondrial malfunctions. Although, amplification of ERBB2 has been reported in different studies; however, the association between changes in mitochondrial DNA content and the ERBB2 gene copy number is poorly understood. The relative mitochondrial DNA content of breast cancer tumor tissues of 70 patients who were referred to Imam Khomeini Hospital Complex was determined using quantitative Real-time PCR. Multiplex ligation-dependent probe amplification (MLPA) was conducted to evaluate the ERBB2 gene copy number variation and finally, two-sample Wilcoxon rank-sum (Mann-Whitney) test was used to investigate the possible association between mitochondrial DNA (mtDNA) content and the ERBB2 gene amplification. Seventeen out of 70 breast cancer tumor tissues were found with ERBB2 gene amplification. Comparison of the mitochondrial DNA content of the aforementioned samples with the rest of the cases showed a significant decrease in the mitochondrial DNA content of the ERBB2-amplified samples (P=0.01). Our data provided evidence that ERBB2 have the potential to have a regulatory role over mitochondrial activity by controlling the mtDNA content.

Integrated bioinformatics analysis reveals key candidate genes and pathways in breast cancer

Breast cancer (BC) is the leading malignancy in women worldwide, yet relatively little is known about the genes and signaling pathways involved in BC tumorigenesis and progression. The present study aimed to elucidate potential key candidate genes and pathways in BC. Five gene expression profile data sets (GSE22035, GSE3744, GSE5764, GSE21422 and GSE26910) were downloaded from the Gene Expression Omnibus (GEO) database, which included data from 113 tumorous and 38 adjacent non-tumorous tissue samples. Differentially expressed genes (DEGs) were identified using t-tests in the limma R package. These DEGs were subsequently investigated by pathway enrichment analysis and a protein-protein interaction (PPI) network was constructed. The most significant module from the PPI network was selected for pathway enrichment analysis. In total, 227 DEGs were identified, of which 82 were upregulated and 145 were downregulated. Pathway enrichment analysis results revealed that the upregulated DEGs were mainly enriched in ‘cell division’, the ‘proteinaceous extracellular matrix (ECM)’, ‘ECM structural constituents’ and ‘ECM-receptor interaction’, whereas downregulated genes were mainly enriched in ‘response to drugs’, ‘extracellular space’, ‘transcriptional activator activity’ and the ‘peroxisome proliferator-activated receptor signaling pathway’. The PPI network contained 174 nodes and 1,257 edges. DNA topoisomerase 2-a, baculoviral inhibitor of apoptosis repeat-containing protein 5, cyclin-dependent kinase 1, G2/mitotic-specific cyclin-B1 and kinetochore protein NDC80 homolog were identified as the top 5 hub genes. Furthermore, the genes in the most significant module were predominantly involved in ‘mitotic nuclear division’, ‘mid-body’, ‘protein binding’ and ‘cell cycle’. In conclusion, the DEGs, relative pathways and hub genes identified in the present study may aid in understanding of the molecular mechanisms underlying BC progression and provide potential molecular targets and biomarkers for BC.

Identification of novel prognosis-related genes associated with cancer using integrative network analysis

Prognosis identifies the seriousness and the chances of survival of a cancer patient. However, it remains a challenge to identify the key cancer genes in prognostic studies. In this study, we collected 2064 genes that were related to prognostic studies by using gene expression measurements curated from published literatures. Among them, 1820 genes were associated with copy number variations (CNVs). The further functional enrichment on 889 genes with frequent copy number gains (CNGs) revealed that these genes were significantly associated with cancer pathways including regulation of cell cycle, cell differentiation and mitogen-activated protein kinase (MAPK) cascade. We further conducted integrative analyses of CNV and their target genes expression using the data from matched tumour samples of The Cancer Genome Atlas (TCGA). Ultimately, 95 key prognosis-related genes were extracted, with concordant CNG events and increased up-regulation in at least 300 tumour samples. These genes, and the number of samples in which they were found, included: ACTL6A (399), ATP6V1C1 (425), EBAG9 (412), FADD (308), MTDH (377), and SENP5 (304). This study provides the first observation of CNV in prognosis-related genes across pan-cancer. The systematic concordance between CNG and up-regulation of gene expression in these novel prognosis-related genes may indicate their prognostic significance.

Isolation and characterization of a new naturally immortalized human breast carcinoma cell line, KAIMRC1

Background

Breast cancer is one of the most common cancer and a leading cause of death in women. Up to date the most commonly used breast cancer cell lines are originating from Caucasians or Afro-Americans but rarely cells are being derived from other ethnic groups. Here we describe for the first time the establishment of a naturally transformed breast cancer cell line, KAIMRC1 from an Arab woman of age 62 suffering from stage IIB breast cancer (T2N1M0). Moreover, we have characterized these cells for the biological and molecular markers, induction of MAPK pathways as well as its response to different commercially available drugs and compounds.

Methods

Breast cancer tissue sections were minced and cultured in media for several weeks. KAIMRC1 cells were successfully isolated from one of the primary breast tumor tissue cultures without any enzymatic digestion. To study the growth characteristics of the cells, wound healing assay, clonogenic assay, cell proliferation assays and live cell time-lapse microscopy was performed. Karyotyping, Immunophenotyping and molecular pathway specific compound treatment was also performed. A selective breast cancer gene expression panel was used to identify genes involved in the signal transduction dysregulation and malfunction of normal biological processes during breast carcinogenesis.

Results

These cells are ER/PR-positive and HER2-negative. The epithelial nature of these cells was confirmed by flow cytometry analysis using epithelial cell markers. They are cuboidal in shape and relatively smaller in size as compared to established cell lines, MCF-7, MDA MB-231 and the normal breast cell line, MCF-10A. In normal cell culture conditions these cells showed the capability of growing both in monolayer as well as in 3-D conformation. They showed a doubling time in vitro of approximately 24 h. They exhibit a modal karyotype of 58–63,X with abnormalities in a couple of chromosomes. KAIMRC1 cells were found to be more responsive to drug treatment in vitro in comparison to the established MDA MB-231 and MCF-7 cell lines.

Conclusions

In conclusion we have isolated and characterized a new naturally immortalized breast cell line, KAIMRC1 with a potential to play a key role in opening up novel avenues towards the understanding of breast carcinoma.

Electronic supplementary material

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

Network analysis revealed aurora kinase dysregulation in five gynecological types of cancer

Gene markers are crucial for cancer prognosis and treatment. Previous studies have placed greater emphasis on individual diagnostic genes, thereby ignoring systemic-level attributes across diseases. Female-specific cells namely, breast, endometrium, cervical, ovarian and vulvar cells are highly susceptible to cancer. To date, a limited number of molecular studies have been performed that evaluate common biological processes across gynecological types of cancer. Differentially expressed genes in breast, cervical, endometrial, vulvar and ovarian cancer were utilized to construct protein-protein interaction networks, and to identify a common module across the five cancer types. A single common module with 8 nodes and 26 edges was mined among the five cancer systems. In total, four hub genes were present across the five cancer gene sets. Genes in the common module were enriched for the common pathways and associated diseases. The aurora kinase pathway was revealed to be conserved across the five cancer types surveyed. The present study, therefore, revealed that the aurora kinase pathway has a crucial function in the pathogenesis of the five aforementioned gynecological types of cancer through cross-tumor conservation.