eRNA profiling uncovers the enhancer landscape of oesophageal adenocarcinoma and reveals new deregulated pathways
Abstract
Cancer is driven by both genetic and epigenetic changes that impact on gene expression profiles and the resulting tumourigenic phenotype. Enhancers are transcriptional regulatory elements that are key to our understanding of how this rewiring of gene expression is achieved in cancer cells. Here we have harnessed the power of RNA-seq data from hundreds of patients with oesophageal adenocarcinoma (OAC) or its precursor state Barrett's oesophagus (BO) coupled with open chromatin maps to identify potential enhancer RNAs (eRNAs) and their associated enhancer regions in this cancer. We identify ~1000 OAC-specific enhancers and use this data to uncover new cellular pathways that are operational in OAC. Among these are enhancers for JUP, MYBL2 and CCNE1, and we show that their activity is required for cancer cell viability. We also demonstrate the clinical utility of our dataset for identifying disease stage and patient prognosis. Our data therefore identify an important set of regulatory elements that enhance our molecular understanding of OAC and point to potential new therapeutic directions.
Data availability
All data have been deposited at ArrayExpress; OE19 KAS-seq and CUT&TAG data (E-MTAB-11357 and E-MTAB-11356, respectively) and OE19 HiC data (E-MTAB-12664).
-
KAS-seq in OE19 cellsEBI ArrayExpress, E-MTAB-11357.
-
CUT&TAG of OE19 cell lineEBI ArrayExpress, E-MTAB-11356.
-
HiC data in OE19 cellsEBI ArrayExpress, E-MTAB-12664.
-
Sequencing data for oesophageal and related samples - Ogden et al releaseEuropean Genome Phenome Archive, EGAD00001007496.
-
ATAC-seq of oesophageal cell lines and tissue samplesEBI ArrayExpress, E-MTAB-5169.
-
ATAC-seq of human Barrett's oesophagus tissueEBI ArrayExpress, E-MTAB-6751.
-
ATAC-seq of oesophageal adenocarcinoma patient samplesEBI ArrayExpress, E-MTAB-8447.
-
ChIP-seq of H3K27Ac in oesophageal adenocarcinoma OE19 cellsEBI ArrayExpress, E-MTAB-10319.
-
ATAC-seq of oesophageal adenocarcinoma patient samplesEBI ArrayExpress, E-MTAB-8447.
-
RNA-seq of OE19 cells treated with siNT or siKLF5 for 72 hoursEBI ArrayExpress, E-MTAB-8568.
-
KLF5 ChIP-seq in CP-A and OE19 cellsEBI ArrayExpress, E-MTAB-8568.
-
OE19 dnFOS RNA-seqEBI ArrayExpress, E-MTAB-10334.
Article and author information
Author details
Funding
Medical Research Council (MR/V010263/1)
- Samuel Ogden
Wellcome Trust (102171/B/13/Z)
- Ibrahim Ahmed
The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.
Copyright
© 2023, Ahmed et al.
This article is distributed under the terms of the Creative Commons Attribution License permitting unrestricted use and redistribution provided that the original author and source are credited.
Metrics
-
- 2,140
- views
-
- 264
- downloads
-
- 8
- citations
Views, downloads and citations are aggregated across all versions of this paper published by eLife.
Download links
Downloads (link to download the article as PDF)
Open citations (links to open the citations from this article in various online reference manager services)
Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)
Further reading
-
- Cancer Biology
- Computational and Systems Biology
Effects from aging in single cells are heterogenous, whereas at the organ- and tissue-levels aging phenotypes tend to appear as stereotypical changes. The mammary epithelium is a bilayer of two major phenotypically and functionally distinct cell lineages: luminal epithelial and myoepithelial cells. Mammary luminal epithelia exhibit substantial stereotypical changes with age that merit attention because these cells are the putative cells-of-origin for breast cancers. We hypothesize that effects from aging that impinge upon maintenance of lineage fidelity increase susceptibility to cancer initiation. We generated and analyzed transcriptomes from primary luminal epithelial and myoepithelial cells from younger <30 (y)ears old and older >55 y women. In addition to age-dependent directional changes in gene expression, we observed increased transcriptional variance with age that contributed to genome-wide loss of lineage fidelity. Age-dependent variant responses were common to both lineages, whereas directional changes were almost exclusively detected in luminal epithelia and involved altered regulation of chromatin and genome organizers such as SATB1. Epithelial expression variance of gap junction protein GJB6 increased with age, and modulation of GJB6 expression in heterochronous co-cultures revealed that it provided a communication conduit from myoepithelial cells that drove directional change in luminal cells. Age-dependent luminal transcriptomes comprised a prominent signal that could be detected in bulk tissue during aging and transition into cancers. A machine learning classifier based on luminal-specific aging distinguished normal from cancer tissue and was highly predictive of breast cancer subtype. We speculate that luminal epithelia are the ultimate site of integration of the variant responses to aging in their surrounding tissue, and that their emergent phenotype both endows cells with the ability to become cancer-cells-of-origin and represents a biosensor that presages cancer susceptibility.
-
- Cancer Biology
Chemotherapy is widely used to treat lung adenocarcinoma (LUAD) patients comprehensively. Considering the limitations of chemotherapy due to drug resistance and other issues, it is crucial to explore the impact of chemotherapy and immunotherapy on these aspects. In this study, tumor samples from nine LUAD patients, of which four only received surgery and five received neoadjuvant chemotherapy, were subjected to scRNA-seq analysis. In vitro and in vivo assays, including flow cytometry, immunofluorescence, Seahorse assay, and tumor xenograft models, were carried out to validate our findings. A total of 83,622 cells were enrolled for subsequent analyses. The composition of cell types exhibited high heterogeneity across different groups. Functional enrichment analysis revealed that chemotherapy drove significant metabolic reprogramming in tumor cells and macrophages. We identified two subtypes of macrophages: Anti-mac cells (CD45+CD11b+CD86+) and Pro-mac cells (CD45+CD11b+ARG +) and sorted them by flow cytometry. The proportion of Pro-mac cells in LUAD tissues increased significantly after neoadjuvant chemotherapy. Pro-mac cells promote tumor growth and angiogenesis and also suppress tumor immunity. Moreover, by analyzing the remodeling of T and B cells induced by neoadjuvant therapy, we noted that chemotherapy ignited a relatively more robust immune cytotoxic response toward tumor cells. Our study demonstrates that chemotherapy induces metabolic reprogramming within the tumor microenvironment of LUAD, particularly affecting the function and composition of immune cells such as macrophages and T cells. We believe our findings will offer insight into the mechanisms of drug resistance and provide novel therapeutic targets for LUAD in the future.