Genome-wide mapping of DNA-binding sites identifies stemness-related genes as directly repressed targets of SNAIL1 in colorectal cancer cells

The epigenetic landscape in intestinal stem cells and its deregulation in colorectal cancer

Epigenetic mechanisms play a pivotal role in controlling gene expression and cellular plasticity in both normal physiology and pathophysiological conditions. These mechanisms are particularly important in the regulation of stem cell self-renewal and differentiation, both in embryonic development and within adult tissues. A prime example of this finely tuned epigenetic control is observed in the gastrointestinal lining, where the small intestine undergoes renewal approximately every 3-5 days. How various epigenetic mechanisms modulate chromatin functions in intestinal stem cells (ISCs) is currently an active area of research. In this review, we discuss the main epigenetic mechanisms that control ISC differentiation under normal homeostasis. Furthermore, we explore the dysregulation of these mechanisms in the context of colorectal cancer (CRC) development. By outlining the main epigenetic mechanisms contributing to CRC, we highlight the recent therapeutics development and future directions for colorectal cancer research.

The Roles of Zinc Finger Proteins in Colorectal Cancer

Despite colorectal cancer remaining a leading worldwide cause of cancer-related death, there remains a paucity of effective treatments for advanced disease. The molecular mechanisms underlying the development of colorectal cancer include altered cell signaling and cell cycle regulation that may result from epigenetic modifications of gene expression and function. Acting as important transcriptional regulators of normal biological processes, zinc finger proteins also play key roles in regulating the cellular mechanisms underlying colorectal neoplasia. These actions impact cell differentiation and proliferation, epithelial-mesenchymal transition, apoptosis, homeostasis, senescence, and maintenance of stemness. With the goal of highlighting promising points of therapeutic intervention, we review the oncogenic and tumor suppressor roles of zinc finger proteins with respect to colorectal cancer tumorigenesis and progression.

Mammalian evolution of human cis-regulatory elements and transcription factor binding sites

Understanding the regulatory landscape of the human genome is a long-standing objective of modern biology. Using the reference-free alignment across 241 mammalian genomes produced by the Zoonomia Consortium, we charted evolutionary trajectories for 0.92 million human candidate cis-regulatory elements (cCREs) and 15.6 million human transcription factor binding sites (TFBSs). We identified 439,461 cCREs and 2,024,062 TFBSs under evolutionary constraint. Genes near constrained elements perform fundamental cellular processes, whereas genes near primate-specific elements are involved in environmental interaction, including odor perception and immune response. About 20% of TFBSs are transposable element-derived and exhibit intricate patterns of gains and losses during primate evolution whereas sequence variants associated with complex traits are enriched in constrained TFBSs. Our annotations illuminate the regulatory functions of the human genome.

Noncanonical Wnt signaling promotes colon tumor growth, chemoresistance and tumor fibroblast activation

Colon tumors of the mesenchymal subtype have the lowest overall survival. Snail1 is essential for the acquisition of this phenotype, characterized by increased tumor stemness and invasion, and high resistance to chemotherapy. Here, we find that Snail1 expression in colon tumor cells is dependent on an autocrine noncanonical Wnt pathway. Accordingly, depletion of Ror2, the co-receptor for noncanonical Wnts such as Wnt5a, potently decreases Snail1 expression. Wnt5a, Ror2, and Snail1 participate in a self-stimulatory feedback loop since Wnt5a increases its own synthesis in a Ror2- and Snail1-dependent fashion. This Wnt5a/Ror2/Snail1 axis controls tumor invasion, chemoresistance, and formation of tumor spheres. It also stimulates TGFβ synthesis; consequently, tumor cells expressing Snail1 are more efficient in activating cancer-associated fibroblasts than the corresponding controls. Ror2 downmodulation or inhibition of the Wnt5a pathway decreases Snail1 expression in primary colon tumor cells and their ability to form tumors and liver metastases. Finally, the expression of SNAI1, ROR2, and WNT5A correlates in human colon and other tumors. These results identify inhibition of the noncanonical Wnt pathway as a putative colon tumor therapy.

TGFβ1-Induced EMT in the MCF10A Mammary Epithelial Cell Line Model Is Executed Independently of SNAIL1 and ZEB1 but Relies on JUNB-Coordinated Transcriptional Regulation

Epithelial-mesenchymal transition (EMT) fosters cancer cell invasion and metastasis, the main cause of cancer-related mortality. Growing evidence that SNAIL and ZEB transcription factors, typically portrayed as master regulators of EMT, may be dispensable for this process, led us to re-investigate its mechanistic underpinnings. For this, we used an unbiased computational approach that integrated time-resolved analyses of chromatin structure and differential gene expression, to predict transcriptional regulators of TGFβ1-inducible EMT in the MCF10A mammary epithelial cell line model. Bioinformatic analyses indicated comparatively minor contributions of SNAIL proteins and ZEB1 to TGFβ1-induced EMT, whereas the AP-1 subunit JUNB was anticipated to have a much larger impact. CRISPR/Cas9-mediated loss-of-function studies confirmed that TGFβ1-induced EMT proceeded independently of SNAIL proteins and ZEB1. In contrast, JUNB was necessary and sufficient for EMT in MCF10A cells, but not in A549 lung cancer cells, indicating cell-type-specificity of JUNB EMT-regulatory capacity. Nonetheless, the JUNB-dependence of EMT-associated transcriptional reprogramming in MCF10A cells allowed to define a gene expression signature which was regulated by TGFβ1 in diverse cellular backgrounds, showed positively correlated expression with TGFβ signaling in multiple cancer transcriptomes, and was predictive of patient survival in several cancer types. Altogether, our findings provide novel mechanistic insights into the context-dependent control of TGFβ1-driven EMT and thereby may lead to improved diagnostic and therapeutic options.

Inside the stemness engine: Mechanistic links between deregulated transcription factors and stemness in cancer

Cell identity is largely determined by its transcriptional profile. In tumour, deregulation of transcription factor expression and/or activity enables cancer cell to acquire a stem-like state characterised by capacity to self-renew, differentiate and form tumours in vivo. These stem-like cancer cells are highly metastatic and therapy resistant, thus warranting a more complete understanding of the molecular mechanisms downstream of the transcription factors that mediate the establishment of stemness state. Here, we review recent research findings that provide a mechanistic link between the commonly deregulated transcription factors and stemness in cancer. In particular, we describe the role of master transcription factors (SOX, OCT4, NANOG, KLF, BRACHYURY, SALL, HOX, FOX and RUNX), signalling-regulated transcription factors (SMAD, β-catenin, YAP, TAZ, AP-1, NOTCH, STAT, GLI, ETS and NF-κB) and unclassified transcription factors (c-MYC, HIF, EMT transcription factors and P53) across diverse tumour types, thereby yielding a comprehensive overview identifying shared downstream targets, highlighting unique mechanisms and discussing complexities.

Loss of SNAI1 induces cellular plasticity in invasive triple-negative breast cancer cells

The transcription factor SNAI1 mediates epithelial-mesenchymal transition, fibroblast activation and controls inter-tissue migration. High SNAI1 expression characterizes metastatic triple-negative breast carcinomas, and its knockout by CRISPR/Cas9 uncovered an epithelio-mesenchymal phenotype accompanied by reduced signaling by the cytokine TGFβ. The SNAI1 knockout cells exhibited plasticity in differentiation, drifting towards the luminal phenotype, gained stemness potential and could differentiate into acinar mammospheres in 3D culture. Loss of SNAI1 de-repressed the transcription factor FOXA1, a pioneering factor of mammary luminal progenitors. FOXA1 induced a specific gene program, including the androgen receptor (AR). Inhibiting AR via a specific antagonist regenerated the basal phenotype and blocked acinar differentiation. Thus, loss of SNAI1 in the context of triple-negative breast carcinoma cells promotes an intermediary luminal progenitor phenotype that gains differentiation plasticity based on the dual transcriptional action of FOXA1 and AR. This function of SNAI1 provides means to separate cell invasiveness from progenitor cell de-differentiation as independent cellular programs.

Snail acetylation by autophagy-derived acetyl-coenzyme A promotes invasion and metastasis of KRAS-LKB1 co-mutated lung cancer cells

Background

Autophagy is elevated in metastatic tumors and is often associated with active epithelial‐to‐mesenchymal transition (EMT). However, the extent to which EMT is dependent on autophagy is largely unknown. This study aimed to identify the mechanisms by which autophagy facilitates EMT.

Methods

We employed a liquid chromatography‐based metabolomic approach with kirsten rat sarcoma viral oncogene (KRAS) and liver kinase B1 (LKB1) gene co‐mutated (KL) cells that represent an autophagy/EMT‐coactivated invasive lung cancer subtype for the identification of metabolites linked to autophagy‐driven EMT activation. Molecular mechanisms of autophagy‐driven EMT activation were further investigated by quantitative real‐time polymerase chain reaction (qRT‐PCR), Western blotting analysis, immunoprecipitation, immunofluorescence staining, and metabolite assays. The effects of chemical and genetic perturbations on autophagic flux were assessed by two orthogonal approaches: microtubule‐associated protein 1A/1B‐light chain 3 (LC3) turnover analysis by Western blotting and monomeric red fluorescent protein‐green fluorescent protein (mRFP‐GFP)‐LC3 tandem fluorescent protein quenching assay. Transcription factor EB (TFEB) activity was measured by coordinated lysosomal expression and regulation (CLEAR) motif‐driven luciferase reporter assay. Experimental metastasis (tail vein injection) mouse models were used to evaluate the impact of calcium/calmodulin‐dependent protein kinase kinase 2 (CAMKK2) or ATP citrate lyase (ACLY) inhibitors on lung metastasis using IVIS luciferase imaging system.

Results

We found that autophagy in KL cancer cells increased acetyl‐coenzyme A (acetyl‐CoA), which facilitated the acetylation and stabilization of the EMT‐inducing transcription factor Snail. The autophagy/acetyl‐CoA/acetyl‐Snail axis was further validated in tumor tissues and in autophagy‐activated pancreatic cancer cells. TFEB acetylation in KL cancer cells sustained pro‐metastatic autophagy in a mammalian target of rapamycin complex 1 (mTORC1)‐independent manner. Pharmacological inhibition of this axis via CAMKK2 inhibitors or ACLY inhibitors consistently reduced the metastatic capacity of KL cancer cells in vivo.

Conclusions

This study demonstrates that autophagy‐derived acetyl‐CoA promotes Snail acetylation and thereby facilitates invasion and metastasis of KRAS‐LKB1 co‐mutated lung cancer cells and that inhibition of the autophagy/acetyl‐CoA/acetyl‐Snail axis using CAMKK2 or ACLY inhibitors could be a potential therapeutic strategy to suppress metastasis of KL lung cancer.

LncRNAs, the Molecules Involved in Communications With Colorectal Cancer Stem Cells

Colorectal cancer stem cells (CRCSCs) can actively self-renew, as well as having multidirectional differentiation and tumor regeneration abilities. Because the high functional activities of CRCSCs are associated with low cure rates in patients with colorectal cancer, efforts have sought to determine the function and regulatory mechanisms of CRCSCs. To date, however, the potential regulatory mechanisms of CRCSCs remain incompletely understood. Many non-coding genes are involved in tumor invasion and spread through their regulation of CRCSCs, with long non-coding RNAs (lncRNAs) being important non-coding RNAs. LncRNAs may be involved in the colorectal cancer development and drug resistance through their regulation of CRCSCs. This review systematically evaluates the latest research on the ability of lncRNAs to regulate CRCSC signaling pathways and the involvement of these lncRNAs in colorectal cancer promotion and suppression. The regulatory network of lncRNAs in the CRCSC signaling pathway has been determined. Further analysis of the potential clinical applications of lncRNAs as novel clinical diagnostic and prognostic biomarkers and therapeutic targets for colorectal cancer may provide new ideas and protocols for the prevention and treatment of colorectal cancer.

Canonical TGFβ signaling induces collective invasion in colorectal carcinogenesis through a Snail1- and Zeb1-independent partial EMT

Local invasion is the initial step towards metastasis, the main cause of cancer mortality. In human colorectal cancer (CRC), malignant cells predominantly invade as cohesive collectives and may undergo partial epithelial-mesenchymal transition (pEMT) at the invasive front. How this particular mode of stromal infiltration is generated is unknown. Here we investigated the impact of oncogenic transformation and the microenvironment on tumor cell invasion using genetically engineered organoids as CRC models. We found that inactivation of the Apc tumor suppressor combined with expression of oncogenic Kras^G12D and dominant-negative Trp53^R172H did not cell-autonomously induce invasion in vitro. However, oncogenic transformation primed organoids for activation of a collective invasion program upon exposure to the prototypical microenvironmental factor TGFβ1. Execution of this program co-depended on a permissive extracellular matrix which was further actively remodeled by invading organoids. Although organoids shed some epithelial properties particularly at the invasive edge, TGFβ1-stimulated organoids largely maintained epithelial gene expression while additionally implementing a mesenchymal transcription pattern, resulting in a pEMT phenotype that did not progress to a fully mesenchymal state. Notably, while TGFβ1 induced pEMT and promoted collective invasion, it abrogated self-renewal capacity of TKA organoids which correlated with the downregulation of intestinal stem cell (ISC) marker genes. Mechanistically, induction of the non-progressive pEMT required canonical TGFβ signaling mediated by Smad transcription factors (TFs), whereas the EMT master regulators Snail1 and Zeb1 were dispensable. Gene expression profiling provided further evidence for pEMT of TGFβ1-treated organoids and showed that their transcriptomes resemble those of human poor prognosis CMS4 cancers which likewise exhibit pEMT features. We propose that collective invasion in colorectal carcinogenesis is triggered by microenvironmental stimuli through activation of a novel, transcription-mediated form of non-progressive pEMT independently of classical EMT regulators.

Transcription factor c-Myb: novel prognostic factor in osteosarcoma

The transcription factor c-Myb is an oncoprotein promoting cell proliferation and survival when aberrantly activated/expressed, thus contributing to malignant transformation. Overexpression of c-Myb has been found in leukemias, breast, colon and adenoid cystic carcinoma. Recent studies revealed its expression also in osteosarcoma cell lines and suggested its functional importance during bone development. However, the relevance of c-Myb in control of osteosarcoma progression remains unknown. A retrospective clinical study was carried out to assess a relationship between c-Myb expression in archival osteosarcoma tissues and prognosis in a cohort of high-grade osteosarcoma patients. In addition, MYB was depleted in metastatic osteosarcoma cell lines SAOS-2 LM5 and 143B and their growth, chemosensitivity, migration and metastatic activity were determined. Immunohistochemical analysis revealed that high c-Myb expression was significantly associated with poor overall survival in the cohort and metastatic progression in young patients. Increased level of c-Myb was detected in metastatic osteosarcoma cell lines and its depletion suppressed their growth, colony-forming capacity, migration and chemoresistance in vitro in a cell line-dependent manner. MYB knock-out resulted in reduced metastatic activity of both SAOS-2 LM5 and 143B cell lines in immunodeficient mice. Transcriptomic analysis revealed the c-Myb-driven functional programs enriched for genes involved in the regulation of cell growth, stress response, cell adhesion and cell differentiation/morphogenesis. Wnt signaling pathway was identified as c-Myb target in osteosarcoma cells. Taken together, we identified c-Myb as a negative prognostic factor in osteosarcoma and showed its involvement in the regulation of osteosarcoma cell growth, chemosensitivity, migration and metastatic activity.

SMAD4 mutations do not preclude epithelial-mesenchymal transition in colorectal cancer

Transforming growth factor beta (TGFβ) superfamily signaling is a prime inducer of epithelial-mesenchymal transitions (EMT) that foster cancer cell invasion and metastasis, a major cause of cancer-related deaths. Yet, TGFβ signaling is frequently inactivated in human tumor entities including colorectal cancer (CRC) and pancreatic adenocarcinoma (PAAD) with a high proportion of mutations incapacitating SMAD4, which codes for a transcription factor (TF) central to canonical TGFβ and bone morphogenetic protein (BMP) signaling. Beyond its role in initiating EMT, SMAD4 was reported to crucially contribute to subsequent gene regulatory events during EMT execution. It is therefore widely assumed that SMAD4-mutant (SMAD4^mut) cancer cells are unable to undergo EMT. Here, we scrutinized this notion and probed for potential SMAD4-independent EMT execution using SMAD4^mut CRC cell lines. We show that SMAD4^mut cells exhibit morphological changes, become invasive, and regulate EMT marker genes upon induction of the EMT-TF SNAIL1. Furthermore, SNAIL1-induced EMT in SMAD4^mut cells was found to be entirely independent of TGFβ/BMP receptor activity. Global assessment of the SNAIL1-dependent transcriptome confirmed the manifestation of an EMT gene regulatory program in SMAD4^mut cells highly related to established EMT signatures. Finally, analyses of human tumor transcriptomes showed that SMAD4 mutations are not underrepresented in mesenchymal tumor samples and that expression patterns of EMT-associated genes are similar in SMAD4^mut and SMAD4 wild-type (SMAD4^wt) cases. Altogether, our findings suggest that alternative TFs take over the gene regulatory functions of SMAD4 downstream of EMT-TFs, arguing for considerable plasticity of gene regulatory networks operating in EMT execution. Further, they establish that EMT is not categorically precluded in SMAD4^mut tumors, which is relevant for their diagnostic and therapeutic evaluation.

Comprehensive analysis of EMT-related genes and lncRNAs in the prognosis, immunity, and drug treatment of colorectal cancer

BACKGROUND

EMT is an important biological process in the mechanism of tumor invasion and metastasis. However, there are still many unknowns about the specific mechanism of EMT in tumor. At present, a comprehensive analysis of EMT-related genes in colorectal cancer (CRC) is still lacking.

METHODS

All the data were downloaded from public databases including TCGA database (488 tumor samples and 52 normal samples) as the training set and the GEO database (GSE40967 including 566 tumor samples and 19 normal samples, GSE12945 including 62 tumor samples, GSE17536 including 177 tumor samples, GSE17537 including 55 tumor samples) as the validation sets. One hundred and sixty-six EMT-related genes (EMT-RDGs) were selected from the Molecular Signatures Database. Bioinformatics methods were used to analyze the correlation between EMT-RDGs and CRC prognosis, metastasis, drug efficacy, and immunity.

RESULTS

We finally obtained nine prognostic-related EMT-RDGs (FGF8, NOG, PHLDB2, SIX2, SNAI1, TBX5, TIAM1, TWIST1, TCF15) through differential expression analysis, Unicox and Lasso regression analysis, and then constructed a risk prognosis model. There were significant differences in clinical characteristics, 22 immune cells, and immune functions between the high-risk and low-risk groups and the different states of the nine prognostic-related EMT-RDGs. The methylation level and mutation status of nine prognostic-related EMT-RDGs all affect their regulation of EMT. The Cox proportional hazards regression model was also constructed by the methylation sites of nine prognostic-related EMT-RDGs. In addition, the expression of FGF8, PHLDB2, SIX2, and SNAIL was higher and the expression level of NOG and TWIST1 was lower in the non-metastasis CRC group. Nine prognostic-related EMT-RDGs also affected the drug treatment response of CRC.

CONCLUSIONS

Targeting these nine prognostic-related EMT-RDGs can regulate CRC metastasis and immune, which is beneficial for the prognosis of CRC patients, improve drug sensitivity in CRC patients.

Coordination of PRKCA/PRKCA-AS1 interplay facilitates DNA methyltransferase 1 recruitment on DNA methylation to affect protein kinase C alpha transcription in mitral valve of rheumatic heart disease

In the present study, mitral valve tissues from three mitral stenosis patients with RHD by valve replacement and two healthy donors were harvested and conducted DNA methylation signature on PRKCA by MeDIP-qPCR. The presence of hypomethylated CpG islands at promoter and 5' terminal of PRKCA was observed in RHD accompanied with highly expressed PRKCA and down-regulated antisense long non-coding RNA (lncRNA) PRKCA-AS1 compared to health control. Furthermore, the enrichments of DNMT1/3A/3B on PRKCA were detected by ChIP-qPCR assay in vivo and in human cardiomyocyte AC16 and RL-14 cells exposed to TNF-α in vitro, and both demonstrated that DNMT1 substantially contributed to DNA methylation. Additionally, PRKCA-AS1 was further determined to bind with promoter of PRKCA via 5' terminal and interact with DNMT1 via 3' terminal. Taken together, our results illuminated a novel regulatory mechanism of DNA methylation on regulating PRKCA transcription through lncRNA PRKCA-AS1, and shed light on the molecular pathogenesis of RHD occurrence.

MYB mediates downregulation of the colorectal cancer metastasis suppressor heterogeneous nuclear ribonucleoprotein L-like during epithelial-mesenchymal transition

Heterogeneous nuclear ribonucleoprotein L-like (HNRNPLL), a suppressor of colorectal cancer (CRC) metastasis, is transcriptionally downregulated when CRC cells undergo epithelial-mesenchymal transition (EMT). Here we show that decrease of MYB mediates the downregulation of HNRNPLL during EMT. The promoter activity was attributed to a region from -273 to -10 base pairs upstream of the transcription start site identified by 5'-RACE analysis, and the region contained potential binding sites for MYB and SP1. Luciferase reporter gene assays and knockdown or knockout experiments for genes encoding the MYB family proteins, MYB, MYBL1, and MYBL2, revealed that MYB was responsible for approximately half of the promoter activity. On the other hand, treatment with mithramycin A, an inhibitor for SP1 and SP3, suppressed the promoter activity and their additive contribution was confirmed by knockout experiments. The expression level of MYB was reduced on EMT while that of SP1 and SP3 was unchanged, suggesting that the downregulation of HNRNPLL during EMT was mediated by the decrease of MYB expression while SP1 and SP3 determine the basal transcription level of HNRNPLL. Histopathological analysis confirmed the accumulation of MYB-downregulated cancer cells at the invasion front of clinical CRC tissues. These results provide an insight into the molecular mechanism underlying CRC progression.

Identification of Cancer Cell Stemness-Associated Long Noncoding RNAs for Predicting Prognosis of Patients with Hepatocellular Carcinoma

Long noncoding RNAs (lncRNAs) are emerging as crucial contributors to the development of hepatocellular carcinoma (HCC) and are involved in the stemness regulation of liver cancer stem cells (LCSCs). However, cancer cell stemness-associated lncRNAs and their relevance in prediction of clinical prognosis remain largely unexplored. In this study, through the transcriptome-wide screen, we identified a total of 136 LCSC-associated lncRNAs. We evaluated the prognostic value of these lncRNAs and optimally established an 11-lncRNA (including AC008622.2, AC015908.3, AC020915.2, AC025176.1, AC026356.2, AC099850.3, CYTOR, DDX11-AS1, HTR2A-AS1, LINC02870, and SNHG3) prognostic risk model. Multivariate analysis revealed that the risk score is an independent prognostic predictor for HCC patients, which outperforms the traditional clinical pathological factors. Gene set enrichment analysis suggested that the high-risk score reflects the alteration of pathways involved in cell cycle, oxidative phosphorylation, and metabolism. Furthermore, functional studies on SNHG12, the leading candidate of the risk lncRNAs, revealed that knockdown of SNHG12 reduces the abilities of HCC cells stemness, proliferation, migration, and invasion. In summary, we constructed a prognostic risk model based on 11 LCSC-associated lncRNAs, which might be a promising prognostic predictor for HCC patients and highlight the involvement of lncRNAs in LCSC-associated treatment strategy in clinical practice.

Snail-Family Proteins: Role in Carcinogenesis and Prospects for Antitumor Therapy

The review analyzes Snail family proteins, which are transcription factors involved in the regulation of the epithelial-mesenchymal transition (EMT) of tumor cells. We describe the structure of these proteins, their post-translational modification, and the mechanisms of Snail-dependent regulation of genes. The role of Snail proteins in carcinogenesis, invasion, and metastasis is analyzed. Furthermore, we focus on EMT signaling mechanisms involving Snail proteins. Next, we dissect Snail signaling in hypoxia, a condition that complicates anticancer treatment. Finally, we offer classes of chemical compounds capable of down-regulating the transcriptional activity of Snails. Given the important role of Snail proteins in cancer biology and the potential for pharmacological inhibition, Snail family proteins may be considered promising as therapeutic targets.

TGFβ promotes widespread enhancer chromatin opening and operates on genomic regulatory domains

The Transforming Growth Factor-β (TGFβ) signaling pathway controls transcription by regulating enhancer activity. How TGFβ-regulated enhancers are selected and what chromatin changes are associated with TGFβ-dependent enhancers regulation are still unclear. Here we report that TGFβ treatment triggers fast and widespread increase in chromatin accessibility in about 80% of the enhancers of normal mouse mammary epithelial-gland cells, irrespective of whether they are activated, repressed or not regulated by TGFβ. This enhancer opening depends on both the canonical and non-canonical TGFβ pathways. Most TGFβ-regulated genes are located around enhancers regulated in the same way, often creating domains of several co-regulated genes that we term TGFβ regulatory domains (TRD). CRISPR-mediated inactivation of enhancers within TRDs impairs TGFβ-dependent regulation of all co-regulated genes, demonstrating that enhancer targeting is more promiscuous than previously anticipated. The area of TRD influence is restricted by topologically associating domains (TADs) borders, causing a bias towards co-regulation within TADs.

Radiosensitivity-Related Genes and Clinical Characteristics of Nasopharyngeal Carcinoma

Materials and Methods

Clinicopathological data of 185 patients with NPC treated at Nanfang Hospital of Southern Medical University between January 2013 and December 2014 were retrospectively analyzed. SPSS statistical software was used to analyze the clinicopathological data related to radiotherapy efficacy. Three patients who achieved complete remission and three with disease progression after CRT were selected. Differentially expressed genes (DEGs) were screened via mRNA microarray analysis of primary diagnostic endoscopy specimens.

Results

The peripheral blood leukocyte count, platelet count, and EBV-DNA copy number in NPC patients who were resistant to radiotherapy were higher than those in NPC patients who were sensitive to radiotherapy. The RobustRankAggreg (RRA) analysis method identified 392 DEGs, and the 66 most closely related genes among the DEGs were identified from the PPI network.

Conclusion

The results of this study indicate that screening for DEGs and pathways in NPC using integrated in silico analyses can help identify a series of genetic and clinical signatures for NPC patients treated with neoadjuvant chemotherapy followed by concurrent chemoradiotherapy.

NGS-Integrator: An efficient tool for combining multiple NGS data tracks using minimum Bayes' factors

BACKGROUND

Next-generation sequencing (NGS) is widely used for genome-wide identification and quantification of DNA elements involved in the regulation of gene transcription. Studies that generate multiple high-throughput NGS datasets require data integration methods for two general tasks: 1) generation of genome-wide data tracks representing an aggregate of multiple replicates of the same experiment; and 2) combination of tracks from different experimental types that provide complementary information regarding the location of genomic features such as enhancers.

RESULTS

NGS-Integrator is a Java-based command line application, facilitating efficient integration of multiple genome-wide NGS datasets. NGS-Integrator first transforms all input data tracks using the complement of the minimum Bayes' factor so that all values are expressed in the range [0,1] representing the probability of a true signal given the background noise. Then, NGS-Integrator calculates the joint probability for every genomic position to create an integrated track. We provide examples using real NGS data generated in our laboratory and from the mouse ENCODE database.

CONCLUSIONS

Our results show that NGS-Integrator is both time- and memory-efficient. Our examples show that NGS-Integrator can integrate information to facilitate downstream analyses that identify functional regulatory domains along the genome.

Long noncoding RNA MAPKAPK5-AS1 promotes colorectal cancer progression by cis-regulating the nearby gene MK5 and acting as a let-7f-1-3p sponge

BACKGROUND

Long noncoding RNAs (lncRNAs) are considered critical regulators in cancers; however, the clinical significance and mechanisms of MAPKAPK5-AS1 (hereinafter referred to as MK5-AS1) in colorectal cancer (CRC) remain mostly unknown.

METHODS

In this study, quantitative real-time PCR (qPCR) and western blotting were utilized to detect the levels of MK5-AS1, let-7f-1-3p and MK5 (MAPK activated protein kinase 5) in CRC tissues and cell lines. The biological functions of MK5-AS1, let-7f-1-3p and MK5 in CRC cells were explored using Cell Counting Kit-8 (CCK8), colony formation and transwell assays. The potential mechanisms of MK5-AS1 were evaluated by RNA pull-down, RNA immunoprecipitation (RIP), dual luciferase reporter assay, chromatin immunoprecipitation (ChIP) and bioinformatics analysis. The effects of MK5-AS1 and MK5 on CRC were investigated by a xenotransplantation model.

RESULTS

We confirmed that MK5-AS1 was significantly increased in CRC tissues. Knockdown of MK5-AS1 suppressed cell migration and invasion in vitro and inhibited lung metastasis in mice. Mechanistically, MK5-AS1 regulated SNAI1 expression by sponging let-7f-1-3p and cis-regulated the adjacent gene MK5. Moreover, MK5-AS1 recruited RBM4 and eIF4A1 to promote the translation of MK5. Our study verified that MK5 promoted the phosphorylation of c-Jun, which activated the transcription of SNAI1 by directly binding to its promoter.

CONCLUSIONS

MK5-AS1 cis-regulated the nearby gene MK5 and acted as a let-7f-1-3p sponge, playing a vital role in CRC tumorigenesis. This study could provide novel insights into molecular therapeutic targets of CRC.

Translation regulatory long non-coding RNA 1 represents a potential prognostic biomarker for colorectal cancer

Long non-coding RNAs (lncRNAs) have attracted a lot of attention for their role in the development, progression and prognosis of colorectal cancer (CRC). However, little is known on the clinical significance of the translation regulatory lncRNA 1 (TRERNA1) in CRC. The present study aimed to explore the clinical value of TRERNA1 in patients with CRC. A total of 89 cancer-associated lncRNA genes were analyzed using the RT² lncRNA PCR array Human Cancer PathwayFinder. Following the PCR array, reverse transcription-quantitative (RT-q)PCR was conducted to identify the differential expression of TRERNA1 between 130 CRC and corresponding non-tumorous adjacent tissues. Additionally, the association between TRERNA1 expression and clinical characteristics was analyzed. Furthermore, TRERNA1 expression was knocked down via small interfering RNAs. The results of the PCR array and RT-qPCR revealed that TRERNA1 expression was significantly upregulated in CRC tissues compared with in adjacent normal tissues. TRERNA1 upregulation was positively associated with distant metastasis, perineural invasion, TNM stage, node metastasis stage and tumor diameter. Multivariate analysis revealed that patients with higher TRERNA1 expression had a shorter overall survival (OS) time and a less favorable prognosis compared with those in the low TRERNA1 expression group. Knockdown of TRERNA1 inhibited invasion and metastasis of CRC cells via regulating Snail expression. In conclusion, TRERNA1 expression was upregulated in CRC tissues. High expression levels of TRERNA1 may be associated with poor OS times, a less favorable prognosis and lymph node metastasis in patients with CRC. TRERNA1 may therefore serve as a useful and novel biomarker for CRC lymph node metastasis and prognosis.

Loss of the nuclear Wnt pathway effector TCF7L2 promotes migration and invasion of human colorectal cancer cells

The transcription factor TCF7L2 is indispensable for intestinal tissue homeostasis where it transmits mitogenic Wnt/β-Catenin signals in stem and progenitor cells, from which intestinal tumors arise. Yet, TCF7L2 belongs to the most frequently mutated genes in colorectal cancer (CRC), and tumor-suppressive functions of TCF7L2 were proposed. This apparent paradox warrants to clarify the role of TCF7L2 in colorectal carcinogenesis. Here, we investigated TCF7L2 dependence/independence of CRC cells and the cellular and molecular consequences of TCF7L2 loss-of-function. By genome editing we achieved complete TCF7L2 inactivation in several CRC cell lines without loss of viability, showing that CRC cells have widely lost the strict requirement for TCF7L2. TCF7L2 deficiency impaired G1/S progression, reminiscent of the physiological role of TCF7L2. In addition, TCF7L2-negative cells exhibited morphological changes, enhanced migration, invasion, and collagen adhesion, albeit the severity of the phenotypic alterations manifested in a cell-line-specific fashion. To provide a molecular framework for the observed cellular changes, we performed global transcriptome profiling and identified gene-regulatory networks in which TCF7L2 positively regulates the proto-oncogene MYC, while repressing the cell cycle inhibitors CDKN2C/CDKN2D. Consistent with its function in curbing cell motility and invasion, TCF7L2 directly suppresses the pro-metastatic transcription factor RUNX2 and impinges on the expression of cell adhesion molecules. Altogether, we conclude that the proliferation-stimulating activity of TCF7L2 persists in CRC cells. In addition, TCF7L2 acts as invasion suppressor. Despite its negative impact on cell cycle progression, TCF7L2 loss-of-function may thereby increase malignancy, which could explain why TCF7L2 is mutated in a sizeable fraction of colorectal tumors.

Canonical BMP Signaling Executes Epithelial-Mesenchymal Transition Downstream of SNAIL1

Epithelial-mesenchymal transition (EMT) is a pivotal process in development and disease. In carcinogenesis, various signaling pathways are known to trigger EMT by inducing the expression of EMT transcription factors (EMT-TFs) like SNAIL1, ultimately promoting invasion, metastasis and chemoresistance. However, how EMT is executed downstream of EMT-TFs is incompletely understood. Here, using human colorectal cancer (CRC) and mammary cell line models of EMT, we demonstrate that SNAIL1 critically relies on bone morphogenetic protein (BMP) signaling for EMT execution. This activity requires the transcription factor SMAD4 common to BMP/TGFβ pathways, but is TGFβ signaling-independent. Further, we define a signature of BMP-dependent genes in the EMT-transcriptome, which orchestrate EMT-induced invasiveness, and are found to be regulated in human CRC transcriptomes and in developmental EMT processes. Collectively, our findings substantially augment the knowledge of mechanistic routes whereby EMT can be effectuated, which is relevant for the conceptual understanding and therapeutic targeting of EMT processes.

Non-coding RNAs and nuclear architecture during epithelial-mesenchymal transition in lung cancer and idiopathic pulmonary fibrosis

Lung cancer (LC) is the leading cause of cancer-related deaths worldwide. On the other hand, idiopathic pulmonary fibrosis (IPF) is the most common interstitial lung disease showing a prevalence of 20 new cases per 100,000 persons per year. Despite differences in cellular origin and pathological phenotypes, LC and IPF are lung diseases that share common features, including hyperproliferation of specific cell types in the lung, involvement of epithelial-mesenchymal transition (EMT) and enhanced activity of signaling pathways, such as tissue growth factor (TGFB), epidermal growth factor (EGF), fibroblast growth factor (FGF), wingless secreted glycoprotein (WNT) signaling, among others. EMT is a process during which epithelial cells lose their cell polarity and cell-cell adhesion, and acquire migratory and invasive properties to become mesenchymal cells. EMT involves numerous morphological hallmarks of hyperproliferative diseases, like cell plasticity, resistance to apoptosis, dedifferentiation and proliferation, thereby playing a central role during organ fibrosis and cancer progression. EMT was considered as an "all-or-none" process. In contrast to these outdated dichotomist interpretations, recent reports suggest that EMT occurs gradually involving different epithelial cell intermediate states with mesenchyme-like characteristics. These cell intermediate states of EMT differ from each other in their cell plasticity, invasiveness and metastatic potential, which in turn are induced by signals from their microenvironment. EMT is regulated by several transcription factors (TFs), which are members of prominent families of master regulators of transcription. In addition, there is increasing evidence for the important contribution of noncoding RNAs (ncRNAs) to EMT. In our review we highlight articles dissecting the function of different ncRNAs subtypes and nuclear architecture in cell intermediate states of EMT, as well as their involvement in LC and IPF.