Molecular subtyping of breast cancer intrinsic taxonomy with oligonucleotide microarray and NanoString nCounter

Homologous recombination deficiency reflects the heterogeneity and monitoring treatment response for patients with breast cancer

BACKGROUND

In breast cancer (BC), homologous recombination defect (HRD) is a common carcinogenic mechanism. It is meaningful to classify BC according to HRD biomarkers and to develop a platform for identifying BC molecular features, pathological features and therapeutic responses.

METHODS

In total, 109 HRD genes were collected and screened by univariate Cox regression analysis to determine the prognostic genes, which were used to construct a consensus matrix to identify BC subtype. Differentially expressed genes (DEGs) were filtered by the Limma package and screened by random forest analysis to build a model to analyze the immunotherapy response and sensitivity and prognosis of patients suffering from BC to different drugs.

RESULTS

Thirteen out of 109 HRD genes were prognostic genes of BC, and BC was classified into two subgroups based on their expression. Cluster 1 had a significantly backward survival outcome and a significantly higher adaptive immunity score relative to cluster 2. Six genes were identified by random forest analysis as factors for developing the model. The model provided a prediction called risk score, which showed a significant stratification effect on BC prognosis, immunotherapy response and IC50 values of 62 drugs.

CONCLUSIONS

In the present study, two HRD subtypes of BC were successfully identified, for which mutation and immunological features were determined. A model based on differential genes of HRD subtypes was established, which was a potential predictor of prognosis, immunotherapy response and drug sensitivity of BC.

Integration of bioinformatics and machine learning strategies identifies APM-related gene signatures to predict clinical outcomes and therapeutic responses for breast cancer patients

BACKGROUND

Tumor antigenicity and efficiency of antigen presentation jointly influence tumor immunogenicity, which largely determines the effectiveness of immune checkpoint blockade (ICB). However, the role of altered antigen processing and presentation machinery (APM) in breast cancer (BRCA) has not been fully elucidated.

METHODS

A series of bioinformatic analyses and machine learning strategies were performed to construct APM-related gene signatures to guide personalized treatment for BRCA patients. A single-sample gene set enrichment analysis (ssGSEA) algorithm and weighted gene co-expression network analysis (WGCNA) were combined to screen for BRCA-specific APM-related genes. The non-negative matrix factorization (NMF) algorithm was used to divide the cohort into different clusters and the fgsea algorithm was applied to investigate the altered signaling pathways. Random survival forest (RSF) and the least absolute shrinkage and selection operator (Lasso) Cox regression analysis were combined to construct an APM-related risk score (APMrs) signature to predict overall survival. Furthermore, a nomogram and decision tree were generated to improve predictive accuracy and risk stratification for individual patients. Based on Tumor Immune Dysfunction and Exclusion (TIDE) method, random forest (RF) and Lasso logistic regression model were combined to establish an APM-related immunotherapeutic response score (APMis). Finally, immune infiltration, immunomodulators, mutational patterns, and potentially applicable drugs were comprehensively analyzed in different APM-related risk groups. IHC staining was used to assess the expression of APM-related genes in clinical samples.

RESULTS

In this study, APMrs and APMis showed favorable performances in risk stratification and therapeutic prediction for BRCA patients. APMrs exhibited more powerful prognostic capacity and accurate survival prediction compared to conventional clinicopathological features. APMrs was closely associated with distinct mutational patterns, immune cell infiltration and immunomodulators expression. Furthermore, the two APM-related gene signatures were independently validated in external cohorts with prognosis or immunotherapeutic responses. Potential applicable drugs and targets were mined in the APMrs-high group. APM-related genes were further validated in our in-house samples.

CONCLUSION

The APM-related gene signatures established in our study could improve the personalized assessment of survival risk and guide ICB decision-making for BRCA patients.

Integrating Single-Cell RNA-Seq and Bulk RNA-Seq Data to Explore the Key Role of Fatty Acid Metabolism in Breast Cancer

Cancer immune escape is associated with the metabolic reprogramming of the various infiltrating cells in the tumor microenvironment (TME), and combining metabolic targets with immunotherapy shows great promise for improving clinical outcomes. Among all metabolic processes, lipid metabolism, especially fatty acid metabolism (FAM), plays a major role in cancer cell survival, migration, and proliferation. However, the mechanisms and functions of FAM in the tumor immune microenvironment remain poorly understood. We screened 309 fatty acid metabolism-related genes (FMGs) for differential expression, identifying 121 differentially expressed genes. Univariate Cox regression models in The Cancer Genome Atlas (TCGA) database were then utilized to identify the 15 FMGs associated with overall survival. We systematically evaluated the correlation between FMGs' modification patterns and the TME, prognosis, and immunotherapy. The FMGsScore was constructed to quantify the FMG modification patterns using principal component analysis. Three clusters based on FMGs were demonstrated in breast cancer, with three patterns of distinct immune cell infiltration and biological behavior. An FMGsScore signature was constructed to reveal that patients with a low FMGsScore had higher immune checkpoint expression, higher immune checkpoint inhibitor (ICI) scores, increased immune microenvironment infiltration, better survival advantage, and were more sensitive to immunotherapy than those with a high FMGsScore. Finally, the expression and function of the signature key gene NDUFAB1 were examined by in vitro experiments. This study significantly demonstrates the substantial impact of FMGs on the immune microenvironment of breast cancer, and that FMGsScores can be used to guide the prediction of immunotherapy efficacy in breast cancer patients. In vitro experiments, knockdown of the NDUFAB1 gene resulted in reduced proliferation and migration of MCF-7 and MDA-MB-231 cell lines.

Molecular Subtyping and Survival Analysis of Osteosarcoma Reveals Prognostic Biomarkers and Key Canonical Pathways

Osteosarcoma (OS) is a common bone malignancy in children and adolescents. Although histological subtyping followed by improved OS treatment regimens have helped achieve favorable outcomes, a lack of understanding of the molecular subtypes remains a challenge to characterize its genetic heterogeneity and subsequently to identify diagnostic and prognostic biomarkers for developing effective treatments. In the present study, global analysis of DNA methylation, and mRNA and miRNA gene expression in OS patient samples were correlated with their clinical characteristics. The mucin family of genes, MUC6, MUC12, and MUC4, were found to be highly mutated in the OS patients. Results revealed the enrichment of molecular pathways including Wnt signaling, Calcium signaling, and PI3K-Akt signaling in the OS tumors. Survival analyses showed that the expression levels of several genes such as RAMP1, CRIP1, CORT, CHST13, and DDX60L, miRNAs and lncRNAs were associated with survival of OS patients. Molecular subtyping using Cluster-Of-Clusters Analysis (COCA) for mRNA, lncRNA, and miRNA expression; DNA methylation; and mutation data from the TARGET dataset revealed two distinct molecular subtypes, each with a distinctive gene expression profile. Between the two subtypes, three upregulated genes, POP4, HEY1, CERKL, and seven downregulated genes, CEACAM1, ABLIM1, LTBP2, ISLR, LRRC32, PTPRF, and GPX3, associated with OS metastasis were found to be differentially regulated. Thus, the molecular subtyping results provide a strong basis for classification of OS patients that could be used to develop better prognostic treatment strategies.

The extended concurrent genes signature for disease-free survival in breast cancer

BACKGROUND/PURPOSE

Previously we had identified concurrent genes, which highlighted the interplay between copy number variation (CNV) and differential gene expression (GE) for Han Chinese breast cancers. The merit of the approach is to discovery biomarkers not identifiable by conventional GE only data, for which phenotype-correlation or gene variability is the criteria of gene selection.

MATERIALS AND METHODS

Thirty-one comparative genomic hybridization (CGH) and 83 GE microarrays were performed, with 29 breast cancers assayed from both platforms. Potential targets were revealed by Genomic Identification of Significant Targets in Cancer (GISTIC) from CGH arrays. Concurrent genes and genes with significant GISTIC scores were used to derive the extended concurrent genes signature, which was consensus from leading edge analysis across all studies and a supervised partial least square (PLS) regression predictive model of disease-free survival was constructed.

RESULTS

There were 1584 concurrent genes from 29 samples with both CGH and GE microarrays. Enriched concurrent genes sets for disease-free survival were identified independently from 83 GE arrays and another one with Han Chinese origin as well as three studies of Western origin. For five studies with disease-free survival follow up, prognostic discrepancy was observed between predicted high-risk and low-risk group patients.

CONCLUSION

We concluded that through parallel analyses of CGH and GE microarrays, the proposed extended concurrent gene expression signature can identify biomarkers with prognostic values.