Yiwen Chen

ChIP-seq guidelines and practices of the ENCODE and modENCODE consortia

Chromatin immunoprecipitation (ChIP) followed by high-throughput DNA sequencing (ChIP-seq) has become a valuable and widely used approach for mapping the genomic location of transcription-factor binding and histone modifications in living cells. Despite its widespread use, there are considerable differences in how these experiments are conducted, how the results are scored and evaluated for quality, and how the data and metadata are archived for public use. These practices affect the quality and utility of any global ChIP experiment. Through our experience in performing ChIP-seq experiments, the ENCODE and modENCODE consortia have developed a set of working standards and guidelines for ChIP experiments that are updated routinely. The current guidelines address antibody validation, experimental replication, sequencing depth, data and metadata reporting, and data quality assessment. We discuss how ChIP quality, assessed in these ways, affects different uses of ChIP-seq data. All data sets used in the analysis have been deposited for public viewing and downloading at the ENCODE (http://encodeproject.org/ENCODE/) and modENCODE (http://www.modencode.org/) portals.

Cistrome: an integrative platform for transcriptional regulation studies

The increasing volume of ChIP-chip and ChIP-seq data being generated creates a challenge for standard, integrative and reproducible bioinformatics data analysis platforms. We developed a web-based application called Cistrome, based on the Galaxy open source framework. In addition to the standard Galaxy functions, Cistrome has 29 ChIP-chip- and ChIP-seq-specific tools in three major categories, from preliminary peak calling and correlation analyses to downstream genome feature association, gene expression analyses, and motif discovery. Cistrome is available at http://cistrome.org/ap/.

REST and stress resistance in ageing and Alzheimer’s disease

Human neurons are functional over an entire lifetime, yet the mechanisms that preserve function and protect against neurodegeneration during ageing are unknown. Here we show that induction of the repressor element 1-silencing transcription factor (REST; also known as neuron-restrictive silencer factor, NRSF) is a universal feature of normal ageing in human cortical and hippocampal neurons. REST is lost, however, in mild cognitive impairment and Alzheimer’s disease. Chromatin immunoprecipitation with deep sequencing and expression analysis show that REST represses genes that promote cell death and Alzheimer’s disease pathology, and induces the expression of stress response genes. Moreover, REST potently protects neurons from oxidative stress and amyloid β-protein toxicity, and conditional deletion of REST in the mouse brain leads to age-related neurodegeneration. A functional orthologue of REST, Caenorhabditis elegans SPR-4, also protects against oxidative stress and amyloid β-protein toxicity. During normal ageing, REST is induced in part by cell non-autonomous Wnt signalling. However, in Alzheimer’s disease, frontotemporal dementia and dementia with Lewy bodies, REST is lost from the nucleus and appears in autophagosomes together with pathological misfolded proteins. Finally, REST levels during ageing are closely correlated with cognitive preservation and longevity. Thus, the activation state of REST may distinguish neuroprotection from neurodegeneration in the ageing brain.

XBP1 promotes triple-negative breast cancer by controlling the HIF1α pathway

Cancer cells induce a set of adaptive response pathways to survive in the face of stressors due to inadequate vascularization. One such adaptive pathway is the unfolded protein (UPR) or endoplasmic reticulum (ER) stress response mediated in part by the ER-localized transmembrane sensor IRE1 (ref. 2) and its substrate XBP1 (ref. 3). Previous studies report UPR activation in various human tumours, but the role of XBP1 in cancer progression in mammary epithelial cells is largely unknown. Triple-negative breast cancer (TNBC)—a form of breast cancer in which tumour cells do not express the genes for oestrogen receptor, progesterone receptor and HER2 (also called ERBB2 or NEU)—is a highly aggressive malignancy with limited treatment options. Here we report that XBP1 is activated in TNBC and has a pivotal role in the tumorigenicity and progression of this human breast cancer subtype. In breast cancer cell line models, depletion of XBP1 inhibited tumour growth and tumour relapse and reduced the CD44highCD24low population. Hypoxia-inducing factor 1α (HIF1α) is known to be hyperactivated in TNBCs. Genome-wide mapping of the XBP1 transcriptional regulatory network revealed that XBP1 drives TNBC tumorigenicity by assembling a transcriptional complex with HIF1α that regulates the expression of HIF1α targets via the recruitment of RNA polymerase II. Analysis of independent cohorts of patients with TNBC revealed a specific XBP1 gene expression signature that was highly correlated with HIF1α and hypoxia-driven signatures and that strongly associated with poor prognosis. Our findings reveal a key function for the XBP1 branch of the UPR in TNBC and indicate that targeting this pathway may offer alternative treatment strategies for this aggressive subtype of breast cancer.

Integrative genomic analyses reveal clinically relevant long noncoding RNAs in human cancer

Despite growing appreciation of the importance of long noncoding RNAs (lncRNAs) in normal physiology and disease, our knowledge of cancer-related lncRNAs remains limited. By repurposing microarray probes, we constructed expression profiles of 10,207 lncRNA genes in approximately 1,300 tumors over four different cancer types. Through integrative analysis of the lncRNA expression profiles with clinical outcome and somatic copy-number alterations, we identified lncRNAs that are associated with cancer subtypes and clinical prognosis and predicted those that are potential drivers of cancer progression. We validated our predictions by experimentally confirming prostate cancer cell growth dependence on two newly identified lncRNAs. Our analysis provides a resource of clinically relevant lncRNAs for the development of lncRNA biomarkers and the identification of lncRNA therapeutic targets. It also demonstrates the power of integrating publically available genomic data sets and clinical information for discovering disease-associated lncRNAs.

Targeting androgen receptor in estrogen receptor-negative breast cancer.

Endocrine therapies for breast cancer that target the estrogen receptor (ER) are ineffective in the 25%-30% of cases that are ER negative (ER-). Androgen receptor (AR) is expressed in 60%-70% of breast tumors, independent of ER status. How androgens and AR regulate breast cancer growth remains largely unknown. We find that AR is enriched in ER- breast tumors that overexpress HER2. Through analysis of the AR cistrome and androgen-regulated gene expression in ER-/HER2+ breast cancers we find that AR mediates ligand-dependent activation of Wnt and HER2 signaling pathways through direct transcriptional induction of WNT7B and HER3. Specific targeting of AR, Wnt or HER2 signaling impairs androgen-stimulated tumor cell growth suggesting potential therapeutic approaches for ER-/HER2+ breast cancers.

Enhancer RNAs participate in androgen receptor-driven looping that selectively enhances gene activation

Significance We report that enhancer RNAs (eRNAs), a class of long noncoding RNAs, participate in the androgen receptor (AR)-dependent looping complex that enhances spatial communication of distal enhancers and target promoters, leading to transcriptional activation events. Furthermore, our data show that KLK3 eRNA (KLK3e) selectively enhances the gene expression of AR-regulated genes, and provide evidence for a positive regulatory loop in which AR-dependent transcription is modulated by an intermediate eRNA. These findings may translate into improved RNA-based therapy (eRNA suppression) to enhance the durability of androgen deprivation therapy (ADT) and prediction of the efficacy of ADT by measuring the enhancer-derived activity (eRNA expression) in prostate tumors. The androgen receptor (AR) is a key factor that regulates the behavior and fate of prostate cancer cells. The AR-regulated network is activated when AR binds enhancer elements and modulates specific enhancer–promoter looping. Kallikrein-related peptidase 3 (KLK3), which codes for prostate-specific antigen (PSA), is a well-known AR-regulated gene and its upstream enhancers produce bidirectional enhancer RNAs (eRNAs), termed KLK3e. Here, we demonstrate that KLK3e facilitates the spatial interaction of the KLK3 enhancer and the KLK2 promoter and enhances long-distance KLK2 transcriptional activation. KLK3e carries the core enhancer element derived from the androgen response element III (ARE III), which is required for the interaction of AR and Mediator 1 (Med1). Furthermore, we show that KLK3e processes RNA-dependent enhancer activity depending on the integrity of core enhancer elements. The transcription of KLK3e was detectable and its expression is significantly correlated with KLK3 (R2 = 0.6213, P < 5 × 10−11) and KLK2 (R2 = 0.5893, P < 5 × 10−10) in human prostate tissues. Interestingly, RNAi silencing of KLK3e resulted in a modest negative effect on prostate cancer cell proliferation. Accordingly, we report that an androgen-induced eRNA scaffolds the AR-associated protein complex that modulates chromosomal architecture and selectively enhances AR-dependent gene expression.

Systematic evaluation of factors influencing ChIP-seq fidelity

We evaluated how variations in sequencing depth and other parameters influence interpretation of chromatin immunoprecipitation–sequencing (ChIP-seq) experiments. Using Drosophila melanogaster S2 cells, we generated ChIP-seq data sets for a site-specific transcription factor (Suppressor of Hairy-wing) and a histone modification (H3K36me3). We detected a chromatin-state bias: open chromatin regions yielded higher coverage, which led to false positives if not corrected. This bias had a greater effect on detection specificity than any base-composition bias. Paired-end sequencing revealed that single-end data underestimated ChIP-library complexity at high coverage. Removal of reads originating at the same base reduced false-positives but had little effect on detection sensitivity. Even at mappable-genome coverage depth of ∼1 read per base pair, ∼1% of the narrow peaks detected on a tiling array were missed by ChIP-seq. Evaluation of widely used ChIP-seq analysis tools suggests that adjustments or algorithm improvements are required to handle data sets with deep coverage.

A Comprehensive View of Nuclear Receptor Cancer Cistromes

Nuclear receptors comprise a superfamily of ligand-activated transcription factors that play important roles in both physiology and diseases including cancer. The technologies of chromatin immunoprecipitation followed by array hybridization (ChIP-chip) or massively parallel sequencing (ChIP-seq) has been used to map, at an unprecedented rate, the in vivo genome-wide binding (cistrome) of nuclear receptors in both normal and cancer cells. We developed a curated database of 88 nuclear receptor cistrome data sets and other associated high-throughput data sets including 121 collaborating factor cistromes, 94 epigenomes, and 319 transcriptomes. Through integrative analysis of the curated nuclear receptor ChIP-chip/seq data sets, we discovered novel factor-specific noncanonical motifs that may have important regulatory roles. We also revealed a common feature of nuclear receptor pioneering factors to recognize relatively short and AT-rich motifs. Most nuclear receptors bind predominantly to introns and distal intergenetic regions, and binding sites closer to transcription start sites were found to be neither stronger nor more evolutionarily conserved. Interestingly, while most nuclear receptors appear to be predominantly transcriptional activators, our analysis suggests that the binding of ESR1, RARA, and RARG has both activating and repressive effects. Through meta-analysis of different omic data of the same cancer cell line model from multiple studies, we generated consensus cistrome and expression profiles. We further made probabilistic predictions of the nuclear receptor target genes by integrating cistrome and transcriptome data and validated the predictions using expression data from tumor samples. The final database, with comprehensive cistrome, epigenome, and transcriptome data sets and downstream analysis results, constitutes a valuable resource for the nuclear receptor and cancer community.

Integrative analyses reveal a long noncoding RNA-mediated sponge regulatory network in prostate cancer.

Mounting evidence suggests that long noncoding RNAs (lncRNAs) can function as microRNA sponges and compete for microRNA binding to protein-coding transcripts. However, the prevalence, functional significance and targets of lncRNA-mediated sponge regulation of cancer are mostly unknown. Here we identify a lncRNA-mediated sponge regulatory network that affects the expression of many protein-coding prostate cancer driver genes, by integrating analysis of sequence features and gene expression profiles of both lncRNAs and protein-coding genes in tumours. We confirm the tumour-suppressive function of two lncRNAs (TUG1 and CTB-89H12.4) and their regulation of PTEN expression in prostate cancer. Surprisingly, one of the two lncRNAs, TUG1, was previously known for its function in polycomb repressive complex 2 (PRC2)-mediated transcriptional regulation, suggesting its sub-cellular localization-dependent function. Our findings not only suggest an important role of lncRNA-mediated sponge regulation in cancer, but also underscore the critical influence of cytoplasmic localization on the efficacy of a sponge lncRNA.