Gene expression associated with lymphovascular invasion and genomic risk in early-stage breast cancer.

559 Background: Lymphovascular invasion (LVI), the passage of carcinoma cells through lymphatic and blood vessels, is an important early step in metastasis; however, LVI is excluded from most breast cancer (BC) clinical risk assessments. Previous studies assessed the prognostic value of LVI to estimate clinical outcomes. To gain understanding of the molecular basis of LVI, we evaluated differentially expressed genes (DEGs) between tumors with LVI versus those without LVI, stratified by the 70-gene signature (MammaPrint/MP) and 80-gene molecular subtyping signature (BluePrint/BP). Methods: The prospective, observational FLEX Study (NCT03053193) includes stage I-III BC patients who receive MP/BP testing and consent to full transcriptome and clinical data collection. Patients with LVI (n=581) and without LVI (n=600, randomly selected), enrolled from 2017 to present, were included. LVI was assessed by local pathology laboratories. Differential gene expression analysis of 44k Agilent microarray data was performed with R limma package. DEGs were compared within all samples, BP Luminal subtype, MP risk groups (Low Risk [LR]/Luminal A and High Risk [HR]/Luminal B), and by lymph node (LN) status. DEGs with FDR<0.05 were considered significant. Results: Of tumors with LVI (LVI+), 66% were MP HR; notably, 51% of tumors without LVI (LVI-) were MP HR. LVI was associated with larger T stage, LN involvement, high grade, negative ER status by IHC, and younger patient age (LVI+ vs. LVI-, p<0.05 for all comparisons). Patient ethnicity, obesity, and tumor type did not differ by LVI status; however, prevalence of type 2 diabetes trended higher in patients with LVI+ HR tumors (21%), compared with LVI- HR (15%, p=0.09) and LVI+ LR (11%, p=0.004). There were significant transcriptomic differences between LVI+ and LVI, with most DEGs evident in the Luminal B subset. DEGs in LVI+, LN-negative (LN-) tumors overlapped substantially with the overall Luminal group analysis. Functional enrichment analysis showed dysregulation of cell cycle, extracellular matrix (ECM) organization, cell adhesion, and cytokine receptor pathways. Gene sets related to insulin growth factor pathways were also enriched in LVI+ tumors. Conclusions: DEGs associated with LVI were primarily found in MP HR Luminal, LN-negative tumors; enrichment analysis suggested dysregulation of ECM organization and cell adhesion pathways, consistent with previous reports. DEGs were not associated with LVI presence in LN+ tumors, suggesting that LVI assessment may be less relevant in LN+ breast cancer. Future studies will assess clinical outcomes, as well as LVI-associated gene expression in BP Basal- and HER2-type tumors. However, the current analysis indicates few DEGs in LVI+ MP LR tumors; thus, the potential prognostic information gained from LVI-associated gene expression is likely already captured by the MP and BP signatures. Clinical trial information: NCT03053193.

Download Full-text

Co-expression of key gene modules and pathways of human breast cancer cell lines

Bioscience Reports ◽

10.1042/bsr20181925 ◽

2019 ◽

Vol 39 (7) ◽

Cited By ~ 2

Author(s):

Yadong Wu ◽

Feng liu ◽

Siyang Luo ◽

Xinhai Yin ◽

Dengqi He ◽

...

Keyword(s):

Breast Cancer ◽

Gene Expression ◽

Network Analysis ◽

Cell Lines ◽

Enrichment Analysis ◽

Functional Enrichment ◽

Correlation Network ◽

Breast Cancer Cell Lines ◽

Hub Genes ◽

Weighted Correlation

Abstract Breast cancer (BC) is the most common leading cause of cancer-related death in women worldwide. Gene expression profiling analysis for human BCs has been studied previously. However, co-expression analysis for BC cell lines is still devoid to date. The aim of the study was to identify key pathways and hub genes that may serve as a biomarker for BC and uncover potential molecular mechanism using weighted correlation network analysis. We analyzed microarray data of BC cell lines (GSE 48213) listed in the Gene Expression Omnibus database. Gene co-expression networks were used to construct and explore the biological function in hub modules using the weighted correlation network analysis algorithm method. Meanwhile, Gene ontology and KEGG pathway analysis were performed using Cytoscape plug-in ClueGo. The network of the key module was also constructed using Cytoscape. A total of 5000 genes were selected, 28 modules of co-expressed genes were identified from the gene co–expression network, one of which was found to be significantly associated with a subtype of BC lines. Functional enrichment analysis revealed that the brown module was mainly involved in the pathway of the autophagy, spliceosome, and mitophagy, the black module was mainly enriched in the pathway of colorectal cancer and pancreatic cancer, and genes in midnightblue module played critical roles in ribosome and regulation of lipolysis in adipocytes pathway. Three hub genes CBR3, SF3B6, and RHPN1 may play an important role in the development and malignancy of the disease. The findings of the present study could improve our understanding of the molecular pathogenesis of breast cancer.

Download Full-text

Identification of Potential Target Genes in HER-2 Positive Breast Cancer by Bioinformatics Analysis.

10.21203/rs.3.rs-34161/v1 ◽

2020 ◽

Author(s):

Song Wang ◽

Yi Quan ◽

Hongying Lyu ◽

Jian Deng

Keyword(s):

Breast Cancer ◽

Gene Expression ◽

Signaling Pathway ◽

Target Genes ◽

Enrichment Analysis ◽

Functional Enrichment Analysis ◽

Functional Enrichment ◽

Ppi Network ◽

Her 2 ◽

Positive Breast Cancer

Abstract Background: HER-2 positive breast cancer has a high risk of for relapse, metastasis and drug resistance, and is correlated with a poor prognosis. Thus, the study objective was to reveal target genes and key pathways in HER-2 subtype breast cancer. Methods: The gene expression dataset (GSE29431) was downloaded from the Gene Expression Omnibus database(GEO), and the differentially expressed genes (DEGs) were determined using LIMMA package in R software. Subsequently, Functional enrichment analysis were performed in ClusterProfiler package of R platform. The Search Tool for the Retrieval of Interacting Genes (STRING) database was used to construct a Protein-Protein Interaction (PPI) network of DEGs. Module analysis and target genes were identified by Cytoscape software. Further more, The influence of target genes on overall survival (OS) was assessed using the Kaplan-Meier plotter database.Results: The differential expression analysis revealed 96 genes were up-regulated while 407 genes were down-regulated in HER-2 positive breast cancer tissue compared to normal breast tissue. Functional enrichment analysis showed that the DEGs were mainly involved in regulation of lipid metabolic process, PPAR signaling pathway and PI3K-Akt signaling pathway. PPI network construction revealed a total of 199 nodes and 560 edges, and 12 target genes were identified by the highest value of degree. In addition, target genes were associated with worse overall prognosis, including NUSAP1, PTTG1, CEP55, TOP2A, CCNB1, CENPF, MELK, AURKA, UBE2C, BUB1B, KIF20A and RRM2.Conclusion: The present study identified 12 target genes associated with the development of HER-2 subtype breast cancer, which may help to provide new biomarkers and therapeutic targets.

Download Full-text

Inflammatory Breast Cancer: Clinical Implications of Genomic Alterations and Mutational Profiling

Cancers ◽

10.3390/cancers12102816 ◽

2020 ◽

Vol 12 (10) ◽

pp. 2816

Author(s):

Flávia L. C. Faldoni ◽

Rolando A. R. Villacis ◽

Luisa M. Canto ◽

Carlos E. Fonseca-Alves ◽

Sarah S. Cury ◽

...

Keyword(s):

Breast Cancer ◽

Inflammatory Breast Cancer ◽

Enrichment Analysis ◽

Microarray Platform ◽

Functional Enrichment Analysis ◽

Functional Enrichment ◽

Tumor Type ◽

Similar Frequency ◽

Targeted Next Generation Sequencing ◽

Mutational Profiling

Inflammatory breast cancer (IBC) is a rare and aggressive type of breast cancer whose molecular basis is poorly understood. We performed a comprehensive molecular analysis of 24 IBC biopsies naïve of treatment, using a high-resolution microarray platform and targeted next-generation sequencing (105 cancer-related genes). The genes more frequently affected by gains were MYC (75%) and MDM4 (71%), while frequent losses encompassed TP53 (71%) and RB1 (58%). Increased MYC and MDM4 protein expression levels were detected in 18 cases. These genes have been related to IBC aggressiveness, and MDM4 is a potential therapeutic target in IBC. Functional enrichment analysis revealed genes associated with inflammatory regulation and immune response. High homologous recombination (HR) deficiency scores were detected in triple-negative and metastatic IBC cases. A high telomeric allelic imbalance score was found in patients having worse overall survival (OS). The mutational profiling was compared with non-IBC (TCGA, n = 250) and IBC (n = 118) from four datasets, validating our findings. Higher frequency of TP53 and BRCA2 variants were detected compared to non-IBC, while PIKC3A showed similar frequency. Variants in mismatch repair and HR genes were associated with worse OS. Our study provided a framework for improved diagnosis and therapeutic alternatives for this aggressive tumor type.

Download Full-text

Identification of potential biomarkers for abdominal pain in IBS patients by bioinformatics approach

BMC Gastroenterology ◽

10.1186/s12876-021-01626-7 ◽

2021 ◽

Vol 21 (1) ◽

Author(s):

Zhongyuan Lin ◽

Yimin Wang ◽

Shiqing Lin ◽

Decheng Liu ◽

Guohui Mo ◽

...

Keyword(s):

Gene Expression ◽

Irritable Bowel Syndrome ◽

Abdominal Pain ◽

Gastrointestinal Disease ◽

Rectal Mucosa ◽

Interaction Network ◽

Enrichment Analysis ◽

Functional Enrichment ◽

Irritable Bowel ◽

Potential Biomarkers

Abstract Background Irritable bowel syndrome (IBS) is the most common functional gastrointestinal disease characterized by chronic abdominal discomfort and pain. The mechanisms of abdominal pain, as a relevant symptom, in IBS are still unclear. We aimed to explore the key genes and neurobiological changes specially involved in abdominal pain in IBS. Methods Gene expression data (GSE36701) was downloaded from Gene Expression Omnibus database. Fifty-three rectal mucosa samples from 27 irritable bowel syndrome with diarrhea (IBS-D) patients and 40 samples from 21 healthy volunteers as controls were included. Differentially expressed genes (DEGs) between two groups were identified using the GEO2R online tool. Functional enrichment analysis of DEGs was performed on the DAVID database. Then a protein–protein interaction network was constructed and visualized using STRING database and Cytoscape. Results The microarray analysis demonstrated a subset of genes (CCKBR, CCL13, ACPP, BDKRB2, GRPR, SLC1A2, NPFF, P2RX4, TRPA1, CCKBR, TLX2, MRGPRX3, PAX2, CXCR1) specially involved in pain transmission. Among these genes, we identified GRPR, NPFF and TRPA1 genes as potential biomarkers for irritating abdominal pain of IBS patients. Conclusions Overexpression of certain pain-related genes (GRPR, NPFF and TRPA1) may contribute to chronic visceral hypersensitivity, therefore be partly responsible for recurrent abdominal pain or discomfort in IBS patients. Several synapses modification and biological process of psychological distress may be risk factors of IBS.

Download Full-text

Integrative transcriptome data mining for identification of core lncRNAs in breast cancer

PeerJ ◽

10.7717/peerj.7821 ◽

2019 ◽

Vol 7 ◽

pp. e7821 ◽

Cited By ~ 7

Author(s):

Xiaoming Zhang ◽

Jing Zhuang ◽

Lijuan Liu ◽

Zhengguo He ◽

Cun Liu ◽

...

Keyword(s):

Breast Cancer ◽

Gene Expression ◽

Prognostic Value ◽

Early Stage ◽

Expression Profiles ◽

Diagnostic Value ◽

The Cancer Genome Atlas ◽

Functional Enrichment ◽

Differentially Expressed ◽

Data Set

Background Cumulative evidence suggests that long non-coding RNAs (lncRNAs) play an important role in tumorigenesis. This study aims to identify lncRNAs that can serve as new biomarkers for breast cancer diagnosis or screening. Methods First, the linear fitting method was used to identify differentially expressed genes from the breast cancer RNA expression profiles in The Cancer Genome Atlas (TCGA). Next, the diagnostic value of all differentially expressed lncRNAs was evaluated using a receiver operating characteristic (ROC) curve. Then, the top ten lncRNAs with the highest diagnostic value were selected as core genes for clinical characteristics and prognosis analysis. Furthermore, core lncRNA-mRNA co-expression networks based on weighted gene co-expression network analysis (WGCNA) were constructed, and functional enrichment analysis was performed using the Database for Annotation, Visualization and Integrated Discovery (DAVID). The differential expression level and diagnostic value of core lncRNAs were further evaluated by using independent data set from Gene Expression Omnibus (GEO). Finally, the expression status and prognostic value of core lncRNAs in various tumors were analyzed based on Gene Expression Profiling Interactive Analysis (GEPIA). Results Seven core lncRNAs (LINC00478, PGM5-AS1, AL035610.1, MIR143HG, RP11-175K6.1, AC005550.4, and MIR497HG) have good single-factor diagnostic value for breast cancer. AC093850.2 has a prognostic value for breast cancer. AC005550.4 and MIR497HG can better distinguish breast cancer patients in early-stage from the advanced-stage. Low expression of MAGI2-AS3, LINC00478, AL035610.1, MIR143HG, and MIR145 may be associated with lymph node metastasis in breast cancer. Conclusion Our study provides candidate biomarkers for the diagnosis and prognosis of breast cancer, as well as a bioinformatics basis for the further elucidation of the molecular pathological mechanism of breast cancer.

Download Full-text

Identification of a prognostic long non-coding RNA signature in breast cancer

10.21203/rs.3.rs-27286/v1 ◽

2020 ◽

Author(s):

Gaochen Lan ◽

Xiaoling Yu ◽

Yanna Zhao ◽

Jinjian Lan ◽

Wan Li ◽

...

Keyword(s):

Breast Cancer ◽

Cox Regression ◽

Prognostic Biomarker ◽

Expression Profiles ◽

Enrichment Analysis ◽

Functional Enrichment Analysis ◽

The Cancer Genome Atlas ◽

Functional Enrichment ◽

Differentially Expressed ◽

Cox Regression Analysis

Abstract Background: Breast cancer is the most common malignant disease among women. At present, more and more attention has been paid to long non-coding RNAs (lncRNAs) in the field of breast cancer research. We aimed to investigate the expression profiles of lncRNAs and construct a prognostic lncRNA for predicting the overall survival (OS) of breast cancer.Methods: The expression profiles of lncRNAs and clinical data with breast cancer were obtained from The Cancer Genome Atlas (TCGA). Differentially expressed lncRNAs were screened out by R package (limma). The survival probability was estimated by the Kaplan‑Meier Test. The Cox Regression Model was performed for univariate and multivariate analysis. The risk score (RS) was established on the basis of the lncRNAs’ expression level (exp) multiplied regression coefficient (β) from the multivariate cox regression analysis with the following formula: RS=exp a1 * β a1 + exp a2 * β a2 +……+ exp an * β an. Functional enrichment analysis was performed by Metascape.Results: A total of 3404 differentially expressed lncRNAs were identified. Among them, CYTOR, MIR4458HG and MAPT-AS1 were significantly associated with the survival of breast cancer. Finally, The RS could predict OS of breast cancer (RS=exp CYTOR * β CYTOR + exp MIR4458HG * β MIR4458HG + exp MAPT-AS1 * β MAPT-AS1). Moreover, it was confirmed that the three-lncRNA signature could be an independent prognostic biomarker for breast cancer (HR=3.040, P=0.000).Conclusions: This study established a three-lncRNA signature, which might be a novel prognostic biomarker for breast cancer.

Download Full-text

Diagnostic model of combined ceRNA and DNA methylation related genes in esophageal carcinoma

PeerJ ◽

10.7717/peerj.8831 ◽

2020 ◽

Vol 8 ◽

pp. e8831 ◽

Cited By ~ 1

Author(s):

Xiaojiao Guan ◽

Yao Yao ◽

Guangyao Bao ◽

Yue Wang ◽

Aimeng Zhang ◽

...

Keyword(s):

Gene Expression ◽

Esophageal Cancer ◽

Esophageal Carcinoma ◽

Differentially Expressed Genes ◽

Enrichment Analysis ◽

Functional Enrichment ◽

Differentially Expressed ◽

Diagnostic Model ◽

Link Type ◽

Key Genes

Esophageal cancer is a common malignant tumor in the world, and the aim of this study was to screen key genes related to the development of esophageal cancer using a variety of bioinformatics analysis tools and analyze their biological functions. The data of esophageal squamous cell carcinoma from the Gene Expression Omnibus (GEO) were selected as the research object, processed and analyzed to screen differentially expressed microRNAs (miRNAs) and differential methylation genes. The competing endogenous RNAs (ceRNAs) interaction network of differentially expressed genes was constructed by bioinformatics tools DAVID, String, and Cytoscape. Biofunctional enrichment analysis was performed using Gene Ontology (GO) and the Kyoto Encyclopedia of Genes and Genomes (KEGG). The expression of the screened genes and the survival of the patients were verified. By analyzing GSE59973 and GSE114110, we found three down-regulated and nine up-regulated miRNAs. The gene expression matrix of GSE120356 was calculated by Pearson correlation coefficient, and the 11696 pairs of ceRNA relation were determined. In the ceRNA network, 643 lncRNAs and 147 mRNAs showed methylation difference. Functional enrichment analysis showed that these differentially expressed genes were mainly concentrated in the FoxO signaling pathway and were involved in the corresponding cascade of calcineurin. By analyzing the clinical data in The Cancer Genome Atlas (TCGA) database, it was found that four lncRNAs had an important impact on the survival and prognosis of esophageal carcinoma patients. QRT-PCR was also conducted to identify the expression of the key lncRNAs (RNF217-AS1, HCP5, ZFPM2-AS1 and HCG22) in ESCC samples. The selected key genes can provide theoretical guidance for further research on the molecular mechanism of esophageal carcinoma and the screening of molecular markers.

Download Full-text

Prognostic Significance of SQSTM1 in Breast Cancer: A Comprehensive Analysis

10.21203/rs.3.rs-64416/v1 ◽

2020 ◽

Author(s):

Yang Liu ◽

Qian Du ◽

Dan Sun ◽

Ruiying Han ◽

Mengmeng Teng ◽

...

Keyword(s):

Breast Cancer ◽

Cancer Patients ◽

Somatic Mutation ◽

Enrichment Analysis ◽

Gene Set Enrichment Analysis ◽

Functional Enrichment ◽

Regulation Mechanism ◽

Breast Cancer Patients ◽

Gene Set ◽

Multiple Datasets

Abstract Background: SQSTM1 (Sequestosome 1, p62) is degraded by activated autophagy and involved in the progression of in various types of cancers. However, the prognostic role and underlying regulation mechanism of SQSTM1 in the progression and development of breast cancer remain unclear.Methods: In this study, 1336 samples with available mRNA data from Molecular Taxonomy of Breast Cancer International Consortium (METABRIC) database and 27 formalin fixation and paraffin embedding (FFPE) tissue samples from the First Affiliated Hospital of Xi’an Jiaotong University were collected to evaluate SQSTM1 expression in mRNA and protein levels. Kaplan–Meier and Cox regression were used for revealing prognostic value in three independent breast cancer independent datasets. Tumor Immune Estimation Resource (TIMER) database and Gene Set Variation Analysis (GSVA) was used to explore the relationship of SQSTM1 mRNA expression and immune infiltration level in breast cancer. Dysregulation mechanisms of SQSTM1 were also explored including copy number variation (CNV), somatic mutation, epigenetic alterations and other transcription and post-transcription level using multiple datasets. Finally, Gene Set Enrichment Analysis (GSEA) was constructed to elucidate functional regulating performance of SQSTM1 in breast cancer.Results: The results showed that mRNA and protein level of SQSTM1 were significantly elevated in breast cancer and receiver operating characteristic (ROC) curve showed that p62 may act as diagnostic biomarker. Lower expression of SQSTM1 predicted better outcome through multiple datasets. It was also found that SQSTM1 correlated with immune infiltrates in breast cancer. Moreover, CNV and methylation of SQSTM1 DNA was correlated with SQSTM1 dysregulation and act as prognostic factors for breast cancer patients. Yet, somatic mutation status of SQSTM1 didn’t show any prognostic relevance. We also identified diverse transcription factors that directly bound to SQSTM1 DNA and the miRNAs which may regulate SQSTM1 mRNA. Finally, functional enrichment analysis revealed that SQSTM1 is related to cell signal transduction, oxidative stress and autophagy in breast cancer.Conclusion: Our findings revealed that SQSTM1 plays a key role in the progression of breast cancer and might be a promising biomarker for the diagnosis and personalized treatment of breast cancer patients.

Download Full-text

Nomogram Models Based on Gene Expression in Prediction of Breast Cancer Bone Metastasis

10.21203/rs.3.rs-1109324/v1 ◽

2021 ◽

Author(s):

Teng-di Fan ◽

Di-kai Bei ◽

Song-wei Li

Keyword(s):

Breast Cancer ◽

Gene Expression ◽

Bone Metastasis ◽

Gene Signature ◽

Functional Enrichment ◽

Training Set ◽

Hub Genes ◽

Kaplan Meier ◽

Patient Prognosis ◽

Kaplan Meier Plotter

Abstract Objective: To design a weighted co-expression network and build gene expression signature-based nomogram (GESBN) models for predicting the likelihood of bone metastasis in breast cancer (BC) patients. Methods: Dataset GSE124647 was used as a training set, and GSE14020 was taken as a validation set. In the training cohort, limma package in R was adopted to obtain differentially expressed genes (DEGs) between BC non-bone metastasis and bone metastasis patients, which were used for functional enrichment analysis. After weighted co-expression network analysis (WGCNA), univariate Cox regression and Kaplan-Meier plotter analyses were performed to screen potential prognosis-related genes. Then, GESBN models were constructed and evaluated. Further, the expression levels of genes in the models were explored in the training set, which was validated in GSE14020. Finally, the prognostic value of hub genes in BC was explored. Results: A total of 1858 DEGs were obtained. WGCNA result showed that the blue module was most significantly related to bone metastasis and prognosis. After survival analyses, GAJ1, SLC24A3, ITGBL1, and SLC44A1 were subjected to construct a GESBN model for overall survival. While GJA1, IGFBP6, MDFI, ITGFBI, ANXA2, and SLC24A3 were subjected to build a GESBN model for progression-free survival. Kaplan-Meier plotter and receiver operating characteristic analyses presented the reliable prediction ability of the models. Besides, GJA1, IGFBP6, ITGBL1, SLC44A1, and TGFBI expressions were significantly different between the two groups in GSE124647 and GSE14020. The hub genes had a significant impact on patient prognosis. Conclusion: Both the four-gene signature and six-gene signature could accurately predict patient prognosis, which may provide novel treatment insights for BC bone metastasis.

Download Full-text