VGLL1 expression is associated with a triple-negative basal-like phenotype in breast cancer

Placental co-transcriptional activator Vestigial-like 1 (VGLL1) drives tumorigenesis via increasing transcription of proliferation and invasion genes

Introduction

Vestigial-like 1 (VGLL1) is a co-transcriptional activator that binds to TEA domain-containing transcription factors (TEADs). Its expression is upregulated in a variety of aggressive cancer types, including pancreatic and basal-like breast cancer, and increased transcription of VGLL1 is strongly correlated with poor prognosis and decreased overall patient survival. In normal tissues, VGLL1 is most highly expressed within placental trophoblast cells, which share the common attributes of rapid cellular proliferation and invasion with tumor cells. The impact of VGLL1 in cancer has not been fully elucidated and no VGLL1-targeted therapy currently exists.

Methods

The aim of this study was to evaluate the cellular function and downstream genomic targets of VGLL1 in placental, pancreatic, and breast cancer cells. Functional assays were employed to assess the role of VGLL1 in cellular invasion and proliferation, and ChIP-seq and RNAseq assays were performed to identify VGLL1 target genes and potential impact using pathway analysis.

Results

ChIP-seq analysis identified eight transcription factors with a VGLL1-binding motif that were common between all three cell types, including TEAD1-4, AP-1, and GATA6, and revealed ~3,000 shared genes with which VGLL1 interacts. Furthermore, increased VGLL1 expression led to an enhancement of cell invasion and proliferation, which was supported by RNAseq analysis showing transcriptional changes in several genes known to be involved in these processes.

Discussion

This work expands our mechanistic understanding of VGLL1 function in tumor cells and provides a strong rationale for developing VGLL1-targeted therapies for treating cancer patients.

VGLL1 stabilization of cytoplasmic TAZ promotes EGFR expression and maintains tumor initiating cells in breast cancer independent of TEAD

Vestigial-like family member 1 (VGLL1) is one of the X-linked genes whose expression is elevated in basal-like breast cancer (BLBC) because of X-chromosome isodisomy. As an approach towards understanding its function, we performed correlation study using transcript data of breast cancer patients from cBioPortal for Cancer Genomics. Our analysis identified EGFR as the most correlated transcript with VGLL1. We demonstrate that VGLL1 promotes EGFR expression and increases the frequency of breast tumor initiating cells (CD44high/+CD24low/-). These findings are crucial because an elevated EGFR expression and high frequency of CD44high/+CD24low/- cells are defining features of BLBC, and we provide a new mechanistic insight into how their expressions are controlled. Importantly, VGLL1 regulation of EGFR and CD44high/+CD24low/- population is mediated by the hippo-transducer TAZ which exerts its oncogenic roles by binding and activating TEAD transcription factors. A crucial finding is that TEAD-binding domain of TAZ is dispensable for its regulation of EGFR and CD44high/+CD24low/- cells. Instead, VGLL1 stabilization of cytoplasmic TAZ is essential for these functions. Also, we show that VGLL1/TAZ restricts the surface expression of CD24 which contributes to the increased number of CD44high/+CD24low/- cells. In addition, we observed that VGLL1 represses AXL expression and suppresses claudin-low phenotype, and that is caused by the VGLL1 mediated nuclear expulsion of TAZ. Therefore, EGFR and AXL seem to represent two different breast tumor subtypes, and their differential expressions is controlled by VGLL1.

Aggressive high-grade NF2 mutant meningiomas downregulate oncogenic YAP signaling via the upregulation of VGLL4 and FAT3/4

Meningiomas are the most common primary brain tumors in adults. Although generally benign, a subset of meningiomas is of higher grade, shows aggressive growth behavior and recurs even after multiple surgeries. Around half of all meningiomas harbor inactivating mutations in NF2. While benign low-grade NF2 mutant meningiomas exhibit few genetic events in addition to NF2 inactivation, aggressive high-grade NF2 mutant meningiomas frequently harbor a highly aberrant genome. We and others have previously shown that NF2 inactivation leads to YAP1 activation and that YAP1 acts as the pivotal oncogenic driver in benign NF2 mutant meningiomas. Using bulk and single-cell RNA-Seq data from a large cohort of human meningiomas, we show that aggressive NF2 mutant meningiomas harbor decreased levels YAP1 activity compared to their benign counterparts. Decreased expression levels of YAP target genes are significantly associated with an increased risk of recurrence. We then identify the increased expression of the YAP1 competitor VGLL4 as well as the YAP1 upstream regulators FAT3/4 as a potential mechanism for the downregulation of YAP activity in aggressive NF2 mutant meningiomas. High expression of these genes is significantly associated with an increased risk of recurrence. In vitro, overexpression of VGLL4 resulted in the downregulation of YAP activity in benign NF2 mutant meningioma cells, confirming the direct link between VGLL4 expression and decreased levels of YAP activity observed in aggressive NF2 mutant meningiomas. Our results shed new insight on the biology of benign and aggressive NF2 mutant meningiomas and may have important implications for the efficacy of therapies targeting oncogenic YAP1 activity in NF2 mutant meningiomas.

N6-methyladenosine-modified VGLL1 promotes ovarian cancer metastasis through high-mobility group AT-hook 1/Wnt/β-catenin signaling

The main causes of death in patients with ovarian cancer (OC) are invasive lesions and the spread of metastasis. The present study aimed to explore the mechanisms that might promote OC metastasis. Here, we identified that VGLL1 expression was remarkably increased in metastatic OC samples. The role of VGLL1 in OC metastasis and tumor growth was examined by cell function assays and mouse models. Mechanistically level, METTL3-mediated N6-methyladenosine (m6A) modification contributed to VGLL1 upregulation in an IGF2BP2 recognition-dependent manner. Furthermore, VGLL1 directly interacts with TEAD4 and co-transcriptionally activates HMGA1. HMGA1 further activates Wnt/β-catenin signaling to enhance OC metastasis by promoting the epithelial-mesenchyme transition traits. Rescue assays indicated that the upregulation of HMGA1 was essential for VGLL1-induced metastasis. Collectively, these findings showed that the m6A-induced VGLL1/HMGA1/β-catenin axis might play a vital role in OC metastasis and tumor growth. VGLL1 might serve as a prognostic marker and therapeutic target against the metastasis of OC.

ER Negative Breast Cancer and miRNA: There Is More to Decipher Than What the Pathologist Can See!

Breast cancer (BC), the most prevalent cancer in women, is a heterogenous disease. Despite advancements in BC diagnosis, prognosis, and therapeutics, survival rates have drastically decreased in the metastatic setting. Therefore, BC still remains a medical challenge. The evolution of high-throughput technology has highlighted gaps in the classification system of BCs. Of particular interest is the notorious triple negative BC, which was recounted as being heterogenous itself and it overlaps with distinct subtypes, namely molecular apocrine (MA) and luminal androgen (LAR) BCs. These subtypes are, even today, still misdiagnosed and poorly treated. As such, researchers and clinicians have been looking for ways through which to refine BC classification in order to properly understand the initiation, development, progression, and the responses to the treatment of BCs. One tool is biomarkers and, specifically, microRNA (miRNA), which are highly reported as associated with BC carcinogenesis. In this review, the diverse roles of miRNA in estrogen receptor negative (ER-) and androgen receptor positive (AR+) BC are depicted. While highlighting their oncogenic and tumor suppressor functions in tumor progression, we will discuss their diagnostic, prognostic, and predictive biomarker potentials, as well as their drug sensitivity/resistance activity. The association of several miRNAs in the KEGG-reported pathways that are related to ER-BC carcinogenesis is presented. The identification and verification of accurate miRNA panels is a cornerstone for tackling BC classification setbacks, as is also the deciphering of the carcinogenesis regulators of ER - AR + BC.

Local augmented graph neural network for multi-omics cancer prognosis prediction and analysis

Cancer prognosis prediction and analysis can help patients understand expected life and help clinicians provide correct therapeutic guidance. Thanks to the development of sequencing technology, multi-omics data, and biological networks have been used for cancer prognosis prediction. Besides, graph neural networks can simultaneously consider multi-omics features and molecular interactions in biological networks, becoming mainstream in cancer prognosis prediction and analysis. However, the limited number of neighboring genes in biological networks restricts the accuracy of graph neural networks. To solve this problem, a local augmented graph convolutional network named LAGProg is proposed in this paper for cancer prognosis prediction and analysis. The process follows: first, given a patient's multi-omics data features and biological network, the corresponding augmented conditional variational autoencoder generates features. Then, the generated augmented features and the original features are fed into a cancer prognosis prediction model to complete the cancer prognosis prediction task. The conditional variational autoencoder consists of two parts: encoder-decoder. In the encoding phase, an encoder learns the conditional distribution of the multi-omics data. As a generative model, a decoder takes the conditional distribution and the original feature as inputs to generate the enhanced features. The cancer prognosis prediction model consists of a two-layer graph convolutional neural network and a Cox proportional risk network. The Cox proportional risk network consists of fully connected layers. Extensive experiments on 15 real-world datasets from TCGA demonstrated the effectiveness and efficiency of the proposed method in predicting cancer prognosis. LAGProg improved the C-index values by an average of 8.5% over the state-of-the-art graph neural network method. Moreover, we confirmed that the local augmentation technique could enhance the model's ability to represent multi-omics features, improve the model's robustness to missing multi-omics features, and prevent the model's over-smoothing during training. Finally, based on genes identified through differential expression analysis, we discovered 13 prognostic markers highly associated with breast cancer, among which ten genes have been proved by literature review.

Bioinformatic miRNA-mRNAs Analysis Revels to miR-934 as a Potential Regulator of the Epithelial-Mesenchymal Transition in Triple-Negative Breast Cancer

Triple-negative breast cancer (TNBC) is one of the most aggressive subtypes of breast cancer and has the worst prognosis. In patients with TNBC tumors, the tumor cells have been reported to have mesenchymal features, which help them migrate and invade. Various studies on cancer have revealed the importance of microRNAs (miRNAs) in different biological processes of the cell in that aberrations, in their expression, lead to alterations and deregulations in said processes, giving rise to tumor progression and aggression. In the present work, we determined the miRNAs that are deregulated in the epithelial-mesenchymal transition process in breast cancer. We discovered that 25 miRNAs that regulate mesenchymal genes are overexpressed in patients with TNBC. We found that miRNA targets modulate different processes and pathways, such as apoptosis, FoxO signaling pathways, and Hippo. We also found that the expression level of miR-934 is specific to the molecular subtype of the triple-negative breast cancer and modulates a set of related epithelial-mesenchymal genes. We determined that miR-934 inhibition in TNBC cell lines inhibits the migratory abilities of tumor cells.

Vestigial-like 1 (VGLL1): An ancient co-transcriptional activator linking wing, placenta, and tumor development

Vestigial-like 1 (VGLL1) is a recently discovered driver of proliferation and invasion that is expressed in many aggressive human malignancies and is strongly associated with poor prognosis. The VGLL1 gene encodes for a co-transcriptional activator that shows intriguing structural similarity to key activators in the hippo pathway, providing important clues to its functional role. VGLL1 binds to TEADs in an analogous fashion to YAP1 but appears to activate a distinct set of downstream gene targets. In mammals, VGLL1 expression is found almost exclusively in placental trophoblasts, cells that share many hallmarks of cancer. Due to its role as a driver of tumor progression, VGLL1 has become a target of interest for potential anticancer therapies. In this review, we discuss VGLL1 from an evolutionary perspective, contrast its role in placental and tumor development, summarize the current knowledge of how signaling pathways can modulate VGLL1 function, and discuss potential approaches for targeting VGLL1 therapeutically.

Construction of a miRNA-Based Nomogram Model to Predict the Prognosis of Endometrial Cancer

Objective: To investigate the differential expression of microRNA (miRNA) in patients with endometrial cancer and its relationship with prognosis and survival. Method: We used The Cancer Genome Atlas (TCGA) database to analyze differentially expressed miRNAs in endometrial cancer tissues and adjacent normal tissues. In addition, we successfully screened out key microRNAs to build nomogram models for predicting prognosis and we performed survival analysis on the key miRNAs as well. Result: We identified 187 differentially expressed miRNAs, which includes 134 up-regulated miRNAs and 53 down-regulated miRNAs. Further univariate Cox regression analysis screened out 47 significantly differentially expressed miRNAs and selected 12 miRNAs from which the prognostic nomogram model for ECA patients by LASSO analysis was constructed. Survival analysis showed that high expression of hsa-mir-138-2, hsa-mir-548f-1, hsa-mir-934, hsa-mir-940, and hsa-mir-4758 as well as low-expression of hsa-mir-146a, hsa-mir-3170, hsa-mir-3614, hsa-mir-3616, and hsa-mir-4687 are associated with poor prognosis in EC patients. However, significant correlations between the expressions levels of has-mir-876 and hsa-mir-1269a and patients' prognosis are not found. Conclusion: Our study found that 12 significantly differentially expressed miRNAs might promote the proliferation, invasion, and metastasis of cancer cells by regulating the expression of upstream target genes, thereby affecting the prognosis of patients with endometrial cancer.

The opposite effects of VGLL1 and VGLL4 genes on granulosa cell proliferation and apoptosis of hen ovarian prehierarchical follicles

Transcription cofactors Vestigial like family (VGLL) members consisting of four homologs (VGLL1-4) are associated with cell growth and metastasis in mammals, among which VGLL1 gene has been documented to possess tumorigenic functions in various types of tumor, and VGLL4 acts as a new tumor suppressor; likewise several studies indicated that they potentially play a role in the regulation of ovary growth and function. However, the biological effects of chicken VGLL1 and VGLL4 on the proliferation, apoptosis, and steroidogenesis of the granulosa cells (GCs) during ovarian follicle development remain unknown now. This study found that VGLL1 and VGLL4 genes present divergent expression patterns of the transcripts in the GCs of various sized prehierarchical follicles (PFs) before follicle selection. Specific small interfering RNA (siRNA) was employed to elucidate the exact roles of VGLL1 and VGLL4 in regulating the PF development of the hen ovary. The results demonstrated that the mRNA expression levels of the steroidogenic-related enzyme steroidogenic acute regulatory protein (STAR) gene and the cell proliferation-related factors B-cell lymphoma-2 (BCL2), and cyclin D1 (CCND1) genes were significantly down-regulated in the cells with VGLL1 silence but remarkably up-regulated in the cells lacking VGLL4. Whereas the expression level of the cell apoptosis biomarker caspase-3 (CASP3) transcript was noticeably enhanced in the GCs without VGLL1 but significantly decreased in the GCs deprived of VGLL4. Further results showed that the siRNA-mediated silence of VGLL1 caused a significant increase in apoptosis with a reduction in the proliferation of GCs. Nevertheless, knockdown of VGLL4 resulted in a remarkable decrement in apoptosis but a memorable augment in proliferation of the GCs. Taken together, this study proved that VGLL1 promotes cell proliferation and steroidogenesis but inhibits apoptosis. In contrast, VGLL4 stimulates GC apoptosis while suppressing the GC proliferation and steroidogenesis in the hen ovarian follicles. We conluded that VGLL1 and VGLL4 affect oppositely the ovarian prehierarchical follicle development by the different regulatory manner in the GC proliferation and apoptosis of chicken ovary.

An Efficient Algorithm for the Detection of Outliers in Mislabeled Omics Data

High dimensionality and noise have made it difficult to detect related biomarkers in omics data. Through previous study, penalized maximum trimmed likelihood estimation is effective in identifying mislabeled samples in high-dimensional data with mislabeled error. However, the algorithm commonly used in these studies is the concentration step (C-step), and the C-step algorithm that is applied to robust penalized regression does not ensure that the criterion function is gradually optimized iteratively, because the regularized parameters change during the iteration. This makes the C-step algorithm runs very slowly, especially when dealing with high-dimensional omics data. The AR-Cstep (C-step combined with an acceptance-rejection scheme) algorithm is proposed. In simulation experiments, the AR-Cstep algorithm converged faster (the average computation time was only 2% of that of the C-step algorithm) and was more accurate in terms of variable selection and outlier identification than the C-step algorithm. The two algorithms were further compared on triple negative breast cancer (TNBC) RNA-seq data. AR-Cstep can solve the problem of the C-step not converging and ensures that the iterative process is in the direction that improves criterion function. As an improvement of the C-step algorithm, the AR-Cstep algorithm can be extended to other robust models with regularized parameters.

Discovery of a cryptic site at the interface 2 of TEAD - Towards a new family of YAP/TAZ-TEAD inhibitors

The Hippo pathway is involved in organ size control and tissue homeostasis by regulating cell growth, proliferation and apoptosis. It controls the phosphorylation of the transcription co-activator YAP (Yes associated protein) and TAZ (Transcriptional coactivator with PDZ-binding motif) in order to control their nuclear import and their interaction with TEAD (Transcriptional Enhanced Associated Domain). YAP, TAZ and TEADs are dysregulated in several cancers making YAP/TAZ-TEAD interaction a new emerging anti-cancer target. We report the synthesis of a set of trisubstituted pyrazoles which bind to hTEAD2 at the interface 2 revealing for the first time a cryptic pocket created by the movement of the phenol ring of Y382. Compound 6 disrupts YAP/TAZ-TEAD interaction in HEK293T cells and inhibits TEAD target genes and cell proliferation in MDA-MB-231 cells. Compound 6 is therefore the first inhibitor of YAP/TAZ-TEAD targeting interface 2. This molecule could serve with other pan-TEAD inhibitors such as interface 3 ligands, for the delineation of the relative importance of VGLL vs YAP/TAZ in a given cellular model.

A panel of miRNAs as prognostic markers for African-American patients with triple negative breast cancer

BACKGROUND

To investigate the global expression profile of miRNAs, their impact on cellular signaling pathways, and their association with poor prognostic parameters in African-American (AA) patients with triple negative breast cancer (TNBC).

METHODS

Twenty-five samples of AA TNBC patients were profiled for global miRNA expression and stratified considering three clinical-pathological parameters: tumor size, lymph node (LN), and recurrence (REC) status. Differential miRNA expression analysis was performed for each parameter, and their discriminatory power was determined by Receiver Operating Characteristic (ROC) curve analysis. KMplotter was assessed to determine the association of the miRNAs with survival, and functional enrichment analysis to determine the main affected pathways and miRNA/mRNA target interactions.

RESULTS

A panel of eight, 23 and 27 miRNAs were associated with tumor size, LN, and REC status, respectively. Combined ROC analysis of two (miR-2117, and miR-378c), seven (let-7f-5p, miR-1255b-5p, miR-1268b, miR-200c-3p, miR-520d, miR-527, and miR-518a-5p), and three (miR-1200, miR-1249-3p, and miR-1271-3p) miRNAs showed a robust discriminatory power based on tumor size (AUC = 0.917), LN (AUC = 0.945) and REC (AUC = 0.981) status, respectively. Enrichment pathway analysis revealed their involvement in proteoglycans and glycan and cancer-associated pathways. Eight miRNAs with deregulated expressions in patients with large tumor size, positive LN metastasis, and recurrence were significantly associated with lower survival rates. Finally, the construction of miRNA/mRNA networks based in experimentally validated mRNA targets, revealed nodes of critical cancer genes, such as AKT1, BCL2, CDKN1A, EZR and PTEN.

CONCLUSIONS

Altogether, our data indicate that miRNA deregulated expression is a relevant biological factor that can be associated with the poor prognosis in TNBC of AA patients, by conferring to their TNBC cells aggressive phenotypes that are reflected in the clinical characteristics evaluated in this study.

Weighted gene co-expression network analysis of the association between upregulated AMD1, EN1 and VGLL1 and the progression and poor prognosis of breast cancer

Breast cancer is the most prevalent malignancy among females, but the molecular mechanisms involved in its pathogenesis and progression have remained to be fully elucidated. The aim of the present study was to identify novel potential therapeutic targets for breast cancer. The dataset GSE76275 was downloaded from the Gene Expression Omnibus database and weighted gene co-expression network analysis (WGCNA) was performed to identify hub genes. Furthermore, the dataset GSE25055, containing gene expression data and clinical information, was downloaded to validate the expression and survival association of these hub genes. In addition, the datasets GSE25065 and GSE42568 were used to validate the association between hub gene expression levels and clinical features. Immunohistochemistry (IHC), reverse transcription-quantitative PCR, as well as proliferation, migration, invasion and apoptosis assays, were used to verify gene expression and function. A total of 4,052 genes were selected for WGCNA and 18 modules were established; the red module was identified as the key module, as it had a strong positive correlation with the tumor grade. Survival analyses of hub genes [S-adenosylmethionine decarboxylase proenzyme (AMD1), homeobox protein engrailed-1 (EN1) and vestigial-like protein (VGLL1)] indicated that higher levels of gene expression were associated with poor prognosis of patients with breast cancer. This association was based on survival analysis of GSE25055 using the Kaplan-Meier plotter tool. Expression validation revealed that the upregulation of hub genes was associated with advanced tumor grade and malignant molecular subtype (basal-like). IHC results from the Human Protein Atlas also demonstrated that protein expression levels of the hub genes were higher in tumor tissues compared with those in adjacent normal tissues. Furthermore, the expression levels of AMD1, EN1 and VGLL1 were strongly correlated with each other. These results demonstrated that AMD1 is highly expressed in breast cancer tissues and cells and AMD1 knockdown decreased the proliferation and metastatic potential, while increasing apoptosis of breast cancer cells. These results suggested that AMD1, EN1 and VGLL1 are likely to contribute to breast cancer progression and unfavorable prognosis.

Genome Instability-Derived Genes Are Novel Prognostic Biomarkers for Triple-Negative Breast Cancer

Background

Triple-negative breast cancer (TNBC) is an aggressive disease. Recent studies have identified genome instability-derived genes for patient outcomes. However, most of the studies mainly focused on only one or a few genome instability-related genes. Prognostic potential and clinical significance of genome instability-associated genes in TNBC have not been well explored.

Methods

In this study, we developed a computational approach to identify TNBC prognostic signature. It consisted of (1) using somatic mutations and copy number variations (CNVs) in TNBC to build a binary matrix and identifying the top and bottom 25% mutated samples, (2) comparing the gene expression between the top and bottom 25% samples to identify genome instability-related genes, and (3) performing univariate Cox proportional hazards regression analysis to identify survival-associated gene signature, and Kaplan–Meier, log-rank test, and multivariate Cox regression analyses to obtain overall survival (OS) information for TNBC outcome prediction.

Results

From the identified 111 genome instability-related genes, we extracted a genome instability-derived gene signature (GIGenSig) of 11 genes. Through survival analysis, we were able to classify TNBC cases into high- and low-risk groups by the signature in the training dataset (log-rank test p = 2.66e−04), validated its prognostic performance in the testing (log-rank test p = 2.45e−02) and Molecular Taxonomy of Breast Cancer International Consortium (METABRIC) (log-rank test p = 2.57e−05) datasets, and further validated the predictive power of the signature in five independent datasets.

Conclusion

The identified novel signature provides a better understanding of genome instability in TNBC and can be applied as prognostic markers for clinical TNBC management.

miR-934 promotes breast cancer metastasis by regulation of PTEN and epithelial-mesenchymal transition

Breast cancer (BC) is the most commonly diagnosed malignancy and the leading cause of cancer-related mortality among females. Over 90 % of the cases of death in BC patients are attributed to tumor cell metastasis. Therefore, it is urgently needed to investigate the molecular mechanisms of BC metastasis. The expression of miRNA in BC was evaluated by qRT-PCR and bioinformatics analysis. Clone formation, EdU assays, and subcutaneous xenograft model were used to test the growth of BC cells. Wound healing, transwell assays, and lung-metastasis model were used to explore the effect of miR-934 knockdown on cell metastasis. The miR-934 targets in BC were identified through bioinformatics analysis and luciferase reporter assays. The expression of protein was tested by western blot. The binding of mRNA and RNA-binding-protein was verified using RIP assays. miR-934 expression was significantly elevated in BC tissues, especially in those with lymph node metastasis and associated with poor patient prognosis. Experiments in vitro and in vivo showed that that upregulated miR-934 was not necessarily required for the growth of BC cells. However, miR-934 knockdown significantly inhibited the migration and invasion abilities of BC cells. Moreover, PTEN as identified as the direct target of miR-934 in BC, and miR-934 could promote BC cell metastasis by regulation of PTEN and epithelial-mesenchymal transition (EMT). Our results suggested that targeting miR-934 may be a practical treatment for BC cell metastasis.

Vestigial suppresses apoptosis and cell migration in a manner dependent on the level of JNK-Caspase signaling in the Drosophila wing disc

The Decapentaplegic (Dpp) and Wingless (Wg) signal pathways play important roles in numerous biological processes in Drosophila. The Drosophila vestigial (vg) gene is selectively required for wing imaginal disc cell proliferation, which is essential for the formation of the adult wing and halter structures, and is regulated by Dpp and Wg signaling. Using a Drosophila invasion model of wing epithelium, we showed herein that inhibition of Dpp or Wg signaling promoted cells to migrate across the cell lineage restrictive anterior/posterior (A/P) compartment boundary. Being downstream of both Dpp and Wg signaling, vg can block cell migration induced by loss of either pathway. In addition, suppression of vg is sufficient to induce cell migration across the A/P boundary. Transcriptomic analysis revealed potential downstream genes involved in the cell migration after suppressing vg in the wing disc. We further demonstrated that the c-Jun N-terminal kinase (JNK) signaling promoted cell migration induced by vg suppression by upregulating Caspase activity. Taken together, our results revealed the requirement of Vg for suppressing cell migration and clarified how developmental signals collaborate to stabilize cells along the compartment boundary.

Modulatory Role of microRNAs in Triple Negative Breast Cancer with Basal-Like Phenotype

Development of new research, classification, and therapeutic options are urgently required due to the fact that TNBC is a heterogeneous malignancy. The expression of high molecular weight cytokeratins identifies a biologically and clinically distinct subgroup of TNBCs with a basal-like phenotype, representing about 75% of TNBCs, while the remaining 25% includes all other intrinsic subtypes. The triple negative phenotype in basal-like breast cancer (BLBC) makes it unresponsive to endocrine therapy, i.e., tamoxifen, aromatase inhibitors, and/or anti-HER2-targeted therapies; for this reason, only chemotherapy can be considered an approach available for systemic treatment even if it shows poor prognosis. Therefore, treatment for these subgroups of patients is a strong challenge for oncologists due to disease heterogeneity and the absence of unambiguous molecular targets. Dysregulation of the cellular miRNAome has been related to huge cellular process deregulations underlying human malignancy. Consequently, epigenetics is a field of great promise in cancer research. Increasing evidence suggests that specific miRNA clusters/signatures might be of clinical utility in TNBCs with basal-like phenotype. The epigenetic mechanisms behind tumorigenesis enable progress in the treatment, diagnosis, and prevention of cancer. This review intends to summarize the epigenetic findings related to miRNAome in TNBCs with basal-like phenotype.

Vestigial-like 1 is a shared targetable cancer-placenta antigen expressed by pancreatic and basal-like breast cancers

Cytotoxic T lymphocyte (CTL)-based cancer immunotherapies have shown great promise for inducing clinical regressions by targeting tumor-associated antigens (TAA). To expand the TAA landscape of pancreatic ductal adenocarcinoma (PDAC), we performed tandem mass spectrometry analysis of HLA class I-bound peptides from 35 PDAC patient tumors. This identified a shared HLA-A*0101 restricted peptide derived from co-transcriptional activator Vestigial-like 1 (VGLL1) as a putative TAA demonstrating overexpression in multiple tumor types and low or absent expression in essential normal tissues. Here we show that VGLL1-specific CTLs expanded from the blood of a PDAC patient could recognize and kill in an antigen-specific manner a majority of HLA-A*0101 allogeneic tumor cell lines derived not only from PDAC, but also bladder, ovarian, gastric, lung, and basal-like breast cancers. Gene expression profiling reveals VGLL1 as a member of a unique group of cancer-placenta antigens (CPA) that may constitute immunotherapeutic targets for patients with multiple cancer types.

Sparse semiparametric canonical correlation analysis for data of mixed types

Canonical correlation analysis investigates linear relationships between two sets of variables, but often works poorly on modern datasets due to high-dimensionality and mixed data types (continuous/binary/zero-inflated). We propose a new approach for sparse canonical correlation analysis of mixed data types that does not require explicit parametric assumptions. Our main contribution is the use of truncated latent Gaussian copula to model the data with excess zeroes, which allows us to derive a rank-based estimator of latent correlation matrix without the estimation of marginal transformation functions. The resulting semiparametric sparse canonical correlation analysis method works well in high-dimensional settings as demonstrated via numerical studies, and application to the analysis of association between gene expression and micro RNA data of breast cancer patients.

Development of prediction models for the sensitivity of oral squamous cell carcinomas to preoperative S-1 administration

S-1 is an anticancer agent that is comprised of tegafur, gimeracil, and oteracil potassium, and is widely used in various carcinomas including oral squamous cell carcinoma (OSCC). Although an established prediction tool is not available, we aimed to develop prediction models for the sensitivity of primary OSCC cases to the preoperative administration of S-1.

We performed DNA microarray analysis of 95 cases with OSCC. Using global gene expression data and the clinical data, we developed two different prediction models, namely, model 1 that comprised the complete response (CR) + the partial response (PR) versus stable disease (SD) + progressive disease (PD), and model 2 that comprised responders versus non-responders. Twelve and 18 genes were designated as feature genes (FGs) in models 1 and 2, respectively, and, of these, six genes were common to both models. The sensitivity was 96.3%, the specificity was 91.2%, and the accuracy was 92.6% for model 1, and the sensitivity was 95.6%, the specificity was 85.2%, and the accuracy was 92.6% for model 2. These models were validated using receiver operating characteristic analysis, and the areas under the curves were 0.967 and 0.949 in models 1 and 2, respectively. The data led to the development of models that can reliably predict the sensitivity of patients with OSCC to the preoperative administration of S-1. The mechanism that regulates S-1 sensitivity remains unclear; however, the prediction models developed provide hope that further functional investigations into the FGs will lead to a greater understanding of drug resistance.

Comparison of methods for the detection of outliers and associated biomarkers in mislabeled omics data

Background

Previous studies have reported that labeling errors are not uncommon in omics data. Potential outliers may severely undermine the correct classification of patients and the identification of reliable biomarkers for a particular disease. Three methods have been proposed to address the problem: sparse label-noise-robust logistic regression (Rlogreg), robust elastic net based on the least trimmed square (enetLTS), and Ensemble. Ensemble is an ensembled classification based on distinct feature selection and modeling strategies. The accuracy of biomarker selection and outlier detection of these methods needs to be evaluated and compared so that the appropriate method can be chosen.

Results

The accuracy of variable selection, outlier identification, and prediction of three methods (Ensemble, enetLTS, Rlogreg) were compared for simulated and an RNA-seq dataset. On simulated datasets, Ensemble had the highest variable selection accuracy, as measured by a comprehensive index, and lowest false discovery rate among the three methods. When the sample size was large and the proportion of outliers was ≤5%, the positive selection rate of Ensemble was similar to that of enetLTS. However, when the proportion of outliers was 10% or 15%, Ensemble missed some variables that affected the response variables.

Overall, enetLTS had the best outlier detection accuracy with false positive rates < 0.05 and high sensitivity, and enetLTS still performed well when the proportion of outliers was relatively large. With 1% or 2% outliers, Ensemble showed high outlier detection accuracy, but with higher proportions of outliers Ensemble missed many mislabeled samples. Rlogreg and Ensemble were less accurate in identifying outliers than enetLTS. The prediction accuracy of enetLTS was better than that of Rlogreg. Running Ensemble on a subset of data after removing the outliers identified by enetLTS improved the variable selection accuracy of Ensemble.

Conclusions

When the proportion of outliers is ≤5%, Ensemble can be used for variable selection. When the proportion of outliers is > 5%, Ensemble can be used for variable selection on a subset after removing outliers identified by enetLTS. For outlier identification, enetLTS is the recommended method. In practice, the proportion of outliers can be estimated according to the inaccuracy of the diagnostic methods used.

Multiple Roles of Vestigial-Like Family Members in Tumor Development

Vestigial-like family (VGLL) members are mammalian orthologs of vestigial gene in Drosophila, and they consist of four homologs (VGLL1–4). VGLL members have TDU motifs that are binding regions to TEA/ATSS-DNA-binding domain transcription factor (TEAD). Through TDU motifs, VGLL members act as transcriptional cofactors for TEAD. VGLL1-3 have single TDU motif, whereas VGLL4 has two tandem TDU motifs, suggesting that VGLL4 has distinct molecular functions among this family. Although molecular and physiological functions of VGLL members are still obscure, emerging evidence has shown that these members are involved in tumor development. Gene alterations and elevated expression of VGLL1-3 were observed in various types of tumors, and VGLL1-3 have been shown to possess tumorigenic functions. In contrast, down-regulation of VGLL4 was detected in various tumors, and the tumor-suppressing role of VGLL4 has been demonstrated. In this review, we summarize the recently identified multiple roles of VGLL members in tumor development and provide important and novel insights regarding tumorigenesis.

The Transcriptional Cofactor VGLL1 Drives Transcription of Human Papillomavirus Early Genes via TEAD1

The TEAD family of transcription factors requires associating cofactors to induce gene expression. TEAD1 is known to activate the early promoter of human papillomavirus (HPV), but the precise mechanisms of TEAD1-mediated transactivation of the HPV promoter, including its relevant cofactors, remain unexplored. Here, we reveal that VGLL1, a TEAD-interacting cofactor, contributes to HPV early gene expression. Knockdown of VGLL1 and/or TEAD1 led to a decrease in viral early gene expression in human cervical keratinocytes and cervical cancer cell lines. We identified 11 TEAD1 target sites in the HPV16 long control region (LCR) by in vitro DNA pulldown assays; 8 of these sites contributed to the transcriptional activation of the early promoter in luciferase reporter assays. VGLL1 bound to the HPV16 LCR via its interaction with TEAD1 both in vitro and in vivo Furthermore, introducing HPV16 and HPV18 whole genomes into primary human keratinocytes led to increased levels of VGLL1, due in part to the upregulation of TEADs. These results suggest that multiple VGLL1/TEAD1 complexes are recruited to the LCR to support the efficient transcription of HPV early genes.IMPORTANCE Although a number of transcription factors have been reported to be involved in HPV gene expression, little is known about the cofactors that support HPV transcription. In this study, we demonstrate that the transcriptional cofactor VGLL1 plays a prominent role in HPV early gene expression, dependent on its association with the transcription factor TEAD1. Whereas TEAD1 is ubiquitously expressed in a variety of tissues, VGLL1 displays tissue-specific expression and is implicated in the development and differentiation of epithelial lineage tissues, where HPV gene expression occurs. Our results suggest that VGLL1 may contribute to the epithelial specificity of HPV gene expression, providing new insights into the mechanisms that regulate HPV infection. Further, VGLL1 is also critical for the growth of cervical cancer cells and may represent a novel therapeutic target for HPV-associated cancers.

miR-934 as a Prognostic Marker Facilitates Cell Proliferation and Migration of Pancreatic Tumor by Targeting PROX1

Background

Pancreatic cancer is an extremely lethal digestive cancer with late diagnosis and poor prognosis. miR-934 has been reported to serve as an oncogene in multiple cancers, such as ovarian cancer and bladder cancer. However, its role in pancreatic cancer remains undiscovered.

Materials and Methods

The expression data of miR-934 were obtained from the Gene Expression Omnibus database and from our own patient samples. The clinicopathological data and corresponding follow-up data were retrieved from The Cancer Genome Atlas database. CCK8 and colony formation assays were conducted to measure cell proliferation capacity in vitro. Wound healing and transwell assays were performed to detect the migration ability of pancreatic cancer cell.

Results

We found that miR-934 was significantly upregulated in pancreatic tumor samples and cell lines. The expression of miR-934 was related to pathological stages. Upregulated miR-934 was associated with poor prognosis in patients with pancreatic cancer. Mir-934 inhibition reduced, while overexpression promoted, cell proliferation and migration. Mechanically, we found  miR-934 could directly bind to 3'-UTR of PROX1 leading to mRNA derogation. Furthermore, increased cell proliferation and migration caused by miR-934 overexpression could be reversed by forced PROX1 expression.

Conclusion

miR-934 is an oncogene in pancreatic cancer and could serve as a prognosis indicator for patients with pancreatic cancer, suggesting that miR-934 is a promising therapeutic target for pancreatic cancer.

USF1-mediated upregulation of lncRNA GAS6-AS2 facilitates osteosarcoma progression through miR-934/BCAT1 axis

Long noncoding RNAs (lncRNAs) have been certified as important regulators in tumorigenesis. LncRNA GAS6-AS2 (GAS6-AS2) was a newly identified tumor-related lncRNA, and its dysregulation and oncogenic effects in melanoma and bladder cancer had been reported in previous studies. However, the expression pattern and potential function of GAS6-AS2 in osteosarcoma (OS) have not been investigated. In this study, we identified a novel OS-related lncRNA GAS6-AS2. We found that GAS6-AS2 was distinctly upregulated in both OS specimens and cell lines. Distinct up-regulation of GAS6-AS2 in OS was correlated with advanced clinical stages and shorter survivals. In addition, USF1 could directly bind to the GAS6-AS2 promoter and contribute to its overexpression. Furthermore, GAS6-AS2 knockdown caused tumor suppressive effects via reducing cellular proliferation, migration and invasion, and promoting OS cell apoptosis. Besides, GAS6-AS2 directly bound to miR-934 and downregulated its expression. Mechanistically, GAS6-AS2 positively regulated the expression of BCAT1 through sponging miR-934. Taken together, our data illustrated how GAS6-AS2 played an oncogenic role in OS and might offer a potential therapeutic target for treating OS.

PI3K/AKT/β-Catenin Signaling Regulates Vestigial-Like 1 Which Predicts Poor Prognosis and Enhances Malignant Phenotype in Gastric Cancer

Although gastric cancer is a common cause of cancer mortality worldwide, its biological heterogeneity limits the available therapeutic options. Therefore, identifying novel therapeutic targets for developing effective targeted therapy of gastric cancer is a pressing need. Here, we investigate molecular function and regulatory mechanisms of Vestigial-like 1 (VGLL1) in gastric cancer. Microarray analysis of 556 gastric cancer tissues revealed that VGLL1 was a prognostic biomarker that correlated with PI3KCA and PI3KCB. VGLL1 regulates the proliferation of gastric cancer cells, as shown in live cell imaging, sphere formation, and in vivo xenograft model. Tail vein injection of NUGC3 cells expressing shVGLL1 resulted in less lung metastasis occurring when compared to the control. In contrast, larger metastatic lesions in lung and liver were detected in the VGLL1-overexpressing NUGC3 cell xenograft excision mouse model. Importantly, VGLL1 expression is transcriptionally regulated by the PI3K-AKT-β-catenin pathway. Subsequently, MMP9, a key molecule in gastric cancer, was explored as one of target genes that were transcribed by VGLL1-TEAD4 complex, a component of the transcription factor. Taken together, PI3K/AKT/β-catenin signaling regulates the transcription of VGLL1, which promotes the proliferation and metastasis in gastric cancer. This finding suggests VGLL1 as a novel prognostic biomarker and a potential therapeutic target.

Zebrafish Vestigial Like Family Member 4b Is Required for Valvulogenesis Through Sequestration of Transcription Factor Myocyte Enhancer Factor 2c

A variety of cardiac transcription factors/cofactors, signaling pathways, and downstream structural genes integrate to form the regulatory hierarchies to ensure proper cardiogenesis in vertebrate. Major interaction proteins of the transcription cofactor vestigial like family member 4 (VGLL4) include myocyte enhancer factor 2 (MEF2) and TEA domain transcription factors (TEAD), both of which play important roles in embryonic cardiac development and in adulthood. In this study, we identified that the deficiency of zebrafish vgll4b paralog, a unique family member expressed in developing heart, led to an impaired valve development. Mechanistically, in vgll4b mutant embryos the disruption of Vgll4b-Mef2c complex, rather than that of Vgll4b-Tead complex, resulted in an aberrant expression of krüppel-like factor 2a (klf2a) in endocardium. Such misexpression of klf2a eventually evoked the valvulogenesis defects. Our findings suggest that zebrafish Vgll4b plays an important role in modulating the transcription activity of Mef2c on klf2a during valve development in a blood-flow-independent manner.

Differences in microRNA expression in breast cancer between women of African and European ancestry

Breast cancer is a heterogeneous disease, characterized by molecularly and phenotypically distinct tumor subtypes, linked to disparate clinical outcomes. American women of African ancestry (AA) are more likely than those of European ancestry (EA) to be diagnosed with aggressive, estrogen receptor negative (ER-) or triple negative breast cancer, and to die of this disease. However, the underlying causes of AA predisposition to ER-/triple negative breast cancer are still largely unknown. In this study, we performed high-throughput whole-genome miRNA expression profiling in breast tissue samples from both AA and EA women. A number of differentially expressed miRNAs, i.e., DEmiRs defined as >2-fold change in expression and false discovery rate <0.05, were identified as up- or downregulated by tumor ER status or by ancestry. We found that among 102 ER-subtype-related DEmiRs identified in breast tumors, the majority of these DEmiRs were race specific, with only 23 DEmiRs shared in tumors from both AAs and EAs; this finding indicates that there are unique subsets of miRNAs differentially expressed between ER- and ER positive tumors within AAs versus EAs. Our overall results support the notion that miRNA expression patterns may differ not only by tumor subtype but by ancestry, indicating differences in tumor biology and heterogeneity of breast cancer between AAs and EAs. These results will provide the basis for further functional analysis to elucidate biological differences between AAs and EAs and to help develop targeted treatment strategies to reduce disparities in breast cancer.

Coupled Genome-Wide DNA Methylation and Transcription Analysis Identified Rich Biomarkers and Drug Targets in Triple-Negative Breast Cancer

Triple-negative breast cancer (TNBC) has poor clinical prognosis. Lack of TNBC-specific biomarkers prevents active clinical intervention. We reasoned that TNBC must have its specific signature due to the lack of three key receptors to distinguish TNBC from other types of breast cancer. We also reasoned that coupling methylation and gene expression as a single unit may increase the specificity for the detected TNBC signatures. We further reasoned that choosing the proper controls may be critical to increasing the sensitivity to identify TNBC-specific signatures. Furthermore, we also considered that specific drugs could target the detected TNBC-specific signatures. We developed a system to identify potential TNBC signatures. It consisted of (1) coupling methylation and expression changes in TNBC to identify the methylation-regulated signature genes for TNBC; (2) using TPBC (triple-positive breast cancer) as the control to detect TNBC-specific signature genes; (3) searching in the drug database to identify those targeting TNBC signature genes. Using this system, we identified 114 genes with both altered methylation and expression, and 356 existing drugs targeting 10 of the 114 genes. Through docking and molecular dynamics simulation, we determined the structural basis between sapropterin, a drug used in the treatment of tetrahydrobiopterin deficiency, and PTGS2, a TNBC signature gene involved in the conversion of arachidonic acid to prostaglandins. Our study reveals the existence of rich TNBC-specific signatures, and many can be drug target and biomarker candidates for clinical applications.

Inhibition of UBE2N-dependent CDK6 protein degradation by miR-934 promotes human bladder cancer cell growth

Because bladder cancer (BC) is one of the most common malignant cancers of the urinary system, identification of BC cell growth-associated effectors is of great significance. Cyclin-dependent kinase (CDK)6 is a member of the CDK family of cell cycle-related proteins and plays an important role in cancer cell growth. This is borne out by the fact that a CDK6 inhibitor had been approved to treat several types of cancers. Nevertheless, underlying molecular mechanisms concerning how to regulate CDK6 expression in BC remains unclear. In the present study, it was observed that miR-934 was much higher in human BCs and human BC cell lines as well. The results also revealed that miR-934 inhibition dramatically decreased human BC cell monolayer growth in vitro and xenograft tumor growth in vivo; the outcomes were accompanied by CDK6 protein down-regulation and G₀-G₁ cell cycle arrest. Moreover, overexpression of CDK6 reversed the inhibition of BC cell growth induced by miR-934. Further studies showed that miR-934 binds to a 3'-UTR of ubiquitin-conjugating enzyme 2N (ube2n) mRNA, down-regulated UBE2N protein expression; this, in turn, attenuated CDK6 protein degradation and led to CDK6 protein accumulation as well as the promotion of BC tumor growth. Collectively, this study not only establishes a novel regulatory axis of miR-934/UBE2N of CDK6 but also provides data suggesting that miR-934 and UBE2N may be potentially promising targets for therapeutic strategies against BC.-Yan, H., Ren, S., Lin, Q., Yu, Y., Chen, C., Hua, X., Jin, H., Lu, Y., Zhang, H., Xie, Q., Huang, C., Huang, H. Inhibition of UBE2N-dependent CDK6 protein degradation by miR-934 promotes human bladder cancer cell growth.

Potential Five-MicroRNA Signature Model for the Prediction of Prognosis in Patients with Wilms Tumor

Background

Wilms tumor (WT) is the most common type of pediatric renal malignancy, and is associated with poor prognosis. The aim of the present study was to identify microRNA (miRNA) signatures which might predict prognosis and categorize WTs into high- and low-risk subgroups.

Material/Methods

The miRNA expression profiles of WT patients and normal samples were obtained from the Therapeutically Applicable Research to Generate Effective Treatment database. Differentially expressed miRNAs between WT patients and normal samples were identified using the EdgeR package. Subsequently, correlations between differentially expressed miRNAs and the prognosis of overall survival were analyzed. Enrichment analyses for the targeted mRNAs were conducted via the Database for Annotation, Visualization, and Integration Discovery.

Results

A total of 154 miRNAs were identified as differentially expressed in WT. Of those, 18 miRNAs were associated with overall survival (P<0.05). A prognostic signature of 5 differentially expressed miRNAs (i.e., has-mir-149, has-mir-7112, has-mir-940, has-mir-1248, and has-mir-490) was constructed to classify the patients into high- and low-risk subgroups. The targeted mRNAs of these prognostic miRNAs were primarily enriched in Gene Ontology terms (i.e., protein autophosphorylation, protein dephosphorylation, and stress-activated MAPK cascade) and the Kyoto Encyclopedia of Genes and Genomes signaling pathways (i.e., MAPK, AMPK, and PI3K-Akt).

Conclusions

The 5-miRNA signature model might be useful in determining the prognosis of WT patients. As a promising prediction tool, this prognosis signature might serve as a potential biomarker for WT patients.

Amplification of the CD24 Gene Is an Independent Predictor for Poor Prognosis of Breast Cancer

CD24 is a glycosyl-phosphatidyl-inositol linked glycoprotein expressed in a broad range of cell types including cancer cells. Although it is overexpressed in nearly 70% of human cancers, copy number variation of the CD24 locus has not been reported for any cancer. Here, we analyzed the genomics, transcriptomics, and clinical data of 1082 breast cancer (BRCA) samples and other cancer samples from the clinically annotated genomic database, The Cancer Genome Atlas (TCGA). The GISTIC2 method was applied to stratify the CD24 copy number, and Cox regression was performed to compare hazard ratio (HR) of CD24 overexpression, amplification and other traditional prognosis features for overall survival (OS). Our data demonstrated that CD24 amplification strongly correlated with its mRNA overexpression as well as TP53 mutant, cancer proliferation and metastasis features. In particular, CD24 amplification was enriched in basal-like subtype samples and associated with poor clinical outcome. Surprisingly, based on the univariate Cox regression analysis, CD24 overexpression (HR = 1.62, P = 0.010) and copy number amplification (HR = 1.79, P = 0.022) was more relevant to OS than TP53 mutant, mutation counts, diagnosis age, and BRCA subtypes. And based on multivariate survival analysis, CD24 amplification remained the most significant and independent predictor for worse OS (HR = 1.88, P = 0.015).

Robust identification of target genes and outliers in triple-negative breast cancer data

Correct classification of breast cancer subtypes is of high importance as it directly affects the therapeutic options. We focus on triple-negative breast cancer which has the worst prognosis among breast cancer types. Using cutting edge methods from the field of robust statistics, we analyze Breast Invasive Carcinoma transcriptomic data publicly available from The Cancer Genome Atlas data portal. Our analysis identifies statistical outliers that may correspond to misdiagnosed patients. Furthermore, it is illustrated that classical statistical methods may fail to identify outliers due to their heavy influence, prompting the need for robust statistics. Using robust sparse logistic regression we obtain 36 relevant genes, of which ca. 60% have been previously reported as biologically relevant to triple-negative breast cancer, reinforcing the validity of the method. The remaining 14 genes identified are new potential biomarkers for triple-negative breast cancer. Out of these, JAM3, SFT2D2, and PAPSS1 were previously associated to breast tumors or other types of cancer. The relevance of these genes is confirmed by the new DetectDeviatingCells outlier detection technique. A comparison of gene networks on the selected genes showed significant differences between triple-negative breast cancer and non-triple-negative breast cancer data. The individual role of FOXA1 in triple-negative breast cancer and non-triple-negative breast cancer, and the strong FOXA1-AGR2 connection in triple-negative breast cancer stand out. The goal of our paper is to contribute to the breast cancer/triple-negative breast cancer understanding and management. At the same time it demonstrates that robust regression and outlier detection constitute key strategies to cope with high-dimensional clinical data such as omics data.

Roles of Grainyhead-like transcription factors in cancer

The mammalian homologs of the D. melanogaster Grainyhead gene, Grainyhead-like 1-3 (GRHL1, GRHL2 and GRHL3), are transcription factors implicated in wound healing, tubulogenesis and cancer. Their induced target genes encode diverse epithelial cell adhesion molecules, while mesenchymal genes involved in cell migration and invasion are repressed. Moreover, GRHL2 suppresses the oncogenic epithelial-mesencyhmal transition, thereby acting as a tumor suppressor. Mechanisms, some involving established cancer-related signaling/transcription factor pathways (for example, Wnt, TGF-β, mir200, ZEB1, OVOL2, p63 and p300) and translational implications of the Grainyhead proteins in cancer are discussed in this review article.

MicroRNA Signature of Lung Adenocarcinoma with EGFR Exon 19 Deletion

The findings of EGFR mutations and the development of targeted therapies have significantly improved the overall survival of lung cancer patients. Still, the prognosis remains poor, so we need to know more about the genetic alterations in lung cancer. MicroRNAs are dysregulated in lung cancer, and some of them can regulate EGFR. So it is very important to predict the candidate microRNAs that target mutated EGFR and to investigate the role of these candidate microRNAs in lung cancer. In this study, we investigated the difference of microRNAs expression between lung adenocarcinoma cell lines with EGFR exon 19 deletion (H1650 and PC9) and wild-type (H1299 and A549) using the Phalanx Human Whole Genome Microarray. Then the expression of individual microRNAs was validated by qRT-PCR assays. Moreover, we detected the microRNAs expression in plasma of lung adenocarcinoma patients with EGFR exon 19 deletion and wild-type. Lastly, we explored the function of the positive microRNA in EGFR tyrosine kinase inhibitors (EGFR-TKIs ) resistance using MTT and Annexin V-APC assays. The expression of 1,732 microRNAs was evaluated, and we found that microRNAs expression was different between these two groups. Hsa-miR-141-3p, hsa-miR-200c-3p, hsa-miR-203, hsa-miR-3182, hsa-miR-934 were up-regulated and hsa-miR-3196 was down-regulated in the EGFR exon 19 deletion group compared with wild-type group. The detection of circulating microRNAs showed that miR-3196 was down-regulated in lung adenocarcinoma patients with EGFR exon 19 deletion compared with wild-type. And then the MTT assay results showed that miR-3196 had no effect on the sensitivity of erlotinib. The results of apoptosis analysis showed that inhibition of miR-3196 and erlotinib induced more apoptosis in H1299 cells than erlotinib alone, and overexpressed miR-3196 and erlotinib induced less apoptosis in PC9 cells than erlotinib alone (P<0.05). It is suggested that microRNAs associate with EGFR exon 19 deletion and miR-3196 may be further explored as a potential predictor and targeted biomarker when it is difficult to get the tumors.

The activity of the Drosophila Vestigial protein is modified by Scalloped-dependent phosphorylation

The Drosophila vestigial gene is required for proliferation and differentiation of the adult wing and for differentiation of larval and adult muscle identity. Vestigial is part of a multi-protein transcription factor complex, which includes Scalloped, a TEAD-class DNA binding protein. Binding Scalloped is necessary for translocation of Vestigial into the nucleus. We show that Vestigial is extensively post-translationally modified and at least one of these modifications is required for proper function during development. We have shown that there is p38-dependent phosphorylation of Serine 215 in the carboxyl-terminal region of Vestigial. Phosphorylation of Serine 215 occurs in the nucleus and requires the presence of Scalloped. Comparison of a phosphomimetic and non-phosphorylatable mutant forms of Vestigial shows differences in the ability to rescue the wing and muscle phenotypes associated with a null vestigial allele.

Expression Profiling of Clinical Specimens Supports the Existence of Neural Progenitor-Like Stem Cells in Basal Breast Cancers

BACKGROUND

We previously characterized in salivary adenoid cystic carcinoma (ACC) a novel population of cancer stem cells (CSCs) marked by coexpression of 2 stemness genes, sex-determining region Y (SRY)-related HMG box-containing factor 10 (SOX10) and CD133. We also reported that in ACC and basal-like breast carcinoma (BBC), a triple-negative breast cancer subtype, expression of SOX10 similarly demarcates a highly conserved gene signature enriched with neural stem cell genes. On the basis of these findings, we hypothesized that BBC might be likewise driven by SOX10-positive (SOX10⁺)/CD133⁺ cells with neural stem cell properties.

MATERIALS AND METHODS

To validate our hypothesis on clinical data, we used a novel approach to meta-analysis that merges gene expression data from independent breast cancer studies and ranks genes according to statistical significance of their coexpression with the gene of interest. Genes that showed strong association with CD133/PROM1 as well as SOX10 were validated across different platforms and data sets and analyzed for enrichment with genes involved in neurogenesis.

RESULTS

We identified in clinical breast cancer data sets a highly conserved SOX10/PROM1 gene signature that contains neural stem cell markers common for Schwann cells, ACC, BBC, and melanoma. Identification of tripartite motif-containing 2 (TRIM2), TRIM29, MPZL2, potassium calcium-activated channel subfamily N member 4 (KCNN4), and V-set domain containing T cell activation inhibitor 1 (VTCN1)/B7 homolog 4 (B7H4) within this signature provides insight into molecular mechanisms of CSC maintenance.

CONCLUSION

Our results suggest that BBC is driven by SOX10⁺/CD133⁺ cells that express neural stem cell-specific markers and share molecular similarities with CSCs of neural crest origin. Our study provides clinically relevant information on possible drivers of these cells that might facilitate development of CSC-targeting therapies against this cancer distinguished with poor prognosis and resistance to conventional therapies.

Exploring the intrinsic differences among breast tumor subtypes defined using immunohistochemistry markers based on the decision tree

Exploring the intrinsic differences among breast cancer subtypes is of crucial importance for precise diagnosis and therapeutic decision-making in diseases of high heterogeneity. The subtypes defined with several layers of information are related but not consistent, especially using immunohistochemistry markers and gene expression profiling. Here, we explored the intrinsic differences among the subtypes defined by the estrogen receptor, progesterone receptor and human epidermal growth factor receptor 2 based on the decision tree. We identified 30 mRNAs and 7 miRNAs differentially expressed along the tree's branches. The final signature panel contained 30 mRNAs, whose performance was validated using two public datasets based on 3 well-known classifiers. The network and pathway analysis were explored for feature genes, from which key molecules including FOXQ1 and SFRP1 were revealed to be densely connected with other molecules and participate in the validated metabolic pathways. Our study uncovered the differences among the four IHC-defined breast tumor subtypes at the mRNA and miRNA levels, presented a novel signature for breast tumor subtyping, and identified several key molecules potentially driving the heterogeneity of such tumors. The results help us further understand breast tumor heterogeneity, which could be availed in clinics.

Peritoneal tumor spread in serous ovarian cancer-epithelial mesenchymal status and outcome

In this study we aimed to analyze the biological mechanisms underlying apparently different modes of peritoneal tumor spread in serous ovarian cancer: miliary (widespread, millet-like lesions) versus non-miliary (bigger, exophytically growing implants). Tumor tissues and ascites from 23 chemotherapy naive patients were analyzed by RNA-sequencing and flow cytometry. On the basis of differential gene expression between miliary and non-miliary, gene signatures were developed. A calculated tumor spread factor revealed a significant independent negative impact of miliary spread on overall survival (HR 3.77; CI₉₅ 1.14–12.39; p = 0.029) in an independent cohort of 165 serous ovarian cancer patients. Comparing previously published epithelial-mesenchymal transition (EMT) gene signatures, non-miliary spread correlated significantly with a reduced epithelial status. We conclude that serous ovarian cancer is a heterogeneous disease with distinct modes of peritoneal tumor spread, differing not only in clinical appearance, but also in molecular characteristics and outcome.. EMT, peritoneal inflammation status, and therapeutic options are discussed.

Construction of an initial microRNA regulation network in breast invasive carcinoma by bioinformatics analysis

INTRODUCTION

microRNAs (miRNA) are involved in many biological processes. They repress target gene expression and play a vital role in breast invasive carcinoma (BRCA). Although many miRNAs are identified to be aberrantly expressed in BRCA and deemed as tumor markers, only sporadic individual studies report their target genes and the pathways involved.

METHODS

miRNA and mRNA expression data were collected from the Cancer Genome Atlas (TCGA) pilot project. Aberrantly expressed miRNAs and mRNAs in BRCA were identified by comparing tumor samples with normal adjacent tissues. Differentially expressed miRNAs and mRNAs in different breast cancer subtypes were also analyzed. miRNA/target correlations were predicted by calculating the spearman correlation coefficients between miRNA and mRNA, and validated by luciferase assay.

RESULTS

31 up-regulated miRNAs, 37 down-regulated miRNAs, 1105 up-regulated mRNAs and 1222 down-regulated mRNAs were identified in BRCA; 125 miRNA/target correlations were predicted, 6 of them were validated. In addition, we also found 9 miRNAs and 143 mRNAs differently expressed between estrogen receptor positive and negative breast cancers, and 4 miRNAs and 46 mRNAs differently expressed between progesterone receptor positive and negative breast cancers. Twelve miRNA/target correlations determined the breast cancer subtypes.

CONCLUSION

We developed a new systematic analytic method for analyzing TCGA database, which took into account both miRNA and mRNA data to dissect the miRNA regulation network in BRCA.

Non-coding RNAs in Mammary Gland Development and Disease

Non-coding RNAs (ncRNAs) are untranslated RNA molecules that function to regulate the expression of numerous genes and associated biochemical pathways and cellular functions. NcRNAs include small interfering RNAs (siRNAs), microRNAs (miRNAs), PIWI-interacting RNAs (piRNAs), small nucleolar RNAs (snoRNAs) and long non-coding RNAs (lncRNAs). They participate in the regulation of all developmental processes and are frequently aberrantly expressed or functionally defective in disease. This Chapter will focus on the role of ncRNAs, in particular miRNAs and lncRNAs, in mammary gland development and disease.

Alcohol-dysregulated miR-30a and miR-934 in head and neck squamous cell carcinoma

BACKGROUND

Alcohol consumption is a well-established risk factor for head and neck squamous cell carcinoma (HNSCC); however, the molecular mechanisms by which alcohol promotes HNSCC pathogenesis and progression remain poorly understood. Our study sought to identify microRNAs that are dysregulated in alcohol-associated HNSCC and investigate their contribution to the malignant phenotype.

METHOD

Using RNA-sequencing data from 136 HNSCC patients, we compared the expression levels of 1,046 microRNAs between drinking and non-drinking cohorts. Dysregulated microRNAs were verified by qRT-PCR in normal oral keratinocytes treated with biologically relevant doses of ethanol and acetaldehyde. The most promising microRNA candidates were investigated for their effects on cellular proliferation and invasion, sensitivity to cisplatin, and expression of cancer stem cell genes. Finally, putative target genes were identified and evaluated in vitro to further establish roles for these miRNAs in alcohol-associated HNSCC.

RESULTS

From RNA-sequencing analysis we identified 8 miRNAs to be significantly upregulated in alcohol-associated HNSCCs. qRT-PCR experiments determined that among these candidates, miR-30a and miR-934 were the most highly upregulated in vitro by alcohol and acetaldehyde. Overexpression of miR-30a and miR-934 in normal and HNSCC cell lines produced up to a 2-fold increase in cellular proliferation, as well as induction of the anti-apoptotic gene BCL-2. Upon inhibition of these miRNAs, HNSCC cell lines exhibited increased sensitivity to cisplatin and reduced matrigel invasion. miRNA knockdown also indicated direct targeting of several tumor suppressor genes by miR-30a and miR-934.

CONCLUSIONS

Alcohol induces the dysregulation of miR-30a and miR-934, which may play crucial roles in HNSCC pathogenesis and progression. Future investigation of the alcohol-mediated pathways effecting these transformations will prove valuable for furthering the understanding and treatment of alcohol-associated HNSCC.

Nuclear TAZ expression associates with the triple-negative phenotype in breast cancer

The Hippo signaling pathway, a conserved regulator of organ size, has emerged as an important regulatory pathway in cancer. The final transducer effectors of this pathway in mammals are the oncoproteins TAZ and YAP1, which are transcriptional coactivators of target genes involved in cell proliferation and survival. TAZ has been previously reported to play a role in tumorigenesis in breast cancer, but detailed analyses of the different breast cancer phenotypes have not been conducted thus far. We analyzed TAZ expression by immunohistochemistry in a retrospective series of 640 invasive breast carcinomas, comprising estrogen/progesterone receptor-positive (ER+/PR+), HER2-positive, and triple-negative (TN) tumors. We found a strong association of TAZ nuclear expression with the TN phenotype (60.5% TAZ-positive, P<0.001), which was strengthened when stratified into the basal-like subtype (70.8% TAZ-positive, P<0.001). Moreover, 90% of metaplastic breast carcinomas with morphological epithelial-mesenchymal transition features were TAZ-positive. We also investigated whether amplification or differential DNA methylation of the TAZ-encoding locus could account for the observed enhanced TAZ protein expression in the TN/basal phenotype. Amplification of the TAZ locus was analyzed by fluorescence in situ hybridization in 30 TN tumors, and we found gene amplification in some cases (6.45%). DNA methylation analysis was performed using the Sequenom MassArray MALDI-TOF platform, and we observed similar low methylation levels both in TN (n=25) and ER+/PR+ (n=26) tumors. These results were further confirmed using a panel of breast cancer cell lines and using the TCGA dataset. Finally, patients with strong TAZ expression showed poorer clinical outcomes with respect to both recurrence and overall survival.