Integrated Chromatin Accessibility and Transcriptome Landscapes of Doxorubicin-Resistant Breast Cancer Cells

Analysis of chromatin accessibility in peripheral blood mononuclear cells from patients with early-stage breast cancer

Objective: This study was aimed at exploring a specific open region of chromatin in the peripheral blood mononuclear cells (PBMCs) of patients with breast cancer and evaluating its feasibility as a biomarker for diagnosing and predicting breast cancer prognosis. Methods: We obtained PBMCs from breast cancer patients and healthy people for the assay for transposase-accessible chromatin (ATAC) sequencing (n = 3) and obtained the GSE27562 chip sequencing data for secondary analyses. Through bioinformatics analysis, we mined the pattern changes for chromatin accessibility in the PBMCs of breast cancer patients. Results: A total of 1,906 differentially accessible regions (DARs) and 1,632 differentially expressed genes (DEGs) were identified via ATAC sequencing. The upregulated DEGs in the disease group were mainly distributed in the cells, organelles, and cell-intima-related structures and were mainly responsible for biological functions such as cell nitrogen complex metabolism, macromolecular metabolism, and cell communication, in addition to functions such as nucleic acid binding, enzyme binding, hydrolase reaction, and transferase activity. Combined with microarray data analysis, the following set of nine DEGs showed intersection between the ATAC and microarray data: JUN, MSL2, CDC42, TRIB1, SERTAD3, RAB14, RHOB, RAB40B, and PRKDC. HOMER predicted and identified five transcription factors that could potentially bind to these peak sites, namely NFY, Sp 2, GFY, NRF, and ELK 1. Conclusion: Chromatin accessibility analysis of the PBMCs from patients with early-stage breast cancer underscores its potential as a significant avenue for biomarker discovery in breast cancer diagnostics and treatment. By screening the transcription factors and DEGs related to breast cancer, this study provides a comprehensive theoretical foundation that is expected to guide future clinical applications and therapeutic developments.

Multi-omic lineage tracing predicts the transcriptional, epigenetic and genetic determinants of cancer evolution

Cancer is a highly heterogeneous disease, where phenotypically distinct subpopulations coexist and can be primed to different fates. Both genetic and epigenetic factors may drive cancer evolution, however little is known about whether and how such a process is pre-encoded in cancer clones. Using single-cell multi-omic lineage tracing and phenotypic assays, we investigate the predictive features of either tumour initiation or drug tolerance within the same cancer population. Clones primed to tumour initiation in vivo display two distinct transcriptional states at baseline. Remarkably, these states share a distinctive DNA accessibility profile, highlighting an epigenetic basis for tumour initiation. The drug tolerant niche is also largely pre-encoded, but only partially overlaps the tumour-initiating one and evolves following two genetically and transcriptionally distinct trajectories. Our study highlights coexisting genetic, epigenetic and transcriptional determinants of cancer evolution, unravelling the molecular complexity of pre-encoded tumour phenotypes.

Drug repositioning for immunotherapy in breast cancer using single-cell analysis

Immunomodulatory peptides, while exhibiting potential antimicrobial, antifungal, and/or antiviral properties, can play a role in stimulating or suppressing the immune system, especially in pathological conditions like breast cancer (BC). Thus, deregulation of these peptides may serve as an immunotherapeutic strategy to enhance the immune response. In this meta-analysis, we utilized single-cell RNA sequencing data and known therapeutic peptides to investigate the deregulation of these peptides in malignant versus normal human breast epithelial cells. We corroborated our findings at the chromatin level using ATAC-seq. Additionally, we assessed the protein levels in various BC cell lines. Moreover, our in-house drug repositioning approach was employed to identify potential drugs that could positively impact the relapse-free survival of BC patients. Considering significantly deregulated therapeutic peptides and their role in BC pathology, our approach aims to downregulate B2M and SLPI, while upregulating PIGR, DEFB1, LTF, CLU, S100A7, and SCGB2A1 in BC epithelial cells through our drug repositioning pipeline. Leveraging the LINCS L1000 database, we propose BRD-A06641369 for B2M downregulation and ST-4070043 and BRD-K97926541 for SLPI downregulation without negatively affecting the MHC complex as a significantly correlated pathway with these two genes. Furthermore, we have compiled a comprehensive list of drugs for the upregulation of other selected immunomodulatory peptides. Employing an immunotherapeutic approach by integrating our drug repositioning pipeline with single-cell analysis, we proposed potential drugs and drug targets to fortify the immune system against BC.

Integrative analysis of chromatin accessibility and transcriptome landscapes in the induction of peritoneal fibrosis by high glucose

BACKGROUND

Peritoneal fibrosis is the prevailing complication induced by prolonged exposure to high glucose in patients undergoing peritoneal dialysis.

METHODS

To elucidate the molecular mechanisms underlying this process, we conducted an integrated analysis of the transcriptome and chromatin accessibility profiles of human peritoneal mesothelial cells (HMrSV5) during high-glucose treatment.

RESULTS

Our study identified 2775 differentially expressed genes (DEGs) related to high glucose-triggered pathological changes, including 1164 upregulated and 1611 downregulated genes. Genome-wide DEGs and network analysis revealed enrichment in the epithelial-mesenchymal transition (EMT), inflammatory response, hypoxia, and TGF-beta pathways. The enriched genes included VEGFA, HIF-1α, TGF-β1, EGF, TWIST2, and SNAI2. Using ATAC-seq, we identified 942 hyper (higher ATAC-seq signal in high glucose-treated HMrSV5 cells than in control cells) and 714 hypo (lower ATAC-seq signal in high glucose-treated HMrSV5 cells versus control cells) peaks with differential accessibility in high glucose-treated HMrSV5 cells versus controls. These differentially accessible regions were positively correlated (R = 0.934) with the nearest DEGs. These genes were associated with 566 up- and 398 downregulated genes, including SNAI2, TGF-β1, HIF-1α, FGF2, VEGFA, and VEGFC, which are involved in critical pathways identified by transcriptome analysis. Integrated ATAC-seq and RNA-seq analysis also revealed key transcription factors (TFs), such as HIF-1α, ARNTL, ELF1, SMAD3 and XBP1. Importantly, we demonstrated that HIF-1α is involved in the regulation of several key genes associated with EMT and the TGF-beta pathway. Notably, we predicted and experimentally validated that HIF-1α can exacerbate the expression of TGF-β1 in a high glucose-dependent manner, revealing a novel role of HIF-1α in high glucose-induced pathological changes in human peritoneal mesothelial cells (HPMCs).

CONCLUSIONS

In summary, our study provides a comprehensive view of the role of transcriptome deregulation and chromosome accessibility alterations in high glucose-induced pathological fibrotic changes in HPMCs. This analysis identified hub genes, signaling pathways, and key transcription factors involved in peritoneal fibrosis and highlighted the novel glucose-dependent regulation of TGF-β1 by HIF-1α. This integrated approach has offered a deeper understanding of the pathogenesis of peritoneal fibrosis and has indicated potential therapeutic targets for intervention.

BHLHE40 Inhibits Ferroptosis in Pancreatic Cancer Cells via Upregulating SREBF1

Pancreatic cancer (PCa) is one of the most fatal human malignancies. The enhanced infiltration of stromal tissue into the PCa tumor microenvironment limits the identification of key tumor-specific transcription factors and epigenomic abnormalities in malignant epithelial cells. Integrated transcriptome and epigenetic multiomics analyses of the paired PCa organoids indicate that the basic helix-loop-helix transcription factor 40 (BHLHE40) is significantly upregulated in tumor samples. Increased chromatin accessibility at the promoter region and enhanced mTOR pathway activity contribute to the elevated expression of BHLHE40. Integrated analysis of chromatin immunoprecipitation-seq, RNA-seq, and high-throughput chromosome conformation capture data, together with chromosome conformation capture assays, indicate that BHLHE40 not only regulates sterol regulatory element-binding factor 1 (SREBF1) transcription as a classic transcription factor but also links the enhancer and promoter regions of SREBF1. It is found that the BHLHE40-SREBF1-stearoyl-CoA desaturase axis protects PCa cells from ferroptosis, resulting in the reduced accumulation of lipid peroxidation. Moreover, fatostatin, an SREBF1 inhibitor, significantly suppresses the growth of PCa tumors with high expressions of BHLHE40. This study highlights the important roles of BHLHE40-mediated lipid peroxidation in inducing ferroptosis in PCa cells and provides a novel mechanism underlying SREBF1 overexpression in PCa.

An integrated multi-omics analysis of identifies distinct molecular characteristics in pulmonary infections of Pseudomonas aeruginosa

Pseudomonas aeruginosa (P. aeruginosa) can cause severe acute infections, including pneumonia and sepsis, and cause chronic infections, commonly in patients with structural respiratory diseases. However, the molecular and pathophysiological mechanisms of P. aeruginosa respiratory infection are largely unknown. Here, we performed assays for transposase-accessible chromatin using sequencing (ATAC-seq), transcriptomics, and quantitative mass spectrometry-based proteomics and ubiquitin-proteomics in P. aeruginosa-infected lung tissues for multi-omics analysis, while ATAC-seq and transcriptomics were also examined in P. aeruginosa-infected mouse macrophages. To identify the pivotal factors that are involved in host immune defense, we integrated chromatin accessibility and gene expression to investigate molecular changes in P. aeruginosa-infected lung tissues combined with proteomics and ubiquitin-proteomics. Our multi-omics investigation discovered a significant concordance for innate immunological and inflammatory responses following P. aeruginosa infection between hosts and alveolar macrophages. Furthermore, we discovered that multi-omics changes in pioneer factors Stat1 and Stat3 play a crucial role in the immunological regulation of P. aeruginosa infection and that their downstream molecules (e.g., Fas) may be implicated in both immunosuppressive and inflammation-promoting processes. Taken together, these findings indicate that transcription factors and their downstream signaling molecules play a critical role in the mobilization and rebalancing of the host immune response against P. aeruginosa infection and may serve as potential targets for bacterial infections and inflammatory diseases, providing insights and resources for omics analyses.

Multi-omic analysis reveals dynamic changes of three-dimensional chromatin architecture during T cell differentiation

Cell differentiation results in widespread changes in transcriptional programs as well as multi-level remodeling of three-dimensional genome architecture. Nonetheless, few synthetically investigate the chromatin higher-order landscapes in different T helper (Th) cells. Using RNA-Seq, ATAC-Seq and Hi-C assays, we characterize dynamic changes in chromatin organization at different levels during Naive CD4⁺ T cells differentiation into T helper 17 (Th17) and T helper 1 (Th1) cells. Upon differentiation, we observe decreased short-range and increased extra-long-range chromatin interactions. Although there is no apparent global switch in the A/B compartments, Th cells display the weaker compartmentalization. A portion of topologically associated domains are rearranged. Furthermore, we identify cell-type specific enhancer-promoter loops, many of which are associated with functional genes in Th cells, such as Rorc facilitating Th17 differentiation and Hif1a responding to intracellular oxygen levels in Th1. Taken together, these results uncover the general patterns of chromatin reorganization and epigenetic landscapes of gene regulation during T helper cell differentiation.

Inhibition of VEGF binding to neuropilin-2 enhances chemosensitivity and inhibits metastasis in triple-negative breast cancer

Although blocking the binding of vascular endothelial growth factor (VEGF) to neuropilin-2 (NRP2) on tumor cells is a potential strategy to treat aggressive carcinomas, a lack of effective reagents that can be used clinically has hampered this potential therapy. Here, we describe the generation of a fully humanized, high-affinity monoclonal antibody (aNRP2-10) that specifically inhibits the binding of VEGF to NRP2, conferring antitumor activity without causing toxicity. Using triple-negative breast cancer as a model, we demonstrated that aNRP2-10 could be used to isolate cancer stem cells (CSCs) from heterogeneous tumor populations and inhibit CSC function and epithelial-to-mesenchymal transition. aNRP2-10 sensitized cell lines, organoids, and xenografts to chemotherapy and inhibited metastasis by promoting the differentiation of CSCs to a state that is more responsive to chemotherapy and less prone to metastasis. These data provide justification for the initiation of clinical trials designed to improve the response of patients with aggressive tumors to chemotherapy using this monoclonal antibody.

Sustained activation of non-canonical NF-κB signalling drives glycolytic reprogramming in doxorubicin-resistant DLBCL

DLBCL is the most common lymphoma with high tumor heterogeneity. Treatment refractoriness and relapse from R-CHOP therapy in patients remain a clinical problem. Activation of the non-canonical NF-κB pathway is associated with R-CHOP resistance. However, downstream targets of non-canonical NF-κB mediating R-CHOP-induced resistance remains uncharacterized. Here, we identify the common mechanisms underlying both intrinsic and acquired resistance that are induced by doxorubicin, the main cytotoxic component of R-CHOP. We performed global transcriptomic analysis of (1) a panel of resistant versus sensitive and (2) isogenic acquired doxorubicin-resistant DLBCL cell lines following short and chronic exposure to doxorubicin respectively. Doxorubicin-induced stress in resistant cells activates a distinct transcriptional signature that is enriched in metabolic reprogramming and oncogenic signalling. Selective and sustained activation of non-canonical NF-κB signalling in these resistant cells exacerbated their survival by augmenting glycolysis. In response to doxorubicin, p52-RelB complexes transcriptionally activated multiple glycolytic regulators with prognostic significance through increased recruitment at their gene promoters. Targeting p52-RelB and their targets in resistant cells increased doxorubicin sensitivity in vitro and in vivo. Collectively, our study uncovered novel molecular drivers of doxorubicin-induced resistance that are regulated by non-canonical NF-κB pathway. We reveal new avenues of therapeutic targeting for R-CHOP-treated refractory/relapsed DLBCL patients.

When cancer drug resistance meets metabolomics (bulk, single-cell and/or spatial): Progress, potential, and perspective

Resistance to drug treatment is a critical barrier in cancer therapy. There is an unmet need to explore cancer hallmarks that can be targeted to overcome this resistance for therapeutic gain. Over time, metabolic reprogramming has been recognised as one hallmark that can be used to prevent therapeutic resistance. With the advent of metabolomics, targeting metabolic alterations in cancer cells and host patients represents an emerging therapeutic strategy for overcoming cancer drug resistance. Driven by technological and methodological advances in mass spectrometry imaging, spatial metabolomics involves the profiling of all the metabolites (metabolomics) so that the spatial information is captured bona fide within the sample. Spatial metabolomics offers an opportunity to demonstrate the drug-resistant tumor profile with metabolic heterogeneity, and also poses a data-mining challenge to reveal meaningful insights from high-dimensional spatial information. In this review, we discuss the latest progress, with the focus on currently available bulk, single-cell and spatial metabolomics technologies and their successful applications in pre-clinical and translational studies on cancer drug resistance. We provide a summary of metabolic mechanisms underlying cancer drug resistance from different aspects; these include the Warburg effect, altered amino acid/lipid/drug metabolism, generation of drug-resistant cancer stem cells, and immunosuppressive metabolism. Furthermore, we propose solutions describing how to overcome cancer drug resistance; these include early detection during cancer initiation, monitoring of clinical drug response, novel anticancer drug and target metabolism, immunotherapy, and the emergence of spatial metabolomics. We conclude by describing the perspectives on how spatial omics approaches (integrating spatial metabolomics) could be further developed to improve the management of drug resistance in cancer patients.

Chromatin Variants Reveal the Genetic Determinants of Oncogenesis in Breast Cancer

Breast cancer presents as multiple distinct disease entities. Each tumor harbors diverse cell populations defining a phenotypic heterogeneity that impinges on our ability to treat patients. To date, efforts mainly focused on genetic variants to find drivers of inter- and intratumor phenotypic heterogeneity. However, these efforts have failed to fully capture the genetic basis of breast cancer. Through recent technological and analytical approaches, the genetic basis of phenotypes can now be decoded by characterizing chromatin variants. These variants correspond to polymorphisms in chromatin states at DNA sequences that serve a distinct role across cell populations. Here, we review the function and causes of chromatin variants as they relate to breast cancer inter- and intratumor heterogeneity and how they can guide the development of treatment alternatives to fulfill the goal of precision cancer medicine.

Reorganization of 3D chromatin architecture in doxorubicin-resistant breast cancer cells

Background: Doxorubicin resistance remains a major therapeutic challenge leading to poor survival prognosis and treatment failure in breast cancer. Although doxorubicin induces massive changes in the transcriptional landscape are well known, potential diagnostic or therapeutic targets associated with the reorganization of three-dimensional (3D) chromatin architecture have not yet been systematically investigated.

Methods: Here we performed in situ high-throughput chromosome conformation capture (Hi-C) on parental and doxorubicin-resistant MCF7 (MCF7-DR) human breast cancer cells, followed by integrative analysis of HiC, ATAC-seq, RNA-seq and TCGA data.

Results: It revealed that A/B compartment switching was positively correlated to genome-wide differential gene expression. The genome of MCF7-DR cells was spatially reorganized into smaller topologically associating domains (TADs) and chromatin loops. We also revealed the contribution of increased chromatin accessibility and potential transcription factor families, including CTCF, AP-1 and bHLH, to gained TADs or loops. Intriguingly, we observed two condensed genomic regions (∼20 kb) with decreased chromatin accessibility flanking TAD boundaries, which might play a critical role in the formation or maintenance of TADs. Finally, combining data from TCGA, we identified a number of gained and lost enhancer-promoter interactions and their corresponding differentially expressed genes involved in chromatin organization and breast cancer signaling pathways, including FA2H, FOXA1 and JRKL, which might serve as potential treatment targets for breast cancer.

Conclusion: These data uncovered a close connection between 3D genome reorganization, chromatin accessibility as well as gene transcription and provide novel insights into the epigenomic mechanisms involving doxorubicin resistance in breast cancer.

Changes in the Transcriptome and Chromatin Landscape in BRAFi-Resistant Melanoma Cells

Metastatic and drug-resistant melanoma are leading causes of skin cancer-associated death. Mitogen-associated protein kinase (MAPK) pathway inhibitors and immunotherapies have provided substantial benefits to patients with melanoma. However, long-term therapeutic efficacy has been limited due to emergence of treatment resistance. Despite the identification of several molecular mechanisms underlying the development of resistant phenotypes, significant progress has still not been made toward the effective treatment of drug-resistant melanoma. Therefore, the identification of new targets and mechanisms driving drug resistance in melanoma represents an unmet medical need. In this study, we performed unbiased RNA-sequencing (RNA-seq) and assay for transposase-accessible chromatin with sequencing (ATAC-seq) to identify new targets and mechanisms that drive resistance to MAPK pathway inhibitors targeting BRAF and MAPK kinase (MEK) in BRAF-mutant melanoma cells. An integrative analysis of ATAC-seq combined with RNA-seq showed that global changes in chromatin accessibility affected the mRNA expression levels of several known and novel genes, which consequently modulated multiple oncogenic signaling pathways to promote resistance to MAPK pathway inhibitors in melanoma cells. Many of these genes were also associated with prognosis predictions in melanoma patients. This study resulted in the identification of new genes and signaling pathways that might be targeted to treat MEK or BRAF inhibitors resistant melanoma patients. The present study applied new and advanced approaches to identify unique changes in chromatin accessibility regions that modulate gene expression associated with pathways to promote the development of resistance to MAPK pathway inhibitors.

Novel Chemicals Derived from Tadalafil Exhibit PRMT5 Inhibition and Promising Activities against Breast Cancer

Breast cancer seriously endangers women’s health worldwide. Protein arginine methyltransferase 5 (PRMT5) is highly expressed in breast cancer and represents a potential druggable target for breast cancer treatment. However, because the currently available clinical PRMT5 inhibitors are relatively limited, there is an urgent need to develop new PRMT5 inhibitors. Our team previously found that the FDA-approved drug tadalafil can act as a PRMT5 inhibitor and enhance the sensitivity of breast cancer patients to doxorubicin treatment. To further improve the binding specificity of tadalafil to PRMT5, we chemically modified tadalafil, and designed three compounds, A, B, and C, based on the PRMT5 protein structure. These three compounds could bind to PRMT5 through different binding modes and inhibit histone arginine methylation. They arrested the proliferation and triggered the apoptosis of breast cancer cells in vitro and also promoted the antitumor effects of the chemotherapy drugs cisplatin, doxorubicin, and olaparib in combination regimens. Among them, compound A possessed the highest potency. Finally, the anti-breast cancer effects of PRMT5 inhibitor A and its ability to enhance chemosensitivity were further verified in a xenograft mouse model. These results indicate that the new PRMT5 inhibitors A, B, and C may be potential candidates for breast cancer treatment.

Discovery of Proteins Responsible for Resistance to Three Chemotherapy Drugs in Breast Cancer Cells Using Proteomics and Bioinformatics Analysis

Chemoresistance is a daunting obstacle to the effective treatment of breast cancer patients receiving chemotherapy. Although the mechanism of chemotherapy drug resistance has been explored broadly, the precise mechanism at the proteome level remains unclear. Especially, comparative studies between widely used anticancer drugs in breast cancer are very limited. In this study, we employed proteomics and bioinformatics approaches on chemoresistant breast cancer cell lines to understand the underlying resistance mechanisms that resulted from doxorubicin (DR), paclitaxel (PR), and tamoxifen (TAR). In total, 10,385 proteins were identified and quantified from three TMT 6-plex and one TMT 10-plex experiments. Bioinformatics analysis showed that Notch signaling, immune response, and protein re-localization processes were uniquely associated with DR, PR, and TAR resistance, respectively. In addition, proteomic signatures related to drug resistance were identified as potential targets of many FDA-approved drugs. Furthermore, we identified potential prognostic proteins with significant effects on overall survival. Representatively, PLXNB2 expression was associated with a highly significant increase in risk, and downregulation of ACOX3 was correlated with a worse overall survival rate. Consequently, our study provides new insights into the proteomic aspects of the distinct mechanisms underlying chemoresistance in breast cancer.