A bioinformatics approach reveals novel interactions of the OVOL transcription factors in the regulation of epithelial - mesenchymal cell reprogramming and cancer progression

Single-cell resolution of MET- and EMT-like programs in osteoblasts during zebrafish fin regeneration

Zebrafish regenerate fin rays following amputation through epimorphic regeneration, a process that has been proposed to involve the epithelial-to-mesenchymal transition (EMT). We performed single-cell RNA sequencing (scRNA-seq) to elucidate osteoblastic transcriptional programs during zebrafish caudal fin regeneration. We show that osteoprogenitors are enriched with components associated with EMT and its reverse, mesenchymal-to-epithelial transition (MET), and provide evidence that the EMT markers cdh11 and twist2 are co-expressed in dedifferentiating cells at the amputation stump at 1 dpa, and in differentiating osteoblastic cells in the regenerate, the latter of which are enriched in EMT signatures. We also show that esrp1, a regulator of alternative splicing in epithelial cells that is associated with MET, is expressed in a subset of osteoprogenitors during outgrowth. This study provides a single cell resource for the study of osteoblastic cells during zebrafish fin regeneration, and supports the contribution of MET- and EMT-associated components to this process.

•

Osteoblasts express EMT/MET signatures during zebrafish fin regeneration

•

De/re-differentiating osteoblasts express cdh11, an EMT marker

•

A subset of osteoprogenitors express the MET marker esrp1

•

Our scRNA-seq data can be explored online

Shavenbaby and Yorkie mediate Hippo signaling to protect adult stem cells from apoptosis

To compensate for accumulating damages and cell death, adult homeostasis (e.g., body fluids and secretion) requires organ regeneration, operated by long-lived stem cells. How stem cells can survive throughout the animal life remains poorly understood. Here we show that the transcription factor Shavenbaby (Svb, OvoL in vertebrates) is expressed in renal/nephric stem cells (RNSCs) of Drosophila and required for their maintenance during adulthood. As recently shown in embryos, Svb function in adult RNSCs further needs a post-translational processing mediated by the Polished rice (Pri) smORF peptides and impairing Svb function leads to RNSC apoptosis. We show that Svb interacts both genetically and physically with Yorkie (YAP/TAZ in vertebrates), a nuclear effector of the Hippo pathway, to activate the expression of the inhibitor of apoptosis DIAP1. These data therefore identify Svb as a nuclear effector in the Hippo pathway, critical for the survival of adult somatic stem cells.

Cellular responses to human cytomegalovirus infection: Induction of a mesenchymal-to-epithelial transition (MET) phenotype

Human cytomegalovirus (HCMV) is the prototypical human β-herpes virus. Here we perform a systems analysis of the HCMV host-cell transcriptome, using gene set enrichment analysis (GSEA) as an engine to globally map the host-pathogen interaction across two cell types. Our analysis identified several previously unknown signatures of infection, such as induction of potassium channels and amino acid transporters, derepression of genes marked with histone H3 lysine 27 trimethylation (H3K27me3), and inhibition of genes related to epithelial-to-mesenchymal transition (EMT). The repression of EMT genes was dependent on early viral gene expression and correlated with induction E-cadherin (CDH1) and mesenchymal-to-epithelial transition (MET) genes. Infection of transformed breast carcinoma and glioma stem cells similarly inhibited EMT and induced MET, arguing that HCMV induces an epithelium-like cellular environment during infection.

Discovery and validation of information theory-based transcription factor and cofactor binding site motifs

Data from ChIP-seq experiments can derive the genome-wide binding specificities of transcription factors (TFs) and other regulatory proteins. We analyzed 765 ENCODE ChIP-seq peak datasets of 207 human TFs with a novel motif discovery pipeline based on recursive, thresholded entropy minimization. This approach, while obviating the need to compensate for skewed nucleotide composition, distinguishes true binding motifs from noise, quantifies the strengths of individual binding sites based on computed affinity and detects adjacent cofactor binding sites that coordinate with the targets of primary, immunoprecipitated TFs. We obtained contiguous and bipartite information theory-based position weight matrices (iPWMs) for 93 sequence-specific TFs, discovered 23 cofactor motifs for 127 TFs and revealed six high-confidence novel motifs. The reliability and accuracy of these iPWMs were determined via four independent validation methods, including the detection of experimentally proven binding sites, explanation of effects of characterized SNPs, comparison with previously published motifs and statistical analyses. We also predict previously unreported TF coregulatory interactions (e.g. TF complexes). These iPWMs constitute a powerful tool for predicting the effects of sequence variants in known binding sites, performing mutation analysis on regulatory SNPs and predicting previously unrecognized binding sites and target genes.

Khan F, Schmitz U, Wolkenhauer O, Gupta SK. MicroRNA and Transcription Factor Gene Regulatory Network Analysis Reveals Key Regulatory Elements Associated with Prostate Cancer Progression

Technological and methodological advances in multi-omics data generation and integration approaches help elucidate genetic features of complex biological traits and diseases such as prostate cancer. Due to its heterogeneity, the identification of key functional components involved in the regulation and progression of prostate cancer is a methodological challenge. In this study, we identified key regulatory interactions responsible for primary to metastasis transitions in prostate cancer using network inference approaches by integrating patient derived transcriptomic and miRomics data into gene/miRNA/transcription factor regulatory networks. One such network was derived for each of the clinical states of prostate cancer based on differentially expressed and significantly correlated gene, miRNA and TF pairs from the patient data. We identified key elements of each network using a network analysis approach and validated our results using patient survival analysis. We observed that HOXD10, BCL2 and PGR are the most important factors affected in primary prostate samples, whereas, in the metastatic state, STAT3, JUN and JUNB are playing a central role. Benefiting integrative networks our analysis suggests that some of these molecules were targeted by several overexpressed miRNAs which may have a major effect on the dysregulation of these molecules. For example, in the metastatic tumors five miRNAs (miR-671-5p, miR-665, miR-663, miR-512-3p and miR-371-5p) are mainly responsible for the dysregulation of STAT3 and hence can provide an opportunity for early detection of metastasis and development of alternative therapeutic approaches. Our findings deliver new details on key functional components in prostate cancer progression and provide opportunities for the development of alternative therapeutic approaches.

Expression of Ovol2 is related to epithelial characteristics and shows a favorable clinical outcome in hepatocellular carcinoma

Metastasis involves epithelial–mesenchymal transition (EMT) and mesenchymal–epithelial transition. Ovol2 belongs to the Ovo-like family (Ovol) of evolutionarily conserved zinc-finger transcription factors that regulate gene expression in various differentiation processes. Recent studies have demonstrated that Ovols affect mesenchymal–epithelial transition by inducing the expression of miR-200 in a range of human cancers. Downregulated Ovol2 expression is involved in the invasion and metastasis of breast and prostate cancers, but little is known about its expression and prognostic value in other cancers, including hepatocellular carcinoma (HCC). This study was designed to explore the clinical and prognostic significance of Ovol2 in patients with HCC. The expression of Ovol2 in tumor samples from patients with HCC and HCC cell lines was examined using Western blotting, real-time polymerase chain reaction, and immunohistochemistry. The expression levels of EMT-related markers, including E-cadherin, N-cadherin, and vimentin, were detected in relation to Ovol2 expression. The prognostic significance of Ovol2 in patients with HCC was statistically analyzed by Kaplan–Meier and Cox regression analyses. Ovol2 expression was significantly lower in HCC tissues than in adjacent noncancerous tissues. Low expression of Ovol2 was detected in HCC tissues with poor histological differentiation, microvascular invasion, and cirrhosis. A significant relationship was observed between Ovol2 and EMT marker expression levels. Kaplan–Meier analysis showed that overall survival was significantly worse in patients with HCC with low Ovol2 expression, indicating that Ovol2 deletion was an independent predictor of unfavorable prognosis in patients with HCC. Elevated Ovol2 expression may suppress HCC cell invasion and metastasis via restricting EMT.

Coupling the modules of EMT and stemness: A tunable 'stemness window' model

Metastasis of carcinoma involves migration of tumor cells to distant organs and initiate secondary tumors. Migration requires a complete or partial Epithelial-to-Mesenchymal Transition (EMT), and tumor-initiation requires cells possessing stemness. Epithelial cells (E) undergoing a complete EMT to become mesenchymal (M) have been suggested to be more likely to possess stemness. However, recent studies suggest that stemness can also be associated with cells undergoing a partial EMT (hybrid E/M phenotype). Therefore, the correlation between EMT and stemness remains elusive. Here, using a theoretical framework that couples the core EMT and stemness modules (miR-200/ZEB and LIN28/let-7), we demonstrate that the positioning of ‘stemness window’ on the ‘EMT axis’ need not be universal; rather it can be fine-tuned. Particularly, we present OVOL as an example of a modulating factor that, due to its coupling with miR-200/ZEB/LIN28/let-7 circuit, fine-tunes the EMT-stemness interplay. Coupling OVOL can inhibit the stemness likelihood of M and elevate that of the hybrid E/M (partial EMT) phenotype, thereby pulling the ‘stemness window’ away from the M end of ‘EMT axis’. Our results unify various apparently contradictory experimental findings regarding the interconnection between EMT and stemness, corroborate the emerging notion that partial EMT associates with stemness, and offer new testable predictions.

OVOL guides the epithelial-hybrid-mesenchymal transition

Metastasis involves multiple cycles of Epithelial-to-Mesenchymal Transition (EMT) and its reverse-MET. Cells can also undergo partial transitions to attain a hybrid epithelial/mesenchymal (E/M) phenotype that has maximum cellular plasticity and allows migration of Circulating Tumor Cells (CTCs) as a cluster. Hence, deciphering the molecular players helping to maintain the hybrid E/M phenotype may inform anti-metastasis strategies. Here, we devised a mechanism-based mathematical model to couple the transcription factor OVOL with the core EMT regulatory network miR-200/ZEB that acts as a three-way switch between the E, E/M and M phenotypes. We show that OVOL can modulate cellular plasticity in multiple ways - restricting EMT, driving MET, expanding the existence of the hybrid E/M phenotype and turning both EMT and MET into two-step processes. Our theoretical framework explains the differences between the observed effects of OVOL in breast and prostate cancer, and provides a platform for investigating additional signals during metastasis.

Loss of signal transducer and activator of transcription 1 is associated with prostate cancer recurrence

STAT1 loss has previously been implicated in cell line studies to modify prostate cancer cell growth and survival, however the clinical significance of this has not previously been established. This study investigated if STAT1 loss was associated with patient outcome measures and the phenotypic consequence of STAT1 silencing. STAT1 expression was assessed in two patient cohorts with localised (n = 78) and advanced prostate cancer at initial diagnosis (n = 39) by immunohistochemistry (IHC). Impact of STAT1 silencing on prostate cancer cells lines was assessed using Cell Death detection ELISA, TLDA gene signature apoptosis arrays, WST-1 assay, xCELLigence system, clonogenic assay, and wound healing assay. In the localised patient cohort, low expression of STAT1 was associated with shorter time to disease recurrence (3.8 vs 7.3 years, P = 0.02) and disease specific survival (6.6 vs 9.3 years, P = 0.05). In the advanced patient cohort, low expression was associated with shorter time to disease recurrence (2.0 vs 3.9 years, P = 0.001). When STAT1 was silenced in PC3 cells (AR negative) and LNCaP cells (AR positive) silencing did not influence levels of apoptosis in either cell line and had little effect on cell viability in the LNCaP cells. In contrast, STAT1 silencing in the PC3 cells resulted in a pronounced increase in cell viability (WST-1 assay: mock silenced vs STAT1 silenced, P < 0.001), clonagenicity (clonogenic assay: mock silenced vs STAT1 silenced, P < 0.001), and migration (wound healing: mock silenced vs STAT1 silenced, P < 0.001). In conclusion, loss of STAT1 may promote prostate cancer recurrence in AR negative patients via increasing cell viability. © 2015 Wiley Periodicals, Inc.

Implications of the Hybrid Epithelial/Mesenchymal Phenotype in Metastasis

Transitions between epithelial and mesenchymal phenotypes – the epithelial to ­mesenchymal transition (EMT) and its reverse the mesenchymal to epithelial transition (MET) – are hallmarks of cancer metastasis. While transitioning between the epithelial and mesenchymal phenotypes, cells can also attain a hybrid epithelial/mesenchymal (E/M) (i.e., partial or intermediate EMT) phenotype. Cells in this phenotype have mixed epithelial (e.g., adhesion) and mesenchymal (e.g., migration) properties, thereby allowing them to move collectively as clusters. If these clusters reach the bloodstream intact, they can give rise to clusters of circulating tumor cells (CTCs), as have often been seen experimentally. Here, we review the operating principles of the core regulatory network for EMT/MET that acts as a “three-way” switch giving rise to three distinct phenotypes – E, M and hybrid E/M – and present a theoretical framework that can elucidate the role of many other players in regulating epithelial plasticity. Furthermore, we highlight recent studies on partial EMT and its association with drug resistance and tumor-initiating potential; and discuss how cell–cell communication between cells in a partial EMT phenotype can enable the formation of clusters of CTCs. These clusters can be more apoptosis-resistant and have more tumor-initiating potential than singly moving CTCs with a wholly mesenchymal (complete EMT) phenotype. Also, more such clusters can be formed under inflammatory conditions that are often generated by various therapies. Finally, we discuss the multiple advantages that the partial EMT or hybrid E/M phenotype have as compared to a complete EMT phenotype and argue that these collectively migrating cells are the primary “bad actors” of metastasis.