Search

The interplay between glycolysis and mitochondrial oxidative phosphorylation (OXPHOS) is central to maintain energy homeostasis. It remains to be determined whether there is a mechanism governing metabolic fluxes based on substrate availability in microenvironments. Here we show that menin is a key transcription factor regulating the expression of OXPHOS and glycolytic genes in cancer cells and primary tumors with poor prognosis. A group of menin-associated proteins (MAPs), including KMT2A, MED12, WAPL, and GATA3, is found to restrain menin&rsquo, s full function in this transcription regulation. shRNA knockdowns of menin and MAPs result in reduced ATP production with proportional alterations of cellular energy generated through glycolysis and OXPHOS. When shRNA knockdown cells are exposed to metabolic stress, the dual functionality can clearly be distinguished among these metabolic regulators. A MAP can negatively counteract the regulatory mode of menin for OXPHOS while the same protein positively influences glycolysis. A close-proximity interaction between menin and MAPs allows transcriptional regulation for metabolic adjustment. This coordinate regulation by menin and MAPs is necessary for cells to rapidly adapt to fluctuating microenvironments and to maintain essential metabolic functions.

Emerging evidence indicates that adipose stromal cells (ASC) are recruited to enhance cancer development. In this study, we examined the role these adipocyte progenitors play relating to intercellular communication in obesity-associated endometrial cancer. This is particularly relevant given that gap junctions have been implicated in tumor suppression. Examining the effects of ASCs on the transcriptome of endometrial epithelial cells (EEC) in an in vitro coculture system revealed transcriptional repression of GJA1 (encoding the gap junction protein Cx43) and other genes related to intercellular communication. This repression was recapitulated in an obesity mouse model of endometrial cancer. Furthermore, inhibition of plasminogen activator inhibitor 1 (PAI-1), which was the most abundant ASC adipokine, led to reversal of cellular distribution associated with the GJA1 repression profile, suggesting that PAI-1 may mediate actions of ASC on transcriptional regulation in EEC. In an endometrial cancer cohort (n = 141), DNA hypermethylation of GJA1 and related loci TJP2 and PRKCA was observed in primary endometrial endometrioid tumors and was associated with obesity. Pharmacologic reversal of DNA methylation enhanced gap-junction intercellular communication and cell–cell interactions in vitro. Restoring Cx43 expression in endometrial cancer cells reduced cellular migration; conversely, depletion of Cx43 increased cell migration in immortalized normal EEC. Our data suggest that persistent repression by ASC adipokines leads to promoter hypermethylation of GJA1 and related genes in the endometrium, triggering long-term silencing of these loci in endometrial tumors of obese patients. Significance: Studies reveal that adipose-derived stem cells in endometrial cancer pathogenesis influence epigenetic repression of gap junction loci, which suggests targeting of gap junction activity as a preventive strategy for obesity-associated endometrial cancer.

Endometrial cancer is the most common gynecological cancer in the United States. We wanted to identify epigenetic aberrations involving microRNAs (miRNAs), whose genes become hypermethylated in endometrial primary tumors. By integrating known miRNA sequences from the miRNA database (miRBase) with DNA methylation data from methyl-CpG-capture sequencing, we identified 111 differentially methylated regions (DMRs) associated with CpG islands (CGIs) and miRNAs. Among them, 22 DMRs related to 29 miRNAs and within 8 kb of CGIs were hypermethylated in endometrial tumors but not in normal endometrium. miR-137 was further validated in additional endometrial primary tumors. Hypermethylation of miR-137 was found in both endometrioid and serous endometrial cancer (P

Radiation therapy (XRT) is a standard treatment for prostate cancer (PCa). Although dose escalation increases local control, toxicity hampers further escalation. Broader improvement will be possible by the addition of adjuvant therapies, which can synergize with radiation and thus improve efficacy. We have identified a natural compound (Nexrutine, Nx) that inhibits the survival and growth of PCa cells in combination with radiation. Combination studies demonstrated strong interaction between Nx and radiation both in vitro in multiple PCa cell lines and in the Transgenic adenocarcinoma of mouse prostate (TRAMP) model. Nx potentiated growth inhibitory effects of IR by down regulating ribosomal protein S6K (RPS6KB1), CyclinD1, Chk1 and HIF-1 α and prolonging G2/M checkpoint block. RPS6KB1 is upregulated in prostate cancers and its expression is correlated with tumor grade. Knockdown of RPS6KB1 in PCa cells increased their sensitivity toward radiation-induced survival inhibition. Overall, we provide scientific evidence (i) in support of Nx as an adjuvant in PCa patients receiving XRT (ii) suggesting that RPS6KB1 is an important player in Nx-mediated combinatorial benefits and emphasizes that RPS6KB1 is a novel target for PCa treatment. These data underscore the need to test the agent in additional preclinical models to validate these observations.

Mucinous type of epithelial ovarian cancer (MuOC) is a unique subtype with a poor survival outcome in recurrent and advanced stages. The role of type-specific epigenomics and its clinical significance remains uncertain. We analyzed the methylomic profiles of 6 benign mucinous adenomas, 24 MuOCs, 103 serous type of epithelial ovarian cancers (SeOCs) and 337 nonepithelial ovarian cancers. MuOC and SeOC exhibited distinct DNA methylation profiles comprising 101 genes, 81 of which exhibited low methylation in MuOC and were associated with the response to glucocorticoid, ATP hydrolysis-coupled proton transport, proteolysis involved in the cellular protein catabolic process and ion transmembrane transport. Hierarchical clustering analysis showed that the profiles of MuOC were similar to colorectal adenocarcinoma and stomach adenocarcinoma. Genetic interaction network analysis of differentially methylated genes in MuOC showed a dominant network module is the proteasome subunit beta (PSMB) family. Combined functional module and methylation analysis identified PSMB8 as a candidate marker for MuOC. Immunohistochemical staining of PSMB8 used to validate in 94 samples of ovarian tumors (mucinous adenoma, MuOC or SeOC) and 62 samples of gastrointestinal cancer. PSMB8 was commonly expressed in MuOC and gastrointestinal cancer samples, predominantly as strong cytoplasmic and occasionally weak nuclei staining, but was not expressed in SeOC samples. Carfilzomib, a second-generation proteasome inhibitor, suppressed MuOC cell growth in vitro. This study unveiled a mucinous-type-specific methylation profile and suggests the potential use of a proteasome inhibitor to treat MuOC.

Increasing evidence suggests the presence of minor cell subpopulations in prostate cancer that are androgen independent and poised for selection as dominant clones after androgen deprivation therapy. In this study, we investigated this phenomenon by stratifying cell subpopulations based on transcriptome profiling of 144 single LNCaP prostate cancer cells treated or untreated with androgen after cell-cycle synchronization. Model-based clustering of 397 differentially expressed genes identified eight potential subpopulations of LNCaP cells, revealing a previously unappreciable level of cellular heterogeneity to androgen stimulation. One subpopulation displayed stem-like features with a slower cell doubling rate, increased sphere formation capability, and resistance to G2–M arrest induced by a mitosis inhibitor. Advanced growth of this subpopulation was associated with enhanced expression of 10 cell-cycle–related genes (CCNB2, DLGAP5, CENPF, CENPE, MKI67, PTTG1, CDC20, PLK1, HMMR, and CCNB1) and decreased dependence upon androgen receptor signaling. In silico analysis of RNA-seq data from The Cancer Genome Atlas further demonstrated that concordant upregulation of these genes was linked to recurrent prostate cancers. Analysis of receiver operating characteristic curves implicates aberrant expression of these genes and could be useful for early identification of tumors that subsequently develop biochemical recurrence. Moreover, this single-cell approach provides a better understanding of how prostate cancer cells respond heterogeneously to androgen deprivation therapies and reveals characteristics of subpopulations resistant to this treatment. Significance: Illustrating the challenge in treating cancers with targeted drugs, which by selecting for drug resistance can drive metastatic progression, this study characterized the plasticity and heterogeneity of prostate cancer cells with regard to androgen dependence, defining the character or minor subpopulations of androgen-independent cells that are poised for clonal selection after androgen-deprivation therapy. Cancer Res; 78(4); 853–64. ©2017 AACR.

Negative valine (V) to phenylalanine (F) switch at the Janus kinase (JAK2) 617 codon (V617F) is the dominant driver mutation in patients with myeloproliferative neoplasms (MPNs). JAK2V617F was proved to be sufficient for cell transformation; however, independent mutations might influence the following epigenomic modifications. To assess the JAK2V617F-induced downstream epigenomic changes without interferences, we profiled the epigenomic changes in ectopically expressed JAK2V617F in Ba/F3 cells. Antibodies against phosphorylated signal transducer and activator of transcription 3 (pSTAT3) and enhancer of zeste homolog 2 (EZH2) were used for chromatin-immunoprecipitation sequencing (ChIP-seq) to detect the downstream epigenomic targets in the JAK2-STAT3 signaling pathway. To confirm the JAK2V617F-induced epigenetic changes in vivo, DNA methylation changes in the target loci in patients with MPNs were detected through methylation-specific polymerase chain reaction and were clustered against the changes within controls. We found that ectopically expressed JAK2V617F in Ba/F3 cells reduced the binding specificity; it was associated with cis-regulatory elements and recognized DNA motifs in both pSTAT3-downstream and EZH2-associated targets. Overlapping target loci between the control and JAK2V617F were3% and 0.4%, respectively, as identified through pSTAT3 and EZH2 ChIP-seq. Furthermore, the methylation changes in the direct target loci (FOXH1, HOXC9, and SRF) were clustered independently from the control locus (L1TD1) and other mutation genes (HMGA2 and Lin28A) in the analyzed MPN samples. Therefore, JAK2V617F influences target binding in both pSTAT3 and EZH2. Without mutations in epigenetic regulators, JAK2V617F can induce downstream epigenomic modifications. Thus, epigenetic changes in JAK2 downstream targets might be trackable in vivo.

Background A genomic signal track is a set of genomic intervals associated with values of various types, such as measurements from high-throughput experiments. Analysis of signal tracks requires complex computational methods, which often make the analysts focus too much on the detailed computational steps rather than on their biological questions. Results Here we propose Signal Track Query Language (STQL) for simple analysis of signal tracks. It is a Structured Query Language (SQL)-like declarative language, which means one only specifies what computations need to be done but not how these computations are to be carried out. STQL provides a rich set of constructs for manipulating genomic intervals and their values. To run STQL queries, we have developed the Signal Track Analytical Research Tool (START, http://yiplab.cse.cuhk.edu.hk/start/), a system that includes a Web-based user interface and a back-end execution system. The user interface helps users select data from our database of around 10,000 commonly-used public signal tracks, manage their own tracks, and construct, store and share STQL queries. The back-end system automatically translates STQL queries into optimized low-level programs and runs them on a computer cluster in parallel. We use STQL to perform 14 representative analytical tasks. By repeating these analyses using bedtools, Galaxy and custom Python scripts, we show that the STQL solution is usually the simplest, and the parallel execution achieves significant speed-up with large data files. Finally, we describe how a biologist with minimal formal training in computer programming self-learned STQL to analyze DNA methylation data we produced from 60 pairs of hepatocellular carcinoma (HCC) samples. Conclusions Overall, STQL and START provide a generic way for analyzing a large number of genomic signal tracks in parallel easily. Electronic supplementary material The online version of this article (doi:10.1186/s12864-017-4071-1) contains supplementary material, which is available to authorized users.

TEA domain (TEAD) transcription factors are key components of the Hippo-YAP1 signaling pathway, but their functional role and regulatory mechanisms remain unclear. This study aims to comprehensively explore the expression pattern and functional role of TEAD family in gastric carcinogenesis and investigate its regulation by microRNAs (miRNAs). The mRNA and protein expression of TEAD family were examined by quantitative reverse transcription-PCR (qRT-PCR) and western blot. Their functional roles were determined by in vitro and in vivo studies. The clinicopathological association of TEAD4 in gastric cancer (GC) was studied using immunohistochemistry on tissue microarray. The prediction of miRNAs, which potentially target TEAD1/4, was performed by TargetScan and miRDB. The regulation of TEAD1/4 by miRNAs was confirmed by qRT-PCR, western blot and luciferase assays. TEAD1/4 were overexpressed in GC cell lines and primary GC tissues. Knockdown of TEAD1/4 induced a significant anticancer effect in vitro and in vivo. TEAD1 was confirmed to be a direct target of miR-377-3p and miR-4269, while TEAD4 was negatively regulated by miR-1343-3p and miR-4269. Among them, miR-4269 was the most effective inhibitor of TEAD1/4. Ectopic expression of these miRNAs substantiated their tumor-suppressive effects. In primary GC tumors, downregulation of miR-4269 was associated with poor disease-specific survival and showed a negative correlation with TEAD4. TEAD1 and TEAD4 are oncogenic factors, whose aberrant activation are, in part, mediated by the silence of miR-377-3p, miR-1343-3p and miR-4269. For the first time, the nuclear accumulated TEAD4 and downregulated miR-4269 are proposed to serve as novel prognostic biomarkers in GC.

Despite numerous studies done on understanding the role of DNA methylation, limited work has focused on systems integration of cell type-specific interplay between DNA methylation and gene transcription. Through a genome-wide analysis of DNA methylation across 19 cell types with T-47D as reference, we identified 106,252 cell type-specific differentially-methylated CpGs categorized into 7,537 differentially (46.6% hyper- and 53.4% hypo-) methylated regions. We found 44% promoter regions and 75% CpG islands were T-47D cell type-specific methylated. Pyrosequencing experiments validated the cell type-specific methylation across three benchmark cell lines. Interestingly, these DMRs overlapped with 1,145 known tumor suppressor genes. We then developed a Bayesian Gaussian Regression model to measure the relationship among DNA methylation, genomic segment distribution, differential gene expression and tumor suppressor gene status. The model uncovered that 3′UTR methylation has much less impact on transcriptional activity than other regions. Integration of DNA methylation and 82 transcription factor binding information across the 19 cell types suggested diverse interplay patterns between the two regulators. Our integrative analysis reveals cell type-specific and genomic region-dependent regulatory patterns and provides a perspective for integrating hundreds of various omics-seq data together.

Cytometry by time-of-flight (CyTOF) simultaneously measures multiple cellular proteins at the single-cell level and is used to assess intertumor and intratumor heterogeneity. This approach may be used to investigate the variability of individual tumor responses to treatments. Herein, we stratified lung tumor subpopulations based on AXL signaling as a potential targeting strategy. Integrative transcriptome analyses were used to investigate how TP-0903, an AXL kinase inhibitor, influences redundant oncogenic pathways in metastatic lung cancer cells. CyTOF profiling revealed that AXL inhibition suppressed SMAD4/TGFβ signaling and induced JAK1–STAT3 signaling to compensate for the loss of AXL. Interestingly, high JAK1–STAT3 was associated with increased levels of AXL in treatment-naïve tumors. Tumors with high AXL, TGFβ, and JAK1 signaling concomitantly displayed CD133-mediated cancer stemness and hybrid epithelial-to-mesenchymal transition features in advanced-stage patients, suggesting greater potential for distant dissemination. Diffusion pseudotime analysis revealed cell-fate trajectories among four different categories that were linked to clinicopathologic features for each patient. Patient-derived organoids (PDO) obtained from tumors with high AXL and JAK1 were sensitive to TP-0903 and ruxolitinib (JAK inhibitor) treatments, supporting the CyTOF findings. This study shows that single-cell proteomic profiling of treatment-naïve lung tumors, coupled with ex vivo testing of PDOs, identifies continuous AXL, TGFβ, and JAK1–STAT3 signal activation in select tumors that may be targeted by combined AXL–JAK1 inhibition. Significance: Single-cell proteomic profiling of clinical samples may facilitate the optimal selection of novel drug targets, interpretation of early-phase clinical trial data, and development of predictive biomarkers valuable for patient stratification.

NAD[P]H:quinone oxidoreductase 1 (NQO1) regulates cell fate decisions in response to stress. Oxidative stress supports cancer maintenance and progression. Previously we showed that knockdown of NQO1 (NQO1low) prostate cancer cells upregulate pro-inflammatory cytokines and survival under hormone-deprived conditions. Here, we tested the ability of NQO1low cells to form tumors. We found NQO1low cells form aggressive tumors compared with NQO1high cells. Biopsy specimens and circulating tumor cells showed biochemical recurrent prostate cancer was associated with low NQO1. NQO1 silencing was sufficient to induce SMAD-mediated TGFβ signaling and mesenchymal markers. TGFβ treatment decreased NQO1 levels and induced molecular changes similar to NQO1 knockdown cells. Functionally, NQO1 depletion increased migration and sensitivity to oxidative stress. Collectively, this work reveals a possible new gatekeeper role for NQO1 in counteracting cellular plasticity in prostate cancer cells. Further, combining NQO1 with TGFβ signaling molecules may serve as a better signature to predict biochemical recurrence., Thapa et al find that depletion of the antioxidant enzyme NAD[P]H:Quinone Oxidoreductase 1 (NQO1) accelerates prostate tumorigenesis and induces the epithelial-to-mesenchymal transition by activating TGFβ signaling. They also find that low NQO1 is associated with mesenchymal signature and biochemical recurrence in clinical samples.

Advanced prostate cancer is a very heterogeneous disease reflecting in diverse regulations of oncogenic signaling pathways. Aberrant spatial dynamics of epidermal growth factor receptor (EGFR) promote their dimerization and clustering, leading to constitutive activation in oncogenesis. The EphB2 and Src signaling pathways are associated with the reorganization of the cytoskeleton leading to malignancy, but their roles in regulating EGFR dynamics and activation are scarcely reported. Using single-particle tracking techniques, we found that highly phosphorylated EGFR in the advanced prostate cancer cell line, PC3, was associated with higher EGFR diffusivity, as compared with LNCaP and less aggressive DU145. The increased EGFR activation and biophysical dynamics were consistent with high proliferation, migration, and invasion. After performing single-cell RNA-seq on prostate cancer cell lines and circulating tumor cells from patients, we identified that upregulated gene expression in the EphB2 and Src pathways are associated with advanced malignancy. After dasatinib treatment or siRNA knockdowns of EphB2 or Src, the PC3 cells exhibited significantly lower EGFR dynamics, cell motility, and invasion. Partial inhibitory effects were also found in DU145 cells. The upregulation of parts of the EphB2 and Src pathways also predicts poor prognosis in the prostate cancer patient cohort of The Cancer Genome Atlas. Our results provide evidence that overexpression of the EphB2 and Src signaling pathways regulate EGFR dynamics and cellular aggressiveness in some advanced prostate cancer cells.

The role of cholesterol content on monocyte biomechanics remains understudied despite the well-established link between cholesterol and monocytes/macrophages in atherosclerosis, and the effect on other cell types. In this work, we have investigated the effect of cholesterol on monocyte deformability and the underlying molecular mechanisms. We altered the baseline cholesterol in human monocytic cell line THP-1, and investigated the changes in monocyte deformability using a custom microfluidic platform and atomic force microscopy. We observed that the cholesterol depletion lowered deformability while enrichment increased deformability compared to untreated cells. As a consequence of altered deformability, cholesterol depleted cells spread more on collagen-coated surfaces with elongated morphology, whereas cholesterol enriched cells had a more rounded morphology. We observed that the decreased deformability in cholesterol depleted cells, despite an increase in the fluidity of the membrane, is due to an increase in phosphorylation of Protein Kinase C (PKC), which translates to a higher degree of actin polymerization. Together, our results highlight the importance of biophysical regulation of monocyte response to cholesterol levels.

Purpose: Endometrial cancer is a common gynecologic cancer whose incidence is increasing annually worldwide. Current methods to detect endometrial cancer are unreliable and biomarkers are unsatisfactory for screening. Cervical scrapings were reported as a potential source of material for molecular testing. DNA methylation is a promising cancer biomarker, but limited use for detecting endometrial cancer. Experimental Design: We analyzed two methylomics databases of endometrioid-type endometrial cancer. Using nonnegative matrix factorization algorithm clustered the methylation pattern and reduced the candidate genes. We verified in pools DNA from endometrial cancer tissues and cervical scrapings, and validated in 146 cervical scrapings from patients with endometrioid-type endometrial cancer (n = 50), uterine myoma (n = 40), and healthy controls (n = 56) using quantitative methylation–specific PCR (QMSP). The logistic regression was used to evaluate the performance of methylation signal and gene combination. Results: We filtered out 180 methylated genes, which constituted four consensus clusters. Serial testing of tissues and cervical scrapings detected 14 genes that are hypermethylated in endometrial cancer. Three genes, BHLHE22, CDO1, and CELF4, had the best performance. Individual genes were sensitivity of 83.7%–96.0% and specificity of 78.7%–96.0%. A panel comprising any two of the three hypermethylated genes reached a sensitivity of 91.8%, specificity of 95.5%, and odds ratio of 236.3 (95% confidence interval, 56.4–989.6). These markers were also applied to cervical scrapings of type II endometrial cancer patients, and detected in 13 of 14 patients. Conclusions: This study demonstrates the potential use of methylated BHLHE22/CDO1/CELF4 panel for endometrial cancer screening of cervical scrapings. Clin Cancer Res; 23(1); 263–72. ©2016 AACR.

The compaction of nucleosomal structures creates a barrier for DNA-binding transcription factors (TFs) to access their cognate cis-regulatory elements. Pioneer factors (PFs) such as FOXA1 are able to directly access these cis-targets within compact chromatin. However, how these PFs interplay with nucleosomes remains to be elucidated, and is critical for us to understand the underlying mechanism of gene regulation. Here, we have conducted a computational analysis on a strand-specific paired-end ChIP-exo (termed as ChIP-ePENS) data of FOXA1 in LNCaP cells by our novel algorithm ePEST. We find that FOXA1 chromatin binding occurs via four distinct border modes (or footprint boundary patterns), with a preferential footprint boundary patterns relative to FOXA1 motif orientation. In addition, from this analysis three fundamental nucleotide positions (oG, oS and oH) emerged as major determinants for blocking exo-digestion and forming these four distinct border modes. By integrating histone MNase-seq data, we found an astonishingly consistent, 'well-positioned' configuration occurs between FOXA1 motifs and dyads of nucleosomes genome-wide. We further performed ChIP-seq of eight chromatin remodelers and found an increased occupancy of these remodelers on FOXA1 motifs for all four border modes (or footprint boundary patterns), indicating the full occupancy of FOXA1 complex on the three blocking sites (oG, oS and oH) likely produces an active regulatory status with well-positioned phasing for protein binding events. Together, our results suggest a positional-nucleosome-oriented accessing model for PFs seeking target motifs, in which FOXA1 can examine each underlying DNA nucleotide and is able to sense all potential motifs regardless of whether they face inward or outward from histone octamers along the DNA helix axis.

The transcription factor nuclear factor kappa B (NF-κB) and the long non-coding RNA (lncRNA) HOTAIR (HOX transcript antisense RNA) play diverse functional roles in cancer. In this study, we show that upregulation of HOTAIR induced platinum resistance in ovarian cancer, and increased HOTAIR levels were observed in recurrent platinum-resistant ovarian tumors vs. primary ovarian tumors. To investigate the role of HOTAIR during DNA damage induced by platinum, we monitored double-strand breaks and show that HOTAIR expression results in sustained activation of DNA damage response after platinum treatment. We demonstrate that ectopic expression of HOTAIR induces NF-κB activation during DNA damage response and MMP-9 and IL-6 expression, both key NF-κB target genes. We show that HOTAIR regulates activation of NF-κB by decreasing Iκ-Bα (NF-κB inhibitor) and establish that by inducing prolonged NF-κB activation and expression of NF-κB target genes during DNA damage, HOTAIR plays a critical role in cellular senescence and platinum sensitivity. Our findings suggest that a NF-κB-HOTAIR axis drives a positive-feedback loop cascade during DNA damage response and contributes to cellular senescence and chemotherapy resistance in ovarian and other cancers.

SUMMARY While plasminogen activator inhibitor-1 (PAI-1) is known to potentiate cellular migration via proteolytic regulation, this adipokine is implicated as an oncogenic ligand in the tumor microenvironment. To understand the underlying paracrine mechanism, here, we conduct transcriptomic analysis of 1,898 endometrial epithelial cells (EECs) exposed and unexposed to PAI-1-secreting adipose stromal cells. The PAI-1-dependent action deregulates crosstalk among tumor-promoting and tumor-repressing pathways, including transforming growth factor β (TGF-β). When PAI-1 is tethered to lipoprotein receptor-related protein 1 (LRP1), the internalized signaling causes downregulation of SMAD4 at the transcriptional and post-translational levels that attenuates TGF-β-related transcription programs. Repression of genes encoding the junction and adhesion complex preferentially occurs in SMAD4-underexpressed EECs of persons with obesity. The findings highlight a role of PAI-1 signaling that renders ineffective intercellular communication for the development of adiposity-associated endometrial cancer., Graphical Abstract, In Brief Lin et al. demonstrate that PAI-1 secreted by adipose stromal cells interacts with LRP1 to repress TGF-β/SMAD4-regulated genes linked to cellular junction and adhesion complexes. This action disrupts cell-cell communication and facilitates the development of obesity-driven endometrial cancer.

Hepatocellular carcinomas (HCC) exhibit distinct promoter hypermethylation patterns, but the epigenetic regulation and function of transcriptional enhancers remain unclear. Here, our affinity- and bisulfite-based whole-genome sequencing analyses reveal global enhancer hypomethylation in human HCCs. Integrative epigenomic characterization further pinpoints a recurrent hypomethylated enhancer of CCAAT/enhancer-binding protein-beta (C/EBPβ) which correlates with C/EBPβ over-expression and poorer prognosis of patients. Demethylation of C/EBPβ enhancer reactivates a self-reinforcing enhancer-target loop via direct transcriptional up-regulation of enhancer RNA. Conversely, deletion of this enhancer via CRISPR/Cas9 reduces C/EBPβ expression and its genome-wide co-occupancy with BRD4 at H3K27ac-marked enhancers and super-enhancers, leading to drastic suppression of driver oncogenes and HCC tumorigenicity. Hepatitis B X protein transgenic mouse model of HCC recapitulates this paradigm, as C/ebpβ enhancer hypomethylation associates with oncogenic activation in early tumorigenesis. These results support a causal link between aberrant enhancer hypomethylation and C/EBPβ over-expression, thereby contributing to hepatocarcinogenesis through global transcriptional reprogramming., There are distinct hypermethylation patterns in gene promoters in hepatocellular carcinomas (HCCs). Here, the authors show that the enhancer of C/EBPβ is recurrently hypomethylated in human HCCs, recapitulating this in a transgenic murine model and linking aberrant enhancer hypomethylation to hepatocarcinogenesis.

The importance of RNA methylation in biological processes is an emerging focus of investigation. We report that altering m6A levels by silencing either N6-adenosine methyltransferase METTL14 (methyltransferase-like 14) or demethylase ALKBH5 (ALKB homolog 5) inhibits cancer growth and invasion. METTL14/ALKBH5 mediate their protumorigenic function by regulating m6A levels of key epithelial-mesenchymal transition and angiogenesis-associated transcripts, including transforming growth factor–β signaling pathway genes. Using MeRIP-seq (methylated RNA immunoprecipitation sequencing) analysis and functional studies, we find that these target genes are particularly sensitive to changes in m6A modifications, as altered m6A status leads to aberrant expression of these genes, resulting in inappropriate cell cycle progression and evasion of apoptosis. Our results reveal that METTL14 and ALKBH5 determine the m6A status of target genes by controlling each other’s expression and by inhibiting m6A reader YTHDF3 (YTH N6-methyladenosine RNA binding protein 3), which blocks RNA demethylase activity. Furthermore, we show that ALKBH5/METTL14 constitute a positive feedback loop with RNA stability factor HuR to regulate the stability of target transcripts. We discover that hypoxia alters the level/activity of writers, erasers, and readers, leading to decreased m6A and consequently increased expression of target transcripts in cancer cells. This study unveils a previously undefined role for m6A in cancer and shows that the collaboration among writers-erasers-readers sets up the m6A threshold to ensure the stability of progrowth/proliferation-specific genes, and protumorigenic stimulus, such as hypoxia, perturbs that m6A threshold, leading to uncontrolled expression/activity of those genes, resulting in tumor growth, angiogenesis, and progression.

// Izhar Batth 1, 8 , Huiyoung Yun 2 , Suleman Hussain 2 , Peng Meng 1, 7 , Pawel Osmulski 3 , Tim Hui-Ming Huang 3, 5 , Roble Bedolla 1 , Amanda Profit 4 , Robert Reddick 4 , Addanki Kumar 1, 2, 3, 5, 6 1 Department of Urology, The University of Texas Health Science Center, San Antonio, TX, USA 2 Department of Pharmacology, The University of Texas Health Science Center, San Antonio, TX, USA 3 Department of Molecular Medicine, The University of Texas Health Science Center, San Antonio, TX, USA 4 Department of Pathology, The University of Texas Health Science Center, San Antonio, TX, USA 5 Cancer Therapy and Research Center, The University of Texas Health Science Center, San Antonio, TX, USA 6 The University of Texas Health Science Center at San Antonio and South Texas Veterans Health Care System, San Antonio, TX, USA 7 Current address: Life Sciences Division, Lawrence Berkley National Laboratory, Berkley, CA, USA 8 Current address: Department of Pediatrics, The University of Texas MD Anderson Cancer Center, Houston, TX, USA Correspondence to: Addanki Kumar, e-mail: kumara3@uthscsa.edu Keywords: castrate resistant prostate cancer, apoptosis, FLIP, RON, MST1R Received: November 18, 2015 Accepted: January 29, 2016 Published: February 9, 2016 ABSTRACT Castrate-resistant prostate cancer (CRPC) is the fatal form of prostate cancer. Although reactivation of androgen receptor (AR) occurs following androgen deprivation, the precise mechanism involved is unclear. Here we show that the receptor tyrosine kinase, RON alters mechanical properties of cells to influence epithelial to mesenchymal transition and functions as a transcription factor to differentially regulate AR signaling. RON inhibits AR activation and subset of AR-regulated transcripts in androgen responsive LNCaP cells. However in C4-2B, a castrate-resistant sub-line of LNCaP and AR-negative androgen independent DU145 cells, RON activates subset of AR-regulated transcripts. Expression of AR in PC-3 cells leads to activation of RON under androgen deprivation but not under androgen proficient conditions implicating a role for RON in androgen independence. Consistently, RON expression is significantly elevated in castrate resistant prostate tumors. Taken together our results suggest that RON activation could aid in promoting androgen independence and that inhibition of RON in combination with AR antagonist(s) merits serious consideration as a therapeutic option during hormone deprivation therapy.

Methylation changes within tumor suppressor (TS) genes or polycomb group target (PcG) genes alter cell fates. Chromatin associated with PcG targets is bivalent in stem cells, while TS genes are not normally bivalent. PcG target methylation changes have been identified in tumor stem cells, and abnormal methylation is found in TS genes in cancers. If the epigenetic states of genes influence DNA methylation, then methylation of PcG targets and TS genes may evolve differently during cancer development. More importantly, methylation changes may be part of a sequence in tumorigenesis. Chromatin and methylation states of 4 PcG targets and 2 TS genes were determined in colon cancer cells. The methylation states were also detected in 100 pairs of colon cancer samples. Principle component analysis (PCA) was used to reveal whether TS methylation or PcG methylation was the main methylation change associated with colon cancers. Chromatin and methylation states differ in colon cancer cell lines. The methylation states within PcG targets clustered independently from the methylation states in TS genes, a finding we previously reported in liver cancers. PCA in colon cancers revealed the strongest association with methylation changes in 2 TS genes, HIC1 and RassF1A. Loss of HIC1 methylation correlated with decreased tumor migration. PcG and TS methylation states cluster independently from each other. The deduced principle component correlated better with TS methylation than PcG methylation in colon cancer. Abnormal methylation changes may represent a sequential biomarker profile to identify developing colon cancer.

Resistance to cancer treatment can be driven by epigenetic reprogramming of specific transcriptomes in favor of the refractory phenotypes. Here we discover that tamoxifen resistance in breast cancer is driven by a regulatory axis consisting of a master transcription factor, its cofactor, and an epigenetic regulator. The oncogenic histone methyltransferase EZH2 conferred tamoxifen resistance by silencing the expression of the estrogen receptor a (ERa) cofactor GREB1. In clinical specimens, induction of DNA methylation of a particular CpG-enriched region at the GREB1 promoter negatively correlated with GREB1 levels and cell sensitivity to endocrine agents. GREB1 also ensured proper cellular reactions to different ligands by recruiting distinct sets of ERa cofactors to cis-regulatory elements, which explains the contradictory biological effects of GREB1 on breast cancer cell growth in response to estrogen or antiestrogen. In refractory cells, EZH2-dependent repression of GREB1 triggered chromatin reallocation of ERa coregulators, converting the antiestrogen into an agonist. In clinical specimens from patients receiving adjuvant tamoxifen treatment, expression levels of EZH2 and GREB1 were correlated negatively, and taken together better predicted patient responses to endocrine therapy. Overall, our work suggests a new strategy to overcome endocrine resistance in metastatic breast cancer by targeting a particular epigenetic program. Significance: This study suggests a new strategy to overcome endocrine resistance in metastatic breast cancer by targeting a particular epigenetic program defined within. Fil: Wu, Yanming. University of Texas at San Antonio; Estados Unidos Fil: Zhang, Zhao. University of Texas at San Antonio; Estados Unidos Fil: Cenciarini, Mauro Ezequiel. Consejo Nacional de Investigaciones Científicas y Técnicas. Instituto de Biología y Medicina Experimental. Fundación de Instituto de Biología y Medicina Experimental. Instituto de Biología y Medicina Experimental; Argentina Fil: Proietti Anastasi, Cecilia Jazmín. Consejo Nacional de Investigaciones Científicas y Técnicas. Instituto de Biología y Medicina Experimental. Fundación de Instituto de Biología y Medicina Experimental. Instituto de Biología y Medicina Experimental; Argentina Fil: Amasino, Matías Federico. Consejo Nacional de Investigaciones Científicas y Técnicas. Instituto de Biología y Medicina Experimental. Fundación de Instituto de Biología y Medicina Experimental. Instituto de Biología y Medicina Experimental; Argentina Fil: Hong, Tao. University of Texas at San Antonio; Estados Unidos. Central South University; China Fil: Yang, Mei. University of Texas at San Antonio; Estados Unidos Fil: Liao, Yiji. University of Texas at San Antonio; Estados Unidos Fil: Chiang, Huai-Chin. University of Texas at San Antonio; Estados Unidos. University Of Texas Health Science Center At San Antonio (ut Health San Antonio) ; University Of Texas At San Antonio; Fil: Kaklamani, Virginia G.. University Of Texas Health Science Center At San Antonio (ut Health San Antonio) ; University Of Texas At San Antonio; Fil: Jeselsohn, Rinath. Dana-farber Cancer Institute; Estados Unidos Fil: Vadlamudi, Ratna K.. University Of Texas Health Science Center At San Antonio (ut Health San Antonio) ; University Of Texas At San Antonio; Fil: Huang, Tim Hui Ming. University Of Texas Health Science Center At San Antonio (ut Health San Antonio) ; University Of Texas At San Antonio; Fil: Li, Rong. University Of Texas Health Science Center At San Antonio (ut Health San Antonio) ; University Of Texas At San Antonio; Fil: De Angelis, Carmine. Baylor College Of Medicine; Estados Unidos Fil: Fu, Xiaoyong. Baylor College Of Medicine; Estados Unidos Fil: Elizalde, Patricia Virginia. Consejo Nacional de Investigaciones Científicas y Técnicas. Instituto de Biología y Medicina Experimental. Fundación de Instituto de Biología y Medicina Experimental. Instituto de Biología y Medicina Experimental; Argentina Fil: Schiff, Rachel. Baylor College Of Medicine; Estados Unidos Fil: Brown, Myles. Dana farber Cancer Institute; Estados Unidos Fil: Xu, Kexin. University Of Texas Health Science Center At San Antonio (ut Health San Antonio) ; University Of Texas At San Antonio

Cholesterol is an important risk factor of atherosclerosis, due to its active uptake by monocytes/macrophages. Monocyte recruitment from flowing blood to atherosclerotic foci is the key first step in the development of atherosclerosis. Cholesterol content alters cell membrane stiffness, and lateral lipid and protein diffusion. We hypothesized that cholesterol content will modulate the recruitment of monocytes to inflamed endothelial surface by altering the dynamics of adhesion receptors. We depleted or enriched the cellular cholesterol levels using methyl-β-cyclodextran in freshly isolated human monocytes. We investigated the effect of these changes on the mechanics of monocyte rolling on E-selectin surfaces at 1 dyn/cm2 in microchannels. Using imaging flow cytometry and atomic force microscopy, we characterized the distribution of lipid rafts and the E-selectin counterreceptor CD44 on the monocyte surface. We observed that lower levels of cholesterol resulted in the uniform, CD44-mediated rolling of monocytes on the E-selectin-coated surfaces. We also observed that cells depleted of cholesterol had higher membrane fluidity, and more uniform distribution of CD44 counterreceptor, which resulted in smooth motion of the cells compared to cells enriched with cholesterol. This work demonstrates that cholesterol can modulate monocyte adhesion by regulating the receptor mobility, and our results provide insights into the biophysical regulation of inflammation for the better understanding of diseases like atherosclerosis and hypercholesterolemia.

Triple-negative breast cancer (TNBC) cells are modulated in reaction to tumor-infiltrating lymphocytes. However, their specific responses to this immune pressure are unknown. In order to address this question, we first used mRNA sequencing to compare the immunophenotype of the TNBC cell line MDA-MB-231 and the luminal breast cancer cell line MCF7 after both were cocultured with activated human T cells. Despite similarities in the cytokine-induced immune signatures of the two cell lines, MDA-MD-231 cells were able to transcribe more IDO1 than MCF7 cells. The two cell lines had similar upstream JAK/STAT1 signaling and IDO1 mRNA stability. However, using a series of breast cancer cell lines, IFNγ stimulated IDO1 protein expression and enzymatic activity only in ER−, not ER+, cell lines. Treatment with 5-aza-deoxycytidine reversed the suppression of IDO1 expression in MCF7 cells, suggesting that DNA methylation was potentially involved in IDO1 induction. By analyzing several breast cancer datasets, we discovered subtype-specific mRNA and promoter methylation differences in IDO1, with TNBC/basal subtypes exhibiting lower methylation/higher expression and ER+/luminal subtypes exhibiting higher methylation/lower expression. We confirmed this trend of IDO1 methylation by bisulfite pyrosequencing breast cancer cell lines and an independent cohort of primary breast tumors. Taken together, these findings suggest that IDO1 promoter methylation regulates anti-immune responses in breast cancer subtypes and could be used as a predictive biomarker for IDO1 inhibitor–based immunotherapy. Cancer Immunol Res; 5(4); 330–44. ©2017 AACR.

The constitutively active androgen receptor (AR) splice variant 7 (AR-V7) plays an important role in the progression of castration-resistant prostate cancer (CRPC). Although biomarker studies established the role of AR-V7 in resistance to AR-targeting therapies, how AR-V7 mediates genomic functions in CRPC remains largely unknown. Using a ChIP-exo approach, we show AR-V7 binds to distinct genomic regions and recognizes a full-length androgen-responsive element in CRPC cells and patient tissues. Remarkably, we find dramatic differences in AR-V7 cistromes across diverse CRPC cells and patient tissues, regulating different target gene sets involved in CRPC progression. Surprisingly, we discover that HoxB13 is universally required for and colocalizes with AR-V7 binding to open chromatin across CRPC genomes. HoxB13 pioneers AR-V7 binding through direct physical interaction, and collaborates with AR-V7 to up-regulate target oncogenes. Transcriptional coregulation by HoxB13 and AR-V7 was further supported by their coexpression in tumors and circulating tumor cells from CRPC patients. Importantly, HoxB13 silencing significantly decreases CRPC growth through inhibition of AR-V7 oncogenic function. These results identify HoxB13 as a pivotal upstream regulator of AR-V7-driven transcriptomes that are often cell context-dependent in CRPC, suggesting that HoxB13 may serve as a therapeutic target for AR-V7-driven prostate tumors.

Alternative splicing is widely recognized for playing roles in regulating genes and creating gene diversity. Consequently the identification and quantification of differentially spliced transcripts are pivotal for transcriptome analysis. However, how these diversified isoforms are spliced during genomic transcription and protein expression and what biological factors might influence the regulation of this are still required for further exploration. The advances in next-generation sequencing of messenger RNA (RNA-seq) have enabled us to survey gene expression and splicing more accurately. We have introduced a novel computational method, graph-based exon-skipping scanner (GESS), for de novo detection of skipping event sites from raw RNA-seq reads without prior knowledge of gene annotations, as well as for determining the dominant isoform generated from such sites. We have applied our method to publicly available RNA-seq data in GM12878 and K562 cells from the ENCODE consortium, and integrated other sequencing-based genomic data to investigate the impact of splicing activities, transcription factors (TFs) and epigenetic histone modifications on splicing outcomes. In a separate study, we also apply this algorithm in prostate cancer in The Cancer Genomics Atlas (TCGA) for de novo skipping event discovery to the understanding of abnormal splicing in each patient and to identify potential markers for prediction and progression of diseases.

RNA-Seq is the leading technology for analyzing gene expression on a global scale across a broad spectrum of sample types. However, due to chemical modifications by fixation or degradation due to collection methods, samples often contain an abundance of RNA that is no longer intact, and the capability of current RNA-Seq protocols to accurately quantify such samples is often limited. We have developed an RNA-Seq protocol to address these key issues as well as quantify gene expression from the whole transcriptome. Furthermore, for compatibility with improved sequencing platforms, we use restructured adapter sequences to generate libraries for Illumina HiSeq, MiSeq, and NextSeq platforms. Our protocol utilizes duplex-specific nuclease (DSN) to remove abundant ribosomal RNA sequences while retaining other types of RNA for superior transcriptome profiling from low quantity input. We employ the Illumina sequencing platform, but this method is described in sufficient detail to adapt to other platforms.

// Guangcun Huang 1 , Pawel A. Osmulski 1 , Hakim Bouamar 2 , Devalingam Mahalingam 3 , Chun-Lin Lin 1 , Michael A. Liss 4 , Addanki Pratap Kumar 4, 5 , Chun-Liang Chen 1 , Ian M. Thompson 4 , Lu-Zhe Sun 2 , Maria E. Gaczynska 1 , Tim H.-M. Huang 1 1 Departments of Molecular Medicine Cancer Research and Therapy Center and School of Medicine, University of Texas Health Science Center at San Antonio, San Antonio, Texas, USA 2 Departments of Cellular and Structural Biology Cancer Research and Therapy Center and School of Medicine, University of Texas Health Science Center at San Antonio, San Antonio, Texas, USA 3 Departments of Medicine Cancer Research and Therapy Center and School of Medicine, University of Texas Health Science Center at San Antonio, San Antonio, Texas, USA 4 Departments of Urology Cancer Research and Therapy Center and School of Medicine, University of Texas Health Science Center at San Antonio, San Antonio, Texas, USA 5 Departments of Radiation Oncology Cancer Research and Therapy Center and School of Medicine, University of Texas Health Science Center at San Antonio, San Antonio, Texas, USA Correspondence to: Tim H.-M. Huang, email: huangt3@uthscsa.edu Maria E. Gaczynska, email: gaczynska@uthscsa.edu Keywords: tumor metastasis, circulating tumor cells, epithelial-mesenchymal transition (EMT), transforming growth factor-β (TGF-β), positive feedback signaling Received: August 08, 2016 Accepted: October 12, 2016 Published: October 21, 2016 ABSTRACT Activation of TGF-β signaling is known to promote epithelial-mesenchymal transition (EMT) for the development of metastatic castration-resistant prostate cancer (mCRPC). To determine whether targeting TGF-β signaling alone is sufficient to mitigate mCRPC, we used the CRISPR/Cas9 genome-editing approach to generate a dominant-negative mutation of the cognate receptor TGFBRII that attenuated TGF-β signaling in mCRPC cells. As a result, the delicate balance of oncogenic homeostasis is perturbed, profoundly uncoupling proliferative and metastatic potential of TGFBRII -edited tumor xenografts. This signaling disturbance triggered feedback rewiring by enhancing ERK signaling known to promote EMT-driven metastasis. Circulating tumor cells displaying upregulated EMT genes had elevated biophysical deformity and an increase in interactions with chaperone macrophages for facilitating metastatic extravasation. Treatment with an ERK inhibitor resulted in decreased aggressive features of CRPC cells in vitro . Therefore, combined targeting of TGF-β and its backup partner ERK represents an attractive strategy for treating mCRPC patients.

Background The Pap smear has remained the foundation for cervical cancer screening for over 70 years. With advancements in molecular diagnostics, primary high-risk human papillomavirus (hrHPV) screening has recently become an accepted stand-alone or co-test with conventional cytology. However, both diagnostic tests have distinct limitations. The aim of this study was to determine the association between HPV genotypes and cellular epigenetic modifications in three grades of cervical cytology for screening biomarker discovery. Methods This prospective, cross-sectional study used residual liquid-based cytology samples for HPV genotyping and epigenetic analysis. Extracted DNA was subjected to parallel polymerase chain reactions using three primer sets (MY09/11, FAP59/64, E6-E7 F/B) for HPV DNA amplification. HPV+ samples were genotyped by DNA sequencing. Promoter methylation of four candidate tumor suppressor genes (adenylate cyclase 8 (ADCY8), cadherin 8, type 2 (CDH8), MGMT, and zinc finger protein 582 (ZNF582)) out of 48 genes screened was quantified by bisulfite-pyrosequencing of genomic DNA. Independent validation of methylation profiles was performed by analyzing data from cervical cancer cell lines and clinical samples from The Cancer Genome Atlas (TCGA). Results Two hundred seventy-seven quality cytology samples were analyzed. HPV was detected in 31/100 (31 %) negative for intraepithelial lesion or malignancy (NILM), 95/100 (95 %) low-grade squamous intraepithelial lesion (LSIL), and 71/77 (92 %) high-grade squamous intraepithelial lesion (HSIL) samples. The proportion of IARC-defined carcinogenic HPV types in sequenced samples correlated with worsening grade: NILM 7/29 (24 %), LSIL 53/92 (58 %), and HSIL 65/70 (93 %). Promoter methylation of ADCY8, CDH8, and ZNF582 was measured in 170 samples: NILM (N = 33), LSIL (N = 70), and HSIL (N = 67) also correlated with worsening grade. Similar hypermethylation patterns were found in cancer cell lines and TCGA samples. The combination of four biomarkers, i.e., HPV genotype and three-gene promoter methylation, predicted HSIL (AUC 0.89) better than HPV alone (AUC 0.74) by logistic regression and probabilistic modeling. Conclusions HPV genotype and DNA methylation of ADCY8, CDH8, and ZNF582 are correlated with cytological grade. Collectively, these biomarkers may serve as a molecular classifier of Pap smears. Electronic supplementary material The online version of this article (doi:10.1186/s13148-016-0263-9) contains supplementary material, which is available to authorized users.

Overexpression of epithelial cell adhesion molecule (EpCAM) has been implicated in advanced endometrial cancer, but its roles in this progression remain to be elucidated. In addition to its structural role in modulating cell-surface adhesion, here we demonstrate that EpCAM is a regulatory molecule in which its internalization into the nucleus turns on a transcription program. Activation of EGF/EGFR signal transduction triggered cell-surface cleavage of EpCAM, leading to nuclear internalization of its cytoplasmic domain EpICD. ChIP-seq analysis identified target genes that are coregulated by EpICD and its transcription partner, LEF-1. Network enrichment analysis further uncovered a group of 105 genes encoding functions for tight junction, adherent, and cell migration. Furthermore, nanomechanical analysis by atomic force microscopy revealed increased softness and decreased adhesiveness of EGF-stimulated cancer cells, implicating acquisition of an epithelial–mesenchymal transition (EMT) phenotype. Thus, genome editing of EpCAM could be associated with altering these nanomechanical properties towards a less aggressive phenotype. Using this integrative genomic–biophysical approach, we demonstrate for the first time an intricate relationship between EpCAM-regulated transcription and altered biophysical properties of cells that promote EMT in advanced endometrial cancer. Cancer Res; 76(21); 6171–82. ©2016 AACR.

Aberrant DNA methylation at promoters is often linked to tumorigenesis. But many aspects of DNA methylation remain unexplored, including the individual roles of distal and gene body methylation, as well as their collaborative roles with promoter methylation. Here we performed a MBD-seq analysis on prostate specimens classified into low, high and very high risk group based on Gleason score and TNM stages. We identified gene sets with differential methylation regions (DMRs) in Distal, TSS, gene body and TES. To understand the collaborative roles, TSS was compared with the other three DMRs, resulted in 12 groups of genes with collaborative differential methylation patterns (CDMPs). We found several groups of genes that show opposite methylation patterns in Distal and Genic regions compared to TSS region and in general they are differentially expressed genes (DEGs) in tumors in TCGA RNA-seq data. IPA (Ingenuity Pathway Analysis) reveals AR/TP53 signaling network to be a major signaling pathway and survival analysis indicates genes subsets significantly associated with prostate cancer recurrence. Our results suggest that DNA methylation in Distal and Genic regions also plays critical roles in contributing to prostate tumorigenesis and may act either positively or negatively with TSSs to alter gene regulation in tumors.

Several mutual information (MI)-based algorithms have been developed to identify dynamic gene-gene and function-function interactions governed by key modulators (genes, proteins, etc.). Due to intensive computation, however, these methods rely heavily on prior knowledge and are limited in genome-wide analysis. We present the modulated gene/gene set interaction (MAGIC) analysis to systematically identify genome-wide modulation of interaction networks. Based on a novel statistical test employing conjugate Fisher transformations of correlation coefficients, MAGIC features fast computation and adaption to variations of clinical cohorts. In simulated datasets MAGIC achieved greatly improved computation efficiency and overall superior performance than the MI-based method. We applied MAGIC to construct the estrogen receptor (ER) modulated gene and gene set (representing biological function) interaction networks in breast cancer. Several novel interaction hubs and functional interactions were discovered. ER+ dependent interaction between TGFβ and NFκB was further shown to be associated with patient survival. The findings were verified in independent datasets. Using MAGIC, we also assessed the essential roles of ER modulation in another hormonal cancer, ovarian cancer. Overall, MAGIC is a systematic framework for comprehensively identifying and constructing the modulated interaction networks in a whole-genome landscape. MATLAB implementation of MAGIC is available for academic uses at https://github.com/chiuyc/MAGIC.

Aberrant DNA methylation is a hallmark of many cancers. Currently, there are four intrinsic molecular subtypes in breast cancer (BC): Luminal A, B, Her2-positive, and triple negative (TNBC). Recently, The Cancer Genome Atlas (TCGA) project has revealed that Luminal subtypes have higher levels of genome-wide methylation that may be a result of Estrogen/Estrogen receptor α (E2/ERα) signaling pathway activation. In this study, we analyze promoter CpG-island (CGIs) of the Reprimo (RPRM) gene in breast cancers (n = 77), cell lines (n = 38), and normal breast tissue (n = 10) using a MBDCap-seq database. Then, a validation cohort (n = 26) was used to confirm the results found in the MBDCap-seq platform. A differential methylation pattern was found between BC and cell lines compared to normal breast tissue. In BC, a higher DNA methylation was observed in tissues that were ERα-positive than in ERα-negative ones; more precisely, subtypes Luminal A compared to TNBC. Also, significant reverse correlation was observed between DNA methylation and RPRM mRNA expression in BC. Our data suggest that ERα expression in BC may affect the DNA methylation of CGIs in the RPRM gene. This approach suggests that DNA methylation status in CGIs of some tumor suppressor genes could be driven by E2 availability, subsequently inducing the activation of the ERα pathway.

Early postnatal exposures to Bisphenol A (BPA) and genistein (GEN) have been reported to predispose for and against mammary cancer, respectively, in adult rats. Since the changes in cancer susceptibility occurs in the absence of the original chemical exposure, we have investigated the potential of epigenetics to account for these changes. DNA methylation studies reveal that prepubertal BPA exposure alters signaling pathways that contribute to carcinogenesis. Prepubertal exposure to GEN and BPA + GEN revealed pathways involved in maintenance of cellular function, indicating that the presence of GEN either reduces or counters some of the alterations caused by the carcinogenic properties of BPA. We subsequently evaluated the potential of epigenetic changes in the rat mammary tissues to predict survival in breast cancer patients via the Cancer Genomic Atlas (TCGA). We identified 12 genes that showed strong predictive values for long-term survival in estrogen receptor positive patients. Importantly, two genes associated with improved long term survival, HPSE and RPS9, were identified to be hypomethylated in mammary glands of rats exposed prepuberally to GEN or to GEN + BPA respectively, reinforcing the suggested cancer suppressive properties of GEN.

Identifying prostate cancer-driving transcription factors (TFs) in addition to the androgen receptor promises to improve our ability to effectively diagnose and treat this disease. We employed an integrative genomics analysis of master TFs CREB1 and FoxA1 in androgen-dependent prostate cancer (ADPC) and castration-resistant prostate cancer (CRPC) cell lines, primary prostate cancer tissues and circulating tumor cells (CTCs) to investigate their role in defining prostate cancer gene expression profiles. Combining genome-wide binding site and gene expression profiles we define CREB1 as a critical driver of pro-survival, cell cycle and metabolic transcription programs. We show that CREB1 and FoxA1 co-localize and mutually influence each other's binding to define disease-driving transcription profiles associated with advanced prostate cancer. Gene expression analysis in human prostate cancer samples found that CREB1/FoxA1 target gene panels predict prostate cancer recurrence. Finally, we showed that this signaling pathway is sensitive to compounds that inhibit the transcription co-regulatory factor MED1. These findings not only reveal a novel, global transcriptional co-regulatory function of CREB1 and FoxA1, but also suggest CREB1/FoxA1 signaling is a targetable driver of prostate cancer progression and serves as a biomarker of poor clinical outcomes.

We would like to thank Dr. Gires for the letter concerning our article ([1][1]). We are very glad that Dr. Gires considers our data on how EGFR activation affects endometrial cancer cell invasiveness via EpCAM activity novel and interesting. The study highlights not only the function of EpCAM as a

Background Human papillomavirus (HPV) is the carcinogen of almost all invasive cervical cancer and a major cause of oral and other anogenital malignancies. HPV genotyping by dideoxy (Sanger) sequencing is currently the reference method of choice for clinical diagnostics. However, for samples with multiple HPV infections, genotype identification is singular and occasionally imprecise or indeterminable due to overlapping chromatograms. Our aim was to explore and compare HPV metagenomes in abnormal cervical cytology by deep sequencing for correlation with disease states. Results Low- and high-grade intraepithelial lesion (LSIL and HSIL) cytology samples were DNA extracted for PCR-amplification of the HPV E6/E7 genes. HPV+ samples were sequenced by dideoxy and deep methods. Deep sequencing revealed ~60% of all samples (n = 72) were multi-HPV infected. Among LSIL samples (n = 43), 27 different genotypes were found. The 3 dominant (most abundant) genotypes were: HPV-39, 11/43 (26%); -16, 9/43 (21%); and -35, 4/43 (9%). Among HSIL (n = 29), 17 HPV genotypes were identified; the 3 dominant genotypes were: HPV-16, 21/29 (72%); -35, 4/29 (14%); and -39, 3/29 (10%). Phylogenetically, type-specific E6/E7 genetic distances correlated with carcinogenic potential. Species diversity analysis between LSIL and HSIL revealed loss of HPV diversity and domination by HPV-16 in HSIL samples. Conclusions Deep sequencing resolves HPV genotype composition within multi-infected cervical cytology. Biodiversity analysis reveals loss of diversity and gain of dominance by carcinogenic genotypes in high-grade cytology. Metagenomic profiles may therefore serve as a biomarker of disease severity and a population surveillance tool for emerging genotypes. Electronic supplementary material The online version of this article (doi:10.1186/s12864-017-3612-y) contains supplementary material, which is available to authorized users.