Molecular Atlas of HER2+ Breast Cancer Cells Stimulated with EGF and HRG: Temporal Insights into Mechanisms in Trastuzumab Resistance

Using this Tool

This RShiny app was developed as an extension of Mukund, et al. (pending). It allows the user to easily generate their own plots using RPPA, RNAseq, and ATACseq data. Click on the question mark symbol next to each plot for more information.

Project Description

Trastuzumab therapy in HER2+ breast cancer patients have mixed success owing to acquired resistance to therapy. In this study, we investigate the cellular mechanisms underlying acquired resistance using -sensitive and -resistant cancer cells (BT474 and BT474R) treated with endogenous ligands EGF and HRG, across time. Our measurements probe early receptor organization through microscopy, signaling events through multi-omic measurements and assess the cellular energetic state through mitochondrial measurements. Our integrative analyses of these multimodal measurements highlight differential mechanisms in both cell lines in response to ligands. In BT474, an active PI3K-AKT-mTORC1 signaling contributes to an active mitochondrial bioenergetic state (glycolysis and lipid metabolism) for both ligands EGF and HRG. A HIF1A mediated increase in cellular bioenergetics is also seen. In BT474R, there is a ligand-dependent activation of signaling cascades. In EGF treated BT474R, an EGFR driven IRF1/STAT1/STAT2 activation with likely impact on cellular bioenergetics is seen, whilst an AR mediated alteration of lipid metabolism is pronounced after HRG treatment.

Study Design and Experimental Setup

We performed qSMLM, RPPA, RNA and ATAC sequencing to dynamically interrogate the cell state as defined by HER2 receptor localization, downstream protein signaling, its transcriptomic landscape and the accessibility of chromatin in response to ligands in BT474 and BT474R cultured cell lines. To capture the temporal dynamics of ligand treatments, samples were sequenced at four timepoints: T0 (baseline, pretreatment only), 1, 12 and 24 hours. Vehicle/water treated cells at each of the same timepoints were used as their respective controls. Protein lysates were additionally collected at 30 minutes for RPPA measurements, to capture the phosphorylated protein dynamics. All samples for the above measurements were obtained in duplicates (n=2) and differential analyses performed with respect to controls at a given timepoint. Given the time complexity of qSMLM experiments, T0 was used as baseline control (n=2-3). Sample preparation, extraction, and imaging/sequencing protocols for each of the aforementioned are presented in detail within the Methods section of the paper. Ligand concentrations were based on published studies and by titration experiments.

Contributors

This app was made by Darya Veraksa and Kavitha Mukund, Subramaniam Lab at UCSD.

RPPA Expression Heatmap

Heatmap showing log-normalized expression for selected antibody targets. Protein lysates from BT474 (Sensitive) and BT474R (Resistant) treated with EGF and HRG were RPPA sequenced at five timepoints: T0 (baseline, pretreatment only), 30 min, 1 hour, 12 hours, and 24 hours. Vehicle/water treated cells at each of the same timepoints were used as their respective controls.


RPPA Fold Change Dotplot

Dotplot of highest differentially expressed antibody targets between user-selected Condition 1 and Condition 2. Protein lysates from BT474 (Sensitive) and BT474R (Resistant) treated with EGF and HRG were RPPA sequenced at five timepoints: T0 (baseline, pretreatment only), 30 min, 1 hour, 12 hours, and 24 hours. Vehicle/water treated cells at each of the same timepoints were used as their respective controls.

Cell Line

Condition 1

Condition 2


RNASeq Volcano Plot

Volcano plot showing differentially expressed genes between user-selected Condition 1 and Condition 2. RNA from BT474 (Sensitive) and BT474R (Resistant) treated with EGF and HRG was sequenced at four timepoints: T0 (baseline, pretreatment only), 1 hour, 12 hours, and 24 hours. Vehicle/water treated cells at each of the same timepoints were used as their respective controls.

Cell Line

Time


RNASeq Boxplot

Boxplot showing rlog-normalized RNAseq expression values across two replicates. RNA from BT474 (Sensitive) and BT474R (Resistant) treated with EGF and HRG was sequenced at four timepoints: T0 (baseline, pretreatment only), 1 hour, 12 hours, and 24 hours. Vehicle/water treated cells at each of the same timepoints were used as their respective controls.


RNASeq Expression Heatmap

Heatmap showing average rlog-normalized expression values for selected genes. RNA from BT474 (Sensitive) and BT474R (Resistant) treated with EGF and HRG was sequenced at four timepoints: T0 (baseline, pretreatment only), 1 hour, 12 hours, and 24 hours. Vehicle/water treated cells at each of the same timepoints were used as their respective controls.


RNASeq Fold Change Heatmap

Heatmap showing significant fold-change (adj p<0.05) in counts between Control and selected conditions, for selected genes. If a gene is selected but not displayed for a certain condition, it did not have a statistically significant fold-change for the comparison with Control. RNA from BT474 (Sensitive) and BT474R (Resistant) treated with EGF and HRG was sequenced at four timepoints: T0 (baseline, pretreatment only), 1 hour, 12 hours, and 24 hours. Vehicle/water treated cells at each of the same timepoints were used as their respective controls.


ATACseq DAR Barplot

Barplot showing proportion of differentially accessible regions that were found in the promoter (1kb upstream or downstream of the TSS) and non-promoter regions of annotated genes.


ATACseq Chromosome Region Barplot

Barplot showing feature distribution based on chromosome region.


ATACseq Distance to TSS Barplot

Barplot showing feature distribution based on distance to transcription start site.

ATACseq Peak Profile Plot

Plot of average peak profiles. DNA from BT474 (Sensitive) and BT474R (Resistant) treated with EGF and HRG was ATAC sequenced at four timepoints: T0 (baseline, pretreatment only), 1 hour, 12 hours, and 24 hours. Vehicle/water treated cells at each of the same timepoints were used as their respective controls.

ATACSeq Enhancer Site Plot

Chromosome map of enhancer sites that regulate the chosen gene. DNA from BT474 (Sensitive) and BT474R (Resistant) treated with EGF and HRG was ATAC sequenced at four timepoints: T0 (baseline, pretreatment only), 1 hour, 12 hours, and 24 hours. Vehicle/water treated cells at each of the same timepoints were used as their respective controls.

ATACSeq Open Binding Site Plot

Map of open binding sites in selected region and chromosome. DNA from BT474 (Sensitive) and BT474R (Resistant) treated with EGF and HRG was ATAC sequenced at four timepoints: T0 (baseline, pretreatment only), 1 hour, 12 hours, and 24 hours. Vehicle/water treated cells at each of the same timepoints were used as their respective controls.