Epigenome-wide association study reveals longitudinally stable DNA methylation differences in CD4+ T cells from children with IgE-mediated food allergy

Food allergy is mediated by a combination of genetic and environmental risk factors, potentially mediated by epigenetic mechanisms. CD4+ T-cells are key drivers of the allergic response, and may therefore harbor epigenetic variation in association with the disease phenotype. Here we retrospectively examined genome-wide DNA methylation profiles (~450?000 CpGs) from CD4+ T-cells on a birth cohort of 12 children with IgE-mediated food allergy diagnosed at 12-months, and 12 non-allergic controls. DNA samples were available at two time points, birth and 12-months. Case:control comparisons of CD4+ methylation profiles identified 179 differentially methylated probes (DMP) at 12-months and 136 DMP at birth (FDR-adjusted P value < 0.05, delta β > 0.1). Approximately 30% of DMPs were coincident with previously annotated SNPs. A total of 96 allergy-associated non-SNP DMPs were present at birth when individuals were initially disease-free, potentially implicating these loci in the causal pathway. Pathway analysis of differentially methylated genes identified several MAP kinase signaling molecules. Mass spectrometry was used to validate 15 CpG sites at 3 candidate genes. Combined analysis of differential methylation with gene expression profiles revealed gene expression differences at some but not all allergy associated differentially methylated genes. Thus, dysregulation of DNA methylation at MAPK signaling-associated genes during early CD4+ T-cell development may contribute to suboptimal T-lymphocyte responses in early childhood associated with the development of food allergy.     

Genome-Wide DNA Methylation Analysis of Systemic Lupus Erythematosus Reveals Persistent Hypomethylation of Interferon Genes and Compositional Changes to CD4+ T-cell Populations

Systemic lupus erythematosus (SLE) is an autoimmune disease with known genetic, epigenetic, and environmental risk factors. To assess the role of DNA methylation in SLE, we collected CD4+ T-cells, CD19+ B-cells, and CD14+ monocytes from 49 SLE patients and 58 controls, and performed genome-wide DNA methylation analysis with Illumina Methylation450 microarrays. We identified 166 CpGs in B-cells, 97 CpGs in monocytes, and 1,033 CpGs in T-cells with highly significant changes in DNA methylation levels (p<1×10⁻⁸) among SLE patients. Common to all three cell-types were widespread and severe hypomethylation events near genes involved in interferon signaling (type I). These interferon-related changes were apparent in patients collected during active and quiescent stages of the disease, suggesting that epigenetically-mediated hypersensitivity to interferon persists beyond acute stages of the disease and is independent of circulating interferon levels. This interferon hypersensitivity was apparent in memory, naïve and regulatory T-cells, suggesting that this epigenetic state in lupus patients is established in progenitor cell populations. We also identified a widespread, but lower amplitude shift in methylation in CD4+ T-cells (>16,000 CpGs at FDR<1%) near genes involved in cell division and MAPK signaling. These cell type-specific effects are consistent with disease-specific changes in the composition of the CD4+ population and suggest that shifts in the proportion of CD4+ subtypes can be monitored at CpGs with subtype-specific DNA methylation patterns.     

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.     

DNA methylation analysis of Homeobox genes implicates HOXB7 hypomethylation as risk factor for neural tube defects

Neural tube defects (NTDs) are common birth defects of complex etiology. Though family- and population-based studies have confirmed a genetic component, the responsible genes for NTDs are still largely unknown. Based on the hypothesis that folic acid prevents NTDs by stimulating methylation reactions, epigenetic factors, such as DNA methylation, are predicted to be involved in NTDs. Homeobox (HOX) genes play a role in spinal cord development and are tightly regulated in a spatiotemporal and collinear manner, partly by epigenetic modifications. We have quantified DNA methylation for the different HOX genes by subtracting values from a genome-wide methylation analysis using leukocyte DNA from 10 myelomeningocele (MMC) patients and 6 healthy controls. From the 1575 CpGs profiled for the 4 HOX clusters, 26 CpGs were differentially methylated (P-value < 0.05; β-difference > 0.05) between MMC patients and controls. Seventy-seven percent of these CpGs were located in the HOXA and HOXB clusters, with the most profound difference for 3 CpGs within the HOXB7 gene body. A validation case-control study including 83 MMC patients and 30 unrelated healthy controls confirmed a significant association between MMC and HOXB7 hypomethylation (-14.4%; 95% CI: 11.9–16.9%; P-value < 0.0001) independent of the MTHFR 667C>T genotype. Significant HOXB7 hypomethylation was also present in 12 unaffected siblings, each related to a MMC patient, suggestive of an epigenetic change induced by the mother. The inclusion of a neural tube formation model using zebrafish showed that Hoxb7a overexpression but not depletion resulted in deformed body axes with dysmorphic neural tube formation. Our results implicate HOXB7 hypomethylation as risk factor for NTDs and highlight the importance for future genome-wide DNA methylation analyses without preselecting candidate pathways.     

DNA methylation abnormalities in congenital heart disease

Congenital heart defects represent the most common malformation at birth, occurring also in ∼50% of individuals with Down syndrome. Congenital heart defects are thought to have multifactorial etiology, but the main causes are largely unknown. We have explored the global methylation profile of fetal heart DNA in comparison to blood DNA from control subjects: an absolute correlation with the type of tissue was detected. Pathway analysis revealed a significant enrichment of differential methylation at genes related to muscle contraction and cardiomyopathies in the developing heart DNA. We have also searched for abnormal methylation profiles on developing heart-tissue DNA of syndromic and non-syndromic congenital heart defects. On average, 3 regions with aberrant methylation were detected per sample and 18 regions were found differentially methylated between groups. Several epimutations were detected in candidate genes involved in growth regulation, apoptosis and folate pathway. A likely pathogenic hypermethylation of several intragenic sites at the MSX1 gene, involved in outflow tract morphogenesis, was found in a fetus with isolated heart malformation. In addition, hypermethylation of the GATA4 gene was present in fetuses with Down syndrome with or without congenital heart defects, as well as in fetuses with isolated heart malformations. Expression deregulation of the abnormally methylated genes was detected. Our data indicate that epigenetic alterations of relevant genes are present in developing heart DNA in fetuses with both isolated and syndromic heart malformations. These epimutations likely contribute to the pathogenesis of the malformation by cis-acting effects on gene expression.     

An Epigenetic Signature in Peripheral Blood Associated with the Haplotype on 17q21.31, a Risk Factor for Neurodegenerative Tauopathy

Little is known about how changes in DNA methylation mediate risk for human diseases including dementia. Analysis of genome-wide methylation patterns in patients with two forms of tau-related dementia – progressive supranuclear palsy (PSP) and frontotemporal dementia (FTD) – revealed significant differentially methylated probes (DMPs) in patients versus unaffected controls. Remarkably, DMPs in PSP were clustered within the 17q21.31 region, previously known to harbor the major genetic risk factor for PSP. We identified and replicated a dose-dependent effect of the risk-associated H1 haplotype on methylation levels within the region in blood and brain. These data reveal that the H1 haplotype increases risk for tauopathy via differential methylation at that locus, indicating a mediating role for methylation in dementia pathophysiology.     

The overexpression of GRASP might inhibit cell proliferation and invasion in hepatocellular carcinoma

This study aimed to validate the methylation of key genes in hepatocellular carcinoma (HCC) screened by bioinformatics analysis and explore whether they affected HCC cell proliferation, migration, and invasion. Using The Cancer Genome Atlas (TCGA) database, HCC-related differentially methylated positions (DMPs) were screened, genes corresponding to DMPs were selected, and prognosis-related genes were identified. A representative DMP was used to divide the DMPs into hyper- and hypomethylated groups. Expression of key genes in cell lines was detected using quantitative real-time polymerase chain reaction and western blot analysis. After treatment of HepG2 cells with 5-Aza-2′-deoxycytidine (5-Aza-DC), gene expression was observed. Bisulfite sequencing PCR assay was used to detect methylation frequency. Overexpressed GRASP lentiviral vectors were constructed to analyze their influence on cell proliferation, migration, and invasion using cell counting kit-8 and transwell assays. Forty-three HCC prognosis-related genes were screened using the TCGA database. cg00249511 (SCT) was used to divide the DMPs into hyper- and hypomethylated groups, distinguishing between high- and low-risk samples. The prognosis survival model constructed using 12 genes revealed the prognosis type. GRASP messenger RNA was downregulated in HepG2 and upregulated after 5-Aza-DC treatment. In HCC tissues, methylation frequency of GRASP was upregulated. GRASP overexpression inhibited HepG2 cell proliferation, invasion, and G-CSFR expression. Thus, GRASP might be a prognosis-related gene controlled by methylation.     

Allele-Specific Transcriptome and Methylome Analysis Reveals Stable Inheritance and Cis-Regulation of DNA Methylation in Nasonia

Gene expression divergence between closely related species could be attributed to both cis- and trans- DNA sequence changes during evolution, but it is unclear how the evolutionary dynamics of epigenetic marks are regulated. In eutherian mammals, biparental DNA methylation marks are erased and reset during gametogenesis, resulting in paternal or maternal imprints, which lead to genomic imprinting. Whether DNA methylation reprogramming exists in insects is not known. Wasps of the genus Nasonia are non-social parasitoids that are emerging as a model for studies of epigenetic processes in insects. In this study, we quantified allele-specific expression and methylation genome-wide in Nasonia vitripennis and Nasonia giraulti and their reciprocal F1 hybrids. No parent-of-origin effect in allelic expression was found for >8,000 covered genes, suggesting a lack of genomic imprinting in adult Nasonia. As we expected, both significant cis- and trans- effects are responsible for the expression divergence between N. vitripennis and N. giraulti. Surprisingly, all 178 differentially methylated genes are also differentially methylated between the two alleles in F1 hybrid offspring, recapitulating the parental methylation status with nearly 100% fidelity, indicating the presence of strong cis-elements driving the target of gene body methylation. In addition, we discovered that total and allele-specific expression are positively correlated with allele-specific methylation in a subset of the differentially methylated genes. The 100% cis-regulation in F1 hybrids suggests the methylation machinery is conserved and DNA methylation is targeted by cis features in Nasonia. The lack of genomic imprinting and parent-of-origin differentially methylated regions in Nasonia, together with the stable inheritance of methylation status between generations, suggests either a cis-regulatory motif for methylation at the DNA level or highly stable inheritance of an epigenetic signal in Nasonia.     

Identification of a DNA Methylome Profile of Esophageal Squamous Cell Carcinoma and Potential Plasma Epigenetic Biomarkers for Early Diagnosis

DNA methylation is a critical epigenetic mechanism involved in key cellular processes. Its deregulation has been linked to many human cancers including esophageal squamous cell carcinoma (ESCC). This study was designed to explore the whole methylation status of ESCC and to identify potential plasma biomarkers for early diagnosis. We used Infinium Methylation 450k array to analyze ESCC tissues (n = 4), paired normal surrounding tissues (n = 4) and normal mucosa from healthy individuals (n = 4), and combined these with gene expression data from the GEO database. One hundred and sixty eight genes had differentially methylated CpG sites in their promoter region and a gene expression pattern inverse to the direction of change in DNA methylation. These genes were involved in several cancer-related pathways. Three genes were validated in additional 42 ESCC tissues and paired normal surrounding tissues. The methylation frequency of EPB41L3, GPX3, and COL14A1 were higher in tumor tissues than in normal surrounding tissues (P<0.017). The higher methylation frequency of EPB41l3 was correlated with large tumor size (P = 0.044) and advanced pT tumor stage (P = 0.001). The higher methylation frequency of GPX3 and COL14A1 were correlated with advanced pN tumor stage (P = 0.001 and P<0.001). The methylation of EPB41L3, GPX3, and COL14A1 genes were only found in ESCC patients'' plasma, but not in normal individuals upon testing 42 ESCC patients and 50 healthy individuals. Diagnostic sensitivity was increased when methylation of any of the 3 genes were counted (64.3% sensitivity and 100% specificity). These differentially methylated genes in plasma may be used as biomarkers for early diagnosis of ESCC.     

Genome-wide DNA methylation profiling of recurrent and non-recurrent chordomas

Chordomas are an aggressive rare type of malignant bone tumors arising from the remnant of the notochord. Chordomas occur mainly in vertebral bones and account for 1–4% of malignant bone tumors. Management and treatment of chordomas are difficult as they are resistant to conventional chemotherapy; therefore, they are mainly treated with surgery and radiation therapy. In this study, we performed DNA methylation profiling of 26 chordomas and normal nucleus pulposus samples plus UCH-1 chordoma cell line using the Illumina Infinium HumanMethylation450 BeadChips. Combined bisulfite restriction analysis and bisulfite sequencing was used to confirm the methylation data. Gene expression was analyzed using RT-PCR before and after 5-aza-2’-deoxycytidine (5-azaDC) treatment of chordoma cell lines. Analysis of the HumanMethylation450 BeadChip data led to the identification of 8,819 loci (2.9%) that were significantly differentially methylated (>0.2 average β-value difference) between chordomas and nucleus pulposus samples (adjusted P < 0.05). Among these, 5,868 probes (66.5%) were hypomethylated, compared to 2,951 (33.5%) loci that were hypermethylated in chordomas compared to controls. From the 2,951 differentially hypermethylated probes, 33.3% were localized in the promoter region (982 probes) and, among these, 104 probes showed cancer-specific hypermethylation. Ingenuity Pathway Analysis indicates that the cancer-specific differentially methylated loci are involved in various networks including cancer disease, nervous system development and function, cell death and survival, cellular growth, cellular development, and proliferation. Furthermore, we identified a subset of probes that were differentially methylated between recurrent and non-recurrent chordomas. BeadChip methylation data was confirmed for these genes and gene expression was shown to be upregulated in methylated chordoma cell lines after treatment with 5-azaDC. Understanding epigenetic changes in chordomas may provide insights into chordoma tumorigenesis and development of epigenetic biomarkers.     

DNA methylation and gene expression changes in monozygotic twins discordant for psoriasis: identification of epigenetically dysregulated genes

Monozygotic (MZ) twins do not show complete concordance for many complex diseases; for example, discordance rates for autoimmune diseases are 20%-80%. MZ discordance indicates a role for epigenetic or environmental factors in disease. We used MZ twins discordant for psoriasis to search for genome-wide differences in DNA methylation and gene expression in CD4(+) and CD8(+) cells using Illumina's HumanMethylation27 and HT-12 expression assays, respectively. Analysis of these data revealed no differentially methylated or expressed genes between co-twins when analyzed separately, although we observed a substantial amount of small differences. However, combined analysis of DNA methylation and gene expression identified genes where differences in DNA methylation between unaffected and affected twins were correlated with differences in gene expression. Several of the top-ranked genes according to significance of the correlation in CD4(+) cells are known to be associated with psoriasis. Further, gene ontology (GO) analysis revealed enrichment of biological processes associated with the immune response and clustering of genes in a biological pathway comprising cytokines and chemokines. These data suggest that DNA methylation is involved in an epigenetic dysregulation of biological pathways involved in the pathogenesis of psoriasis. This is the first study based on data from MZ twins discordant for psoriasis to detect epigenetic alterations that potentially contribute to development of the disease.     

Whole-Genome Saliva and Blood DNA Methylation Profiling in Individuals with a Respiratory Allergy

The etiology of respiratory allergies (RA) can be partly explained by DNA methylation changes caused by adverse environmental and lifestyle factors experienced early in life. Longitudinal, prospective studies can aid in the unravelment of the epigenetic mechanisms involved in the disease development. High compliance rates can be expected in these studies when data is collected using non-invasive and convenient procedures. Saliva is an attractive biofluid to analyze changes in DNA methylation patterns. We investigated in a pilot study the differential methylation in saliva of RA (n = 5) compared to healthy controls (n = 5) using the Illumina Methylation 450K BeadChip platform. We evaluated the results against the results obtained in mononuclear blood cells from the same individuals. Differences in methylation patterns from saliva and mononuclear blood cells were clearly distinguishable (P_Adj<0.001 and |Δβ|>0.2), though the methylation status of about 96% of the cg-sites was comparable between peripheral blood mononuclear cells and saliva. When comparing RA cases with healthy controls, the number of differentially methylated sites (DMS) in saliva and blood were 485 and 437 (P<0.05 and |Δβ|>0.1), respectively, of which 216 were in common. The methylation levels of these sites were significantly correlated between blood and saliva. The absolute levels of methylation in blood and saliva were confirmed for 3 selected DMS in the PM20D1, STK32C, and FGFR2 genes using pyrosequencing analysis. The differential methylation could only be confirmed for DMS in PM20D1 and STK32C genes in saliva. We show that saliva can be used for genome-wide methylation analysis and that it is possible to identify DMS when comparing RA cases and healthy controls. The results were replicated in blood cells of the same individuals and confirmed by pyrosequencing analysis. This study provides proof-of-concept for the applicability of saliva-based whole-genome methylation analysis in the field of respiratory allergy.     

Methylomic analysis of monozygotic twins discordant for childhood psychotic symptoms

Childhood psychotic symptoms are associated with increased rates of schizophrenia, other psychiatric disorders, and suicide attempts in adulthood; thus, elucidating early risk indicators is crucial to target prevention efforts. There is considerable discordance for psychotic symptoms between monozygotic twins, indicating that child-specific non-genetic factors must be involved. Epigenetic processes may constitute one of these factors and have not yet been investigated in relation to childhood psychotic symptoms. Therefore, this study explored whether differences in DNA methylation at age 10 were associated with monozygotic twin discordance for psychotic symptoms at age 12. The Environmental Risk (E-Risk) Longitudinal Twin Study cohort of 2,232 children (1,116 twin pairs) was assessed for age-12 psychotic symptoms and 24 monozygotic twin pairs discordant for symptoms were identified for methylomic comparison. Children provided buccal samples at ages 5 and 10. DNA was bisulfite modified and DNA methylation was quantified using the Infinium HumanMethylation450 array. Differentially methylated positions (DMPs) associated with psychotic symptoms were subsequently tested in post-mortem prefrontal cortex tissue from adult schizophrenia patients and age-matched controls. Site-specific DNA methylation differences were observed at age 10 between monozygotic twins discordant for age-12 psychotic symptoms. Similar DMPs were not found at age 5. The top-ranked psychosis-associated DMP (cg23933044), located in the promoter of the C5ORF42 gene, was also hypomethylated in post-mortem prefrontal cortex brain tissue from schizophrenia patients compared to unaffected controls. These data tentatively suggest that epigenetic variation in peripheral tissue is associated with childhood psychotic symptoms and may indicate susceptibility to schizophrenia and other mental health problems.