CARM1 drives triple-negative breast cancer progression by coordinating with HIF1A.

Coactivator-associated arginine methyltransferase 1 (CARM1) promotes the development and metastasis of estrogen receptor alpha (ERα)-positive breast cancer. The function of CARM1 in triple-negative breast cancer (TNBC) is still unclear and requires further exploration. Here, we report that CARM1 promotes proliferation, epithelial-mesenchymal transition (EMT), and stemness in TNBC. CARM1 is upregulated in multiple cancers and its expression correlates with breast cancer progression. Genome-wide analysis of CARM1 showed that CARM1 is recruited by hypoxia-inducible factor 1 subunit alpha (HIF1A) and occupy the promoters of CDK4, Cyclin D1, ß-catenin, HIF1A, MALAT1, and SIX1 critically involved in cell cycle, HIF-1 signaling pathway, Wnt signaling pathway, VEGF signaling pathway, thereby modulating the proliferation and invasion of TNBC cells. We demonstrated that CARM1 is physically associated with and directly interacts with HIF1A. Moreover, we found that ellagic acid, an inhibitor of CARM1, can suppress the proliferation and metastasis of TNBC by directly inhibiting CDK4 expression. Our research has determined the molecular basis of CARM1 carcinogenesis in TNBC and its effective natural inhibitor, which may provide new ideas and drugs for cancer therapy.

JARID2 coordinates with the NuRD complex to facilitate breast tumorigenesis through response to adipocyte-derived leptin.

BACKGROUND: Proteins containing the Jumonji C (JmjC) domain participated in tumorigenesis and cancer progression. However, the mechanisms underlying this effect are still poorly understood. Our objective was to investigate the role of Jumonji and the AT-rich interaction domain-containing 2 (JARID2) - a JmjC family protein - in breast cancer, as well as its latent association with obesity. METHODS: Immunohistochemistry, The Cancer Genome Atlas, Gene Expression Omnibus, and other databases were used to analyze the expression of JARID2 in breast cancer cells. Growth curve, 5-ethynyl-2-deoxyuridine (EdU), colony formation, and cell invasion experiments were used to detect whether JARID2 affected breast cancer cell proliferation and invasion. Spheroidization-based experiments and xenotumor transplantation in NOD/SCID mice were used to examine the association between JARID2 and breast cancer stemness. RNA-sequencing, Kyoto Encyclopedia of Genes and Genomes, and Gene Set Enrichment Analysis were used to identify the cell processes in which JARID2 participates. Immunoaffinity purification and silver staining mass spectrometry were conducted to search for proteins that might interact with JARID2. The results were further verified using co-immunoprecipitation and glutathione S-transferase (GST) pull-down experiments. Using chromatin immunoprecipitation (ChIP) sequencing, we sought the target genes that JARID2 and metastasis-associated protein 1 (MTA1) jointly regulated; the results were validated by ChIP-PCR, quantitative ChIP (qChIP) and ChIP-reChIP assays. A coculture experiment was used to explore the interactions between breast cancer cells and adipocytes. RESULTS: In this study, we found that JARID2 was highly expressed in multiple types of cancer including breast cancer. JARID2 promoted glycolysis, lipid metabolism, proliferation, invasion, and stemness of breast cancer cells. Furthermore, JARID2 physically interacted with the nucleosome remodeling and deacetylase (NuRD) complex, transcriptionally repressing a series of tumor suppressor genes such as BRCA2 DNA repair associated (BRCA2), RB transcriptional corepressor 1 (RB1), and inositol polyphosphate-4-phosphatase type II B (INPP4B). Additionally, JARID2 expression was regulated by the obesity-associated adipokine leptin via Janus kinase 2/signal transducer and activator of transcription 3 (JAK2/STAT3) pathway in the breast cancer microenvironment. Analysis of various online databases also indicated that JARID2/MTA1 was associated with a poor prognosis of breast cancer. CONCLUSION: Our data indicated that JARID2 promoted breast tumorigenesis and development, confirming JARID2 as a target for cancer treatment.

Prognostic analysis of patients with breast cancer based on tumor mutational burden and DNA damage repair genes.

Background: Breast cancer has a high tumor-specific death rate and poor prognosis. In this study, we aimed to provide a basis for the prognostic risk in patients with breast cancer using significant gene sets selected by analyzing tumor mutational burden (TMB) and DNA damage repair (DDR). Methods: Breast cancer genomic and transcriptomic data were obtained from The Cancer Genome Atlas (TCGA). Breast cancer samples were dichotomized into high- and low-TMB groups according to TMB values. Differentially expressed DDR genes between high- and low-TMB groups were incorporated into univariate and multivariate cox regression model to build prognosis model. Performance of the prognosis model was validated in an independently new GEO dataset and evaluated by time-dependent ROC curves. Results: Between high- and low-TMB groups, there were 6,424 differentially expressed genes, including 67 DDR genes. Ten genes associated with prognosis were selected by univariate cox regression analysis, among which seven genes constituted a panel to predict breast cancer prognosis. The seven-gene prognostic model, as well as the gene copy numbers are closely associated with tumor-infiltrating immune cells. Conclusion: We established a seven-gene prognostic model comprising MDC1, PARP3, PSMB1, PSMB9, PSMD2, PSMD7, and PSMD14 genes, which provides a basis for further exploration of a population-based prediction of prognosis and immunotherapy response in patients with breast cancer.

Identifying key transcription factors and immune infiltration in non-small-cell lung cancer using weighted correlation network and Cox regression analyses.

Introduction: Lung cancer is one of the most common cancers and a significant cause of cancer-related deaths. Non-small cell lung cancer (NSCLC) accounts for about 85% of all lung cancer cases. Therefore, it is crucial to identify effective diagnostic and therapeutic methods. In addition, transcription factors are essential for eukaryotic cells to regulate their gene expression, and aberrant expression transcription factors are an important step in the process of oncogenesis in NSCLC. Methods: Differentially expressed transcription factors between NSCLC and normal tissues by analyzing mRNA profiling from The Cancer Genome Atlas (TCGA) database program were identified. Weighted correlation network analysis (WGCNA) and line plot of least absolute shrinkage and selection operator (LASSO) were performed to find prognosis-related transcription factors. The cellular functions of transcription factors were performed by 5-ethynyl-2'-deoxyuridine (EdU) assay, wound healing assay, cell invasion assay in lung cancer cells. Results: We identified 725 differentially expressed transcription factors between NSCLC and normal tissues. Three highly related modules for survival were discovered, and transcription factors highly associated with survival were obtained by using WGCNA. Then line plot of LASSO was applied to screen transcription factors related to prognosis and build a prognostic model. Consequently, SETDB2, SNAI3, SCML4, and ZNF540 were identified as prognosis-related transcription factors and validated in multiple databases. The low expression of these hub genes in NSCLC was associated with poor prognosis. The deletions of both SETDB2 and SNAI3 were found to promote proliferation, invasion, and stemness in lung cancer cells. Furthermore, there were significant differences in the proportions of 22 immune cells between the high- and low-score groups. Discussion: Therefore, our study identified the transcription factors involved in regulating NSCLC, and we constructed a panel for the prediction of prognosis and immune infiltration to inform the clinical application of transcription factor analysis in the prevention and treatment of NSCLC.

The PRMT6/PARP1/CRL4B Complex Regulates the Circadian Clock and Promotes Breast Tumorigenesis.

Circadian rhythms, as physiological systems with self-regulatory functions in living organisms, are controlled by core clock genes and are involved in tumor development. The protein arginine methyltransferase 6 (PRMT6) serves as an oncogene in a myriad of solid tumors, including breast cancer. Hence, the primary aim of the current study is to investigate the molecular mechanisms by which the PRMT6 complex promotes breast cancer progression. The results show that PRMT6, poly(ADP-ribose) polymerase 1 (PARP1), and the cullin 4 B (CUL4B)-Ring E3 ligase (CRL4B) complex interact to form a transcription-repressive complex that co-occupies the core clock gene PER3 promoter. Moreover, genome-wide analysis of PRMT6/PARP1/CUL4B targets identifies a cohort of genes that is principally involved in circadian rhythms. This transcriptional-repression complex promotes the proliferation and metastasis of breast cancer by interfering with circadian rhythm oscillation. Meanwhile, the PARP1 inhibitor Olaparib enhances clock gene expression, thus, reducing breast carcinogenesis, indicating that PARP1 inhibitors have potential antitumor effects in high-PRMT6 expression breast cancer.

Ehm2 transcript variant 1 inhibits breast cancer progression and increases E-cadherin stability.

Ehm2/1, an Ehm2 transcript variant, regulates the cytoskeleton by binding to plasma membrane proteins. However, the role of Ehm2/1 in breast cancer development remains poorly understood. This study shows that, the expression of Ehm2/1 was decreased in breast cancer and that patients with low Ehm2/1 expression had a significantly poorer prognosis than those with high expression of Ehm2/1. Overexpression of Ehm2/1 in MCF-7 breast cancer cells inhibited cell migration and invasion. Ehm2/1 markedly increased the stability and half-life of E-cadherin. Moreover, Ehm2/1 was collocated with E-cadherin in the plasma membrane of MCF-7 cells. Furthermore, downregulation of Ehm2/1 promoted ubiquitination of E-cadherin, whereas overexpression of Ehm2/1 inhibited ubiquitination of E-cadherin. These results suggest that Ehm2/1 could suppress the migration and invasion of breast cancer cells by increasing E-cadherin stability.

RUNX2 recruits the NuRD(MTA1)/CRL4B complex to promote breast cancer progression and bone metastasis.

Runt-related transcription factor 2 (RUNX2) is an osteogenesis-related transcription factor that has emerged as a prominent transcription repressing factor in carcinogenesis. However, the role of RUNX2 in breast cancer metastasis remains poorly understood. Here, we show that RUNX2 recruits the metastasis-associated 1 (MTA1)/NuRD and the Cullin 4B (CUL4B)-Ring E3 ligase (CRL4B) complex to form a transcriptional-repressive complex, which catalyzes the histone deacetylation and ubiquitylation. Genome-wide analysis of the RUNX2/NuRD(MTA1)/CRL4B complex targets identified a cohort of genes including peroxisome proliferator-activated receptor alpha (PPARα) and superoxide dismutase 2 (SOD2), which are critically involved in cell growth, epithelial-to-mesenchymal transition (EMT) and invasion. We demonstrate that the RUNX2/NuRD(MTA1)/CRL4B complex promotes the proliferation, invasion, tumorigenesis, bone metastasis, cancer stemness of breast cancer in vitro and in vivo. Strikingly, RUNX2 expression is upregulated in multiple human carcinomas, including breast cancer. Our study suggests that RUNX2 is a promising potential target for the future treatment strategies of breast cancer.

The hub ten gene-based risk score system using RNA m6A methylation regulator features and tumor immune microenvironment in breast cancer.

BACKGROUND: RNA N6-methyladenosine (m6A) modification is primarily regulated by m6A regulators, which play significant epigenetic regulatory roles in tumorigenesis, tumor development, and tumor immune microenvironment. However, the correlation between m6A regulators and immune cell infiltration in breast cancer remains unclear. METHODS: In this study, m6A modification patterns were evaluated based on 31 m6A modification regulators. m6A clusters were determined by consensus clustering. Immune landscape and immune cell infiltration subgroups were characterized by m6A clusters. Key module and hub genes related to m6A regulators and immune infiltration cells were identified by WGCNA. LASSO algorithm was applied to select prognostic signatures. Multivariate Cox regression analysis was applied to assess the prognostic value of gene signatures. RESULTS: Two distinct m6A clusters were determined based on the expression of 31 m6A modification regulators and characterized by two tumor immune microenvironment (TIME) immune cell infiltration subgroups. Further, a total of 1971 differentially expressed genes between breast cancer patients and healthy controls were screened, nine modules associated with clinical characteristics of breast cancer patients were identified. Later, one key module and 13 hub genes correlated with m6A regulators and immune infiltration cells were identified. LASSO Cox regression analysis selected and constructed a ten-gene prognostic model to build a risk score system for individual breast cancer patient prognosis. The performance of the ten-gene-based risk score system was further validated in an independent dataset with an AUC of 0.659. CONCLUSIONS: This study revealed that m6A modification regulators played a significant role in the TIME regulation of breast cancer. The hub ten gene-based risk score system is valuable in predicting the prognosis of breast cancer patients, which may provide potential significance for breast cancer diagnosis, prognosis, and immunotherapy in the future.

m6A regulators featured by tumor immune microenvironment landscapes and correlated with immunotherapy in non-small cell lung cancer (NSCLC).

Introduction: Recent research has confirmed the critical role that epigenetic factors play in regulating the immune response. Nonetheless, what role m6A methylation modification might play in the immune response of non-small cell lung cancer (NSCLC) remains vague. Methods: Herein, the gene expression, copy number variations (CNVs), and somatic mutations of 31 m6A regulators in NSCLC and adjacent control samples from the GEO and TCGA databases were comprehensively explored. Using consensus clustering, m6A modification patterns were identified. Correlations between m6A modification patterns and immune cell infiltration traits in the tumor immune microenvironment (TME) were systematically analyzed. Differentially expressed genes were verified and screened by random forest and cox regression analysis by comparing different m6A modification patterns. Based on the retained gene panel, a risk model was built, and m6Ascore for each sample was calculated. The function of m6Ascore in NSCLC prognosis, tumor somatic mutations, and chemotherapy/immunotherapy response prediction were evaluated. Results: Consensus clustering classified all NSCLC samples into two m6A clusters (m6A_clusterA and m6A_clusterB) according to the expression levels of 25 m6A regulator genes. Hierarchical clustering further divides the NSCLC samples into two m6A gene clusters: m6AgeneclusterA and m6AgeneclusterB. A panel of 83 genes was screened from the 194 differentially expressed genes between m6A gene clusters. Based on this, a risk score model was established. m6A modification clusters, m6A gene clusters, and m6Ascore calculated from the risk model were able to predict tumor stages, immune cell infiltration, clinical prognosis, and tumor somatic mutations. NSCLC patients with high m6Ascore have poor drug resistance to chemotherapy drugs (Cisplatin and Gemcitabine) and exhibit considerable therapeutic benefits and favorable clinical responses to anti-PD1 or anti-CTLA4 immunotherapy. Discussion: In conclusion, methylation modification patterns mediated by the m6A regulators in individuals play a non-negligible role in prognosis prediction and immunotherapy response, which will facilitate personalized treatment and immunotherapeutic strategies for NSCLC patients in the future.

BCAS3 exhibits oncogenic properties by promoting CRL4A-mediated ubiquitination of p53 in breast cancer.

OBJECTIVES: Breast cancer-amplified sequence 3 (BCAS3) was initially found to be amplified in human breast cancer (BRCA); however, there has been little consensus on the functions of BCAS3 in breast tumours. MATERIALS AND METHODS: We analysed BCAS3 expression in BRCA using bio-information tools. Affinity purification and mass spectrometry were employed to identify BCAS3-associated proteins. GST pull-down and ubiquitination assays were performed to analyse the interaction mechanism between BCAS3/p53 and CUL4A-RING E3 ubiquitin ligase (CRL4A) complex. BCAS3 was knocked down individually or in combination with p53 in MCF-7 cells to further explore the biological functions of the BCAS3/p53 axis. The clinical values of BCAS3 for BRCA progression were evaluated via semiquantitative immunohistochemistry (IHC) analysis and Cox regression. RESULTS: We reported that the expression level of BCAS3 in BRCA was higher than that in adjacent normal tissues. High BCAS3 expression promoted growth, inhibited apoptosis and conferred chemoresistance in breast cancer cells. Mechanistically, BCAS3 overexpression fostered BRCA cell growth by interacting with the CRL4A complex and promoting ubiquitination and proteasomal degradation of p53. Furthermore, BCAS3 could regulate cell growth, apoptosis and chemoresistance through a p53-mediated mechanism. Clinically, BCAS3 overexpression was significantly correlated with a malignant phenotype. Moreover, higher expression of BCAS3 correlates with shorter overall survival (OS) in BRCA. CONCLUSIONS: The functional characterization of BCAS3 offers new insights into the oncogenic properties and chemotherapy resistance in breast cancer.

Snail/PRMT5/NuRD complex contributes to DNA hypermethylation in cervical cancer by TET1 inhibition.

The biological function of PRMT5 remains poorly understood in cervical cancer metastasis. Here, we report that PRMT5 physically associates with the transcription factor Snail and the NuRD(MTA1) complex to form a transcriptional-repressive complex that catalyzes the symmetrical histone dimethylation and deacetylation. This study shows that the Snail/PRMT5/NuRD(MTA1) complex targets genes, such as TET1 and E-cadherin, which are critical for epithelial-mesenchymal transition (EMT). This complex also affects the conversion of 5mC to 5hmC. This study demonstrates that the Snail/PRMT5/NuRD(MTA1) complex promotes the invasion and metastasis of cervical cancer in vitro and in vivo. This study also shows that PRMT5 expression is upregulated in cervical cancer and various human cancers, and the PRMT5 inhibitor EPZ015666 suppresses EMT and the invasion potential of cervical cancer cells by disinhibiting the expression of TET1 and increasing 5hmC, suggesting that PRMT5 is a potential target for cancer therapy.

The histone deacetylase inhibitor PCI-24781 impairs calcium influx and inhibits proliferation and metastasis in breast cancer.

Histone deacetylases (HDACs) are involved in key cellular processes and have been implicated in cancer. As such, compounds that target HDACs or drugs that target epigenetic markers may be potential candidates for cancer therapy. This study was therefore aimed to identify a potential epidrug with low toxicity and high efficiency as anti-tumor agents. Methods: We first screened an epigenetic small molecule inhibitor library to screen for an epidrug for breast cancer. The candidate was identified as PCI-24781 and was characterized for half maximal inhibitory concentration (IC50), for specificity to breast cancer cells, and for effects on carcinogenesis and metastatic properties of breast cancer cell lines in vitro. A series of in silico and in vitro analyses were further performed of PCI-24781 to identify and understand its target. Results: Screening of an epigenetic inhibitor library in MDA-MB-231 cells, a malignant cancer cell line, showed that PCI-24781 is a potential anti-tumor drug specific to breast cancer. Ca2+ related pathways were identified as a potential target of PCI-24781. Further analyses showed that PCI-24781 inhibited Gαq-PLCß3-mediated calcium signaling by activating the expression of regulator of G-protein signaling 2 (RGS2) to reduce cell proliferation, metastasis, and differentiation, resulting in cell death in breast cancer. In addition, RGS2 depletion reversed anti-tumor effect and inhibition of calcium influx induced by PCI-24781 treatment in breast cancer cells. Conclusions: We have demonstrated that PCI-24781 is an effective anti-tumor therapeutic agent that targets calcium signaling by activating RGS2. This study also provides a novel perspective into the use of HDAC inhibitors for cancer therapy.

Identification of Key Transcription Factors and Immune Infiltration Patterns Associated With Breast Cancer Prognosis Using WGCNA and Cox Regression Analysis.

Breast cancer is the most frequently diagnosed cancer and the second leading cause of cancer death among women worldwide. Therefore, the need for effective breast cancer treatment is urgent. Transcription factors (TFs) directly participate in gene transcription, and their dysregulation plays a key role in breast cancer. Our study identified 459 differentially expressed TFs between tumor and normal samples from The Cancer Genome Atlas database. Based on gene expression analysis and weighted gene co-expression network analysis, the co-expression yellow module was found to be integral for breast cancer progression. A total of 121 genes in the yellow module were used for function enrichment. To further confirm prognosis-related TFs, COX regression and LASSO analyses were performed; consequently, a prognostic risk model was constructed, and its validity was verified. Ten prognosis-related TFs were identified according to their expression profile, survival probability, and target genes. COPS5, HDAC2, and NONO were recognized as hub TFs in breast cancer. These TFs were highly expressed in human breast cancer cell lines and clinical breast cancer samples; this result was consistent with the information from multiple databases. Immune infiltration analysis revealed that the proportions of resting dendritic and mast cells were greater in the low-risk group than those in the high-risk group. Thus, in this study, we identified three hub biomarkers related to breast cancer prognosis. The results provide a framework for the co-expression of TF modules and immune infiltration in breast cancer.

Identification of key modules and genes associated with breast cancer prognosis using WGCNA and ceRNA network analysis.

Breast cancer is one of the leading causes of cancer-associated mortality in women worldwide and has become a major public health problem. Although the definitive cause of breast cancer is not known, many genes sensitive to breast cancer have been detected using advanced technologies. Our study identified 3301 differentially expressed lncRNAs and mRNAs between tumor and normal samples from The Cancer Genome Atlas database. Based on the gene expression analysis and clinical traits as well as weighted gene co-expression network analysis, the co-expression Brown module was found to be key for breast cancer prognosis. A total of 453 genes in the Brown module were used for functional enrichment, protein-protein interaction analysis, lncRNA-miRNA-mRNA ceRNA network, and lncRNA-RNA binding protein-mRNA network construction. GRM4, SSTR2, PARD6B, PRR15, COX6C, and lncRNA DSCAM-AS1 were the hub genes according to protein-protein interaction, lncRNA-miRNA-mRNA and lncRNA-RNA binding protein-mRNA network. Their high expression was found to be correlated with breast cancer development, according to multiple databases. In conclusion, this study provides a framework of the co-expression gene modules of breast cancer and identifies several important biomarkers in breast cancer development and prognosis.

PHF20L1 as a H3K27me2 reader coordinates with transcriptional repressors to promote breast tumorigenesis.

TUDOR domain-containing proteins (TDRDs) are chiefly responsible for recognizing methyl-lysine/arginine residue. However, how TDRD dysregulation contributes to breast tumorigenesis is poorly understood. Here, we report that TUDOR domain-containing PHF20L1 as a H3K27me2 reader exerts transcriptional repression by recruiting polycomb repressive complex 2 (PRC2) and Mi-2/nucleosome remodeling and deacetylase (NuRD) complex, linking PRC2-mediated methylation and NuRD-mediated deacetylation of H3K27. Furthermore, PHF20L1 was found to serve as a potential MYC and hypoxia-driven oncogene, promoting glycolysis, proliferation, and metastasis of breast cancer cells by directly inhibiting tumor suppressors such as HIC1, KISS1, and BRCA1. PHF20L1 expression was also strongly correlated with higher histologic grades of breast cancer and markedly up-regulated in several cancers. Meanwhile, Phf20l1 deletion not only induces growth retardation and mammary ductal outgrowth delay but also inhibits tumorigenesis in vivo. Our data indicate that PHF20L1 promotes tumorigenesis, supporting the pursuit of PHF20L1 as a target for cancer therapy.

Exploring TCGA database for identification of potential prognostic genes in stomach adenocarcinoma.

BACKGROUND: Stomach adenocarcinoma (STAD) is the fifth most prevalent cancer in the world and ranks third among cancer-related deaths worldwide. The tumour microenvironment (TME) plays an important role in tumorigenesis, development, and metastasis. Hence, we calculated the immune and stromal scores to find the potential prognosis-related genes in STAD using bioinformatics analysis. METHODS: The ESTIMATE algorithm was used to calculate the immune/stromal scores of the STAD samples. Functional enrichment analysis, protein-protein interaction (PPI) network analysis, and overall survival analysis were then performed on differential genes. And we validated these genes using data from the Gene Expression Omnibus database. Finally, we used the Human Protein Atlas (HPA) databases to verify these genes at the protein levels by IHC. RESULTS: Data analysis revealed correlation between stromal/immune scores and the TNM staging system. The top 10 core genes extracted from the PPI network, and primarily involved in immune responses, extracellular matrix, and cell adhesion. There are 31 genes have been validated with poor prognosis and 16 genes were upregulated in tumour tissues compared with normal tissues at the protein level. CONCLUSIONS: In summary, we identified genes associated with the tumour microenvironment with prognostic implications in STAD, which may become potential therapeutic markers leading to better clinical outcomes.

A Ten-N6-Methyladenosine (m6A)-Modified Gene Signature Based on a Risk Score System Predicts Patient Prognosis in Rectum Adenocarcinoma.

OBJECTIVES: The study aims to analyze the expression of N6-methyladenosine (m6A)-modified genes in rectum adenocarcinoma (READ) and identify reliable prognostic biomarkers to predict the prognosis of READ. MATERIALS AND METHODS: RNA sequence data of READ and corresponding clinical survival data were obtained from The Cancer Genome Atlas (TCGA) database. N6-methyladenosine (m6A)-modified genes in READ were downloaded from the "m6Avar" database. Differentially expressed m6A-modified genes in READ stratified by different clinicopathological characteristics were identified using the "limma" package in R. Protein-protein interaction (PPI) network and co-expression analysis of differentially expressed genes (DEGs) were performed using "STRING" and Cytoscape, respectively. Principal component analysis (PCA) was done using R. In addition, Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways were used to functionally annotate the differentially expressed genes in different subgroups. Univariate Cox regression analyses were conducted to identify the powerful independent prognostic factors in READ associated with overall survival (OS). A robust likelihood-based survival model was built using the "rbsurv" package to screen for survival-associated signature genes. The Support Vector Machine (SVM) was used to predict the prognosis of READ through the risk score of survival-associated signature genes. Correlation analysis were carried out using GraphPad prism 8. RESULTS: We screened 974 differentially expressed m6A-modified genes among four types of READ samples. Two READ subgroups (group 1 and group 2) were identified by K means clustering according to the expression of DEGs. The two subgroups were significantly different in overall survival and pathological stages. Next, 118 differentially expressed genes between the two subgroups were screened and the expression of 112 genes was found to be related to the prognosis of READ. Next, a panel of 10 survival-associated signature genes including adamtsl1, csmd2, fam13c, fam184a, klhl4, olfml2b, pdzd4, sec14l5, setbp1, tmem132b was constructed. The signature performed very well for prognosis prediction, time-dependent receiver-operating characteristic (ROC) analysis displaying an area under the curve (AUC) of 0.863, 0.8721, and 0.8752 for 3-year survival rate, prognostic status, and pathological stage prediction, respectively. Correlation analysis showed that the expression levels of the 10 m6A-modified genes were positively correlated with that of m6A demethylase FTO and ALKBH5. CONCLUSION: This study identified potential m6A-modified genes that may be involved in the pathophysiology of READ and constructed a novel gene expression panel for READ risk stratification and prognosis prediction.

GATA3 recruits UTX for gene transcriptional activation to suppress metastasis of breast cancer.

GATA3 has emerged as a prominent transcription factor required for maintaining mammary-gland homeostasis. GATA3 loss is associated with aggressive breast cancer development, but the mechanism by which breast cancer is affected by the loss of GATA3 function remains unclear. Here, we report that GATA3 expression is positively correlated with the expression of UTX, a histone H3K27 demethylase contained in the MLL4 methyltransferase complex, and that GATA3 recruits the chromatin-remodeling MLL4 complex and interacts directly with UTX, ASH2L, and RBBP5. Using RNA sequencing and chromatin immunoprecipitation and sequencing, we demonstrate that the GATA3/UTX complex synergistically regulates a cohort of genes including Dicer and UTX, which are critically involved in the epithelial-to-mesenchymal transition (EMT). Our results further show that the GATA3-UTX-Dicer axis inhibits EMT, invasion, and metastasis of breast cancer cells in vitro and the dissemination of breast cancer in vivo. Our study implicates the GATA3-UTX-Dicer axis in breast cancer metastasis and provides new mechanistic insights into the pathophysiological function of GATA3.

MicroRNA-455-3p mediates GATA3 tumor suppression in mammary epithelial cells by inhibiting TGF-ß signaling.

GATA3 is a basic and essential transcription factor that regulates many pathophysiological processes and is required for the development of mammary luminal epithelial cells. Loss-of-function GATA3 alterations in breast cancer are associated with poor prognosis. Here, we sought to understand the tumor-suppressive functions GATA3 normally performs. We discovered a role for GATA3 in suppressing epithelial-to-mesenchymal transition (EMT) in breast cancer by activating miR-455-3p expression. Enforced expression of miR-455-3p alone partially prevented EMT induced by transforming growth factor ß (TGF-ß) both in cells and tumor xenografts by directly inhibiting key components of TGF-ß signaling. Pathway and biochemical analyses showed that one miRNA-455-3p target, the TGF-ß-induced protein ZEB1, recruits the Mi-2/nucleosome remodeling and deacetylase (NuRD) complex to the promotor region of miR-455 to strictly repress the GATA3-induced transcription of this microRNA. Considering that ZEB1 enhances TGF-ß signaling, we delineated a double-feedback interaction between ZEB1 and miR-455-3p, in addition to the repressive effect of miR-455-3p on TGF-ß signaling. Our study revealed that a feedback loop between these two axes, specifically GATA3-induced miR-455-3p expression, could repress ZEB1 and its recruitment of NuRD (MTA1) to suppress miR-455, which ultimately regulates TGF-ß signaling. In conclusion, we identified that miR-455-3p plays a pivotal role in inhibiting the EMT and TGF-ß signaling pathway and maintaining cell differentiation. This forms the basis of that miR-455-3p might be a promising therapeutic intervention for breast cancer.

Differential expression and functions of Ehm2 transcript variants in lung adenocarcinoma.

Ehm2 [also known as erythrocyte membrane protein band 4.1­like protein 4B (EPB41L4B)] is a member of the NF2/ERM/4.1 superfamily. The overexpression of Ehm2 has been observed in metastatic cancer cells. Through alternative splicing, the Ehm2 gene produces two transcript variants that encode the two different isoforms, Ehm2/1 and Ehm2/2. The biological functions of these different Ehm2 transcript variants remain unclear. The present study aimed to determine the expression of the Ehm2 variants in lung adenocarcinoma and their involvement in the disease progression of the patients. The expression of Ehm2 transcript variants in human lung adenocarcinoma tissues was analyzed using immunohistochemistry and western blot analysis. Ehm2 variants were overexpressed or knocked down in A549 human lung adenocarcinoma cells. The consequent effects of the genetic modifications on the cellular functions of lung cancer cells were then examined using in vitro cell viability, invasion and migration assays. The expression of epithelial­mesenchymal transition (EMT)­related markers was evaluated by western blot analysis in the cell models. The association of Ehm2 variant expression with patient survival was analyzed using Kaplan­Meier survival analysis. The expression of Ehm2/1 was significantly decreased in lung cancers compared with the paired normal lung tissues (P<0.05), while the Ehm2/2 protein levels were higher in the tumors than in the paired normal lung tissues, although this was not statistically significant. The overexpression of Ehm2/1 exerted inhibitory effects, while the knockdown of Ehm2/1 promoted the growth, invasion and migration of A549 cells in vitro. Ehm2/2 was expressed at low levels in the A549 cells and the enforced expression of Ehm2/2 significantly increased the invasiveness and migration of the A549 cells. Immunofluorescence staining revealed that Ehm2/1 was confined to the plasma membrane, while Ehm2/2 was observed at both the plasma membrane and cytoplasm. The overexpression of Ehm2/1 resulted in the upregulation of the epithelial marker, E­cadherin, and in the decreased expression of the mesenchymal markers, N­cadherin and Snail1, while the knockdown of Ehm2/1 and the enforced expression of Ehm2/2 had the opposite effects on the protein levels of EMT­related markers. Kaplan­Meier survival analysis revealed that higher Ehm2/1 transcript levels were associated with the longer survival of patients with lung adenocarcinoma, while the lower expression of Ehm2/2 exhibited a similar association with patient survival. Taken together, the two Ehm2 variants appear to be differentially expressed in lung adenocarcinoma. Ehm2/1 may function as a putative tumor suppressor in the disease progression of lung adenocarcinoma, while Ehm2/2 may have an opposite function.