ChIP-chip Data Research Articles

Global repression by tailless during segmentation

The orphan nuclear receptor Tailless (Tll) exhibits conserved roles in brain formation and maintenance that are shared, for example, with vertebrate orthologous forms (Tlx). However, the early expression of tll in two gap domains in the segmentation cascade of Drosophila is unusual even for most other insects. Here we investigate tll regulation on pair-rule stripes. With ectopic misexpression of tll we detected unexpected repression of almost all pair-rule stripes of hairy (h), even-skipped (eve), runt (run), and fushi-tarazu (ftz). Examining Tll embryonic ChIP-chip data with regions mapped as Cis-Regulatory Modules (CRMs) of pair-rule stripes we verified Tll interactions to these regions. With the ChIP-chip data we also verified Tll interactions to the CRMs of gap domains and in the misexpression assay, Tll-mediated repression on Kruppel (Kr), kni (kni) and giant (gt) according to their differential sensitivity to Tll. These results with gap genes confirmed previous data from the literature and argue against indirect repression roles of Tll in the striped pattern. Moreover, the prediction of Tll binding sites in the CRMs of eve stripes and the mathematical modeling of their removal using an experimentally validated theoretical framework shows effects on eve stripes compatible with the absence of a repressor binding to the CRMs. In addition, modeling increased tll levels in the embryo results in the differential repression of eve stripes, agreeing well with the results of the misexpression assay. In genetic assays we investigated eve 5, that is strongly repressed by the ectopic domain and representative of more central stripes not previously implied to be under direct regulation of tll. While this stripe is little affected in tll-, its posterior border is expanded in gt- but detected with even greater expansion in gt-;tll-. We end up by discussing tll with key roles in combinatorial repression mechanisms to contain the expression of medial patterns of the segmentation cascade in the extremities of the embryo.

Changes in Hox Gene Chromatin Organization during Odontogenic Lineage Specification.

Craniofacial tissues comprise highly evolved organs characterized by a relative lack of expression in the HOX family transcription factors. In the present study, we sought to define the epigenetic events that limit HOX gene expression from undifferentiated neural crest cells to semi-differentiated odontogenic progenitors and to explore the effects of elevated levels of HOX. The ChIP-chip data demonstrated high levels of repressive H3K27me3 marks on the HOX gene promoters in ES and cranial neural crest cells when compared to the H3K4me3 marks, while the K4/K27 ratio was less repressive in the odontogenic progenitors, dental follicle, dental pulp, periodontal ligament fibroblasts, alveolar bone osteoblasts, and cementoblasts. The gene expression of multiple HOX genes, especially those from the HOXA and HOXB clusters, was significantly elevated and many times higher in alveolar bone cells than in the dental follicle cells. In addition, the HOX levels in the skeletal osteoblasts were many times higher in the trunk osteoblasts compared to the alveolar bone osteoblasts, and the repressive mark H3K27me3 promoter occupancy was substantially and significantly elevated in the alveolar bone osteoblasts when compared to the trunk osteoblasts. To explore the effect of elevated HOX levels in craniofacial neural crest cells, HOX expression was induced by transfecting cells with the Cdx4 transcription factor, resulting in a significant decrease in the mineralization markers, RUNX2, OSX, and OCN upon HOX elevation. Promoting HOX gene expression in developing teeth using the small molecule EZH2 inhibitor GSK126 resulted in an increased number of patterning events, supernumerary cusp formation, and increased Hoxa4 and Hoxb6 gene expression when compared to the controls. Together, these studies illustrate the profound effects of epigenetic regulatory events at all stages of the differentiation of craniofacial peripheral tissues from the neural crest, including lineage specification, tissue differentiation, and patterning.

Abstract 777: Deciphering the role of MSI2 as a regulator of mantle cell lymphoma stem-like properties

Abstract Mantle cell lymphoma (MCL) is a mature B cell neoplasm with two distinct biological subtypes regarding clinical, molecular and pathological features. Conventional MCL (cMCL) subtype is characterized by high aggressiveness and poor outcome with frequent relapses after initial response to treatment, suggesting the presence of a population of cancer stem cells (CSC) able to self-renew and responsible for drug resistance. SOX11 embryonic transcription factor is aberrantly overexpressed mostly in cMCL and has an oncogenic role in MCL. MCL-CSC have been isolated from MCL cases by different groups. However, its relationship to SOX11 expression and contribution to MCL clinical features and aggressive behavior remains unknown. We hypothesized that SOX11 may promote stem cell-like properties through the activation of stem cell-related genes, leading to an aggressive and incurable tumor in MCL patients. Here, we integrated the differential gene expression profile (GEP) between SOX11+ and SOX11- MCL primary cases with stem cell-related genes, and specific SOX11 ChIP-chip and epigenetic data in MCL cell lines and primary cases, to identify potential mediators of CSC directly regulated by SOX11 in MCL. We observed an enrichment of hematopoietic (HSC) and leukemic stem cells (LSC) gene signatures in SOX11+ compared to SOX11- MCL primary cases. MSI2, an RNA binding protein that maintains self-renewal and prevents differentiation in HSC, emerged as one of the most significant upregulated CSC-related gene in SOX11+ compared to SOX11- MCLs, positively correlating with SOX11 expression in MCL primary cases. Our in vitro experiments showed that SOX11 binds to MSI2 regulatory regions increasing its expression in MCL cells. MSI2 enhancers located in intronic regions were activated in SOX11+ but not in SOX11- MCL cell lines and primary cases. MSI2 upregulation significantly associated with poor overall survival, independently of common risk factors (SOX11 expression, high copy number alterations, 17p/TP53, 11q/ATM and 9p/CDKN2A alterations), in MCL patients. MSI2 knockdown MCL cell lines had a gene expression profile enriched in proapoptotic-related genes; whereas HSC and LSC gene signatures were downregulated compared to its control cells. MSI2 knockdown or MSI2 inhibition with Ro 08-2750 treatment decreased clonogenic growth, tumor survival and chemoresistance in MCL cells. SOX11+/MSI2 high MCL primary leukemic samples showed higher percentage of ALDH+ cells compared to SOX11-/MSI2 low samples, which decreased after Ro 08-2750 treatment. Besides, MSI2 silencing delayed tumor growth in vivo in MCL xenograft mice models. Overall, our results suggest that SOX11+ MCL cases acquire stemness features through the upregulation of MSI2 expression, which promotes chemoresistance, self-renewal and tumor survival. MSI2 represents a novel biomarker for MCL-CSC and therapeutic target for relapsed MCL. Citation Format: Marta Sureda-Gómez, Patricia Balsas, Marta Leonor Rodríguez, Ferran Nadeu, Anna De Bolòs, Álvaro Eguileor, Marta Kulis, Giancarlo Castellano, José Ignacio Martin-Subero, Santiago Demajo, Pedro Jares, Eva Giné, Elias Campo, Virginia Amador. Deciphering the role of MSI2 as a regulator of mantle cell lymphoma stem-like properties [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2022; 2022 Apr 8-13. Philadelphia (PA): AACR; Cancer Res 2022;82(12_Suppl):Abstract nr 777.

Roles of the MYST Family in the Pathogenesis of Alzheimer's Disease via Histone or Non-histone Acetylation.

Alzheimer's disease (AD) is one of the most common neurodegenerative diseases and a major cause of death among elderly individuals. The etiology of AD involves a combination of genetic, environmental, and lifestyle factors. A number of epigenetic alterations in AD have recently been reported; for example, studies have found an increase in histone acetylation in patients with AD and the protective function of histone deacetylase inhibitors. The histone acetylases in the MYST family are involved in a number of key nuclear processes, such as gene-specific transcriptional regulation, DNA replication, and DNA damage response. Therefore, it is not surprising that they contribute to epigenetic regulation as an intermediary between genetic and environmental factors. MYST proteins also exert acetylation activity on non-histone proteins that are closely associated with the pathogenesis of AD. In this review, we summarized the current understanding of the roles of MYST acetyltransferases in physiological functions and pathological processes related to AD. Additionally, using published RNA-seq, ChIP-seq, and ChIP-chip data, we identified enriched pathways to further evaluate the correlation between MYST and AD. The recent research described in this review supports the importance of epigenetic modifications and the MYST family in AD, providing a basis for future functional studies.

Genome Wide Approaches to Identify Protein-DNA Interactions.

Transcription factors are DNA-binding proteins that play key roles in many fundamental biological processes. Unraveling their interactions with DNA is essential to identify their target genes and understand the regulatory network. Genome-wide identification of their binding sites became feasible thanks to recent progress in experimental and computational approaches. ChIP-chip, ChIP-seq, and ChIP-exo are three widely used techniques to demarcate genome-wide transcription factor binding sites. This review aims to provide an overview of these three techniques including their experiment procedures, computational approaches, and popular analytic tools. ChIP-chip, ChIP-seq, and ChIP-exo have been the major techniques to study genome- wide in vivo protein-DNA interaction. Due to the rapid development of next-generation sequencing technology, array-based ChIP-chip is deprecated and ChIP-seq has become the most widely used technique to identify transcription factor binding sites in genome-wide. The newly developed ChIP-exo further improves the spatial resolution to single nucleotide. Numerous tools have been developed to analyze ChIP-chip, ChIP-seq and ChIP-exo data. However, different programs may employ different mechanisms or underlying algorithms thus each will inherently include its own set of statistical assumption and bias. So choosing the most appropriate analytic program for a given experiment needs careful considerations. Moreover, most programs only have command line interface so their installation and usage will require basic computation expertise in Unix/Linux.

RNA polymerase II (RNAP II)-associated factors are recruited to tRNA loci, revealing that RNAP II– and RNAP III–mediated transcriptions overlap in yeast

In yeast (Saccharomyces cerevisiae), the synthesis of tRNAs by RNA polymerase III (RNAP III) down-regulates the transcription of the nearby RNAP II-transcribed genes by a mechanism that is poorly understood. To clarify the basis of this tRNA gene-mediated (TGM) silencing, here, conducting a bioinformatics analysis of available ChIP-chip and ChIP-sequencing genomic data from yeast, we investigated whether the RNAP III transcriptional machinery can recruit protein factors required for RNAP II transcription. An analysis of 46 genome-wide protein-density profiles revealed that 12 factors normally implicated in RNAP II-mediated gene transcription are more enriched at tRNA than at mRNA loci. These 12 factors typically have RNA-binding properties, participate in the termination stage of the RNAP II transcription, and preferentially localize to the tRNA loci by a mechanism that apparently is based on the RNAP III transcription level. The factors included two kinases of RNAP II (Bur1 and Ctk1), a histone demethylase (Jhd2), and a mutated form of a nucleosome-remodeling factor (Spt6) that have never been reported to be recruited to tRNA loci. Moreover, we show that the expression levels of RNAP II-transcribed genes downstream of tRNA loci correlate with the distance from the tRNA gene by a mechanism that depends on their orientation. These results are consistent with the notion that pre-tRNAs recruit RNAP II-associated factors, thereby reducing the availability of these factors for RNAP II transcription and contributing, at least in part, to the TGM-silencing mechanism.

A Multi-Laplacian Prior and Augmented Lagrangian Approach to the Exploratory Analysis of Time-Varying Gene and Transcriptional Regulatory Networks for Gene Microarray Data.

This paper proposes a novel multi-Laplacian prior (MLP) and augmented Lagrangian method (ALM) approach for gene interactions and putative transcription factors (TFs) identification from time-course gene microarray data. It employs a non-linear time-varying auto-regressive (N-TVAR) model and the Maximum-A-Posteriori-Probability method for incorporating the multi-Laplacian prior and the continuity constraint. The MLP allows connections to/from a gene to be better preserved for putative TF identification in non-stationarity gene regulatory network as compared with conventional L1-based penalties. Moreover, the ALM allows the resultant non-smooth L1-based penalties to be decoupled from the remaining smooth terms, so that the former and latter can be efficiently solved using a low-complexity proximity operator and smooth optimization technique, respectively. Synthetic and real time-course gene microarray datasets are tested to evaluate the performance of the proposed method. Experimental results show that the proposed method gives better accuracy and higher computational speed than our previous work using smoothed approximation. Moreover, its performance, without the use of ChIP-chip data, is found to be highly comparable with other state-of-the-art methods integrating both ChIP-chip and gene microarray data. It suggests that the proposed method may serve as a useful exploratory tool for putative TF identification with reduced experimental cost.

Genome-Wide Analysis of the Distinct Types of Chromatin Interactions in Arabidopsis thaliana.

The three-dimensional shapes of chromosomes regulate gene expression and genome function. Our knowledge of the role of chromatin interaction is evolving rapidly. Here, we present a study of global chromatin interaction patterns in Arabidopsis thaliana. High-throughput experimental techniques have been developed to map long-range interactions within chromatin. We have integrated data from multiple experimental sources including Hi-C, BS-seq, ChIP-chip and ChIP-seq data for 17 epigenetic marks and 35 transcription factors. We identified seven groups of interacting loci, which can be distinguished by their epigenetic profiles. Furthermore, the seven groups of interacting loci can be divided into three types of chromatin linkages based on expression status. We observed that two interacting loci sometimes share common epigenetic and transcription factor-binding profiles. Different groups of loci display very different relationships between epigenetic marks and the binding of transcription factors. Distinctive types of chromatin linkages exhibit different gene expression profiles. Our study unveils an entirely unexplored regulatory interaction, linking epigenetic profiles, transcription factor binding and the three-dimensional spatial organization of the Arabidopsis nuclear genome.

Hepatitis B virus X protein influences enrichment profiles of H3K9me3 on promoter regions in human hepatoma cell lines.

We previously showed that hepatitis B virus (HBV) X protein (HBx) could promote the trimethylation of histone H3 lysine 9 (H3K9me3) to repress tumor suppressor genes in hepatocellular carcinoma (HCC). In this work, we analyze 23,148 human promoters using ChIP-chip to determine the effects of HBx on H3K9me3 enrichments in hepatoma cells with transfection of HBx-expressing plasmid. Immunohistochemistry for HBx and H3K9me3 was performed in 21 cases of HBV-associated HCC tissues. We identified that H3K9me3 immunoreactivity was significantly correlated with HBx staining in HCC tissues. ChIP-chip data indicated that HBx remarkably altered promoter enrichments of H3K9me3 in hepatoma cells. We identified 25 gene promoters, whose H3K9me3 enrichments are significantly altered in hepatoma cells transfected HBx-expressing plasmid, including 19 gaining H3K9m3, and six losing this mark. Most of these genes have not been previously reported in HCC, and BTBD17, MIR6089, ZNF205-AS1 and ZP1 have not previously been linked to cancer; only two genes (DAB2IP and ZNF185) have been reported in HCC. Genomic analyses suggested that genes with the differential H3K9me3 enrichments function in diverse cellular pathways and many are involved in cancer development and progression.

Activation of Ftz-F1-Responsive Genes through Ftz/Ftz-F1 Dependent Enhancers.

The orphan nuclear receptor Ftz-F1 is expressed in all somatic nuclei in Drosophila embryos, but mutations result in a pair-rule phenotype. This was explained by the interaction of Ftz-F1 with the homeodomain protein Ftz that is expressed in stripes in the primordia of segments missing in either ftz-f1 or ftz mutants. Ftz-F1 and Ftz were shown to physically interact and coordinately activate the expression of ftz itself and engrailed by synergistic binding to composite Ftz-F1/Ftz binding sites. However, attempts to identify additional target genes on the basis of Ftz-F1/ Ftz binding alone has met with only limited success. To discern rules for Ftz-F1 target site selection in vivo and to identify additional target genes, a microarray analysis was performed comparing wildtype and ftz-f1 mutant embryos. Ftz-F1-responsive genes most highly regulated included engrailed and nine additional genes expressed in patterns dependent on both ftz and ftz-f1. Candidate enhancers for these genes were identified by combining BDTNP Ftz ChIP-chip data with a computational search for Ftz-F1 binding sites. Of eight enhancer reporter genes tested in transgenic embryos, six generated expression patterns similar to the corresponding endogenous gene and expression was lost in ftz mutants. These studies identified a new set of Ftz-F1 targets, all of which are co-regulated by Ftz. Comparative analysis of enhancers containing Ftz/Ftz-F1 binding sites that were or were not bona fide targets in vivo suggested that GAF negatively regulates enhancers that contain Ftz/Ftz-F1 binding sites but are not actually utilized. These targets include other regulatory factors as well as genes involved directly in morphogenesis, providing insight into how pair-rule genes establish the body pattern.

The Nuclear Receptor, RORγ, Regulates Pathways Necessary for Breast Cancer Metastasis

We have previously reported that RORγ expression was decreased in ER−ve breast cancer, and increased expression improves clinical outcomes. However, the underlying RORγ dependent mechanisms that repress breast carcinogenesis have not been elucidated. Here we report that RORγ negatively regulates the oncogenic TGF-β/EMT and mammary stem cell (MaSC) pathways, whereas RORγ positively regulates DNA-repair. We demonstrate that RORγ expression is: (i) decreased in basal-like subtype cancers, and (ii) inversely correlated with histological grade and drivers of carcinogenesis in breast cancer cohorts. Furthermore, integration of RNA-seq and ChIP-chip data reveals that RORγ regulates the expression of many genes involved in TGF-β/EMT-signaling, DNA-repair and MaSC pathways (including the non-coding RNA, LINC00511). In accordance, pharmacological studies demonstrate that an RORγ agonist suppresses breast cancer cell viability, migration, the EMT transition (microsphere outgrowth) and mammosphere-growth. In contrast, RNA-seq demonstrates an RORγ inverse agonist induces TGF-β/EMT-signaling. These findings suggest pharmacological modulation of RORγ activity may have utility in breast cancer.

Abstract 1868A: Progesterone receptor promotes inflammatory gene programs in breast cancer

Abstract Mounting clinical data implicates progesterone and its receptor (progesterone receptor; PR) in breast cancer progression. Upon diagnosis, nearly 70% of breast cancers express PR and the estrogen receptor (ER). In contrast, only 7-10% of normal luminal epithelial cells express ER and PR. Mitogenic protein kinases (MAPK, ck2) known to modify and activate PR have been shown to be dysregulated in breast cancer. This creates a context wherein PR is inappropriately activated by aberrant kinases in the absence of hormones. We hypothesize that ck2-dependent phosphorylation of PR Serine 81 (Ser81) “redirects” the modified receptor to a distinct subset of target genes, including STAT5A, a gene required for mammary gland development and key to mediating inflammatory/interferon responses. Microarray data sets were analyzed using Gene Set Enrichment Analysis (GSEA) from progestin (R5020)-treated T47D breast cancer cells stably expressing either wt or Ser81 phospho-mutant PR (S79/81A PR). GSEA revealed that genes from multiple JAK/STAT, inflammatory and interferon-regulated gene sets were significantly enriched in our T47D-wt PR dataset. Enrichment of these gene sets was lost in cells lacking PR Ser81 phosphorylation (S79/81A PR), suggesting that Ser81 phosphorylation is necessary for their PR-dependent regulation. Leading Edge (LE) analysis revealed a significant enrichment in inflammatory/interferon-regulated pathways, again dependent upon PR Ser81 phosphorylation. Regulation of individual genes identified as being highly associated with the overlapping gene sets (LE gene list) was impaired in cells expressing S79/81A PR relative to cells expressing wt PR. Of note, many of these genes are transcriptional targets of interferon signaling pathways and key mediators of inflammation (i.e. IFITs, STATs, CXCL10, CXCL2, OAS1/2). Moreover, CEAS (cis-regulatory element annotation system) analysis of a publically available PR ChIP-chip data set (obtained from T47D breast cancer cells treated with progesterone) revealed significant enrichment of STAT5 binding motifs (GAS sites) within PR binding sites. Finally, using ChIP analysis we identified coordinate PR and STAT5 binding at key PR Ser81-regulated target genes; recruitment failed to occur in PR S79/81A-expressing cells. Cumulatively, these data support a model of coordinate regulation between PR and STAT5 on a select subset of Ser81 phosphorylation-dependent PR target genes. Many any of these genes are involved in cellular immune/inflammatory responses. We predict that direct interaction between STAT5 and PR at enhancer regions drives inflammatory gene programs in breast cancer. Understanding cross-talk between PR and inflammation is critical to developing better treatments, and could expose a novel mechanism through which PR drives early events in breast cancer development and progression that may be exploited using novel PR-based endocrine therapies. Citation Format: Christy R. Hagan, Todd P. Knutson, Carol A. Lange. Progesterone receptor promotes inflammatory gene programs in breast cancer. [abstract]. In: Proceedings of the 106th Annual Meeting of the American Association for Cancer Research; 2015 Apr 18-22; Philadelphia, PA. Philadelphia (PA): AACR; Cancer Res 2015;75(15 Suppl):Abstract nr 1868A. doi:10.1158/1538-7445.AM2015-1868A

ChIP-on-chip analysis methods for Affymetrix tiling arrays.

Although the ChIP-sequencing has gained significant attraction recently, ChIP analysis using microarrays is still an attractive option due to the low cost, ease of analysis, and access to legacy and public data sets. The analysis of ChIP-Chip data entails a multistep approach that requires several different applications to progress from the initial stages of raw data analysis to the identification and characterization of ChIP binding sites. There are multiple approaches to data analysis and there are several applications available for each stage of the analysis pipeline. Each application must be evaluated for its suitability for the particular experiment as well as the investigator's background with computational tools. This chapter is a review of the commonly available applications for Affymetrix ChIP-Chip data analysis, as well as the general workflow of a ChIP-Chip analysis approach. The purpose of the chapter is to allow the researcher to better select the appropriate applications and provide them with the direction necessary to proceed with a ChIP-Chip analysis.

Metaseq: a Python package for integrative genome-wide analysis reveals relationships between chromatin insulators and associated nuclear mRNA.

Here we introduce metaseq, a software library written in Python, which enables loading multiple genomic data formats into standard Python data structures and allows flexible, customized manipulation and visualization of data from high-throughput sequencing studies. We demonstrate its practical use by analyzing multiple datasets related to chromatin insulators, which are DNA-protein complexes proposed to organize the genome into distinct transcriptional domains. Recent studies in Drosophila and mammals have implicated RNA in the regulation of chromatin insulator activities. Moreover, the Drosophila RNA-binding protein Shep has been shown to antagonize gypsy insulator activity in a tissue-specific manner, but the precise role of RNA in this process remains unclear. Better understanding of chromatin insulator regulation requires integration of multiple datasets, including those from chromatin-binding, RNA-binding, and gene expression experiments. We use metaseq to integrate RIP- and ChIP-seq data for Shep and the core gypsy insulator protein Su(Hw) in two different cell types, along with publicly available ChIP-chip and RNA-seq data. Based on the metaseq-enabled analysis presented here, we propose a model where Shep associates with chromatin cotranscriptionally, then is recruited to insulator complexes in trans where it plays a negative role in insulator activity.

Direct and Indirect Targets of the E2A-PBX1 Leukemia-Specific Fusion Protein

E2A-PBX1 is expressed as a result of the t(1;19) chromosomal translocation in nearly 5% of cases of childhood acute lymphoblastic leukemia. The E2A-PBX1 chimeric transcription factor contains the N-terminal transactivation domain of E2A (TCF3) fused to the C-terminal DNA-binding homeodomain of PBX1. While there is no doubt of its oncogenic potential, the mechanisms of E2A-PBX1-mediated pre-B cell transformation and the nature of direct E2A-PBX1 target genes and pathways remain largely unknown. Herein we used chromatin immunoprecipitation assays (ChIP-chip) to identify direct targets of E2A-PBX1, and we used gene expression arrays of siRNA E2A-PBX1-silenced cells to evaluate changes in expression induced by the fusion protein. Combined ChIP-chip and expression data analysis gave rise to direct and functional targets of E2A-PBX1. Further we observe that the set of ChIP-chip identified E2A-PBX1 targets show a collective down-regulation trend in the E2A-PBX1 silenced samples compared to controls suggesting an activating role of this fusion transcription factor. Our data suggest that the expression of the E2A-PBX1 fusion gene interferes with key regulatory pathways and functions of hematopoietic biology. Among these are members of the WNT and apoptosis/cell cycle control pathways, and thus may comprise an essential driving force for the propagation and maintenance of the leukemic phenotype. These findings may also provide evidence of potentially attractive therapeutic targets.

Lysine Acetyltransferase GCN5b Interacts with AP2 Factors and Is Required for Toxoplasma gondii Proliferation

Histone acetylation has been linked to developmental changes in gene expression and is a validated drug target of apicomplexan parasites, but little is known about the roles of individual histone modifying enzymes and how they are recruited to target genes. The protozoan parasite Toxoplasma gondii (phylum Apicomplexa) is unusual among invertebrates in possessing two GCN5-family lysine acetyltransferases (KATs). While GCN5a is required for gene expression in response to alkaline stress, this KAT is dispensable for parasite proliferation in normal culture conditions. In contrast, GCN5b cannot be disrupted, suggesting it is essential for Toxoplasma viability. To further explore the function of GCN5b, we generated clonal parasites expressing an inducible HA-tagged dominant-negative form of GCN5b containing a point mutation that ablates enzymatic activity (E703G). Stabilization of this dominant-negative GCN5b was mediated through ligand-binding to a destabilization domain (dd) fused to the protein. Induced accumulation of the ddHAGCN5b(E703G) protein led to a rapid arrest in parasite replication. Growth arrest was accompanied by a decrease in histone H3 acetylation at specific lysine residues as well as reduced expression of GCN5b target genes in GCN5b(E703G) parasites, which were identified using chromatin immunoprecipitation coupled with microarray hybridization (ChIP-chip). Proteomics studies revealed that GCN5b interacts with AP2-domain proteins, apicomplexan plant-like transcription factors, as well as a “core complex” that includes the co-activator ADA2-A, TFIID subunits, LEO1 polymerase-associated factor (Paf1) subunit, and RRM proteins. The dominant-negative phenotype of ddHAGCN5b(E703G) parasites, considered with the proteomics and ChIP-chip data, indicate that GCN5b plays a central role in transcriptional and chromatin remodeling complexes. We conclude that GCN5b has a non-redundant and indispensable role in regulating gene expression required during the Toxoplasma lytic cycle.

A review of three different studies on hidden markov models for epigenetic problems: a computational perspective.

Recent technical advances, such as chromatin immunoprecipitation combined with DNA microarrays (ChIp-chip) and chromatin immunoprecipitation-sequencing (ChIP-seq), have generated large quantities of high-throughput data. Considering that epigenomic datasets are arranged over chromosomes, their analysis must account for spatial or temporal characteristics. In that sense, simple clustering or classification methodologies are inadequate for the analysis of multi-track ChIP-chip or ChIP-seq data. Approaches that are based on hidden Markov models (HMMs) can integrate dependencies between directly adjacent measurements in the genome. Here, we review three HMM-based studies that have contributed to epigenetic research, from a computational perspective. We also give a brief tutorial on HMM modelling-targeted at bioinformaticians who are new to the field.

Systematic discovery and characterization of regulatory motifs in ENCODE TF binding experiments

Recent advances in technology have led to a dramatic increase in the number of available transcription factor ChIP-seq and ChIP-chip data sets. Understanding the motif content of these data sets is an important step in understanding the underlying mechanisms of regulation. Here we provide a systematic motif analysis for 427 human ChIP-seq data sets using motifs curated from the literature and also discovered de novo using five established motif discovery tools. We use a systematic pipeline for calculating motif enrichment in each data set, providing a principled way for choosing between motif variants found in the literature and for flagging potentially problematic data sets. Our analysis confirms the known specificity of 41 of the 56 analyzed factor groups and reveals motifs of potential cofactors. We also use cell type-specific binding to find factors active in specific conditions. The resource we provide is accessible both for browsing a small number of factors and for performing large-scale systematic analyses. We provide motif matrices, instances and enrichments in each of the ENCODE data sets. The motifs discovered here have been used in parallel studies to validate the specificity of antibodies, understand cooperativity between data sets and measure the variation of motif binding across individuals and species.

SOX11 Regulates Angiogenesis In Mantle Cell Lymphoma

Mantle cell lymphoma (MCL) is an aggressive tumor but recent studies have identified a subtype of this lymphoma with an indolent clinical behavior and long survival of the patients even without chemotherapy that may correspond to a different clinical and biological subtype of the disease. These indolent MCLs differ from conventional MCLs in a particular gene signature that lacks the expression of some transcription factors of the High Mobility Group (HMG). SOX11 is one of the best discriminatory genes between these conventional and indolent MCL tumors. We have recently demonstrated the oncogenic implication of SOX11 expression in the aggressive behavior of MCL as SOX11-knockdown derived tumors display a significant reduction on tumor growth compared to SOX11-control tumors in subcutaneous MCL xenograft experiments (Vegliante MC, Palomero J et al. Blood. 2013; 121:2175-85). To uncover SOX11 regulated genes and transcriptional programs in MCL, we performed an integrative analysis coupling data from human genome-wide promoter analysis by SOX11 ChIP-chip experiments and gene expression profiling (GEP) upon SOX11 silencing in MCL cell lines. In these studies we initially identified the role of SOX11 in blocking the terminal B-cell differentiation by the direct positive regulation of PAX5. However, how SOX11 may enhance tumor growth is still not well known.To further characterize the potential oncogenic mechanisms regulated by SOX11, we have integrated our ChIP-chip data with further analyses including GEP derived from the xenografted SOX11-positive and silenced tumors. The gene ontology term analysis of the genome-wide promoter study revealed “blood vessel development” as one of the most significant biological processes overrepresented among the SOX11-bound genes. Concordantly, a gene set enrichment analysis (GSEA) of the differentially expressed genes in the xenografted tumors showed a significant higher enrichment of gene signatures related to tumor angiogenesis in the SOX11-positive tumors. To validate these results in vivo, we first investigated protein extracts of these tumors using an angiogenesis proteome profiler antibody array and discovered a significant higher expression of 20 pro-angiogenic factors in the SOX11-positive tumors. We have then investigated the angiogenic development in tissue sections from four Z-138 SOX11-positive and nine SOX11-silenced xenografted tumors. The microvascular density area recognized by CD31 staining was significantly larger in SOX11-positive than knockdown tumors (90%±µ2m vs 15%±µ2m, respectively, p<1x10-4). SOX11-silenced tumors had larger necrotic areas than positive tumors (20%±µ2m vs 1%±µ2m, respectively, p<1x10-4), suggesting that vascular paucity could contribute to the significant smaller volume of these tumors.To determine whether angiogenic development was also differentially represented in primary human MCL we have studied the GEP of 16 SOX11-positive and 22 SOX11-negative MCL. A GSEA revealed that SOX11 expressing tumors were enriched in signatures related to tumor vasculature, angiogenesis and vasculature development. Concordantly, the CD34 staining of an independent series of primary SOX11-positive (n=8) and negative (n=9) human MCL showed higher microvascular density in SOX11 expressing tumors (30%±µ2m vs 5%±µ2m, respectively, p<1x10-4).In conclusion, these findings indicate that SOX11 promotes angiogenesis in experimental and primary human MCL and this may be an oncogenic mechanism contributing to the aggressiveness of these tumors. Disclosures:No relevant conflicts of interest to declare.

E2F2 and CREB cooperatively regulate transcriptional activity of cell cycle genes

E2F2 is essential for the maintenance of T lymphocyte quiescence. To identify the full set of E2F2 target genes, and to gain further understanding of the role of E2F2 in transcriptional regulation, we have performed ChIP-chip analyses across the genome of lymph node–derived T lymphocytes. Here we show that during quiescence, E2F2 binds the promoters of a large number of genes involved in DNA metabolism and cell cycle regulation, concomitant with their transcriptional silencing. A comparison of ChIP-chip data with expression profiling data on resting E2f2−/− T lymphocytes identified a subset of 51 E2F2-specific target genes, most of which are upregulated on E2F2 loss. Luciferase reporter assays showed a retinoblastoma-independent role for E2F2 in the negative regulation of these target genes. Importantly, we show that the DNA binding activity of the transcription factor CREB contributes to E2F2-mediated repression of Mcm5 and Chk1 promoters. siRNA-mediated CREB knockdown, expression of a dominant negative KCREB mutant or disruption of CREB binding by mutating a CRE motif on Mcm5 promoter, relieved E2F2-mediated transcriptional repression. Taken together, our data uncover a new regulatory mechanism for E2F-mediated transcriptional control, whereby E2F2 and CREB cooperate in the transcriptional repression of a subset of E2F2 target genes.