In recent years, the development of new genomic technologies has allowed for the investigation of many regulatory epigenetic marks besides expression levels, on a genome-wide scale. As the price for these technologies continues to decrease, study sizes will not only increase, but several different assays are beginning to be used for the same samples. It is therefore desirable to develop statistical methods to integrate multiple data types that can handle the increased computational burden of incorporating large data sets. Furthermore, it is important to develop sound quality control and normalization methods as technical errors can compound when integrating multiple genomic assays. DNA methylation is a commonly studied epigenetic mark, and the Infinium HumanMethylation450 BeadChip has become a popular microarray that provides genome-wide coverage and is affordable enough to scale to larger study sizes. It employs a complex array design that has complicated efforts to develop normalization methods. We propose a novel normalization method that uses a set of stable methylation sites from housekeeping genes as empirical controls to fit a local regression hypersurface to signal intensities. We demonstrate that our method performs favorably compared to other popular methods for the array. We also discuss an approach to estimating cell-type admixtures, which is a frequent biological confound in these studies. For data integration we propose a gene-centric procedure that uses canonical correlation and subsequent permutation testing to examine correlation or other measures of association and co-localization of epigenetic marks on the genome. Specifically, a likelihood ratio test for general association between data modalities is performed after an initial dimension reduction step. Canonical scores are then regressed against covariates of interest using linear mixed effects models. Lastly, permutation testing is performed on weighted correlation matrices to test for co-localization of relationships to physical locations in the genome. We demonstrate these methods on a set of developmental brain samples from the BrainSpan consortium and find substantial relationships between DNA methylation, gene expression, and alternative promoter usage primarily in genes related to axon guidance. We perform a second integrative analysis on another set of brain samples from the Stanley Medical Research Institute.

The genes of the human β-globin locus are located on chromosome 11 in the order of their expression during development: 5' ε, γ, β 3'. During development, silencing of the 5' gene occurs with activation of the immediate 3' gene. This process occurs twice and is termed hemoglobin switching. The exact mechanism(s) of this process have not been fully described. Herein, we describe a role for DNA methylation and methyl binding domain protein 2 in the transcriptional regulation of the human embryonic and fetal beta type globin genes. Adult mice containing the entire human β-globin locus as a yeast artificial chromosome (βYAC) express very low levels of the fetal γ-globin gene. However, treatment of adult βYAC transgenic mice with the DNA methyltransferase inhibitor, 5-azacytidine, induces a >10-fold increase γ-globin mRNA levels. In addition, βYAC transgenic mice null for methyl binding domain protein 2 (MBD2) express a similar level of γ-globin mRNA. DNA methylation and MBD2 appear to induce γ-globin expression via the same pathway(s), as treatment of MBD2 null βYAC transgenic mice do not show an additive boost in γ-globin expression. MBD2 does not bind to the γ-globin promoter region in vivo indicating MBD2 mediated transcriptional silencing does not occur by recruitment of transcriptional repression complexes to the γ-globin gene promoter. Additionally, these transgenic mice contain only the 5' portion of the β-globin locus through the ε-globin, and do not express the ε-globin genes as adults. However, treatment with 5-azacytidine or loss of MBD2 induces expression of the ε-globin gene in adult transgenic mice. A similar induction of ε-globin is seen in βYAC transgenic mice under the same conditions. The level of expression of the ε-globin gene is much lower than the γ-globin gene, indicating the powerful effect of the cis elements mediating transcriptional repression of the ε-globin gene. These studies indicate DNA methylation and MBD2 contribute to the transcriptional repression of the human embryonic and fetal β-type globin genes. Additionally, MBD2 has been identified as a potential target for the therapeutic induction of fetal hemoglobin for the treatment of hemoglobinopathies.

DNA methylation plays a pivotal role in governing cellular processes including genomic imprinting, gene expression, and development. Recently, the Tet family of methylcytosine dioxygenases(Tet1, Tet2 and Tet3) was found to catalyze the oxidation of 5-methylcytosine (5mC) to 5-hydroxymethylcytosine (5hmC), an intermediate in the pathway of DNA demethylation. Tet enzymes catalyze this hydroxylation in a 2-oxoglutarate and Fe2+ dependent manner. We have recently reported significant levels of 5mC and 5hmC modification in immunoprecipitates of mammalian mitochondrial DNA(mtDNA). We provide the first evidence that a DNA Methyltransferase-1 isoform (mtDNMT1) translocates to the mitochondria using an N-terminal mitochondrial targeting sequence. mtDNMT1 expression is upregulated by NRF1 and PGC1α, master regulators of mitochondrial biogenesis and function, as well as by loss of p53. Altered mtDNMT1 expression asymmetrically affects mtDNA transcription. We are now pursuing the role of Tet proteins in generating 5hmC in mtDNA. Using an in vitro enzyme assay, we have successfully detected Tet activity in crude and percoll purified mitochondrial fractions of HCT116 cells. Mitoprot analysis on Tet family predicts that Tet1 may be translocated to the mitochondria. Immunoblot analysis indicates that a band of expected size(235kDa) is present on immunoblots of mitochondrial fraction from mouse embryonic stem cells with an antibody directed against Tet1. This band, however, is not protected from trypsin treatment of mitochondria indicating that Tet1 may not be transported to the mitochondrial matrix. The putative Tet1 mitochondrial targeting sequence (MTS) fails to carry heterologous protein to the mitochondria. Knock out of Tet1 in mouse ES cells also does not alter 5hmC signal in hydroxyMeDIP assay. We now seek to determine if Tet2/Tet3 may be involved in 5hmC generation. In the nucleus, 5hmC serves as an intermediate in the process of DNA demethylation through the combined action of cytidine deaminases and the base excision repair pathway. We plan to investigate if 5hmC holds the same functional significance in the mitochondria as it does in the nucleus. Our overall goal is to understand epigenetic regulation of normal mitochondrial function and changes that occur in diseases involving mitochondrial dysfunction such as ischemic heart disease, neurodegenerative diseases like Parkinsons disease, and cancer.

Understanding the impact of epigenetic mechanisms on tumorigenesis is essential, as epigenetic alterations are associated with tumor initiation and progression. Because epigenetic changes are reversible, they are potential targets for cancer therapy. Nucleosome Remodeling Factor (NURF) is a chromatin-remodeling complex that regulates gene expression by changing nucleosome positioning along the DNA sequence. Previous studies have shown a role for NURF in embryonic development as well as regulating genes involved in tumor progression. In this work we investigated the impact of eliminating NURF function in tumorigenesis in vivo. BALB/c mice challenged with syngeneic 67NR breast cancer cell lines, injected into the mammary fat pad, lacking NURF, due to knockdown of its essential subunits Bptf, showed reduction in tumor growth comparing to control tumors. The observed reduction in tumor growth was abrogated in immunodeficient mice lacking a functional immune system. Bptf KD and control 67NR cells grew at similar rates in vitro. Similar findings were observed in our lab using 66cl4 breast cancer cell lines. Using immunofluorescence staining, no significant difference in CD8+, CD4+, NK and MDSC cells infiltrations into the tumor microenvironment was observed in 66cl4 tumors. Preliminary results from 67NR tumors suggested more CD4+ and CD8+ cells. Gene expression profile of tumor tissues from BALB/c mice injected with 67NR and 66cl4 cell lines showed enrichment of genes associated with immune response. Our findings suggested a role of the immune system in targeting tumor cells lacking Bptf in vivo.

DNA methylation is the most recognizable epigenetic mechanism. In general, DNA hypomethylation is associated with increased gene expression whereas DNA hypermethylation is associated with decreased gene expression. To date, little is known about the role of DNA methylation in the pathophysiology of preeclampsia. In this study, we examined the differences in DNA methylation in omental arteries of normal pregnant and preeclamptic women using the high throughput Illumina HumanMethylation27 BeadChip assay. We found 1,685 genes with a significant difference in DNA methylation at a false discovery rate of < 10% with many inflammatory genes having reduced methylation. The thromboxane synthase gene was the most hypomethylated gene in preeclamptic women as compared to normal pregnant women. When we examined the expression of thromboxane synthase in omental arteries of normal pregnant and preeclamptic women we found it to be significantly increased in preeclamptic women. The increased expression was observed in vascular smooth muscle cells, endothelial cells and infiltrating neutrophils. Experimentally induced DNA hypomethylation increased the expression of thromboxane synthase in the neutrophil-like HL-60 cell line, whereas tumor necrosis factor α (TNFα), a neutrophil product, increased its expression in cultured human vascular smooth muscle cells (VSMC). These finding suggest that DNA methylation and release of TNFα by infiltrating neutrophils could contribute to the increased expression of thromboxane synthase in systemic blood vessels of preeclamptic women, contributing to the hypertension and coagulation abnormalities. We also explored the possible contribution of DNA methylation to the altered expression of genes involved in collagen metabolism in preeclampsia. Several matrix metalloproteinase (MMP) genes, including MMP1 and MMP8, were significantly less methylated in preeclamptic women, whereas TIMP and COL genes were either significantly more methylated or had no significant change in their DNA methylation status. Experimentally induced DNA hypomethylation increased the expression of MMP-1, but not TIMP-1 or COL1A1, in cultured VSMCs and increased the expression of MMP-1 and MMP-8 in HL-60 cells. These findings suggest that DNA methylation contributes to the imbalance in genes involved in collagen metabolism in blood vessels of preeclamptic women.

Chromodomain Helicase DNA-Binding Protein 4 (CHD4) is an ATPase that alters the phasing of nucleosomes on DNA and has recently been implicated in DNA double stranded break (DSB) repair. Here, we show that depletion of CHD4 in Acute Myeloid Leukemia (AML) blasts induces a global relaxation of chromatin that renders cells more susceptible to DSB formation, while concurrently impeding their repair. Furthermore, CHD4 depletion renders AML blasts more sensitive both in vitro and in vivo to genotoxic agents used in clinical therapy: daunorubicin (DNR) and cytarabine (ara-C). Sensitization to DNR and ara-C is mediated in part by activation of the ATM pathway, which is preliminarily activated by a Tip60-dependent mechanism in response to chromatin relaxation and further activated by genotoxic-agent induced DSBs. This sensitization preferentially affects AML cells, as CHD4 depletion in normal CD34+ hematopoetic progenitors does not increase their susceptibility to DNR or ara-C. Unexpectedly, we found that CHD4 is necessary for maintaining the tumor formatting behavior of AML cells, as CHD4 depletion severely restricted the ability of AML cells to form xenografts in mice and colonies in soft agar. Taken together, these results provide evidence for CHD4 as a novel therapeutic target whose inhibition has the potential to enhance the effectiveness of genotoxic agents used in AML therapy.

The myeloid cell lineage is a fundamental element of the immune system and it can give rise to a diverse set of terminally differentiated cells, such as macrophages or osteoclasts among many others. Mutations or misregulation of gene expression may lead to severe clinical conditions, such as arthritis, osteoporosis or cancers. Epigenetics, the regulation of gene expression and chromatin remodelling, is implicated in cell differentiation, function and disease, and hence it is a promising new area to explore in order to explain underlying cellular mechanisms. Firstly, human macrophage subtypes were studied. Chemokine (C-C motif) ligand (CCL) 1 and mannose receptor were validated to be granulocyte macrophage (GM) colony stimulating factor (CSF) induced macrophage markers, while CCL₂ was specifically expressed in macrophage CSF (MCSF) macrophage population. By utilising publicly available high-throughput sequencing data, new biomarkers dehydrogenase/reductase (SDR family) member 2 and CCL₂₆ were discovered to be MCSF-macrophage specific while guanylate binding protein 5 and apolipoprotein B mRNA editing enzyme, catalytic polypeptide-like 3A were highly up-regulated in GMCSF cells. Secondly, a range of gene knock-down techniques for the myeloid cell lineage were optimised and established. Lentiviral short-hairpin RNA (shRNA) delivery methods were shown to induce an undesirable pro-inflammatory response in macrophages. Furthermore, the frequently utilised cytomegalovirus promoter for gene expression was shown to be completely silenced in macrophage populations. Locked nucleic acids were selected as a suitable alternative to shRNA knock-down and by employing this new tool it was shown that a histone demethylase lysine (K)-specific demethylase (KDM) 6B is fundamental for macrophage differentiation. Finally, a small molecule GSK-J₄, a potent inhibitor of histone demethylases KDM6A, KDM6B and KDM₅B specific for H₃K_27me3 and H₃K_4me3, respectively, was used to dissect epigenetic signalling in osteoclasts and multiple myeloma. In osteoclasts it was shown to act mainly by inhibiting transcriptional changes required for osteoclastogenesis when MCSF-macrophages are stimulated with Receptor Activator Of Nuclear Factor Kappa-B Ligand (RANKL), as indicated by the differential increase in H₃K_27me3 marks, leading to inhibition of c-Jun and potentially abolition of transcription factor AP-1, required for the transcriptional initiation of nuclear factor of activated T-cells 1 (NFATc1). In multiple myeloma cells, GSK-J₄ causes a dramatic increase in expression, further supported by the build-up of global H₃K_4me3 marks, which results in the upregulation of the unfolded protein response pathway. In both cell systems, there is an early upregulation of metallothionein genes, which in multiple myeloma was shown to increase potentially due to rapid influx of zinc ions within the first 30 minutes, and as such may cause induction of apoptosis in multiple myeloma and may inhibit differentiation of osteoclasts.

Tout au long de la vie d’un individu, il existe un nombre optimal de cellules à produire et de progéniteurs à conserver en réserve. On parle de maintien de l’homéostasie tissulaire. De façon générale, l’organisme a cinq possibilités pour réguler l’homéostasie : l’autorenouvellement et la quiescence, souvent utilisés pour maintenir un ‘pool’ fonctionnel de progéniteurs, la différenciation qui permet de produire des cellules effectrices, l’apoptose et la sénescence, qui permettent de limiter la production de cellules ou encore d’en faire diminuer le nombre quand elles sont en excès. La régulation de ces quatre mécanismes peut se faire de façon extrinsèque en passant par différentes voies de signalisation combinées à l’action intrinsèque de facteurs de transcription comme Ikaros et GATA1. Le facteur de transcription Ikaros joue un rôle critique dans le devenir des cellules progénitrices et la différenciation des lignages hématopoïétiques. Cependant, il demeure surtout connu pour son influence sur la voie Notch dans les cellules lymphoïdes, notamment les lymphocytes T. Les cellules érythroïdes sont hautement sensibles à l’environnement et donc, particulièrement adaptées à l’étude des régulations de l’homéostasie. Les résultats de différentes études ont permis de démontrer qu’Ikaros et la voie Notch influencent l’érythropoïèse. Cependant le détail de leurs actions demeure en grande partie inconnu à ce jour. Au cours de notre étude nous avons voulu déterminer l’action d’Ikaros dans le maintien de l’homéostasie des cellules érythroïdes et si son rôle passe par un dialogue avec la voie Notch. Nous avons voulu décrypter les mécanismes de régulation transcriptionnelle utilisés par Ikaros et par Notch au cours de l’érythropoïèse et leurs effets. Notre étude montre qu’Ikaros réprime à l’aide de GATA1 le gène Hes1, une cible importante de la voie Notch, en recrutant un complexe de la famille Polycomb, le PRC2 ii (Polycomb Repressive Complex 2). Cette répression permet la promotion de la différenciation des cellules érythroïdes. Au niveau du maintien de l’homéostasie par régulation de l’apoptose, Ikaros est connu pour cibler l’anti-apoptotique Bcl2l1 dans les lymphocytes. Puisque Gata-1, partenaire préférentiel d’Ikaros cible Bcl2l1 dans les cellules érythroïdes, nous avons caractérisé leur effet sur l’expression de Bcl2l1. Nous avons découvert qu’Ikaros active de façon directe Bcl2l1 et qu’il recrute sur le gène deux complexes partenaires d’élongation : un de la famille SET1/MLL, et le complexe P-TEFb-NuRD. En l’absence d’Ikaros, le fragment intracellulaire de Notch (NICD) et son cofacteur RBP-J remplacent Ikaros et favorisent l’hyper-activation de l’expression de Bcl2l1. Ceci est associé à la modification du complexe d’élongation recruté, ainsi qu’à la mise en place de modifications épigénétiques distinctes de celles observées avec Ikaros ce qui modifie l’élongation transcriptionnelle du gène. Ikaros et Notch sont fréquemment mutés ou présentent des fonctions altérées dans les leucémies. Notre étude montre un dialogue Ikaros/Notch influençant aussi bien la différenciation que l’apoptose et met en évidence l’existence d’un circuit génétique dont le dérèglement pourrait favoriser l’apparition d’une hématopoïèse maligne. / Throughout the life of an individual, there is an optimal count of cells to produce and progenitors to conserve in stock. This is the tissue homeostasis maintenance. In a general fashion the organism has five means to regulate the homeostasis. Self-renewal and quiescence, often used in order to maintain a functional progenitors pool. Differentiation enhances effector cells production. Apoptosis and senescence can limit cell production and reduce cell number in case of excess. These regulation mechanisms can be performed in an extrinsic fashion using different signaling pathways combined with the action of transcription factors like Ikaros and GATA1. The transcription factor Ikaros is critical for progenitor cells fate and hematopoietic lineages differentiation. However, Ikaros is mostly known for its influence on Notch signaling in lymphoid cells, notably T lymphocytes. Erythroid cells are highly sensitive to the environment thus, particularly adapted to study homeostasis maintenance regulation. Results obtained in different studies showed Ikaros and Notch signaling influencing erythropoiesis. However, the detail of their effect remains mainly unknown to day. Our aim was to determine Ikaros effect on erythroid cells homeostasis maintenance and if its role involved a cross-talk with Notch signaling. We will decipher transcription regulation mechanisms used by Ikaros and Notch during erythropoiesis and their effects. We show Ikaros uses GATA1 to repress Hes1, a major Notch target by recruiting a Polycomb family complex, the PRC2 (Polycomb Repressive Complex 2). This repression promotes erythroid cells differentiation. At the apoptosis mediated control of homeostasis level, Ikaros is known to target Bcl2l1 in lymphocytes. As GATA1, Ikaros preferential partner, targets Bcl2l1 in erythroid cells, we assessed their effect on Bcl2l1 expression. We discovered Ikaros directly activates Bcl2l1 iv and recruits two elongation associated complexes: one from SET1/MLL complex family, and the P-TEFb-NuRD complex. In the absence of Ikaros, the intracellular fragment of Notch (NICD) and its cofactor RBP-J replace Ikaros and favors Bcl2l1 overactivation. This is associated with a switch of recruited elongation associated complex and the establishment of distinct epigenetic modifications from those observed with Ikaros, which modifies the gene transcriptional elongation. Ikaros and Notch are frequently mutated or present altered functions in leukemia. Our works present an Ikaros/Notch cross-talk influencing as well differentiation as apoptosis and reveal the existence of a genetic circuit for which a malfunction could favor hematologic disorders. Keywords : transcription, homeostasis, erythropoiesis

La betterave sucrière est une plante bisannuelle dont le besoin de vernalisation est absolu. Ce processus correspond à l’acquisition de l’aptitude à la montaison et à la floraison et résulte d’une exposition prolongée à de basses températures. La durée de froid requise pour induire la montaison puis la floraison varie suivant les génotypes et reflète leur tolérance à la montaison, qui constitue donc un caractère agronomique essentiel. Cette thèse visait à (i) mettre en évidence un éventuel contrôle épigénétique (méthylation ADN) de l‘induction de la montaison chez des génotypes de betterave sucrière résistants ou sensibles à la montaison, (ii) identifier les séquences ciblées par des remaniements de méthylation de l’ADN et d’expression associés, et (iii) caractériser certaines séquences candidates en vue de leur utilisation comme marqueurs de la montaison. Nos travaux ont montré que l’amplitude et la cinétique des variations de méthylation de l’ADN observées au cours de la vernalisation semblent être des éléments critiques de l’induction et de la tolérance à la montaison. Par une approche ciblée, des séquences dont la méthylation de l’ADN est remaniée ont été identifiées. L’implication dans la transition florale de deux BvRNMT (RNA METHYLTRANSFERASES) et de la méthylation des ARN, tels que l’ARNm de BvFL1, un répresseur floral, a ainsi été mise en évidence chez la betterave sucrière. Enfin, grâce à une approche génomique, un réseau de gènes intégratif incluant la réponse à l’environnement, la signalisation hormonale et l’induction de la floraison a été identifié. La cinétique d’activation de ces gènes définirait le niveau de tolérance à la montaison des différents génotypes de betterave sucrière. / Sugar beet is a biennial plant with an absolute requirement of vernalization, corresponding to the acquisition of the competence to bolt and flower after a prolonged exposure to low temperatures. The cold duration needed to induce bolting and flowering varies depending on the genotypes, reflecting their bolting tolerance, which is an essential agronomic trait. This work aimed at (i) investigating a possible epigenetic control of bolting induction in bolting sensitive and bolting resistant sugar beet genotypes, (ii) identifying sequences targeted by DNA methylation and expression remodeling, and (iii) characterizing candidate sequences which could be used in marked-assisted selection for plant breeding. Our data suggest that the time course and amplitude of DNA methylation variations are critical points for the induction of sugar beet bolting and represent an epigenetic component of the genotypic bolting tolerance. In addition, we identified differentially methylated sequences exhibiting variations of gene-body DNA methylation and expression during cold exposure and/or between genotypes. Among them, two RNA METHYLCYTOSINE TRANSFERASES, in association with RNA methylation such as BvFL1 mRNA, a floral repressor, were shown to play a role in floral transition. Finally, using microarrays we identified an integrative network of genes including response to environment, phytohormone signalling and flowering induction. The activation kinetics of these genes could define the bolting tolerance level of sugar beet genotypes.

Obesity is widely considered a primary risk factor for type 2 diabetes (T2D). However, less is known about the early adaptive responses to short-term periods of high-fat energy excess (HFEE). Previous reports detailing whole-body adaptation to fat and energy oversupply are equivocal, perhaps, in part, owing to use of different experimental protocols, varying durations of dietary manipulation and participant cohorts with individuals of varying characteristics. In addition to use of different dietary protocols between studies, alterations in functional end-point measures due to the type of dietary fat consumed warrants consideration. Daily n-3 PUFA intake, commonly obtained from pelagic fish oil (FO) consumption, has been shown to positively associate with insulin sensitivity in epidemiological studies and thus may be a useful dietary strategy for slowing insulin resistance development. Chapter 2 of this thesis extends previous literature by demonstrating that 6 d HFEE (150 % habitual energy intake; 60 % of energy from fat) does not clearly alter whole- body insulin sensitivity, irrespective of FO consumption. However, investigation of metabolism at the tissue level, as presented in Chapter 3 of this thesis, offers insight into a potential tissue-specific level of regulation that precedes whole-body regulation. Skeletal muscle insulin signalling protein (e.g. protein kinase B (PKB)) activity, levels of certain ceramide species, and AMPK α2 activity were altered following HFEE and may explain the early maladaptive responses to short-term HFEE. Moreover, FO intake as 10 % of total fats mediated some of these molecular responses, including PKB and AMPK α2 activity, reflecting possible functional effects of FO at the subcellular level. Regulation of these metabolic / molecular responses at both the tissue and whole- body level can be explained, in part, by genetic predisposition, environmental influence and more recently epigenetics, including microRNAs (miRNAs). In Chapter 4, we characterised the plasma and skeletal muscle miRNA responses to HFEE and oral glucose ingestion. We demonstrate transient changes in levels of certain miRNAs following oral glucose ingestion in both tissue types and in response to HFEE in skeletal muscle. However, no significant correlations between basal plasma and skeletal muscle miRNA levels were observed, suggesting that our candidate plasma miRNAs may be co-ordinating functional changes in other tissue types. Plasma miR- 145-5p and skeletal muscle miR-204-5p predicted a significant proportion of the variance in mean whole-body insulin sensitivity change in response to HFEE. These data indicate that these miRNAs may be useful biomarkers of insulin resistance development following HFEE. A constraint of this thesis is that all conclusions are made within the context of statistically unaltered insulin sensitivity. Therefore, future investigations of diet-induced maladaptation should consider establishing a time course of insulin resistance development in response to HFEE, or use different study populations. Populations that are more susceptible to T2D development, e.g., overweight, sedentary individuals would be of particular interest. These data would aid development of a working model of diet-induced insulin resistance that has more direct application to T2D progression and extends the data presented herein.