ZFHX3 is indispensable for ERβ to inhibit cell proliferation via MYC downregulation in prostate cancer cells

Somatic mutational landscape across Indian breast cancer cases by whole exome sequencing

Breast cancer (BC) has emerged as the most common malignancy among females. The genomic profile of BC is diverse in nature and complex due to heterogeneity among various geographically different ethnic groups. The primary objective of this study was to carry out a comprehensive mutational analysis of Indian BC cases by performing whole exome sequencing. The cohort included patients with a median age of 48 years. TTN, TP53, MUC16, SYNE1, and OBSCN were the frequently altered genes found in our cohort. The PIK3CA and KLC3 genes are driver genes implicated in various cellular functions and cargo transportation through microtubules, respectively. Except for CCDC168 and PIK3CA, several gene pairings were found to be significantly linked with co-occurrence. Irrespective of their hormonal receptor status, RTK/RAS was observed with frequently altered signaling pathways. Further analysis of the mutational signature revealed that SBS13, SBS6, and SBS29 were mainly observed in our cohort. This study supplements the discovery of diagnostic biomarkers and provides new therapeutic options for the improved management of BC.

Evaluation of Potential Roles of Zinc Finger Homeobox 3 (Zfhx3) Expressed in Chondrocytes and Osteoblasts on Skeletal Growth in Mice

Bone formation is tightly modulated by genetically encoded molecular proteins that interact to regulate cellular differentiation and secretion of bony matrix. Many transcription factors are known to coordinate these events by controlling gene transcription within networks. However, not all factors involved are known. Here, we identified a novel function for Zinc Finger Homeobox 3 (Zfhx3), a gene encoding a transcription factor, as a regulator of bone metabolism. We knocked out Zfhx3 conditionally in mice in either chondrocytes or osteoblasts and characterized their bones by micro-CT in 12-week-old mice. We observed a negative effect in linear bone growth in both knockout mice but reduced bone mass only in mice with Zfhx3 deleted in osteoblasts. Loss of Zfhx3 expression in osteoblasts affected trabecular bone mass in femurs and vertebrae in both sexes but influenced cortical bone volume fraction only in females. Moreover, transcriptional analysis of femoral bones in osteoblast Zfhx3 conditional knockout mice revealed a reduced expression of osteoblast genes, and histological evaluation of trabecular bones suggests that Zfhx3 causes changes in bone formation and not resorption. The loss of Zfhx3 causes reductions in trabecular bone area and osteoid volume, but no changes in the expression of osteoclast differentiation markers or number of TRAP stained osteoclasts. These studies introduce Zfhx3 as a relevant factor toward understanding gene regulatory networks that control bone formation and development of peak bone mass.

ZFHX3 acts as a tumor suppressor in prostate cancer by targeting FTO-mediated m6A demethylation

Zinc-finger homeobox 3 (ZFHX3, also known as ATBF1) suppresses prostatic tumorigenesis. ZFHX3 is frequently found to have numerous deletions in human prostate cancer (PCa). However, the underlying molecular function of ZFHX3 during prostatic tumorigenesis is not well understood. N6-methyladenosine (m6A) modification in RNA plays a critical role in the development of cancers; however, the relationship between ZFHX3 and m6A modification is largely unknown in PCa. In this study, we found that ZFHX3 knockdown decreased total m6A levels through enhancing the transcriptional activity of FTO in PCa cells. Importantly, FTO inhibition suppressed cell proliferation and rescued the promoting function of ZFHX3 knockdown on cell proliferation. In vivo, we verified that FTO was upregulated and ZFHX3 was decreased in PCa patients and that a high level of ZFHX3 is indispensable for low FTO expression and is correlated with better patient survival. Through transcriptome sequencing and MeRIP sequencing, we revealed that E2F2 and CDKN2C were the direct targets of FTO-mediated m6A modification and ZFXH3 was required for the regulation of FTO on E2F2 and CDKN2C expression. Unexpectedly, we uncovered that ZFHX3 expression was in return regulated by FTO in an m6A-dependent way. These findings establish a novel crosstalk mechanism between ZFHX3 and FTO in prostatic tumorigenesis.

LncRNA-SERB promotes vasculogenic mimicry (VM) formation and tumor metastasis in renal cell carcinoma

A growing body of evidence shows that vasculogenic mimicry (VM) is closely related to the invasion and metastasis of many tumor cells. Although the estrogen receptor (ER) can promote initiation and progression of renal cell carcinoma (RCC), how the downstream biomolecules are involved, and the detailed mechanisms of how ER expression is elevated in RCC remain to be further elucidated. Here, we discovered that long noncoding RNA (LncRNA)-SERB is highly expressed in tumor cells of RCC patients. We used multiple RCC cells and an in vivo mouse model for our study, and results indicated that LncRNA-SERB could boost RCC VM formation and cell invasion in vitro and in vivo. Although a previous report showed that ERβ can affect the VM formation in RCC, it is unclear which factor could upregulate ERβ. This is the first study to show LncRNA-SERB can be the upstream regulator of ERβ to control RCC progression. Mechanistically, LncRNA-SERB may increase ERβ via binding to the promoter area, and ERβ functions through transcriptional regulation of zinc finger E-box binding homeobox 1 (ZEB1) to regulate VM formation. These results suggest that LncRNA-SERB promotes RCC cell VM formation and invasion by upregulating the ERβ/ZEB1 axis and that therapeutic targeting of this newly identified pathway may better inhibit RCC progression.

Follicular lymphoma B cells exhibit heterogeneous transcriptional states with associated somatic alterations and tumor microenvironments

Follicular lymphoma (FL) is an indolent non-Hodgkin lymphoma of germinal center origin, which presents with significant biologic and clinical heterogeneity. Using RNA-seq on B cells sorted from 87 FL biopsies, combined with machine-learning approaches, we identify 3 transcriptional states that divide the biological ontology of FL B cells into inflamed, proliferative, and chromatin-modifying states, with relationship to prior GC B cell phenotypes. When integrated with whole-exome sequencing and immune profiling, we find that each state was associated with a combination of mutations in chromatin modifiers, copy-number alterations to TNFAIP3, and T follicular helper cells (Tfh) cell interactions, or primarily by a microenvironment rich in activated T cells. Altogether, these data define FL B cell transcriptional states across a large cohort of patients, contribute to our understanding of FL heterogeneity at the tumor cell level, and provide a foundation for guiding therapeutic intervention.

S100A2 induces epithelial-mesenchymal transition and metastasis in pancreatic cancer by coordinating transforming growth factor β signaling in SMAD4-dependent manner

Pancreatic ductal adenocarcinoma (PDAC) is a highly aggressive tumor and is associated with a poor prognosis. Treatment strategies for PDAC are largely ineffective primarily because of delay in its diagnosis and limited efficacy of systematic treatment. S100A2 is associated with the proliferation, migration, and differentiation of several tumors; however, its effects on PDAC and the associated molecular mechanisms remain to be explored. We studied the mechanisms underlying the effect of S100A2 on epithelial-mesenchymal transition (EMT) and metastasis in PDAC cells. We found that the level of S100A2 remarkably increased and was associated with poor PDAC prognosis. The overexpression of S100A2 in PANC-1 cells also induced EMT, in addition to increasing the invasion and migration of PDAC cells, whereas the knockdown of S100A2 markedly inhibited cell metastasis. Furthermore, S100A2 was found to enhance metastatic abilities in vivo. The overexpression of S100A2 increased SMAD4 expression, whereas the knockdown of S100A2 reduced SMAD4 expression. SMAD4 overexpression could effectively rescue the effects of S100A2 knockdown on EMT. S100A2 mechanistically activated the transforming growth factor (TGF)-β/Smad2/3 signaling pathway, upregulated SMAD4 expression, induced EMT, and increased PANC-1 cell metastasis. In conclusion, the S100A2/SMAD4 axis modulates EMT to accelerate PDAC development. Our results supplement and enrich the understanding of the pathogenesis underlying PDAC and provide a new theoretical basis and strategy targeting S100A2 for the diagnosis and treatment of PDAC.

To be, or not to be: the dilemma of immunotherapy for non-small cell lung cancer harboring various driver mutations

INTRODUCTION

Lung cancer is one of primary cancer type with high incidence and mortality, non-small cell lung cancer (NSCLC) is the most common type of lung cncer. For advanced lung cancer, traditional chemotherapy and targeted therapy become difficult to solve the dilemma of further progress. In recent years, with the clinical application of immunotherapy, the therapeutic strategy of lung cancer has changed dramatically. At present, immunotherapy has shown conspicuous efficacy in NSCLC patients with high expression of programmed death-ligand 1 (PD-L1) and high tumor mutational burden (TMB). The discovery of driver mutations brings delightful hope for targeted cancer therapy. However, it remains controversial whether immunotherapy can be used in NSCLC patients with these specific driver mutations.

METHOD

This article summarized the latest research progresses of immunotherapy in advanced NSCLC. We paid close attention to the relevance of various driver mutations and immunotherapy in NSCLC patients, and summarized the predictive effects of several driver mutations and immunotherapy.

RESULTS

The mutations of KRAS, KRAS+TP53, EPHA (especially EPHA5), ZFHX3, ZFHX3+TP53, NOTCH, BRAF and LRP1B+FAT3 have potential to be used as biomarkers to predict the positive effectiveness of immunotherapy. ZFHX3, ZFHX3+TP53, STKII/LKB1+KEAP1+SMARCA4+PBRM1 mutations in LUAD patients get more positive effect in immunotherapy. While the mutations of EGFR, KEAP1, STKII/LKB1+KRAS, EML4-ALK, MET exon 14 skipping mutation, PBRM1, STKII/LKB1+KEAP1+SMARCA4+PBRM1, ERBB2, PIK3CA and RET often indicate poor benefit from immunotherapy.

CONCLUSION

Many gene mutations have been shown to be associated with immunotherapy efficacy. Gene mutations should be combined with PD-L1, TMB, etc. to predict the effect of immunotherapy.

Using evolutionary constraint to define novel candidate driver genes in medulloblastoma

Current knowledge of cancer genomics remains biased against noncoding mutations. To systematically search for regulatory noncoding mutations, we assessed mutations in conserved positions in the genome under the assumption that these are more likely to be functional than mutations in positions with low conservation. To this end, we use whole-genome sequencing data from the International Cancer Genome Consortium and combined it with evolutionary constraint inferred from 240 mammals, to identify genes enriched in noncoding constraint mutations (NCCMs), mutations likely to be regulatory in nature. We compare medulloblastoma (MB), which is malignant, to pilocytic astrocytoma (PA), a primarily benign tumor, and find highly different NCCM frequencies between the two, in agreement with the fact that malignant cancers tend to have more mutations. In PA, a high NCCM frequency only affects the BRAF locus, which is the most commonly mutated gene in PA. In contrast, in MB, >500 genes have high levels of NCCMs. Intriguingly, several loci with NCCMs in MB are associated with different ages of onset, such as the HOXB cluster in young MB patients. In adult patients, NCCMs occurred in, e.g., the WASF-2/AHDC1/FGR locus. One of these NCCMs led to increased expression of the SRC kinase FGR and augmented responsiveness of MB cells to dasatinib, a SRC kinase inhibitor. Our analysis thus points to different molecular pathways in different patient groups. These newly identified putative candidate driver mutations may aid in patient stratification in MB and could be valuable for future selection of personalized treatment options.

WFDC3 inhibits tumor metastasis by promoting the ERβ-mediated transcriptional repression of TGFBR1 in colorectal cancer

Estrogen plays a protective role in colorectal cancer (CRC) and primarily functions through estrogen receptor β (ERβ). However, clinical strategies for CRC therapy associated with ERβ are still under investigation. Our discoveries identified WFDC3 as a tumor suppressor that facilitates estrogen-induced inhibition of metastasis through the ERβ/TGFBR1 signaling axis. WFDC3 interacts with ERβ and increases its protein stability by inhibiting its proteasome-dependent degradation. WFDC3 represses TGFBR1 expression through ERβ-mediated transcription. Blocking TGFβ signaling with galunisertib, a drug used in clinical trials that targets TGFBR1, impaired the migration of CRC cells induced by WFDC3 depletion. Moreover, there was clinical significance to WFDC3 in CRC, as CRC patients with high WFDC3 expression in tumor cells had favorable prognoses. Therefore, this work suggests that WFDC3 could be an indicator for therapies targeting the estrogen/ERβ pathway in CRC patients.

Prediction of prostate cancer biochemical recurrence by using discretization supports the critical contribution of the extra-cellular matrix genes

Due to its complexity, much effort has been devoted to the development of biomarkers for prostate cancer that have acquired the utmost clinical relevance for diagnosis and grading. However, all of these advances are limited due to the relatively large percentage of biochemical recurrence (BCR) and the limited strategies for follow up. This work proposes a methodology that uses discretization to predict prostate cancer BCR while optimizing the necessary variables. We used discretization of RNA-seq data to increase the prediction of biochemical recurrence and retrieve a subset of ten genes functionally known to be related to the tissue structure. Equal width and equal frequency data discretization methods were compared to isolate the contribution of the genes and their interval of action, simultaneously. Adding a robust clinical biomarker such as prostate specific antigen (PSA) improved the prediction of BCR. Discretization allowed classifying the cancer patients with an accuracy of 82% on testing datasets, and 75% on a validation dataset when a five-bin discretization by equal width was used. After data pre-processing, feature selection and classification, our predictions had a precision of 71% (testing dataset: MSKCC and GSE54460) and 69% (Validation dataset: GSE70769) should the patients present BCR up to 24 months after their final treatment. These results emphasize the use of equal width discretization as a pre-processing step to improve classification for a limited number of genes in the signature. Functionally, many of these genes have a direct or expected role in tissue structure and extracellular matrix organization. The processing steps presented in this study are also applicable to other cancer types to increase the speed and accuracy of the models in diverse datasets.

Synthesis and evaluation of radioiodinated estrogens for diagnosis and therapy of male urogenital tumours

The preparation of 24 estrogens, their estrogen receptor (ER) affinity and studies of radioiodinated estrogen binding to ER-positive male bladder tumor cells (HTB9) are described. The estrogens with the highest affinity were selected using fluorescence anisotropy assays. A 2,2,2-trifluoroethyl group at the 11β-position caused particularly promising affinity. (Radio)iodination was performed on the 17α-vinyl group. Binding studies on HTB9 cells revealed picomolar affinities of radioconjugates 19 and 31, indicating promising ability for targeting of urogenital tumors.

Genomic, epigenomic, and transcriptomic signatures of prostate cancer between African American and European American patients

Prostate cancer is the second most common cancer in men in the United States, and racial disparities are greatly observed in the disease. Specifically, African American (AA) patients have 60% higher incidence and mortality rates, in addition to higher grade and stage prostate tumors, than European American (EA) patients. In order to narrow the gap between clinical outcomes for these two populations, genetic and molecular signatures contributing to this disparity have been characterized. Over the past decade, profiles of prostate tumor samples from different ethnic groups have been developed using molecular and functional assays coupled with next generation sequencing or microarrays. Comparative genome-wide analyses of genomic, epigenomic, and transcriptomic profiles from prostate tumor samples have uncovered potential race-specific mutations, copy number alterations, DNA methylation, and gene expression patterns. In this study, we reviewed over 20 published studies that examined the aforementioned molecular contributions to racial disparities in AA and EA prostate cancer patients. The reviewed genomic studies revealed mutations, deletions, amplifications, duplications, or fusion genes differentially enriched in AA patients relative to EA patients. Commonly reported genomic alterations included mutations or copy number alterations of FOXA1, KMT2D, SPOP, MYC, PTEN, TP53, ZFHX3, and the TMPRSS2-ERG fusion. The reviewed epigenomic studies identified that CpG sites near the promoters of PMEPA1, RARB, SNRPN, and TIMP3 genes were differentially methylated between AA and EA patients. Lastly, the reviewed transcriptomic studies identified genes (e.g. CCL4, CHRM3, CRYBB2, CXCR4, GALR1, GSTM3, SPINK1) and signaling pathways dysregulated between AA and EA patients. The most frequently found dysregulated pathways were involved in immune and inflammatory responses and neuroactive ligand signaling. Overall, we observed that the genomic, epigenomic, and transcriptomic alterations evaluated between AA and EA prostate cancer patients varied between studies, highlighting the impact of using different methods and sample sizes. The reported genomic, epigenomic, and transcriptomic alterations do not only uncover molecular mechanisms of tumorigenesis but also provide researchers and clinicians valuable resources to identify novel biomarkers and treatment modalities to improve the disparity of clinical outcomes between AA and EA patients.

ATBF1 is a potential diagnostic marker of histological grade and functions via WNT5A in breast cancer

BACKGROUND

Histological grade has been demonstrated to be an important factor of breast cancer outcome and is associated with cell differentiation and is currently being evaluated via H&E-stained sections. Molecular biomarkers are essential to improve the accuracy of histological grading. ATBF1, a large transcription factor, has been considered a tumor suppressor gene with frequent mutations or deletions in multiple cancers. In breast cancer, ATBF1 was reported to function in cell differentiation and mammary development. However, its role in the clinic has rarely been reported.

METHODS

Breast cancer tissues (BCTs) and adjacent noncancerous tissues (ANCTs) were collected to analyze the expression of ATBF1 at the mRNA and protein levels. Three anti-ATBF1 antibodies recognizing independent peptides of ATBF1 (N-terminal end, middle region and C-terminal end) were applied for IHC staining. Small interfering RNA (siRNA) was used to silence ATBF1 expression and to investigate the roles of ATBF1 in MCF7 cells. Microarrays were introduced to analyze the differentially expressed genes, enriched GO terms and KEGG terms regulated by ATBF1 and its potential downstream genes, which were further confirmed in vitro and in clinical samples.

RESULTS

The expression of ATBF1 was reduced in BCTs at both the mRNA and protein levels compared with that in ANCTs. ATBF1 protein was predominantly localized in the nucleus of ANCTs but in the cytoplasm of BCTs. Both the mRNA and protein levels of ATBF1 were significantly correlated with histological grade. Consistently, knockdown of ATBF1 increased stemness marker expression and reduced differentiation markers in vitro. Further analysis identified WNT5A as an essential downstream gene of ATBF1 in breast cancer cells. Treatment of WNT5A disrupted cell proliferation induced by ATBF1 silencing. In BCTs, a significant correlation was observed between the expression of WNT5A and ATBF1.

CONCLUSION

The results indicated that ATBF1 expression might be a useful diagnostic marker associated with histological grade and breast cancer malignancy. WNT5A and its signaling pathway are novel mechanisms by which ATBF1 contributes to breast cancer tumorigenesis.

A six-gene prognostic signature for both adult and pediatric acute myeloid leukemia identified with machine learning

BACKGROUND

Although it is well-known that adult and pediatric acute myeloid leukemias (AMLs) are genetically distinct diseases, they still share certain gene expression profiles. The age-related genetic heterogeneities of AMLs have been well-studied, but the common prognostic signatures and molecular mechanisms of adult and pediatric AMLs are less investigated.

AIM

To identify genes and pathways that are associated with both pediatric and adult AMLs and discover a gene signature for overall survival (OS) prediction.

METHODS

Through mining the transcriptome profiles of The Cancer Genome Atlas (TCGA) data sets of adult cancers and The Therapeutically Applicable Research to Generate Effective Treatments (TARGET) data of pediatric cancers, we identified genes that are commonly dysregulated in both pediatric and adult AMLs, further discovered a common gene signature, and built two risk score models for TCGA and TARGET cohorts, respectively with L ₀ regularized global AUC (area under the receiver operating characteristic curve) summary maximization.

RESULTS

We identified 57 genes that are differentially expressed and prognostically significant in both adult and childhood AMLs. The top 4 Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways enriched with those 57 genes include transcriptional misregulation, focal adhesion, PI3K-Akt signaling pathway, and signaling pathways regulating pluripotency of stem cells. We further identified a 6-gene signature including genes of ADAMTS3, DNMT3B, NYNRIN, SORT1, ZFHX3, and ZG16B for risk prediction. We constructed a risk score model with one dataset (either TCGA or TARGET) and evaluated its performance with the other. The test AUCs for the risk prediction of TCGA data with a 2-year and 5-year OS cutoffs are 0.762 (P = 2.33e-13, 95% CI: 0.69-0.83) and 0.759 (P = 7.26e-08, 95% CI: 0.66-0.85), respectively, while the test AUCs of TARGET data with the same cutoffs are 0.71 (P = 3.3e-07, 95% CI: 0.62-0.79) and 0.72 (P= 5.25e-09, 95% CI: 0.65-0.80), respectively. We further stratified patients into 3 equal sized prognostic subtypes with the 6-gene risk scores. The P-values of the tertile partitions are 1.74e-07 and 3.28e-08 for the TARGET and TCGA cohorts, respectively, which are significantly better than the standard cytogenetic risk stratification of both cohorts (TARGET: P = 1.64e-06; TCGA: P = 1.79e-05). When validated with two other independent cohorts, the 6-gene risk score models remain a significant predictor for OS. Investigating the common gene expression program is significant in that we may extrapolate the findings from adults to children and avoid unnecessary pediatric clinical trials.

Chromatin accessibility landscape and active transcription factors in primary human invasive lobular and ductal breast carcinomas

BACKGROUND

Invasive lobular breast carcinoma (ILC), the second most prevalent histological subtype of breast cancer, exhibits unique molecular features compared with the more common invasive ductal carcinoma (IDC). While genomic and transcriptomic features of ILC and IDC have been characterized, genome-wide chromatin accessibility pattern differences between ILC and IDC remain largely unexplored.

METHODS

Here, we characterized tumor-intrinsic chromatin accessibility differences between ILC and IDC using primary tumors from The Cancer Genome Atlas (TCGA) breast cancer assay for transposase-accessible chromatin with sequencing (ATAC-seq) dataset.

RESULTS

We identified distinct patterns of genome-wide chromatin accessibility in ILC and IDC. Inferred patient-specific transcription factor (TF) motif activities revealed regulatory differences between and within ILC and IDC tumors. EGR1, RUNX3, TP63, STAT6, SOX family, and TEAD family TFs were higher in ILC, while ATF4, PBX3, SPDEF, PITX family, and FOX family TFs were higher in IDC.

CONCLUSIONS

This study reveals the distinct epigenomic features of ILC and IDC and the active TFs driving cancer progression that may provide valuable information on patient prognosis.

Comprehensive Clinical Analysis of Gallbladder Neuroendocrine Neoplasms: A Large-Volume Multicenter Study During One Decade

BACKGROUND

This study aimed to comprehensively investigate the clinicopathologic characteristics and therapeutic situations of gallbladder neuroendocrine neoplasms (GB-NENs) in the real world via a multicenter, large-scale cohort study.

METHODS

The study searched for patients in 143 hospitals in China and enrolled 154 patients with GB-NENs diagnosed in 40 hospitals between 2004 and 2021. Clinicopathologic characteristics and therapeutic approaches were analyzed retrospectively.

RESULTS

The median age at the initial diagnosis of the patients with GB-NENs was 63 years (range 33-83 years), and 61.7% of the patients were women. Tumor-node-metastasis staging classified 92 patients as stage 3 or above. Based on the 2019 World Health Organization classification, 96 cases (62.3%) were confirmed pathologically as poorly differentiated neuroendocrine carcinomas, 13 cases (8.4%) as well-differentiated neuroendocrine tumors, and 45 cases as mixed neuroendocrine-non-neuroendocrine neoplasms. The liver was the most frequent metastatic site. Immunohistochemistry showed that synaptophysin was most frequently positive (80.4%), followed by chromogranin A (61.7%), and CD56 (58.4%). Computed tomography and magnetic resonance imaging showed more common clear boundaries (25/39 cases) and invasive growth features (27 cases). None of these cases had an accurate diagnosis before surgery, with a misdiagnosis rate of 100%. Surgical resection is the main treatment, and platinum-based chemotherapeutic regimens were preferred as adjuvant therapies for patients with GB-NENs. The available survival data for 74 patients showed an overall survival rate of 59% at 1 year, 33% at 3 years, and 29% at 5 years. No significant difference was found between the patients treated with and those treated without adjuvant chemotherapy.

CONCLUSIONS

Gallbladder neuroendocrine neoplasms have high malignancy and a poor prognosis. Importantly, this large-scale cohort study significantly improves our understanding of GB-NENs and will benefit the exploration of its mechanism and treatment modes. Further investigation is necessary to explore the management of this disease.

Transcriptome-wide prediction of prostate cancer gene expression from histopathology images using co-expression based convolutional neural networks

Motivation

Molecular phenotyping by gene expression profiling is central in contemporary cancer research and in molecular diagnostics but remains resource intense to implement. Changes in gene expression occurring in tumours cause morphological changes in tissue, which can be observed on the microscopic level. The relationship between morphological patterns and some of the molecular phenotypes can be exploited to predict molecular phenotypes from routine haematoxylin and eosin-stained whole slide images (WSIs) using convolutional neural networks (CNNs). In this study, we propose a new, computationally efficient approach to model relationships between morphology and gene expression.

Results

We conducted the first transcriptome-wide analysis in prostate cancer, using CNNs to predict bulk RNA-sequencing estimates from WSIs for 370 patients from the TCGA PRAD study. Out of 15 586 protein coding transcripts, 6618 had predicted expression significantly associated with RNA-seq estimates (FDR-adjusted P-value <1×10⁻⁴) in a cross-validation and 5419 (81.9%) of these associations were subsequently validated in a held-out test set. We furthermore predicted the prognostic cell-cycle progression score directly from WSIs. These findings suggest that contemporary computer vision models offer an inexpensive and scalable solution for prediction of gene expression phenotypes directly from WSIs, providing opportunity for cost-effective large-scale research studies and molecular diagnostics.

Availability and implementation

A self-contained example is available from http://github.com/phiwei/prostate_coexpression. Model predictions and metrics are available from doi.org/10.5281/zenodo.4739097.

Supplementary information

Supplementary data are available at Bioinformatics online.

Proteogenomic characterization of 2002 human cancers reveals pan-cancer molecular subtypes and associated pathways

Mass-spectrometry-based proteomic data on human tumors—combined with corresponding multi-omics data—present opportunities for systematic and pan-cancer proteogenomic analyses. Here, we assemble a compendium dataset of proteomics data of 2002 primary tumors from 14 cancer types and 17 studies. Protein expression of genes broadly correlates with corresponding mRNA levels or copy number alterations (CNAs) across tumors, but with notable exceptions. Based on unsupervised clustering, tumors separate into 11 distinct proteome-based subtypes spanning multiple tissue-based cancer types. Two subtypes are enriched for brain tumors, one subtype associating with MYC, Wnt, and Hippo pathways and high CNA burden, and another subtype associating with metabolic pathways and low CNA burden. Somatic alteration of genes in a pathway associates with higher pathway activity as inferred by proteome or transcriptome data. A substantial fraction of cancers shows high MYC pathway activity without MYC copy gain but with mutations in genes with noncanonical roles in MYC. Our proteogenomics survey reveals the interplay between genome and proteome across tumor lineages.

Pan-cancer proteomics analysis enables the analysis of protein expression across multiple cancer types. Here, the authors compare proteomics from 14 cancer types and show 11 distinct subtypes across multiple cancer types. Proteome data could link higher pathway activity levels with somatic alteration of specific genes in the pathway.

Clinical utility of PDX cohorts to reveal biomarkers of intrinsic resistance and clonal architecture changes underlying acquired resistance to cetuximab in HNSCC

Cetuximab is a widely used drug for treating head and neck squamous cell carcinomas (HNSCCs); however, it provides restricted clinical benefits, and its response duration is limited by drug resistance. Here, we conducted randomized “Phase II-like clinical trials” of 49 HNSCC PDX models and reveal multiple informative biomarkers for intrinsic resistance to cetuximab (e.g., amplification of ANKH, up-regulation of PARP3). After validating these intrinsic resistance biomarkers in another HNSCC PDX cohort (61 PDX models), we generated acquired cetuximab resistance PDX models and analyzed them to uncover resistance mechanisms. Whole exome sequencing and transcriptome sequencing revealed diverse patterns of clonal selection in acquired resistant PDXs, including the emergence of subclones with strongly activated RAS/MAPK. Extending these insights, we show that a combination of a RAC1/RAC3 dual-target inhibitor and cetuximab could overcome acquired cetuximab resistance in vitro and in vivo. Beyond revealing intrinsic resistance biomarkers, our PDX-based study shows how clonal architecture changes underlying acquired resistance can be targeted to expand the therapeutic utility of this important drug to more HNSCC patients.

Comprehensive Genomic and Epigenomic Analyses on Transcriptomic Regulation in Stomach Adenocarcinoma

Background: DNA methylation (MET)–mediated transcriptomic disturbance and copy number variations (CNVs) exert a significant influence in stimulating the heterogeneous progression of stomach adenocarcinoma (STAD). Nevertheless, the relation of DNA MET with CNVs, together with its impact on tumor occurrence, is still unclear.

Methods: The messenger RNA (mRNA) expression (EXP) profiles, DNA MET, and DNA copy numbers, together with STAD mutation data, were collected from the TCGA official data portal. We employed circular binary segmentation algorithm in “DNAcopy.” library of R package for mapping DNA CNV data at genetic level for all samples based on the segmented CNV data. Stable clusters of samples were recognized using negative matrix factorization cluster analysis based on 50 iterations and the “brunet” method using the MET-correlated (METcor) and CNV-correlated (CNVcor) genes. The R package “iCluster” method was utilized to comprehensively analyze the EXP, MET, and DNA CNV profiles.

Results: A total of 313 STAD samples were isolated for checking DNA copy numbers and MET and for measuring EXP. In accordance with our results, we discovered obvious co-regulation of CNVcor genes and METcor counterparts. Apart from that, these genes were subject to multi-omics integration. Meanwhile, three subtypes of STAD were detected and confirmed based on independent data. Among them, the subtype with increased aggressiveness was related to decreased mutation frequencies of ARID1A, PIK3CA, ZFHX3, SPECC1, OBSCN, KMT2D, FSIP2, ZBTB20, TTN, and RANBP2, together with the abnormal levels of JPH3, KCNB1, and PLCXD3.

Conclusion: According to the results, these aforementioned genes exerted crucial roles in the development of invasive STAD. Our findings on transcriptomic regulation genomically and epigenetically facilitate the understanding of the STAD pathology from different aspects, which help to develop efficient anti-STAD therapy.

Identification of Tissue Types and Gene Mutations From Histopathology Images for Advancing Colorectal Cancer Biology

Objective: Colorectal cancer (CRC) patients respond differently to treatments and are sub-classified by different approaches. We evaluated a deep learning model, which adopted endoscopic knowledge learnt from AI-doscopist, to characterise CRC patients by histopathological features. Results: Data of 461 patients were collected from TCGA-COAD database. The proposed framework was able to 1) differentiate tumour from normal tissues with an Area Under Receiver Operating Characteristic curve (AUROC) of 0.97; 2) identify certain gene mutations (MYH9, TP53) with an AUROC > 0.75; 3) classify CMS2 and CMS4 better than the other subtypes; and 4) demonstrate the generalizability of predicting KRAS mutants in an external cohort. Conclusions: Artificial intelligent can be used for on-site patient classification. Although KRAS mutants were commonly associated with therapeutic resistance and poor prognosis, subjects with predicted KRAS mutants in this study have a higher survival rate in 30 months after diagnoses.

Exonic variants in multiple myeloma patients associated with relapsed/ refractory and response to bortezomib regimens

Novel treatment in multiple myeloma represented by proteasome inhibitors, immunomodulatory drugs and monoclonal antibodies have produced a deep response. However, relapses are possible, and all classes of drugs are refractory to patients. Next-generation sequencing has improved our understanding of the multiple myeloma genome related to drug resistance and has discovered many genomic variants. Therefore, this study was conducted to investigate new variants associated with drug resistance in MM patients who relapsed and refractory to bortezomib regimen and daratumumab treatment using next-generation sequencing for whole-exome sequencing. Peripheral blood samples were collected in EDTA tubes from six patients; four were in relapsed and refractory to bortezomib regimens and daratumumab; two patients responded to bortezomib regimens. Whole-exome sequencing was performed by the MGI-DNBSEQ-G400 instrument. We identified 21 variants in multiple myeloma patients. Seventeen variants were found in relapsed and refractory multiple myeloma in 11 genes (GNAQ, PMS1, CREB1, NSUNS2, PIK3CG, ROS1, PMS2, FIT4, KDM5A, STK11 and ZFHX3). And four variants were identified in two patients with response to bortezomib regimens in 4 genes (RAF1, CREB1, ZFHX3 and INSR). We have observed several genetic variants in many genes that may have been associated with the poor prognosis and poor response to treatment in these patients. These values should be further confirmed in large sample studies using the RNA-seq technique to identify genome expression.

AR imposes different effects on ZFHX3 transcription depending on androgen status in prostate cancer cells

Both androgen receptor (AR) and the ZFHX3 transcription factor modulate prostate development. While AR drives prostatic carcinogenesis, ZFHX3 is a tumour suppressor whose loss activates the PI3K/AKT signalling in advanced prostate cancer (PCa). However, it is unknown whether ZFHX3 and AR are functionally related in PCa cells and, if so, how. Here, we report that in AR-positive LNCaP and C4-2B PCa cells, androgen upregulates ZFHX3 transcription via androgen-induced AR binding to the androgen-responsive elements (AREs) of the ZFHX3 promoter. Androgen also upregulated ZFHX3 transcription in vivo, as castration dramatically reduced Zfhx3 mRNA and protein levels in mouse prostates, and ZFHX3 mRNA levels correlated with AR activities in human PCa. Interestingly, the binding of AR to one ARE occurred in the absence of androgen, and the binding repressed ZFHX3 transcription as this repressive binding was interrupted by androgen treatment. The enzalutamide antiandrogen prevented androgen from inducing ZFHX3 transcription and caused excess ZFHX3 protein degradation. In human PCa, ZFHX3 was downregulated and the downregulation correlated with worse patient survival. These findings establish a regulatory relationship between AR and ZFHX3, suggest a role of ZFHX3 in AR function and implicate ZFHX3 loss in the antiandrogen therapies of PCa.

Clinicopathological features of tumor mutation burden, Epstein-Barr virus infection, microsatellite instability and PD-L1 status in Chinese patients with gastric cancer

Objectives

Gastric cancer (GC) is the 4th most common type of cancer worldwide. Different GC subtypes have unique molecular features that may have different therapeutic methods. The aim of the present study was to investigate Epstein-Barr virus (EBV) infection, microsatellite instability (MSI) status, the expression of programmed death-ligand 1 (PD-L1) and gene mutations in GC patients.

Methods

The data of 2504 GC patients, who underwent curative gastrectomy with lymphadenectomy at Peking University Cancer Hospital between 2013 and 2018, were reviewed. We analyzed the clinicopathological factors associated with the immunohistochemistry (IHC) profiles of these patients, and genetic alterations were analyzed using next generation sequencing (NGS).

Results

Mismatch repair-deficient (d-MMR) GC patients were found to have a higher probability of expressing PD-L1 (p = 0.000, PD-L1 cutoff value = 1%). In addition, 4 and 6.9% of the 2504 gastric cancer patients were EBV-positive and d-MMR, respectively. The number of MLH1/PMS2-negative cases was 126 (6%), and the number of MSH2/MSH6-negative cases was 14 (0.9%). d-MMR status was associated with a intestinal group (p = 0.012), but not with tumor differentiation. Furthermore, MSI and d-MMR GC status (detected by NGS and IHC, respectively) were consistently high, and the rate of MSI was higher in patients with d-MMR GC. A number of genes associated with DNA damage repair were detected in GC patients with MSI, including POLE, ETV6, BRCA and RNF43. In patients with a high tumor mutation burden, the most significantly mutated genes were LRP1B (79.07%), ARID1A (74.42%), RNF43 (69.77%), ZFHX3 (65.12%), TP53 (58.14%), GANS (51.16%), BRCA2 (51.16%), PIK3CA (51.16%), NOTCH1 (51.16%), SMARCA4 (48.84%), ATR (46.51%), POLE (41.86%) and ATM (39.53%).

Conclusions

Using IHC and NGS, MSI status, protein expression, tumor mutation burden (TMB) and genetic alterations were identified in patients with GC, which provides a theoretical basis for the future clinical treatment of GC.

Reprogramming enriches for somatic cell clones with small-scale mutations in cancer-associated genes

Cellular therapies based on induced pluripotent stem cells (iPSCs) come out of age and an increasing number of clinical trials applying iPSC-based transplants are ongoing or in preparation. Recent studies, however, demonstrated a high number of small-scale mutations in iPSCs. Although the mutational load in iPSCs seems to be largely derived from their parental cells, it is still unknown whether reprogramming may enrich for individual mutations that could lead to loss of functionality and tumor formation from iPSC derivatives. 30 hiPSC lines were analyzed by whole exome sequencing. High accuracy amplicon sequencing showed that all analyzed small-scale variants pre-existed in their parental cells and that individual mutations present in small subpopulations of parental cells become enriched among hiPSC clones during reprogramming. Among those, putatively actionable driver mutations affect genes related to cell-cycle control, cell death, and pluripotency and may confer a selective advantage during reprogramming. Finally, a short hairpin RNA (shRNA)-based experimental approach was applied to provide additional evidence for the individual impact of such genes on the reprogramming efficiency. In conclusion, we show that enriched mutations in curated onco- and tumor suppressor genes may account for an increased tumor risk and impact the clinical value of patient-derived hiPSCs.

Zfh-2 facilitates Notch-induced apoptosis in the CNS and appendages of Drosophila melanogaster

Apoptosis is a fundamental remodeling process for most tissues during development. In this manuscript we examine a pro-apoptotic function for the Drosophila DNA binding protein Zfh-2 during development of the central nervous system (CNS) and appendages. In the CNS we find that a loss-of-function zfh-2 allele gives an overall reduction of apoptotic cells in the CNS, and an altered pattern of expression for the axonal markers 22C10 and FasII. This same loss-of-function zfh-2 allele causes specific cells in the NB7-3 lineage of the CNS that would normally undergo apoptosis to be inappropriately maintained, whereas a gain-of-function zfh-2 allele has the opposite effect, resulting in a loss of normal NB 7-3 progeny. We also demonstrate that Zfh-2 and Hunchback reciprocally repress each other's gene expression which limits apoptosis to later born progeny of the NB7-3 lineage. Apoptosis is also required for proper segmentation of the fly appendages. We find that Zfh-2 co-localizes with apoptotic cells in the folds of the imaginal discs and presumptive cuticular joints. A reduction of Zfh-2 levels with RNAi inhibits expression of the pro-apoptotic gene reaper, and produces abnormal joints in the leg, antenna and haltere. Apoptosis has previously been shown to be activated by Notch signaling in both the NB7-3 CNS lineage and the appendage joints. Our results indicate that Zfh-2 facilitates Notch-induced apoptosis in these structures.

Whole-exome mutational landscape of neuroendocrine carcinomas of the gallbladder

Neuroendocrine carcinoma (NEC) of the gallbladder (GB-NEC) is a rare but extremely malignant subtype of gallbladder cancer (GBC). The genetic and molecular signatures of GB-NEC are poorly understood; thus, molecular targeting is currently unavailable. In the present study, we applied whole-exome sequencing (WES) technology to detect gene mutations and predicted somatic single-nucleotide variants (SNVs) in 15 cases of GB-NEC and 22 cases of general GBC. In 15 GB-NECs, the C > T mutation was predominant among the 6 types of SNVs. TP53 showed the highest mutation frequency (73%, 11/15). Compared with neuroendocrine carcinomas of other organs, significantly mutated genes (SMGs) in GB-NECs were more similar to those in pulmonary large-cell neuroendocrine carcinomas (LCNECs), with driver roles for TP53 and RB1. In the COSMIC database of cancer-related genes, 211 genes were mutated. Strikingly, RB1 (4/15, 27%) and NAB2 (3/15, 20%) mutations were found specifically in GB-NECs; in contrast, mutations in 29 genes, including ERBB2 and ERBB3, were identified exclusively in GBC. Mutations in RB1 and NAB2 were significantly related to downregulation of the RB1 and NAB2 proteins, respectively, according to immunohistochemical (IHC) data (p values = 0.0453 and 0.0303). Clinically actionable genes indicated 23 mutated genes, including ALK, BRCA1, and BRCA2. In addition, potential somatic SNVs predicted by ISOWN and SomVarIUS constituted 6 primary COSMIC mutation signatures (1, 3, 30, 6, 7, and 13) in GB-NEC. Genes carrying somatic SNVs were enriched mainly in oncogenic signaling pathways involving the Notch, WNT, Hippo, and RTK-RAS pathways. In summary, we have systematically identified the mutation landscape of GB-NEC, and these findings may provide mechanistic insights into the specific pathogenesis of this deadly disease.

Gene networks and transcriptional regulators associated with liver cancer development and progression

BACKGROUND

Treatment options for hepatocellular carcinoma (HCC) are limited, and overall survival is poor. Despite the high frequency of this malignoma, its basic disease mechanisms are poorly understood. Therefore, the aim of this study was to use different methodological approaches and combine the results to improve our knowledge on the development and progression of HCC.

METHODS

Twenty-three HCC samples were characterized by histological, morphometric and cytogenetic analyses, as well as comparative genomic hybridization (aCGH) and genome-wide gene expression followed by a bioinformatic search for potential transcriptional regulators and master regulatory molecules of gene networks.

RESULTS

Histological evaluation revealed low, intermediate and high-grade HCCs, and gene expression analysis split them into two main sets: GE1-HCC and GE2-HCC, with a low and high proliferation gene expression signature, respectively. Array-based comparative genomic hybridization demonstrated a high level of chromosomal instability, with recurrent chromosomal gains of 1q, 6p, 7q, 8q, 11q, 17q, 19p/q and 20q in both HCC groups and losses of 1p, 4q, 6q, 13q and 18q characteristic for GE2-HCC. Gene expression and bioinformatics analyses revealed that different genes and gene regulatory networks underlie the distinct biological features observed in GE1-HCC and GE2-HCC. Besides previously reported dysregulated genes, the current study identified new candidate genes with a putative role in liver cancer, e.g. C1orf35, PAFAH1B3, ZNF219 and others.

CONCLUSION

Analysis of our findings, in accordance with the available published data, argues in favour of the notion that the activated E2F1 signalling pathway, which can be responsible for both inappropriate cell proliferation and initial chromosomal instability, plays a pivotal role in HCC development and progression. A dedifferentiation switch that manifests in exaggerated gene expression changes might be due to turning on transcriptional co-regulators with broad impact on gene expression, e.g. POU2F1 (OCT1) and NFY, as a response to accumulating cell stress during malignant development. Our findings point towards the necessity of different approaches for the treatment of HCC forms with low and high proliferation signatures and provide new candidates for developing appropriate HCC therapies.

Genomic profiles and tumor immune microenvironment of primary lung carcinoma and brain oligo-metastasis

Brain metastasis (BM) is a common malignant event in lung cancer. Here, we recruited 33 lung cancer patients with brain oligo-metastasis to explore the genomic features and tumor immune microenvironment (TIME) of the lung and BM independently. For genomic profiling, targeted sequencing was performed. We found that high-frequent ZFHX3 occurred in the lung (40%) and brain tumor (28%), which might relate to brain metastasis event; the vast majority of patients had lesions-shared mutations in primary tumor and BM, confirming the common clonal events; and EGFR was the most frequently clonal gene in both lung and BM, indicating its driver capability. To characterize TIME status, we also sequenced the T cell receptor (TCR) repertoires and performed immunohistochemistry (IHC) on CD8+ tumor-infiltrating lymphocytes (TILs) and PD-L1 expression in 28 patients who had paired samples. Through the comparison, the TCR clonality of BM was higher than lung tumor, indicating the distinct pattern of the stronger oligoclonal T cell expansion in BM; the primary tumor had a higher TMB than oligo-BM (13.9 vs 8.7 mutations, p = 0.019); CD8 + TILs of BM were significantly lower than lung tumor (10% vs 30%, p = 0.015), revealing the lower level of cytotoxic T cell infiltration; BM showed statistically equivalent level of PD-L1 compared with lung tumor (p = 0.722). We further investigated the potential biomarkers associated with overall survival (OS) after brain surgery. We found that higher TCR clonality was related to prolonged OS in EGFR-treated patients (HR 0.175, p < 0.001) but the worse outcomes in non-EGFR-treated (HR 2.623, p = 0.034). More CD8+ TILs were an independently positive indicator for OS, in EGFR-treated (HR 0.160, p = 0.001) and non-EGFR-treated patients (HR 0.308, p = 0.009). These findings provide a meaningful molecular and clinical understanding of lung carcinoma and brain oligo-metastasis.

LncRNA H19 Suppresses Osteosarcomagenesis by Regulating snoRNAs and DNA Repair Protein Complexes

Osteosarcoma is one of the most frequent common primary malignant tumors in childhood and adolescence. Long non-coding RNAs (lncRNAs) have been reported to regulate the initiation and progression of tumors. However, the exact molecular mechanisms involving lncRNA in osteosarcomagenesis remain largely unknown. Li-Fraumeni syndrome (LFS) is a familial cancer syndrome caused by germline p53 mutation. We investigated the tumor suppressor function of lncRNA H19 in LFS-associated osteosarcoma. Analyzing H19-induced transcriptome alterations in LFS induced pluripotent stem cell (iPSC)-derived osteoblasts, we unexpectedly discovered a large group of snoRNAs whose expression was significantly affected by H19. We identified SNORA7A among the H19-suppressed snoRNAs. SNORA7A restoration impairs H19-mediated osteogenesis and tumor suppression, indicating an oncogenic role of SNORA7A. TCGA analysis indicated that SNORA7A expression is associated with activation of oncogenic signaling and poor survival in cancer patients. Using an optimized streptavidin-binding RNA aptamer designed from H19 lncRNA, we revealed that H19-tethered protein complexes include proteins critical for DNA damage response and repair, confirming H19's tumor suppressor role. In summary, our findings demonstrate a critical role of H19-modulated SNORA7A expression in LFS-associated osteosarcomas.

ZFHX3 Promotes the Proliferation and Tumor Growth of ER-Positive Breast Cancer Cells Likely by Enhancing Stem-Like Features and MYC and TBX3 Transcription

Simple Summary

Breast cancer is a common malignancy, but the understanding of its cellular and molecular mechanisms is limited. The ZFHX3 transcription factor regulates mammary epithelial cells’ proliferation and differentiation by interacting with estrogen and progesterone receptors. Both these receptors play crucial roles in breast cancer development, but whether ZFHX3 also impacts breast cancer is unknown. In this study, the authors aim to determine if ZFHX3 promotes breast cancer cells’ proliferation and tumor growth and explore the underlying cellular and molecular mechanisms. Higher ZFHX3 expression is associated with worse patient survival in breast cancer, ZFHX3 promotes the proliferation and tumor growth of breast cancer cells, and several breast cancer stem cell factors appear to be involved in the role of ZFHX3 in breast cancer growth. The findings suggest that ZFHX3 is a novel oncogenic molecule promoting breast cancer development. Such a molecule could provide novel opportunities for the treatment of breast cancer.

Abstract

Breast cancer is a common malignancy, but the understanding of its cellular and molecular mechanisms is limited. ZFHX3, a transcription factor with many homeodomains and zinc fingers, suppresses prostatic carcinogenesis but promotes tumor growth of liver cancer cells. ZFHX3 regulates mammary epithelial cells’ proliferation and differentiation by interacting with estrogen and progesterone receptors, potent breast cancer regulators. However, whether ZFHX3 plays a role in breast carcinogenesis is unknown. Here, we found that ZFHX3 promoted the proliferation and tumor growth of breast cancer cells in culture and nude mice; and higher expression of ZFHX3 in human breast cancer specimens was associated with poorer prognosis. The knockdown of ZFHX3 in ZFHX3-high MCF-7 cells decreased, and ZFHX3 overexpression in ZFHX3-low T-47D cells increased the proportion of breast cancer stem cells (BCSCs) defined by mammosphere formation and the expression of CD44, CD24, and/or aldehyde dehydrogenase 1. Among several transcription factors that have been implicated in BCSCs, MYC and TBX3 were transcriptionally activated by ZFHX3 via promoter binding, as demonstrated by luciferase-reporter and ChIP assays. These findings suggest that ZFHX3 promotes breast cancer cells’ proliferation and tumor growth likely by enhancing BCSC features and upregulating MYC, TBX3, and others.

S100A2 promotes glycolysis and proliferation via GLUT1 regulation in colorectal cancer

The deregulation of S100A2 has been implicated in the pathogenesis of several types of cancers. However, the molecular mechanisms underlying the protumorigenic capacities of S100A2 have not been fully elucidated. Here, we demonstrated the molecular mechanisms underlying the roles of S100A2 in glycolysis reprogramming and proliferation of colorectal cancer (CRC) cells. The results indicated that S100A2 overexpression raises glucose metabolism and proliferation. Mechanistically, S100A2 activated the PI3K/AKT signaling pathway, upregulated GLUT1 expression, induced glycolytic reprogramming, and consequently increased proliferation. Clinical data showed significantly increased S100A2 levels in CRC tissues and the Oncomine database. In addition, analysis revealed a positive correlation between S100A2 and GLUT1 mRNA expression in CRC tissues. Together, these results demonstrate that the S100A2/GLUT1 axis can promote the progression of CRC by modulating glycolytic reprogramming. Our results further suggest that targeting S100A2 could present a promising therapeutic avenue for the prevention of colorectal cancer progression.

Efficient mutation screening for cervical cancers from circulating tumor DNA in blood

Background

Early diagnosis and continuous monitoring are necessary for an efficient management of cervical cancers (CC). Liquid biopsy, such as detecting circulating tumor DNA (ctDNA) from blood, is a simple, non-invasive method for testing and monitoring cancer markers. However, tumor-specific alterations in ctDNA have not been extensively investigated or compared to other circulating biomarkers in the diagnosis and monitoring of the CC. Therfore, Next-generation sequencing (NGS) analysis with blood samples can be a new approach for highly accurate diagnosis and monitoring of the CC.

Method

Using a bioinformatics approach, we designed a panel of 24 genes associated with CC to detect and characterize patterns of somatic single-nucleotide variations, indels, and copy number variations. Our NGS CC panel covers most of the genes in The Cancer Genome Atlas (TCGA) as well as additional cancer driver and tumor suppressor genes. We profiled the variants in ctDNA from 24 CC patients who were being treated with systemic chemotherapy and local radiotherapy at the Jeonbuk National University Hospital, Korea.

Result

Eighteen out of 24 genes in our NGS CC panel had mutations across the 24 CC patients, including somatic alterations of mutated genes (ZFHX3–83%, KMT2C-79%, KMT2D-79%, NSD1–67%, ATM-38% and RNF213–27%). We demonstrated that the RNF213 mutation could be used potentially used as a monitoring marker for response to chemo- and radiotherapy.

Conclusion

We developed our NGS CC panel and demostrated that our NGS panel can be useful for the diagnosis and monitoring of the CC, since the panel detected the common somatic variations in CC patients and we observed how these genetic variations change according to the treatment pattern of the patient.

Establishment of a novel cell cycle-related prognostic signature predicting prognosis in patients with endometrial cancer

Background

Endometrial cancer (EnCa) ranks fourth in menace within women’s malignant tumors. Large numbers of studies have proven that functional genes can change the process of tumors by regulating the cell cycle, thereby achieving the goal of targeted therapy.

Methods

The transcriptional data of EnCa samples obtained from the TCGA database was analyzed. A battery of bioinformatics strategies, which included GSEA, Cox and LASSO regression analysis, establishment of a prognostic signature and a nomogram for overall survival (OS) assessment. The GEPIA and CPTAC analysis were applied to validate the dysregulation of hub genes. For mutation analysis, the “maftools” package was used.

Results

GSEA identified that cell cycle was the most associated pathway to EnCa. Five cell cycle-related genes including HMGB3, EZH2, NOTCH2, UCK2 and ODF2 were identified as prognosis-related genes to build a prognostic signature. Based on this model, the EnCa patients could be divided into low- and high-risk groups, and patients with high-risk score exhibited poorer OS. Time-dependent ROC and Cox regression analyses revealed that the 5-gene signature could predict EnCa prognosis exactly and independently. GEPIA and CPTAC validation exhibited that these genes were notably dysregulated between EnCa and normal tissues. Lower mutation rates of PTEN, TTN, ARID1A, and etc. were found in samples with high-risk score compared with that with low-risk score. GSEA analysis suggested that the samples of the low- and high-risk groups were concentrated on various pathways, which accounted for the different oncogenic mechanisms in patients in two groups.

Conclusion

The current research construct a 5-gene signature to evaluate prognosis of EnCa patients, which may innovative clinical application of prognostic assessment.

ZFHX3 mutation as a protective biomarker for immune checkpoint blockade in non-small cell lung cancer

To date, immunotherapy has opened a new chapter in the treatment of lung cancer. Precise biomarkers can help to screen subpopulations of lung cancer to provide the best treatment. Multiple studies suggest that specific gene mutations may be predictive markers in guiding non-small cell lung cancer (NSCLC) immune checkpoint inhibitor (ICI) treatment. A published immunotherapy cohort with mutational and survival data for 350 NSCLC patients was used. First, the mutational data of the immunotherapy cohort were used to identify gene mutations related to the prognosis of ICI therapy. The immunotherapy cohort and TCGA-NSCLC cohort were further studied to elucidate the relationships between specific gene mutations and tumor immunogenicity, antitumor immune response capabilities, and immune cell and mutation counts in the DNA damage response (DDR) pathway. In the immunotherapy cohort (N = 350), ZFHX3 mutations were an independent predictive biomarker for NSCLC patients receiving ICI treatment. Significant differences were observed between ZFHX3-mutant (ZFHX3-MT) and ZFHX3-wild type (ZFHX3-WT) patients regarding the overall survival (OS) time (P < 0.001, HR = 0.26, 95% Cl 0.17-0.41). ZFHX3-MT is significantly associated with higher tumor mutation burden (TMB) and neoantigen load (NAL), and ZFHX3-MT positively correlates with known immunotherapy response biomarkers, including T-cell infiltration, immune-related gene expression, and mutation counts in the DDR pathway in NSCLC. ZFHX3-MT is closely related to longer OS in NSCLC patients treated with ICIs, suggesting that ZFHX3 mutations be used as a novel predictive marker in guiding NSCLC ICI treatment.

The transcription factor ZFHX3 is crucial for the angiogenic function of hypoxia-inducible factor 1α in liver cancer cells

Angiogenesis is a hallmark of tumorigenesis, and hepatocellular carcinoma (HCC) is hypervascular and therefore very dependent on angiogenesis for tumor development and progression. Findings from previous studies suggest that in HCC cells, hypoxia-induced factor 1α (HIF1A) and zinc finger homeobox 3 (ZFHX3) transcription factors functionally interact in the regulation of genes in HCC cells. Here, we report that hypoxia increases the transcription of the ZFHX3 gene and enhances the binding of HIF1A to the ZFHX3 promoter in the HCC cell lines HepG2 and Huh-7. Moreover, ZFHX3, in turn, physically associated with and was functionally indispensable for HIF1A to exert its angiogenic activity, as indicated by in vitro migration and tube formation assays of human umbilical vein endothelial cells (HUVECs) and microvessel formation in xenograft tumors of HCC cells. Mechanistically, ZFHX3 was required for HIF1A to transcriptionally activate the vascular endothelial growth factor A (VEGFA) gene by binding to its promoter. Functionally, down-regulation of ZFHX3 in HCC cells slowed their tumor growth, and addition of VEGFA to conditioned medium from ZFHX3-silenced HCC cells partially rescued the inhibitory effect of this medium on HUVEC tube formation. In human HCC, ZFHX3 expression was up-regulated, and this up-regulation correlated with both HIF1A up-regulation and worse patient survival, confirming a functional association between ZFHX3 and HIF1A in human HCC. We conclude that ZFHX3 is an angiogenic transcription factor that is integral to the HIF1A/VEGFA signaling axis in HCC cells.

SUMOylation of the transcription factor ZFHX3 at Lys-2806 requires SAE1, UBC9, and PIAS2 and enhances its stability and function in cell proliferation

SUMOylation is a posttranslational modification (PTM) at a lysine residue and is crucial for the proper functions of many proteins, particularly of transcription factors, in various biological processes. Zinc finger homeobox 3 (ZFHX3), also known as AT motif-binding factor 1 (ATBF1), is a large transcription factor that is active in multiple pathological processes, including atrial fibrillation and carcinogenesis, and in circadian regulation and development. We have previously demonstrated that ZFHX3 is SUMOylated at three or more lysine residues. Here, we investigated which enzymes regulate ZFHX3 SUMOylation and whether SUMOylation modulates ZFHX3 stability and function. We found that SUMO1, SUMO2, and SUMO3 each are conjugated to ZFHX3. Multiple lysine residues in ZFHX3 were SUMOylated, but Lys-2806 was the major SUMOylation site, and we also found that it is highly conserved among ZFHX3 orthologs from different animal species. Using molecular analyses, we identified the enzymes that mediate ZFHX3 SUMOylation; these included SUMO1-activating enzyme subunit 1 (SAE1), an E1-activating enzyme; SUMO-conjugating enzyme UBC9 (UBC9), an E2-conjugating enzyme; and protein inhibitor of activated STAT2 (PIAS2), an E3 ligase. Multiple analyses established that both SUMO-specific peptidase 1 (SENP1) and SENP2 deSUMOylate ZFHX3. SUMOylation at Lys-2806 enhanced ZFHX3 stability by interfering with its ubiquitination and proteasomal degradation. Functionally, Lys-2806 SUMOylation enabled ZFHX3-mediated cell proliferation and xenograft tumor growth of the MDA-MB-231 breast cancer cell line. These findings reveal the enzymes involved in, and the functional consequences of, ZFHX3 SUMOylation, insights that may help shed light on ZFHX3's roles in various cellular and pathophysiological processes.

A Systematic Study of the Impact of Estrogens and Selective Estrogen Receptor Modulators on Prostate Cancer Cell Proliferation

The estrogen signaling pathway has been reported to modulate prostate cancer (PCa) progression through the activity of estrogen receptors α and β (ERα and ERβ). Given that selective estrogen receptor modulators (SERMs) are used to treat breast cancer, ERs have been proposed as attractive therapeutic targets in PCa. However, many inconsistencies regarding the expression of ERs and the efficacy of SERMs for PCa treatment exist, notably due to the use of ERβ antibodies lacking specificity and treatments with high SERM concentrations leading to off-target effects. To end this confusion, our objective was to study the impact of estrogenic and anti-estrogenic ligands in well-studied in vitro PCa models with appropriate controls, dosages, and ER subtype-specific antibodies. When using physiologically relevant concentrations of nine estrogenic/anti-estrogenic compounds, including five SERMs, we observed no significant modulation of PCa cell proliferation. Using RNA-seq and validated antibodies, we demonstrate that these PCa models do not express ERs. In contrast, RNA-seq from PCa samples from patients have detectable expression of ERα. Overall, our study reveals that commonly used PCa models are inappropriate to study ERs and indicate that usage of alternative models is essential to properly assess the roles of the estrogen signaling pathway in PCa.

Combination of Rapamycin and MK-2206 Induced Cell Death via Autophagy and Necroptosis in MYCN-Amplified Neuroblastoma Cell Lines

Neuroblastoma (NB) is the most common pediatric malignant extracranial solid tumor. Despite multi-modality therapies, the emergence of drug resistance is an obstacle in the treatment of high-risk NB patients (with MYCN amplification). In our previous study, we found that rapamycin and MK-2206 synergistically induced cell death in MYCN-amplified cell lines but the mechanisms remained unclear. In our present study, either 3-MA or necroatatin-1 blocked the cell death induced by rapamycin and MK-2206, but z-VAD-fmk did not block this cell death. The expressions of autophagy markers (ATG5, ATG7, Beclin-1, LC3 B) and the necroptosis marker RIPK3 increased and another necroptosis marker RIPK1 decreased after the combination treatment of rapamycin and MK-2206, and were accompanied by the morphological characteristics of autophagy and necroptosis. In NB xenograft tumor tissues, the expressions of autophagy and necroptosis markers were consistent with observations in vitro. These data suggested that autophagy and necroptosis contributed to the cell death induced by rapamycin and MK-2206 in NB cells. To understand the role of MYCN in this process, MYCN expression was downregulated in MYCN-amplified cell lines (NGP, BE2) using siRNAs and was upregulated in MYCN non-amplified cell lines (AS, SY5Y) using plasmid. We found the cell death induced by rapamycin and MK-2206 was MYCN-dependent. We also found that the metabolic activity in NB cells was correlated with the expression level of MYCN. This study delineates the role of MYCN in the cell death induced by combination treatment of rapamycin and MK-2206 in MYCN-amplified NB cells.

6-Thioguanine and Its Analogs Promote Apoptosis of Castration-Resistant Prostate Cancer Cells in a BRCA2-Dependent Manner

Background: Mutations in the oncosuppressor gene BReast CAncer susceptibility gene 2 (BRCA2) predispose to aggressive forms of prostate cancer which show poor response to taxane-based therapy, the standard treatment for castration-resistant, aggressive prostate cancer. Herein, we addressed the question whether changes in BRCA2 expression, a potential surrogate marker for BRCA2 activity, may affect the response of castration-resistant prostate cancer cells to 6-thioguanine (6-TG), a thiopurine used in the treatment of haematological malignancies. Methods: Yeast, normal prostate cells and castration-resistant prostate cancer cells were treated with 6-TG or its analogues, in presence or absence of paclitaxel, or with olaparib, a poly-(ADP-ribose) polymerase (PARP) inhibitor currently in clinical trials for treatment of metastatic castration-resistant prostate cancer, and cell proliferation, apoptosis and androgen receptor (AR) levels were measured. Results: 6-TG inhibited cell proliferation in yeast, normal and castration-resistant prostate cancer cells but promoted apoptosis only in cancer cells. Suppression of BRCA2 expression by siRNA or shRNA increased the sensitivity to 6-TG- and olaparib-induced apoptosis but did not affect cancer cell response to taxane. Intriguingly, 6-TG reduced AR expression levels independently on BRCA2 expression. Instead, olaparib decreased AR levels only in BRCA2-knockdown prostate cancer cells. Notably, overexpression of BRCA2 resulted in resistance of castration-resistant prostate cancer cells to 6-TG-, taxane- and olaparib-based treatment but promoted sensitivity to apoptosis induced by 2-amino-6-bromopurine and 2,6-dithiopurine, two 6-TG analogues. Conclusions: Our results provide a pre-clinical rationale for the use of 6-TG in the treatment of BRCA2-deficient castration-resistant prostate cancers, and of certain 6-TG analogues for treatment of BRCA2-proficient prostate cancers.

Zfhx3 is essential for progesterone/progesterone receptor signaling to drive ductal side-branching and alveologenesis in mouse mammary glands

Progesterone (Pg)/progesterone receptor (PR) signaling drives mammary gland side-branching and alveologenesis, but the mechanisms through which Pg/PR signaling functions remain to be clarified. Using in vitro and in vivo models and histological and molecular analyses, we determined the role of Zfhx3 transcription factor in mammary gland development driven by Pg/PR signaling. Postnatal deletion of Zfhx3 in mouse mammary epithelial cells attenuated side-branching morphogenesis and alveologenesis. These effects were undetectable in the absence of Pg/PR signaling. During the estrus cycle, Zfhx3 expression corresponded to that of Pg, being at the highest level at the diestrus stage; Zfhx3 deletion inhibited mammary gland branching more potently at diestrus than estrus stage. Loss of Zfhx3 not only attenuated the expansion of stem/progenitor cells driven by Pg/PR signaling, but also impaired the function of Pg/PR signaling in the transcriptional activation of multiple genes. In addition, Pg/PR signaling significantly expanded PR- and Zfhx3-positive epithelial cells, and induced the physical association of ZFHX3 with PR. These findings establish Zfhx3 as an integral transcription factor of Pg/PR signaling in driving side-branching and alveologenesis during mammary gland development.