ABSTRACT
Impaired telomere length (TL) maintenance in ovarian tissue may play a pivotal role in the onset of epithelial ovarian cancer (OvC). TL in either target or surrogate tissue (blood) is currently being investigated for use as a predictor in anti-OvC therapy or as a biomarker of the disease progression, respectively. There is currently an urgent need for an appropriate approach to chemotherapy response prediction. We performed a monochrome multiplex qPCR measurement of TL in peripheral blood leukocytes (PBL) and tumor tissues of 209 OvC patients. The methylation status and gene expression of the shelterin complex and telomerase catalytic subunit (hTERT) were determined within tumor tissues by High-Throughput DNA methylation profiling and RNA sequencing (RNA-Seq) analysis, respectively. The patients sensitive to cancer treatment (n = 46) had shorter telomeres in PBL compared to treatment-resistant patients (n = 93; P = 0.037). In the patients with a different therapy response, transcriptomic analysis showed alterations in the peroxisome proliferator-activated receptor (PPAR) signaling pathway (q = 0.001). Moreover, tumor TL shorter than the median corresponded to better overall survival (OS) (P = 0.006). TPP1 gene expression was positively associated with TL in tumor tissue (P = 0.026). TL measured in PBL could serve as a marker of platinum therapy response in OvC patients. Additionally, TL determined in tumor tissue provides information on OvC patients' OS.
ABSTRACT
Chemotherapy is used as neoadjuvant therapy for all subgroups of breast cancer, including ER-positive, and HER2-negative cases. However, studies have suggested that using aromatase inhibitors combined with CDK4/6-inhibitors might be an appropriate alternative in selected patients. Thus, the NEOLETRIB trial evaluates the response of ER-positive, HER2-negative luminal A/B breast cancer to the combination of letrozole and ribociclib in the neoadjuvant setting. Comprehensive molecular biology procedures, including sequential single-cell RNA-sequencing of tumor biopsies, are performed during 6 months of treatment with extensive biobanking of blood samples, tumor biopsies and gut microbiome specimens. Our findings will hopefully contribute to an improved selection of patients who may benefit from this drug combination and give new insights into the intra-tumoral changes during this treatment.Trial registration number: NCT05163106 (ClinicalTrials.gov).
[Box: see text].
ABSTRACT
PURPOSE: Development of a computational biomarker to predict, prior to treatment, the response to CDK4/6 inhibition (CDK4/6i) in combination with endocrine therapy in patients with breast cancer. EXPERIMENTAL DESIGN: A mechanistic mathematical model that accounts for protein signaling and drug mechanisms of action was developed and trained on extensive, publicly available data from breast cancer cell lines. The model was built to provide a patient-specific response score based on the expression of six genes (CCND1, CCNE1, ESR1, RB1, MYC and CDKN1A). The model was validated in five independent cohorts of 148 patients in total with early-stage or advanced breast cancer treated with endocrine therapy and CDK4/6i. Response was measured either by evaluating Ki67 levels and PAM50 risk of relapse (ROR) after neoadjuvant treatment or by evaluating progression-free survival (PFS). RESULTS: The model showed significant association with patient´s outcomes in all five cohorts. The model predicted high Ki67 (area under the curve; AUC (95% confidence interval) of 0.80 (0.64 - 0.92), 0.81 (0.60 - 1.00) and 0.80 (0.65 - 0.93)) and high PAM50 ROR (AUC of 0.78 (0.64 - 0.89)). This observation was not obtained in patients treated with chemotherapy. In the other cohorts, patient stratification based on the model prediction was significantly associated with PFS (hazard ratio=2.92 (95% CI 1.08 - 7.86), p=0.034 and HR=2.16 (1.02 4.55), p=0.043). CONCLUSION: A mathematical modeling approach accurately predicts patient outcome following CDK4/6i plus endocrine therapy, which marks a step towards more personalized treatments in patients with Luminal B breast cancer.
ABSTRACT
Breast cancer (BCa) is a prevalent malignancy that predominantly affects women around the world. Somatic copy number alterations (CNAs) are tumor-specific amplifications or deletions of DNA segments that often drive BCa development and therapy resistance. Hence, the complex patterns of CNAs complement BCa classification systems. In addition, understanding the precise contributions of CNAs is essential for tailoring personalized treatment approaches. This review highlights how tumor evolution drives the acquisition of CNAs, which in turn shape the genomic landscapes of BCas. It also discusses advanced methodologies for identifying recurrent CNAs, studying CNAs in BCa and their clinical impact.
Subject(s)
Breast Neoplasms , DNA Copy Number Variations , Humans , Breast Neoplasms/genetics , Female , Genome, Human , Genomics/methodsABSTRACT
Treatment with the anti-angiogenic drug bevacizumab in addition to chemotherapy has shown efficacy for breast cancer in some clinical trials, but better biomarkers are needed to optimally select patients for treatment. Here, we present an omics approach where DNA methylation profiles are integrated with gene expression and results from proteomic data in breast cancer patients to predict response to therapy and pinpoint response-related epigenetic events. Fresh-frozen tumor biopsies taken before, during, and after treatment from human epidermal growth factor receptor 2 negative non-metastatic patients receiving neoadjuvant chemotherapy with or without bevacizumab were subjected to molecular profiling. Here, we report that DNA methylation at enhancer CpGs related to cell cycle regulation can predict response to chemotherapy and bevacizumab for the estrogen receptor positive subset of patients (AUC = 0.874), and we validated this observation in an independent patient cohort with a similar treatment regimen (AUC = 0.762). Combining the DNA methylation scores with the scores from a previously published protein signature resulted in a slight increase in the prediction performance (AUC = 0.784). We also show that tumors receiving the combination treatment underwent more extensive epigenetic alterations. Finally, we performed an integrative expression-methylation quantitative trait loci analysis on alterations in DNA methylation and gene expression levels, showing that the epigenetic alterations that occur during treatment are different between responders and non-responders and that these differences may be explained by the proliferation-epithelial-to-mesenchymal transition axis through the activity of grainyhead like transcription factor 2. Using tumor purity computed from copy number data, we developed a method for estimating cancer cell-specific methylation to confirm that the association to response reflects DNA methylation in cancer cells. Taken together, these results support the potential for clinical benefit of the addition of bevacizumab to chemotherapy when administered to the correct patients.
Subject(s)
Bevacizumab , Breast Neoplasms , DNA Methylation , Epigenesis, Genetic , Neoadjuvant Therapy , Humans , Bevacizumab/therapeutic use , Breast Neoplasms/drug therapy , Breast Neoplasms/genetics , Breast Neoplasms/pathology , Breast Neoplasms/metabolism , Female , Epigenesis, Genetic/drug effects , DNA Methylation/drug effects , Gene Expression Regulation, Neoplastic/drug effects , Proteomics , Middle AgedABSTRACT
Estrogen receptor positive breast cancer is frequently treated with anti-hormonal treatment such as aromatase inhibitors (AI). Interestingly, a high body mass index has been shown to have a negative impact on AI efficacy, most likely due to disturbances in steroid metabolism and adipokine production. Here, we propose a mathematical model based on a system of ordinary differential equations to investigate the effect of high-fat diet on tumor growth. We inform the model with data from mouse experiments, where the animals are fed with high-fat or control (normal) diet. By incorporating AI treatment with drug resistance into the model and by solving optimal control problems we found differential responses for control and high-fat diet. To the best of our knowledge, this is the first attempt to model optimal anti-hormonal treatment for breast cancer in the presence of drug resistance. Our results underline the importance of considering high-fat diet and obesity as factors influencing clinical outcomes during anti-hormonal therapies in breast cancer patients.
Subject(s)
Breast Neoplasms , Humans , Animals , Mice , Female , Breast Neoplasms/pathology , Drug Resistance, Neoplasm , Models, Biological , Mathematical Concepts , Aromatase Inhibitors/therapeutic use , Aromatase Inhibitors/pharmacology , DietABSTRACT
Anti-VEGF (vascular endothelial growth factor) treatment improves response rates, but not progression-free or overall survival in advanced breast cancer. It has been suggested that subgroups of patients may benefit from this treatment; however, the effects of adding anti-VEGF treatment to a standard chemotherapy regimen in breast cancer patients are not well studied. Understanding the effects of the anti-vascular treatment on tumor vasculature may provide a selection of patients that can benefit. The aim of this study was to study the vascular effect of bevacizumab using clinical dynamic contrast-enhanced MRI (DCE-MRI). A total of 70 women were randomized to receive either chemotherapy alone or chemotherapy with bevacizumab for 25 weeks. DCE-MRI was performed at baseline and at 12 and 25 weeks, and in addition 25 of 70 patients agreed to participate in an early MRI after one week. Voxel-wise pharmacokinetic analysis was performed using semi-quantitative methods and the extended Tofts model. Vascular architecture was assessed by calculating the fractal dimension of the contrast-enhanced images. Changes during treatment were compared with baseline and between the treatment groups. There was no significant difference in tumor volume at any point; however, DCE-MRI parameters revealed differences in vascular function and vessel architecture. Adding bevacizumab to chemotherapy led to a pronounced reduction in vascular DCE-MRI parameters, indicating decreased vascularity. At 12 and 25 weeks, the difference between the treatment groups is severely reduced.
ABSTRACT
FANCM germline protein truncating variants (PTVs) are moderate-risk factors for ER-negative breast cancer. We previously described the spectrum of FANCM PTVs in 114 European breast cancer cases. In the present, larger cohort, we report the spectrum and frequency of four common and 62 rare FANCM PTVs found in 274 carriers detected among 44,803 breast cancer cases. We confirmed that p.Gln1701* was the most common PTV in Northern Europe with lower frequencies in Southern Europe. In contrast, p.Gly1906Alafs*12 was the most common PTV in Southern Europe with decreasing frequencies in Central and Northern Europe. We verified that p.Arg658* was prevalent in Central Europe and had highest frequencies in Eastern Europe. We also confirmed that the fourth most common PTV, p.Gln498Thrfs*7, might be a founder variant from Lithuania. Based on the frequency distribution of the carriers of rare PTVs, we showed that the FANCM PTVs spectra in Southwestern and Central Europe were much more heterogeneous than those from Northeastern Europe. These findings will inform the development of more efficient FANCM genetic testing strategies for breast cancer cases from specific European populations.
ABSTRACT
We assessed the PREDICT v 2.2 for prognosis of breast cancer patients with pathogenic germline BRCA1 and BRCA2 variants, using follow-up data from 5453 BRCA1/2 carriers from the Consortium of Investigators of Modifiers of BRCA1/2 (CIMBA) and the Breast Cancer Association Consortium (BCAC). PREDICT for estrogen receptor (ER)-negative breast cancer had modest discrimination for BRCA1 carrier patients overall (Gönen & Heller unbiased concordance 0.65 in CIMBA, 0.64 in BCAC), but it distinguished clearly the high-mortality group from lower risk categories. In an analysis of low to high risk categories by PREDICT score percentiles, the observed mortality was consistently lower than the expected mortality, but the confidence intervals always included the calibration slope. Altogether, our results encourage the use of the PREDICT ER-negative model in management of breast cancer patients with germline BRCA1 variants. For the PREDICT ER-positive model, the discrimination was slightly lower in BRCA2 variant carriers (concordance 0.60 in CIMBA, 0.65 in BCAC). Especially, inclusion of the tumor grade distorted the prognostic estimates. The breast cancer mortality of BRCA2 carriers was underestimated at the low end of the PREDICT score distribution, whereas at the high end, the mortality was overestimated. These data suggest that BRCA2 status should also be taken into consideration with tumor characteristics, when estimating the prognosis of ER-positive breast cancer patients.
ABSTRACT
A limited number of studies are devoted to regulating TRIP6 expression in cancer. Hence, we aimed to unveil the regulation of TRIP6 expression in MCF-7 breast cancer cells (with high TRIP6 expression) and taxane-resistant MCF-7 sublines (manifesting even higher TRIP6 expression). We found that TRIP6 transcription is regulated primarily by the cyclic AMP response element (CRE) in hypomethylated proximal promoters in both taxane-sensitive and taxane-resistant MCF-7 cells. Furthermore, in taxane-resistant MCF-7 sublines, TRIP6 co-amplification with the neighboring ABCB1 gene, as witnessed by fluorescence in situ hybridization (FISH), led to TRIP6 overexpression. Ultimately, we found high TRIP6 mRNA levels in progesterone receptor-positive breast cancer and samples resected from premenopausal women.
Subject(s)
ATP Binding Cassette Transporter, Subfamily B , Drug Resistance, Neoplasm , LIM Domain Proteins , Neoplasms , Female , Humans , Adaptor Proteins, Signal Transducing/genetics , ATP Binding Cassette Transporter, Subfamily B/genetics , Cyclic AMP , Drug Resistance, Neoplasm/genetics , In Situ Hybridization, Fluorescence , LIM Domain Proteins/genetics , MCF-7 Cells , Neoplasms/genetics , Response Elements , Taxoids , Transcription Factors/geneticsABSTRACT
Human epidermal growth factor receptor 2-positive (HER2+) breast cancer is an aggressive subtype of this disease. Targeted treatment has improved outcome, but there is still a need for new therapeutic strategies as some patients respond poorly to treatment. Our aim was to identify compounds that substantially affect viability in HER2+ breast cancer cells in response to combinatorial treatment. We performed a high-throughput drug screen of 278 compounds in combination with trastuzumab and lapatinib using two HER2+ breast cancer cell lines (KPL4 and SUM190PT). The most promising drugs were validated in vitro and in vivo, and downstream molecular changes of the treatments were analyzed. The screen revealed multiple drugs that could be used in combination with lapatinib and/or trastuzumab. The Src-inhibitor dasatinib showed the largest combinatorial effect together with lapatinib in the KPL4 cell line compared to treatment with dasatinib alone (p < 0.01). In vivo, only lapatinib significantly reduced tumor growth (p < 0.05), whereas dasatinib alone, or in combination with lapatinib, did not show significant effects. Protein analyses of the treated xenografts showed significant alterations in protein levels compared to untreated controls, suggesting that all drugs reached the tumor and exerted a measurable effect. In silico analyses suggested activation of apoptosis and reduced activity of survival pathways by all treatments, but the opposite pattern was observed for the combinatorial treatment compared to lapatinib alone.
Subject(s)
Antineoplastic Agents , Breast Neoplasms , Humans , Female , Breast Neoplasms/pathology , Lapatinib/pharmacology , Lapatinib/therapeutic use , Dasatinib/pharmacology , Dasatinib/therapeutic use , Antineoplastic Agents/pharmacology , Antineoplastic Agents/therapeutic use , Quinazolines/pharmacology , Quinazolines/therapeutic use , Receptor, ErbB-2/metabolism , Antibodies, Monoclonal, Humanized/therapeutic use , Cell Line, Tumor , Trastuzumab/therapeutic useABSTRACT
Evidence from literature, including the BRIDGES study, indicates that germline protein truncating variants (PTVs) in FANCM confer moderately increased risk of ER-negative and triple-negative breast cancer (TNBC), especially for women with a family history of the disease. Association between FANCM missense variants (MVs) and breast cancer risk has been postulated. In this study, we further used the BRIDGES study to test 689 FANCM MVs for association with breast cancer risk, overall and in ER-negative and TNBC subtypes, in 39,885 cases (7566 selected for family history) and 35,271 controls of European ancestry. Sixteen common MVs were tested individually; the remaining rare 673 MVs were tested by burden analyses considering their position and pathogenicity score. We also conducted a meta-analysis of our results and those from published studies. We did not find evidence for association for any of the 16 variants individually tested. The rare MVs were significantly associated with increased risk of ER-negative breast cancer by burden analysis comparing familial cases to controls (OR = 1.48; 95% CI 1.07-2.04; P = 0.017). Higher ORs were found for the subgroup of MVs located in functional domains or predicted to be pathogenic. The meta-analysis indicated that FANCM MVs overall are associated with breast cancer risk (OR = 1.22; 95% CI 1.08-1.38; P = 0.002). Our results support the definition from previous analyses of FANCM as a moderate-risk breast cancer gene and provide evidence that FANCM MVs could be low/moderate risk factors for ER-negative and TNBC subtypes. Further genetic and functional analyses are necessary to clarify better the increased risks due to FANCM MVs.
Subject(s)
Breast Neoplasms , DNA Helicases , Humans , Female , Breast Neoplasms/genetics , DNA Helicases/genetics , Triple Negative Breast Neoplasms/genetics , Genetic Predisposition to DiseaseABSTRACT
Digital analysis of pathology whole-slide images has been recently gaining interest in the context of cancer diagnosis and treatment. In particular, deep learning methods have demonstrated significant potential in supporting pathology analysis, recently detecting molecular traits never before recognized in pathology H&E whole-slide images (WSIs). Alongside these advancements in the digital analysis of WSIs, it is becoming increasingly evident that both spatial and overall tumor heterogeneity may be significant determinants of cancer prognosis and treatment outcome. In this chapter, we describe methods that aim to support these two elements. We describe both an end-to-end deep learning pipeline for producing limited spatial transcriptomics from WSIs with associated bulk gene expression data, as well as an algorithm for quantifying spatial tumor heterogeneity based on the results of this pipeline.
Subject(s)
Neoplasms , Humans , Neoplasms/diagnosis , Neoplasms/genetics , Phenotype , Algorithms , Microscopy/methodsABSTRACT
The analysis of whole genomes of pan-cancer data sets provides a challenge for researchers, and we contribute to the literature concerning the identification of robust subgroups with clear biological interpretation. Specifically, we tackle this unsupervised problem via a novel rank-based Bayesian clustering method. The advantages of our method are the integration and quantification of all uncertainties related to both the input data and the model, the probabilistic interpretation of final results to allow straightforward assessment of the stability of clusters leading to reliable conclusions, and the transparent biological interpretation of the identified clusters since each cluster is characterized by its top-ranked genomic features. We applied our method to RNA-seq data from cancer samples from 12 tumor types from the Cancer Genome Atlas. We identified a robust clustering that mostly reflects tissue of origin but also includes pan-cancer clusters. Importantly, we identified three pan-squamous clusters composed of a mix of lung squamous cell carcinoma, head and neck squamous carcinoma, and bladder cancer, with different biological functions over-represented in the top genes that characterize the three clusters. We also found two novel subtypes of kidney cancer that show different prognosis, and we reproduced known subtypes of breast cancer. Taken together, our method allows the identification of robust and biologically meaningful clusters of pan-cancer samples.
Subject(s)
Breast Neoplasms , Carcinoma, Squamous Cell , Head and Neck Neoplasms , Humans , Female , Transcriptome , Bayes Theorem , Carcinoma, Squamous Cell/genetics , Breast Neoplasms/genetics , Cluster AnalysisABSTRACT
Most cancer alterations occur in the noncoding portion of the human genome, where regulatory regions control gene expression. The discovery of noncoding mutations altering the cells' regulatory programs has been limited to few examples with high recurrence or high functional impact. Here, we show that transcription factor binding sites (TFBSs) have similar mutation loads to those in protein-coding exons. By combining cancer somatic mutations in TFBSs and expression data for protein-coding and miRNA genes, we evaluate the combined effects of transcriptional and post-transcriptional alterations on the regulatory programs in cancers. The analysis of seven TCGA cohorts culminates with the identification of protein-coding and miRNA genes linked to mutations at TFBSs that are associated with a cascading trans-effect deregulation on the cells' regulatory programs. Our analyses of cis-regulatory mutations associated with miRNAs recurrently predict 12 mature miRNAs (derived from 7 precursors) associated with the deregulation of their target gene networks. The predictions are enriched for cancer-associated protein-coding and miRNA genes and highlight cis-regulatory mutations associated with the dysregulation of key pathways associated with carcinogenesis. By combining transcriptional and post-transcriptional regulation of gene expression, our method predicts cis-regulatory mutations related to the dysregulation of key gene regulatory networks in cancer patients.
Subject(s)
MicroRNAs , Neoplasms , Humans , Gene Expression Regulation , Neoplasms/genetics , Mutation , MicroRNAs/physiology , Gene Regulatory NetworksABSTRACT
Epithelial ovarian carcinoma (EOC) is known for high mortality due to diagnosis at advanced stages and frequent therapy resistance. Previous findings suggested that the DNA repair system is involved in the therapeutic response of cancer patients and DNA repair genes are promising targets for novel therapies. This study aimed to address complex inter-relations among gene expression levels, methylation profiles, and somatic mutations in DNA repair genes and EOC prognosis and therapy resistance status. We found significant associations of DUT expression with the presence of peritoneal metastases in EOC patients. The high-grade serous EOC subtype was enriched with TP53 mutations compared to other subtypes. Furthermore, somatic mutations in XPC and PRKDC were significantly associated with worse overall survival of EOC patients, and higher FAAP20 expression in platinum-resistant than platinum-sensitive patients was observed. We found higher methylation of RAD50 in platinum-resistant than in platinum-sensitive patients. Somatic mutations in BRCA1 and RAD9A were significantly associated with higher RBBP8 methylation in platinum-sensitive compared to platinum-resistant EOC patients. In conclusion, we discovered associations of several candidate genes from the DNA repair pathway with the prognosis and platinum resistance status of EOC patients, which deserve further validation as potential predictive biomarkers.
ABSTRACT
Kinase activity is frequently altered in renal cell carcinoma (RCC), and tyrosine kinase inhibitors (TKIs) are part of the standard treatment strategy in patients with metastatic disease. However, there are still no established biomarkers to predict clinical benefits of a specific TKI. Here, we performed protein tyrosine kinase (PTK) profiling using PamChip® technology. The aim of this study was to identify differences in PTK activity between normal and malignant kidney tissue obtained from the same patient, and to investigate the inhibitory effects of TKIs frequently used in the clinics: sunitinib, pazopanib, cabozantinib and tivozanib. Briefly, our results showed that 36 kinase substrates differs (FDR < 0.05) between normal and cancer kidney tissue, where members of the Src family kinases and the phosphoinositide-3-kinase (PI3K) pathway exhibit high activity in renal cancer. Furthermore, ex vivo treatment of clear cell RCC with TKIs revealed that pathways such as Rap1, Ras and PI3K pathways were strongly inhibited, whereas the neurotrophin pathway had increased activity upon TKI addition. In our assay, tivozanib and cabozantinib exhibited greater inhibitory effects on PTK activity compared to sunitinib and pazopanib, implying they might be better suitable as TKIs for selected RCC patients.
Subject(s)
Carcinoma, Renal Cell , Protein Kinase Inhibitors , Protein-Tyrosine Kinases , Anilides , Carcinoma, Renal Cell/drug therapy , Carcinoma, Renal Cell/metabolism , Carcinoma, Renal Cell/pathology , Humans , Indazoles , Kidney Neoplasms/drug therapy , Kidney Neoplasms/metabolism , Kidney Neoplasms/pathology , Nerve Growth Factors , Phenylurea Compounds , Phosphatidylinositol 3-Kinases , Phosphatidylinositols/therapeutic use , Protein Kinase Inhibitors/pharmacology , Protein Kinase Inhibitors/therapeutic use , Protein-Tyrosine Kinases/metabolism , Pyridines , Pyrimidines , Quinolines , Sulfonamides , Sunitinib/therapeutic use , src-Family KinasesABSTRACT
Single-strand selective uracil-DNA glycosylase 1 (SMUG1) initiates base excision repair (BER) of uracil and oxidized pyrimidines. SMUG1 status has been associated with cancer risk and therapeutic response in breast carcinomas and other cancer types. However, SMUG1 is a multifunctional protein involved, not only, in BER but also in RNA quality control, and its function in cancer cells is unclear. Here we identify several novel SMUG1 interaction partners that functions in many biological processes relevant for cancer development and treatment response. Based on this, we hypothesized that the dominating function of SMUG1 in cancer might be ascribed to functions other than BER. We define a bad prognosis signature for SMUG1 by mapping out the SMUG1 interaction network and found that high expression of genes in the bad prognosis network correlated with lower survival probability in ER+ breast cancer. Interestingly, we identified hsa-let-7b-5p microRNA as an upstream regulator of the SMUG1 interactome. Expression of SMUG1 and hsa-let-7b-5p were negatively correlated in breast cancer and we found an inhibitory auto-regulatory loop between SMUG1 and hsa-let-7b-5p in the MCF7 breast cancer cells. We conclude that SMUG1 functions in a gene regulatory network that influence the survival and treatment response in several cancers.
Subject(s)
Breast Neoplasms , MicroRNAs , Breast Neoplasms/genetics , Breast Neoplasms/pathology , Female , Humans , MicroRNAs/genetics , Prognosis , Uracil/metabolism , Uracil-DNA Glycosidase/geneticsABSTRACT
Long non-coding RNAs (lncRNAs) are involved in breast cancer pathogenesis through chromatin remodeling, transcriptional and post-transcriptional gene regulation. We report robust associations between lncRNA expression and breast cancer clinicopathological features in two population-based cohorts: SCAN-B and TCGA. Using co-expression analysis of lncRNAs with protein coding genes, we discovered three distinct clusters of lncRNAs. In silico cell type deconvolution coupled with single-cell RNA-seq analyses revealed that these three clusters were driven by cell type specific expression of lncRNAs. In one cluster lncRNAs were expressed by cancer cells and were mostly associated with the estrogen signaling pathways. In the two other clusters, lncRNAs were expressed either by immune cells or fibroblasts of the tumor microenvironment. To further investigate the cis-regulatory regions driving lncRNA expression in breast cancer, we identified subtype-specific transcription factor (TF) occupancy at lncRNA promoters. We also integrated lncRNA expression with DNA methylation data to identify long-range regulatory regions for lncRNA which were validated using ChiA-Pet-Pol2 loops. lncRNAs play an important role in shaping the gene regulatory landscape in breast cancer. We provide a detailed subtype and cell type-specific expression of lncRNA, which improves the understanding of underlying transcriptional regulation in breast cancer.
Subject(s)
Breast Neoplasms , RNA, Long Noncoding , Breast Neoplasms/pathology , DNA Methylation , Female , Gene Expression Regulation , Humans , RNA, Long Noncoding/genetics , RNA, Long Noncoding/metabolism , Tumor MicroenvironmentABSTRACT
BACKGROUND: Locally advanced breast cancer is a heterogeneous disease with respect to response to neoadjuvant chemotherapy (NACT) and survival. It is currently not possible to accurately predict who will benefit from the specific types of NACT. DNA methylation is an epigenetic mechanism known to play an important role in regulating gene expression and may serve as a biomarker for treatment response and survival. We investigated the potential role of DNA methylation as a prognostic marker for long-term survival (> 5 years) after NACT in breast cancer. METHODS: DNA methylation profiles of pre-treatment (n = 55) and post-treatment (n = 75) biopsies from 83 women with locally advanced breast cancer were investigated using the Illumina HumanMethylation450 BeadChip. The patients received neoadjuvant treatment with epirubicin and/or paclitaxel. Linear mixed models were used to associate DNA methylation to treatment response and survival based on clinical response to NACT (partial response or stable disease) and 5-year survival, respectively. LASSO regression was performed to identify a risk score based on the statistically significant methylation sites and Kaplan-Meier curve analysis was used to estimate survival probabilities using ten years of survival follow-up data. The risk score developed in our discovery cohort was validated in an independent validation cohort consisting of paired pre-treatment and post-treatment biopsies from 85 women with locally advanced breast cancer. Patients included in the validation cohort were treated with either doxorubicin or 5-FU and mitomycin NACT. RESULTS: DNA methylation patterns changed from before to after NACT in 5-year survivors, while no significant changes were observed in non-survivors or related to treatment response. DNA methylation changes included an overall loss of methylation at CpG islands and gain of methylation in non-CpG islands, and these changes affected genes linked to transcription factor activity, cell adhesion and immune functions. A risk score was developed based on four methylation sites which successfully predicted long-term survival in our cohort (p = 0.0034) and in an independent validation cohort (p = 0.049). CONCLUSION: Our results demonstrate that DNA methylation patterns in breast tumors change in response to NACT. These changes in DNA methylation show potential as prognostic biomarkers for breast cancer survival.