Weighted gene co-expression network analysis of the association between upregulated AMD1, EN1 and VGLL1 and the progression and poor prognosis of breast cancer

Secondary Analysis of Human Bulk RNA-Seq Dataset Suggests Potential Mechanisms for Letrozole Resistance in Estrogen-Positive (ER+) Breast Cancer

Estrogen receptor-positive (ER+) breast cancer is common among postmenopausal women and is frequently treated with Letrozole, which inhibits aromatase from synthesizing estrogen from androgens. Decreased estrogen slows the growth of tumors and can be an effective treatment. The increase in Letrozole resistance poses a unique problem for patients. To better understand the underlying molecular mechanism(s) of Letrozole resistance, we reanalyzed transcriptomic data by comparing individuals who responded to Letrozole therapy (responders) to those who were resistant to treatment (non-responders). We identified SOX11 and S100A9 as two significant differentially expressed genes (DEGs) between these patient cohorts, with "PLK1 signaling events" being the most significant signaling pathway. We also identified PRDX4 and E2F8 gene products as being the top mechanistic transcriptional markers for ER+ treatment resistance. Many of the significant DEGs that we identified play a known role in ER+ breast cancer or other types of cancer, which partially validate our results. Several of the gene products we identified are novel in the context of ER+ breast cancer. Many of the genes that we identified warrant further research to elucidate the more specific molecular mechanisms of Letrozole resistance in this patient population and could potentially be used as prognostic markers with further wet lab validation. We anticipate that these findings could contribute to improved detection and therapeutic outcomes in aromatase-resistant ER+ breast cancer patients.

VGLL1 stabilization of cytoplasmic TAZ promotes EGFR expression and maintains tumor initiating cells in breast cancer independent of TEAD

Vestigial-like family member 1 (VGLL1) is one of the X-linked genes whose expression is elevated in basal-like breast cancer (BLBC) because of X-chromosome isodisomy. As an approach towards understanding its function, we performed correlation study using transcript data of breast cancer patients from cBioPortal for Cancer Genomics. Our analysis identified EGFR as the most correlated transcript with VGLL1. We demonstrate that VGLL1 promotes EGFR expression and increases the frequency of breast tumor initiating cells (CD44high/+CD24low/-). These findings are crucial because an elevated EGFR expression and high frequency of CD44high/+CD24low/- cells are defining features of BLBC, and we provide a new mechanistic insight into how their expressions are controlled. Importantly, VGLL1 regulation of EGFR and CD44high/+CD24low/- population is mediated by the hippo-transducer TAZ which exerts its oncogenic roles by binding and activating TEAD transcription factors. A crucial finding is that TEAD-binding domain of TAZ is dispensable for its regulation of EGFR and CD44high/+CD24low/- cells. Instead, VGLL1 stabilization of cytoplasmic TAZ is essential for these functions. Also, we show that VGLL1/TAZ restricts the surface expression of CD24 which contributes to the increased number of CD44high/+CD24low/- cells. In addition, we observed that VGLL1 represses AXL expression and suppresses claudin-low phenotype, and that is caused by the VGLL1 mediated nuclear expulsion of TAZ. Therefore, EGFR and AXL seem to represent two different breast tumor subtypes, and their differential expressions is controlled by VGLL1.

Development of a polyamine gene expression score for predicting prognosis and treatment response in clear cell renal cell carcinoma

Backgrounds

Polyamine metabolism (PM) is closely related to the tumor microenvironment (TME) and is involved in antitumor immunity. Clear cell renal cell carcinoma (ccRCC) not only has high immunogenicity but also has significant metabolic changes. However, the role of PM in the immune microenvironment of ccRCC remains unclear. This study aimed to reveal the prognostic value of PM-related genes (PMRGs) expression in ccRCC and their correlation with the TME.

Methods

The expression levels PMRGs in different cells were characterized with single-cell sequencing analysis. The PMRG expression pattern of 777 ccRCC patients was evaluated based on PMRGs. Unsupervised clustering analysis was used in identifying PMRG expression subtypes, and Lasso regression analysis was used in developing polyamine gene expression score (PGES), which was validated in external and internal data sets. The predictive value of PGES for immunotherapy was validated in the IMvigor210 cohort. Multiple algorithms were used in analyzing the correlation between PGES and immune cells. The sensitivity of PGES to chemotherapeutic drugs was analyzed with the "pRRophetic" package. We validated the genes that develop PGES in tissue samples. Finally, weighted gene co-expression network analysis was used in identifying the key PMRGs closely related to ccRCC, and cell function experiments were carried out.

Results

PMRGs were abundantly expressed on tumor cells, and PMRG expression was active in CD8⁺ T cells and fibroblasts. We identified three PMRG expression subtypes. Cancer and immune related pathways were active in PMRG expression cluster A, which had better prognosis. PGES exhibited excellent predictive value. The high-PGES group was characterized by high immune cell infiltration, high expression of T cell depletion markers, high tumor mutation burden and tumor immune dysfunction and exclusion, was insensitive to immunotherapy but sensitive to sunitinib, temsirolimus, and rapamycin, and had poor prognosis. Spermidine synthetase (SRM) has been identified as a key gene and is highly expressed in ccRCC at RNA and protein levels. SRM knockdown can inhibit ccRCC cell proliferation, migration, and invasion.

Conclusions

We revealed the biological characteristics of PMRG expression subtypes and developed PGES to accurately predict the prognosis of patients and response to immunotherapy.

Polygenic and Network-Based Studies in Risk Identification and Demystification of cancer

INTRODUCTION

Diseases were initially thought to be the consequence of a single gene mutation. Advances in DNA sequencing tools and our understanding of gene behavior have revealed that complex diseases, such as cancer, are the product of genes cooperating with each other and with their environment in orchestrated communication networks. Seeing that the function of individual genes is still used to analyze cancer, the shift to using functionally interacting groups of genes as a new unit of study holds promise for demystifying cancer.

AREAS COVERED

The literature search focused on three types of cancer, namely breast, lung, and prostate, but arguments from other cancers were also included. The aim was to prove that multigene analyses can accurately predict and prognosticate cancer risk, subtype cancer for more personalized and effective treatments, and discover anti-cancer therapies. Computational intelligence is being harnessed to analyze this type of data and is proving indispensable to scientific progress.

EXPERT OPINION

In the future, comprehensive profiling of all kinds of patient data (e.g., serum molecules, environmental exposures) can be used to build universal networks that should help us elucidate the molecular mechanisms underlying diseases and provide appropriate preventive measures, ensuring lifelong health and longevity.

Identification and Validation of Key Genes of Differential Correlations in Gastric Cancer

Background: Gastric cancer (GC) is aggressive cancer with a poor prognosis. Previously bulk transcriptome analysis was utilized to identify key genes correlated with the development, progression and prognosis of GC. However, due to the complexity of the genetic mutations, there is still an urgent need to recognize core genes in the regulatory network of GC.

Methods: Gene expression profiles (GSE66229) were retrieved from the GEO database. Weighted correlation network analysis (WGCNA) was employed to identify gene modules mostly correlated with GC carcinogenesis. R package ‘DiffCorr’ was applied to identify differentially correlated gene pairs in tumor and normal tissues. Cytoscape was adopted to construct and visualize the gene regulatory network.

Results: A total of 15 modules were detected in WGCNA analysis, among which three modules were significantly correlated with GC. Then genes in these modules were analyzed separately by “DiffCorr”. Multiple differentially correlated gene pairs were recognized and the network was visualized by the software Cytoscape. Moreover, GEMIN5 and PFDN2, which were rarely discussed in GC, were identified as key genes in the regulatory network and the differential expression was validated by real-time qPCR, WB and IHC in cell lines and GC patient tissues.

Conclusions: Our research has shed light on the carcinogenesis mechanism by revealing differentially correlated gene pairs during transition from normal to tumor. We believe the application of this network-based algorithm holds great potential in inferring relationships and detecting candidate biomarkers.

Estimation of Gene Regulatory Networks from Cancer Transcriptomics Data

Cancer is a genetic disease in which multiple genes are perturbed. Thus, information about the regulatory relationships between genes is necessary for the identification of biomarkers and therapeutic targets. In this review, methods for inference of gene regulatory networks (GRNs) from transcriptomics data that are used in cancer research are introduced. The methods are classified into three categories according to the analysis model. The first category includes methods that use pair-wise measures between genes, including correlation coefficient and mutual information. The second category includes methods that determine the genetic regulatory relationship using multivariate measures, which consider the expression profiles of all genes concurrently. The third category includes methods using supervised and integrative approaches. The supervised approach estimates the regulatory relationship using a supervised learning method that constructs a regression or classification model for predicting whether there is a regulatory relationship between genes with input data of gene expression profiles and class labels of prior biological knowledge. The integrative method is an expansion of the supervised method and uses more data and biological knowledge for predicting the regulatory relationship. Furthermore, simulation and experimental validation of the estimated GRNs are also discussed in this review. This review identified that most GRN inference methods are not specific for cancer transcriptome data, and such methods are required for better understanding of cancer pathophysiology. In addition, more systematic methods for validation of the estimated GRNs need to be developed in the context of cancer biology.