Integrated Analysis of RNA and DNA from the Phase III Trial CALGB 40601 Identifies Predictors of Response to Trastuzumab-Based Neoadjuvant Chemotherapy in HER2-Positive Breast Cancer

iGenSig-Rx: an integral genomic signature based white-box tool for modeling cancer therapeutic responses using multi-omics data

Multi-omics sequencing is poised to revolutionize clinical care in the coming decade. However, there is a lack of effective and interpretable genome-wide modeling methods for the rational selection of patients for personalized interventions. To address this, we present iGenSig-Rx, an integral genomic signature-based approach, as a transparent tool for modeling therapeutic response using clinical trial datasets. This method adeptly addresses challenges related to cross-dataset modeling by capitalizing on high-dimensional redundant genomic features, analogous to reinforcing building pillars with redundant steel rods. Moreover, it integrates adaptive penalization of feature redundancy on a per-sample basis to prevent score flattening and mitigate overfitting. We then developed a purpose-built R package to implement this method for modeling clinical trial datasets. When applied to genomic datasets for HER2 targeted therapies, iGenSig-Rx model demonstrates consistent and reliable predictive power across four independent clinical trials. More importantly, the iGenSig-Rx model offers the level of transparency much needed for clinical application, allowing for clear explanations as to how the predictions are produced, how the features contribute to the prediction, and what are the key underlying pathways. We anticipate that iGenSig-Rx, as an interpretable class of multi-omics modeling methods, will find broad applications in big-data based precision oncology. The R package is available: https://github.com/wangxlab/iGenSig-Rx .

PERCEPTION predicts patient response and resistance to treatment using single-cell transcriptomics of their tumors

Tailoring optimal treatment for individual cancer patients remains a significant challenge. To address this issue, we developed PERCEPTION (PERsonalized Single-Cell Expression-Based Planning for Treatments In ONcology), a precision oncology computational pipeline. Our approach uses publicly available matched bulk and single-cell (sc) expression profiles from large-scale cell-line drug screens. These profiles help build treatment response models based on patients' sc-tumor transcriptomics. PERCEPTION demonstrates success in predicting responses to targeted therapies in cultured and patient-tumor-derived primary cells, as well as in two clinical trials for multiple myeloma and breast cancer. It also captures the resistance development in patients with lung cancer treated with tyrosine kinase inhibitors. PERCEPTION outperforms published state-of-the-art sc-based and bulk-based predictors in all clinical cohorts. PERCEPTION is accessible at https://github.com/ruppinlab/PERCEPTION . Our work, showcasing patient stratification using sc-expression profiles of their tumors, will encourage the adoption of sc-omics profiling in clinical settings, enhancing precision oncology tools based on sc-omics.

Gene-Expression Profiling to Decipher Breast Cancer Inter- and Intratumor Heterogeneity

Breast cancer is heterogeneous and differs substantially across different tumors (intertumor heterogeneity) and even within an individual tumor (intratumor heterogeneity). Gene-expression profiling has considerably impacted our understanding of breast cancer biology. Four main "intrinsic subtypes" of breast cancer (i.e., luminal A, luminal B, HER2-enriched, and basal-like) have been consistently identified by gene expression, showing significant prognostic and predictive value in multiple clinical scenarios. Thanks to the molecular profiling of breast tumors, breast cancer is a paradigm of treatment personalization. Several standardized prognostic gene-expression assays are presently being used in the clinic to guide treatment decisions. Moreover, the development of single-cell-level resolution molecular profiling has allowed us to appreciate that breast cancer is also heterogeneous within a single tumor. There is an evident functional heterogeneity within the neoplastic and tumor microenvironment cells. Finally, emerging insights from these studies suggest a substantial cellular organization of neoplastic and tumor microenvironment cells, thus defining breast cancer ecosystems and highlighting the importance of spatial localizations.

Association of PIK3CA Mutation With Pathologic Complete Response and Outcome by Hormone Receptor Status and Intrinsic Subtype in Early-Stage ERBB2/HER2-Positive Breast Cancer

Importance

PIK3CA mutations may be associated with outcomes of patients with ERBB2/HER2-positive early breast cancer (EBC).

Objectives

To assess if PIK3CA mutations among patients with ERBB2/HER2-positive EBC are associated with treatment response and outcome, and if these associations vary by hormone receptor (HR) status or intrinsic molecular subtype (IMS).

Design, Setting, and Participants

This cohort study derived data on 184 patients from the phase 3 neoadjuvant Cancer and Leukemia Group B (CALGB) 40601 trial that enrolled patients with ERBB2/HER2-positive EBC in North America between January 1, 2008, and December 31, 2012. Participants received neoadjuvant paclitaxel with trastuzumab, lapatinib, or both. Statistical analysis was performed from March 23, 2022, to March 9, 2023.

Exposures

Gene expression profiling by RNA sequencing with Prediction Analysis of Microarray 50-determined IMS and PIK3CA mutations from whole-exome sequencing were obtained from pretreatment biopsies from 184 of 305 trial participants.

Main Outcomes and Measures

The primary end point was pathologic complete response (pCR) and the secondary end point of event-free survival (EFS). The association of PIK3CA mutations with pCR and EFS by HR status and IMS was estimated using logistic and Cox proportional hazards regression models.

Results

All 184 participants were women, with a median age of 49 years (range 24-75 years). A total of 121 participants (66%) had clinical stage II tumors; 32 (17%) had PIK3CA mutations, most frequently H1047R (38% [12 of 32]) and E545K (22% [7 of 32]). PIK3CA mutations were present in 20 of 102 cases of HR-positive EBC (20%) and 12 of 82 cases HR-negative EBC (15%) and varied by IMS (luminal B, 9 of 25 [36%]; luminal A, 2 of 21 [10%]; and ERBB2/HER2-enriched tumors, 19 of 102 [19%]). Pathologic complete response rates were lower in PIK3CA mutated than PIK3CA wild type in the overall population (34% [11 of 32] vs 49% [74 of 152]; P = .14) and were significantly different among those receiving trastuzumab (30% [7 of 23] vs 54% [63 of 117]; P = .045). At a median follow-up of 9 years, PIK3CA mutations were significantly associated with worse EFS in the overall cohort (hazard ratio, 2.58 [95% CI, 1.24-5.35]; P = .01), which persisted in a multivariable model including pCR, HR status, stage, and IMS (hazard ratio, 2.52 [95% CI, 1.16-5.47]; P = .02). The negative association of PIK3CA mutation was significant in HR-positive (hazard ratio, 3.60 [95% CI, 1.45-8.96]; P = .006) and luminal subtypes (hazard ratio, 4.84 [95% CI, 1.08-21.70]; P = .04), but not in nonluminal and HR-negative tumors.

Conclusions and Relevance

In ERBB2/HER2-positive EBC, PIK3CA mutations were associated with lower pCR rates and independently associated with worse long-term EFS. These findings appear to be associated with PIK3CA mutations in HR-positive and luminal EBC.

iGenSig-Rx: an integral genomic signature based white-box tool for modeling cancer therapeutic responses using multi-omics data

Multi-omics sequencing is expected to become clinically routine within the next decade and transform clinical care. However, there is a paucity of viable and interpretable genome-wide modeling methods that can facilitate rational selection of patients for tailored intervention. Here we develop an integral genomic signature-based method called iGenSig-Rx as a white-box tool for modeling therapeutic response based on clinical trial datasets with improved cross-dataset applicability and tolerance to sequencing bias. This method leverages high-dimensional redundant genomic features to address the challenges of cross-dataset modeling, a concept similar to the use of redundant steel rods to reinforce the pillars of a building. Using genomic datasets for HER2 targeted therapies, the iGenSig-Rx model demonstrates stable predictive power across four independent clinical trials. More importantly, the iGenSig-Rx model offers the level of transparency much needed for clinical application, allowing for clear explanations as to how the predictions are produced, how the features contribute to the prediction, and what are the key underlying pathways. We expect that iGenSig-Rx as a class of biologically interpretable multi-omics modeling methods will have broad applications in big-data based precision oncology. The R package is available: https://github.com/wangxlab/iGenSig-Rx. NOTE: the Github website will be released upon publication and the R package is available for review through google drive: https://drive.google.com/drive/folders/1KgecmUoon9-h2Dg1rPCyEGFPOp28Ols3?usp=sharing.

The Add-On Effect of Fluorouracil, Epirubicin, and Cyclophosphamide Regimens for Neoadjuvant Chemotherapy in Human Epidermal Receptor 2 (HER2)-Positive Breast Cancer: A Single-Center Retrospective Study

BACKGROUND

The addition of pertuzumab to trastuzumab in neoadjuvant chemotherapy (NAC) for anti-human epidermal receptor 2 (HER2) positive breast cancer has shown a significant improvement in the pathologic complete response (pCR) rate. However, the add-on effect of an anthracycline-based regimen (standard-of-care regimen) remains unclear. In this retrospective, observational study, participants received pertuzumab combination therapy as NAC for HER2-positive primary breast cancer.

METHODS

This study was conducted from January 1, 2020, to December 31, 2022. Patients who had not received at least three courses of pertuzumab owing to adverse events or those who had received preoperative radiotherapy were excluded.

RESULTS

The pCR rate was 35.3% (12/34 patients). The pCR group had a significantly higher percentage of histopathologic grade III (1/11 patients, p=0.030) and a significantly higher percentage of hormone receptor-negative patients (7/12 patients, p=0.015) than the non-pCR group. The non-pCR group had a significantly higher incidence of vascular invasion than the pCR group (7/22 patients, p=0.036). Menopausal status, stage, and ki-67 values were not significantly different between the two groups.

CONCLUSIONS

This study suggests an unlikely add-on effect of an anthracycline-based regimen for NAC in HER2-positive breast cancer. Moreover, our results support that the pCR rate is high in patients with hormone receptor-negative, HER2-positive breast cancer.

Deep (phospho)proteomics profiling of pre- treatment needle biopsies identifies signatures of treatment resistance in HER2+ breast cancer

Patients with early-stage HER2-overexpressing breast cancer struggle with treatment resistance in 20%-40% of cases. More information is needed to predict HER2 therapy response and resistance in vivo. In this study, we perform (phospho)proteomics analysis of pre-treatment HER2+ needle biopsies of early-stage invasive breast cancer to identify molecular signatures predictive of treatment response to trastuzumab, pertuzumab, and chemotherapy. Our data show that accurate quantification of the estrogen receptor (ER) and HER2 biomarkers, combined with the assessment of associated biological features, has the potential to enable better treatment outcome prediction. In addition, we identify cellular mechanisms that potentially precondition tumors to resist therapy. We find proteins with expression changes that correlate with resistance and constitute to a strong predictive signature for treatment success in our patient cohort. Our results highlight the multifactorial nature of drug resistance in vivo and demonstrate the necessity of deep tumor profiling.

Alterations in immune cell phenotype and cytotoxic capacity in HER2+ breast cancer patients receiving HER2-targeted neo-adjuvant therapy

BACKGROUND

The phase II neo-adjuvant clinical trial ICORG10-05 (NCT01485926) compared chemotherapy in combination with trastuzumab, lapatinib or both in patients with HER2+ breast cancer. We studied circulating immune cells looking for alterations in phenotype, genotype and cytotoxic capacity (direct and antibody-dependent cell-mediated cytotoxicity (ADCC)) in the context of treatment response.

METHODS

Peripheral blood mononuclear cells (PBMCs) were isolated from pre- (n = 41) and post- (n = 25) neo-adjuvant treatment blood samples. Direct/trastuzumab-ADCC cytotoxicity of patient-derived PBMCs against K562/SKBR3 cell lines was determined ex vivo. Pembrolizumab was interrogated in 21 pre-treatment PBMC ADCC assays. Thirty-nine pre-treatment and 21 post-treatment PBMC samples were immunophenotyped. Fc receptor genotype, tumour infiltrating lymphocyte (TIL) levels and oestrogen receptor (ER) status were quantified.

RESULTS

Treatment attenuated the cytotoxicity/ADCC of PBMCs. CD3+/CD4+/CD8+ T cells increased following therapy, while CD56+ NK cells/CD14+ monocytes/CD19+ B cells decreased with significant post-treatment immune cell changes confined to patients with residual disease. Pembrolizumab-augmented ex vivo PBMC ADCC activity was associated with residual disease, but not pathological complete response. Pembrolizumab-responsive PBMCs were associated with lower baseline TIL levels and ER+ tumours.

CONCLUSIONS

PBMCs display altered phenotype and function following completion of neo-adjuvant treatment. Anti-PD-1-responsive PBMCs in ex vivo ADCC assays may be a biomarker of treatment response.

Deciphering breast cancer: from biology to the clinic

Breast cancer remains a leading cause of cancer-related mortality in women, reflecting profound disease heterogeneity, metastasis, and therapeutic resistance. Over the last decade, genomic and transcriptomic data have been integrated on an unprecedented scale and revealed distinct cancer subtypes, critical molecular drivers, clonal evolutionary trajectories, and prognostic signatures. Furthermore, multi-dimensional integration of high-resolution single-cell and spatial technologies has highlighted the importance of the entire breast cancer ecosystem and the presence of distinct cellular "neighborhoods." Clinically, a plethora of new targeted therapies has emerged, now being rapidly incorporated into routine care. Resistance to therapy, however, remains a crucial challenge for the field.

Immune repertoire and responses to neoadjuvant TCHP therapy in HER2-positive breast cancer

Background

Despite the introduction of trastuzumab, pathologic complete response (pCR) is not attained in approximately 30-40% of Human epithelial growth factor receptor-2-positive breast cancer. Tumor-infiltrating lymphocytes (TIL) have been suggested as a predictive marker of treatment response, albeit not always effective. We investigated the relationship between trastuzumab, docetaxel, carboplatin, and pertuzumab (TCHP) treatment and immune repertoire as a treatment response predictor.

Design

In all, 35 cases were divided into two experimental groups: 10 and 25 cases in the preliminary and main experiments, respectively. In the preliminary experiment, the biopsy tissues before TCHP treatment and the surgical tissues after TCHP treatment were compared. In the main experiment, the biopsy tissues before TCHP treatment were compared according to the TCHP treatment response.

Methods

The T-cell repertoire for TRA, TRB, TRG, and TRD, and B-cell repertoire for immunoglobulin heavy, immunoglobulin kappa, and immunoglobulin lambda were evaluated. Whole transcriptome sequencing was also performed.

Results

In the preliminary experiment, the density and richness of the T-cell receptor (TCR) and B-cell receptor (BCR) repertoires decreased after treatment, regardless of TCHP response. In the main experiment, the Shannon's entropy index, density, and length of CDR3 of the TCR and BCR repertoires did not differ significantly in patients who did and did not achieve pCR. The pCR and non-pCR subgroups according to the level of TILs revealed that the non-pCR/lowTIL group had a higher proportion of low-frequency clones than the pCR/lowTIL group in TRA (non-pCR/lowTIL versus pCR/lowTIL, 0.01-0.1%, 63% versus 45.3%; <0.01%, 32.9% versus 51.8%, p < 0.001) and TRB (non-pCR/lowTIL versus pCR/lowTIL, 0.01-0.1%, 26.5% versus 14.7%; <0.01%, 72.0% versus 84.1%, p < 0.001).

Conclusions

The role of the diversity, richness, and density of the TCR and BCR repertoires as predictive markers for TCHP response was not identified. Compositions of low-frequency clones could be candidates for predictive factors of TCHP response; however, validation studies and further research are necessary.

Impact of RNA Signatures on pCR and Survival after 12-Week Neoadjuvant Pertuzumab plus Trastuzumab with or without Paclitaxel in the WSG-ADAPT HER2+/HR- Trial

PURPOSE

To identify associations of biological signatures and stromal tumor-infiltrating lymphocytes (sTIL) with pathological complete response (pCR; ypT0 ypN0) and survival in the Phase II WSG-ADAPT HER2+/HR- trial (NCT01817452).

EXPERIMENTAL DESIGN

Patients with cT1-cT4c, cN0-3 HER2+/HR- early breast cancer (EBC) were randomized to pertuzumab+trastuzumab (P+T, n = 92) or P+T+paclitaxel (n = 42). Gene expression signatures were analyzed in baseline biopsies using NanoString Breast Cancer 360 panel (n = 117); baseline and on-treatment (week 3) sTIL levels were available in 119 and 76 patients, respectively. Impacts of standardized gene expression signatures on pCR and invasive disease-free survival (iDFS) were estimated by logistic and Cox regression.

RESULTS

In all patients, ERBB2 [OR, 1.70; 95% confidence interval (CI), 1.08-2.67] and estrogen receptor (ER) signaling (OR, 1.72; 95% CI, 1.13-2.61) were favorable, whereas PTEN (OR, 0.57; 95% CI, 0.38-0.87) was unfavorable for pCR. After 60 months median follow-up, 13 invasive events occurred (P+T: n = 11, P+T+paclitaxel: n = 2), none following pCR. Gene signatures related to immune response (IR) and ER signaling were favorable for iDFS, all with similar HR about 0.43-0.55. These patterns were even more prominent in the neoadjuvant chemotherapy-free group, where additionally BRCAness signature was unfavorable (HR, 2.00; 95% CI, 1.04-3.84). IR signatures were strongly intercorrelated. sTILs (baseline/week 3/change) were not associated with pCR or iDFS, though baseline sTILs correlated positively with IR signatures.

CONCLUSIONS

Distinct gene signatures were associated with pCR versus iDFS in HER2+/HR- EBC. The potential role of IR in preventing recurrence suggests that patients with upregulated IR signatures could be candidates for de-escalation concepts in HER2+ EBC.

Translating transcriptomic findings from cancer model systems to humans through joint dimension reduction

Model systems are an essential resource in cancer research. They simulate effects that we can infer into humans, but come at a risk of inaccurately representing human biology. This inaccuracy can lead to inconclusive experiments or misleading results, urging the need for an improved process for translating model system findings into human-relevant data. We present a process for applying joint dimension reduction (jDR) to horizontally integrate gene expression data across model systems and human tumor cohorts. We then use this approach to combine human TCGA gene expression data with data from human cancer cell lines and mouse model tumors. By identifying the aspects of genomic variation joint-acting across cohorts, we demonstrate how predictive modeling and clinical biomarkers from model systems can be improved.

Personalized targeted therapy prescription in colorectal cancer using algorithmic analysis of RNA sequencing data

Background:

Overall survival of advanced colorectal cancer (CRC) patients remains poor, and gene expression analysis could potentially complement detection of clinically relevant mutations to personalize CRC treatments.

Methods:

We performed RNA sequencing of formalin-fixed, paraffin-embedded (FFPE) cancer tissue samples of 23 CRC patients and interpreted the data obtained using bioinformatic method Oncobox for expression-based rating of targeted therapeutics. Oncobox ranks cancer drugs according to the efficiency score calculated using target genes expression and molecular pathway activation data. The patients had primary and metastatic CRC with metastases in liver, peritoneum, brain, adrenal gland, lymph nodes and ovary. Two patients had mutations in NRAS, seven others had mutated KRAS gene. Patients were treated by aflibercept, bevacizumab, bortezomib, cabozantinib, cetuximab, crizotinib, denosumab, panitumumab and regorafenib as monotherapy or in combination with chemotherapy, and information on the success of totally 39 lines of therapy was collected.

Results:

Oncobox drug efficiency score was effective biomarker that could predict treatment outcomes in the experimental cohort (AUC 0.77 for all lines of therapy and 0.91 for the first line after tumor sampling). Separately for bevacizumab, it was effective in the experimental cohort (AUC 0.87) and in 3 independent literature CRC datasets, n = 107 (AUC 0.84–0.94). It also predicted progression-free survival in univariate (Hazard ratio 0.14) and multivariate (Hazard ratio 0.066) analyses. Difference in AUC scores evidences importance of using recent biosamples for the prediction quality.

Conclusion:

Our results suggest that RNA sequencing analysis of tumor FFPE materials may be helpful for personalizing prescriptions of targeted therapeutics in CRC.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12885-022-10177-3.

An integral genomic signature approach for tailored cancer therapy using genome-wide sequencing data

Low-cost multi-omics sequencing is expected to become clinical routine and transform precision oncology. Viable computational methods that can facilitate tailored intervention while tolerating sequencing biases are in high demand. Here we propose a class of transparent and interpretable computational methods called integral genomic signature (iGenSig) analyses, that address the challenges of cross-dataset modeling through leveraging information redundancies within high-dimensional genomic features, averaging feature weights to prevent overweighing, and extracting unbiased genomic information from large tumor cohorts. Using genomic dataset of chemical perturbations, we develop a battery of iGenSig models for predicting cancer drug responses, and validate the models using independent cell-line and clinical datasets. The iGenSig models for five drugs demonstrate predictive values in six clinical studies, among which the Erlotinib and 5-FU models significantly predict therapeutic responses in three studies, offering clinically relevant insights into their inverse predictive signature pathways. Together, iGenSig provides a computational framework to facilitate tailored cancer therapy based on multi-omics data.

N6-methyladenosine regulated FGFR4 attenuates ferroptotic cell death in recalcitrant HER2-positive breast cancer

Intrinsic and acquired anti-HER2 resistance remains a major hurdle for treating HER2-positive breast cancer. Using genome-wide CRISPR/Cas9 screening in vitro and in vivo, we identify FGFR4 as an essential gene following anti-HER2 treatment. FGFR4 inhibition enhances susceptibility to anti-HER2 therapy in resistant breast cancer. Mechanistically, m6A-hypomethylation regulated FGFR4 phosphorylates GSK-3β and activates β-catenin/TCF4 signaling to drive anti-HER2 resistance. Notably, suppression of FGFR4 dramatically diminishes glutathione synthesis and Fe²⁺ efflux efficiency via the β-catenin/TCF4-SLC7A11/FPN1 axis, resulting in excessive ROS production and labile iron pool accumulation. Ferroptosis, a unique iron-dependent form of oxidative cell death, is triggered after FGFR4 inhibition. Experiments involving patient-derived xenografts and organoids reveals a synergistic effect of anti-FGFR4 with anti-HER2 therapy in breast cancer with either intrinsic or acquired resistance. Together, these results pinpoint a mechanism of anti-HER2 resistance and provide a strategy for overcoming resistance via FGFR4 inhibition in recalcitrant HER2-positive breast cancer.

Anti-HER2 resistance causes treatment failure in HER2-positive breast cancers. Here the authors identify FGFR4 as one of the vulnerabilities of anti-HER2 resistant breast cancer and show that FGRR4 inhibition enhances sensitivity to anti-HER2 treatment in these resistant cells by triggering ferroptosis.

Gene Expression-Based Signature Can Predict Sorafenib Response in Kidney Cancer

Sorafenib is a tyrosine kinase inhibitory drug with multiple molecular specificities that is approved for clinical use in second-line treatments of metastatic and advanced renal cell carcinomas (RCCs). However, only 10–40% of RCC patients respond on sorafenib-containing therapies, and personalization of its prescription may help in finding an adequate balance of clinical efficiency, cost-effectiveness, and side effects. We investigated whether expression levels of known molecular targets of sorafenib in RCC can serve as prognostic biomarker of treatment response. We used Illumina microarrays to profile RNA expression in pre-treatment formalin-fixed paraffin-embedded (FFPE) samples of 22 metastatic or advanced RCC cases with known responses on next-line sorafenib monotherapy. Among them, nine patients showed partial response (PR), three patients—stable disease (SD), and 10 patients—progressive disease (PD) according to Response Evaluation Criteria In Solid Tumors (RECIST) criteria. We then classified PR + SD patients as “responders” and PD patients as “poor responders”. We found that gene signature including eight sorafenib target genes was congruent with the drug response characteristics and enabled high-quality separation of the responders and poor responders [area under a receiver operating characteristic curve (AUC) 0.89]. We validated these findings on another set of 13 experimental annotated FFPE RCC samples (for 2 PR, 1 SD, and 10 PD patients) that were profiled by RNA sequencing and observed AUC 0.97 for 8-gene signature as the response classifier. We further validated these results in a series of qRT-PCR experiments on the third experimental set of 12 annotated RCC biosamples (for 4 PR, 3 SD, and 5 PD patients), where 8-gene signature showed AUC 0.83.

Circulating Proteins Associated with Response and Resistance to Neoadjuvant Chemotherapy in HER2-Positive Breast Cancer

Simple Summary

The goal of this study was to find circulating proteins that can be easily sampled and incorporated into a clinical setting to improve predictive treatment response in HER2-positive breast cancer patients receiving neoadjuvant chemotherapy. We looked for potential biomarkers in serum, which we identified using two proteomics techniques: qualitative LC-MS/MS and a quantitative assay that assessed protein expression between responders and non-responders HER2-positive breast cancer patients to neoadjuvant chemotherapy.

Abstract

Despite the increasing use of neoadjuvant chemotherapy (NAC) in HER2-positive breast cancer (BC) patients, the clinical problem of predicting individual treatment response remains unanswered. Furthermore, the use of ineffective chemotherapeutic regimens should be avoided. Serum biomarker levels are being studied more and more for their ability to predict therapy response and aid in the development of personalized treatment regimens. This study aims to identify effective protein networks and biomarkers to predict response to NAC in HER2-positive BC patients through an exhaustive large-scale LC-MS/MS-based qualitative and quantitative proteomic profiling of serum samples from responders and non-responders. Serum samples from HER2-positive BC patients were collected before NAC and were processed by three methods (with and without nanoparticles). The qualitative analysis revealed differences in the proteomic profiles between responders and non-responders, mainly in proteins implicated in the complement and coagulation cascades and apolipoproteins. Qualitative analysis confirmed that three proteins (AFM, SERPINA1, APOD) were correlated with NAC resistance. In this study, we show that serum biomarker profiles can predict treatment response and outcome in the neoadjuvant setting. If these findings are further developed, they will be of significant clinical utility in the design of treatment regimens for individual BC patients.

Experimental and Meta-Analytic Validation of RNA Sequencing Signatures for Predicting Status of Microsatellite Instability

Microsatellite instability (MSI) is an important diagnostic and prognostic cancer biomarker. In colorectal, cervical, ovarian, and gastric cancers, it can guide the prescription of chemotherapy and immunotherapy. In laboratory diagnostics of susceptible tumors, MSI is routinely detected by the size of marker polymerase chain reaction products encompassing frequent microsatellite expansion regions. Alternatively, MSI status is screened indirectly by immunohistochemical interrogation of microsatellite binding proteins. RNA sequencing (RNAseq) profiling is an emerging source of data for a wide spectrum of cancer biomarkers. Recently, three RNAseq-based gene signatures were deduced for establishing MSI status in tumor samples. They had 25, 15, and 14 gene products with only one common gene. However, they were developed and tested on the incomplete literature of The Cancer Genome Atlas (TCGA) sampling and never validated experimentally on independent RNAseq samples. In this study, we, for the first time, systematically validated these three RNAseq MSI signatures on the literature colorectal cancer (CRC) (n = 619), endometrial carcinoma (n = 533), gastric cancer (n = 380), uterine carcinosarcoma (n = 55), and esophageal cancer (n = 83) samples and on the set of experimental CRC RNAseq samples (n = 23) for tumors with known MSI status. We found that all three signatures performed well with area under the curve (AUC) ranges of 0.94–1 for the experimental CRCs and 0.94–1 for the TCGA CRC, esophageal cancer, and uterine carcinosarcoma samples. However, for the TCGA endometrial carcinoma and gastric cancer samples, only two signatures were effective with AUC 0.91–0.97, whereas the third signature showed a significantly lower AUC of 0.69–0.88. Software for calculating these MSI signatures using RNAseq data is included.

RNA Sequencing Data for FFPE Tumor Blocks Can Be Used for Robust Estimation of Tumor Mutation Burden in Individual Biosamples

Tumor mutation burden (TMB) is a well-known efficacy predictor for checkpoint inhibitor immunotherapies. Currently, TMB assessment relies on DNA sequencing data. Gene expression profiling by RNA sequencing (RNAseq) is another type of analysis that can inform clinical decision-making and including TMB estimation may strongly benefit this approach, especially for the formalin-fixed, paraffin-embedded (FFPE) tissue samples. Here, we for the first time compared TMB levels deduced from whole exome sequencing (WES) and RNAseq profiles of the same FFPE biosamples in single-sample mode. We took TCGA project data with mean sequencing depth 23 million gene-mapped reads (MGMRs) and found 0.46 (Pearson)–0.59 (Spearman) correlation with standard mutation calling pipelines. This was converted into low (<10) and high (>10) TMB per megabase classifier with area under the curve (AUC) 0.757, and application of machine learning increased AUC till 0.854. We then compared 73 experimental pairs of WES and RNAseq profiles with lower (mean 11 MGMRs) and higher (mean 68 MGMRs) RNA sequencing depths. For higher depth, we observed ~1 AUC for the high/low TMB classifier and 0.85 (Pearson)–0.95 (Spearman) correlation with standard mutation calling pipelines. For the lower depth, the AUC was below the high-quality threshold of 0.7. Thus, we conclude that using RNA sequencing of tumor materials from FFPE blocks with enough coverage can afford for high-quality discrimination of tumors with high and low TMB levels in a single-sample mode.

Metabolic Stress Adaptations Underlie Mammary Gland Morphogenesis and Breast Cancer Progression

Breast cancers display dynamic reprogrammed metabolic activities as cancers develop from premalignant lesions to primary tumors, and then metastasize. Numerous advances focus on how tumors develop pro-proliferative metabolic signaling that differs them from adjacent, non-transformed epithelial tissues. This leads to targetable oncogene-driven liabilities among breast cancer subtypes. Other advances demonstrate how microenvironments trigger stress-response at single-cell resolution. Microenvironmental heterogeneities give rise to cell regulatory states in cancer cell spheroids in three-dimensional cultures and at stratified terminal end buds during mammary gland morphogenesis, where stress and survival signaling juxtapose. The cell-state specificity in stress signaling networks recapture metabolic evolution during cancer progression. Understanding lineage-specific metabolic phenotypes in experimental models is useful for gaining a deeper understanding of subtype-selective breast cancer metabolism.

MultiK: an automated tool to determine optimal cluster numbers in single-cell RNA sequencing data

Single-cell RNA sequencing (scRNA-seq) provides new opportunities to characterize cell populations, typically accomplished through some type of clustering analysis. Estimation of the optimal cluster number (K) is a crucial step but often ignored. Our approach improves most current scRNA-seq cluster methods by providing an objective estimation of the number of groups using a multi-resolution perspective. MultiK is a tool for objective selection of insightful Ks and achieves high robustness through a consensus clustering approach. We demonstrate that MultiK identifies reproducible groups in scRNA-seq data, thus providing an objective means to estimating the number of possible groups or cell-type populations present.

Algorithmically Deduced FREM2 Molecular Pathway Is a Potent Grade and Survival Biomarker of Human Gliomas

Gliomas are the most common malignant brain tumors with high mortality rates. Recently we showed that the FREM2 gene has a role in glioblastoma progression. Here we reconstructed the FREM2 molecular pathway using the human interactome model. We assessed the biomarker capacity of FREM2 expression and its pathway as the overall survival (OS) and progression-free survival (PFS) biomarkers. To this end, we used three literature and one experimental RNA sequencing datasets collectively covering 566 glioblastomas (GBM) and 1097 low-grade gliomas (LGG). The activation level of deduced FREM2 pathway showed strong biomarker characteristics and significantly outperformed the FREM2 expression level itself. For all relevant datasets, it could robustly discriminate GBM and LGG (p < 1.63 × 10-13, AUC > 0.74). High FREM2 pathway activation level was associated with poor OS in LGG (p < 0.001), and low PFS in LGG (p < 0.001) and GBM (p < 0.05). FREM2 pathway activation level was poor prognosis biomarker for OS (p < 0.05) and PFS (p < 0.05) in LGG with IDH mutation, for PFS in LGG with wild type IDH (p < 0.001) and mutant IDH with 1p/19q codeletion(p < 0.05), in GBM with unmethylated MGMT (p < 0.05), and in GBM with wild type IDH (p < 0.05). Thus, we conclude that the activation level of the FREM2 pathway is a potent new-generation diagnostic and prognostic biomarker for multiple molecular subtypes of GBM and LGG.

Transcriptome Analysis Identifies GATA3-AS1 as a Long Noncoding RNA Associated with Resistance to Neoadjuvant Chemotherapy in Locally Advanced Breast Cancer Patients

Breast cancer is one of the leading causes of mortality in women worldwide, and neoadjuvant chemotherapy has emerged as an option for the management of locally advanced breast cancer. Extensive efforts have been made to identify new molecular markers to predict the response to neoadjuvant chemotherapy. Transcripts that do not encode proteins, termed long noncoding RNAs (lncRNAs), have been shown to display abnormal expression profiles in different types of cancer, but their role as biomarkers in response to neoadjuvant chemotherapy has not been extensively studied. Herein, lncRNA expression was profiled using RNA sequencing in biopsies from patients who subsequently showed either response or no response to treatment. The GATA3-AS1 transcript was overexpressed in the nonresponder group and was the most stable feature when performing selection in multiple random forest models. GATA3-AS1 was experimentally validated by RT-qPCR in an extended group of 68 patients. Expression analysis confirmed that GATA3-AS1 is overexpressed primarily in patients who were nonresponsive to neoadjuvant chemotherapy, with a sensitivity of 92.9%, a specificity of 75.0%, and an area under the curve of approximately 0.90, as measured by receiver operating characteristic curve analysis. The statistical model was based on luminal B-like patients and adjusted by menopausal status and phenotype (odds ratio, 37.49; 95% CI, 6.74-208.42; P = 0.001); GATA3-AS1 was established as an independent predictor of response. Thus, lncRNA GATA3-AS1 is proposed as a potential predictive biomarker of nonresponse to neoadjuvant chemotherapy.

Sensitivity to targeted therapy differs between HER2-amplified breast cancer cells harboring kinase and helical domain mutations in PIK3CA

BACKGROUND

HER2-amplified breast cancer is a clinically defined subtype of breast cancer for which there are multiple viable targeted therapies. Resistance to these targeted therapies is a common problem, but the mechanisms by which resistance occurs remain incompletely defined. One mechanism that has been proposed is through mutation of genes in the PI3-kinase pathway. Intracellular signaling from the HER2 pathway can occur through PI3-kinase, and mutations of the encoding gene PIK3CA are known to be oncogenic. Mutations in PIK3CA co-occur with HER2-amplification in ~ 20% of cases within the HER2-amplified subtype.

METHODS

We generated isogenic knockin mutants of each PIK3CA hotspot mutation in HER2-amplified breast cancer cells using adeno-associated virus-mediated gene targeting. Isogenic clones were analyzed using a combinatorial drug screen to determine differential responses to HER2-targeted therapy. Western blot analysis and immunofluorescence uncovered unique intracellular signaling dynamics in cells resistant to HER2-targeted therapy. Subsequent combinatorial drug screens were used to explore neuregulin-1-mediated resistance to HER2-targeted therapy. Finally, results from in vitro experiments were extrapolated to publicly available datasets.

RESULTS

Treatment with HER2-targeted therapy reveals that mutations in the kinase domain (H1047R) but not the helical domain (E545K) increase resistance to lapatinib. Mechanistically, sustained AKT signaling drives lapatinib resistance in cells with the kinase domain mutation, as demonstrated by staining for the intracellular product of PI3-kinase, PIP3. This resistance can be overcome by co-treatment with an inhibitor to the downstream kinase AKT. Additionally, knockout of the PIP3 phosphatase, PTEN, phenocopies this result. We also show that neuregulin-1, a ligand for HER-family receptors, confers resistance to cells harboring either hotspot mutation and modulates response to combinatorial therapy. Finally, we show clinical evidence that the hotspot mutations have distinct expression profiles related to therapeutic resistance through analysis of TCGA and METABRIC data cohorts.

CONCLUSION

Our results demonstrate unique intracellular signaling differences depending on which mutation in PIK3CA the cell harbors. Only mutations in the kinase domain fully activate the PI3-kinase signaling pathway and maintain downstream signaling in the presence of HER2 inhibition. Moreover, we show there is potentially clinical importance in understanding both the PIK3CA mutational status and levels of neuregulin-1 expression in patients with HER2-amplified breast cancer treated with targeted therapy and that these problems warrant further pre-clinical and clinical testing.

Copy Number Aberration Analysis to Predict Response to Neoadjuvant Anti-HER2 Therapy: Results from the NeoALTTO Phase III Clinical Trial

PURPOSE

The heterogeneity of response to anti-HER2 agents represents a major challenge in patients with HER2-positive breast cancer. To better understand the sensitivity and resistance to trastuzumab and lapatinib, we investigated the role of copy number aberrations (CNA) in predicting pathologic complete response (pCR) and survival outcomes in the NeoALTTO trial.

EXPERIMENTAL DESIGN

The neoadjuvant phase III NeoALTTO trial enrolled 455 patients with HER2-positive early-stage breast cancer. DNA samples from 269 patients were assessed for genome-wide copy number profiling. Recurrent CNAs were found with GISTIC2.0.

RESULTS

CNA estimates were obtained for 184 patients included in NeoALTTO. Among those, matched transcriptome and whole-exome data were available for 154 and 181 patients, respectively. A significant association between gene copy number and pCR was demonstrated for ERBB2 amplification. Nevertheless, ERBB2 amplification ceased to be predictive once ERBB2 expression level was considered. GISTIC2.0 analysis revealed 159 recurrent CNA regions. Lower copy number levels of the 6q23-24 locus predicted absence of pCR in the whole cohort and in the estrogen receptor-positive subgroup. 6q23-24 deletion was significantly more frequent in TP53 wild-type (WT) compared with TP53-mutated, resulting in copy number levels significantly associated with lack of pCR only in the TP53 WT subgroup. Interestingly, a gene-ontology analysis highlighted several immune processes correlated to 6q23-24 copy number.

CONCLUSIONS

Our analysis identified ERBB2 copy number as well as 6q23-24 CNAs as predictors of response to anti-HER2-based treatment. ERBB2 expression outperformed ERBB2 amplification. The complexity of the 6q23-24 region warrants further investigation.

Clinical and Genetic Predictive Models for the Prediction of Pathological Complete Response to Optimize the Effectiveness for Trastuzumab Based Chemotherapy

Background

Trastuzumab shows excellent benefits for HER2+ breast cancer patients, although 20% treated remain unresponsive. We conducted a retrospective cohort study to optimize neoadjuvant chemotherapy and trastuzumab treatment in HER2+ breast cancer patients.

Methods

Six hundred patients were analyzed to identify clinical characteristics of those not achieving a pathological complete response (pCR) to develop a clinical predictive model. Available RNA sequence data was also reviewed to develop a genetic model for pCR.

Results

The pCR rate was 39.8% and pCR was associated with superior disease free survival and overall survival. ER negativity and PR negativity, higher HER2 IHC scores, higher Ki-67, and trastuzumab use were associated with improved pCR. Weekly paclitaxel and carboplatin had the highest pCR rate (46.70%) and the anthracycline+taxanes regimen had the lowest rate (11.11%). Four published GEO datasets were analyzed and a 10-gene model and immune signature for pCR were developed. Non-pCR patients were ER⁺PR⁺ and had a lower immune signature and gene model score. Hormone receptor status and immune signatures were independent predictive factors of pCR.

Conclusion

Hormone receptor status and a 10-gene model could predict pCR independently and may be applied for patient selection and drug effectiveness optimization.

Clinical and genomic assessment of PD-L1 SP142 expression in triple-negative breast cancer

PURPOSE

The SP142 PD-L1 assay is a companion diagnostic for atezolizumab in metastatic triple-negative breast cancer (TNBC). We strove to understand the biological, genomic, and clinical characteristics associated with SP142 PD-L1 positivity in TNBC patients.

METHODS

Using 149 TNBC formalin-fixed paraffin-embedded tumor samples, tissue microarray (TMA) and gene expression microarrays were performed in parallel. The VENTANA SP142 assay was used to identify PD-L1 expression from TMA slides. We next generated a gene signature reflective of SP142 status and evaluated signature distribution according to TNBCtype and PAM50 subtypes. A SP142 gene expression signature was identified and was biologically and clinically evaluated on the TNBCs of TCGA, other cohorts, and on other malignancies treated with immune checkpoint inhibitors (ICI).

RESULTS

Using SP142, 28.9% of samples were PD-L1 protein positive. The SP142 PD-L1-positive TNBC had higher CD8+ T cell percentage, stromal tumor-infiltrating lymphocyte levels, and higher rate of the immunomodulatory TNBCtype compared to PD-L1-negative samples. The recurrence-free survival was prolonged in PD-L1-positive TNBC. The SP142-guided gene expression signature consisted of 94 immune-related genes. The SP142 signature was associated with a higher pathologic complete response rate and better survival in multiple TNBC cohorts. In the TNBC of TCGA, this signature was correlated with lymphocyte-infiltrating signature scores, but not with tumor mutational burden or total neoantigen count. In other malignancies treated with ICIs, the SP142 genomic signature was associated with improved response and survival.

CONCLUSIONS

We provide multi-faceted evidence that SP142 PDL1-positive TNBC have immuno-genomic features characterized as highly lymphocyte-infiltrated and a relatively favorable survival.

LncRNA BCAR4 expression and predicts the clinical response to neoadjuvant chemotherapy in patients with locally advanced breast cancer

BACKGROUND

Neoadjuvant chemotherapy (NAC) is an important treatment for locally advanced breast cancer (LABC). However, there are no effective biomarkers to predict the efficacy. Therefore, there is an urgent need for new biomarkers to predict the response of LABC to NAC. LncRNA BCAR4 has been detected in a variety of malignant tumor tissues and used as a new biomarker for diagnosis and prognosis. However, LncRNA BCAR4 predicts the response of LABC to NAC is unclear.

OBJECTIVE

Explore the predictive effect of LncRNA BCAR4 on the efficacy of NAC for LABC in three different evaluation systems.

METHODS

First, the TCGA database was used to analyze the expression of LncRNA BCAR4 in 33 kinds of malignant tumors, and further explore its expression in breast cancer and its impact on the survival and prognosis of breast cancer. Furthermore, quantitative methods were used to measure the expression level of LncRNA BCAR4 in cancer tissues of 48 LABC patients, and the correlation between LncRNA BCAR4 and clinicopathological status and response to NAC under the evaluation system of 3, RECIST1.1, Miller-Payne (MP) score and whether it reaches pCR,was analyzed.

RESULTS

TCGA data analysis found that LncRNA is highly expressed in a variety of malignant tumor tissues, including breast cancer. And relatively low expression, the shorter the overall survival time of high expression patients. The high expression of LncRNA BCAR4 is related to the size of the tumor, and there are differences in expression between stage I and other stages, but there is no obvious correlation with the positive lymph node and hormone receptor status. Among the three evaluation systems, only in the RECIST 1.1 evaluation system LncRNA BCAR4 has a predictive effect on NAC for LABC. The expression of LncRNA BCAR4 has no significant correlation with clinical stage, Ki-67% and hormone receptor status, and has no significant correlation with whether patients with locally advanced breast cancer obtain pCR during neoadjuvant chemotherapy.

CONCLUSION

LncRNA BCAR4 is highly expressed in LABC tissues and may be an effective marker for predicting the efficacy of NAC for LABC.

Elastic Net Models Based on DNA Copy Number Variations Predicts Clinical Features, Expression Signatures, and Mutations in Lung Adenocarcinoma

In the precision medicine of lung adenocarcinoma, the identification and prediction of tumor phenotypes for specific biomolecular events are still not studied in depth. Various earlier researches sheds light on the close correlation between genetic expression signatures and DNA copy number variations (CNVs), for which analysis of CNVs provides valuable information about molecular and phenotypic changes in tumorigenesis. In this study, we propose a comprehensive analysis combining genome-wide association analysis and an Elastic Net Regression predictive model, focus on predicting the levels of many gene expression signatures in lung adenocarcinoma, based upon DNA copy number features alone. Additionally, we predicted many other key phenotypes, including clinical features (pathological stage), gene mutations, and protein expressions. These Elastic Net prediction methods can also be applied to other gene sets, thereby facilitating their use as biomarkers in monitoring therapy.

The tumor immune microenvironment and immune-related signature predict the chemotherapy response in patients with osteosarcoma

BACKGROUND

Genome-wide expression profiles have been shown to predict the response to chemotherapy. The purpose of this study was to develop a novel predictive signature for chemotherapy in patients with osteosarcoma.

METHODS

We analysed the relevance of immune cell infiltration and gene expression profiles of the tumor samples of good responders with those of poor responders from the TARGET and GEO databases. Immune cell infiltration was evaluated using a single-sample gene set enrichment analysis (ssGSEA) and the CIBERSORT algorithm between good and poor chemotherapy responders. Differentially expressed genes were identified based on the chemotherapy response. LASSO regression and binary logistic regression analyses were applied to select the differentially expressed immune-related genes (IRGs) and developed a predictive signature in the training cohort. A receiver operating characteristic (ROC) curve analysis was employed to assess and validate the predictive accuracy of the predictive signature in the validation cohort.

RESULTS

The analysis of immune infiltration showed a positive relationship between high-level immune infiltration and good responders, and T follicular helper cells and CD8 T cells were significantly more abundant in good responders with osteosarcoma. Two hundred eighteen differentially expressed genes were detected between good and poor responders, and a five IRGs panel comprising TNFRSF9, CD70, EGFR, PDGFD and S100A6 was determined to show predictive power for the chemotherapy response. A chemotherapy-associated predictive signature was developed based on these five IRGs. The accuracy of the predictive signature was 0.832 for the training cohort and 0.720 for the validation cohort according to ROC analysis.

CONCLUSIONS

The novel predictive signature constructed with five IRGs can be effectively utilized to predict chemotherapy responsiveness and help improve the efficacy of chemotherapy in patients with osteosarcoma.

Machine Learning Applicability for Classification of PAD/VCD Chemotherapy Response Using 53 Multiple Myeloma RNA Sequencing Profiles

Multiple myeloma (MM) affects ~500,000 people and results in ~100,000 deaths annually, being currently considered treatable but incurable. There are several MM chemotherapy treatment regimens, among which eleven include bortezomib, a proteasome-targeted drug. MM patients respond differently to bortezomib, and new prognostic biomarkers are needed to personalize treatments. However, there is a shortage of clinically annotated MM molecular data that could be used to establish novel molecular diagnostics. We report new RNA sequencing profiles for 53 MM patients annotated with responses on two similar chemotherapy regimens: bortezomib, doxorubicin, dexamethasone (PAD), and bortezomib, cyclophosphamide, dexamethasone (VCD), or with responses to their combinations. Fourteen patients received both PAD and VCD; six received only PAD, and 33 received only VCD. We compared profiles for the good and poor responders and found five genes commonly regulated here and in the previous datasets for other bortezomib regimens (all upregulated in the good responders): FGFR3, MAF, IGHA2, IGHV1-69, and GRB14. Four of these genes are linked with known immunoglobulin locus rearrangements. We then used five machine learning (ML) methods to build a classifier distinguishing good and poor responders for two cohorts: PAD + VCD (53 patients), and separately VCD (47 patients). We showed that the application of FloWPS dynamic data trimming was beneficial for all ML methods tested in both cohorts, and also in the previous MM bortezomib datasets. However, the ML models build for the different datasets did not allow cross-transferring, which can be due to different treatment regimens, experimental profiling methods, and MM heterogeneity.

Synthetic lethality-mediated precision oncology via the tumor transcriptome

Precision oncology has made significant advances, mainly by targeting actionable mutations in cancer driver genes. Aiming to expand treatment opportunities, recent studies have begun to explore the utility of tumor transcriptome to guide patient treatment. Here, we introduce SELECT (synthetic lethality and rescue-mediated precision oncology via the transcriptome), a precision oncology framework harnessing genetic interactions to predict patient response to cancer therapy from the tumor transcriptome. SELECT is tested on a broad collection of 35 published targeted and immunotherapy clinical trials from 10 different cancer types. It is predictive of patients' response in 80% of these clinical trials and in the recent multi-arm WINTHER trial. The predictive signatures and the code are made publicly available for academic use, laying a basis for future prospective clinical studies.

Spatial proteomic characterization of HER2-positive breast tumors through neoadjuvant therapy predicts response

The addition of HER2-targeted agents to neoadjuvant chemotherapy has dramatically improved pathological complete response (pCR) rates in early-stage, HER2-positive breast cancer. Nonetheless, up to 50% of patients have residual disease after treatment, while others are likely overtreated. Here, we performed multiplex spatial proteomic characterization of 122 samples from 57 HER2-positive breast tumors from the neoadjuvant TRIO-US B07 clinical trial sampled pre-treatment, after 14-21 d of HER2-targeted therapy and at surgery. We demonstrated that proteomic changes after a single cycle of HER2-targeted therapy aids the identification of tumors that ultimately undergo pCR, outperforming pre-treatment measures or transcriptomic changes. We further developed and validated a classifier that robustly predicted pCR using a single marker, CD45, measured on treatment, and showed that CD45-positive cell counts measured via conventional immunohistochemistry perform comparably. These results demonstrate robust biomarkers that can be used to enable the stratification of sensitive tumors early during neoadjuvant HER2-targeted therapy, with implications for tailoring subsequent therapy.

The Way of the Future: Personalizing Treatment Plans Through Technology

Advances in tissue analysis methods, image analysis, high-throughput molecular profiling, and computational tools increasingly allow us to capture and quantify patient-to patient variations that impact cancer risk, prognosis, and treatment response. Statistical models that integrate patient-specific information from multiple sources (e.g., family history, demographics, germline variants, imaging features) can provide individualized cancer risk predictions that can guide screening and prevention strategies. The precision, quality, and standardization of diagnostic imaging are improving through computer-aided solutions, and multigene prognostic and predictive tests improved predictions of prognosis and treatment response in various cancer types. A common theme across many of these advances is that individually moderately informative variables are combined into more accurate multivariable prediction models. Advances in machine learning and the availability of large data sets fuel rapid progress in this field. Molecular dissection of the cancer genome has become a reality in the clinic, and molecular target profiling is now routinely used to select patients for various targeted therapies. These technology-driven increasingly more precise and quantitative estimates of benefit versus risk from a given intervention empower patients and physicians to tailor treatment strategies that match patient values and expectations.

Pathologic and molecular responses to neoadjuvant trastuzumab and/or lapatinib from a phase II randomized trial in HER2-positive breast cancer (TRIO-US B07)

In this multicenter, open-label, randomized phase II investigator-sponsored neoadjuvant trial with funding provided by Sanofi and GlaxoSmithKline (TRIO-US B07, Clinical Trials NCT00769470), participants with early-stage HER2-positive breast cancer (N = 128) were recruited from 13 United States oncology centers throughout the Translational Research in Oncology network. Participants were randomized to receive trastuzumab (T; N = 34), lapatinib (L; N = 36), or both (TL; N = 58) as HER2-targeted therapy, with each participant given one cycle of this designated anti-HER2 therapy alone followed by six cycles of standard combination chemotherapy with the same anti-HER2 therapy. The primary objective was to estimate the rate of pathologic complete response (pCR) at the time of surgery in each of the three arms. In the intent-to-treat population, we observed similar pCR rates between T (47%, 95% confidence interval [CI] 30-65%) and TL (52%, 95% CI 38-65%), and a lower pCR rate with L (25%, 95% CI 13-43%). In the T arm, 100% of participants completed all protocol-specified treatment prior to surgery, as compared to 69% in the L arm and 74% in the TL arm. Tumor or tumor bed tissue was collected whenever possible pre-treatment (N = 110), after one cycle of HER2-targeted therapy alone (N = 89), and at time of surgery (N = 59). Higher-level amplification of HER2 and hormone receptor (HR)-negative status were associated with a higher pCR rate. Large shifts in the tumor, immune, and stromal gene expression occurred after one cycle of HER2-targeted therapy. In contrast to pCR rates, the L-containing arms exhibited greater proliferation reduction than T at this timepoint. Immune expression signatures increased in all arms after one cycle of HER2-targeted therapy, decreasing again by the time of surgery. Our results inform approaches to early assessment of sensitivity to anti-HER2 therapy and shed light on the role of the immune microenvironment in response to HER2-targeted agents.

Therapeutic Strategies for the Management of Hormone Receptor-Positive, Human Epidermal Growth Factor Receptor 2-Positive (HR+/HER2+) Breast Cancer: A Review of the Current Literature

Enormous advances have been made in the understanding and treatment of human epidermal growth factor receptor 2-positive breast cancer (HER2+ BC) in the last 30 years that have resulted in survival gains for affected patients. A growing body of evidence suggests that hormone receptor-positive (HR+)/HER2+ BC and HR-negative (HR-)/HER2+ BC are biologically different, with complex molecular bidirectional crosstalk between the estrogen receptor and HER2 pathway potentially affecting sensitivity to both HER2-targeted and endocrine therapy in patients with HR+/HER2+ BC. Subgroup analyses from trials enrolling patients with HER2+ BC and the results of clinical trials specifically designed to evaluate therapy in patients with HR+/HER2+ BC are helping to guide treatment decisions. In this context, encouraging results with strategies aimed at delaying or reversing drug resistance, including extended adjuvant therapy and the addition of drugs targeting alternative pathways, such as cyclin-dependent kinase (CDK) 4 and 6 inhibitors, have recently emerged. We have reached the point where tailoring the treatment according to risk and biology has become the paradigm in early BC. However, further clinical trials are needed that integrate translational research principles and identify and consider specific patient subgroups and biomarkers.

Survival, Pathologic Response, and Genomics in CALGB 40601 (Alliance), a Neoadjuvant Phase III Trial of Paclitaxel-Trastuzumab With or Without Lapatinib in HER2-Positive Breast Cancer

PURPOSE

CALGB 40601 assessed whether dual versus single human epidermal growth factor receptor 2 (HER2) -targeting drugs added to neoadjuvant chemotherapy increased pathologic complete response (pCR). Here, we report relapse-free survival (RFS), overall survival (OS), and gene expression signatures that predict pCR and survival.

PATIENTS AND METHODS

Three hundred five women with untreated stage II and III HER2-positive breast cancer were randomly assigned to receive weekly paclitaxel combined with trastuzumab plus lapatinib (THL), trastuzumab (TH), or lapatinib (TL). The primary end point was pCR, and secondary end points included RFS, OS, and gene expression analyses. mRNA sequencing was performed on 264 pretreatment samples.

RESULTS

One hundred eighteen patients were randomly allocated to THL, 120 to TH, and 67 to TL. At more than 7 years of follow-up, THL had significantly better RFS and OS than did TH (RFS hazard ratio, 0.32; 95% CI, 0.14 to 0.71; P = .005; OS hazard ratio, 0.34; 95% CI, 0.12 to 0.94; P = .037), with no difference between TH and TL. Of 688 previously described gene expression signatures, significant associations were found in 215 with pCR, 45 with RFS, and only 22 with both pCR and RFS (3.2%). Specifically, eight immune signatures were significantly correlated with a higher pCR rate and better RFS. Among patients with residual disease, the immunoglobulin G signature was an independent, good prognostic factor, whereas the HER2-enriched signature, which was associated with a higher pCR rate, showed a significantly shorter RFS.

CONCLUSION

In CALGB 40601, dual HER2-targeting resulted in significant RFS and OS benefits. Integration of intrinsic subtype and immune signatures allowed for the prediction of pCR and RFS, both overall and within the residual disease group. These approaches may provide means for rational escalation and de-escalation treatment strategies in HER2-positive breast cancer.

Cancer gene expression profiles associated with clinical outcomes to chemotherapy treatments

Background

Machine learning (ML) methods still have limited applicability in personalized oncology due to low numbers of available clinically annotated molecular profiles. This doesn’t allow sufficient training of ML classifiers that could be used for improving molecular diagnostics.

Methods

We reviewed published datasets of high throughput gene expression profiles corresponding to cancer patients with known responses on chemotherapy treatments. We browsed Gene Expression Omnibus (GEO), The Cancer Genome Atlas (TCGA) and Tumor Alterations Relevant for GEnomics-driven Therapy (TARGET) repositories.

Results

We identified data collections suitable to build ML models for predicting responses on certain chemotherapeutic schemes. We identified 26 datasets, ranging from 41 till 508 cases per dataset. All the datasets identified were checked for ML applicability and robustness with leave-one-out cross validation. Twenty-three datasets were found suitable for using ML that had balanced numbers of treatment responder and non-responder cases.

Conclusions

We collected a database of gene expression profiles associated with clinical responses on chemotherapy for 2786 individual cancer cases. Among them seven datasets included RNA sequencing data (for 645 cases) and the others – microarray expression profiles. The cases represented breast cancer, lung cancer, low-grade glioma, endothelial carcinoma, multiple myeloma, adult leukemia, pediatric leukemia and kidney tumors. Chemotherapeutics included taxanes, bortezomib, vincristine, trastuzumab, letrozole, tipifarnib, temozolomide, busulfan and cyclophosphamide.

NSABP B-41, a Randomized Neoadjuvant Trial: Genes and Signatures Associated with Pathologic Complete Response

PURPOSE

In NSABP B-41, pathologic complete response (pCR) was associated with prolonged survival among women with HER2-positive operable breast cancer treated with neoadjuvant chemotherapy and lapatinib, trastuzumab, or the combination. We used a large human breast cancer gene expression panel to select candidate prognostic biomarkers for pCR among women treated with trastuzumab in NSABP B-41.

PATIENTS AND METHODS

Eligible patients had a baseline preadjuvant treatment core biopsy sample, known pCR status, and no withdrawal of consent. We analyzed extracted RNA using the human nCounter Breast Cancer 360 gene expression panel. Gene counts were normalized to housekeeping genes and transformed into logarithmic scale with base 2. To screen for candidate genes and metagene signatures prognostic of pCR, we used univariate logistic regression. Variable selection was done by multivariable logistic regression with lasso regularization.

RESULTS

Analyses of data from 130 patients revealed that a composite of gene expression from 19 genes and one gene signature appeared to predict pCR in women with HER2-positive early-stage breast cancer undergoing neoadjuvant chemotherapy with trastuzumab-containing regimens. The identified genes are involved in important pathways such as epithelial-mesenchymal transition, adhesion and migration, estrogen receptor signaling, DNA damage and repair, apoptosis, and proliferation. The AUC from a 10-fold cross-validation on predicting pCR, with these 20 genomic markers in a logistic regression model, was 0.73.

CONCLUSIONS

The expression level of ERBB2, ESR1, and a few other genomic markers was highly predictive of pCR after trastuzumab-containing regimens. These findings need to be validated and calibrated in future studies.

B Cells and T Follicular Helper Cells Mediate Response to Checkpoint Inhibitors in High Mutation Burden Mouse Models of Breast Cancer

This study identifies mechanisms mediating responses to immune checkpoint inhibitors using mouse models of triple-negative breast cancer. By creating new mammary tumor models, we find that tumor mutation burden and specific immune cells are associated with response. Further, we developed a rich resource of single-cell RNA-seq and bulk mRNA-seq data of immunotherapy-treated and non-treated tumors from sensitive and resistant murine models. Using this, we uncover that immune checkpoint therapy induces T follicular helper cell activation of B cells to facilitate the anti-tumor response in these models. We also show that B cell activation of T cells and the generation of antibody are key to immunotherapy response and propose a new biomarker for immune checkpoint therapy. In total, this work presents resources of new preclinical models of breast cancer with large mRNA-seq and single-cell RNA-seq datasets annotated for sensitivity to therapy and uncovers new components of response to immune checkpoint inhibitors.

Serum miR-222-3p as a Double-Edged Sword in Predicting Efficacy and Trastuzumab-Induced Cardiotoxicity for HER2-Positive Breast Cancer Patients Receiving Neoadjuvant Target Therapy

Background: We aimed to explore whether the expression of serum miR-222-3p might contribute to early prediction of therapeutic response, clinical outcomes, and adverse events for HER2-positive breast cancer patients receiving neoadjuvant therapy (NAT).

Methods: A total of 65 HER2-positive breast cancer patients receiving NAT were analyzed. The concentration of serum miR-222-3p was detected by quantitative real-time PCR. Logistic regression analysis was used to identify the association of serum miR-222-3p with pathological complete response (pCR). The relationship of serum miR-222-3p with disease-free survival (DFS) and overall survival (OS) was examined via log-rank test and Cox proportional hazards analysis. The ordered logistic regression was applied to evaluate the association between serum miR-222-3p and adverse events.

Results: The miR-222-3p low group was more likely to achieve pCR [odds ratio (OR) = 0.258, P = 0.043]. The interaction between miR-222-3p and presenting Ki67 level was also detected for pCR (OR = 49.230, P_interaction = 0.025). The miR-222-3p low group was correlated with superior DFS (P = 0.029) and OS (P = 0.0037). The expression of serum miR-222-3p was the independent protective factor for trastuzumab-induced cardiotoxicity (P < 0.05) and anemia (P = 0.013).

Conclusions: Serum miR-222-3p is the potential factor to predict pCR, survival benefit and trastuzumab-induced cardiotoxicity for HER2-positive breast cancer patients receiving NAT.

RNA Sequencing in Comparison to Immunohistochemistry for Measuring Cancer Biomarkers in Breast Cancer and Lung Cancer Specimens

RNA sequencing is considered the gold standard for high-throughput profiling of gene expression at the transcriptional level. Its increasing importance in cancer research and molecular diagnostics is reflected in the growing number of its mentions in scientific literature and clinical trial reports. However, the use of different reagents and protocols for RNA sequencing often produces incompatible results. Recently, we published the Oncobox Atlas of RNA sequencing profiles for normal human tissues obtained from healthy donors killed in road accidents. This is a database of molecular profiles obtained using uniform protocol and reagents settings that can be broadly used in biomedicine for data normalization in pathology, including cancer. Here, we publish new original 39 breast cancer (BC) and 19 lung cancer (LC) RNA sequencing profiles obtained for formalin-fixed paraffin-embedded (FFPE) tissue samples, fully compatible with the Oncobox Atlas. We performed the first correlation study of RNA sequencing and immunohistochemistry-measured expression profiles for the clinically actionable biomarker genes in FFPE cancer tissue samples. We demonstrated high (Spearman's rho 0.65-0.798) and statistically significant (p < 0.00004) correlations between the RNA sequencing (Oncobox protocol) and immunohistochemical measurements for HER2/ERBB2, ER/ESR1 and PGR genes in BC, and for PDL1 gene in LC; AUC: 0.963 for HER2, 0.921 for ESR1, 0.912 for PGR, and 0.922 for PDL1. To our knowledge, this is the first validation that total RNA sequencing of archived FFPE materials provides a reliable estimation of marker protein levels. These results show that in the future, RNA sequencing can complement immunohistochemistry for reliable measurements of the expression biomarkers in FFPE cancer samples.

Sporadic activation of an oxidative stress-dependent NRF2-p53 signaling network in breast epithelial spheroids and premalignancies

Breast and mammary epithelial cells experience different local environments during tissue development and tumorigenesis. Microenvironmental heterogeneity gives rise to distinct cell regulatory states whose identity and importance are just beginning to be appreciated. Cellular states diversify when clonal three-dimensional (3D) spheroids are cultured in basement membrane, and one such state is associated with stress tolerance and poor response to anticancer therapeutics. Here, we found that this state was jointly coordinated by the NRF2 and p53 pathways, which were costabilized by spontaneous oxidative stress within 3D cultures. Inhibition of NRF2 or p53 individually disrupted some of the transcripts defining the regulatory state but did not yield a notable phenotype in nontransformed breast epithelial cells. In contrast, combined perturbation prevented 3D growth in an oxidative stress-dependent manner. By integrating systems models of NRF2 and p53 signaling in a single oxidative stress network, we recapitulated these observations and made predictions about oxidative stress profiles during 3D growth. NRF2 and p53 signaling were similarly coordinated in normal breast epithelial tissue and hormone-negative ductal carcinoma in situ lesions but were uncoupled in triple-negative breast cancer (TNBC), a subtype in which p53 is usually mutated. Using the integrated model, we correlated the extent of this uncoupling in TNBC cell lines with the importance of NRF2 in the 3D growth of these cell lines and their predicted handling of oxidative stress. Our results point to an oxidative stress tolerance network that is important for single cells during glandular development and the early stages of breast cancer.

RNA sequencing profiles and diagnostic signatures linked with response to ramucirumab in gastric cancer

Gastric cancer (GC) is the fifth-ranked cancer type by associated mortality. The proportion of early diagnosis is low, and most patients are diagnosed at the advanced stages. First-line therapy standardly includes fluoropyrimidines and platinum compounds with trastuzumab for HER2-positive cases. For recurrent disease, there are several alternative options including ramucirumab, a monoclonal therapeutic antibody that inhibits VEGF-mediated tumor angiogenesis by binding with VEGFR2, alone or in combination with other cancer drugs. However, overall response rate following ramucirumab or its combinations is 30%-80% of the patients, suggesting that personalization of drug prescription is needed to increase efficacy of treatment. We report here original tumor RNA sequencing profiles for 15 advanced GC patients linked with data on clinical response to ramucirumab or its combinations. Three genes showed differential expression in the tumors for responders versus nonresponders: CHRM3, LRFN1, and TEX15 Of them, CHRM3 was up-regulated in the responders. Using the bioinformatic platform Oncobox we simulated ramucirumab efficiency and compared output model results with actual tumor response data. An agreement was observed between predicted and real clinical outcomes (AUC ≥ 0.7). These results suggest that RNA sequencing may be used to personalize the prescription of ramucirumab for GC and indicate potential molecular mechanisms underlying ramucirumab resistance. The RNA sequencing profiles obtained here are fully compatible with the previously published Oncobox Atlas of Normal Tissue Expression (ANTE) data.

Towards personalized treatment for early stage HER2-positive breast cancer

Advances in HER2-targeted therapies have improved the survival of patients with HER2-positive breast cancer. The standard-of-care treatment for localized disease has been chemotherapy and 1 year of adjuvant HER2-targeted therapy, typically with the anti-HER2 antibody trastuzumab. Despite the effectiveness of this treatment, disease relapse occurs in a subset of patients; thus, focus has been placed on escalating treatment by either combining different HER2-targeted agents or extending the duration of HER2-targeted therapy. Indeed, dual HER2-targeted therapies and extended-duration anti-HER2 therapy, as well as adjuvant therapy with the anti-HER2 antibody-drug conjugate T-DM1, have all been approved for clinical use. Emerging evidence suggests, however, that some patients do not derive sufficient benefit from these additional therapies to offset the associated toxicities and/or costs. Similarly, the universal use of chemotherapy might not benefit all patients, and treatment de-escalation through omission of chemotherapy has shown promise in clinical trials and is currently being explored further. The future of precision medicine should therefore involve tailoring of therapy based on the genetics and biology of each tumour and the clinical characteristics of each patient. Predictive biomarkers that enable the identification of patients who will benefit from either escalated or de-escalated treatment will be crucial to this approach. In this Review, we summarize the available HER2-targeted agents and associated mechanisms of resistance, and describe the current therapeutic landscape of early stage HER2-positive breast cancer, focusing on strategies for treatment escalation or de-escalation.

Breast Cancer: A Molecularly Heterogenous Disease Needing Subtype-Specific Treatments

Breast cancer is the most commonly occurring cancer in women. There were over two-million new cases in world in 2018. It is the second leading cause of death from cancer in western countries. At the molecular level, breast cancer is a heterogeneous disease, which is characterized by high genomic instability evidenced by somatic gene mutations, copy number alterations, and chromosome structural rearrangements. The genomic instability is caused by defects in DNA damage repair, transcription, DNA replication, telomere maintenance and mitotic chromosome segregation. According to molecular features, breast cancers are subdivided in subtypes, according to activation of hormone receptors (estrogen receptor and progesterone receptor), of human epidermal growth factors receptor 2 (HER2), and or BRCA mutations. In-depth analyses of the molecular features of primary and metastatic breast cancer have shown the great heterogeneity of genetic alterations and their clonal evolution during disease development. These studies have contributed to identify a repertoire of numerous disease-causing genes that are altered through different mutational processes. While early-stage breast cancer is a curable disease in about 70% of patients, advanced breast cancer is largely incurable. However, molecular studies have contributed to develop new therapeutic approaches targeting HER2, CDK4/6, PI3K, or involving poly(ADP-ribose) polymerase inhibitors for BRCA mutation carriers and immunotherapy.

Genetic determinants of the molecular portraits of epithelial cancers

The ability to characterize and predict tumor phenotypes is crucial to precision medicine. In this study, we present an integrative computational approach using a genome-wide association analysis and an Elastic Net prediction method to analyze the relationship between DNA copy number alterations and an archive of gene expression signatures. Across breast cancers, we are able to quantitatively predict many gene signatures levels within individual tumors with high accuracy based upon DNA copy number features alone, including proliferation status and Estrogen-signaling pathway activity. We can also predict many other key phenotypes, including intrinsic molecular subtypes, estrogen receptor status, and TP53 mutation. This approach is also applied to TCGA Pan-Cancer, which identify repeatedly predictable signatures across tumor types including immune features in lung squamous and basal-like breast cancers. These Elastic Net DNA predictors could also be called from DNA-based gene panels, thus facilitating their use as biomarkers to guide therapeutic decision making.

Effective precision medicine strategies rely on the ability to predict tumour behaviour based on molecular characteristics. Here, the authors build models to predict multiple distinct gene expression patterns using DNA copy number alterations

Prospective Decision Analysis Study of Clinical Genomic Testing in Metastatic Breast Cancer: Impact on Outcomes and Patient Perceptions

PURPOSE

To evaluate the impact of targeted DNA sequencing on selection of cancer therapy for patients with metastatic breast cancer (MBC).

PATIENTS AND METHODS

In this prospective, single-center, single-arm trial, patients with MBC were enrolled within 10 weeks of starting a new therapy. At enrollment, tumor samples underwent next-generation sequencing for any of 315 cancer-related genes to high depth (> 500×) using FoundationOne CDx. Sequencing results were released to providers at the time of disease progression, and physician treatment recommendations were assessed via questionnaire. We evaluated three prespecified questions to assess patients' perceptions of genomic testing.

RESULTS

In all, 100 patients underwent genomic testing, with a median of five mutations (range, 0 to 13 mutations) detected per patient. Genomic testing revealed one or more potential therapies in 98% of patients (98 of 100), and 60% of patients (60 of 100) had one or more recommended treatments with level I/II evidence for actionability. Among the 94 genomic text reports that were released, there was physician questionnaire data for 87 patients (response rate, 92.6%) and 31.0% of patients (27 of 87) had treatment change recommended by their physician. Of these, 37.0% (10 of 27) received the treatment supported by genomic testing. We did not detect a statistically significant difference in time-to-treatment failure (log-rank P = .87) or overall survival (P = .71) among patients who had treatment change supported by genomic testing versus those who had no treatment change. For patients who completed surveys before and after genomic testing, there was a significant decrease in confidence of treatment success, specifically among patients who did not have treatment change supported by genomic testing (McNemar's test of agreement P = .001).

CONCLUSION

In this prospective study, genomic profiling of tumors in patients with MBC frequently identified potential treatments and resulted in treatment change in a minority of patients. Patients whose therapy was not changed on the basis of genomic testing seemed to have a decrease in confidence of treatment success.