DNA methylation of cord blood cell types: Applications for mixed cell birth studies

Epigenome-wide association studies of disease widely use DNA methylation measured in blood as a surrogate tissue. Cell proportions can vary between people and confound associations of exposure or outcome. An adequate reference panel for estimating cell proportions from adult whole blood for DNA methylation studies is available, but an analogous cord blood cell reference panel is not yet available. Cord blood has unique cell types and the epigenetic signatures of standard cell types may not be consistent throughout the life course. Using magnetic bead sorting, we isolated cord blood cell types (nucleated red blood cells, granulocytes, monocytes, natural killer cells, B cells, CD4⁺T cells, and CD8⁺T cells) from 17 live births at Johns Hopkins Hospital. We confirmed enrichment of the cell types using fluorescence assisted cell sorting and ran DNA from the separated cell types on the Illumina Infinium HumanMethylation450 BeadChip array. After filtering, the final analysis was on 104 samples at 429,794 probes. We compared cell type specific signatures in cord to each other and methylation at 49.2% of CpG sites on the array differed by cell type (F-test P < 10^?8). Differences between nucleated red blood cells and the remainder of the cell types were most pronounced (36.9% of CpG sites at P < 10^?8) and 99.5% of these sites were hypomethylated relative to the other cell types. We also compared the mean-centered sorted cord profiles to the available adult reference panel and observed high correlation between the overlapping cell types for granulocytes and monocytes (both r=0.74), and poor correlation for CD8⁺T cells and NK cells (both r=0.08). We further provide an algorithm for estimating cell proportions in cord blood using the newly developed cord reference panel, which estimates biologically plausible cell proportions in whole cord blood samples.     

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.     

Orthology between genomes of Brachypodium, wheat and rice

Background

Differential expression of perforin (PRF1), a gene with a pivotal role in immune surveillance, can be attributed to differential methylation of CpG sites in its promoter region. A reproducible method for quantitative and CpG site-specific determination of perforin methylation is required for molecular epidemiologic studies of chronic diseases with immune dysfunction.

Findings

We developed a pyrosequencing based method to quantify site-specific methylation levels in 32 out of 34 CpG sites in the PRF1 promoter, and also compared methylation pattern in DNAs extracted from whole blood drawn into PAXgene blood DNA tubes (whole blood DNA) or DNA extracted from peripheral blood mononuclear cells (PBMC DNA) from the same normal subjects. Sodium bisulfite treatment of DNA and touchdown PCR were highly reproducible (coefficient of variation 1.63 to 2.18%) to preserve methylation information. Application of optimized pyrosequencing protocol to whole blood DNA revealed that methylation level varied along the promoter in normal subjects with extremely high methylation (mean 86%; range 82–92%) in the distal enhancer region (CpG sites 1–10), a variable methylation (range 49%–83%) in the methylation sensitive region (CpG sites 11–17), and a progressively declining methylation level (range 12%–80%) in the proximal promoter region (CpG sites 18–32) of PRF1. This pattern of methylation remained the same between whole blood and PBMC DNAs, but the absolute values of methylation in 30 out of 32 CpG sites differed significantly, with higher values for all CpG sites in the whole blood DNA.

Conclusion

This reproducible, site-specific and quantitative method for methylation determination of PRF1 based on pyrosequencing without cloning is well suited for large-scale molecular epidemiologic studies of diseases with immune dysfunction. PBMC DNA may be better suited than whole blood DNA for examining methylation levels in genes associated with immune function.     

Heterogeneity in White Blood Cells Has Potential to Confound DNA Methylation Measurements

Epigenetic studies are commonly conducted on DNA from tissue samples. However, tissues are ensembles of cells that may each have their own epigenetic profile, and therefore inter-individual cellular heterogeneity may compromise these studies. Here, we explore the potential for such confounding on DNA methylation measurement outcomes when using DNA from whole blood. DNA methylation was measured using pyrosequencing-based methodology in whole blood (n = 50–179) and in two white blood cell fractions (n = 20), isolated using density gradient centrifugation, in four CGIs (CpG Islands) located in genes HHEX (10 CpG sites assayed), KCNJ11 (8 CpGs), KCNQ1 (4 CpGs) and PM20D1 (7 CpGs). Cellular heterogeneity (variation in proportional white blood cell counts of neutrophils, lymphocytes, monocytes, eosinophils and basophils, counted by an automated cell counter) explained up to 40% (p<0.0001) of the inter-individual variation in whole blood DNA methylation levels in the HHEX CGI, but not a significant proportion of the variation in the other three CGIs tested. DNA methylation levels in the two cell fractions, polymorphonuclear and mononuclear cells, differed significantly in the HHEX CGI; specifically the average absolute difference ranged between 3.4–15.7 percentage points per CpG site. In the other three CGIs tested, methylation levels in the two fractions did not differ significantly, and/or the difference was more moderate. In the examined CGIs, methylation levels were highly correlated between cell fractions. In summary, our analysis detects region-specific differential DNA methylation between white blood cell subtypes, which can confound the outcome of whole blood DNA methylation measurements. Finally, by demonstrating the high correlation between methylation levels in cell fractions, our results suggest a possibility to use a proportional number of a single white blood cell type to correct for this confounding effect in analyses.     

DNA methylation analysis of bone marrow cells at diagnosis of acute lymphoblastic leukemia and at remission

To detect genes with CpG sites that display methylation patterns that are characteristic of acute lymphoblastic leukemia (ALL) cells, we compared the methylation patterns of cells taken at diagnosis from 20 patients with pediatric ALL to the methylation patterns in mononuclear cells from bone marrow of the same patients during remission and in non-leukemic control cells from bone marrow or blood. Using a custom-designed assay, we measured the methylation levels of 1,320 CpG sites in regulatory regions of 413 genes that were analyzed because they display allele-specific gene expression (ASE) in ALL cells. The rationale for our selection of CpG sites was that ASE could be the result of allele-specific methylation in the promoter regions of the genes. We found that the ALL cells had methylation profiles that allowed distinction between ALL cells and control cells. Using stringent criteria for calling differential methylation, we identified 28 CpG sites in 24 genes with recurrent differences in their methylation levels between ALL cells and control cells. Twenty of the differentially methylated genes were hypermethylated in the ALL cells, and as many as nine of them (AMICA1, CPNE7, CR1, DBC1, EYA4, LGALS8, RYR3, UQCRFS1, WDR35) have functions in cell signaling and/or apoptosis. The methylation levels of a subset of the genes were consistent with an inverse relationship with the mRNA expression levels in a large number of ALL cells from published data sets, supporting a potential biological effect of the methylation signatures and their application for diagnostic purposes.     

High density DNA methylation array with single CpG site resolution

We have developed a new generation of genome-wide DNA methylation BeadChip which allows high-throughput methylation profiling of the human genome. The new high density BeadChip can assay over 480K CpG sites and analyze twelve samples in parallel. The innovative content includes coverage of 99% of RefSeq genes with multiple probes per gene, 96% of CpG islands from the UCSC database, CpG island shores and additional content selected from whole-genome bisulfite sequencing data and input from DNA methylation experts. The well-characterized Infinium® Assay is used for analysis of CpG methylation using bisulfite-converted genomic DNA. We applied this technology to analyze DNA methylation in normal and tumor DNA samples and compared results with whole-genome bisulfite sequencing (WGBS) data obtained for the same samples. Highly comparable DNA methylation profiles were generated by the array and sequencing methods (average R² of 0.95). The ability to determine genome-wide methylation patterns will rapidly advance methylation research.     

Individual variation and longitudinal pattern of genome-wide DNA methylation from birth to the first two years of life

Prenatal development and early childhood are critical periods for establishing the tissue-specific epigenome, and may have a profound impact on health and disease in later life. However, epigenomic profiles at birth and in early childhood remain largely unexplored. The focus of this report is to examine the individual variation and longitudinal pattern of genome-wide DNA methylation levels from birth through the first two years of life in 105 Black children (59 males and 46 females) enrolled at the Boston Medical Center. We performed epigenomic mapping of cord blood at birth and venous blood samples from the same set of children within the first two years of life using Illumina Infinium Humanmethylation27 BeadChip. We observed a wide range of inter-individual variations in genome-wide methylation at each time point including lower levels at CpG islands, TSS200, 5′UTR and 1st Exon locations, but significantly higher levels in CpG shores, shelves, TSS1500, gene body and 3′UTR. We identified CpG sites with significant intra-individual longitudinal changes in the first two years of life throughout the genome. Specifically, we identified 159 CpG sites in males and 149 CpG sites in females with significant longitudinal changes defined by both statistical significance and magnitude of changes. These significant CpG sites appeared to be located within genes with important biological functions including immunity and inflammation. Further studies are needed to replicate our findings, including analysis by specific cell types, and link those individual variations and longitudinal changes with specific health outcomes in early childhood and later life.     

Individual variation and longitudinal pattern of genome-wide DNA methylation from birth to the first two years of life

Prenatal development and early childhood are critical periods for establishing the tissue-specific epigenome, and may have a profound impact on health and disease in later life. However, epigenomic profiles at birth and in early childhood remain largely unexplored. The focus of this report is to examine the individual variation and longitudinal pattern of genome-wide DNA methylation levels from birth through the first two years of life in 105 Black children (59 males and 46 females) enrolled at the Boston Medical Center. We performed epigenomic mapping of cord blood at birth and venous blood samples from the same set of children within the first two years of life using Illumina Infinium Humanmethylation27 BeadChip. We observed a wide range of inter-individual variations in genome-wide methylation at each time point including lower levels at CpG islands, TSS200, 5′UTR and 1st Exon locations, but significantly higher levels in CpG shores, shelves, TSS1500, gene body and 3′UTR. We identified CpG sites with significant intra-individual longitudinal changes in the first two years of life throughout the genome. Specifically, we identified 159 CpG sites in males and 149 CpG sites in females with significant longitudinal changes defined by both statistical significance and magnitude of changes. These significant CpG sites appeared to be located within genes with important biological functions including immunity and inflammation. Further studies are needed to replicate our findings, including analysis by specific cell types, and link those individual variations and longitudinal changes with specific health outcomes in early childhood and later life.     

Differential DNA methylation identified in the blood and retina of AMD patients

Age-related macular degeneration (AMD) is a major cause of blindness in the western world. While genetic studies have linked both common and rare variants in genes involved in regulation of the complement system to increased risk of development of AMD, environmental factors, such as smoking and nutrition, can also significantly affect the risk of developing the disease and the rate of disease progression. Since epigenetics has been implicated in mediating, in part, the disease risk associated with some environmental factors, we investigated a possible epigenetic contribution to AMD. We performed genome-wide DNA methylation profiling of blood from AMD patients and controls. No differential methylation site reached genome-wide significance; however, when epigenetic changes in and around known GWAS-defined AMD risk loci were explored, we found small but significant DNA methylation differences in the blood of neovascular AMD patients near age-related maculopathy susceptibility 2 (ARMS2), a top-ranked GWAS locus preferentially associated with neovascular AMD. The methylation level of one of the CpG sites significantly correlated with the genotype of the risk SNP rs10490924, suggesting a possible epigenetic mechanism of risk. Integrating genome-wide DNA methylation analysis of retina samples with and without AMD together with blood samples, we further identified a consistent, replicable change in DNA methylation in the promoter region of protease serine 50 (PRSS50). These methylation changes may identify sites in novel genes that are susceptible to non-genetic factors known to contribute to AMD development and progression.     

Exploring genome-wide DNA methylation profiles altered in hepatocellular carcinoma using Infinium HumanMethylation 450 BeadChips

Hepatocellular carcinoma (HCC) incidence has increased in the US and also has one of the fastest growing death rates of any cancer. The purpose of the current study was to discover novel genome-wide aberrant DNA methylation patterns in HCC tumors that are predominantly HCV-related. Infinium HumanMethylation 450K BeadChip arrays were used to examine genome-wide DNA methylation profiles in 66 pairs of HCC tumor and adjacent non-tumor tissues. After Bonferroni adjustment, a total of 130,512 CpG sites significantly differed in methylation level in tumor compared with non-tumor tissues, with 28,017 CpG sites hypermethylated and 102,495 hypomethylated in tumor tissues. Absolute tumor/non-tumor methylation differences ≥ 20% were found in 24.9% of the hypermethylated and 43.1% of the hypomethylated CpG sites; almost 10,000 CpG sites have ≥ 30% DNA methylation differences. Most (60.1%) significantly hypermethylated CpG sites are located in CpG islands, with 21.6% in CpG shores and 3.6% in shelves. In contrast, only a small proportion (8.2%) of significantly hypomethylated CpG sites are situated in islands, while most are found in open sea (60.2%), shore (17.3%) or shelf (14.3%) regions. A total of 2,568 significant CpG sites (2,441 hypermethylated and 127 hypomethylated) covering 589 genes are located within 684 differentially methylated regions defined as regions with at least two significant CpG sites displaying > 20% methylation differences in the same direction within 250-bp. The top 500 significant CpG sites can significantly distinguish HCC tumor from adjacent tissues with one misclassification. Within adjacent non-tumor tissues, we also identified 75 CpG sites significantly associated with gender, 228 with HCV infection, 17,207 with cirrhosis, and 56 with both HCV infection and cirrhosis after multiple comparisons adjustment. Aberrant DNA methylation profiles across the genome were identified in tumor tissues from US HCC cases that are predominantly related to HCV infection. These results demonstrate the significance of aberrant DNA methylation in HCC tumorigenesis.     

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.     

Whole blood DNA aberrant methylation in pancreatic adenocarcinoma shows association with the course of the disease: a pilot study

Pancreatic tumors are usually diagnosed at an advanced stage in the progression of the disease, thus reducing the survival chances of the patients. Non-invasive early detection would greatly enhance therapy and survival rates. Toward this aim, we investigated in a pilot study the power of methylation changes in whole blood as predictive markers for the detection of pancreatic tumors. We investigated methylation levels at selected CpG sites in the CpG rich regions at the promoter regions of p16, RARbeta, TNFRSF10C, APC, ACIN1, DAPK1, 3OST2, BCL2 and CD44 in the blood of 30 pancreatic tumor patients and in the blood of 49 matching controls. In addition, we studied LINE-1 and Alu repeats using degenerate amplification approach as a surrogate marker for genome-wide methylation. The site-specific methylation measurements at selected CpG sites were done by the SIRPH method. Our results show that in the patient's blood, tumor suppressor genes were slightly but significantly higher methylated at several CpG sites, while repeats were slightly less methylated compared to control blood. This was found to be significantly associated with higher risk for pancreatic ductal adenocarcinoma. Additionally, high methylation levels at TNFRSCF10C were associated with positive perineural spread of tumor cells, while higher methylation levels of TNFRSF10C and ACIN1 were significantly associated with shorter survival. This pilot study shows that methylation changes in blood could provide a promising method for early detection of pancreatic tumors. However, larger studies must be carried out to explore the clinical usefulness of a whole blood methylation based test for non-invasive early detection of pancreatic tumors.