Figures
Abstract
Introduction and Results
Immunodeficiency, centromeric instability and facial anomalies syndrome (ICF) is a rare autosomal recessive disease, characterized by severe hypomethylation in pericentromeric regions of chromosomes (1, 16 and 9), marked immunodeficiency and facial anomalies. The majority of ICF patients present mutations in the DNMT3B gene, affecting the DNA methyltransferase activity of the protein. In the present study, we have used the Infinium 450K DNA methylation array to evaluate the methylation level of 450,000 CpGs in lymphoblastoid cell lines and untrasformed fibroblasts derived from ICF patients and healthy donors. Our results demonstrate that ICF-specific DNMT3B variants A603T/STP807ins and V699G/R54X cause global DNA hypomethylation compared to wild-type protein. We identified 181 novel differentially methylated positions (DMPs) including subtelomeric and intrachromosomic regions, outside the classical ICF-related pericentromeric hypomethylated positions. Interestingly, these sites were mainly located in intergenic regions and inside the CpG islands. Among the identified hypomethylated CpG-island associated genes, we confirmed the overexpression of three selected genes, BOLL, SYCP2 and NCRNA00221, in ICF compared to healthy controls, which are supposed to be expressed in germ line and silenced in somatic tissues.
Conclusions
In conclusion, this study contributes in clarifying the direct relationship between DNA methylation defect and gene expression impairment in ICF syndrome, identifying novel direct target genes of DNMT3B. A high percentage of the DMPs are located in the subtelomeric regions, indicating a specific role of DNMT3B in methylating these chromosomal sites. Therefore, we provide further evidence that hypomethylation in specific non-pericentromeric regions of chromosomes might be involved in the molecular pathogenesis of ICF syndrome. The detection of DNA hypomethylation at BOLL, SYCP2 and NCRNA00221 may pave the way for the development of specific clinical biomarkers with the aim to facilitate the identification of ICF patients.
Citation: Simo-Riudalbas L, Diaz-Lagares A, Gatto S, Gagliardi M, Crujeiras AB, Matarazzo MR, et al. (2015) Genome-Wide DNA Methylation Analysis Identifies Novel Hypomethylated Non-Pericentromeric Genes with Potential Clinical Implications in ICF Syndrome. PLoS ONE 10(7): e0132517. https://doi.org/10.1371/journal.pone.0132517
Editor: Robert Dante, Institut National de la Santé et de la Recherche Médicale, FRANCE
Received: December 16, 2014; Accepted: June 15, 2015; Published: July 10, 2015
Copyright: © 2015 Simo-Riudalbas et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited
Data Availability: The data reported in this manuscript has been deposited in the Gene Expression Omnibus (GEO) database, www.ncbi.nlm.nih.gov/geo [accession nos. GSE30870] (DNA methylation microarray data from peripheral mononuclear cells obtained from cord blood samples (CB10, CB13, CB20 and CB76) and the following link for fibroblasts and lymphoblastoid cell lines http://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?token=ahcrykaunjaxpcn&acc=GSE68747.
Funding: The research leading to these results has received funding from the Cellex Foundation, Instituto de Sanidad Carlos III (ISCIII) and the Health and Science Departments of the Generalitat de Catalunya. ADL and ABC were supported by the “Sara Borrell” postdoctoral contract (CD12/00738) and (C09/00365) from the ISCIII at the Spanish Ministry of Economy and Competitiveness. MRM is supported by the Epigenomics Flagship Project EPIGEN, MIUR-CNR. ME is an ICREA Research Professor. JS is a “Miguel Servet” researcher and this work was supported by CP13/00055 ISCIII grant. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.
Competing interests: The authors have declared that no competing interests exist.
Introduction
The immunodeficiency, centromeric instability and facial anomalies syndrome (ICF) is a rare recessive disorder, with less than 60 cases reported worldwide. ICF syndrome is characterized by two peculiar signs: a variable immunodeficiency and a recurrent instability of pericentromeric heterochromatin, which usually leads to chromosome breakage in mitogen-stimulated lymphocytes. The chromosomal abnormalities are found exclusively in hypomethylated pericentromeric regions of chromosome 1, 16 and less frequently in 9. Other ICF symptoms count in facial anomalies, psychomotor and mental retardation and developmental delay [1].
The importance of ICF pathology, at the molecular level, relies on the fact that it is the only human disease showing mendelian inheritance of aberrant DNA methylation, caused by mutations in one of the three main DNA-methyltransferase genes, DNMT3B. Approximately 50% of the ICF cases, defined as ICF type1, present biallelic DNMT3B mutations located mainly in the catalytic domain of the protein, often leading to the impairment of its methyltransferase activity [2]. Among the rest of patients some carry nonsense mutations in zinc-finger and BTB domain-containing 24 gene (ZBTB24), designated as ICF2 patients, while a small group of them has still unknown etiology, and are designated as ICFX [2].
The biochemical defects in DNMT3B-mediated de novo DNA methylation have been recently assessed by in vitro studies of the ICF-associated DNMT3B variants [3]. These results reveal that catalysis by DNMT3B is much more complex than expected. In that context, ICF mutations cause a broad spectrum of biochemical defects in DNMT3B function, including defects in homo-oligomerization, SAM binding, SAM utilization and DNA binding [3].
Although it is well established in literature that all ICF1 patient derived cells exhibit targeted DNA hypomethylation of the pericentromeric heterochromatin at chromosomes 1, 16 and sometimes 9, the molecular defects in global DNA methylation caused by ICF-specific DNMT3B mutants remain relatively uncharacterized. Subtle differences in local DNA methylation patterns have been undetected by less sensitive assays previously employed, such as two-dimensional gel electrophoresis of methylation-sensitive digested genomic DNA and COBRA analysis [4],[5]. Recently, several studies suggested that DNA hypomethylation might be more extended than previously thought [6][7]. In this sense, previous genome-wide work from our laboratory, detected DNA hypomethylation in the inactive heterochromatic regions, satellite repeats and transposons associated to two heterozygous DNMT3B mutations in one ICF1 patient [8].
In this study, a state-of-the-art DNA-methylation high-resolution microarray (HumanMethylation450k BeadChip) from Illumina (San Diego, CA) that interrogates 450,000 CpGs sites in the human genome was used to find out the impact of four different ICF-specific DNMT3B mutant alleles (A603T/STP807ins and V699G/R54X) on DNA methylation at a global level. This will help us elucidate and understand the potential pivotal role of this heterogeneous epigenetic mechanism in ICF syndrome pathogenesis. The identification of potential biomarkers for ICF patients after validation in peripheral blood samples of ICF patients might allow us to design more effective strategies to address the diagnosis or treatment of this disease.
Materials and Methods
Sample preparation
The current analysis was performed by evaluating one lymphoblastoid (LCL) cell line and untrasformed fibroblasts derived from two different ICF patients with compound and different heterozygous mutations in DNMT3B gene (PT5 and GM08747) compared to LCLs (XX and MS) and fibroblasts (3674) derived from 3 unrelated controls. PT5 sample comes from an ICF male with heterozygous DNMT3B mutations (V699G/R54X). GM08747 comes from a female diagnosed for ICF syndrome with heterozygous DNMT3B mutations (A603T/STP807ins). The healthy donors samples belonged to two females (XX, GM03674) and one male (MS). All of these cell lines were obtained from Dr. RS Hansen laboratory (University of Washington, USA) and the Coriell Cell Repository and have been used in a previous published article [9]. Additionally, DNA from mononuclear cells isolated from four cord blood donors (three males and one female) was kindly provided by Dr. D. Monk from PEBC (IDIBELL, Spain). For validation steps, we extended the cohort with three lymphoblastoid cell lines: GM08714 from the ICF female patient with mutations A603T/STP807ins (see above) and its related healthy donors (GM08728 and GM08729), identified as the GM08714’s mother (GM08728) and father (GM08729) and another primary fibroblast PT3 from the ICF male with heterozygous DNMT3B mutations (V699G/R54X). Another unrelated LCL control was used (LDA). All these cell lines were also obtained and purchased from the Coriell Cell Repository. The project has been approved by the local Ethical Committee of the IDIBELL Institution.
Genome-wide DNA methylation analysis
Genome-wide DNA methylation analysis was performed using the Infinium HumanMethylation450 BeadChip from Illumina. The 450K DNA methylation array by Illumina is an established, highly reproducible method for DNA methylation detection and has been validated in two independent laboratories [10],[11].
DNA from ICF patients and healthy donors were isolated using Phenol:Chloroform:Isoamylalcohol (Sigma) and quantified by Quant-iT PicoGreen dsDNA Reagent (Invitrogen). The integrity was analyzed in a 1.3% agarose gel. Bisulfite conversion of 600 ng of each DNA sample was performed according to the manufacturer's recommendation for Illumina Infinium Assay. Effective bisulfite conversion was checked for three controls that were converted simultaneously with the samples. Four μl of bisulfite converted DNA were used to hybridize on Infinium Human Methylation 450 BeadChip, following Illumina Infinium HD Methylation protocol. Chip analysis was performed using Illumina HiScan SQ fluorescent scanner. The intensities of the images were extracted using GenomeStudio (2010.3) Methylation module (1.8.5) software (San Diego, California). Methylation score of each CpG is represented as beta (β) value. The 450K DNA Methylation array includes 485,764 cytosine positions of the human genome that were filtered by sex chromosomes CpGs (avoiding sex link alterations) and non-valid CpGs (p-value<0.001). The intensities of the images were extracted and normalized using GenomeStudio (2011.1) Methylation module (1.9.0) software. Unsupervised (using the 5000 random CpGs) and supervised heatmaps were obtained using hierarchical clustering analysis with Manhattan metrics.
For determining differentially methylated CpGs a parametric analysis using an absolute difference in beta values of 0.65 and standard deviation <0.15 in ICF1 patients compared to controls were used for selecting the most relevant positions.
Bisulfite genomic DNA sequencing
The Methyl Primer Express v1.0 software (Applied Biosystems, Life technologies, Grand Island, New York) was used to identify the CpG islands in gene promoter regions and to design specific primers for the methylation analysis (S1 Table). DNA methylation status was established by bisulfite genomic sequencing. DNA was extracted from samples using the DNeasy tissue kit (Qiagen, Milan, Italy) and 1 μg of DNA was modified with sodium bisulfite using the EZ DNA methylation-gold kit (Zymo Research, CA USA) according to manufacturer´s instructions. Multiple clones were analyzed for each sample and the methylation frequency was calculated in each case.
Gene expression analysis
Total RNA from lymphoblastoid cell lines and untrasformed fibroblasts was extracted using TRIzol reagent (Life Technologies. Grand Island, New York) and was reverse-transcribed using iScript cDNA Synthesis kit (Bio-rad). Quantitative real-time PCR (qRT-PCR) was performed using iQ Supermix SYBR Green 2X (Bio-rad. San Diego, California) on the Bio-Rad iCycler according to the manufacturer’s protocols. The ΔΔCt method [12] was used to determine relative quantitative levels. GAPDH gene was used to normalize the data. Primer sequences for gene expression analysis are shown in (S2 Table). Statistical analysis was performed using Student t test.
Results and Discussion
ICF-specific DNMT3B variants cause a global decrease in DNA methylation profile
Seminal studies regarding DNA methylation in ICF1 patients have reported significant hypomethylation at pericentromeric satellite DNA sequences in ICF cells of chromosomes 1, 9, 16, Alu sequences, D4Z4 and NBL2 repeats. Heterogeneous hypomethylation has been described at few single copy loci, X-linked and imprinted genes, while the genomic hypomethylation in ICF1 patients has been thought to be involved only in a rather small percentage of the 5- methylcytosine residues [4],[13–17]. Recently, a genome-wide DNA methylation analysis has been performed in our laboratory with the limitation of using a unique patient sample [8].
Since the majority of ICF patients deal with mutations within the DNMT3B catalytic domain expected to variably interfere with the methyltransferase activity of the protein, the primary goal of this study was to describe the global DNA methylation profile affected by ICF specific DNMT3B mutant alleles. Here, we characterized the methylome of one lymphoblastoid cell line and of untrasformed fibroblasts derived from two different compound heterozygous ICF patients with the DNMT3B mutations V699G/R54X and A603T/STP807ins (PT5 and GM08747, respectively) compared to three control LCLs (XX, MS) and fibroblasts 3674 derived from healthy donors. Using this strategy that includes ICF patients derived from different tissue types, the variability and interference due to tissue-specific genes, will be reduced. The analysis, by calculating first the averaged Beta values for each CpG from the three controls and average Beta values from the two ICF patients and later the delta values (average ICF-average controls), shows that ICFs globally contain more poorly methylated (βvalue<0.33) and less highly methylated CpGs (βvalue>0.66) compared to controls (Fig 1A). In this sense, the accumulated number of poorly methylated CpGs ranging Beta values from 0 to 0.33 of ICF patients is 215,227; while for controls decrease to 202,003 (Delta ICF-Control = +13,224). However, an opposite pattern is obtained for highly methylated CpGs ranging from Beta values 0.66 to 1. In this case, control donor showed 143,108 highly methylated CpGs compared to lower number for ICF patients 138,873 CpGs (Delta ICF-Control = -4,235) (Fig 1B). A more comprehensive representation is the scatter plot of the DNA methylation levels (βvalue) of ICF patients compared to controls showing a higher accumulation of hypomethylated CpGs in ICFs than in controls, see triangle area in Fig 1C. Confirming these results we observed, using a non-parametic Mann-withney U test after testing normality with the Shapiro-Wilk test, a significant decrease in methylation level in ICF samples compared to controls (Fig 1D). We provide individual histograms, scatter and box plots for all the hybridized samples. Individual methylation levels were consistent, although control 2 (XX) presented lower global levels than the other two controls (MS and GM03674) (S1 Fig). Therefore, our results are in agreement with previous studies reporting that ICF syndrome is a disease characterized by DNA hypomethylation and we further demonstrate that the combination of the specific DNMT3B variants A603T/STP807ins and V699G/R54X derives in a global loss of DNA methylation levels.
(A) Histograms shows bimodal distribution pattern of DNA methylation profiles in ICF patients and normal donors. The frequency of CpGs according to DNA methylation levels are depicted in the graph. (B) Table showing number of average poorly methylated (methylation levels beta<0.33) and average highly methylated (methylation levels Beta>0.66). (C) Scatter plot represents comparison of DNA methylation levels of total CpG sites using the Infinium 450K DNA methylation assay. Green triangle selects hypomethylated area for ICF patients compared to controls. (D) Box plot displaying the distribution of Beta-values of total CpG sites of ICF versus healthy control donors. Normality was tested using the Shapiro-Wilk test and significance was evaluated with the Mann-Whitney U test and is indicated by three asterisks *** (p<0.001).
Identification of differentially DNA methylated genes in ICF
Previous studies, using candidate-gene approaches, have been searching for differentially methylated regions in ICF patients that would account for the severe clinical features that characterize the ICF patients. In this sense, Jin et al reported subtle but significant changes in DNA methylation levels associated with transcript level variations in a few genes involved in development, neurogenesis and immunological function of ICF patients using expression microarrays [6]. Moreover, high degree correlation between DNA methylation changes and gene expression patterns for a number of heterochromatic genes located at the pericentromeric region of chromosome 21 in ICF patients has been also reported [7]. Therefore, a genome-scale approach would contribute to identify new differentially methylated genes and increase the knowledge in the etiology and development of the disease.
Since we observed that ICF patients with DNMT3B mutations show a global reduction of DNA methylation, we focused on studying the CpG positions with loss of methylation. To gain robustness and reliability, we added to the previous set of ICF1 samples four new controls of peripheral mononuclear cells obtained from cord blood samples (CB10, CB13, CB20 and CB76), being aware of the limitation of the cell type heterogeneity in these samples. The unsupervised hierarchical clustering, using 5000 random CpGs mimicking the global methylome, shows that samples are grouped based on their tissue type. The methylome of the four peripheral mononuclear derived cells samples is homogeneous and clustered together. The immortalized cell lines derived from controls and ICF patients clusterized in the same group, but in a separate subgroup. Finally, fibroblast cells clustered in an independent group (Fig 2A). These results indicate that the distance between samples is mainly due to tissue type.
(A) Unsupervised hierarchical clustering and heatmap of four cord blood donors (purple), three unrelated healthy donors (blue) and two ICF patients (orange) using 5000 random selected CpGs. DNA Methylation levels scale is shown. Each column represents patients and each row represents the different CpGs. (B) Supervised cluster and heatmap representing the distinctive 181 CpGs corresponding to the comparison between ICF patients (orange) and all control samples (dark blue).
We focused on ICF-associated regions which had loss of methylation, in order to obtain reliable and ICF-specific differentially methylated candidates. To achieve this, we used a restrictive threshold to overcome the limited number of samples and the heterogeneous effect of the mutations. Then, we performed a parametric analysis comparing average Beta values from ICF1 patients versus controls selecting those with differences in methylation levels higher than 65% (delta<-0.65) and a standard deviation value lower than 15% (Desvest<0.15). Using this strategy, with restrictive cut-offs, we identified 181 hypomethylated CpGs in ICF1 patients compared to control group. A table with the 181 differentially methylated positions (DMPs) with their gene characteristics including target ID, name and Beta values for ICF1 and controls groups, chromosome location, genomic distribution and CpG context is provided (Table 1). The hierarchical clustering heatmap for all these selected 181 CpG sites clearly segregates both groups (Fig 2B). It is worth emphasizing that both ICF1 patients present different tissue types and DNMT3B mutations, thereby producing subtle distinct levels of DNA methylation defects. However, the selected approach allowed us to identify the epigenetic differences that are common to both ICF1 patients leading to the establishment of potential markers for the ICF1 disease. Furthermore, we aimed to compare our selected 181 DMPs with the previous published methylome in our laboratory by Heyn et al. of the LCL GM08714 from the same ICF patient from which the GM08747 fibroblasts derive [8]. Using whole genome bisulfite sequencing (WGBS) we previously obtained 296,964 differentially methylated regions (DMRs). Consistently, we observed that a high number of 138 out of 181 (76.2%) of DMPs were inside a previously defined DMRs and therefore are common to both studies (S2A Fig). It is important to emphasize that although we are using different ICF1 samples from the previously published work, there is a high percentage of concordance that favors the reliability of our results.
Characterization of the genomic localization and gene features of obtained DMPs
ICF-specific DNA methylation changes have been previously reported mainly in pericentromeric regions of chromosomes 1, 9 and 16 [18]. However, when we analyzed the subchromosomal localization of our DMPs we observed that the majority of the DMPs are located in the intra and subtelomeric regions of the chromosomes (Figs 3 and 4A). These results support previous studies where hypomethylation of subtelomeric regions was associated in the ICF1 cells with advanced telomere replication timing and elevated levels of transcripts emanating from telomeric regions. These findings may explain the abnormal telomeric phenotype observed in ICF syndrome [19].
Gene names and intergenic CpGs are represented and localized by blue lines.
(A) Chromosomal sub-localization classified in different groups: subtelomeric, pericentromeric and intrachromosomal. (B) Associated RNA transcription classified in: coding, non-coding and intergenic. (C) CpG context and neighborhood classified in: island, shore, shelf and open sea/others. (D) Functional genomic distribution classified in: promoter (TSS1500, TSS200, 5´UTR), genic (1stexon, body and 3´UTR) and intergenic.
We further evaluated the association of the RNA transcripts with significant DMPs, finding that 39.8% and 2.2% were associated to coding and non-coding genes, respectively. Interestingly, there were 4 DMPs associated to a long non-coding RNA (Fig 4B). These results indicate that specific non coding RNAs are target of DNMT3B and possibly regulated by DNA methylation. From the CpG content and neighborhood context standpoint (Fig 4C), the CpG island, which are regions with high dense number of CpGs, are the most extensively screened regions (48.6%) and are over-represented according to the design of the 450K Infinium array where CGs in islands represent only a 31% [10]. It is shown the classification according to functional genome distribution, indicating that the majority of the DMPs are located in intergenic regions (58%) leading to an over-representation of these regions regarding the design of the 450K infinium array where the percentage of probes in intergenic region is 24.6% (Fig 4D). Meanwhile, 24.3% of the DMPs correspond to promoter regions, defined as CpGs located at TSS1500, TSS200 and UTR regions. Interestingly this promoter group is under-represented based on the expected 38.9% from the 450K array [20]. Moreover, we took advantage of results published in [8] and compared the functional genomic distribution between our 181 DMPs and the 296,964 DMRs previously reported. We observed overall a discrepancy in promoter region (24.3% vs 4.5% in our study and Heyn et al, respectively) (S2B Fig, left panel). However, this striking result may be explained by the distribution of the CpGs in the array and in the whole genome that mimics our obtained pattern (S2B Fig, right panel). Finally, analyzing a subset of DMPs located in promoter regions only, we observed that 13 (87%) are associated with CpG islands or shores. Interestingly, these regions have been reported for being important regulatory regions for disease [21],[22]. As a global conclusion, the majority of the hypomethylated DMPs in ICF patients are located at intergenic and CpG islands regions.
DNA hypomethylated candidate validation
To confirm the methylation results obtained by genome-scale techniques, we used a small scale, site-specific technique as targeted bisulfite sequencing. This technique takes advantage of the activity of sodium bisulfite that converts non-methylated cytosines into uraciles by deamination, while methyl-cytosines remain unaltered. We performed a technical validation using the same samples that were hybridized; the ICF fibroblasts and LCL (GM08747 and PT5) and 3 unrelated controls (XX, MS and 3674). Moreover, to gain robustness we also performed a biological validation using non-hybridized LCLs from the same ICF patient (GM08714) and ICF fibroblast compared with related controls (GM08728, GM08729 respectively mother and father of GM08714/GM08747) and the unrelated control (LDA). GM8714 was the ICF patient sample used in our previously paper [8].
Based on previous knowledge that epigenetic disruption of germ line function in somatic tissues has been associated with ICF [23][24], we selected four representative genes related to germ line specific pathways and/or are expressed exclusively in germ cells [25–27]. These four genes (BOLL, SYCP2, LDHAL6A and NCRNA00221) presented DNA methylation differences at crucial regulatory elements such as promoter with CpG islands. It is worth to mention that more than one DMCpGs in those important regulatory regions were identified in our analysis for BOLL and NCRNA00221, 6 and 4 respectively. Thus, coding genes such as, BOLL, SYCP2 and LDHAL6A, might be additional examples of DNMT3b target genes functionally regulated by DNA methylation in somatic tissues. Moreover, NCRNA00221 (Linc00221) is a long intergenic non-coding RNA (lincRNA). These non-coding genes are important regulators of gene expression that have been described in several diseases [28].
The pyrosequencing analysis confirmed, by comparing ICFs patients with controls, the different methylation levels in the targeted promoter CpGs of the four genes (BOLL, SYCP2, LDHAL6A and NCRNA00221) (Fig 5, left panel). Both DNA methylation assays, bisulfite genomic sequencing and 450K infinium array, showed similar levels of methylation. Right panel of Fig 5 depicts CpGs DNA methylation values (including the CpGs obtained by 450K array analysis marked with a red arrow) located in the regions analyzed for the four genes in two representative samples: an ICF patient and a healthy control. DNA methylation values for all samples (ICF patients and control donors) and a schematic representation of the four selected gene regions are shown in (S3, S4, S5 and S6 Figs). Globally, we can conclude that not only the infinium-targeted CpGs are unmethylated in ICF patients (red arrow), but the surrounding CpGs also show the same pattern favoring the idea that global DNA demethylation landscape is maintained. This effect could indicate that the entire CpG island, and not only a single base demethylation, is altered as consequence of the DNMT3B deficiency in the ICF disease and might cause gene expression deregulation.
Left panel depicts DNA methylation values (Beta) from ICF and healthy donors using both methodologies bisulfite genomic sequencing (BSG) and 450K array. Right panel shows bisulfite genomic sequencing analysis of genes in one representative control and ICF patient. CpG dinucleotides are shown in vertical lines. Multiple single clones are represented for each sample. Presence of unmethylated or methylated CpGs is indicated by white or black squares, respectively. Red arrows mark the localization of the differentially methylated CpGs by 450K array. The distance to transcription start site (bp) is also indicated.
Association of gene expression and DNA methylation changes
The existence of aberrations in the DNA methylation patterns of ICF cells, particularly the hypermethylation of the CpG island sequences located in the promoter regions of key regulatory genes, that lead to gene silencing, have been extensively described in literature [29]. Conversely, DNA hypomethylation has been associated mainly with DNA methylation loss at genome-wide level, although it also occurs locally. In this light, the effect of DNA demethylation makes accessible the transcription machinery and hence facilitating gene activation, which have been mainly described in cancer, involving the role of oncogenes [22]. Although some hypomethylated genes have been found associated to ICF1 [23], the number of disrupted epigenetic genes identified is very limited and little is known about their association with the etiology of this disease. In line with this finding, we aimed to evaluate the impact of DNA methylation changes on the transcriptional activity and detect gene promoters that are highly methylated in healthy controls a suffer de novo established loss-of-methylation in ICF patients concordant with gene up-regulation.
Therefore, based on the established idea that promoter CpG islands (CGI) are the prominent and crucial regulatory regions for gene expression and with high variable methylation level between normal and diseased tissue [30],[31], we sought to evaluate gene expression of those previously validated genes. To elucidate the impact of DNA methylation on the transcriptional activity of SYCP2, BOLL, LDHAL6A and NCRNA00221, we analyzed by qPCR their gene expression in ICF patient cells (GM08714 and PT5) and control samples including related and unrelated healthy donors (GM08728, GM08729, MS and LDA). We also analyzed in parallel the untransformed fibroblasts form ICF patient (GM08747) and the healthy donor (GM03674). The results clearly showed a significant up-regulation, using the unpaired Student t test, for three of four genes in ICF cells (except LDHAL6A that only show significance in GM08714, but not in PT5) compared to controls, indicating that the impaired DNA methylation at the identified DMPs is critical for controlling their gene activity (Fig 6). Differences in expression levels between both ICF samples in SYCP2 and BOLL could be due to DNMT3 type of mutation and the idea that, the other factors together with DNA methylation changes may be involved in the complex regulation of gene expression. In line, consistent results were obtained in fibroblasts expression analysis, except for BOLL that we were not able to detect expression in any fibroblast. This could be explained due to a tissue-specific regulation beyond DNA methylation. Interestingly, even if at different extents overall comparing levels in LCLs with untransformed fibroblasts, three genes were epigenetically regulated in ICF patients with four different DNMT3B mutations (GM8714, PT5). This suggests an important role in the DNMT3B-mediated regulatory pathway that could contribute to explain some aspects of the characteristic phenotype of these ICF patients. These results suggest a potential interaction between the DNMT3B type of mutation and the epigenetic regulation and intensity of gene activation of the studied genes.
Fold change values of the differentially DNA methylated genes in lymphoblastoid ICF patients and healthy donors were evaluated by qRT-PCR. In parallel, Fold change values were also tested in untransformed fibroblast form an ICF patient and a healthy donor. Values were determined at least in triplicate. Statistic analysis was evaluated using student t test and significance symbols correspond to (* p<0.05; ** p<0.01 and *** p<0.001).
The role of DNMT3b in protecting somatic cells against the aberrant expression of the germ line program has recently been suggested [23]. Moreover, the DNMT3b-mediated silencing of a subset of germ line genes in somatic cells occurs through the recruitment of the E2F6 transcriptional repressor at their promoter region [32]. In this light, by identifying novel germ line genes, which are hypomethylated and inappropriately expressed, our results suggest that this phenomenon in the context of DNMT3B deficiency might be rather widespread. How this specific deregulation may contribute to ICF molecular pathogenesis remains to be established. Notably, the ectopic expression of meiotic genes in cancer cells has been functionally related to abnormal chromosome segregation and aneuploidy [33]. Because chromosomal instability is a hallmark of ICF syndrome, this raises the possibility that loss of silencing at particular germ line genes drives the typical cytological abnormalities seen in ICF patient lymphocytes.
It is known that early diagnosis of ICF syndrome is crucial since early treatment can improve the course of disease. However, ICF is probably underdiagnosed, especially in patients that present incomplete phenotype or born to families with no affected relatives [23]. Therefore, the DNA hypomethylation profile of NCRNA00221 especially, and in a minor extent SYCP2 and BOLL, could be further investigated and validated in peripheral blood in order to develop specific clinical biomarkers to facilitate the identification of ICF patients.
Conclusion
Our results contribute to elucidate how different mutations in DNMT3b result in deficiency of DNA methyltransferase activity, eventually causing ICF1 syndrome. It seems accepted that is the DNA methylation deficiency, and not other aspects of impaired DNMT3b activity, responsible for the ICF syndrome.
The regions with aberrant methylation in ICF patients were almost exclusive of pericentromeric regions of chromosome 1, 16, sometimes 9 and associated to repeated DNA sequences or heterochromatic genes. Although in vitro studies identified a spectrum of biochemical defects in the catalytic function associated to ICF-specific DNMT3B mutations, their impact on the genome-wide DNA methylation level in patient-derived cells is unsolved. Our results contributed to characterize the global defects of DNA methylation pattern in two heterozygous DNMT3B-mutant backgrounds, uncovering novel ICF-specific hypomethylated sites, outside the pericentromeric regions and in other chromosomes compared to those previously mentioned. Interestingly, we identified additional DNMT3B target loci whose expression must be restricted to germ cells. Establishment and maintenance of promoter DNA methylation in somatic tissues by DNMT3B is critical for their transcriptional repression. In addition, to provide further evidence on DNMT3B role in silencing germ line genes, these findings are of particular interest in the context of other human disease, like cancer. It is remarkable that the expression of catalytically inactive DNMT3b splice variants, the aberrant transcription of germ line genes and chromosomal instability are shared features. We thus believe that these genome-wide studies will help to elucidate the relationship between DNA hypomethylation and pathological phenotypes.
Finally, from the ICF syndrome point of view, our results contribute to further evaluate the utility of these potential biomarkers as diagnostic markers.
Supporting Information
S1 Fig. Individual genome-wide DNA methylation profiles in ICF patients and control samples.
(A) Histograms showing bimodal distribution pattern of DNA methylation profiles in the two ICF patients (PT5 and GM8748) and normal donors (XX, MS and GM3674). (B) Table showing number of poorly methylated and highly methylated for each sample hybridized (C) Individual scatter plots combining ICF patients and normal donors. (D) Individual box plots for ICF patients compared to normal donors. Normality was tested using the Shapiro-Wilk test and significance was evaluated with the Mann-Whitney U test and is indicated by three asterisks *** (p<0.001).
https://doi.org/10.1371/journal.pone.0132517.s001
(TIF)
S2 Fig. Comparison of Hypomethylated regions reported by Heyn et al. and current differentially methylated positions.
(A) Venn diagram illustrating the number of hypomethylated features. The intersecting region represents those CpGs that are common to both analyses. (B) Left panel: graph depicting percentages of differentially methylated features based on functional genomic distribution comparing both analyses (dark blue for our analysis and dark red for Heyn et al.). Right panel: graph showing percentages based on the 450K array design and the whole genome (light blue for 450K array and red for whole genome). Data obtained from [20].
https://doi.org/10.1371/journal.pone.0132517.s002
(TIF)
S3 Fig. Bisulfite genomic sequencing for technical and biological validation of differentially DNA methylated gene BOLL.
For technical validation two ICF (PT5 and GM08747) and two unrelated controls (XX, MS and GM03674) were used. For biological validation two ICF patient samples (GM8714 and PT3) and three controls (two related controls GM8728, GM8729 and one unrelated control LDA) were analyzed. CpG dinucleotides is shown in vertical lines. Multiple single clones are represented for each sample. Presence of unmethylated or methylated CpGs is indicated by white or black squares, respectively. Red arrows mark the localization of the differentially methylated CpGs by 450K array. The distance to Trasncription Start Site (bp) is also indicated.
https://doi.org/10.1371/journal.pone.0132517.s003
(TIF)
S4 Fig. Bisulfite genomic sequencing for technical and biological validation of differentially DNA methylated gene SYCP2.
The design is similar to S3 Fig.
https://doi.org/10.1371/journal.pone.0132517.s004
(TIF)
S5 Fig. Bisulfite genomic sequencing for technical and biological validation of differentially DNA methylated gene LDHAL6A.
The design is similar to S3 Fig.
https://doi.org/10.1371/journal.pone.0132517.s005
(TIF)
S6 Fig. Bisulfite genomic sequencing for technical and biological validation of differentially DNA methylated gene NCRNA00221.
The design is similar to S3 Fig.
https://doi.org/10.1371/journal.pone.0132517.s006
(TIF)
S1 Table. Primers for bisulfite genomic sequencing.
https://doi.org/10.1371/journal.pone.0132517.s007
(DOCX)
Acknowledgments
We thank Diana Garcia, Carles Arribas and Sebastian Moran for their technical support and Dr. Sylwia Leppert for her help during the early phase of the work.
Author Contributions
Conceived and designed the experiments: JS ME. Performed the experiments: LSR SG ADL MG ABC MM JS. Analyzed the data: JS ME. Contributed reagents/materials/analysis tools: JS ME MM. Wrote the paper: JS ME MM.
References
- 1. Ehrlich M, Jackson K, Weemaes C. (2006) Immunodeficiency, centromeric region instability, facial anomalies syndrome (ICF). Orphanet J Rare Dis 1: 2. pmid:16722602
- 2. Weemaes CM, van Tol MJ, Wang J, van Ostaijen-ten Dam MM, van Eggermond MC, Thijssen PE et al. (2013) Heterogeneous clinical presentation in ICF syndrome: correlation with underlying gene defects. Eur J Hum Genet. 21(11):1219–25. pmid:23486536
- 3. Moarefi AH1, Chédin F. (2011) ICF syndrome mutations cause a broad spectrum of biochemical defects in DNMT3B-mediated de novo DNA methylation. J Mol Biol. Jun 24;409(5):758–72. pmid:21549127
- 4. Kondo T, Bobek MP, Kuick R, Lamb B, Zhu X, Narayan A et al. (2000) Whole-genome methylation scan in ICF syndrome: hypomethylation of non-satellite DNA repeats D4Z4 and NBL2. Hum Mol Genet. 9(4):597–604. pmid:10699183
- 5. Ehrlich M, Buchanan KL, Tsien F, Jiang G, Sun B, Uicker W et al. (2001) DNA methyltransferase 3B mutations linked to the ICF syndrome cause dysregulation of lymphogenesis genes. Hum Mol Genet 10(25):2917–31. pmid:11741835
- 6. Jin B, Tao Q, Peng J, Soo HM, Wu W, Ying J et al. (2008) DNA methyltransferase 3B (DNMT3B) mutations in ICF syndrome lead to altered epigenetic modifications and aberrant expression of genes regulating development, neurogenesis and immune function. Hum Mol Genet 17(5):690–709. pmid:18029387
- 7. Brun ME, Lana E, Rivals I, Lefranc G, Sarda P, Claustres M et al. (2011). Heterochromatic genes undergo epigenetic changes and escape silencing in immunodeficiency, centromeric instability, facial anomalies (ICF) syndrome. PLoS One 6(4):e19464. pmid:21559330
- 8. Heyn H, Vidal E, Sayols S, Sanchez-Mut JV, Moran S, Medina I et al. (2012) Whole-genome bisulfite DNA sequencing of a DNMT3B mutant patient. Epigenetics 7(6):542–50. pmid:22595875
- 9. Gatto S, Della Ragione F, Cimmino A, Strazzullo M, Fabbri M, Mutarelli M et al. (2010) Epigenetic alteration of microRNAs in DNMT3B-mutated patients of ICF syndrome. Epigenetics 5(5):427–43. pmid:20448464
- 10. Sandoval J, Heyn HA, Moran S, Serra-Musach J, Pujana MA, Bibikova M et al. (2011) Validation of a DNA methylation microarray for 450,000 CpG sites in the human genome. Epigenetics 6(6):692–702. pmid:21593595
- 11. Bibikova M, Barnes B, Tsan C, Ho V, Klotzle B, Le JM et al. (2011) High density DNA methylation array with single CpG site resolution. Genomics 98(4):288–95. pmid:21839163
- 12. Livak KJ and Schmittgen TD. (2001) Analysis of relative gene expression data using real-time quantitative PCR and the 2(-Delta Delta C(T)) Method. Methods 25(4):402–8. pmid:11846609
- 13. Miniou P, Jeanpierre M, Blanquet V, Sibella V, Bonneau D, Herbelin C et al. (1994) Abnormal methylation pattern in constitutive and facultative (X inactive chromosome) heterochromatin of ICF patients. Hum Mol Genet 3(12):2093–102. pmid:7881405
- 14. Miniou P, Jeanpierre M, Bourc'his D, Coutinho Barbosa AC, Blanquet V and Viegas-Péquignot E. (1997) Alpha-satellite DNA methylation in normal individuals and in ICF patients: heterogeneous methylation of constitutive heterochromatin in adult and fetal tissues. Hum Genet 99(6):738–45. pmid:9187666
- 15. Hansen RS, Stöger R, Wijmenga C, Stanek AM, Canfield TK, Luo P et al. (2000) Escape from gene silencing in ICF syndrome: evidence for advanced replication time as a major determinant. Hum Mol Genet. 9(18):2575–87. pmid:11063717
- 16. Matarazzo MR, De Bonis ML, Gregory RI, Vacca M, Hansen RS, Mercadante G et al. (2002) Allelic inactivation of the pseudoautosomal gene SYBL1 is controlled by epigenetic mechanisms common to the X and Y chromosomes. Hum Mol Genet. 11(25):3191–8. pmid:12444103
- 17. Matarazzo MR, De Bonis ML, Strazzullo M, Cerase A, Ferraro M, Vastarelli P et al. (2007). Multiple binding of methyl-CpG and polycomb proteins in long-term gene silencing events. J Cell Physiol. 210(3):711–9. pmid:17133344
- 18. Jeanpierre M, Turleau C, Aurias A, Prieur M, Ledeist F, Fischer A et al. (1993) An embryonic-like methylation pattern of classical satellite DNA is observed in ICF syndrome. Hum Mol Genet. 2(6):731–5. pmid:8102570
- 19. Yehezkel S, Segev Y, Viegas-Péquignot E, Skorecki K and Selig S. (2008) Hypomethylation of subtelomeric regions in ICF syndrome is associated with abnormally short telomeres and enhanced transcription from telomeric regions. Hum Mol Genet. 17(18):2776–89. pmid:18558631
- 20. Stirzaker C, Taberlay PC, Statham AL, Clark SJ. (2014) Mining cancer methylomes: prospects and challenges. Trends Genet. 30(2):75–84. pmid:24368016
- 21. Irizarry RA, Ladd-Acosta C, Wen B, Wu Z, Montano C, Onyango P et al. (2009) The human colon cancer methylome shows similar hypo- and hypermethylation at conserved tissue-specific CpG island shores. Nat Genet 41(2):178–86. pmid:19151715
- 22. Sandoval J and Esteller M. (2012) Cancer epigenomics: beyond genomics. Curr Opin Genet Dev. 22(1):50–5. pmid:22402447
- 23. Velasco G, Walton EL, Sterlin D, Hédouin S, Nitta H, Ito Y et al. (2014) Germline genes hypomethylation and expression define a molecular signature in peripheral blood of ICF patients: implications for diagnosis and etiology. Orphanet J Rare Dis. 9:56. pmid:24742017
- 24. Walton EL, Francastel C and Velasco G. (2014) Dnmt3b Prefers Germ Line Genes and Centromeric Regions: Lessons from the ICF Syndrome and Cancer and Implications for Diseases. Biology (Basel) 3(3): 578–605.
- 25. Li M, Liu C, Zhu H, Sun J, Yu M, Song W et al. (2013) Expression pattern of Boule in dairy goat testis and its function in promoting the meiosis in male germline stem cells (mGSCs). J Cell Biochem. 114(2):294–302. pmid:22930651
- 26. Zheng YH, Rengaraj D, Choi JW, Park KJ, Lee SI and Han JY. (2009) Expression pattern of meiosis associated SYCP family members during germline development in chickens. Reproduction 138(3):483–92. pmid:19525366
- 27. Wang H, Zhou Z, Lu L, Xu Z and Sha J. (2005) Cloning and characterization of a novel intronless lactate dehydrogenase gene in human testis. Int J Mol Med 15(6):949–53. pmid:15870898
- 28. Della Ragione F, Gagliardi M, D'Esposito M and Matarazzo MR. (2014) Non-coding RNAs in chromatin disease involving neurological defects Front Cell Neurosci. 8:54. pmid:24616662
- 29. Esteller M. (2008) Epigenetics in cancer. N Engl J Med 358(11):1148–59. pmid:18337604
- 30. Lopez-Serra P, Marcilla M, Villanueva A, Ramos-Fernandez A, Palau A, leal L et al. (2014) DERL3-associated defect in the degradation of SLC2A1 mediates the Warburg effect. Nat Commun. 5:3608. pmid:24699711
- 31. Quintavalle C, Mangani D, Roscigno G, Romano G, Diaz-Lagares A, Iaboni M et al. (2013) MiR-221/222 target the DNA methyltransferase MGMT in glioma cells. Plos One 8(9):e74466. pmid:24147153
- 32. Velasco G, Hube F, Rollin J, Neuillet D, Philippe C, Bouzinba-Segard H et al. (2010). Dnmt3b recruitment through E2F6 transcriptional repressor mediates germ-line gene silencing in murine somatic tissues. Proc Natl Acad Sci U S A 107:9281–9286. pmid:20439742
- 33. Simpson AJ, Caballero OL, Jungbluth A, Chen YT, Old LJ. (2005) Cancer/testis antigens, gametogenesis and cancer. Nat. Rev. Cancer 5: 615–625. pmid:16034368