Multi-Omic Regulation of the PAM50 Gene Signature in Breast Cancer Molecular Subtypes

Comprehensive Transcriptome Analysis Expands lncRNA Functional Profiles in Breast Cancer

Breast cancer is a heterogeneous disease that arises as a multi-stage process involving multiple cell types. Patients diagnosed with the same clinical stage and pathological classification may have different prognoses and therapeutic responses due to alterations in molecular genetics. As an essential marker for the molecular subtyping of breast cancer, long non-coding RNAs (lncRNAs) play a crucial role in gene expression regulation, cell differentiation, and the maintenance of genomic stability. Here, we developed a modular framework for lncRNA identification and applied it to a breast cancer cohort to identify novel lncRNAs not previously annotated. To investigate the potential biological function, regulatory mechanisms, and clinical relevance of the novel lncRNAs, we elucidated the genomic and chromatin features of these lncRNAs, along with the associated protein-coding genes and putative enhancers involved in the breast cancer regulatory networks. Furthermore, we uncovered that the expression patterns of novel and annotated lncRNAs identified in breast cancer were related to the hormone response in the PAM50 subtyping criterion, as well as the immune response and progression states of breast cancer across different immune cells and immune checkpoint genes. Collectively, the comprehensive identification and functional analysis of lncRNAs revealed that these lncRNAs play an essential role in breast cancer by altering gene expression and participating in the regulatory networks, contributing to a better insight into breast cancer heterogeneity and potential avenues for therapeutic intervention.

Few-shot genes selection: subset of PAM50 genes for breast cancer subtypes classification

BACKGROUND

In recent years, researchers have made significant strides in understanding the heterogeneity of breast cancer and its various subtypes. However, the wealth of genomic and proteomic data available today necessitates efficient frameworks, instruments, and computational tools for meaningful analysis. Despite its success as a prognostic tool, the PAM50 gene signature's reliance on many genes presents challenges in terms of cost and complexity. Consequently, there is a need for more efficient methods to classify breast cancer subtypes using a reduced gene set accurately.

RESULTS

This study explores the potential of achieving precise breast cancer subtype categorization using a reduced gene set derived from the PAM50 gene signature. By employing a "Few-Shot Genes Selection" method, we randomly select smaller subsets from PAM50 and evaluate their performance using metrics and a linear model, specifically the Support Vector Machine (SVM) classifier. In addition, we aim to assess whether a more compact gene set can maintain performance while simplifying the classification process. Our findings demonstrate that certain reduced gene subsets can perform comparable or superior to the full PAM50 gene signature.

CONCLUSIONS

The identified gene subsets, with 36 genes, have the potential to contribute to the development of more cost-effective and streamlined diagnostic tools in breast cancer research and clinical settings.

Drug repurposing for Basal breast cancer subpopulations using modular network signatures

Breast cancer is characterized as being a heterogeneous pathology with a broad phenotype variability. Breast cancer subtypes have been developed in order to capture some of this heterogeneity. Each of these breast cancer subtypes, in turns retains varied characteristic features impacting diagnostic, prognostic and therapeutics. Basal breast tumors, in particular have been challenging in these regards. Basal breast cancer is often more aggressive, of rapid evolution and no tailor-made targeted therapies are available yet to treat it. Arguably, epigenetic variability is behind some of these intricacies. It is possible to further classify basal breast tumor in groups based on their non-coding transcriptome and methylome profiles. It is expected that these groups will have differences in survival as well as in sensitivity to certain classes of drugs. With this in mind, we implemented a computational learning approach to infer different subpopulations of basal breast cancer (from TCGA multi-omic data) based on their epigenetic signatures. Such epigenomic signatures were associated with different survival profiles; we then identified their associated gene co-expression network structure, extracted a signature based on modules within these networks, and use these signatures to find and prioritize drugs (in the LINCS dataset) that may be used to target these types of cancer. In this way we are introducing the analytical workflow for an epigenomic signature-based drug repurposing structure.

Transcription factor E2F8 is a therapeutic target in the basal-like subtype of breast cancer

Introduction

Tumorigenesis in breast cancers usually accompanied by the dysregulation of transcription factors (TFs). Abnormal amplification of TFs leads aberrant expression of its downstream target genes. However, breast cancers are heterogeneous disease with different subtypes that have distinguished clinical behaviours, and the identification of prognostic TFs may enable to provide diagnosis and treatment of breast cancer based on subtypes, especially in Basal-like breast cancer.

Methods

The RNA-sequencing was performed to screen differential TFs in breast cancer subtypes. The GEPIA dataset analysis was used to analyze the genes expression in invasive breast carcinoma. The expression of MYBL2, HOXC13, and E2F8 was verified by qRT-PCR assay in breast cancers. The depiction analysis of co-expressed proteins was revealed using the STRING datasets. The cellular infiltration level analysis by the TISIDB and TIMER databases. The transwell assay was performed to analyze cellular migration and invasion. CCK-8 assay was used to evaluate cellular drug susceptibility for docetaxel treatment. Predicted targeted drugs in breast cancers by GSCA Lite database online.

Results

Kaplan-Meier plotter suggested that high expression of both E2F8 and MYBL2 in Basal-like subtype had a poor relapse-free survival. Functional enrichment results identified that apoptosis, cell cycle, and hormone ER pathway were represented the crucial regulation pathways by both E2F8 and MYBL2. In the meantime, database analysis indicated that high expression of E2F8 responded to chemotherapy, while those patients of high expression of MYBL2 responded to endocrinotherapy, and a positive correlation between the expression of E2F8 and PD-L1/CTLA4. Our cell line experiments confirmed the importance of E2F8 and MYBL2 in proliferation and chemotherapy sensitivity, possibly, the relationship with PD-L1. Additionally, we also observed that the up-regulation of E2F8 was accompanied with higher enrichments of CD4+ T cells and CD8+ T cells in breast cancers.

Conclusion

Taken together, our findings elucidated a prospective target in Basal-like breast cancer, providing underlying molecular biomarkers for the development of breast cancer treatment.

Relationship of micro-RNA, mRNA and eIF Expression in Tamoxifen-Adapted MCF-7 Breast Cancer Cells: Impact of miR-1972 on Gene Expression, Proliferation and Migration

Background: Tamoxifen-adapted MCF-7-Tam cells represent an in-vitro model for acquired tamoxifen resistance, which is still a problem in clinics. We here investigated the correlation of microRNA-, mRNA- and eukaryotic initiation factors (eIFs) expression in this model. Methods: MicroRNA- and gene expression were analyzed by nCounter and qRT-PCR technology; eIFs by Western blotting. Protein translation mode was determined using a reporter gene assay. Cells were transfected with a miR-1972-mimic. Results: miR-181b-5p,-3p and miR-455-5p were up-, miR-375, and miR-1972 down-regulated and are significant in survival analysis. About 5% of the predicted target genes were significantly altered. Pathway enrichment analysis suggested a contribution of the FoxO1 pathway. The ratio of polio-IRES driven to cap-dependent protein translation shifted towards cap-dependent initiation. Protein expression of eIF2A, -4G, -4H and -6 decreased, whereas eIF3H was higher in MCF-7-Tam. Significant correlations between tamoxifen-regulated miRNAs and eIFs were found in representative breast cancer cell lines. Transfection with a miR-1972-mimic reverses tamoxifen-induced expression for a subset of genes and increased proliferation in MCF-7, but reduced proliferation in MCF-7-Tam, especially in the presence of 4OH-tamoxifen. Migration was inhibited in MCF-7-Tam cells. Translation mode remained unaffected. Conclusions: miR-1972 contributes to the orchestration of gene-expression and physiological consequences of tamoxifen adaption.

The Breast Cancer Protein Co-Expression Landscape

Simple Summary

Proteins are among the most fundamental building blocks and molecular players behind the functions of cells and tissues. Their abundance and interaction patterns shape, to a large extent, what happens at the cellular and organ levels. This is also true regarding tumor tissues. In this work, we explored the patterns of abundance and co-occurrence of a large number of proteins in breast cancer cells and their healthy counterparts. We discovered the main differences and tried to see whether those differences may be associated with relevant aspects of the biology of these tumors. Our final goal is to provide information to empower cancer clinicians and pharmacologists to develop better diagnostic, prognostic, and therapeutic tools.

Abstract

Breast cancer is a complex phenotype (or better yet, several complex phenotypes) characterized by the interplay of a large number of cellular and biomolecular entities. Biological networks have been successfully used to capture some of the heterogeneity of intricate pathophenotypes, including cancer. Gene coexpression networks, in particular, have been used to study large-scale regulatory patterns. Ultimately, biological processes are carried out by proteins and their complexes. However, to date, most of the tumor profiling research has focused on the genomic and transcriptomic information. Here, we tried to expand this profiling through the analysis of open proteomic data via mutual information co-expression networks’ analysis. We could observe that there are distinctive biological processes associated with communities of these networks and how some transcriptional co-expression phenomena are lost at the protein level. These kinds of data and network analyses are a broad resource to explore cellular behavior and cancer research.

Introduction to Genomic Network Reconstruction for Cancer Research

High-throughput genomic technologies have revolutionized the study of cancer. Current research in oncology is now limited more for the capacity of analyzing and interpreting data, rather than the availability of data itself. Integrative approaches to obtain functional information from data are at the core of the disciplines gathered under the systems biology banner. In this context, network models have been used to study cancer, from the identification of key molecules involved in the disease to the discovery of functional alterations associated with specific manifestations of the disease.In this chapter, we describe the state of the art of network reconstruction from genomic data, with an emphasis in gene expression experiments. We explore the strengths and limitations of correlation, Bayesian, and information theoretic approaches to network reconstruction. We present tools that leverage the flexibility of network science to gain a deeper understanding of cancer biology.

Non-Coding RNAs Are Brokers in Breast Cancer Interactome Networks and Add Discrimination Power between Subtypes

Despite the power of high-throughput genomics, most non-coding RNA (ncRNA) biotypes remain hard to identify, characterize, and validate. This is a clear indication that intensive next-generation sequencing research has led to great efficiency and accuracy in detecting ncRNAs, but not in their functionalization. Computational scientists continue to support the discovery process by spotting significant data features (expression or mutational profiles), elucidating phenotype uncertainty, and delineating complex regulation landscapes for biological pathways and pathophysiological processes. With reference to transcriptome regulation dynamics in cancer, this work introduces a novel network-driven inference approach designed to reveal the potential role of computationally identified ncRNAs in discriminating between breast cancer (BC) subtypes beyond the traditional gene expression signatures. As heterogeneity cast in the subtypes is a characteristic of most cancers, the proposed approach is generalizable beyond BC. Expression profiles of a wide transcriptome spectrum were obtained for a number of BC patients (and controls) listed in TCGA and processed with RNA-Seq. The well-known PAM50 subtype signature was available for the samples and used to move from differentially expressed transcript profiles to subtype-specific biclusters associating gene patterns with patients. Co-expressed gene networks were then generated and annotations were provided, focusing on the biclusters with basal and luminal signatures. These were used to build template maps, i.e., networks in which to embed the ncRNAs and contextually functionalize them based on their interactors. This inference approach is able to assess the influence of ncRNAs at the level of BC subtype. Network topology was considered through the brokerage measure to account for disruptiveness effects induced by the removal of nodes corresponding to ncRNAs. Equivalently, it is shown that ncRNAs can act as brokers of network interactome dynamics, and removing them allows the refinement of subtype-related characteristics previously obtained by gene signatures only. The results of the study elucidate the role of pseudogenes in two major BC subtypes, considering the contextual annotations. Put into a wider perspective, ncRNA brokers may help predictive functionalization studies targeted to new disease phenotypes, for instance those linked to the tumor microenvironment or metabolism, or those specifically involving metastasis. Overall, the approach may represent an in silico prioritization strategy toward the systems identification of new diagnostic and prognostic biomarkers.

Dissecting the Role of N6-Methylandenosine-Related Long Non-coding RNAs Signature in Prognosis and Immune Microenvironment of Breast Cancer

Breast cancer (BC) represents a molecularly and clinically heterogeneous disease. Recent progress in immunotherapy has provided a glimmer of hope for several BC subtypes. The relationship between N6-methyladenosine (m6A) modification and long non-coding RNAs (LncRNAs) is still largely unexplored in BC. Here, with the intention to dissect the landscape of m6A-related lncRNAs and explore the immunotherapeutic value of the m6A-related lncRNA signature, we identified m6A-related lncRNAs by co-expression analysis from The Cancer Genome Atlas (TCGA) and stratified BC patients into different subgroups. Furthermore, we generated an m6A-related lncRNA prognostic signature. Four molecular subtypes were identified by consensus clustering. Cluster 3 preferentially had favorable prognosis, upregulated immune checkpoint expression, and high level of immune cell infiltration. Twenty-one m6A-related lncRNAs were applied to construct the m6A-related lncRNA model (m6A-LncRM). Survival analysis and receiver operating characteristic (ROC) curves further confirmed the prognostic value and prediction performance of m6A-LncRM. Finally, high- and low-risk BC subgroups displayed significantly different clinical features and immune cell infiltration status. Overall, our study systematically explored the prognostic value of the m6A-related LncRNAs and identified a high immunogenicity BC subtype. The proposed m6A-related LncRNA model might serve as a robust prognostic signature and attractive immunotherapeutic targets for BC treatment.

Kinectin 1 promotes the growth of triple-negative breast cancer via directly co-activating NF-kappaB/p65 and enhancing its transcriptional activity

Triple-negative breast cancer (TNBC) is the most challenging subtype of breast cancer. Various endeavor has been made to explore the molecular biology basis of TNBC. Herein, we reported a novel function of factor Kinectin 1 (KTN1) as a carcinogenic promoter in TNBC. KTN1 expression in TNBC was increased compared with adjacent tissues or luminal or Her2 subtypes of breast cancer, and TNBC patients with high KTN1 expression have poor prognosis. In functional studies, knockdown of KTN1 inhibited the proliferation and invasiveness of TNBC both in vitro and in vivo, while overexpression of KTN1 promoted cancer cell proliferation and invasiveness. RNA-seq analysis revealed that the interaction of cytokine-cytokine receptor, particularly CXCL8 gene, was upregulated by KTN1, which was supported by the further experiments. CXCL8 depletion inhibited the tumorigenesis and progression of TNBC. Additionally, rescue experiments validated that KTN1-mediated cell growth acceleration in TNBC was dependent on CXCL8 both in vitro and in vivo. Furthermore, it was found that KTN1 enhanced the phosphorylation of NF-κB/p65 protein at Ser536 site, and specifically bound to NF-κB/p65 protein in the nucleus and cytoplasm of cells. Moreover, the transcription of CXCL8 gene was directly upregulated by the complex of KTN1 and NF-κB/p65 protein. Taken together, our results elucidated a novel mechanism of KTN1 gene in TNBC tumorigenesis and progression. KTN1 may be a potential molecular target for the development of TNBC treatment.

Analysis of immune subtypes across the epithelial-mesenchymal plasticity spectrum

Epithelial-mesenchymal plasticity plays a critical role in many solid tumor types as a mediator of metastatic dissemination and treatment resistance. In addition, there is also a growing appreciation that the epithelial/mesenchymal status of a tumor plays a role in immune evasion and immune suppression. A deeper understanding of the immunological features of different tumor types has been facilitated by the availability of large gene expression datasets and the development of methods to deconvolute bulk RNA-Seq data. These resources have generated powerful new ways of characterizing tumors, including classification of immune subtypes based on differential expression of immunological genes. In the present work, we combine scoring algorithms to quantify epithelial-mesenchymal plasticity with immune subtype analysis to understand the relationship between epithelial plasticity and immune subtype across cancers. We find heterogeneity of epithelial-mesenchymal transition (EMT) status both within and between cancer types, with greater heterogeneity in the expression of EMT-related factors than of MET-related factors. We also find that specific immune subtypes have associated EMT scores and differential expression of immune checkpoint markers.

An Information Theoretical Multilayer Network Approach to Breast Cancer Transcriptional Regulation

Breast cancer is a complex, highly heterogeneous disease at multiple levels ranging from its genetic origins and molecular processes to clinical manifestations. This heterogeneity has given rise to the so-called intrinsic or molecular breast cancer subtypes. Aside from classification, these subtypes have set a basis for differential prognosis and treatment. Multiple regulatory mechanisms-involving a variety of biomolecular entities-suffer from alterations leading to the diseased phenotypes. Information theoretical approaches have been found to be useful in the description of these complex regulatory programs. In this work, we identified the interactions occurring between three main mechanisms of regulation of the gene expression program: transcription factor regulation, regulation via noncoding RNA, and epigenetic regulation through DNA methylation. Using data from The Cancer Genome Atlas, we inferred probabilistic multilayer networks, identifying key regulatory circuits able to (partially) explain the alterations that lead from a healthy phenotype to different manifestations of breast cancer, as captured by its molecular subtype classification. We also found some general trends in the topology of the multi-omic regulatory networks: Tumor subtype networks present longer shortest paths than their normal tissue counterpart; epigenomic regulation has frequently focused on genes enriched for certain biological processes; CpG methylation and miRNA interactions are often part of a regulatory core of conserved interactions. The use of probabilistic measures to infer information regarding theoretical-derived multilayer networks based on multi-omic high-throughput data is hence presented as a useful methodological approach to capture some of the molecular heterogeneity behind regulatory phenomena in breast cancer, and potentially other diseases.

Precision Oncology, Signaling and Anticancer Agents in Cancer Therapeutics

BACKGROUND

The global alliance for genomics and healthcare facilities provides innovational solutions to expedite research and clinical practices for complex and incurable health conditions. Precision oncology is an emerging field explicitly tailored to facilitate cancer diagnosis, prevention and treatment based on patients' genetic profile. Advancements in "omics" techniques, next-generation sequencing, artificial intelligence and clinical trial designs provide a platform for assessing the efficacy and safety of combination therapies and diagnostic procedures.

METHOD

Data were collected from Pubmed and Google scholar using keywords: "Precision medicine", "precision medicine and cancer", "anticancer agents in precision medicine" and reviewed comprehensively.

RESULTS

Personalized therapeutics including immunotherapy, cancer vaccines, serve as a groundbreaking solution for cancer treatment. Herein, we take a measurable view of precision therapies and novel diagnostic approaches targeting cancer treatment. The contemporary applications of precision medicine have also been described along with various hurdles identified in the successful establishment of precision therapeutics.

CONCLUSION

This review highlights the key breakthroughs related to immunotherapies, targeted anticancer agents, and target interventions related to cancer signaling mechanisms. The success story of this field in context to drug resistance, safety, patient survival and in improving quality of life is yet to be elucidated. We conclude that, in the near future, the field of individualized treatments may truly revolutionize the nature of cancer patient care.

Single-cell RNA-seq dissects the intratumoral heterogeneity of triple-negative breast cancer based on gene regulatory networks

Triple-negative breast cancer (TNBC) is a subtype of breast cancer with high intratumoral heterogeneity. Recent studies revealed that TNBC patients might comprise cells with distinct molecular subtypes. In addition, gene regulatory networks (GRNs) constructed based on single-cell RNA sequencing (scRNA-seq) data have demonstrated the significance for decoding the key regulators. We performed a comprehensive analysis of the GRNs for the intrinsic subtypes of TNBC patients using scRNA-seq. The copy number variations (CNVs) were inferred from scRNA-seq data and identified 545 malignant cells. The subtypes of the malignant cells were assigned based on the PAM50 model. The cell-cell communication analysis revealed that the macrophage plays a dominant role in the tumor microenvironment. Next, the GRN for each subtype was constructed through integrating gene co-expression and enrichment of transcription-binding motifs. Then, we identified the critical genes based on the centrality metrics of genes. Importantly, the critical gene ETV6 was ubiquitously upregulated in all subtypes, but it exerted diverse roles in each subtype through regulating different target genes. In conclusion, the construction of GRNs based on scRNA-seq data could help us to dissect the intratumoral heterogeneity and identify the critical genes of TNBC.