Zhenqing Ye

Cell type-specific binding patterns reveal that TCF7L2 can be tethered to the genome by association with GATA3

BackgroundThe TCF7L2 transcription factor is linked to a variety of human diseases, including type 2 diabetes and cancer. One mechanism by which TCF7L2 could influence expression of genes involved in diverse diseases is by binding to distinct regulatory regions in different tissues. To test this hypothesis, we performed ChIP-seq for TCF7L2 in six human cell lines.ResultsWe identified 116,000 non-redundant TCF7L2 binding sites, with only 1,864 sites common to the six cell lines. Using ChIP-seq, we showed that many genomic regions that are marked by both H3K4me1 and H3K27Ac are also bound by TCF7L2, suggesting that TCF7L2 plays a critical role in enhancer activity. Bioinformatic analysis of the cell type-specific TCF7L2 binding sites revealed enrichment for multiple transcription factors, including HNF4alpha and FOXA2 motifs in HepG2 cells and the GATA3 motif in MCF7 cells. ChIP-seq analysis revealed that TCF7L2 co-localizes with HNF4alpha and FOXA2 in HepG2 cells and with GATA3 in MCF7 cells. Interestingly, in MCF7 cells the TCF7L2 motif is enriched in most TCF7L2 sites but is not enriched in the sites bound by both GATA3 and TCF7L2. This analysis suggested that GATA3 might tether TCF7L2 to the genome at these sites. To test this hypothesis, we depleted GATA3 in MCF7 cells and showed that TCF7L2 binding was lost at a subset of sites. RNA-seq analysis suggested that TCF7L2 represses transcription when tethered to the genome via GATA3.ConclusionsOur studies demonstrate a novel relationship between GATA3 and TCF7L2, and reveal important insights into TCF7L2-mediated gene regulation.

MicroRNA-31 predicts the presence of lymph node metastases and survival in patients with lung adenocarcinoma

Purpose: We conducted genome-wide miRNA-sequencing (miRNA-seq) in primary cancer tissue from patients of lung adenocarcinoma to identify markers for the presence of lymph node metastasis. Experimental Design: Markers for lymph node metastasis identified by sequencing were validated in a separate cohort using quantitative PCR. After additional validation in the The Cancer Genome Atlas (TCGA) dataset, functional characterization studies were conducted in vitro. Results: MiR-31 was upregulated in lung adenocarcinoma tissues from patients with lymph node metastases compared with those without lymph node metastases. We confirmed miR-31 to be upregulated in lymph node-positive patients in a separate patient cohort (P = 0.009, t test), and to be expressed at higher levels in adenocarcinoma tissue than in matched normal adjacent lung tissues (P < 0.0001, paired t test). MiR-31 was then validated as a marker for lymph node metastasis in an external validation cohort of 233 lung adenocarcinoma cases of the TCGA (P = 0.031, t test). In vitro functional assays showed that miR-31 increases cell migration, invasion, and proliferation in an ERK1/2 signaling-dependent manner. Notably, miR-31 was a significant predictor of survival in a multivariate cox regression model even when controlling for cancer staging. Exploratory in silico analysis showed that low expression of miR-31 is associated with excellent survival for T2N0 patients. Conclusions: We applied miRNA-seq to study microRNomes in lung adenocarcinoma tissue samples for the first time and potentially identified a miRNA predicting the presence of lymph node metastasis and survival outcomes in patients of lung adenocarcinoma. Clin Cancer Res; 19(19); 5423–33. ©2013 AACR.

MicroRNA-224 promotes tumor progression in nonsmall cell lung cancer

Significance Aberrant microRNA (miRNA) expression is involved in tumorigenesis, and miR-224 was observed to be up-regulated in certain tumor types. However, the role of miR-224 in the pathogenesis of lung cancer remains poorly understood. Here, we comprehensively analyzed and revealed mechanisms of miR-224 up-regulation and its oncogenic role in nonsmall cell lung cancer (NSCLC). We showed that miR-224 promotes cellular migratory, invasive, and proliferative capacity and tumor growth both in vitro and in vivo. Furthermore, we identified TNFα-induced protein 1 and SMAD4 as targets of miR-224. In addition, up-regulated miR-224 expression in NSCLC is partially controlled by its promoter region’s hypomethylation and activated ERK signaling. Our finding suggests that targeting miR-224 might be a promising therapeutic strategy in the treatment of NSCLC. Lung cancer is the leading cause of cancer-related deaths worldwide. Despite advancements and improvements in surgical and medical treatments, the survival rate of lung cancer patients remains frustratingly poor. Local control for early-stage nonsmall cell lung cancer (NSCLC) has dramatically improved over the last decades for both operable and inoperable patients. However, the molecular mechanisms of NSCLC invasion leading to regional and distant disease spread remain poorly understood. Here, we identify microRNA-224 (miR-224) to be significantly up-regulated in NSCLC tissues, particularly in resected NSCLC metastasis. Increased miR-224 expression promotes cell migration, invasion, and proliferation by directly targeting the tumor suppressors TNFα-induced protein 1 (TNFAIP1) and SMAD4. In concordance with in vitro studies, mouse xenograft studies validated that miR-224 functions as a potent oncogenic miRNA in NSCLC in vivo. Moreover, we found promoter hypomethylation and activated ERK signaling to be involved in the regulation of miR-224 expression in NSCLC. Up-regulated miR-224, thus, facilitates tumor progression by shifting the equilibrium of the partially antagonist functions of SMAD4 and TNFAIP1 toward enhanced invasion and growth in NSCLC. Our findings indicate that targeting miR-224 could be effective in the treatment of certain lung cancer patients.

Three-tiered role of the pioneer factor GATA2 in promoting androgen-dependent gene expression in prostate cancer

In prostate cancer, androgen receptor (AR) binding and androgen-responsive gene expression are defined by hormone-independent binding patterns of the pioneer factors FoxA1 and GATA2. Insufficient evidence of the mechanisms by which GATA2 contributes to this process precludes complete understanding of a key determinant of tissue-specific AR activity. Our observations suggest that GATA2 facilitates androgen-responsive gene expression by three distinct modes of action. By occupying novel binding sites within the AR gene locus, GATA2 positively regulates AR expression before and after androgen stimulation. Additionally, GATA2 engages AR target gene enhancers prior to hormone stimulation, producing an active and accessible chromatin environment via recruitment of the histone acetyltransferase p300. Finally, GATA2 functions in establishing and/or sustaining basal locus looping by recruiting the Mediator subunit MED1 in the absence of androgen. These mechanisms may contribute to the generally positive role of GATA2 in defining AR genome-wide binding patterns that determine androgen-responsive gene expression profiles. We also find that GATA2 and FoxA1 exhibit both independent and codependent co-occupancy of AR target gene enhancers. Identifying these determinants of AR transcriptional activity may provide a foundation for the development of future prostate cancer therapeutics that target pioneer factor function.

ZBTB33 binds unmethylated regions of the genome associated with actively expressed genes

BackgroundDNA methylation and repressive histone modifications cooperate to silence promoters. One mechanism by which regions of methylated DNA could acquire repressive histone modifications is via methyl DNA-binding transcription factors. The zinc finger protein ZBTB33 (also known as Kaiso) has been shown in vitro to bind preferentially to methylated DNA and to interact with the SMRT/NCoR histone deacetylase complexes. We have performed bioinformatic analyses of Kaiso ChIP-seq and DNA methylation datasets to test a model whereby binding of Kaiso to methylated CpGs leads to loss of acetylated histones at target promoters.ResultsOur results suggest that, contrary to expectations, Kaiso does not bind to methylated DNA in vivo but instead binds to highly active promoters that are marked with high levels of acetylated histones. In addition, our studies suggest that DNA methylation and nucleosome occupancy patterns restrict access of Kaiso to potential binding sites and influence cell type-specific binding.ConclusionsWe propose a new model for the genome-wide binding and function of Kaiso whereby Kaiso binds to unmethylated regulatory regions and contributes to the active state of target promoters.

Agonist and antagonist switch DNA motifs recognized by human androgen receptor in prostate cancer

Human transcription factors recognize specific DNA sequence motifs to regulate transcription. It is unknown whether a single transcription factor is able to bind to distinctly different motifs on chromatin, and if so, what determines the usage of specific motifs. By using a motif‐resolution chromatin immunoprecipitation‐exonuclease (ChIP‐exo) approach, we find that agonist‐liganded human androgen receptor (AR) and antagonist‐liganded AR bind to two distinctly different motifs, leading to distinct transcriptional outcomes in prostate cancer cells. Further analysis on clinical prostate tissues reveals that the binding of AR to these two distinct motifs is involved in prostate carcinogenesis. Together, these results suggest that unique ligands may switch DNA motifs recognized by ligand‐dependent transcription factors in vivo. Our findings also provide a broad mechanistic foundation for understanding ligand‐specific induction of gene expression profiles.

Ligand-dependent genomic function of glucocorticoid receptor in triple-negative breast cancer.

Glucocorticoids (GCs) have been widely used as coadjuvants in the treatment of solid tumours, but GC treatment may be associated with poor pharmacotherapeutic response or prognosis. The genomic action of GC in these tumours is largely unknown. Here we find that dexamethasone (Dex, a synthetic GC)-regulated genes in triple-negative breast cancer (TNBC) cells are associated with drug resistance. Importantly, these GC-regulated genes are aberrantly expressed in TNBC patients and are associated with unfavourable clinical outcomes. Interestingly, in TNBC cells, Compound A (CpdA, a selective GR modulator) only regulates a small number of genes not involved in carcinogenesis and therapy resistance. Mechanistic studies using a ChIP-exo approach reveal that Dex- but not CpdA-liganded glucocorticoid receptor (GR) binds to a single glucocorticoid response element (GRE), which drives the expression of pro-tumorigenic genes. Our data suggest that development of safe coadjuvant therapy should consider the distinct genomic function between Dex- and CpdA-liganded GR.

Computational analysis reveals a correlation of exon-skipping events with splicing, transcription and epigenetic factors

Alternative splicing (AS), in higher eukaryotes, is one of the mechanisms of post-transcriptional regulation that generate multiple transcripts from the same gene. One particular mode of AS is the skipping event where an exon may be alternatively excluded or constitutively included in the resulting mature mRNA. Both transcript isoforms from this skipping event site, i.e. in which the exon is either included (inclusion isoform) or excluded (skipping isoform), are typically present in one cell, and maintain a subtle balance that is vital to cellular function and dynamics. However, how the prevailing conditions dictate which isoform is expressed and what biological factors might influence the regulation of this process remain areas requiring further exploration. In this study, we have developed 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 experimentally validated several skipping site predictions by RT-PCR. Furthermore, we integrated other sequencing-based genomic data to investigate the impact of splicing activities, transcription factors (TFs) and epigenetic histone modifications on splicing outcomes. Our computational analysis found that splice sites within the skipping-isoform-dominated group (SIDG) tended to exhibit weaker MaxEntScan-calculated splice site strength around middle, ‘skipping’, exons compared to those in the inclusion-isoform-dominated group (IIDG). We further showed the positional preference pattern of splicing factors, characterized by enrichment in the intronic splice sites immediately bordering middle exons. Finally, our analysis suggested that different epigenetic factors may introduce a variable obstacle in the process of exon–intron boundary establishment leading to skipping events.

Computational Methods and Correlation of Exon-skipping Events with Splicing, Transcription, and Epigenetic Factors

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.

A survey of computational methods in transcriptome-wide alternative splicing analysis.

Abstract Alternative splicing is widely recognized for its roles in regulating genes and creating gene diversity. Consequently the identification and quantification of differentially spliced transcripts is pivotal for transcriptome analysis. Here, we review the currently available computational approaches for the analysis of RNA-sequencing data with a focus on exon-skipping events of alternative splicing and discuss the novelties as well as challenges faced to perform differential splicing analyses. In accordance with operational needs we have classified the software tools, which may be instrumental for a specific analysis based on the experimental objectives and expected outcomes. In addition, we also propose a framework for future directions by pinpointing more extensive experimental validation to assess the accuracy of the software predictions and improvements that would facilitate visualizations, data processing, and downstream analyses along with their associated software implementations.