Search

To explore the feature of cancer cells and tumor subclones, we analyzed 101,065 single-cell transcriptomes from 12 colorectal cancer (CRC) patients and 92 single cell genomes from one of these patients. We found cancer cells, endothelial cells and stromal cells in tumor tissue expressed much more genes and had stronger cell-cell interactions than their counterparts in normal tissue. We identified copy number variations (CNVs) in each cancer cell and found correlation between gene copy number and expression level in cancer cells at single cell resolution. Analysis of tumor subclones inferred by CNVs showed accumulation of mutations in each tumor subclone along lineage trajectories. We found differentially expressed genes (DEGs) between tumor subclones had two populations: DEGCNVand DEGreg. DEGCNV, showing high CNV-expression correlation and whose expression differences depend on the differences of CNV level, enriched in housekeeping genes and cell adhesion associated genes. DEGreg, showing low CNV-expression correlation and mainly in low CNV variation regions and regions without CNVs, enriched in cytokine signaling genes. Furthermore, cell-cell communication analyses showed that DEGCNVtends to involve in cell-cell contact while DEGregtends to involve in secreted signaling, which further support that DEGCNVand DEGregare two regulatorily and functionally distinct categories.

Primers used in the study.

Lysine to methionine mutations at position 27 (K27M) in the histone H3 (H3.3 and H3.1) are highly prevalent in pediatric high-grade gliomas (HGG) that arise in the midline of the central nervous system. H3K27M perturbs the activity of polycomb repressor complex 2 and correlates with DNA hypomethylation; however, the pathways whereby H3K27M drives the development of pediatric HGG remain poorly understood. To understand the mechanism of pediatric HGG development driven by H3.3K27M and discover potential therapeutic targets or biomarkers, we established pediatric glioma cell model systems harboring H3.3K27M and performed microarray analysis. H3.3K27M caused the upregulation of multiple cancer/testis (CT) antigens, such as ADAMTS1, ADAM23, SPANXA1, SPANXB1/2, IL13RA2, VCY, and VCX3A, in pediatric glioma cells. Chromatin immunoprecipitation analysis from H3.3K27M cells revealed decreased H3K27me3 levels and increased H3K4me3 levels on the VCX3A promoter. Knockdown of VCX3A by siRNA significantly inhibited the growth of pediatric glioma cells harboring H3.3K27M. Overexpression of VCX3A/B genes stimulated the expression of several HLA genes, including HLA-A, HLA-B, HLA-E, HLA-F, and HLA-G. The expression of VCX3A in pediatric HGG was confirmed using a tissue microarray. Gene set enrichment analysis revealed that CT antigens are enriched in pediatric HGG clinical specimens with H3.3K27M, with the upregulation of IL13RA2 contributing to the enrichment significantly. These results indicate that the upregulation of CT antigens, such as VCX3A and IL13RA2, correlates with pediatric gliomagenesis. Mol Cancer Res; 16(4); 623–33. ©2018 AACR.

Supplementary Figure S1. H3.3K27M activates cancer/testis (CT) antigens in SF-188 cell line; Supplementary Figure S2. The putative CpGislands in the human VCX3A locus; Supplementary Figure S3. The genomic sequence of human glyceraldehyde 3-phosphate dehydrogenase (GAPDH), and an intragenic CpG rich region of GAPDH (S3A). No change in intragenic CpG methylation of GAPDH in Res259 cells carrying H3.3K27M (S3B); Supplementary Figure S4. Subcellular localization of VCX3A/B in Res259 cells; Supplementary Figure S5. The different isoforms of VCX3A/B cause similar gene expression changes; Supplementary Figure S6. The mRNA expression of HLA-A, -B, E, F and G in Res259 cells stably expressing H3.3K27M or empty vector.

A list of cancer-testis antigen genes: the custom Gene Matrix Transposed (GMT) file containing the list of cancer testis antigens used for GSEA in the study.

S1. The relative expression of KLK3/PSA. S2. Kaplan-Meier analysis of prostate cancer outcome was obtained from GEPIA tool. S3. LNCaP cells with the indicated treatments were seeded in the upper chamber of the Boyden's transwell chamber. S4. LNCaP cells with the indicated treatments were seeded in the upper chamber of the Boyden's transwell chamber which were coated with 1:20 diluted Matrigel. S5. LNCaP cells with the indicated treatments were seeded in the upper chamber of the Boyden's transwell chamber. S6. Snapshots form the AR ChIP-seq showing AR binding to the ARBS. S7. Screenshot from the AR ChIP-seq shows marked reduction in AR binding to the ARBS located at ~ 30kb from the NDRG1 after LINC00844 knock-down in LNCaP cells. S8. LNCaP cells with the indicated treatments were seeded in the upper chamber of the Boyden's transwell chamber. S9. LNCaP cells we co-transfected with the with the indicated siRNA and plasmids and seeded in the upper chamber of the Boyden's transwell chamber. S10. LNCaP cells we co-transfected with the with the indicated siRNA and plasmids cells were harvested and the required amount of cells were used for the migration and invasion assay and rest were lysed for western blot analysis.

The human genome is mostly transcribed, yielding a rich repository of noncoding transcripts that are involved in a myriad of biological processes including cancer. However, how many noncoding transcripts such as long noncoding RNAs (lncRNA) function in cancer is still unclear. This study identified a novel set of clinically relevant androgen-regulated lncRNAs in prostate cancer. Among this group, LINC00844 was demonstrated to be a direct androgen-regulated target that is actively transcribed in androgen receptor (AR)–dependent prostate cancer cells. The expression of LINC00844 is higher in normal prostate compared with malignant and metastatic prostate cancer clinical specimens, and patients with low expression had a poor prognosis and significantly increased biochemical recurrence, suggesting LINC00844 functions in suppressing tumor progression and metastasis. Indeed, in vitro loss-of-function studies revealed that LINC00844 prevents prostate cancer cell migration and invasion. Moreover, findings from gene expression profiling analysis indicated that LINC00844 functions in trans, affecting global androgen-regulated gene transcription. Mechanistic evidence reveals that LINC00844 is important in facilitating AR binding to the chromatin. Finally, LINC00844 mediates its phenotypic effects in part by activating the expression of NDRG1, a crucial cancer metastasis suppressor. Collectively, LINC00844 is a novel coregulator of AR that plays a central role in the androgen transcriptional network and the development and progression of prostate cancer.Implications:This study highlights the function of the lncRNA, LINC00844, in regulating global AR-regulated genes in prostate cancer by modulating AR binding to chromatin.

S1. Differentially regulated lncRNAs FDR

The sequence analysis of different variants of VCX3A and VCX3B.

Activator Protein 2 gamma (AP-2γ) is a master transcription factor that plays a critical role in the development and progression of breast cancer. However, the underlying mechanism is still unclear. Herein, using a proteomics approach, we identified Tripartite motif-containing 37 (TRIM37) as a novel coactivator of AP-2γ-mediated transcription in breast cancer cells. We demonstrate that TRIM37 facilitates AP-2γ chromatin binding to directly regulate the AP-2γ mediated transcriptional program. We also show that TRIM37 achieves this by stimulating K63 chain-linked ubiquitination of AP-2γ, promoting protein localization from the cytoplasm to the nucleus. In clinical analyses, we find TRIM37 is upregulated in multiple breast cancer datasets, supporting our findings that the TRIM37-AP-2γ interaction is essential for breast cancer tumor growth. Overall, our work reveals that TRIM37 is an oncogenic coactivator of AP-2γ in breast cancer and provides a novel therapeutic target for treating the disease.

Androgen receptor (AR) plays a central role in driving prostate cancer (PCa) progression. How AR promotes this process is still not completely clear. Herein, we used single-cell transcriptome analysis to reconstruct the transcriptional network of AR in PCa. Our work shows AR directly regulates a set of signature genes in the ER-to-Golgi protein vesicle-mediated transport pathway. The expression of these genes is required for maximum androgen-dependent ER-to-Golgi trafficking, cell growth, and survival. Our analyses also reveal the signature genes are associated with PCa progression and prognosis. Moreover, we find inhibition of the ER-to-Golgi transport process with a small molecule enhanced antiandrogen-mediated tumor suppression of hormone-sensitive and insensitive PCa. Finally, we demonstrate AR collaborates with CREB3L2 in mediating ER-to-Golgi trafficking in PCa. In summary, our findings uncover a critical role for dysregulation of ER-to-Golgi trafficking expression and function in PCa progression, provide detailed mechanistic insights for how AR tightly controls this process, and highlight the prospect of targeting the ER-to-Golgi pathway as a therapeutic strategy for advanced PCa.

Spatial omics is a rapidly evolving approach for exploring tissue microenvironment and cellular networks by integrating spatial knowledge with transcript or protein expression information. However, there is a lack of databases for users to access and analyze spatial omics data. To address this limitation, we developed Aquila, a comprehensive platform for managing and analyzing spatial omics data. Aquila contains 107 datasets from 30 diseases, including 6500+ regions of interest, and 15.7 million cells. The database covers studies from spatial transcriptome and proteome analyses, 2D and 3D experiments, and different technologies. Aquila provides visualization of spatial omics data in multiple formats such as spatial cell distribution, spatial expression and co-localization of markers. Aquila also lets users perform many basic and advanced spatial analyses on any dataset. In addition, users can submit their own spatial omics data for visualization and analysis in a safe and secure environment. Finally, Aquila can be installed as an individual app on a desktop and offers the RESTful API service for power users to access the database. Overall, Aquila provides a detailed insight into transcript and protein expression in tissues from a spatial perspective. Aquila is available at https://aquila.cheunglab.org.

Animal interphase chromosomes are organized into topologically associating domains (TADs). How TADs are formed is not fully understood. Here, we combined high-throughput chromosome conformation capture and gene silencing to obtain insights into TAD dynamics in Xenopus tropicalis embryos. First, TAD establishment in X. tropicalis is similar to that in mice and flies and does not depend on zygotic genome transcriptional activation. This process is followed by further refinements in active and repressive chromatin compartments and the appearance of loops and stripes. Second, within TADs, higher self-interaction frequencies at one end of the boundary are associated with higher DNA occupancy of the architectural proteins CTCF and Rad21. Third, the chromatin remodeling factor ISWI is required for de novo TAD formation. Finally, TAD structures are variable in different tissues. Our work shows that X. tropicalis is a powerful model for chromosome architecture analysis and suggests that chromatin remodeling plays an essential role in de novo TAD establishment., Hi-C and single-molecule sequencing analysis provide an improved assembly of the Xenopus tropicalis genome and insights into three-dimensional genome dynamics throughout embryogenesis.

RNA N

Androgen receptor (AR) is a master transcription factor that drives prostate cancer (PCa) development and progression. Alterations in the expression or activity of AR coregulators significantly impact the outcome of the disease. Using a proteomics approach, we identified the tripartite motif‐containing 33 (TRIM33) as a novel transcriptional coactivator of AR. We demonstrate that TRIM33 facilitates AR chromatin binding to directly regulate a transcription program that promotes PCa progression. TRIM33 further stabilizes AR by protecting it from Skp2‐mediated ubiquitination and proteasomal degradation. We also show that TRIM33 is essential for PCa tumor growth by avoiding cell‐cycle arrest and apoptosis, and TRIM33 knockdown sensitizes PCa cells to AR antagonists. In clinical analyses, we find TRIM33 upregulated in multiple PCa patient cohorts. Finally, we uncover an AR‐TRIM33‐coactivated gene signature highly expressed in PCa tumors and predict disease recurrence. Overall, our results reveal that TRIM33 is an oncogenic AR coactivator in PCa and a potential therapeutic target for PCa treatment.

Earlier detection of biomarkers responsible for cancer relapse facilitates more rational cancer treatment regimens to be designed. Herein, we develop a mass cytometry-based strategy for unbiased mining of cell subsets that potentially contribute to cancer recurrence through panoramic examination of the immunophenotypic features and multidrug resistance characteristics. The incorporation of metal tags enables multiplexed information of single cells to be interrogated based on metal fingerprint. Using acute lymphoblastic leukemia (B-ALL) as a showcase, we show overexpressed multidrug resistance biomarkers, i. e., BCRP, Bcl-2, MRP1, and P-gp in B-ALL cells compared with healthy control, and a positive correlation among different multidrug resistance biomarkers. Different cell subsets with multidrug resistance are well-defined, featured with CD34+ CD38+ CD10- and CD34+ CD38+/int CD10+ . Importantly, we uncovered that CD34 expression level is positively correlated to multidrug resistance, indicative of a higher potential of immature cells to induce B-ALL relapse. In addition, the cell subsets positively expressing CD73 and CD304 (CD34+ CD10+ CD304+ ; CD34+ CD38+/int CD10+ CD73+ ) also overexpress multidrug resistance biomarkers, suggesting that they may serve as additional new biomarkers for B-ALL stratification and prognosis. Our data provide the first evidence that highly expressed multidrug resistance biomarkers in certain cell subpopulations with specific immunophenotypes may potentially induce B-ALL recurrence. The incorporation of multidrug resistance features with cell phenotypes using mass cytometry proposed in this study provides a general strategy for risk assessment and the prediction of recurrence of different types of cancers.

We elucidated the anti-inflammatory mechanisms of IL-38 in allergic asthma. Human bronchial epithelial cells and eosinophils were cocultured upon stimulation with the viral RLR ligand poly (I:C)/LyoVec or infection-related cytokine TNF-α to induce expression of cytokines/chemokines/adhesion molecules. House dust mite (HDM)-induced allergic asthma and humanized allergic asthma NOD/SCID murine models were established to assess anti-inflammatory mechanisms in vivo. IL-38 significantly inhibited induced proinflammatory IL-6, IL-1β, CCL5, and CXCL10 production, and antiviral interferon-β and intercellular adhesion molecule-1 expression in the coculture system. Mass cytometry and RNA-sequencing analysis revealed that IL-38 could antagonize the activation of the intracellular STAT1, STAT3, p38 MAPK, ERK1/2, and NF-κB pathways, and upregulate the expression of the host defense-related gene POU2AF1 and anti-allergic response gene RGS13. Intraperitoneal injection of IL-38 into HDM-induced allergic asthma mice could ameliorate airway hyperreactivity by decreasing the accumulation of eosinophils in the lungs and inhibiting the expression of the Th2-related cytokines IL-4, IL-5, and IL-13 in the bronchoalveolar lavage fluid (BALF) and lung homogenates. Histological examination indicated lung inflammation was alleviated by reductions in cell infiltration and goblet cell hyperplasia, together with reduced Th2, Th17, and innate lymphoid type 2 cell numbers but increased proportions of regulatory T cells in the lungs, spleen, and lymph nodes. IL-38 administration suppressed airway hyperreactivity and asthma-related IL-4 and IL-5 expression in humanized mice, together with significantly decreased CCR3+ eosinophil numbers in the BALF and lungs, and a reduced percentage of human CD4+CRTH2+ Th2 cells in the lungs and mediastinal lymph nodes. Together, our results demonstrated the anti-inflammatory mechanisms of IL-38 and provided a basis for the development of a regulatory cytokine-based treatment for allergic asthma.

Summary Pyruvate is a key metabolite in glycolysis and the tricarboxylic acid (TCA) cycle. Exogenous pyruvate modulates metabolism, provides cellular protection, and is essential for the maintenance of human preimplantation embryos and human embryonic stem cells (hESCs). However, little is known about how pyruvate contributes to cell-fate determination during epiblast stage. In this study, we used hESCs as a model to demonstrate that elevated exogenous pyruvate shifts metabolic balance toward oxidative phosphorylation in both maintenance and differentiation conditions. During differentiation, pyruvate potentiates mesoderm and endoderm lineage specification. Pyruvate production and its mitochondrial metabolism are required in BMP4-induced mesoderm differentiation. However, the TCA-cycle metabolites do not have the same effect as pyruvate on differentiation. Further study shows that pyruvate increases AMP/ATP ratio, activates AMPK, and modulates the mTOR pathway to enhance mesoderm differentiation. This study reveals that exogenous pyruvate not only controls metabolism but also modulates signaling pathways in hESC differentiation., Graphical Abstract, Highlights • Exogenous pyruvate enhances TCA-cycle metabolism in hESCs • Exogenous pyruvate promotes BMP4-induced mesoderm differentiation • Exogenous pyruvate influences differentiation through AMPK/mTOR pathways, In this report, Chen and colleagues demonstrate that pyruvate promotes mesoderm differentiation in human embryonic stem cells. Pyruvate elevates TCA-cycle metabolism and modulates AMPK and mTOR pathways to promote BMP4-induced mesoderm differentiation. This discovery highlights that pyruvate as a common metabolite can regulate signaling pathways in cell-fate determination.

RNA N6-Methyladenosine (m6A) is the most abundant mRNA modification, and forms part of an epitranscriptomic system that modulates RNA function. RNA modifications can be reversibly catalyzed by several specific enzymes, and those modifications can be recognized by RNA binding proteins that in turn regulate biological processes. Although there are many reports demonstrating m6A participation in critical biological functions, this exploration has mainly been conducted through the global knockout or knockdown of the writers, erasers, or readers of m6A. Consequently, there is a lack of information about the role of m6A on single transcripts in biological processes, posing a challenge in understanding the biological functions of m6A. Here, we demonstrate a CRISPR/dCas13a-based RNA m6A-editor which can target mRNAs using single crRNA or multiple crRNAs array to methylate or demethylate m6A. We systematically assay its capabilities to enable the targeted rewriting of m6A dynamics, including modulation of circular RNA translation and transcript half-life. Finally, we demonstrate the utility of the system by specifically modulating XIST m6A levels, which can control X chromosome silencing and activation. Based on our editors, m6A on single and multiple transcripts can be modified to allow the exploration of the role of m6A on in biological processes.

Metazoan genomes are folded into 3D structures in interphase nuclei. However, the molecular mechanism remains unknown. Here, we show that topologically associating domains (TADs) form in two waves during Xenopus tropicalis embryogenesis, first at zygotic genome activation and then as the expression of CTCF and Rad21 is elevated. We also found TAD structures continually change for at least three times during development. Surprisingly, the directionality index is preferentially stronger on one side of TADs where orientation-biased CTCF and Rad21 binding are observed, a conserved pattern that is found in human cells as well. Depletion analysis revealed CTCF, Rad21, and RPB1, a component of RNAPII, are required for the establishment of TADs. Overall, our work shows that Xenopus is a powerful model for chromosome architecture analysis. Furthermore, our findings indicate that cohesin-mediated extrusion may anchor at orientation-biased CTCF binding sites, supporting a CTCF-anchored extrusion model as the mechanism for TAD establishment.

The human genome is mostly transcribed, yielding a rich repository of noncoding transcripts that are involved in a myriad of biological processes including cancer. However, how many noncoding transcripts such as long noncoding RNAs (lncRNA) function in cancer is still unclear. This study identified a novel set of clinically relevant androgen-regulated lncRNAs in prostate cancer. Among this group, LINC00844 was demonstrated to be a direct androgen-regulated target that is actively transcribed in androgen receptor (AR)–dependent prostate cancer cells. The expression of LINC00844 is higher in normal prostate compared with malignant and metastatic prostate cancer clinical specimens, and patients with low expression had a poor prognosis and significantly increased biochemical recurrence, suggesting LINC00844 functions in suppressing tumor progression and metastasis. Indeed, in vitro loss-of-function studies revealed that LINC00844 prevents prostate cancer cell migration and invasion. Moreover, findings from gene expression profiling analysis indicated that LINC00844 functions in trans, affecting global androgen-regulated gene transcription. Mechanistic evidence reveals that LINC00844 is important in facilitating AR binding to the chromatin. Finally, LINC00844 mediates its phenotypic effects in part by activating the expression of NDRG1, a crucial cancer metastasis suppressor. Collectively, LINC00844 is a novel coregulator of AR that plays a central role in the androgen transcriptional network and the development and progression of prostate cancer. Implications: This study highlights the function of the lncRNA, LINC00844, in regulating global AR-regulated genes in prostate cancer by modulating AR binding to chromatin.

Cell density has profound impacts on the cell culture practices of human pluripotent stem cells. The regulation of cell growth, cell death, pluripotency and differentiation converge at high density, but it is largely unknown how different regulatory mechanisms act at this stage. We use a chemically defined medium to systemically examine cellular activities and the impact of medium components in high-density culture. We show that medium acidosis is the main factor that alters cell cycle, gene expression and cellular metabolism at high cell density. The low medium pH leads to inhibition of glucose consumption, cell cycle arrest, and subsequent cell death. At high cell density, the suppression of medium acidosis with sodium bicarbonate (NaHCO3) significantly increases culture capacity for stem cell survival, derivation, maintenance and differentiation. Our study provides a simple and effective tool to improve stem cell maintenance and applications.

Identification of prognostic biomarkers helps facilitate the prediction of patient outcomes as well as guide treatments. Accumulating evidence now suggests that long non-coding RNAs (lncRNAs) play key roles in tumor progression with diagnostic and prognostic values. However, little is known about the biological functions of lncRNAs and how they contribute to the pathogenesis of cancer. Herein, we performed weighted correlation network analysis (WGCNA) on 380 RNA-seq samples from prostate cancer patients to create networks comprising of microRNAs, lncRNAs, and protein-coding genes. Our analysis revealed expression modules that associated with pathological parameters. More importantly, we identified a gene module that is involved in protein translation and is associated with patient survival. In this gene module, we explored the regulation axis involving GAS5, ZFAS1, and miR-940. We show that GAS5, ZFAS1, and miR-940 are up-regulated in tumors relative to normal prostate tissues, and high expression of either lncRNA is an indicator of poor patient outcome. Finally, we constructed a co-expression network involving GAS5, ZFAS1, and miR-940, as well as the targets of miR-940. Our results show that GAS5 and ZFAS1 are targeted by miR-940 via NAA10 and RPL28. Taken together, co-expression analysis of gene expression profiling from RNA-seq can accelerate the identification and functional characterization of novel prognostic markers in prostate cancer.

The aberrant activities of transcription factors such as the androgen receptor (AR) underpin prostate cancer development. While the AR cis-regulation has been extensively studied in prostate cancer, information pertaining to the spatial architecture of the AR transcriptional circuitry remains limited. In this paper, we propose a novel framework to profile long-range chromatin interactions associated with AR and its collaborative transcription factor, erythroblast transformation-specific related gene (ERG), using chromatin interaction analysis by paired-end tag (ChIA-PET). We identified ERG-associated long-range chromatin interactions as a cooperative component in the AR-associated chromatin interactome, acting in concert to achieve coordinated regulation of a subset of AR target genes. Through multifaceted functional data analysis, we found that AR-ERG interaction hub regions are characterized by distinct functional signatures, including bidirectional transcription and cotranscription factor binding. In addition, cancer-associated long noncoding RNAs were found to be connected near protein-coding genes through AR-ERG looping. Finally, we found strong enrichment of prostate cancer genome-wide association study (GWAS) single nucleotide polymorphisms (SNPs) at AR-ERG co-binding sites participating in chromatin interactions and gene regulation, suggesting GWAS target genes identified from chromatin looping data provide more biologically relevant findings than using the nearest gene approach. Taken together, our results revealed an AR-ERG-centric higher-order chromatin structure that drives coordinated gene expression in prostate cancer progression and the identification of potential target genes for therapeutic intervention.

How deregulation of chromatin modifiers causes malignancies is of general interest. Here, we show that histone H2A T120 is phosphorylated in human cancer cell lines and demonstrate that this phosphorylation is catalyzed by hVRK1. Cyclin D1 was one of ten genes downregulated upon VRK1 knockdown in two different cell lines and showed loss of H2A T120 phosphorylation and increased H2A K119 ubiquitylation of its promoter region, resulting in impaired cell growth. In vitro, H2A T120 phosphorylation and H2A K119 ubiquitylation are mutually inhibitory, suggesting that histone phosphorylation indirectly activates chromatin. Furthermore, expression of a phosphomimetic H2A T120D increased H3 K4 methylation. Finally, both VRK1 and the H2A T120D mutant histone transformed NIH/3T3 cells. These results suggest that histone H2A T120 phosphorylation by hVRK1 causes inappropriate gene expression, including upregulated cyclin D1, which promotes oncogenic transformation.

Assay for Transposase Accessible Chromatin with high-throughput sequencing (ATAC-seq) is a powerful genomic technology that is used for the global mapping and analysis of open chromatin regions. However, for users to process and analyze such data they either have to use a number of complicated bioinformatic tools or attempt to use the currently available ATAC-seq analysis software, which are not very user friendly and lack visualization of the ATAC-seq results. Because of these issues, biologists with minimal bioinformatics background who wish to process and analyze their own ATAC-seq data by themselves will find these tasks difficult and ultimately will need to seek help from bioinformatics experts. Moreover, none of the available tools provide complete solution for ATAC-seq data analysis. Therefore, to enable non-programming researchers to analyze ATAC-seq data on their own, we developed a tool called Graphical User interface for the Analysis and Visualization of ATAC-seq data (GUAVA). GUAVA is a standalone software that provides users with a seamless solution from beginning to end including adapter trimming, read mapping, the identification and differential analysis of ATAC-seq peaks, functional annotation, and the visualization of ATAC-seq results. We believe GUAVA will be a highly useful and time-saving tool for analyzing ATAC-seq data for biologists with minimal or no bioinformatics background. Since GUAVA can also operate through command-line, it can easily be integrated into existing pipelines, thus providing flexibility to users with computational experience.

The identification and functional characterization of novel biomarkers in cancer requires survival analysis and gene expression analysis of both patient samples and cell line models. To help facilitate this process, we have developed KM-Express. KM-Express holds an extensive manually curated transcriptomic profile of 45 different datasets for prostate and breast cancer with phenotype and pathoclinical information, spanning from clinical samples to cell lines. KM-Express also contains The Cancer Genome Atlas datasets for 30 other cancer types with matching cell line expression data for 23 of them. We present KM-Express as a hypothesis generation tool for researchers to identify potential new prognostic RNA biomarkers as well as targets for further downstream functional cell-based studies. Specifically, KM-Express allows users to compare the expression level of genes in different groups of patients based on molecular, genetic, clinical and pathological status. Moreover, KM-Express aids the design of biological experiments based on the expression profile of the genes in different cell lines. Thus, KM-Express provides a one-stop analysis from bench work to clinical prospects. We have used this tool to successfully evaluate the prognostic potential of previously published biomarkers for prostate cancer and breast cancer. We believe KM-Express will accelerate the translation of biomedical research from bench to bed. Database URL: http://ec2-52-201-246-161.compute-1.amazonaws.com/kmexpress/index.php

Objectives Drug combination in cancer therapy aims to achieve synergistic therapeutic effect, reduced drug dosage, reduced drug toxicity and minimizes or delays the induction of drug resistance. In the present study, we investigated the anticancer effects of the combination of two metabolic modulators, dichloroacetate (DCA) and bacillus caldovelox arginase (BCA) (or pegyated human arginase (HA)). Methods The combination treatments were evaluated in MCF-7 and MDA-MB 231 cells as well as in MDA-MB 231 breast cancer xenograft model. Key findings Dichloroacetate and BCA combination exhibited anti-proliferative effects on MCF-7 cells, which were found to be synergistic. Analysis of the gene expression upon drug treatments revealed that the synergistic anti-proliferative effect on MCF-7 cells was possibly in part due to the activation of the p53 pathway. A similar synergistic anti-proliferative effect was observed in the combined use of DCA and HA on MCF-7 and MDA-MB231 cells, which was due to induction of cell cycle arrest at G2/M phase. Moreover, the combination enhanced anti-tumour activity in a MDA-MB 231 xenograft mouse model. Conclusions Our results suggested that dichloroacetate and arginase combination exhibited enhanced anti-cancer effects in preclinical breast cancer models which may offer an additional treatment option for breast cancer.

Despite altered metabolism being an accepted hallmark of cancer, it is still not completely understood which signaling pathways regulate these processes. Given the central role of androgen receptor (AR) signaling in prostate cancer, we hypothesized that AR could promote prostate cancer cell growth in part through increasing glucose uptake via the expression of distinct glucose transporters. Here, we determined that AR directly increased the expression ofSLC2A12, the gene that encodes the glucose transporter GLUT12. In support of these findings, gene signatures of AR activity correlated withSLC2A12expression in multiple clinical cohorts. Functionally, GLUT12 was required for maximal androgen-mediated glucose uptake and cell growth in LNCaP and VCaP cells. Knockdown of GLUT12 also decreased the growth of C4-2, 22Rv1 and AR-negative PC-3 cells. This latter observation corresponded with a significant reduction in glucose uptake, indicating that additional signaling mechanisms could augment GLUT12 function in an AR-independent manner. Interestingly, GLUT12 trafficking to the plasma membrane was modulated by calcium/calmodulin-dependent protein kinase kinase 2 (CaMKK2)-5′-AMP-activated protein kinase (AMPK) signaling, a pathway we previously demonstrated to be a downstream effector of AR. Inhibition of CaMKK2-AMPK signaling decreased GLUT12 translocation to the plasma membrane by inhibiting the phosphorylation of TBC1D4, a known regulator of glucose transport. Further, AR increasedTBC1D4expression. Correspondingly, expression ofTBC1D4correlated with AR activity in prostate cancer patient samples. Taken together, these data demonstrate that prostate cancer cells can increase the functional levels of GLUT12 through multiple mechanisms to promote glucose uptake and subsequent cell growth.

The human genome is mostly transcribed, yielding a rich repository of noncoding transcripts that are involved in a myriad of biological processes including cancer. However, how many noncoding transcripts such as long noncoding RNAs (lncRNA) function in cancer is still unclear. This study identified a novel set of clinically relevant androgen-regulated lncRNAs in prostate cancer. Among this group

Human embryonic stem cell (hESC) derived mesenchymal stem cells (EMSC) are efficacious in treating a series of autoimmune, inflammatory, and degenerative diseases in animal models. However, all the EMSC derivation methods reported so far rely on two-dimensional (2D) culture systems, which are inefficient, costive and difficult for large-scale production. HESC, as an unlimited source, can be successively propagated in spheroids. Here, we demonstrate that hESC spheroids can directly differentiate into MSC spheroids (EMSCSp) within 20 days in one vessel without passaging and the system is scalable to any desired size. EMSCSp can further differentiate into osteocytes and chondrocytes in spheres or demineralized bone matrix. EMSCSp also retains immune-modulatory effects in vitro and therapeutic effects on two mouse models of colitis after dissociation. Compared to EMSC differentiated in monolayer, EMSCSp-derived cells have faster proliferation and higher yield and develop less apoptosis and slower senescence. Thus, the 3D differentiation system allows simple, cost-effective, and scalable production of high-quality EMSC and subsequently bone and cartilage tissues for therapeutic application.

Estrogen regulates diverse physiological effects and drives breast tumor progression by directly activating estrogen receptor α (ERα). However, due to the stochastic nature of gene transcription and the resulting heterogeneous cellular response, it is important to investigate estrogen-stimulated gene expression profiles at the single-cell level in order to fully understand how ERα regulates transcription in breast cancer cells. In this study, we performed single-cell transcriptome analysis on ERα-positive breast cancer cell lines following 17β-estradiol stimulation. Overall, we observed robust gene expression diversity between individual cells. Moreover, we found over two thirds of the genes in breast cancer cells displayed a bimodal expression pattern, which caused averaging artifacts and masked the identification of potential estrogen-regulated genes. We overcame this issue by reconstructing a dynamic estrogen-responsive transcriptional network from discrete time points into a pseudotemporal continuum. Pathway analysis of the differentially expressed genes derived from the pseudotemporal analysis showed an estrogen-stimulated metabolic switch that favored biosynthesis and cell proliferation but reduced estrogen degradation. In addition, we identified folate-mediated one-carbon metabolism as a novel estrogen-regulated pathway in breast cancer cells. Notably, estrogen stimulation reprogramed this pathway through the mitochondrial folate pathway to coordinately fuel polyamine and de novo purine synthesis. Finally, we showed AZIN1 and PPAT, key regulators in the above pathways, are direct ERα target genes and essential for breast cancer cell survival and growth. In summary, our single-cell study illustrated a dynamic transcriptional heterogeneity in ERα-positive breast cancer cells in response to estrogen stimulation and uncovered a novel mechanism of an estrogen-mediated metabolic switch.

The majority of the human genome is transcribed, yielding a rich repository of non-coding transcripts that are involved in a myriad of biological processes including cancer. However, how non-coding transcripts such as long non-coding RNAs (lncRNAs) function in prostate cancer is still unclear. In this study, we have identified a novel set of clinically relevant androgen-regulated lncRNAs in prostate cancer. Among this group, we showed LINC00844 is a direct androgen-regulated target that is actively transcribed in androgen receptor (AR)-dependent prostate cancer cells. The expression of LINC00844 is higher in normal prostate compared to malignant and metastatic prostate cancer samples and patients with low expression demonstrate poor prognosis and significantly increased biochemical recurrence, suggesting LINC00844 may function in suppressing tumor progression and metastasis. Indeed,in-vitroloss-of-function studies revealed that LINC00844 prevents prostate cancer cell migration and invasion. Moreover, findings from gene expression analysis indicated that LINC00844 functions intrans, affecting global androgen-regulated gene transcription. Mechanistically, we provide evidence to show LINC00844 is important in facilitating AR binding to the chromatin. Finally, we demonstrated LINC00844 mediates its phenotypic effects in part by activating the expression of NDRG1, a crucial cancer metastasis suppressor. Collectively, our findings suggest LINC00844 is a novel coregulator of AR that plays a central role in the androgen transcriptional network and the development and progression of prostate cancer.