CRISPR/dCAS9-mediated DNA demethylation screen identifies functional epigenetic determinants of colorectal cancer.

BACKGROUND: Promoter hypermethylation of tumour suppressor genes is frequently observed during the malignant transformation of colorectal cancer (CRC). However, whether this epigenetic mechanism is functional in cancer or is a mere consequence of the carcinogenic process remains to be elucidated. RESULTS: In this work, we performed an integrative multi-omic approach to identify gene candidates with strong correlations between DNA methylation and gene expression in human CRC samples and a set of 8 colon cancer cell lines. As a proof of concept, we combined recent CRISPR-Cas9 epigenome editing tools (dCas9-TET1, dCas9-TET-IM) with a customized arrayed gRNA library to modulate the DNA methylation status of 56 promoters previously linked with strong epigenetic repression in CRC, and we monitored the potential functional consequences of this DNA methylation loss by means of a high-content cell proliferation screen. Overall, the epigenetic modulation of most of these DNA methylated regions had a mild impact on the reactivation of gene expression and on the viability of cancer cells. Interestingly, we found that epigenetic reactivation of RSPO2 in the tumour context was associated with a significant impairment in cell proliferation in p53-/- cancer cell lines, and further validation with human samples demonstrated that the epigenetic silencing of RSPO2 is a mid-late event in the adenoma to carcinoma sequence. CONCLUSIONS: These results highlight the potential role of DNA methylation as a driver mechanism of CRC and paves the way for the identification of novel therapeutic windows based on the epigenetic reactivation of certain tumour suppressor genes.

Nicotinamide N-methyltransferase (NNMT) regulates the glucocorticoid signaling pathway during the early phase of adipogenesis.

Obesity is associated with adipose tissue dysfunction through the differentiation and expansion of pre-adipocytes to adipocytes (hyperplasia) and/or increases in size of pre-existing adipocytes (hypertrophy). A cascade of transcriptional events coordinates the differentiation of pre-adipocytes into fully differentiated adipocytes; the process of adipogenesis. Although nicotinamide N-methyltransferase (NNMT) has been associated with obesity, how NNMT is regulated during adipogenesis, and the underlying regulatory mechanisms, remain undefined. In present study we used genetic and pharmacological approaches to elucidate the molecular signals driving NNMT activation and its role during adipogenesis. Firstly, we demonstrated that during the early phase of adipocyte differentiation NNMT is transactivated by CCAAT/Enhancer Binding Protein beta (CEBPB) in response to glucocorticoid (GC) induction. We found that Nnmt knockout, using CRISPR/Cas9 approach, impaired terminal adipogenesis by influencing the timing of cellular commitment and cell cycle exit during mitotic clonal expansion, as demonstrated by cell cycle analysis and RNA sequencing experiments. Biochemical and computational methods showed that a novel small molecule, called CC-410, stably binds to and highly specifically inhibits NNMT. CC-410 was, therefore, used to modulate protein activity during pre-adipocyte differentiation stages, demonstrating that, in line with the genetic approach, chemical inhibition of NNMT at the early stages of adipogenesis impairs terminal differentiation by deregulating the GC network. These congruent results conclusively demonstrate that NNMT is a key component of the GC-CEBP axis during the early stages of adipogenesis and could be a potential therapeutic target for both early-onset obesity and glucocorticoid-induced obesity.

The outcome of boosting mitochondrial activity in alcohol-associated liver disease is organ-dependent.

BACKGROUND AND AIMS: Alcohol-associated liver disease (ALD) accounts for 70% of liver-related deaths in Europe, with no effective approved therapies. Although mitochondrial dysfunction is one of the earliest manifestations of alcohol-induced injury, restoring mitochondrial activity remains a problematic strategy due to oxidative stress. Here, we identify methylation-controlled J protein (MCJ) as a mediator for ALD progression and hypothesize that targeting MCJ may help in recovering mitochondrial fitness without collateral oxidative damage. APPROACH AND RESULTS: C57BL/6 mice [wild-type (Wt)] Mcj knockout and Mcj liver-specific silencing (MCJ-LSS) underwent the NIAAA dietary protocol (Lieber-DeCarli diet containing 5% (vol/vol) ethanol for 10 days, plus a single binge ethanol feeding at day 11). To evaluate the impact of a restored mitochondrial activity in ALD, the liver, gut, and pancreas were characterized, focusing on lipid metabolism, glucose homeostasis, intestinal permeability, and microbiota composition. MCJ, a protein acting as an endogenous negative regulator of mitochondrial respiration, is downregulated in the early stages of ALD and increases with the severity of the disease. Whole-body deficiency of MCJ is detrimental during ALD because it exacerbates the systemic effects of alcohol abuse through altered intestinal permeability, increased endotoxemia, and dysregulation of pancreatic function, which overall worsens liver injury. On the other hand, liver-specific Mcj silencing prevents main ALD hallmarks, that is, mitochondrial dysfunction, steatosis, inflammation, and oxidative stress, as it restores the NAD + /NADH ratio and SIRT1 function, hence preventing de novo lipogenesis and improving lipid oxidation. CONCLUSIONS: Improving mitochondrial respiration by liver-specific Mcj silencing might become a novel therapeutic approach for treating ALD.

Multi-omic rejuvenation of naturally aged tissues by a single cycle of transient reprogramming.

The expression of the pluripotency factors OCT4, SOX2, KLF4, and MYC (OSKM) can convert somatic differentiated cells into pluripotent stem cells in a process known as reprogramming. Notably, partial and reversible reprogramming does not change cell identity but can reverse markers of aging in cells, improve the capacity of aged mice to repair tissue injuries, and extend longevity in progeroid mice. However, little is known about the mechanisms involved. Here, we have studied changes in the DNA methylome, transcriptome, and metabolome in naturally aged mice subject to a single period of transient OSKM expression. We found that this is sufficient to reverse DNA methylation changes that occur upon aging in the pancreas, liver, spleen, and blood. Similarly, we observed reversion of transcriptional changes, especially regarding biological processes known to change during aging. Finally, some serum metabolites and biomarkers altered with aging were also restored to young levels upon transient reprogramming. These observations indicate that a single period of OSKM expression can drive epigenetic, transcriptomic, and metabolomic changes toward a younger configuration in multiple tissues and in the serum.

Epigenetic Profiling and Response to CD19 Chimeric Antigen Receptor T-Cell Therapy in B-Cell Malignancies.

BACKGROUND: Chimeric antigen receptor (CAR) T cells directed against CD19 (CART19) are effective in B-cell malignancies, but little is known about the molecular factors predicting clinical outcome of CART19 therapy. The increasingly recognized relevance of epigenetic changes in cancer immunology prompted us to determine the impact of the DNA methylation profiles of CART19 cells on the clinical course. METHODS: We recruited 114 patients with B-cell malignancies, comprising 77 patients with acute lymphoblastic leukemia and 37 patients with non-Hodgkin lymphoma who were treated with CART19 cells. Using a comprehensive DNA methylation microarray, we determined the epigenomic changes that occur in the patient T cells upon transduction of the CAR vector. The effects of the identified DNA methylation sites on clinical response, cytokine release syndrome, immune effector cell-associated neurotoxicity syndrome, event-free survival, and overall survival were assessed. All statistical tests were 2-sided. RESULTS: We identified 984 genomic sites with differential DNA methylation between CAR-untransduced and CAR-transduced T cells before infusion into the patient. Eighteen of these distinct epigenetic loci were associated with complete response (CR), adjusting by multiple testing. Using the sites linked to CR, an epigenetic signature, referred to hereafter as the EPICART signature, was established in the initial discovery cohort (n = 79), which was associated with CR (Fisher exact test, P < .001) and enhanced event-free survival (hazard ratio [HR] = 0.36; 95% confidence interval [CI] = 0.19 to 0.70; P = .002; log-rank P = .003) and overall survival (HR = 0.45; 95% CI = 0.20 to 0.99; P = .047; log-rank P = .04;). Most important, the EPICART profile maintained its clinical course predictive value in the validation cohort (n = 35), where it was associated with CR (Fisher exact test, P < .001) and enhanced overall survival (HR = 0.31; 95% CI = 0.11 to 0.84; P = .02; log-rank P = .02). CONCLUSIONS: We show that the DNA methylation landscape of patient CART19 cells influences the efficacy of the cellular immunotherapy treatment in patients with B-cell malignancy.

Physical exercise shapes the mouse brain epigenome.

OBJECTIVE: To analyze the genome-wide epigenomic and transcriptomic changes induced by long term resistance or endurance training in the hippocampus of wild-type mice. METHODS: We performed whole-genome bisulfite sequencing (WGBS) and RNA sequencing (RNA-seq) of mice hippocampus after 4 weeks of specific training. In addition, we used a novel object recognition test before and after the intervention to determine whether the exercise led to an improvement in cognitive function. RESULTS: Although the majority of DNA methylation changes identified in this study were training-model specific, most were associated with hypomethylation and were enriched in similar histone marks, chromatin states, and transcription factor biding sites. It is worth highlighting the significant association found between the loss of DNA methylation in Tet1 binding sites and gene expression changes, indicating the importance of these epigenomic changes in transcriptional regulation. However, endurance and resistance training activate different gene pathways, those being associated with neuroplasticity in the case of endurance exercise, and interferon response pathways in the case of resistance exercise, which also appears to be associated with improved learning and memory functions. CONCLUSIONS: Our results help both understand the molecular mechanisms by which different exercise models exert beneficial effects for brain health and provide new potential therapeutic targets for future research.

Epigenetic Deregulation of the Histone Methyltransferase KMT5B Contributes to Malignant Transformation in Glioblastoma.

Glioblastoma multiforme (GBM) is the most common and aggressive type of brain tumor in adulthood. Epigenetic mechanisms are known to play a key role in GBM although the involvement of histone methyltransferase KMT5B and its mark H4K20me2 has remained largely unexplored. The present study shows that DNA hypermethylation and loss of DNA hydroxymethylation is associated with KMT5B downregulation and genome-wide reduction of H4K20me2 levels in a set of human GBM samples and cell lines as compared with non-tumoral specimens. Ectopic overexpression of KMT5B induced tumor suppressor-like features in vitro and in a mouse tumor xenograft model, as well as changes in the expression of several glioblastoma-related genes. H4K20me2 enrichment was found immediately upstream of the promoter regions of a subset of deregulated genes, thus suggesting a possible role for KMT5B in GBM through the epigenetic modulation of key target cancer genes.

Stability of Imprinting and Differentiation Capacity in Naïve Human Cells Induced by Chemical Inhibition of CDK8 and CDK19.

Pluripotent stem cells can be stabilized in vitro at different developmental states by the use of specific chemicals and soluble factors. The naïve and primed states are the best characterized pluripotency states. Naïve pluripotent stem cells (PSCs) correspond to the early pre-implantation blastocyst and, in mice, constitute the optimal starting state for subsequent developmental applications. However, the stabilization of human naïve PSCs remains challenging because, after short-term culture, most current methods result in karyotypic abnormalities, aberrant DNA methylation patterns, loss of imprinting and severely compromised developmental potency. We have recently developed a novel method to induce and stabilize naïve human PSCs that consists in the simple addition of a chemical inhibitor for the closely related CDK8 and CDK19 kinases (CDK8/19i). Long-term cultured CDK8/19i-naïve human PSCs preserve their normal karyotype and do not show widespread DNA demethylation. Here, we investigate the long-term stability of allele-specific methylation at imprinted loci and the differentiation potency of CDK8/19i-naïve human PSCs. We report that long-term cultured CDK8/19i-naïve human PSCs retain the imprinting profile of their parental primed cells, and imprints are further retained upon differentiation in the context of teratoma formation. We have also tested the capacity of long-term cultured CDK8/19i-naïve human PSCs to differentiate into primordial germ cell (PGC)-like cells (PGCLCs) and trophoblast stem cells (TSCs), two cell types that are accessible from the naïve state. Interestingly, long-term cultured CDK8/19i-naïve human PSCs differentiated into PGCLCs with a similar efficiency to their primed counterparts. Also, long-term cultured CDK8/19i-naïve human PSCs were able to differentiate into TSCs, a transition that was not possible for primed PSCs. We conclude that inhibition of CDK8/19 stabilizes human PSCs in a functional naïve state that preserves imprinting and potency over long-term culture.

Conservation of Aging and Cancer Epigenetic Signatures across Human and Mouse.

Aging and cancer are two interrelated processes, with aging being a major risk factor for the development of cancer. Parallel epigenetic alterations have been described for both, although differences, especially within the DNA hypomethylation scenario, have also been recently reported. Although many of these observations arise from the use of mouse models, there is a lack of systematic comparisons of human and mouse epigenetic patterns in the context of disease. However, such comparisons are significant as they allow to establish the extent to which some of the observed similarities or differences arise from pre-existing species-specific epigenetic traits. Here, we have used reduced representation bisulfite sequencing to profile the brain methylomes of young and old, tumoral and nontumoral brain samples from human and mouse. We first characterized the baseline epigenomic patterns of the species and subsequently focused on the DNA methylation alterations associated with cancer and aging. Next, we described the functional genomic and epigenomic context associated with the alterations, and finally, we integrated our data to study interspecies DNA methylation levels at orthologous CpG sites. Globally, we found considerable differences between the characteristics of DNA methylation alterations in cancer and aging in both species. Moreover, we describe robust evidence for the conservation of the specific cancer and aging epigenomic signatures in human and mouse. Our observations point toward the preservation of the functional consequences of these alterations at multiple levels of genomic regulation. Finally, our analyses reveal a role for the genomic context in explaining disease- and species-specific epigenetic traits.

Epigenetic downregulation of TET3 reduces genome-wide 5hmC levels and promotes glioblastoma tumorigenesis.

Loss of 5-hydroxymethylcytosine (5hmC) has been associated with mutations of the ten-eleven translocation (TET) enzymes in several types of cancer. However, tumors with wild-type TET genes can also display low 5hmC levels, suggesting that other mechanisms involved in gene regulation might be implicated in the decline of this epigenetic mark. Here we show that DNA hypermethylation and loss of DNA hydroxymethylation, as well as a marked reduction of activating histone marks in the TET3 gene, impair TET3 expression and lead to a genome-wide reduction in 5hmC levels in glioma samples and cancer cell lines. Epigenetic drugs increased expression of TET3 in glioblastoma cells and ectopic overexpression of TET3 impaired in vitro cell growth and markedly reduced tumor formation in immunodeficient mice models. TET3 overexpression partially restored the genome-wide patterns of 5hmC characteristic of control brain samples in glioblastoma cell lines, while elevated TET3 mRNA levels were correlated with better prognosis in glioma samples. Our results suggest that epigenetic repression of TET3 might promote glioblastoma tumorigenesis through the genome-wide alteration of 5hmC.

Chromatin regulation by Histone H4 acetylation at Lysine 16 during cell death and differentiation in the myeloid compartment.

Histone H4 acetylation at Lysine 16 (H4K16ac) is a key epigenetic mark involved in gene regulation, DNA repair and chromatin remodeling, and though it is known to be essential for embryonic development, its role during adult life is still poorly understood. Here we show that this lysine is massively hyperacetylated in peripheral neutrophils. Genome-wide mapping of H4K16ac in terminally differentiated blood cells, along with functional experiments, supported a role for this histone post-translational modification in the regulation of cell differentiation and apoptosis in the hematopoietic system. Furthermore, in neutrophils, H4K16ac was enriched at specific DNA repeats. These DNA regions presented an accessible chromatin conformation and were associated with the cleavage sites that generate the 50 kb DNA fragments during the first stages of programmed cell death. Our results thus suggest that H4K16ac plays a dual role in myeloid cells as it not only regulates differentiation and apoptosis, but it also exhibits a non-canonical structural role in poising chromatin for cleavage at an early stage of neutrophil cell death.

Longitudinal genome-wide DNA methylation analysis uncovers persistent early-life DNA methylation changes.

BACKGROUND: Early life is a period of drastic epigenetic remodeling in which the epigenome is especially sensitive to extrinsic and intrinsic influence. However, the epigenome-wide dynamics of the DNA methylation changes that occur during this period have not been sufficiently characterized in longitudinal studies. METHODS: To this end, we studied the DNA methylation status of more than 750,000 CpG sites using Illumina MethylationEPIC arrays on 33 paired blood samples from 11 subjects at birth and at 5 and 10 years of age, then characterized the chromatin context associated with these loci by integrating our data with histone, chromatin-state and enhancer-element external datasets, and, finally, validated our results through bisulfite pyrosequencing in two independent longitudinal cohorts of 18 additional subjects. RESULTS: We found abundant DNA methylation changes (110,726 CpG sites) during the first lustrum of life, while far fewer alterations were observed in the subsequent 5 years (460 CpG sites). However, our analysis revealed persistent DNA methylation changes at 240 CpG sites, indicating that there are genomic locations of considerable epigenetic change beyond immediate birth. The chromatin context of hypermethylation changes was associated with repressive genomic locations and genes with developmental and cell signaling functions, while hypomethylation changes were linked to enhancer regions and genes with immunological and mRNA and protein metabolism functions. Significantly, our results show that genes that suffer simultaneous hyper- and hypomethylation are functionally distinct from exclusively hyper- or hypomethylated genes, and that enhancer-associated methylation is different in hyper- and hypomethylation scenarios, with hypomethylation being more associated to epigenetic changes at blood tissue-specific enhancer elements. CONCLUSIONS: These data show that epigenetic remodeling is dramatically reduced after the first 5 years of life. However, there are certain loci which continue to manifest DNA methylation changes, pointing towards a possible functionality beyond early development. Furthermore, our results deepen the understanding of the genomic context associated to hyper- or hypomethylation alterations during time, suggesting that hypomethylation of blood tissue-specific enhancer elements could be of importance in the establishment of functional states in blood tissue during early-life.

Deregulation of the imprinted DLK1-DIO3 locus ncRNAs is associated with replicative senescence of human adipose-derived stem cells.

BACKGROUND: Human adult adipose-derived stem cells (hADSCs) have become the most promising cell source for regenerative medicine. However the prolonged ex vivo expansion periods required to obtain the necessary therapeutic dose promotes progressive senescence, with the concomitant reduction of their therapeutic potential. AIM AND SCOPE: A better understanding of the determinants of hADSC senescence is needed to improve biosafety while preserving therapeutic efficiency. Here, we investigated the association between deregulation of the imprinted DLK1-DIO3 region and replicative senescence in hADSC cultures. METHODS: We compared hADSC cultures at short (PS) and prolonged (PL) passages, both in standard and low [O2] (21 and 3%, respectively), in relation to replicative senescence. hADSCs were evaluated for expression alterations in the DLK1-DIO3 region on chromosome 14q32, and particularly in its main miRNA cluster. RESULTS: Comparison of hADSCs cultured at PL or PS surprisingly showed a quite significant fraction (69%) of upregulated miRNAs in PL cultures mapping to the imprinted 14q32 locus, the largest miRNA cluster described in the genome. In agreement, expression of the lncRNA MEG3 (Maternally Expressed 3; Meg3/Gtl2), cultured at 21 and 3% [O2], was also significantly higher in PL than in PS passages. During hADSC replicative senescence the AcK16H4 activating mark was found to be significantly associated with the deregulation of the entire DLK1-DIO3 locus, with a secondary regulatory role for the methylation of DMR regions. CONCLUSION: A direct relationship between DLK1-DIO3 deregulation and replicative senescence of hADSCs is reported, involving upregulation of a very significant fraction of its largest miRNA cluster (14q32.31), paralleled by the progressive overexpression of the lncRNA MEG3, which plays a central role in the regulation of Dlk1/Dio3 activation status in mice.

Epigenetic dysregulation of TET2 in human glioblastoma.

Ten-eleven translocation (TET) enzymes are frequently deregulated in cancer, but the underlying molecular mechanisms are still poorly understood. Here we report that TET2 shows frequent epigenetic alterations in human glioblastoma including DNA hypermethylation and hypo-hydroxymethylation, as well as loss of histone acetylation. Ectopic overexpression of TET2 regulated neural differentiation in glioblastoma cell lines and impaired tumor growth. Our results suggest that epigenetic dysregulation of TET2 plays a role in human glioblastoma.

Loss of 5hmC identifies a new type of aberrant DNA hypermethylation in glioma.

Aberrant DNA hypermethylation is a hallmark of cancer although the underlying molecular mechanisms are still poorly understood. To study the possible role of 5-hydroxymethylcytosine (5hmC) in this process we analyzed the global and locus-specific genome-wide levels of 5hmC and 5-methylcytosine (5mC) in human primary samples from 12 non-tumoral brains and 53 gliomas. We found that the levels of 5hmC identified in non-tumoral samples were significantly reduced in gliomas. Strikingly, hypo-hydroxymethylation at 4627 (9.3%) CpG sites was associated with aberrant DNA hypermethylation and was strongly enriched in CpG island shores. The DNA regions containing these CpG sites were enriched in H3K4me2 and presented a different genuine chromatin signature to that characteristic of the genes classically aberrantly hypermethylated in cancer. As this 5mC gain is inversely correlated with loss of 5hmC and has not been identified with classical sodium bisulfite-based technologies, we conclude that our data identifies a novel 5hmC-dependent type of aberrant DNA hypermethylation in glioma.

Epigenome-wide analysis reveals specific DNA hypermethylation of T cells during human hematopoietic differentiation.

AIM: Epigenetic regulation plays an important role in cellular development and differentiation. A detailed map of the DNA methylation dynamics that occur during cell differentiation would contribute to decipher the molecular networks governing cell fate commitment. METHODS: Illumina MethylationEPIC BeadChip platform was used to describe the genome-wide DNA methylation changes observed throughout hematopoietic maturation by analyzing multiple myeloid and lymphoid hematopoietic cell types. RESULTS: We identified a plethora of DNA methylation changes that occur during human hematopoietic differentiation. We observed that T lymphocytes display substantial enhancement of de novo CpG hypermethylation as compared with other hematopoietic cell populations. T-cell-specific hypermethylated regions were strongly associated with open chromatin marks and enhancer elements, as well as binding sites of specific key transcription factors involved in hematopoietic differentiation, such as PU.1 and TAL1. CONCLUSION: These results provide novel insights into the role of DNA methylation at enhancer elements in T-cell development.

DNA methylation changes in human lung epithelia cells exposed to multi-walled carbon nanotubes.

Humans are increasingly exposed to nanoparticles and, although many of their physiological effects have been described, the molecular mechanisms underlying them are still largely unknown. The present study aimed to determine the possible role of certain epigenetic mechanisms in the cellular response of human lung epithelial cells that are triggered by long-term exposure to titanium dioxide nanoparticles (TiO2NPs) and multi-walled carbon nanotubes (MWCNTs). The results showed that exposure to TiO2NPs had only minor effects on genome-wide DNA methylation. However, we identified 755 CpG sites showing consistent DNA hypomethylation in cells exposed to MWCNTs. These sites were mainly located at low density CpG regions and enhancers, and very frequently on the X chromosome. Our results thus suggest that long-term MWCNT exposure may have important effects on the epigenome.

Age-associated hydroxymethylation in human bone-marrow mesenchymal stem cells.

BACKGROUND: Age-associated changes in genomic DNA methylation have been primarily attributed to 5-methylcytosine (5mC). However, the recent discovery of 5-hydroxymethylcytosine (5hmC) suggests that this epigenetic mark might also play a role in the process. METHODS: Here, we analyzed the genome-wide profile of 5hmc in mesenchymal stem cells (MSCs) obtained from bone-marrow donors, aged 2-89 years. RESULTS: We identified 10,685 frequently hydroxymethylated CpG sites in MSCs that were, as in other cell types, significantly associated with low density CpG regions, introns, the histone posttranslational modification H3k4me1 and enhancers. Study of the age-associated changes to 5hmC identified 785 hyper- and 846 hypo-hydroxymethylated CpG sites in the MSCs obtained from older individuals. CONCLUSIONS: DNA hyper-hydroxymethylation in the advanced-age group was associated with loss of 5mC, which suggests that, at specific CpG sites, this epigenetic modification might play a role in DNA methylation changes during lifetime. Since bone-marrow MSCs have many clinical applications, and the fact that the epigenomic alterations in this cell type associated with aging identified in this study could have associated functional effects, the age of donors should be taken into account in clinical settings.

Longitudinal study of DNA methylation during the first 5 years of life.

BACKGROUND: Early life epigenetic programming influences adult health outcomes. Moreover, DNA methylation levels have been found to change more rapidly during the first years of life. Our aim was the identification and characterization of the CpG sites that are modified with time during the first years of life. We hypothesize that these DNA methylation changes would lead to the detection of genes that might be epigenetically modulated by environmental factors during early childhood and which, if disturbed, might contribute to susceptibility to diseases later in life. METHODS: The study of the DNA methylation pattern of 485577 CpG sites was performed on 30 blood samples from 15 subjects, collected both at birth and at 5 years old, using Illumina(®) Infinium 450 k array. To identify differentially methylated CpG (dmCpG) sites, the methylation status of each probe was examined using linear models and the Empirical Bayes Moderated t test implemented in the limma package of R/Bioconductor. Surogate variable analysis was used to account for batch effects. RESULTS: DNA methylation levels significantly changed from birth to 5 years of age in 6641 CpG sites. Of these, 36.79 % were hypermethylated and were associated with genes related mainly to developmental ontology terms, while 63.21 % were hypomethylated probes and associated with genes related to immune function. CONCLUSIONS: Our results suggest that DNA methylation alterations with age during the first years of life might play a significant role in development and the regulation of leukocyte-specific functions. This supports the idea that blood leukocytes experience genome remodeling related to their interaction with environmental factors, underlining the importance of environmental exposures during the first years of life and suggesting that new strategies should be take into consideration for disease prevention.