Genome-wide DNA methylation profiling in rheumatoid arthritis identifies disease-associated methylation changes that are distinct to individual T- and B-lymphocyte populations

Changes to the DNA methylome have been described in patients with rheumatoid arthritis (RA). In previous work, we reported genome-wide methylation differences in T-lymphocyte and B-lymphocyte populations from healthy individuals. Now, using HumanMethylation450 BeadChips to interrogate genome-wide DNA methylation, we have determined disease-associated methylation changes in blood-derived T- and B-lymphocyte populations from 12 female patients with seropositive established RA, relative to 12 matched healthy individuals. Array data were analyzed using NIMBL software and bisulfite pyrosequencing was used to validate array candidates. Genome-wide DNA methylation, determined by analysis of LINE-1 sequences, revealed higher methylation in B-lymphocytes compared with T-lymphocytes (P ≤ 0.01), which is consistent with our findings in healthy individuals. Moreover, loci-specific methylation differences that distinguished T-lymphocytes from B-lymphocytes in healthy individuals were also apparent in RA patients. However, disease-associated methylation differences were also identified in RA. In these cases, we identified 509 and 252 CpGs in RA-derived T- and B-lymphocytes, respectively, that showed significant changes in methylation compared with their cognate healthy counterparts. Moreover, this included a restricted set of 32 CpGs in T-lymphocytes and 20 CpGs in B-lymphocytes (representing 15 and 10 genes, respectively, and including two, MGMT and CCS, that were common to both cell types) that displayed more substantial changes in methylation. These changes, apparent as hyper- or hypo-methylation, were independently confirmed by pyrosequencing analysis. Validation by pyrosequencing also revealed additional sites in some candidate genes that also displayed altered methylation in RA. In this first study of genome-wide DNA methylation in individual T- and B-lymphocyte populations in RA patients, we report disease-associated methylation changes that are distinct to each cell type and which support a role for discrete epigenetic regulation in this disease.     

Genome-wide DNA methylation profiling in rheumatoid arthritis identifies disease-associated methylation changes that are distinct to individual T- and B-lymphocyte populations

Changes to the DNA methylome have been described in patients with rheumatoid arthritis (RA). In previous work, we reported genome-wide methylation differences in T-lymphocyte and B-lymphocyte populations from healthy individuals. Now, using HumanMethylation450 BeadChips to interrogate genome-wide DNA methylation, we have determined disease-associated methylation changes in blood-derived T- and B-lymphocyte populations from 12 female patients with seropositive established RA, relative to 12 matched healthy individuals. Array data were analyzed using NIMBL software and bisulfite pyrosequencing was used to validate array candidates. Genome-wide DNA methylation, determined by analysis of LINE-1 sequences, revealed higher methylation in B-lymphocytes compared with T-lymphocytes (P ≤ 0.01), which is consistent with our findings in healthy individuals. Moreover, loci-specific methylation differences that distinguished T-lymphocytes from B-lymphocytes in healthy individuals were also apparent in RA patients. However, disease-associated methylation differences were also identified in RA. In these cases, we identified 509 and 252 CpGs in RA-derived T- and B-lymphocytes, respectively, that showed significant changes in methylation compared with their cognate healthy counterparts. Moreover, this included a restricted set of 32 CpGs in T-lymphocytes and 20 CpGs in B-lymphocytes (representing 15 and 10 genes, respectively, and including two, MGMT and CCS, that were common to both cell types) that displayed more substantial changes in methylation. These changes, apparent as hyper- or hypo-methylation, were independently confirmed by pyrosequencing analysis. Validation by pyrosequencing also revealed additional sites in some candidate genes that also displayed altered methylation in RA. In this first study of genome-wide DNA methylation in individual T- and B-lymphocyte populations in RA patients, we report disease-associated methylation changes that are distinct to each cell type and which support a role for discrete epigenetic regulation in this disease.     

Genome-Wide DNA Methylation Profiles Indicate CD8+ T Cell Hypermethylation in Multiple Sclerosis

Objective

Determine whether MS-specific DNA methylation profiles can be identified in whole blood or purified immune cells from untreated MS patients.

Methods

Whole blood, CD4+ and CD8+ T cell DNA from 16 female, treatment naïve MS patients and 14 matched controls was profiled using the HumanMethylation450K BeadChip. Genotype data were used to assess genetic homogeneity of our sample and to exclude potential SNP-induced DNA methylation measurement errors.

Results

As expected, significant differences between CD4+ T cells, CD8+ T cells and whole blood DNA methylation profiles were observed, regardless of disease status. Strong evidence for hypermethylation of CD8+ T cell, but not CD4+ T cell or whole blood DNA in MS patients compared to controls was observed. Genome-wide significant individual CpG-site DNA methylation differences were not identified. Furthermore, significant differences in gene DNA methylation of 148 established MS-associated risk genes were not observed.

Conclusion

While genome-wide significant DNA methylation differences were not detected for individual CpG-sites, strong evidence for DNA hypermethylation of CD8+ T cells for MS patients was observed, indicating a role for DNA methylation in MS. Further, our results suggest that large DNA methylation differences for CpG-sites tested here do not contribute to MS susceptibility. In particular, large DNA methylation differences for CpG-sites within 148 established MS candidate genes tested in our study cannot explain missing heritability. Larger studies of homogenous MS patients and matched controls are warranted to further elucidate the impact of CD8+ T cell and more subtle DNA methylation changes in MS development and pathogenesis.     

EBV transformation and cell culturing destabilizes DNA methylation in human lymphoblastoid cell lines

Recent research suggests that epigenetic alterations involving DNA methylation can be causative for neurodevelopmental, growth and metabolic disorders. Although lymphoblastoid cell lines have been an invaluable resource for the study of both genetic and epigenetic disorders, the impact of EBV transformation, cell culturing and freezing on epigenetic patterns is unknown. We compared genome-wide DNA methylation patterns of four white blood cell samples, four low-passage lymphoblastoid cell lines pre and post freezing and four high-passage lymphobastoid cell lines, using two microarray platforms: Illumina HumanMethylation27 platform containing 27,578 CpG sites and Agilent Human CpG island Array containing 27,800 CpG islands. Comparison of genome-wide methylation profiles between white blood cells and lymphoblastoid cell lines demonstrated methylation alterations in lymphoblastoid cell lines occurring at random genomic locations. These changes were more profound in high-passage cells. Freezing at low-passages did not have a significant effect on DNA methylation. Methylation changes were observed in several imprinted differentially methylated regions, including DIRAS3, NNAT, H19, MEG3, NDN and MKRN3, but not in known imprinting centers. Our results suggest that lymphoblastoid cell lines should be used with caution for the identification of disease-associated DNA methylation changes or for discovery of new imprinted genes, as the methylation patterns seen in these cell lines may not always be representative of DNA methylation present in the original B-lymphocytes of the patient.     

Genome-wide profiling of DNA methylation reveals a class of normally methylated CpG island promoters

The role of CpG island methylation in normal development and cell differentiation is of keen interest, but remains poorly understood. We performed comprehensive DNA methylation profiling of promoter regions in normal peripheral blood by methylated CpG island amplification in combination with microarrays. This technique allowed us to simultaneously determine the methylation status of 6,177 genes, 92% of which include dense CpG islands. Among these 5,549 autosomal genes with dense CpG island promoters, we have identified 4.0% genes that are nearly completely methylated in normal blood, providing another exception to the general rule that CpG island methylation in normal tissue is limited to X inactivation and imprinted genes. We examined seven genes in detail, including ANKRD30A, FLJ40201, INSL6, SOHLH2, FTMT, C12orf12, and DPPA5. Dense promoter CpG island methylation and gene silencing were found in normal tissues studied except testis and sperm. In both tissues, bisulfite cloning and sequencing identified cells carrying unmethylated alleles. Interestingly, hypomethylation of several genes was associated with gene activation in cancer. Furthermore, reactivation of silenced genes could be induced after treatment with a DNA demethylating agent or in a cell line lacking DNMT1 and/or DNMT3b. Sequence analysis identified five motifs significantly enriched in this class of genes, suggesting that cis-regulatory elements may facilitate preferential methylation at these promoter CpG islands. We have identified a group of non-X-linked bona fide promoter CpG islands that are densely methylated in normal somatic tissues, escape methylation in germline cells, and for which DNA methylation is a primary mechanism of tissue-specific gene silencing.     

Identification of a DNA methylome signature of esophageal squamous cell carcinoma and potential epigenetic biomarkers

Esophageal squamous cell carcinoma (ESCC) is believed to arise from esophageal mucosa through accumulation of both genetic and epigenetic changes. DNA methylation is a critical epigenetic mechanism involved in key cellular processes and its deregulation has been linked to many human cancers, including ESCC. The aim of this study is to examine the global deregulation of methylation states in ESCC and identify potential early biomarkers. With this purpose, we performed a bead array analysis of more than 800 cancer-related genes in ten ESCC samples, ten matched surrounding tissues and four esophageal mucosa from healthy individuals. Pyrosequencing was used for validation of DNA methylation changes in up to 106 cases and 27 controls. A total of 37 CpG sites were found to be differentially methylated between tumors and surrounding tissues. These CpG sites were significantly enriched in genes related to several pathways including IL-10 anti-inflammatory signaling pathway and cell communication pathway. In addition, by comparing with healthy esophageal mucosa, we identified TFF1 gene as a potential early marker of ESCC. This is the first study to address methylation changes in ESCC in a large set of genes. Methylome analysis is shown as a sensitive and powerful tool to identify molecular players in ESCC. These data should prove to be the reference for future studies identifying potential biomarkers and molecular targets in ESCC.     

Identification of New Differentially Methylated Genes That Have Potential Functional Consequences in Prostate Cancer

Many differentially methylated genes have been identified in prostate cancer (PCa), primarily using candidate gene-based assays. Recently, several global DNA methylation profiles have been reported in PCa, however, each of these has weaknesses in terms of ability to observe global DNA methylation alterations in PCa. We hypothesize that there remains unidentified aberrant DNA methylation in PCa, which may be identified using higher resolution assay methods. We used the newly developed Illumina HumanMethylation450 BeadChip in PCa (n = 19) and adjacent normal tissues (n = 4) and combined these with gene expression data for identifying new DNA methylation that may have functional consequences in PCa development and progression. We also confirmed our methylation results in an independent data set. Two aberrant DNA methylation genes were validated among an additional 56 PCa samples and 55 adjacent normal tissues. A total 28,735 CpG sites showed significant differences in DNA methylation (FDR adjusted P<0.05), defined as a mean methylation difference of at least 20% between PCa and normal samples. Furthermore, a total of 122 genes had more than one differentially methylated CpG site in their promoter region and a gene expression pattern that was inverse to the direction of change in DNA methylation (e.g. decreased expression with increased methylation, and vice-versa). Aberrant DNA methylation of two genes, AOX1 and SPON2, were confirmed via bisulfate sequencing, with most of the respective CpG sites showing significant differences between tumor samples and normal tissues. The AOX1 promoter region showed hypermethylation in 92.6% of 54 tested PCa samples in contrast to only three out of 53 tested normal tissues. This study used a new BeadChip combined with gene expression data in PCa to identify novel differentially methylated CpG sites located within genes. The newly identified differentially methylated genes may be used as biomarkers for PCa diagnosis.     

Maternal genome-wide DNA methylation profiling in gestational diabetes shows distinctive disease-associated changes relative to matched healthy pregnancies

Several recent reports have described associations between gestational diabetes (GDM) and changes to the epigenomic landscape where the DNA samples were derived from either cord or placental sources. We employed genome-wide 450K array analysis to determine changes to the epigenome in a unique cohort of maternal blood DNA from 11 pregnant women prior to GDM development relative to matched controls. Hierarchical clustering segregated the samples into 2 distinct clusters comprising GDM and healthy pregnancies. Screening identified 100 CpGs with a mean β-value difference of ≥0.2 between cases and controls. Using stringent criteria, 5 CpGs (within COPS8, PIK3R5, HAAO, CCDC124, and C5orf34 genes) demonstrated potentials to be clinical biomarkers as revealed by differential methylation in 8 of 11 women who developed GDM relative to matched controls. We identified, for the first time, maternal methylation changes prior to the onset of GDM that may prove useful as biomarkers for early therapeutic intervention.     

Differential DNA methylation in purified human blood cells: implications for cell lineage and studies on disease susceptibility

Methylation of cytosines at CpG sites is a common epigenetic DNA modification that can be measured by a large number of methods, now even in a genome-wide manner for hundreds of thousands of sites. The application of DNA methylation analysis is becoming widely popular in complex disorders, for example, to understand part of the "missing heritability". The DNA samples most readily available for methylation studies are derived from whole blood. However, blood consists of many functionally and developmentally distinct cell populations in varying proportions. We studied whether such variation might affect the interpretation of methylation studies based on whole blood DNA. We found in healthy male blood donors there is important variation in the methylation profiles of whole blood, mononuclear cells, granulocytes, and cells from seven selected purified lineages. CpG methylation between mononuclear cells and granulocytes differed for 22% of the 8252 probes covering the selected 343 genes implicated in immune-related disorders by genome-wide association studies, and at least one probe was differentially methylated for 85% of the genes, indicating that whole blood methylation results might be unintelligible. For individual genes, even if the overall methylation patterns might appear similar, a few CpG sites in the regulatory regions may have opposite methylation patterns (i.e., hypo/hyper) in the main blood cell types. We conclude that interpretation of whole blood methylation profiles should be performed with great caution and for any differences implicated in a disorder, the differences resulting from varying proportions of white blood cell types should be considered.     

Symmetrical Dose-Dependent DNA-Methylation Profiles in Children with Deletion or Duplication of 7q11.23

Epigenetic dysfunction has been implicated in a growing list of disorders that include cancer, neurodevelopmental disorders, and neurodegeneration. Williams syndrome (WS) and 7q11.23 duplication syndrome (Dup7) are rare neurodevelopmental disorders with broad phenotypic spectra caused by deletion and duplication, respectively, of a 1.5-Mb region that includes several genes with a role in epigenetic regulation. We have identified striking differences in DNA methylation across the genome between blood cells from children with WS or Dup7 and blood cells from typically developing (TD) children. Notably, regions that were differentially methylated in both WS and Dup7 displayed a significant and symmetrical gene-dose-dependent effect, such that WS typically showed increased and Dup7 showed decreased DNA methylation. Differentially methylated genes were significantly enriched with genes in pathways involved in neurodevelopment, autism spectrum disorder (ASD) candidate genes, and imprinted genes. Using alignment with ENCODE data, we also found the differentially methylated regions to be enriched with CCCTC-binding factor (CTCF) binding sites. These findings suggest that gene(s) within 7q11.23 alter DNA methylation at specific sites across the genome and result in dose-dependent DNA-methylation profiles in WS and Dup7. Given the extent of DNA-methylation changes and the potential impact on CTCF binding and chromatin regulation, epigenetic mechanisms most likely contribute to the complex neurological phenotypes of WS and Dup7. Our findings highlight the importance of DNA methylation in the pathogenesis of WS and Dup7 and provide molecular mechanisms that are potentially shared by WS, Dup7, and ASD.     

Identification of a DNA methylome signature of esophageal squamous cell carcinoma and potential epigenetic biomarkers

Esophageal squamous cell carcinoma (ESCC) is believed to arise from esophageal mucosa through accumulation of both genetic and epigenetic changes. DNA methylation is a critical epigenetic mechanism involved in key cellular processes and its deregulation has been linked to many human cancers, including ESCC. The aim of this study is to examine the global deregulation of methylation states in ESCC and identify potential early biomarkers. With this purpose, we performed a bead array analysis of more than 800 cancer-related genes in ten ESCC samples, ten matched surrounding tissues and four esophageal mucosa from healthy individuals. Pyrosequencing was used for validation of DNA methylation changes in up to 106 cases and 27 controls. A total of 37 CpG sites were found to be differentially methylated between tumors and surrounding tissues. These CpG sites were significantly enriched in genes related to several pathways including IL-10 anti-inflammatory signaling pathway and cell communication pathway. In addition, by comparing with healthy esophageal mucosa, we identified TFF1 gene as a potential early marker of ESCC. This is the first study to address methylation changes in ESCC in a large set of genes. Methylome analysis is shown as a sensitive and powerful tool to identify molecular players in ESCC. These data should prove to be the reference for future studies identifying potential biomarkers and molecular targets in ESCC.     

Variations in DNA Methylation Patterns During the Cell Cycle of HeLa Cells

DNA methylation has been viewed as a stable component of the epigenome, which is established during development and fixed thereafter. We show here using nearest neighbor analysis, immunocytochemistry, and high performance capillary electrophoresis that the DNA methylation pattern varies in HeLa cells during a single cell cycle. Immunocytochemical analysis in primary human fibroblasts shows similar variations. The global levels of DNA methylation decreased in G1 and increase during the S phase of the cell cycle. Since there was little change in the DNA methylation levels in repetitive sequences throughout the cell cycle, we examined the DNA methylation pattern of unique sequences using a human CpG island microarray. Hybridization with methylated DNA from G1 and S phase of the cell cycle revealed that 174 CG-containing sequences were differentially methylated between G1 and S. 75% of all the variations in DNA methylation detected in unique sequences represented hypomethylation at G0, with changes occurring in both CpG islands and non-CpG islands. Bisulfite mapping confirmed these changes in methylation in the regions identified by the microarray. This is the first demonstration of a dynamic DNA methylation pattern within a single cell cycle of a mature somatic cell. These data are important for our understanding of the stability of DNA methylation patterns in somatic cells.     

Altered DNA methylation in leukocytes with trisomy 21

The primary abnormality in Down syndrome (DS), trisomy 21, is well known; but how this chromosomal gain produces the complex DS phenotype, including immune system defects, is not well understood. We profiled DNA methylation in total peripheral blood leukocytes (PBL) and T-lymphocytes from adults with DS and normal controls and found gene-specific abnormalities of CpG methylation in DS, with many of the differentially methylated genes having known or predicted roles in lymphocyte development and function. Validation of the microarray data by bisulfite sequencing and methylation-sensitive Pyrosequencing (MS-Pyroseq) confirmed strong differences in methylation (p<0.0001) for each of 8 genes tested: TMEM131, TCF7, CD3Z/CD247, SH3BP2, EIF4E, PLD6, SUMO3, and CPT1B, in DS versus control PBL. In addition, we validated differential methylation of NOD2/CARD15 by bisulfite sequencing in DS versus control T-cells. The differentially methylated genes were found on various autosomes, with no enrichment on chromosome 21. Differences in methylation were generally stable in a given individual, remained significant after adjusting for age, and were not due to altered cell counts. Some but not all of the differentially methylated genes showed different mean mRNA expression in DS versus control PBL; and the altered expression of 5 of these genes, TMEM131, TCF7, CD3Z, NOD2, and NPDC1, was recapitulated by exposing normal lymphocytes to the demethylating drug 5-aza-2'deoxycytidine (5aza-dC) plus mitogens. We conclude that altered gene-specific DNA methylation is a recurrent and functionally relevant downstream response to trisomy 21 in human cells.     

Differential DNA methylation with age displays both common and dynamic features across human tissues that are influenced by CpG landscape

Background

DNA methylation is an epigenetic modification that changes with age in human tissues, although the mechanisms and specificity of this process are still poorly understood. We compared CpG methylation changes with age across 283 human blood, brain, kidney, and skeletal muscle samples using methylation arrays to identify tissue-specific age effects.

Results

We found age-associated CpGs (ageCGs) that are both tissue-specific and common across tissues. Tissue-specific ageCGs are frequently located outside CpG islands with decreased methylation, and common ageCGs show the opposite trend. AgeCGs are significantly associated with poorly expressed genes, but those with decreasing methylation are linked with higher tissue-specific expression levels compared with increasing methylation. Therefore, tissue-specific gene expression may protect against common age-dependent methylation. Distinguished from other tissues, skeletal muscle ageCGs are more associated with expression, enriched near genes related to myofiber contraction, and closer to muscle-specific CTCF binding sites. Kidney-specific ageCGs are more increasingly methylated compared to other tissues as measured by affiliation with kidney-specific expressed genes. Underlying chromatin features also mark common and tissue-specific age effects reflective of poised and active chromatin states, respectively. In contrast with decreasingly methylated ageCGs, increasingly methylated ageCGs are also generally further from CTCF binding sites and enriched within lamina associated domains.

Conclusions

Our data identified common and tissue-specific DNA methylation changes with age that are reflective of CpG landscape and suggests both common and unique alterations within human tissues. Our findings also indicate that a simple epigenetic drift model is insufficient to explain all age-related changes in DNA methylation.     

Maternal and paternal chromosomes 7 show differential methylation of many genes in lymphoblast DNA

Genomic imprinting, the differential expression of paternal and maternal alleles, involves many chromosomal regions and plays a role in development and growth. Differential methylation of maternal and paternal alleles is a hallmark of imprinted genes, and thus methylation assays are widely used to support the identification of novel imprinted genes. Either blood or lymphoblast DNAs are most often used in these assays, even though methylation levels may change in cell culture. We undertook a systematic survey of parent-of-origin-specific methylation of chromosome 7 genes and ESTs by comparing DNA samples from cases of maternal and paternal uniparental disomy for chromosome 7 using DNA from fresh blood and lymphoblast cell lines. Our results revealed that up to 41% of genes and ESTs show parent-of-origin-specific methylation differences in lymphoblast DNA after only a short time in culture, whereas methylation differences were not seen in blood DNA. The methylation changes occurred most commonly on paternal chromosome 7, whereas alterations on maternal chromosome 7 were more infrequent and weaker. These findings indicate that methylation patterns may change significantly during cell culture in a parent-of-origin-dependent manner and suggest that methylation is maintained differently on maternal and paternal chromosomes 7.