HDAC inhibitors cause site-specific chromatin remodeling at PU.1-bound enhancers in K562 cells

BACKGROUND: Small molecule inhibitors of histone deacetylases (HDACi) hold promise as anticancer agents for particular malignancies. However, clinical use is often confounded by toxicity, perhaps due to indiscriminate hyperacetylation of cellular proteins. Therefore, elucidating the mechanisms by which HDACi trigger differentiation, cell cycle arrest, or apoptosis of cancer cells could inform development of more targeted therapies. We used the myelogenous leukemia line K562 as a model of HDACi-induced differentiation to investigate chromatin accessibility (DNase-seq) and expression (RNA-seq) changes associated with this process. RESULTS: We identified several thousand specific regulatory elements [~10 % of total DNase I-hypersensitive (DHS) sites] that become significantly more or less accessible with sodium butyrate or suberanilohydroxamic acid treatment. Most of the differential DHS sites display hallmarks of enhancers, including being enriched for non-promoter regions, associating with nearby gene expression changes, and increasing luciferase reporter expression in K562 cells. Differential DHS sites were enriched for key hematopoietic lineage transcription factor motifs, including SPI1 (PU.1), a known pioneer factor. We found PU.1 increases binding at opened DHS sites with HDACi treatment by ChIP-seq, but PU.1 knockdown by shRNA fails to block the chromatin accessibility and expression changes. A machine-learning approach indicates H3K27me3 initially marks PU.1-bound sites that open with HDACi treatment, suggesting these sites are epigenetically poised. CONCLUSIONS: We find HDACi treatment of K562 cells results in site-specific chromatin remodeling at epigenetically poised regulatory elements. PU.1 shows evidence of a pioneer role in this process by marking poised enhancers but is not required for transcriptional activation.

Genome-wide analysis of p63 binding sites identifies AP-2 factors as co-regulators of epidermal differentiation

The p63 transcription factor (TP63) is critical in development, growth and differentiation of stratifying epithelia. This is highlighted by the severity of congenital abnormalities caused by TP63 mutations in humans, the dramatic phenotypes in knockout mice and de-regulation of TP63 expression in neoplasia altering the tumour suppressive roles of the TP53 family. In order to define the normal role played by TP63 and provide the basis for better understanding how this network is perturbed in disease, we used chromatin immunoprecipitation combined with massively parallel sequencing (ChIP-seq) to identify >7500 high-confidence TP63-binding regions across the entire genome, in primary human neonatal foreskin keratinocytes (HFKs). Using integrative strategies, we demonstrate that only a subset of these sites are bound by TP53 in response to DNA damage. We identify a role for TP63 in transcriptional regulation of multiple genes genetically linked to cleft palate and identify AP-2alpha (TFAP2A) as a co-regulator of a subset of these genes. We further demonstrate that AP-2gamma (TFAP2C) can bind a subset of these regions and that acute depletion of either TFAP2A or TFAP2C alone is sufficient to reduce terminal differentiation of organotypic epidermal skin equivalents, indicating overlapping physiological functions with TP63.

Global identification of androgen response elements

Chromatin immunoprecipitation (ChIP) is an invaluable tool in the study of transcriptional regulation. ChIP methods require both a priori knowledge of the transcriptional regulators which are important for a given biological system and high-quality specific antibodies for these targets. The androgen receptor (AR) is known to play essential roles in male sexual development, in prostate cancer and in the function of many other AR-expressing cell types (e.g. neurons and myocytes). As a ligand-activated transcription factor the AR also represents an endogenous, inducible system to study transcriptional biology. Therefore, ChIP studies of the AR can make use of treatment contrast experiments to define its transcriptional targets. To date several studies have mapped AR binding sites using ChIP in combination with genome tiling microarrays (ChIP-chip) or direct sequencing (ChIP-seq), mainly in prostate cancer cell lines and with varying degrees of genomic coverage. These studies have provided new insights into the DNA sequences to which the AR can bind, identified AR cooperating transcription factors, mapped thousands of potential AR regulated genes and provided insights into the biological processes regulated by the AR. However, further ChIP studies will be required to fully characterise the dynamics of the AR-regulated transcriptional programme, to map the occupancy of different AR transcriptional complexes which result in different transcriptional output and to delineate the transcriptional networks downstream of the AR.

Novel role of androgens in mitochondrial fission and apoptosis

Androgen and androgen receptors (AR) play critical roles in the proliferation of prostate cancer through transcriptional regulation of target genes. Here, we found that androgens upregulated the expression of dynamin-related protein 1 (Drp1), which is involved in the induction of mitochondrial fission, a common event in mitosis and apoptosis. Clinical tissue samples and various prostate cancer cell lines revealed a positive correlation between Drp1 and AR levels. Treatment of androgen-sensitive cells with an AR agonist, R1881, and antagonist, bicalutamide, showed that Drp1 is transcriptionally regulated by androgens, as confirmed by an AR ChIP-seq assay. Live imaging experiments using pAcGFP1-Mito stably transfected LNCaP (mito-green) cells revealed that androgen did not induce significant mitochondrial fission by itself, although Drp1 was upregulated. However, when treated with CGP37157 (CGP), an inhibitor of mitochondrial Ca²⁺ efflux, these cells exhibited mitochondrial fission, which was further enhanced by pretreatment with R1881, suggesting that androgen-induced Drp1 expression facilitated CGP-induced mitochondrial fission. This enhanced mitochondrial fission was correlated with increased apoptosis. Transfection with dominant-negative (DN-Drp1, K38A) rescued cells from increased apoptosis, confirming the role of androgen-induced Drp1 in the observed apoptosis with combination treatment. Furthermore, we found that CGP reduced the expression of Mfn1, a protein that promotes mitochondrial fusion, a process which opposes fission. We suggest that androgen-increased Drp1 enhanced mitochondrial fission leading to apoptosis. The present study shows a novel role for androgens in the regulation of mitochondrial morphology that could potentially be utilized in prostate cancer therapy.

CLOCK-controlled polyphonic regulation of circadian rhythms through canonical and noncanonical E-boxes.

In mammalian circadian clockwork, the CLOCK-BMAL1 complex binds to DNA enhancers of target genes and drives circadian oscillation of transcription. Here we identified 7,978 CLOCK-binding sites in mouse liver by chromatin immunoprecipitation-sequencing (ChIP-Seq), and a newly developed bioinformatics method, motif centrality analysis of ChIP-Seq (MOCCS), revealed a genome-wide distribution of previously unappreciated noncanonical E-boxes targeted by CLOCK. In vitro promoter assays showed that CACGNG, CACGTT, and CATG(T/C)G are functional CLOCK-binding motifs. Furthermore, we extensively revealed rhythmically expressed genes by poly(A)-tailed RNA-Seq and identified 1,629 CLOCK target genes within 11,926 genes expressed in the liver. Our analysis also revealed rhythmically expressed genes that have no apparent CLOCK-binding site, indicating the importance of indirect transcriptional and posttranscriptional regulations. Indirect transcriptional regulation is represented by rhythmic expression of CLOCK-regulated transcription factors, such as Krüppel-like factors (KLFs). Indirect posttranscriptional regulation involves rhythmic microRNAs that were identified by small-RNA-Seq. Collectively, CLOCK-dependent direct transactivation through multiple E-boxes and indirect regulations polyphonically orchestrate dynamic circadian outputs.

Génomique fonctionnelle des cellules corticotropes hypophysaires : contrôle génétique de la gestion systémique des stress

L'axe hypothalamo-hypophyso-surrénalien (HPA) permet de maintenir l'homéostasie de l'organisme face à divers stress. Qu'ils soient de nature psychologique, physique ou inflammatoire/infectieux, les stress provoquent la synthèse et la libération de CRH par l'hypothalamus. Les cellules corticotropes hypophysaires perçoivent ce signal et en réaction, produisent et sécrètent l'ACTH. Ceci induit la synthèse des glucocorticoïdes (Gc) par le cortex surrénalien; ces stéroïdes mettent le système métabolique en état d’alerte pour la réponse au stress et à l’agression. Les Gc ont le rôle essentiel de contrôler les défenses de l'organisme, en plus d'exercer une rétro-inhibition sur l'axe HPA. L'ACTH est une petite hormone peptidique produite par le clivage d'un précurseur: la pro-opiomélanocortine (POMC). À cause de sa position critique dans la normalisation de l'homéostasie, le contrôle transcriptionnel du gène Pomc a fait l'objet d'études approfondies au cours des dernières décennies. Nous savons maintenant que la région promotrice du gène Pomc permet une expression ciblée dans les cellules POMC hypophysaires. L'étude du locus Pomc par des technologies génomiques m'a permis de découvrir un nouvel élément de régulation qui est conservé à travers l'évolution des mammifères. La caractérisation de cet enhancer a démontré qu'il dirige une expression restreinte à l'hypophyse, et plus particulièrement dans les cellules corticotropes. De façon intéressante, l'activité de cet élément dépend d'un nouveau site de liaison recrutant un homodimère du facteur de transcription Tpit, dont l'expression est également limitée aux cellules POMC de l'hypophyse. La découverte de cet enhancer ajoute une toute nouvelle dimension à la régulation de l'expression de POMC. Les cytokines pro-inflammatoires IL6/LIF et les Gc sont connus pour leur antagonisme sur la réaction inflammatoire et sur le promoteur Pomc via l'action des facteurs de transcription Stat3 et GR respectivement. L'analyse génomique des sites liés ii par ces deux facteurs nous a révélé une interrelation complexe et a permis de définir un code transcriptionnel entre ces voies de signalisation. En plus de leur action par interaction directe avec l’ADN au niveau des séquences régulatrices, ces facteurs interagissent directement entre eux avec des résultats transcriptionnels différents. Ainsi, le recrutement de GR par contact protéine:protéine (tethering) sur Stat3 étant lié à l'ADN provoque un antagonisme transcriptionnel. Inversement, le tethering de Stat3 sur GR supporte une action synergique, tout comme leur co-recrutement à l'ADN sur des sites contigus ou composites. Lors d'une activation soutenue, ce synergisme entre les voies IL6/LIF et Gc induit une réponse innée de défense cellulaire. Ainsi lors d'un stress majeur, ce mécanisme de défense est mis en branle dans toutes les cellules et tissus. En somme, les travaux présentés dans cette thèse définissent les mécanismes transcriptionnels engagés dans le combat de l'organisme contre les stress. Plus particulièrement, ces mécanismes ont été décrits au niveau de la réponse globale des corticotropes et du gène Pomc. Il est essentiel pour l'organisme d'induire adéquatement ces mécanismes afin de faire face aux stress et d'éviter des dérèglements comme les maladies inflammatoires et métaboliques.

Understanding the transcriptional control of EIF4E and its dysregulation in acute myeloid leukemia: role of NF-κB

EIF4E, le facteur d’initiation de la traduction chez les eucaryotes est un oncogène puissant et qui se trouve induit dans plusieurs types de cancers, parmi lesquels les sous-types M4 et M5 de la leucémie aiguë myéloblastique (LAM). EIF4E est régulé à plusieurs niveaux cependant, la régulation transcriptionnelle de ce gène est peu connue. Mes résultats montrent que EIF4E est une cible transcriptionnelle directe du facteur nucléaire « kappa-light- chain- enhancer of activated B cells » (NF-κB).Dans les cellules hématopoïétiques primaires et les lignées cellulaires, les niveaux de EIF4E sont induits par des inducteurs de NF-κB. En effet, l’inactivation pharmaceutique ou génétique de NF-κB réprime l’activation de EIF4E. En effet, suite à l’activation de NF-κB chez l’humain, le promoteur endogène de EIF4E recrute p65 (RelA) et c-Rel aux sites évolutionnaires conservés κB in vitro et in vivo en même temps que p300 ainsi que la forme phosphorylée de Pol II. De plus, p65 est sélectivement associé au promoteur de EIF4E dans les sous-types LAM M4/M5 mais non pas dans les autres sous-types LAM ou dans les cellules hématopoïétiques primaires normales. Ceci indique que ce processus représente un facteur essentiel qui détermine l’expression différentielle de EIF4E dans la LAM. Les analyses de données d’expressions par séquençage de l’ARN provenant du « Cancer Genome Atlas » (TCGA) suggèrent que les niveaux d’ARNm de EIF4E et RELA se trouvent augmentés dans les cas LAM à pronostic intermédiaire ou faible mais non pas dans les groupes cytogénétiquement favorables. De plus, des niveaux élevés d’ARNm de EIF4E et RELA sont significativement associés avec un taux de survie relativement bas chez les patients. En effet, les sites uniques κB se trouvant dans le promoteur de EIF4E recrutent le régulateur de transcription NF-κB p65 dans 47 nouvelles cibles prévues. Finalement, 6 nouveaux facteurs de transcription potentiellement impliqués dans la régulation du gène EIF4E ont été prédits par des analyses de données ChIP-Seq provenant de l’encyclopédie des éléments d’ADN (ENCODE). Collectivement, ces résultats fournissent de nouveaux aperçus sur le control transcriptionnel de EIF4E et offrent une nouvelle base moléculaire pour sa dérégulation dans au moins un sous-groupe de spécimens de LAM. L’étude et la compréhension de ce niveau de régulation dans le contexte de spécimens de patients s’avère important pour le développement de nouvelles stratégies thérapeutiques ciblant l’expression du gène EIF4E moyennant des inhibiteurs de NF-κB en combinaison avec la ribavirine.

Identification de nouvelles cibles transcriptionnelles d’ETV6, un facteur de transcription fréquemment altéré dans la leucémie aiguë lymphoblastique

La leucémie aiguë lymphoblastique (LAL) est le cancer pédiatrique le plus fréquent. Plusieurs réarrangements chromosomiques ont été associés à cette maladie, dont la translocation t(12;21), qui est observée dans 25% des cas de LAL de type pré-B. Cette translocation engendre l’expression de la protéine de fusion ETV6-AML1. Toutefois, celle-ci n’est pas suffisante pour initier seule une leucémie, ce qui suggère que des mutations additionnelles sont nécessaires à la transformation oncogénique. Or, on observe que l’allèle non-réarrangé d’ETV6 est perdu dans 75% des cas de t(12;21). Cette délétion entraîne l’inactivation complète du facteur de transcription ETV6 et l’abolition de sa fonction biologique. Puisqu’ETV6 semble jouer un rôle de suppresseur de tumeurs, nous croyons que son inactivation favoriserait le développement de la leucémie via la dérégulation de ses gènes cibles. Ce projet visait donc à identifier de nouvelles cibles transcriptionnelles d’ETV6, afin d’élucider son implication dans la leucémie. Une expérience de RNA-Seq a permis d’identifier plus de 200 gènes dont l’expression est corrélée avec celle d’ETV6 dans des cellules souches hématopoïétiques CD34+. Parmi ceux-ci, plusieurs gènes sont impliqués dans la réponse immunitaire et inflammatoire, la migration cellulaire, l’homéostasie ionique et la signalisation intracellulaire. Nous avons également mis en place une approche d’immunoprécipitation de la chromatine afin d’identifier les régions auxquelles le facteur de transcription ETV6 peut se lier. À l’aide de cette méthode, nous avons démontré une interaction entre ETV6 et SLCO2B1, un gène dont l’expression est également co-régulée avec ETV6. Finalement, notre étude suggère qu’ETV6 contribuerait à la leucémogenèse en dérégulant l’expression de certains gènes ayant des propriétés oncogéniques.

Trans-Ancestral Studies Fine Map the SLE-Susceptibility Locus TNFSF4

We previously established an 80 kb haplotype upstream of TNFSF4 as a susceptibility locus in the autoimmune disease SLE. SLE-associated alleles at this locus are associated with inflammatory disorders, including atherosclerosis and ischaemic stroke. In Europeans, the TNFSF4 causal variants have remained elusive due to strong linkage disequilibrium exhibited by alleles spanning the region. Using a trans-ancestral approach to fine-map the locus, utilising 17,900 SLE and control subjects including Amerindian/Hispanics (1348 cases, 717 controls), African-Americans (AA) (1529, 2048) and better powered cohorts of Europeans and East Asians, we find strong association of risk alleles in all ethnicities; the AA association replicates in African-American Gullah (152,122). The best evidence of association comes from two adjacent markers: rs2205960-T (P = 1.71×10-34, OR = 1.43[1.26-1.60]) and rs1234317-T (P = 1.16×10-28, OR = 1.38[1.24-1.54]). Inference of fine-scale recombination rates for all populations tested finds the 80 kb risk and non-risk haplotypes in all except African-Americans. In this population the decay of recombination equates to an 11 kb risk haplotype, anchored in the 5′ region proximal to TNFSF4 and tagged by rs2205960-T after 1000 Genomes phase 1 (v3) imputation. Conditional regression analyses delineate the 5′ risk signal to rs2205960-T and the independent non-risk signal to rs1234314-C. Our case-only and SLE-control cohorts demonstrate robust association of rs2205960-T with autoantibody production. The rs2205960-T is predicted to form part of a decameric motif which binds NF-κBp65 with increased affinity compared to rs2205960-G. ChIP-seq data also indicate NF-κB interaction with the DNA sequence at this position in LCL cells. Our research suggests association of rs2205960-T with SLE across multiple groups and an independent non-risk signal at rs1234314-C. rs2205960-T is associated with autoantibody production and lymphopenia. Our data confirm a global signal at TNFSF4 and a role for the expressed product at multiple stages of lymphocyte dysregulation during SLE pathogenesis. We confirm the validity of trans-ancestral mapping in a complex trait. © 2013 Manku et al.

Non-referenced genome assembly from epigenomic short-read data

Current computational methods used to analyze changes in DNA methylation and chromatin modification rely on sequenced genomes. Here we describe a pipeline for the detection of these changes from short-read sequence data that does not require a reference genome. Open source software packages were used for sequence assembly, alignment, and measurement of differential enrichment. The method was evaluated by comparing results with reference-based results showing a strong correlation between chromatin modification and gene expression. We then used our de novo sequence assembly to build the DNA methylation profile for the non-referenced Psammomys obesus genome. The pipeline described uses open source software for fast annotation and visualization of unreferenced genomic regions from short-read data.

Paediatric leukaemia DNA methylation profiling using MBD enrichment and SOLiD sequencing on archival bone marrow smears

BACKGROUND: Acute Lymphoblastic Leukaemia (ALL) is the most common cancer in children. Over the past four decades, research has advanced the treatment of this cancer from a less than 60% chance of survival to over 85% today. The causal molecular mechanisms remain unclear. Here, we performed sequencing-based genomic DNA methylation profiling of eight paediatric ALL patients using archived bone marrow smear microscope slides. FINDINGS: SOLiD™ sequencing data was collected from Methyl-Binding Domain (MBD) enriched fractions of genomic DNA. The primary tumour and remission bone marrow sample was analysed from eight patients. Four patients relapsed and the relapsed tumour was analysed. Input and MBD-enriched DNA from each sample was sequenced, aligned to the hg19 reference genome and analysed for enrichment peaks using MACS (Model-based Analysis for ChIP-Seq) and HOMER (Hypergeometric Optimization of Motif EnRichment). In total, 3.67 gigabases (Gb) were sequenced, 2.74 Gb were aligned to the reference genome (average 74.66% alignment efficiency). This dataset enables the interrogation of differential DNA methylation associated with paediatric ALL. Preliminary results reveal concordant regions of enrichment indicative of a DNA methylation signature. CONCLUSION: Our dataset represents one of the first SOLiD™MBD-Seq studies performed on paediatric ALL and is the first to utilise archival bone marrow smears. Differential DNA methylation between cancer and equivalent disease-free tissue can be identified and correlated with existing and published genomic studies. Given the rarity of paediatric haematopoietic malignancies, relative to adult counterparts, our demonstration of the utility of archived bone marrow smear samples to high-throughput methylation sequencing approaches offers tremendous potential to explore the role of DNA methylation in the aetiology of cancer.

The SUCEST-FUN regulatory network database: Designing an energy grass

Modern sugarcane cultivars are complex hybrids resulting from crosses among several Saccharum species. Traditional breeding methods have been employed extensively in different countries over the past decades to develop varieties with increased sucrose yield and resistance to pests and diseases. Conventional variety improvement, however, may be limited by the narrow pool of suitable genes. Thus, molecular genetics is seen as a promising tool to assist in the process of developing improved varieties. The SUCEST-FUN Project (http://sucest-fun.org) aims to associate function with sugarcane genes using a variety of tools, in particular those that enable the study of the sugarcane transcriptome. An extensive analysis has been conducted to characterise, phenotypically, sugarcane genotypes with regard to their sucrose content, biomass and drought responses. Through the analysis of different cultivars, genes associated with sucrose content, yield, lignin and drought have been identified. Currently, tools are being developed to determine signalling and regulatory networks in grasses, and to sequence the sugarcane genome, as well as to identify sugarcane promoters. This is being implemented through the SUCEST-FUN (http://sucest-fun.org) and GRASSIUS databases (http://grassius.org), the cloning of sugarcane promoters, the identification of cis-regulatory elements (CRE) using Chromatin Immunoprecipitation-sequencing (ChIP-Seq) and the generation of a comprehensive Signal Transduction and Transcription gene catalogue (SUCAST Catalogue).

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.

Hsf1 and Hsp90 orchestrate temperature-dependent global transcriptional remodelling and chromatin architecture in Candida albicans

We thank Karim Gharbi and Urmi Trivedi for their assistance with RNA sequencing, carried out in the GenePool genomics facility (University of Edinburgh). We also thank Susan Fairley and Eduardo De Paiva Alves (Centre for Genome Enabled Biology and Medicine, University of Aberdeen) for help with the initial bioinformatics analysis. We thank Aaron Mitchell for kindly providing the ALS3 mutant, Julian Naglik for the gift of TR146 cells, and Jon Richardson for technical assistance. We thank the Genomics and Bioinformatics core of the Faculty of Health Sciences for Next Generation Sequencing and Bioinformatics support, the Information and Communication Technology Office at the University of Macau for providing access to a High Performance Computer and Jacky Chan and William Pang for their expert support on the High Performance Computer. Finally, we thank Amanda Veri for generating CaLC2928. M.D.L. is supported by a Sir Henry Wellcome Postdoctoral Fellowship (Wellcome Trust 096072), R.A.F. by a Wellcome Trust-Massachusetts Institute of Technology (MIT) Postdoctoral Fellowship, L.E.C. by a Canada Research Chair in Microbial Genomics and Infectious Disease and by Canadian Institutes of Health Research Grants MOP-119520 and MOP-86452, A.J. P.B. was supported by the UK Biotechnology and Biological Sciences Research Council (BB/F00513X/1) and by the European Research Council (ERC-2009-AdG-249793-STRIFE), KHW is supported by the Science and Technology Development Fund of Macau S.A.R (FDCT) (085/2014/A2) and the Research and Development Administrative Office of the University of Macau (SRG2014-00003-FHS) and R.T.W. by the Burroughs Wellcome fund and NIH R15AO094406. Data availability RNA-sequencing data sets are available at ArrayExpress (www.ebi.ac.uk) under accession code E-MTAB-4075. ChIP-seq data sets are available at the NCBI SRA database (http://www.ncbi.nlm.nih.gov) under accession code SRP071687. The authors declare that all other data supporting the findings of this study are available within the article and its supplementary information files, or from the corresponding author upon request.

Understanding human genome regulation through entropic graph-based regularization and submodular optimization

Thesis (Ph.D.)--University of Washington, 2016-06