THE p63 ISOFORM ∆Np63α INHIBITS EPITHELIAL – MESENCHYMAL TRANSITION BY PROMOTING THE EXPRESSION OF MIR-205 IN HUMAN BLADDER CANCER CELLS

p63, a p53 family member, is a transcription factor that has complex roles in cancer. This study focuses on the role of the ∆Np63α isoform in bladder cancer (BC). Epithelial – mesenchymal transition (EMT) is a physiological process that plays an important part in metastasis and drug resistance. At the molecular level, EMT is characterized by the loss of the epithelial marker E-cadherin, and the acquisition of the transcriptional repressors of E-cadherin (ZEB1, ZEB2, TWIST, SNAI1 and SNAI2). Recent publications highlight the role of microRNAs belonging to the miR-200 family and miR-205 in preventing EMT through suppression of ZEB1 and ZEB2. p53, the homologue of p63, is implicated in regulating EMT by modulating the expression of miR-200c; however, the mechanisms underlying miR-205 control remain unclear. Here we show that ∆Np63α regulates the transcription of miR-205 and controls EMT in human BC cells. We observed a strong correlation between the expression of ∆Np63α, miR-205 and E-cadherin in a panel of BC cell lines (n=28) and also in bladder primary tumors from a cohort of patients (n=98). A remarkably inverse correlation is observed between ∆Np63α and ZEB1/2 in cell lines. Stable knockdown (KD) ∆Np63α in UC6, an “epithelial” BC cell line, decreased the expression of miR-205 and induced ZEB1/2 expression, the effects that were reversed by expression of exogenous miR-205. Moreover, overexpressing ∆Np63α in UC3, a “messenchymal” BC cell line, brought about opposite results, an increase in miR-205 expression and a reduction in ZEB1/2 expression. Modulation of ∆Np63α expression resulted in a parallel change in the expression of miR-205 and miR-205 “host” gene (miR-205HG). Nuclear run-on and chromatin immunoprecipitation experiments demonstrated that ∆Np63α regulates the transcription of miR-205 through controlling the recruitment of RNA Polymerase II to the promoter of miR-205HG. Interestingly, high miR-205 expression correlated with poor clinical outcome in BC patients, consistent with our recent publication highlighting the enrichment of ∆Np63 in a lethal subset of muscle invasive BC. In summary, our data present the important roles of ∆Np63α in preventing EMT mediated by miR-205. Our study also identifies miR-205 as a potential molecular marker to predict clinical outcome in BC patients.

GENOME-WIDE PROFILING UNVEILS CRITICIAL FUNCTIONS OF p53 IN HUMAN EMBRYONIC STEM CELLS

Embryonic stem cells (ESCs) possess two unique characteristics: infinite self-renewal and the potential to differentiate into almost every cell type (pluripotency). Recently, global expression analyses of metastatic breast and lung cancers revealed an ESC-like expression program or signature, specifically for cancers that are mutant for p53 function. Surprisingly, although p53 is widely recognized as the guardian of the genome, due to its roles in cell cycle checkpoints, programmed cell death or senescence, relatively little is known about p53 functions in normal cells, especially in ESCs. My hypothesis is that p53 has specific transcription regulatory functions in human ESCs (hESCs) that a) oppose pluripotency and b) protect the stem cell genome in response to DNA damage and stress signaling. In mouse ESCs, these roles are believed to coincide, as p53 promotes differentiation in response to DNA damage, but this is unexplored in hESCs. To determine the biological roles of p53, specifically in hESCs, we mapped genome-wide chromatin interactions of p53 by chromatin immunoprecipitation and massively parallel tag sequencing (ChIP-Seq), and did so under three VIdifferent conditions of hESC status: pluripotency, differentiation-initiated and DNA-damage-induced. ChIP-Seq showed that p53 is enriched at distinct, induction-specific gene loci during each of these different conditions. Microarray gene expression analysis and functional annotation of the distinct p53-target genes revealed that p53 regulates specific genes encoding developmental regulators, which are expressed in differentiation-initiated but not DNA- damaged hESCs. We further discovered that, in response to differentiation signaling, p53 binds regions of chromatin that are repressed but also poised for rapid activation by core pluripotency factors OCT4 and NANOG in pluripotent hESCs. In response to DNA damage, genes associated with migration and motility are targeted by p53; whereas, the prime targets of p53 in control of cell death are conserved for p53 regulation in both differentiation and DNA damage. Our genome-wide profiling and bioinformatics analyses show that p53 occupies a special set of developmental regulatory genes during early differentiation of hESCs and functions in an induction-specific manner. In conclusion, our research unveiled previously unknown functions of p53 in ESC biology, which augments our understanding of one of the most deregulated proteins in human cancers.

Large-scale identification of gibberellin-related transcription factor defines group VII ERFs as functional of DELLA partners

DELLA proteins are the master negative regulators in gibberellin (GA) signaling acting in the nucleus as transcriptional regulators. The current view of DELLA action indicates that their activity relies on the physical interaction with transcription factors (TFs). Therefore, the identification of TFs through which DELLAs regulate GA responses is key to understanding these responses from a mechanistic point of view. Here, we have determined the TF interactome of the Arabidopsis (Arabidopsis thaliana) DELLA protein GIBBERELLIN INSENSITIVE and screened a collection of conditional TF overexpressors in search of those that alter GA sensitivity. As a result, we have found RELATED TO APETALA2.3, an ethylene-induced TF belonging to the group VII ETHYLENE RESPONSE FACTOR of the APETALA2/ethylene responsive element binding protein superfamily, as a DELLA interactor with physiological relevance in the context of apical hook development. The combination of transactivation assays and chromatin immunoprecipitation indicates that the interaction with GIBBERELLIN INSENSITIVE impairs the activity of RELATED TO APETALA2.3 on the target promoters. This mechanism represents a unique node in the cross regulation between the GA and ethylene signaling pathways controlling differential growth during apical hook development.

HIV-1 Tat transactivator recruits p300 and CREB-binding protein histone acetyltransferases to the viral promoter

In cells infected with HIV type 1 (HIV-1), the integrated viral promoter is present in a chromatin-bound conformation and is transcriptionally silent in the absence of stimulation. The HIV-1 Tat protein binds to a stem-loop structure at the 5′ end of viral mRNA and relieves this inhibition by inducing a remodeling of the nucleosome arrangement downstream of the transcription-initiation site. Here we show that Tat performs this activity by recruiting to the viral long terminal repeat (LTR) the transcriptional coactivator p300 and the closely related CREB-binding protein (CBP), having histone acetyltransferase (HAT) activity. Tat associates with HAT activity in human nuclear extracts and binds to p300 and CBP both in vitro and in vivo. Integrity of the basic domain of Tat is essential for this interaction. By a quantitative chromatin immunoprecipitation assay we show that the delivery of recombinant Tat induces the association of p300 and CBP with the chromosomally integrated LTR promoter. Expression of human p300 in both human and rodent cells increases the levels of Tat transactivation of the integrated LTR. These results reinforce the evidence that p300 and CBP have a pivotal function at both cellular and viral promoters and demonstrate that they also can be recruited by an RNA-targeted activator. Additionally, these findings have important implications for the understanding of the mechanisms of HIV-1 latency and reactivation.

Switch from Myc/Max to Mad1/Max binding and decrease in histone acetylation at the telomerase reverse transcriptase promoter during differentiation of HL60 cells

Recent evidence suggests that the Myc and Mad1 proteins are implicated in the regulation of the gene encoding the human telomerase reverse transcriptase (hTERT), the catalytic subunit of telomerase. We have analyzed the in vivo interaction between endogenous c-Myc and Mad1 proteins and the hTERT promoter in HL60 cells with the use of the chromatin immunoprecipitation assay. The E-boxes at the hTERT proximal promoter were occupied in vivo by c-Myc in exponentially proliferating HL60 cells but not in cells induced to differentiate by DMSO. In contrast, Mad1 protein was induced and bound to the hTERT promoter in differentiated HL60 cells. Concomitantly, the acetylation of the histones at the promoter was significantly reduced. These data suggest that the reciprocal E-box occupancy by c-Myc and Mad1 is responsible for activation and repression of the hTERT gene in proliferating and differentiated HL60 cells, respectively. Furthermore, the histone deacetylase inhibitor trichostatin A inhibited deacetylation of histones at the hTERT promoter and attenuated the repression of hTERT transcription during HL60 cell differentiation. In addition, trichostatin A treatment activated hTERT transcription in resting human lymphocytes and fibroblasts. Taken together, these results indicate that acetylation/deacetylation of histones is operative in the regulation of hTERT expression.

Selective association of the methyl-CpG binding protein MBD2 with the silent p14/p16 locus in human neoplasia

DNA methylation of tumor suppressor genes is a common feature of human cancer. The cyclin-dependent kinase inhibitor gene p16/Ink4A is hypermethylated in a wide range of malignant tissues and the p14/ARF gene located 20 kb upstream on chromosome 9p21 is also methylated in carcinomas. p14/ARF (ARF, alternative reading frame) does not inhibit the activities of cyclins or cyclin-dependent kinase complexes; however, the importance of the two gene products in the etiology of cancer resides in their involvement in two major cell cycle regulatory pathways: p53 and the retinoblastoma protein, Rb, respectively. Distinct first exons driven from separate promoters are spliced onto the common exons 2 and 3 and the resulting proteins are translated in different reading frames. Both genes are expressed in normal cells but can be alternatively or coordinately silenced when their CpG islands are hypermethylated. Herein, we examined the presence of methyl-CpG binding proteins associated with aberrantly methylated promoters, the distribution of acetylated histones H3 and H4 by chromatin immunoprecipitation assays, and the effect of chemical treatment with 5-aza-2′-deoxycytidine (5aza-dC) and trichostatin A on gene induction in colon cell lines by quantitative reverse transcriptase–PCR. We observed that the methyl-CpG binding protein MBD2 is targeted to methylated regulatory regions and excludes the acetylated histones H3 and H4, resulting in a localized inactive chromatin configuration. When methylated, the genes can be induced by 5aza-dC but the combined action of 5aza-dC and trichostatin A results in robust gene expression. Thus, methyl-CpG binding proteins and histone deacetylases appear to cooperate in vivo, with a dominant effect of DNA methylation toward histone acetylation, and repress expression of tumor suppressor genes hypermethylated in cancers.

Epigenetic engineering shows H3K4me2 is required for HJURP targeting and CENP-A assembly on a synthetic human kinetochore

Kinetochores assemble on distinct 'centrochromatin' containing the histone H3 variant CENP-A and interspersed nucleosomes dimethylated on H3K4 (H3K4me2). Little is known about how the chromatin environment at active centromeres governs centromeric structure and function. Here, we report that centrochromatin resembles K4-K36 domains found in the body of some actively transcribed housekeeping genes. By tethering the lysine-specific demethylase 1 (LSD1), we specifically depleted H3K4me2, a modification thought to have a role in transcriptional memory, from the kinetochore of a synthetic human artificial chromosome (HAC). H3K4me2 depletion caused kinetochores to suffer a rapid loss of transcription of the underlying α-satellite DNA and to no longer efficiently recruit HJURP, the CENP-A chaperone. Kinetochores depleted of H3K4me2 remained functional in the short term, but were defective in incorporation of CENP-A, and were gradually inactivated. Our data provide a functional link between the centromeric chromatin, α-satellite transcription, maintenance of CENP-A levels and kinetochore stability.

Transient Activation of Meox1 Is an Early Component of the Gene Regulatory Network Downstream of Hoxa2

Hox genes encode transcription factors that regulate morphogenesis in all animals with bilateral symmetry. Although Hox genes have been extensively studied, their molecular function is not clear in vertebrates, and only a limited number of genes regulated by Hox transcription factors have been identified. Hoxa2 is required for correct development of the second branchial arch, its major domain of expression. We now show that Meox1 is genetically downstream from Hoxa2 and is a direct target. Meox1 expression is downregulated in the second arch of Hoxa2 mouse mutant embryos. In chromatin immunoprecipitation (ChIP), Hoxa2 binds to the Meox1 proximal promoter. Two highly conserved binding sites contained in this sequence are required for Hoxa2-dependent activation of the Meox1 promoter. Remarkably, in the absence of Meox1 and its close homolog Meox2, the second branchial arch develops abnormally and two of the three skeletal elements patterned by Hoxa2 are malformed. Finally, we show that Meox1 can specifically bind the DNA sequences recognized by Hoxa2 on its functional target genes. These results provide new insight into the Hoxa2 regulatory network that controls branchial arch identity.

Menin associates with a trithorax family histone methyltransferase complex and with the Hoxc8 locus

The cellular function of the menin tumor suppressor protein, product of the MEN1 gene mutated in familial multiple endocrine neoplasia type 1, has not been defined. We now show that menin is associated with a histone methyltransferase complex containing two trithorax family proteins, MLL2 and Ash2L, and other homologs of the yeast Set1 assembly. This menin-associated complex methylates histone H3 on lysine 4. A subset of tumor-derived menin mutants lacks the associated histone methyltransferase activity. In addition, menin is associated with RNA polymerase II whose large subunit carboxyl-terminal domain is phosphorylated on Ser5. Men1 knockout embryos and cells show decreased expression of the homeobox genes Hoxc6 and Hoxc8. Chromatin immunoprecipitation experiments reveal that menin is bound to the Hoxc8 locus. These results suggest that menin activates the transcription of differentiation-regulating genes by covalent histone modification, and that this activity is related to tumor suppression by MEN1.

Knowledge integration for analyzing ChIP-seq

To capture the genomic profiles for histone modification, chromatin immunoprecipitation (ChIP) is combined with next generation sequencing, which is called ChIP-seq. However, enriched regions generated from the ChIP-seq data are only evaluated on the limited knowledge acquired from manually examining the relevant biological literature. This paper proposes a novel framework, which integrates multiple knowledge sources such as biological literature, Gene Ontology, and microarray data. In order to precisely analyze ChIP-seq data for histone modification, knowledge integration is based on a unified probabilistic model. The model is employed to re-rank the enriched regions generated from peak finding algorithms. Through filtering the reranked enriched regions using some predefined threshold, more reliable and precise results could be generated. The combination of the multiple knowledge sources with the peaking finding algorithm produces a new paradigm for ChIP-seq data analysis. © (2012) Trans Tech Publications, Switzerland.

REST is a hypoxia-responsive transcriptional repressor

Cellular exposure to hypoxia results in altered gene expression in a range of physiologic and pathophysiologic states. Discrete cohorts of genes can be either up- or down-regulated in response to hypoxia. While the Hypoxia-Inducible Factor (HIF) is the primary driver of hypoxia-induced adaptive gene expression, less is known about the signalling mechanisms regulating hypoxiadependent gene repression. Using RNA-seq, we demonstrate that equivalent numbers of genes are induced and repressed in human embryonic kidney (HEK293) cells. We demonstrate that nuclear localization of the Repressor Element 1-Silencing Transcription factor (REST) is induced in hypoxia and that REST is responsible for regulating approximately 20% of the hypoxia-repressed genes. Using chromatin immunoprecipitation assays we demonstrate that REST-dependent gene repression is at least in part mediated by direct binding to the promoters of target genes. Based on these data, we propose that REST is a key mediator of gene repression in hypoxia.

Understanding Ten-Eleven Translocation-2 in Hematological and Nervous Systems

I proposed the study of two distinct aspects of Ten-Eleven Translocation 2 (TET2) protein for understanding specific functions in different body systems. In Part I, I characterized the molecular mechanisms of Tet2 in the hematological system. As the second member of Ten-Eleven Translocation protein family, TET2 is frequently mutated in leukemic patients. Previous studies have shown that the TET2 mutations frequently occur in 20% myelodysplastic syndrome/myeloproliferative neoplasm (MDS/MPN), 10% T-cell lymphoma leukemia and 2% B-cell lymphoma leukemia. Genetic mouse models also display distinct phenotypes of various types of hematological malignancies. I performed 5-hydroxymethylcytosine (5hmC) chromatin immunoprecipitation sequencing (ChIP-Seq) and RNA sequencing (RNA-Seq) of hematopoietic stem/progenitor cells to determine whether the deletion of Tet2 can affect the abundance of 5hmC at myeloid, T-cell and B-cell specific gene transcription start sites, which ultimately result in various hematological malignancies. Subsequent Exome sequencing (Exome-Seq) showed that disease-specific genes are mutated in different types of tumors, which suggests that TET2 may protect the genome from being mutated. The direct interaction between TET2 and Mutator S Homolog 6 (MSH6) protein suggests TET2 is involved in DNA mismatch repair. Finally, in vivo mismatch repair studies show that the loss of Tet2 causes a mutator phenotype. Taken together, my data indicate that TET2 binds to MSH6 to protect genome integrity. In Part II, I intended to better understand the role of Tet2 in the nervous system. 5-hydroxymethylcytosine regulates epigenetic modification during neurodevelopment and aging. Thus, Tet2 may play a critical role in regulating adult neurogenesis. To examine the physiological significance of Tet2 in the nervous system, I first showed that the deletion of Tet2 reduces the 5hmC levels in neural stem cells. Mice lacking Tet2 show abnormal hippocampal neurogenesis along with 5hmC alternations at different gene promoters and corresponding gene expression downregulation. Through the luciferase reporter assay, two neural factors Neurogenic differentiation 1 (NeuroD1) and Glial fibrillary acidic protein (Gfap) were down-regulated in Tet2 knockout cells. My results suggest that Tet2 regulates neural stem/progenitor cell proliferation and differentiation in adult brain.

Characterizing Genetic Drivers of Lymphoma through High-Throughput Sequencing

The advent of next-generation sequencing, now nearing a decade in age, has enabled, among other capabilities, measurement of genome-wide sequence features at unprecedented scale and resolution.

In this dissertation, I describe work to understand the genetic underpinnings of non-Hodgkin’s lymphoma through exploration of the epigenetics of its cell of origin, initial characterization and interpretation of driver mutations, and finally, a larger-scale, population-level study that incorporates mutation interpretation with clinical outcome.

In the first research chapter, I describe genomic characteristics of lymphomas through the lens of their cells of origin. Just as many other cancers, such as breast cancer or lung cancer, are categorized based on their cell of origin, lymphoma subtypes can be examined through the context of their normal B Cells of origin, Naïve, Germinal Center, and post-Germinal Center. By applying integrative analysis of the epigenetics of normal B Cells of origin through chromatin-immunoprecipitation sequencing, we find that differences in normal B Cell subtypes are reflected in the mutational landscapes of the cancers that arise from them, namely Mantle Cell, Burkitt, and Diffuse Large B-Cell Lymphoma.

In the next research chapter, I describe our first endeavor into understanding the genetic heterogeneity of Diffuse Large B Cell Lymphoma, the most common form of non-Hodgkin’s lymphoma, which affects 100,000 patients in the world. Through whole-genome sequencing of 1 case as well as whole-exome sequencing of 94 cases, we characterize the most recurrent genetic features of DLBCL and lay the groundwork for a larger study.

In the last research chapter, I describe work to characterize and interpret the whole exomes of 1001 cases of DLBCL in the largest single-cancer study to date. This highly-powered study enabled sub-gene, gene-level, and gene-network level understanding of driver mutations within DLBCL. Moreover, matched genomic and clinical data enabled the connection of these driver mutations to clinical features such as treatment response or overall survival. As sequencing costs continue to drop, whole-exome sequencing will become a routine clinical assay, and another diagnostic dimension in addition to existing methods such as histology. However, to unlock the full utility of sequencing data, we must be able to interpret it. This study undertakes a first step in developing the understanding necessary to uncover the genomic signals of DLBCL hidden within its exomes. However, beyond the scope of this one disease, the experimental and analytical methods can be readily applied to other cancer sequencing studies.

Thus, this dissertation leverages next-generation sequencing analysis to understand the genetic underpinnings of lymphoma, both by examining its normal cells of origin as well as through a large-scale study to sensitively identify recurrently mutated genes and their relationship to clinical outcome.

Mapping Protein-DNA Interactions Using ChIP-exo and Illumina-Based Sequencing

Chromatin immunoprecipitation (ChIP) provides a means of enriching DNA associated with transcription factors, histone modifications, and indeed any other proteins for which suitably characterized antibodies are available. Over the years, sequence detection has progressed from quantitative real-time PCR and Southern blotting to microarrays (ChIP-chip) and now high-throughput sequencing (ChIP-seq). This progression has vastly increased the sequence coverage and data volumes generated. This in turn has enabled informaticians to predict the identity of multi-protein complexes on DNA based on the overrepresentation of sequence motifs in DNA enriched by ChIP with a single antibody against a single protein. In the course of the development of high-throughput sequencing, little has changed in the ChIP methodology until recently. In the last three years, a number of modifications have been made to the ChIP protocol with the goal of enhancing the sensitivity of the method and further reducing the levels of nonspecific background sequences in ChIPped samples. In this chapter, we provide a brief commentary on these methodological changes and describe a detailed ChIP-exo method able to generate narrower peaks and greater peak coverage from ChIPped material.

FACS-based purification of Arabidopsis microspores, sperm cells and vegetative nuclei

Background: The male germline in flowering plants differentiates by asymmetric division of haploid uninucleated microspores, giving rise to a vegetative cell enclosing a smaller generative cell, which eventually undergoes a second mitosis to originate two sperm cells. The vegetative cell and the sperm cells activate distinct genetic and epigenetic mechanisms to control pollen tube growth and germ cell specification, respectively. Therefore, a comprehensive characterization of these processes relies on efficient methods to isolate each of the different cell types throughout male gametogenesis. Results: We developed stable transgenic Arabidopsis lines and reliable purification tools based on Fluorescence-Activated Cell Sorting (FACS) in order to isolate highly pure and viable fractions of each cell/nuclei type before and after pollen mitosis. In the case of mature pollen, this was accomplished by expressing GFP and RFP in the sperm and vegetative nuclei, respectively, resulting in 99% pure sorted populations. Microspores were also purified by FACS taking advantage of their characteristic small size and autofluorescent properties, and were confirmed to be 98% pure. Conclusions: We provide simple and efficient FACS-based purification protocols for Arabidopsis microspores, vegetative nuclei and sperm cells. This paves the way for subsequent molecular analysis such as transcriptomics, DNA methylation analysis and chromatin immunoprecipitation, in the developmental context of microgametogenesis in Arabidopsis.