Person:

Liu, Xiaole

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 ((R^2 = 0.6213, P < 5 × 10^{−11})) and KLK2 ((R^2 = 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.

The CRISPR/Cas9 system has revolutionized mammalian somatic cell genetics. Genome-wide functional screens using CRISPR/Cas9-mediated knockout or dCas9 fusion-mediated inhibition/activation (CRISPRi/a) are powerful techniques for discovering phenotype-associated gene function. We systematically assessed the DNA sequence features that contribute to single guide RNA (sgRNA) efficiency in CRISPR-based screens. Leveraging the information from multiple designs, we derived a new sequence model for predicting sgRNA efficiency in CRISPR/Cas9 knockout experiments. Our model confirmed known features and suggested new features including a preference for cytosine at the cleavage site. The model was experimentally validated for sgRNA-mediated mutation rate and protein knockout efficiency. Tested on independent data sets, the model achieved significant results in both positive and negative selection conditions and outperformed existing models. We also found that the sequence preference for CRISPRi/a is substantially different from that for CRISPR/Cas9 knockout and propose a new model for predicting sgRNA efficiency in CRISPRi/a experiments. These results facilitate the genome-wide design of improved sgRNA for both knockout and CRISPRi/a studies.

Cell differentiation requires remodeling of tissue-specific gene loci and activities of key transcriptional regulators, which are recognized for their dominant control over cellular programs. Using epigenomic methods, we characterized enhancer elements specifically modified in differentiating intestinal epithelial cells and found enrichment of transcription factor-binding motifs corresponding to CDX2, a critical regulator of the intestine. Directed investigation revealed surprising lability in CDX2 occupancy of the genome, with redistribution from hundreds of sites occupied only in proliferating cells to thousands of new sites in differentiated cells. Knockout mice confirmed distinct Cdx2 requirements in dividing and mature adult intestinal cells, including responsibility for the active enhancer configuration associated with maturity. Dynamic CDX2 occupancy corresponds with condition-specific gene expression and, importantly, to differential co-occupancy with other tissue-restricted transcription factors such as GATA6 and HNF4A. These results reveal dynamic, context-specific functions and mechanisms of a prominent transcriptional regulator within a cell lineage.

Selective activity of a specific set of enhancers defines tissue-specific gene transcription. The pioneer factor FOXA1 has been shown to induce functional enhancer competency through chromatin openings. We have previously found that FOXA1 is recruited to thousands of regions across the genome of a given cell type. Here, we monitored the chromatin structure at FOXA1 binding sites on a chromosome-wide scale using formaldehyde assisted isolation of regulatory elements (FAIRE). Surprisingly, we find that a significant fraction of FOXA1-bound sites have a relatively closed chromatin conformation linked to a shift of the epigenetic signature toward repressive histone marks. Importantly, these sites are not correlated with gene expression in a given cell type suggesting that FOXA1 is required, but not sufficient, for the functional activity of bound enhancers. Interestingly, we find that a significant proportion of the inactive FOXA1-bound regulatory sites in one cell type are actually functional in another cellular context. We found that at least half of the FOXA1 binding sites from a given cell type are shared with another cell lineage. Mechanisms that restrict the activity of shared FOXA1-bound enhancers likely play a significant role in defining the cell-type-specific functions of FOXA1.

Surprisingly few pathways signal between cells, raising questions about mechanisms for tissue-specific responses. In particular, Wnt ligands signal in many mammalian tissues, including the intestinal epithelium, where constitutive signaling causes cancer. Genome-wide analysis of DNA cis-regulatory regions bound by the intestine-restricted transcription factor CDX2 in colonic cells uncovered highly significant overrepresentation of sequences that bind TCF4, a transcriptional effector of intestinal Wnt signaling. Chromatin immunoprecipitation confirmed TCF4 occupancy at most such sites and co-occupancy of CDX2 and TCF4 across short distances. A region spanning the single nucleotide polymorphism rs6983267, which lies within a MYC enhancer and confers colorectal cancer risk in humans, represented one of many co-occupied sites. Co-occupancy correlated with intestine-specific gene expression and CDX2 loss reduced TCF4 binding. These results implicate CDX2 in directing TCF4 binding in intestinal cells. Co-occupancy of regulatory regions by signal-effector and tissue-restricted transcription factors may represent a general mechanism for ubiquitous signaling pathways to achieve tissue-specific outcomes.

Genome-wide screening using CRISPR coupled with nuclease Cas9 (CRISPR/Cas9) is a powerful technology for the systematic evaluation of gene function. Statistically principled analysis is needed for the accurate identification of gene hits and associated pathways. Here, we describe how to perform computational analysis of CRISPR screens using the MAGeCKFlute pipeline. MAGeCKFlute combines the MAGeCK and MAGeCK-VISPR algorithms and incorporates additional downstream analysis functionalities. MAGeCKFlute is distinguished from other currently available tools by being a comprehensive pipeline that contains a series of functions for analyzing CRISPR screen data. This protocol explains how to use MAGeCKFlute to perform quality control, normalization, batch effect removal, copy number bias correction, gene hit identification, and downstream functional enrichment analysis for CRISPR screens. We also describe gene identification and data analysis in CRISPR screens involving drug treatment. Completing the entire MAGeCKFlute pipeline requires approximately two hours on a desktop computer running Linux or Mac OS and with R support. The MAGeCKFlute package is available at http://www.bioconductor.org/packages/release/bioc/html/MAGeCKFlute.html.