Analysis of the transcription regulator, CNOT7, as a candidate chromosome 8 tumor suppressor gene in colorectal cancer

Oncogenic Gene CNOT7 Promotes Progression and Induces Poor Prognosis of Glioma

Glioma is the most common malignant brain tumor in the central nervous system with the poor prognosis of patients. The CNOT7 (CCR4-NOT Transcription Complex Subunit 7) is an important functional subunit of CCR4-NOT protein complex that has not been reported in glioma. In this study, we aimed to explore the function of CNOT7 in glioma. The TCGA (The Cancer Genome Atlas) and CGGA (Chinese Glioma Genome Atlas) databases were used for investigating the expression and survival condition of CNOT7 in glioma. The cellular function experiments of qRT-PCR, CCK-8 assays, wound healing assays, and Transwell assays were conducted to verify the function of knockdown CNOT7 in the glioma cell lines DBTRG and U251. Enrichment analysis was used to explore the molecular mechanism of CONT7 in glioma. What is more, the upstream regulation transcription factors of CNOT7 were analyzed based on the ChIP-Atlas and cBioportal (provisional) databases, and verified by the qRT-PCR and luciferase reporter assay. The CNOT7 was highly expressed in glioma and presented the poorer prognosis. The knockdown of CNOT7 inhibited the proliferation, migration, and invasion of glioma cell line, compared to control group. The enrichment analysis revealed that the CNOT7 participated in the development of glioma via G2M checkpoint, E2F targets, IL6-JAK-STAT3, and TNF-α signaling pathways via NF-κB. Besides, it was found that the HDAC2 (Human histone deacetylase-2) contributes to increased CNOT7 expression in glioma. The high-expressed CNOT7 is an oncogene with poor prognosis and participate the progression of glioma.

Transcriptomic correlates of cell cycle checkpoints with distinct prognosis, molecular characteristics, immunological regulation, and therapeutic response in colorectal adenocarcinoma

Backgrounds

Colorectal adenocarcinoma (COAD), accounting for the most common subtype of colorectal cancer (CRC), is a kind of malignant digestive tumor. Some cell cycle checkpoints (CCCs) have been found to contribute to CRC progression, whereas the functional roles of a lot of CCCs, especially the integrated role of checkpoint mechanism in the cell cycle, remain unclear.

Materials and methods

The Genomic Data Commons (GDC) The Cancer Genome Atlas (TCGA) COAD cohort was retrieved as the training dataset, and GSE24551 and GSE29623 were downloaded from Gene Expression Omnibus (GEO) as the validation datasets. A total of 209 CCC-related genes were derived from the Gene Ontology Consortium and were subsequently enrolled in the univariate, multivariate, and least absolute shrinkage and selection operator (LASSO) Cox regression analyses, finally defining a CCC signature. Cell proliferation and Transwell assay analyses were utilized to evaluate the functional roles of signature-related CCCs. The underlying CCC signature, molecular characteristics, immune-related features, and therapeutic response were finally estimated. The Genomics of Drug Sensitivity in Cancer (GDSC) database was employed for the evaluation of chemotherapeutic responses.

Results

The aberrant gene expression of CCCs greatly contributed to COAD development and progression. Univariate Cox regression analysis identified 27 CCC-related genes significantly affecting the overall survival (OS) of COAD patients; subsequently, LASSO analysis determined a novel CCC signature. Noticeably, CDK5RAP2, MAD1L1, NBN, RGCC, and ZNF207 were first identified to be correlated with the prognosis of COAD, and it was proven that all of them were significantly correlated with the proliferation and invasion of HCT116 and SW480 cells. In TCGA COAD cohort, CCC signature robustly stratified COAD patients into high and low CCC score groups (median OS: 57.24 months vs. unreached, p< 0.0001), simultaneously, with the good AUC values for OS prediction at 1, 2, and 3 years were 0.74, 0.78, and 0.77. Furthermore, the prognostic capacity of the CCC signature was verified in the GSE24551 and GSE29623 datasets, and the CCC signature was independent of clinical features. Moreover, a higher CCC score always indicated worse OS, regardless of clinical features, histological subtypes, or molecular subgroups. Intriguingly, functional enrichment analysis confirmed the CCC score was markedly associated with extracellular, matrix and immune (chemokine)-related signaling, cell cycle-related signaling, and metabolisms. Impressively, a higher CCC score was positively correlated with a majority of chemokines, receptors, immunostimulators, and anticancer immunity, indicating a relatively immune-promoting microenvironment. In addition, GSE173839, GSE25066, GSE41998, and GSE194040 dataset analyses of the underlying CCC signature suggested that durvalumab with olaparib and paclitaxel, taxane-anthracycline chemotherapy, neoadjuvant cyclophosphamide/doxorubicin with ixabepilone or paclitaxel, and immunotherapeutic strategies might be suitable for COAD patients with higher CCC score. Eventually, the GDSC database analysis showed that lower CCC scores were likely to be more sensitive to 5-fluorouracil, bosutinib, gemcitabine, gefitinib, methotrexate, mitomycin C, and temozolomide, while patients with higher CCC score seemed to have a higher level of sensitivity to bortezomib and elesclomol.

Conclusion

The novel CCC signature exhibited a good ability for prognosis prediction for COAD patients, and the CCC score was found to be highly correlated with molecular features, immune-related characteristics, and therapeutic responses, which would greatly promote clinical management and precision medicine for COAD.

Methylation recognition protein YTH N6-methyladenosine RNA binding protein 1 (YTHDF1) regulates the proliferation, migration and invasion of osteosarcoma by regulating m6A level of CCR4-NOT transcription complex subunit 7 (CNOT7)

N6-methyladenosine (m6A) is one of the most significant modifications in human mRNAs. Emerging evidence indicates that m6A participates in the initiation and development of malignant tumors. Nevertheless, the biological roles and mechanism of m6A in osteosarcoma (OS) remain unclear. The purpose of this study was to investigate the role and mechanism of the methylation recognition protein-YTH N6-methyladenosine RNA binding protein 1 (YTHDF1) in OS. The YTHDF1 expression in OS was detected by qRT-PCR and Western blot assay. M6A quantification was utilized to measure the methylation level of OS. Cell counting kit-8 (CCK8), 5-Ethynyl-2’-deoxyuridine (EdU) assay and transwell experiments were conducted to confirm the biological effects of YTHDF1 on OS cells. The bioinformatics websites and in vitro assays were conducted to analyze the downstream targets of YTHDF1 was upregulated in OS tissues at mRNA and protein level. The results showed that the expression level of YTHDF1 might be closely associated with the poor prognosis for OS patients. Inhibition of YTHDF1 could suppress the proliferation, migration and invasion of the OS cells. Moreover, we found that CCR4-NOT transcription complex subunit 7 (CNOT7) might be the potential target of YTHDF1, which was upregulated in OS tissues. YTHDF1 could recognize the m6A sites of CONT7 and promote its expression in an m6A manner. Moreover, methyltransferase-like 3 (METTL3) could promote the m6A level of CONT7. YTHDF1 was upregulated in OS and could promote cell proliferation, migration and invasion. The METTL3-CONT7-YTHDF1 regulatory axis might be the potential target for the prognosis and therapy of OS.

Human Ccr4 and Caf1 Deadenylases Regulate Proliferation and Tumorigenicity of Human Gastric Cancer Cells via Modulating Cell Cycle Progression

Simple Summary

Cancer cells generally reprogram their gene expression profiles to satisfy continuous growth, proliferation, and metastasis. Most eukaryotic mRNAs are degraded in a deadenylation-dependent pathway, in which deadenylases are the key enzymes. We found that human Ccr4 (hCcr4a/b) and Caf1 (hCaf1a/b), the dominant cytosolic deadenylases, were dysregulated in several types of cancers including stomach adenocarcinoma. Stably knocking down hCaf1a/b or hCcr4a/b blocks cell cycle progression by enhancing the levels of cell cycle inhibitors and by inhibiting the formation of processing bodies, which are cytosolic foci involved in mRNA metabolism. More importantly, depletion of hCaf1a/b or hCcr4a/b dramatically inhibits cell proliferation and tumorigenicity. Our results suggest that perturbating global RNA metabolism may provide a potential novel strategy for cancer treatment.

Abstract

Cancer cells generally have reprogrammed gene expression profiles to meet the requirements of survival, continuous division, and metastasis. An interesting question is whether the cancer cells will be affected by interfering their global RNA metabolism. In this research, we found that human Ccr4a/b (hCcr4a/b) and Caf1a/b (hCaf1a/b) deadenylases, the catalytic components of the Ccr4-Not complex, were dysregulated in several types of cancers including stomach adenocarcinoma. The impacts of the four deadenylases on cancer cell growth were studied by the establishment of four stable MKN28 cell lines with the knockdown of hCcr4a/b or hCaf1a/b or transient knockdown in several cell lines. Depletion of hCcr4a/b or hCaf1a/b significantly inhibited cell proliferation and tumorigenicity. Mechanistic studies indicated that the cells were arrested at the G2/M phase by knocking down hCaf1a, while arrested at the G0/G1 phase by depleting hCaf1b or hCcr4a/b. The four enzymes did not affect the levels of CDKs and cyclins but modulated the levels of CDK–cyclin inhibitors. We identified that hCcr4a/b, but not hCaf1a/b, targeted the p21 mRNA in the MKN28 cells. Furthermore, depletion of any one of the four deadenylases dramatically impaired processing-body formation in the MKN28 and HEK-293T cells. Our results highlight that perturbating global RNA metabolism may severely affect cancer cell proliferation, which provides a potential novel strategy for cancer treatment.

Diverse functions of deadenylases in DNA damage response and genomic integrity

DNA damage response (DDR) is a coordinated network of diverse cellular processes including the detection, signaling, and repair of DNA lesions, the adjustment of metabolic network and cell fate determination. To deal with the unavoidable DNA damage caused by either endogenous or exogenous stresses, the cells need to reshape the gene expression profile to allow efficient transcription and translation of DDR-responsive messenger RNAs (mRNAs) and to repress the nonessential mRNAs. A predominant method to adjust RNA fate is achieved by modulating the 3'-end oligo(A) or poly(A) length via the opposing actions of polyadenylation and deadenylation. Poly(A)-specific ribonuclease (PARN) and the carbon catabolite repressor 4 (CCR4)-Not complex, the major executors of deadenylation, are indispensable to DDR and genomic integrity in eukaryotic cells. PARN modulates cell cycle progression by regulating the stabilities of mRNAs and microRNA (miRNAs) involved in the p53 pathway and contributes to genomic stability by affecting the biogenesis of noncoding RNAs including miRNAs and telomeric RNA. The CCR4-Not complex is involved in diverse pathways of DDR including transcriptional regulation, signaling pathways, mRNA stabilities, translation regulation, and protein degradation. The RNA targets of deadenylases are tuned by the DDR signaling pathways, while in turn the deadenylases can regulate the levels of DNA damage-responsive proteins. The mutual feedback between deadenylases and the DDR signaling pathways allows the cells to precisely control DDR by dynamically adjusting the levels of sensors and effectors of the DDR signaling pathways. Here, the diverse functions of deadenylases in DDR are summarized and the underlying mechanisms are proposed according to recent findings. This article is categorized under: RNA Processing > 3' End Processing RNA in Disease and Development > RNA in Disease RNA Turnover and Surveillance > Turnover/Surveillance Mechanisms.

Identification of Crucial Candidate Genes and Pathways in Glioblastoma Multiform by Bioinformatics Analysis

The present study aimed to investigate the molecular mechanisms underlying glioblastoma multiform (GBM) and its biomarkers. The differentially expressed genes (DEGs) were diagnosed using the limma software package. The ToppGene (ToppFun) was used to perform pathway and Gene Ontology (GO) enrichment analysis of the DEGs. Protein-protein interaction (PPI) networks, extracted modules, miRNA-target genes regulatory network and TF-target genes regulatory network were used to obtain insight into the actions of DEGs. Survival analysis for DEGs was carried out. A total of 590 DEGs, including 243 up regulated and 347 down regulated genes, were diagnosed between scrambled shRNA expression and Lin7A knock down. The up-regulated genes were enriched in ribosome, mitochondrial translation termination, translation, and peptide biosynthetic process. The down-regulated genes were enriched in focal adhesion, VEGFR3 signaling in lymphatic endothelium, extracellular matrix organization, and extracellular matrix. The current study screened the genes in the PPI network, extracted modules, miRNA-target genes regulatory network, and TF-target genes regulatory network with higher degrees as hub genes, which included NPM1, CUL4A, YIPF1, SHC1, AKT1, VLDLR, RPL14, P3H2, DTNA, FAM126B, RPL34, and MYL5. Survival analysis indicated that the high expression of RPL36A and MRPL35 were predicting longer survival of GBM, while high expression of AP1S1 and AKAP12 were predicting shorter survival of GBM. High expression of RPL36A and AP1S1 were associated with pathogenesis of GBM, while low expression of ALPL was associated with pathogenesis of GBM. In conclusion, the current study diagnosed DEGs between scrambled shRNA expression and Lin7A knock down samples, which could improve our understanding of the molecular mechanisms in the progression of GBM, and these crucial as well as new diagnostic markers might be used as therapeutic targets for GBM.

Deadenylation: enzymes, regulation, and functional implications

Lengths of the eukaryotic messenger RNA (mRNA) poly(A) tails are dynamically changed by the opposing effects of poly(A) polymerases and deadenylases. Modulating poly(A) tail length provides a highly regulated means to control almost every stage of mRNA lifecycle including transcription, processing, quality control, transport, translation, silence, and decay. The existence of diverse deadenylases with distinct properties highlights the importance of regulating poly(A) tail length in cellular functions. The deadenylation activity can be modulated by subcellular locations of the deadenylases, cis-acting elements in the target mRNAs, trans-acting RNA-binding proteins, posttranslational modifications of deadenylase and associated factors, as well as transcriptional and posttranscriptional regulation of the deadenylase genes. Among these regulators, the physiological functions of deadenylases are largely dependent on the interactions with the trans-acting RNA-binding proteins, which recruit deadenylases to the target mRNAs. The task of these RNA-binding proteins is to find and mark the target mRNAs based on their sequence features. Regulation of the regulators can switch on or switch off deadenylation and thereby destabilize or stabilize the targeted mRNAs, respectively. The distinct domain compositions and cofactors provide various deadenylases the structural basis for the recruitments by distinct RNA-binding protein subsets to meet dissimilar cellular demands. The diverse deadenylases, the numerous types of regulators, and the reversible posttranslational modifications together make up a complicated network to precisely regulate intracellular mRNA homeostasis. This review will focus on the diverse regulators of various deadenylases and will discuss their functional implications, remaining problems, and future challenges.

A metastasis modifier locus on human chromosome 8p in uveal melanoma identified by integrative genomic analysis

PURPOSE

To identify genes that modify metastatic risk in uveal melanoma, a type of cancer that is valuable for studying metastasis because of its remarkably consistent metastatic pattern and well-characterized gene expression signature associated with metastasis.

EXPERIMENTAL DESIGN

We analyzed 53 primary uveal melanomas by gene expression profiling, array-based comparative genomic hybridization, array-based global DNA methylation profiling, and single nucleotide polymorphism-based detection of loss of heterozygosity to identify modifiers of metastatic risk. A candidate gene, leucine zipper tumor suppressor-1 (LZTS1), was examined for its effect on proliferation, migration, and motility in cultured uveal melanoma cells.

RESULTS

In metastasizing primary uveal melanomas, deletion of chromosome 8p12-22 and DNA hypermethylation of the corresponding region of the retained hemizygous 8p allele were associated with more rapid metastasis. Among the 11 genes located within the deleted region, LZTS1 was most strongly linked to rapid metastasis. LZTS1 was silenced in rapidly metastasizing and metastatic uveal melanomas but not in slowly metastasizing and nonmetastasizing uveal melanomas. Forced expression of LZTS1 in metastasizing uveal melanoma cells inhibited their motility and invasion, whereas depletion of LZTS1 increased their motility.

CONCLUSIONS

We have described a metastatic modifier locus on chromosome 8p and identified LZTS1 as a potential metastasis suppressor within this region. This study shows the utility of integrative genomic methods for identifying modifiers of metastatic risk in human cancers and may suggest new therapeutic targets in metastasizing tumor cells.

Human MicroRNA targets

MicroRNAs (miRNAs) interact with target mRNAs at specific sites to induce cleavage of the message or inhibit translation. The specific function of most mammalian miRNAs is unknown. We have predicted target sites on the 3' untranslated regions of human gene transcripts for all currently known 218 mammalian miRNAs to facilitate focused experiments. We report about 2,000 human genes with miRNA target sites conserved in mammals and about 250 human genes conserved as targets between mammals and fish. The prediction algorithm optimizes sequence complementarity using position-specific rules and relies on strict requirements of interspecies conservation. Experimental support for the validity of the method comes from known targets and from strong enrichment of predicted targets in mRNAs associated with the fragile X mental retardation protein in mammals. This is consistent with the hypothesis that miRNAs act as sequence-specific adaptors in the interaction of ribonuclear particles with translationally regulated messages. Overrepresented groups of targets include mRNAs coding for transcription factors, components of the miRNA machinery, and other proteins involved in translational regulation, as well as components of the ubiquitin machinery, representing novel feedback loops in gene regulation. Detailed information about target genes, target processes, and open-source software for target prediction (miRanda) is available at http://www.microrna.org. Our analysis suggests that miRNA genes, which are about 1% of all human genes, regulate protein production for 10% or more of all human genes.

Mapping of a candidate colorectal cancer tumor-suppressor gene to a 900-kilobase region on the short arm of chromosome 8

Loss of heterozygosity (LOH) on 8p occurs at high frequencies in many tumor types, including colorectal carcinoma (CRC). We previously used microcell-mediated chromosome transfer (MMCT) into the CRC cell line SW620 to map a approximately 7.7-Mb colorectal cancer-suppressor region (CRCSR) at 8p22-23.1. In the current study, we transferred small fragments of this CRCSR into SW620 to refine the region further. Two microcell hybrids containing a 321- to 898-kb region around the D8S552 marker at 8p23.1 showed suppression of soft agar clonicity and tumorigenicity in athymic mice when compared to control cell lines. These data suggest that the putative colorectal tumor-suppressor gene is within this small region. We analyzed two candidate genes within this region: FLJ23749 and KIAA1456. Expression of both genes was detected in normal colonic crypt cells and in mucosa. Quantitative reverse transcriptase polymerase chain reaction showed downregulation of KIAA1456 in 9 of 12 primary colorectal tumors compared to matching normal mucosa, but normal or increased expression of FLJ23749. FLJ23749 was expressed in all CRC cell lines tested; however, KIAA1456 was downregulated in three cell lines, including SW620, and was restored in the suppressed microcell hybrids. 5'aza-2'Deoxycytidine treatment of the downregulated cell lines restored expression of KIAA1456, but bisulfite genomic sequencing did not show a correlation between promoter methylation and expression. Forty percent of the primary tumors showed LOH at this CRCSR locus, and mutation analysis revealed somatic mutations in 1 of 88 primary colorectal tumors for both KIAA1456 and FLJ23749. Despite the rarity of somatic mutations, the expression data suggest that KIAA1456 is still a candidate for the putative 8p colorectal cancer tumor-suppressor gene.