Overexpressed HSF1 cancer signature genes cluster in human chromosome 8q

Low-Pass Genomic Sequencing Reveals Novel Subtypes of Pancreatic Cystic Neoplasms

BACKGROUND

Over the years, the detection rate of pancreatic cystic neoplasms (PCNs) has significantly increased; however, the differential diagnosis and identification of high-risk PCNs remain challenging. We sought to investigate whether chromosomal instability (CIN) features in cell-free DNA in the cystic fluid of PCNs could help to identify high-risk PCNs.

METHODS

Pancreatic cystic fluid samples from 102 patients with PCNs were intraoperatively collected for detection of CIN using an ultrasensitive chromosomal aneuploidy detector. Clinical and imaging data were retrospectively collected, and statistical analysis was performed to assess the potential role of CIN in clinical practice.

RESULTS

CIN was investigated in a total of 100 patients. Sixteen of 26 serous cystic cystadenomas (SCAs) harbored deletions of chr3p and/or chr6p, whereas low rates of CIN were detected in mucinous cystic neoplasms. Most malignant PCNs presented with more than one type of CIN; amplification of chr1q and chr8q found in nine and seven of 11 malignant PCNs (81.8% and 63.6%), respectively, could aid in distinguishing high-risk IPMNs from low-risk ones, with a higher sensitivity than imaging. A combination of the mural nodule imaging feature and amplification of chr1q and chr8q achieved a sensitivity of 70.0% and a specificity of 82.4% in identifying high-risk IPMNs.

CONCLUSIONS

Our work revealed the distinct CIN signature of different types of PCNs. Deletions of chr3p and chr6p defined a subtype of SCAs. Gains of chr1q and chr8q were associated with insidious malignant PCNs and helped identify high-risk IPMNs.

Prognostic Values of Gene Copy Number Alterations in Prostate Cancer

Whilst risk prediction for individual prostate cancer (PCa) cases is of a high priority, the current risk stratification indices for PCa management have severe limitations. This study aimed to identify gene copy number alterations (CNAs) with prognostic values and to determine if any combination of gene CNAs could have risk stratification potentials. Clinical and genomic data of 500 PCa cases from the Cancer Genome Atlas stable were retrieved from the Genomic Data Commons and cBioPortal databases. The CNA statuses of a total of 52 genetic markers, including 21 novel markers and 31 previously identified potential prognostic markers, were tested for prognostic significance. The CNA statuses of a total of 51/52 genetic markers were significantly associated with advanced disease at an odds ratio threshold of ≥1.5 or ≤0.667. Moreover, a Kaplan-Meier test identified 27/52 marker CNAs which correlated with disease progression. A Cox Regression analysis showed that the amplification of MIR602 and deletions of MIR602, ZNF267, MROH1, PARP8, and HCN1 correlated with a progression-free survival independent of the disease stage and Gleason prognostic group grade. Furthermore, a binary logistic regression analysis identified twenty-two panels of markers with risk stratification potentials. The best model of 7/52 genetic CNAs, which included the SPOP alteration, SPP1 alteration, CCND1 amplification, PTEN deletion, CDKN1B deletion, PARP8 deletion, and NKX3.1 deletion, stratified the PCa cases into a localised and advanced disease with an accuracy of 70.0%, sensitivity of 85.4%, specificity of 44.9%, positive predictive value of 71.67%, and negative predictive value of 65.35%. This study validated prognostic gene level CNAs identified in previous studies, as well as identified new genetic markers with CNAs that could potentially impact risk stratification in PCa.

The Revelation of Continuously Organized, Co-Overexpressed Protein-Coding Genes with Roles in Cellular Communications in Breast Cancer

Many human cancers, including breast cancer, are polygenic and involve the co-dysregulation of multiple regulatory molecules and pathways. Though the overexpression of genes and amplified chromosomal regions have been closely linked in breast cancer, the notion of the co-upregulation of genes at a single locus remains poorly described. Here, we describe the co-overexpression of 34 continuously organized protein-coding genes with diverse functions at 8q.24.3(143437655-144326919) in breast and other cancer types, the CanCord34 genes. In total, 10 out of 34 genes have not been reported to be overexpressed in breast cancer. Interestingly, the overexpression of CanCord34 genes is not necessarily associated with genomic amplification and is independent of hormonal or HER2 status in breast cancer. CanCord34 genes exhibit diverse known and predicted functions, including enzymatic activities, cell viability, multipotency, cancer stem cells, and secretory activities, including extracellular vesicles. The co-overexpression of 33 of the CanCord34 genes in a multivariant analysis was correlated with poor survival among patients with breast cancer. The analysis of the genome-wide RNAi functional screening, cell dependency fitness, and breast cancer stem cell databases indicated that three diverse overexpressed CanCord34 genes, including a component of spliceosome PUF60, a component of exosome complex EXOSC4, and a ribosomal biogenesis factor BOP1, shared roles in cell viability, cell fitness, and stem cell phenotypes. In addition, 17 of the CanCord34 genes were found in the microvesicles (MVs) secreted from the mesenchymal stem cells that were primed with MDA-MB-231 breast cancer cells. Since these MVs were important in the chemoresistance and dedifferentiation of breast cancer cells into cancer stem cells, these findings highlight the significance of the CanCord34 genes in cellular communications. In brief, the persistent co-overexpression of CanCord34 genes with diverse functions can lead to the dysregulation of complementary functions in breast cancer. In brief, the present study provides new insights into the polygenic nature of breast cancer and opens new research avenues for basic, preclinical, and therapeutic studies in human cancer.

NSMCE2, a novel super-enhancer-regulated gene, is linked to poor prognosis and therapy resistance in breast cancer

BACKGROUND

Despite today's advances in the treatment of cancer, breast cancer-related mortality remains high, in part due to the lack of effective targeted therapies against breast tumor types that do not respond to standard treatments. Therefore, identifying additional breast cancer molecular targets is urgently needed. Super-enhancers are large regions of open chromatin involved in the overactivation of oncogenes. Thus, inhibition of super-enhancers has become a focus in clinical trials for its therapeutic potential. Here, we aimed to identify novel super-enhancer dysregulated genes highly associated with breast cancer patients' poor prognosis and negative response to treatment.

METHODS

Using existing datasets containing super-enhancer-associated genes identified in breast tumors and public databases comprising genomic and clinical information for breast cancer patients, we investigated whether highly expressed super-enhancer-associated genes correlate to breast cancer patients' poor prognosis and to patients' poor response to therapy. Our computational findings were experimentally confirmed in breast cancer cells by pharmacological SE disruption and gene silencing techniques.

RESULTS

We bioinformatically identified two novel super-enhancer-associated genes - NSMCE2 and MAL2 - highly upregulated in breast tumors, for which high RNA levels significantly and specifically correlate with breast cancer patients' poor prognosis. Through in-vitro pharmacological super-enhancer disruption assays, we confirmed that super-enhancers upregulate NSMCE2 and MAL2 transcriptionally, and, through bioinformatics, we found that high levels of NSMCE2 strongly associate with patients' poor response to chemotherapy, especially for patients diagnosed with aggressive triple negative and HER2 positive tumor types. Finally, we showed that decreasing NSMCE2 gene expression increases breast cancer cells' sensitivity to chemotherapy treatment.

CONCLUSIONS

Our results indicate that moderating the transcript levels of NSMCE2 could improve patients' response to standard chemotherapy consequently improving disease outcome. Our approach offers a new avenue to identify a signature of tumor specific genes that are not frequently mutated but dysregulated by super-enhancers. As a result, this strategy can lead to the discovery of potential and novel pharmacological targets for improving targeted therapy and the treatment of breast cancer.

Copy Number Variation of Circulating Tumor DNA (ctDNA) Detected Using NIPT in Neoadjuvant Chemotherapy-Treated Ovarian Cancer Patients

Analysis of circulating tumor DNA (ctDNA) can be used to characterize and monitor cancers. Recently, non-invasive prenatal testing (NIPT) as a new next-generation sequencing (NGS)-based approach has been applied for detecting ctDNA. This study aimed to investigate the copy number variations (CNVs) utilizing the non-invasive prenatal testing in plasma ctDNA from ovarian cancer (OC) patients who were treated with neoadjuvant chemotherapy (NAC). The plasma samples of six patients, including stages II–IV, were collected during the pre- and post-NAC treatment that were divided into NAC-sensitive and NAC-resistant groups during the follow-up time. CNV analysis was performed using the NIPT via two methods “an open-source algorithm WISECONDORX and NextGENe software.” Results of these methods were compared in pre- and post-NAC of OC patients. Finally, bioinformatics tools were used for data mining from The Cancer Genome Atlas (TCGA) to investigate CNVs in OC patients. WISECONDORX analysis indicated fewer CNV changes on chromosomes before treatment in the NAC-sensitive rather than NAC-resistant patients. NextGENe data indicated that CNVs are not only observed in the coding genes but also in non-coding genes. CNVs in six genes were identified, including HSF1, TMEM249, MROH1, GSTT2B, ABR, and NOMO2, only in NAC-resistant patients. The comparison of these six genes in NAC-resistant patients with The Cancer Genome Atlas data illustrated that the total alteration frequency is amplification, and the highest incidence of the CNVs (≥35% based on TCGA data) is found in MROH1, TMEM249, and HSF1 genes on the chromosome (Chr) 8. Based on TCGA data, survival analysis showed a significant reduction in the overall survival among chemotherapy-resistant patients as well as a high expression level of these three genes compared to that of sensitive samples (all, p < 0.0001). The continued Chr8 study using WISECONDORX revealed CNV modifications in NAC-resistant patients prior to NAC therapy, but no CNV changes were observed in NAC-sensitive individuals. Our findings showed that low coverage whole-genome sequencing analysis used for NIPT could identify CNVs in ctDNA of OC patients before and after chemotherapy. These CNVs are different in NAC-sensitive and -resistant patients highlighting the potential application of this approach in cancer patient management.

Heat Shock Proteins and HSF1 in Cancer

Fitness of cells is dependent on protein homeostasis which is maintained by cooperative activities of protein chaperones and proteolytic machinery. Upon encountering protein-damaging conditions, cells activate the heat-shock response (HSR) which involves HSF1-mediated transcriptional upregulation of a group of chaperones - the heat shock proteins (HSPs). Cancer cells experience high levels of proteotoxic stress due to the production of mutated proteins, aneuploidy-induced excess of components of multiprotein complexes, increased translation rates, and dysregulated metabolism. To cope with this chronic state of proteotoxic stress, cancers almost invariably upregulate major components of HSR, including HSF1 and individual HSPs. Some oncogenic programs show dependence or coupling with a particular HSR factor (such as frequent coamplification of HSF1 and MYC genes). Elevated levels of HSPs and HSF1 are typically associated with drug resistance and poor clinical outcomes in various malignancies. The non-oncogene dependence ("addiction") on protein quality controls represents a pancancer target in treating human malignancies, offering a potential to enhance efficacy of standard and targeted chemotherapy and immune checkpoint inhibitors. In cancers with specific dependencies, HSR components can serve as alternative targets to poorly druggable oncogenic drivers.

Targeting therapy-resistant prostate cancer via a direct inhibitor of the human heat shock transcription factor 1

Heat shock factor 1 (HSF1) is a cellular stress-protective transcription factor exploited by a wide range of cancers to drive proliferation, survival, invasion, and metastasis. Nuclear HSF1 abundance is a prognostic indicator for cancer severity, therapy resistance, and shortened patient survival. The HSF1 gene was amplified, and nuclear HSF1 abundance was markedly increased in prostate cancers and particularly in neuroendocrine prostate cancer (NEPC), for which there are no available treatment options. Despite genetic validation of HSF1 as a therapeutic target in a range of cancers, a direct and selective small-molecule HSF1 inhibitor has not been validated or developed for use in the clinic. We described the identification of a direct HSF1 inhibitor, Direct Targeted HSF1 InhiBitor (DTHIB), which physically engages HSF1 and selectively stimulates degradation of nuclear HSF1. DTHIB robustly inhibited the HSF1 cancer gene signature and prostate cancer cell proliferation. In addition, it potently attenuated tumor progression in four therapy-resistant prostate cancer animal models, including an NEPC model, where it caused profound tumor regression. This study reports the identification and validation of a direct HSF1 inhibitor and provides a path for the development of a small-molecule HSF1-targeted therapy for prostate cancers and other therapy-resistant cancers.

Molecular Mechanisms of Heat Shock Factors in Cancer

Malignant transformation is accompanied by alterations in the key cellular pathways that regulate development, metabolism, proliferation and motility as well as stress resilience. The members of the transcription factor family, called heat shock factors (HSFs), have been shown to play important roles in all of these biological processes, and in the past decade it has become evident that their activities are rewired during tumorigenesis. This review focuses on the expression patterns and functions of HSF1, HSF2, and HSF4 in specific cancer types, highlighting the mechanisms by which the regulatory functions of these transcription factors are modulated. Recently developed therapeutic approaches that target HSFs are also discussed.

HSF1: Primary Factor in Molecular Chaperone Expression and a Major Contributor to Cancer Morbidity

Heat shock factor 1 (HSF1) is the primary component for initiation of the powerful heat shock response (HSR) in eukaryotes. The HSR is an evolutionarily conserved mechanism for responding to proteotoxic stress and involves the rapid expression of heat shock protein (HSP) molecular chaperones that promote cell viability by facilitating proteostasis. HSF1 activity is amplified in many tumor contexts in a manner that resembles a chronic state of stress, characterized by high levels of HSP gene expression as well as HSF1-mediated non-HSP gene regulation. HSF1 and its gene targets are essential for tumorigenesis across several experimental tumor models, and facilitate metastatic and resistant properties within cancer cells. Recent studies have suggested the significant potential of HSF1 as a therapeutic target and have motivated research efforts to understand the mechanisms of HSF1 regulation and develop methods for pharmacological intervention. We review what is currently known regarding the contribution of HSF1 activity to cancer pathology, its regulation and expression across human cancers, and strategies to target HSF1 for cancer therapy.

The Multifaceted Role of HSF1 in Tumorigenesis

Heat Shock Factor 1 (HSF1), the master transcriptional regulator of the heat shock response (HSR), was first cloned more than 30 years ago. Most early research interrogating the role that HSF1 plays in biology focused on its cytoprotective functions, as a factor that promotes the survival of organisms by protecting against the proteotoxicity associated with neurodegeneration and other pathological conditions. However, recent studies have revealed a deleterious role of HSF1, as a factor that is co-opted by cancer cells to promote their own survival to the detriment of the organism. In cancer, HSF1 operates in a multifaceted manner to promote oncogenic transformation, proliferation, metastatic dissemination, and anti-cancer drug resistance. Here we review our current understanding of HSF1 activation and function in malignant progression and discuss the potential for HSF1 inhibition as a novel anticancer strategy. Collectively, this ever-growing body of work points to a prominent role of HSF1 in nearly every aspect of carcinogenesis.

Copy number of 8q24.3 drives HSF1 expression and patient outcome in cancer: an individual patient data meta-analysis

Background

The heat-shock transcription factor 1 (HSF1) has been linked to cell proliferation and survival in cancer and has been proposed as a biomarker for poor prognosis. Here, we assessed the role of HSF1 expression in relation to copy number alteration (CNA) and cancer prognosis.

Methods

Using 10,287 cancer genomes from The Cancer Genome Atlas and Cbioportal databases, we assessed the association of HSF1 expression with CNA and cancer prognosis. CNA of 8q24.3 was categorized as diploid (reference), deletion (fewer copies), gain (+ 1 copy) and amplification (≥ + 2 copies). Multivariate logistic regression modeling was used to assess 5-year survival among those with a first cancer diagnosis and complete follow-up data (N = 9568), categorized per anatomical location and histology, assessing interaction with tumor stage, and expressed as odds ratios and 95% confidence intervals.

Results

We found that only 54.1% of all tumors have a normal predicted 8q24.3 copy number and that 8q24.3 located genes including HSF1 are mainly overexpressed due to increased copies number of 8q24.3 in different cancers. The tumor of patients having respectively gain (+ 1 copy) and amplification (≥ + 2 copies) of 8q24.3 display a global increase of 5-year mortality (odds ratio = 1.98, 95% CI 1.22–3.21) and (OR = 2.19, 1.13–4.26) after full adjustment. For separate cancer types, tumor patients with 8q24.3 deletion showed a marked increase of 5-year mortality in uterine (OR = 4.84, [2.75–8.51]), colorectal (OR = 4.12, [1.15–14.82]), and ovarian (OR = 1.83, [1.39–2.41]) cancers; and decreased mortality in kidney cancer (OR = 0.41, [0.21–0.82]). Gain of 8q24.3 resulted in significant mortality changes in 5-year mortality for cancer of the uterus (OR = 3.67, [2.03–6.66]), lung (OR = 1.76, [1.24–2.51]), colorectal (OR = 1.75, [1.32–2.31]) cancers; and amplification for uterine (OR = 4.58, [1.43–14.65]), prostate (OR = 4.41 [3.41–5.71]), head and neck (OR = 2.68, [2.17–3.30]), and stomach (OR = 0.56, [0.36–0.87]) cancers.

Conclusions

Here, we show that CNAs of 8q24.3 genes, including HSF1, are tightly linked to 8q24.3 copy number in tumor patients and can affect patient outcome. Our results indicate that the integration of 8q24.3 CNA detection may be a useful predictor for cancer prognosis.

Heat shock factor 1 (HSF1)-targeted anticancer therapeutics: overview of current preclinical progress

INTRODUCTION

The heat shock factor 1 (HSF1) plays a pivotal role in guarding proteome stability or proteostasis by induction of heat shock proteins (HSPs). While HSF1 remains mostly latent in unstressed normal cells, it is constitutively active in malignant cells, rendering them addicted to HSF1 for their growth and survival. HSF1 affects tumorigenesis, cancer progression, and treatment resistance by preserving cancer proteostasis, thus suggesting disruption of HSF1 activity as a potential anticancer strategy. Areas covered: In this review, we focus on the HSF1 activation cycle and its interaction with HSPs, the role of HSF1 in oncogenesis, and development of HSF1-targeted drugs as a potential anticancer therapy for disrupting cancer proteostasis. Expert opinion: HSF1 systematically maintains proteostasis in malignant cancer cells. Although genomic instability is widely accepted as a hallmark of cancer, little is known about the role of proteostasis in cancer. Unveiling the complicated mechanism of HSF1 regulation, particularly in cancer cells, will enable further development of proteostasis-targeted anticancer therapy.

ABBREVIATIONS

AMPK: AMP-activated protein kinase; DBD: DNA-binding domain; HR-A/B; HR-C: heptad repeats; HSE: heat shock elements; HSF1: heat shock factor; HSPs: heat shock proteins; HSR: heat shock response; MEK: mitogen-activated protein kinase kinase; mTOR: mammalian target of rapamycin; NF1: neurofibromatosis type 1; P-TEFb: positive transcription elongation factor b; RD: regulatory domain; RNAi: RNA interference; TAD: transactivation domain; TRiC: TCP-1 ring complex.

Functionally Related Genes Cluster into Genomic Regions That Coordinate Transcription at a Distance in Saccharomyces cerevisiae

The two-dimensional, physical positioning of genes along a chromosome can impact proper transcriptional regulation throughout a genomic region. The transcription of neighboring genes is correlated in a genome-wide manner, which is a characteristic of eukaryotes. Many coregulated gene families can be found clustered with another member of the same set—which can result in adjacent gene coregulation of the pair. Due to the myriad gene families that exhibit a nonrandom genomic distribution, there are likely multiple mechanisms working in concert to properly regulate transcriptional coordination of functionally clustered genes. In this study, we utilized budding yeast in an attempt to elucidate mechanisms that underlie this coregulation: testing and empirically validating the enhancer-promoter hypothesis in this species and reporting that functionally related genes cluster to genomic regions that are more conducive to transcriptional regulation at a distance. These clusters rely, in part, on chromatin maintenance and remodelers to maintain proper transcriptional coordination. Our work provides insight into the mechanisms underlying adjacent gene coregulation.

Balancing gene expression is a fundamental challenge of all cell types. To properly regulate transcription on a genome-wide level, there are myriad mechanisms employed by the cell. One layer to this regulation is through spatial positioning, with particular chromosomal loci exerting an influence on transcription throughout a region. Many coregulated gene families utilize spatial positioning to coordinate transcription, with functionally related genes clustering together which can allow coordinated expression via adjacent gene coregulation. The mechanisms underlying this process have not been elucidated, though there are many coregulated gene families that exhibit this genomic distribution. In the present study, we tested for a role for the enhancer-promoter (EP) hypothesis, which demonstrates that regulatory elements can exert transcriptional effects over a broad distance, in coordinating transcriptional coregulation using budding yeast, Saccharomyces cerevisiae. We empirically validated the EP model, finding that the genomic distance a promoter can affect varies by locus, which can profoundly affect levels of transcription, phenotype, and the extent of transcriptional disruption throughout a genomic region. Using the nitrogen metabolism, ribosomal protein, toxin response, and heat shock gene families as our test case, we report functionally clustered genes localize to genomic loci that are more conducive to transcriptional regulation at a distance compared to the unpaired members of the same families. Furthermore, we report that the coregulation of functional clusters is dependent, in part, on chromatin maintenance and remodeling, providing one mechanism underlying adjacent gene coregulation.

IMPORTANCE The two-dimensional, physical positioning of genes along a chromosome can impact proper transcriptional regulation throughout a genomic region. The transcription of neighboring genes is correlated in a genome-wide manner, which is a characteristic of eukaryotes. Many coregulated gene families can be found clustered with another member of the same set—which can result in adjacent gene coregulation of the pair. Due to the myriad gene families that exhibit a nonrandom genomic distribution, there are likely multiple mechanisms working in concert to properly regulate transcriptional coordination of functionally clustered genes. In this study, we utilized budding yeast in an attempt to elucidate mechanisms that underlie this coregulation: testing and empirically validating the enhancer-promoter hypothesis in this species and reporting that functionally related genes cluster to genomic regions that are more conducive to transcriptional regulation at a distance. These clusters rely, in part, on chromatin maintenance and remodelers to maintain proper transcriptional coordination. Our work provides insight into the mechanisms underlying adjacent gene coregulation.

Regional perturbation of gene transcription is associated with intrachromosomal rearrangements and gene fusion transcripts in high grade ovarian cancer

Genomic rearrangements are a hallmark of cancer biology and progression, allowing cells to rapidly transform through alterations in regulatory structures, changes in expression patterns, reprogramming of signaling pathways, and creation of novel transcripts via gene fusion events. Though functional gene fusions encoding oncogenic proteins are the most dramatic outcomes of genomic rearrangements, we investigated the relationship between rearrangements evidenced by fusion transcripts and local expression changes in cancer using transcriptome data alone. 9,953 gene fusion predictions from 418 primary serious ovarian cancer tumors were analyzed, identifying depletions of gene fusion breakpoints within coding regions of fused genes as well as an N-terminal enrichment of breakpoints within fused genes. We identified 48 genes with significant fusion-associated upregulation and furthermore demonstrate that significant regional overexpression of intact genes in patient transcriptomes occurs within 1 megabase of 78 novel gene fusions that function as central markers of these regions. We reveal that cancer transcriptomes select for gene fusions that preserve protein and protein domain coding potential. The association of gene fusion transcripts with neighboring gene overexpression supports rearrangements as mechanism through which cancer cells remodel their transcriptomes and identifies a new way to utilize gene fusions as indicators of regional expression changes in diseased cells with only transcriptomic data.