Tissue-specific patterns of allelically-skewed DNA methylation

Metastable epialleles in humans

First identified in isogenic mice, metastable epialleles (MEs) are loci where the extent of DNA methylation (DNAm) is variable between individuals but correlates across tissues derived from different germ layers within a given individual. This property, termed systemic interindividual variation (SIV), is attributed to stochastic methylation establishment before germ layer differentiation. Evidence suggests that some putative human MEs are sensitive to environmental exposures in early development. In this review we introduce key concepts pertaining to human MEs, describe methods used to identify MEs in humans, and review their genomic features. We also highlight studies linking DNAm at putative human MEs to early environmental exposures and postnatal (including disease) phenotypes.

Association between methylation in the promoter region of the GAD2 gene and opioid use disorder

DNA methylation is one of the epigenetic mechanisms involved in opioid use disorder. GAD2 is a key catalyticase in gamma amino butyric acid (GABA) synthesis from glutamate, that is implicated in opioid-induced rewarding effect. To reveal the relationship and the underlying mechanism between GAD2 gene methylation and opioid use disorder, we first examined and compared the methylation levels in the promoter region of the GAD2 gene in peripheral blood between 120 patients with opioid use disorder and 110 healthy controls by using a targeted approach. A diagnostic model with methylation biomarkers was established to distinguish opioid use disorder and healthy control groups. Correlations between methylation levels in the promoter region of the GAD2 gene and the duration and dosage of opioid use were then determined. Finally, the transcription factors that potentially bind to the target sequences including the detected CpG sites were predicted with the JASPAR database. Our results demonstrated that hypermethylation in the promoter region of the GAD2 gene was associated with opioid use disorder. A diagnostic model based on 10 methylation biomarkers could distinguish the opioid use disorder and healthy control groups. Several correlations between methylation levels in the GAD2 gene promoter and the duration and dosage of opioid use were observed. Transcription factors TFAP2A, Arnt and Runx1 were predicted to bind to the target sequences including several CpG sites detected in the present study in the GAD2 gene promoter. Our findings highlight and extend the role of DNA methylation in the GAD2 gene in opioid use disorder.

A comparative methylome analysis reveals conservation and divergence of DNA methylation patterns and functions in vertebrates

Background

Cytosine DNA methylation is a heritable epigenetic mark present in most eukaryotic groups. While the patterns and functions of DNA methylation have been extensively studied in mouse and human, their conservation in other vertebrates remains poorly explored. In this study, we interrogated the distribution and function of DNA methylation in primary fibroblasts of seven vertebrate species including bio-medical models and livestock species (human, mouse, rabbit, dog, cow, pig, and chicken).

Results

Our data highlight both divergence and conservation of DNA methylation patterns and functions. We show that the chicken genome is hypomethylated compared to other vertebrates. Furthermore, compared to mouse, other species show a higher frequency of methylation of CpG-rich DNA. We reveal the conservation of large unmethylated valleys and patterns of DNA methylation associated with X-chromosome inactivation through vertebrate evolution and make predictions of conserved sets of imprinted genes across mammals. Finally, using chemical inhibition of DNA methylation, we show that the silencing of germline genes and endogenous retroviruses (ERVs) are conserved functions of DNA methylation in vertebrates.

Conclusions

Our data highlight conserved properties of DNA methylation in vertebrate genomes but at the same time point to differences between mouse and other vertebrate species.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12915-022-01270-x.

Allele-specific DNA methylation maps in monozygotic twins discordant for psychiatric disorders reveal that disease-associated switching at the EIPR1 regulatory loci modulates neural function

The non-Mendelian features of phenotypic variations within monozygotic twins are likely complicated by environmental modifiers of genetic effects that have yet to be elucidated. Here, we performed methylome and genome analyses of blood DNA from psychiatric disorder-discordant monozygotic twins to study how allele-specific methylation (ASM) mediates phenotypic variations. We identified that thousands of genetic variants with ASM imbalances exhibit phenotypic variation-associated switching at regulatory loci. These ASMs have plausible causal associations with psychiatric disorders through effects on interactions between transcription factors, DNA methylations, and other epigenomic markers and then contribute to dysregulated gene expression, which eventually increases disease susceptibility. Moreover, we also experimentally validated the model that the rs4854158 alternative C allele at an ASM switching regulatory locus of EIPR1 encoding endosome-associated recycling protein-interacting protein 1, is associated with demethylation and higher RNA expression and shows lower TF binding affinities in unaffected controls. An epigenetic ASM switching induces C allele hypermethylation and then recruits repressive Polycomb repressive complex 2 (PRC2), reinforces trimethylation of lysine 27 on histone 3 and inhibits its transcriptional activity, thus leading to downregulation of EIPR1 in schizophrenia. Moreover, disruption of rs4854158 induces gain of EIPR1 function and promotes neural development and vesicle trafficking. Our study provides a powerful framework for identifying regulatory risk variants and contributes to our understanding of the interplay between genetic and epigenetic variants in mediating psychiatric disorder susceptibility.

DAMEfinder: a method to detect differential allele-specific methylation

Background

DNA methylation is a highly studied epigenetic signature that is associated with regulation of gene expression, whereby genes with high levels of promoter methylation are generally repressed. Genomic imprinting occurs when one of the parental alleles is methylated, i.e., when there is inherited allele-specific methylation (ASM). A special case of imprinting occurs during X chromosome inactivation in females, where one of the two X chromosomes is silenced, to achieve dosage compensation between the sexes. Another more widespread form of ASM is sequence dependent (SD-ASM), where ASM is linked to a nearby heterozygous single nucleotide polymorphism (SNP).

Results

We developed a method to screen for genomic regions that exhibit loss or gain of ASM in samples from two conditions (treatments, diseases, etc.). The method relies on the availability of bisulfite sequencing data from multiple samples of the two conditions. We leverage other established computational methods to screen for these regions within a new R package called DAMEfinder. It calculates an ASM score for all CpG sites or pairs in the genome of each sample, and then quantifies the change in ASM between conditions. It then clusters nearby CpG sites with consistent change into regions. In the absence of SNP information, our method relies only on reads to quantify ASM. This novel ASM score compares favorably to current methods that also screen for ASM. Not only does it easily discern between imprinted and non-imprinted regions, but also females from males based on X chromosome inactivation. We also applied DAMEfinder to a colorectal cancer dataset and observed that colorectal cancer subtypes are distinguishable according to their ASM signature. We also re-discover known cases of loss of imprinting.

Conclusion

We have designed DAMEfinder to detect regions of differential ASM (DAMEs), which is a more refined definition of differential methylation, and can therefore help in breaking down the complexity of DNA methylation and its influence in development and disease.

Influence of pre-pregnancy body mass index (p-BMI) and gestational weight gain (GWG) on DNA methylation and protein expression of obesogenic genes in umbilical vein

Understanding the regulatory mechanisms that affect obesogenic genes expression in newborns is essential for early prevention efforts, but they remain unclear. Our study aimed to explore whether the maternal p-BMI and GWG were associated with regulatory single-locus DNA methylation in selected obesogenic genes. For this purpose, DNA methylation was assayed by Methylation-Sensitive High Resolution Melting (MS-HRM) technique and Sanger allele-bisulfite sequencing in fifty samples of umbilical vein to evaluate glucosamine-6-phosphate deaminase 2 (GNPDA2), peroxisome proliferator-activated receptor gamma coactivator 1 alpha (PGC1α), and leptin receptor (LEPR) genes. Correlations between DNA methylation levels and indicators of maternal nutritional status were carried out. Western blotting was used to evaluate protein expression in extracts of the same samples. Results indicated that GNPDA2 and PGC1α genes have the same level of DNA methylation in all samples; however, a differential DNA methylation of LEPR gene promoter was found, correlating it with GWG and this correlation is unaffected by maternal age or unhealthy habits. Furthermore, leptin receptor (Lep-Rb) was upregulated in samples that showed the lowest levels of DNA methylation. This study highlights the association between poor GWG and adjustments on obesogenic genes expression in newborn tissues with potential consequences for development of obesity in the future.

Genome-wide allele-specific methylation is enriched at gene regulatory regions in a multi-generation pedigree from the Norfolk Island isolate

Background

Allele-specific methylation (ASM) occurs when DNA methylation patterns exhibit asymmetry among alleles. ASM occurs at imprinted loci, but its presence elsewhere across the human genome is indicative of wider importance in terms of gene regulation and disease risk. Here, we studied ASM by focusing on blood-based DNA collected from 24 subjects comprising a 3-generation pedigree from the Norfolk Island genetic isolate. We applied a genome-wide bisulphite sequencing approach with a genotype-independent ASM calling method to map ASM across the genome. Regions of ASM were then tested for enrichment at gene regulatory regions using Genomic Association Test (GAT) tool.

Results

In total, we identified 1.12 M CpGs of which 147,170 (13%) exhibited ASM (P ≤ 0.05). When including contiguous ASM signal spanning ≥ 2 CpGs, this condensed to 12,761 ASM regions (AMRs). These AMRs tagged 79% of known imprinting regions and most (98.1%) co-localised with known single nucleotide variants. Notably, miRNA and lncRNA showed a 3.3- and 1.8-fold enrichment of AMRs, respectively (P < 0.005). Also, the 5′ UTR and start codons each showed a 3.5-fold enrichment of AMRs (P < 0.005). There was also enrichment of AMRs observed at subtelomeric regions of many chromosomes. Five out of 11 large AMRs localised to the protocadherin cluster on chromosome 5.

Conclusions

This study shows ASM extends far beyond genomic imprinting in humans and that gene regulatory regions are hotspots for ASM. Future studies of ASM in pedigrees should help to clarify transgenerational inheritance patterns in relation to genotype and disease phenotypes.

A genomic atlas of systemic interindividual epigenetic variation in humans

BACKGROUND

DNA methylation is thought to be an important determinant of human phenotypic variation, but its inherent cell type specificity has impeded progress on this question. At exceptional genomic regions, interindividual variation in DNA methylation occurs systemically. Like genetic variants, systemic interindividual epigenetic variants are stable, can influence phenotype, and can be assessed in any easily biopsiable DNA sample. We describe an unbiased screen for human genomic regions at which interindividual variation in DNA methylation is not tissue-specific.

RESULTS

For each of 10 donors from the NIH Genotype-Tissue Expression (GTEx) program, CpG methylation is measured by deep whole-genome bisulfite sequencing of genomic DNA from tissues representing the three germ layer lineages: thyroid (endoderm), heart (mesoderm), and brain (ectoderm). We develop a computational algorithm to identify genomic regions at which interindividual variation in DNA methylation is consistent across all three lineages. This approach identifies 9926 correlated regions of systemic interindividual variation (CoRSIVs). These regions, comprising just 0.1% of the human genome, are inter-correlated over long genomic distances, associated with transposable elements and subtelomeric regions, conserved across diverse human ethnic groups, sensitive to periconceptional environment, and associated with genes implicated in a broad range of human disorders and phenotypes. CoRSIV methylation in one tissue can predict expression of associated genes in other tissues.

CONCLUSIONS

In addition to charting a previously unexplored molecular level of human individuality, this atlas of human CoRSIVs provides a resource for future population-based investigations into how interindividual epigenetic variation modulates risk of disease.

Epigenetics in eating disorders: a systematic review

Eating disorders are complex heritable conditions influenced by both genetic and environmental factors. Given the progress of genomic discovery in anorexia nervosa, with the identification of the first genome-wide significant locus, as well as animated discussion of epigenetic mechanisms in linking environmental factors with disease onset, our goal was to conduct a systematic review of the current body of evidence on epigenetic factors in eating disorders to inform future directions in this area. Following PRISMA guidelines, two independent authors conducted a search within PubMed and Web of Science and identified 18 journal articles and conference abstracts addressing anorexia nervosa (n = 13), bulimia nervosa (n = 6), and binge-eating disorder (n = 1), published between January 2003 and October 2017. We reviewed all articles and included a critical discussion of field-specific methodological considerations. The majority of epigenetic analyses of eating disorders investigated methylation at candidate genes (n = 13), focusing on anorexia and bulimia nervosa in very small samples with considerable sample overlap across published studies. Three studies used microarray-based technologies to examine DNA methylation across the genome of anorexia nervosa and binge-eating disorder patients. Overall, results were inconclusive and were primarily exploratory in nature. The field of epigenetics in eating disorders remains in its infancy. We encourage the scientific community to apply methodologically sound approaches using genome-wide designs including epigenome-wide association studies (EWAS), to increase sample sizes, and to broaden the focus to include all eating disorder types.

Crosstalk of Genetic Variants, Allele-Specific DNA Methylation, and Environmental Factors for Complex Disease Risk

Over the past decades, genome-wide association studies (GWAS) have identified thousands of phenotype-associated DNA sequence variants for potential explanations of inter-individual phenotypic differences and disease susceptibility. However, it remains a challenge for translating the associations into causative mechanisms for complex diseases, partially due to the involved variants in the noncoding regions and the inconvenience of functional studies in human population samples. So far, accumulating evidence has suggested a complex crosstalk among genetic variants, allele-specific binding of transcription factors (ABTF), and allele-specific DNA methylation patterns (ASM), as well as environmental factors for disease risk. This review aims to summarize the current studies regarding the interactions of the aforementioned factors with a focus on epigenetic insights. We present two scenarios of single nucleotide polymorphisms (SNPs) in coding regions and non-coding regions for disease risk, via potentially impacting epigenetic patterns. While a SNP in a coding region may confer disease risk via altering protein functions, a SNP in non-coding region may cause diseases, via SNP-altering ABTF, ASM, and allele-specific gene expression (ASE). The allelic increases or decreases of gene expression are key for disease risk during development. Such ASE can be achieved via either a "SNP-introduced ABTF to ASM" or a "SNP-introduced ASM to ABTF." Together with our additional in-depth review on insulator CTCF, we are convinced to propose a working model that the small effect of a SNP acts through altered ABTF and/or ASM, for ASE and eventual disease outcome (named as a "SNP intensifier" model). In summary, the significance of complex crosstalk among genetic factors, epigenetic patterns, and environmental factors requires further investigations for disease susceptibility.

Neuroepigenetics of Schizophrenia

Schizophrenia is a complex disorder of the brain, where genetic variants explain only a portion of risk. Neuroepigenetic mechanisms may explain the remaining share of risk, as well as the transition from susceptibility to the actual disease. Here, we discuss the most recent findings in the field of brain epigenetics applied to the study of schizophrenia. Methylome studies have found several candidates exhibiting methylation modifications in association with the disorder, but genes affected do not always overlap. Notably, these studies converge in that genes within the schizophrenia risk loci or genes differentially methylated in patients affected with the disorder are dynamically regulated during early life. They also imply that schizophrenia-associated genetic variation may affect DNA methylation in fetal and adult brains. Histone modifications may help mediating the effect of genetic risk variants associated with schizophrenia, and regulating chromatin higher-order structure. The 3D-organization of chromatin in the brain creates physical interactions within chromosomes, so that schizophrenia-associated genetic variants can be linked with genes distant from their loci; this suggests that chromatin conformation matters in the mechanism of risk for the disorder. Non-coding RNAs provide a novel and complex mechanism of gene regulation potentially significant for schizophrenia, as proposed by research on specific microRNAs and long non-coding RNAs (lncRNAs). Finally, a recent study in epitranscriptomics identifies RNA methylation as a further epigenetic mechanism active in human brain and specifically in a portion of the transcriptome associated with schizophrenia susceptibility. These findings indicate that, as expected from the complexity of the brain and its development, several epigenetic mechanisms may intervene in the etiopathogenesis of schizophrenia. An understanding of their roles calls for research approaches integrating the investigation of different epigenetic mechanisms and of environmental and genetic risk, in the context of development.

Maternal 5mCpG Imprints at the PARD6G-AS1 and GCSAML Differentially Methylated Regions Are Decoupled From Parent-of-Origin Expression Effects in Multiple Human Tissues

A hallmark of imprinted genes in mammals is the occurrence of parent-of-origin-dependent asymmetry of DNA cytosine methylation (5^mC) of alleles at CpG islands (CGIs) in their promoter regions. This 5^mCpG asymmetry between the parental alleles creates allele-specific imprinted differentially methylated regions (iDMRs). iDMRs are often coupled to the transcriptional repression of the methylated allele and the activation of the unmethylated allele in a tissue-specific, developmental-stage-specific and/or isoform-specific fashion. iDMRs function as regulatory platforms, built through the recruitment of chemical modifications to histones to achieve differential, parent-of-origin-dependent chromatin segmentation states. Here, we used a comparative computational data mining approach to identify 125 novel constitutive candidate iDMRs that integrate the maximal number of allele-specific methylation region records overlapping CGIs in human methylomes. Twenty-nine candidate iDMRs display gametic 5^mCpG asymmetry, and another 96 are candidate secondary iDMRs. We established the maternal origin of the 5^mCpG imprints of one gametic (PARD6G-AS1) and one secondary (GCSAML) iDMRs. We also found a constitutively hemimethylated, nonimprinted domain at the PWWP2AP1 promoter CGI with oocyte-derived methylation asymmetry. Given that the 5^mCpG level at the iDMRs is not a sufficient criterion to predict active or silent locus states and that iDMRs can regulate genes from a distance of more than 1 Mb, we used RNA-Seq experiments from the Genotype-Tissue Expression project and public archives to assess the transcriptional expression profiles of SNPs across 4.6 Mb spans around the novel maternal iDMRs. We showed that PARD6G-AS1 and GCSAML are expressed biallelically in multiple tissues. We found evidence of tissue-specific monoallelic expression of ZNF124 and OR2L13, located 363 kb upstream and 419 kb downstream, respectively, of the GCSAML iDMR. We hypothesize that the GCSAML iDMR regulates the tissue-specific, monoallelic expression of ZNF124 but not of OR2L13. We annotated the non-coding epigenomic marks in the two maternal iDMRs using data from the Roadmap Epigenomics project and showed that the PARD6G-AS1 and GCSAML iDMRs achieve contrasting activation and repression chromatin segmentations. Lastly, we found that the maternal 5^mCpG imprints are perturbed in several hematopoietic cancers. We conclude that the maternal 5^mCpG imprints at PARD6G-AS1 and GCSAML iDMRs are decoupled from parent-of-origin transcriptional expression effects in multiple tissues.

Reelin (RELN) DNA methylation in the peripheral blood of schizophrenia

The epigenetic changes of RELN that are involved in the development of dopaminergic neurons may fit the developmental theory of schizophrenia. However, evidence regarding the association of RELN DNA methylation with schizophrenia is far from sufficient, as studies have only been conducted on a few limited brain samples. As DNA methylation in the peripheral blood may mirror the changes taking place in the brain, the use of peripheral blood for a DNA methylation study in schizophrenia is feasible due to the scarcity of brain samples. Therefore, the aim of our study was to examine the relationship of DNA methylation levels of RELN promoters with schizophrenia using genomic DNA derived from the peripheral blood of patients with the disorder. The case control studies consisted of 110 schizophrenia participants and 122 healthy controls who had been recruited from the same district. After bisufhite conversion, the methylation levels of the DNA samples were calculated based on their differences of the Cq values assayed using the highly sensitive real-time MethyLight TaqMan^® procedure. A significantly higher level of methylation of the RELN promoter was found in patients with schizophrenia compared to controls (p = 0.005) and also in males compared with females (p = 0.004). Subsequently, the RELN expression of the methylated group was 25 fold less than that of the non-methylated group. Based upon the assumption of parallel methylation changes in the brain and peripheral blood, we concluded that RELN DNA methylation might contribute to the pathogenesis of schizophrenia. However, the definite effects of methylation on RELN function during development and also in adult life still require further elaboration.

MeDIP-seq reveals the features of mitochondrial genomic methylation in immature testis of Chinese mitten crab Eriocheir sinensis

In this study, the methylation of mitochondrial genome in the immature testis of Chinese mitten crab Eriocheir sinensis of the Yangtze River system was determined for the first time using MeDIP-seq. Our methylated DNA fragments covered more than 99% of the mitochondrial genome in E. sinensis loaded from GenBank. There were 8 mutated bases and 42 SNPs in the crab mitochondrial genome. The methylation presented in all genes as well as in an A + T region, but less in intergenic regions in the mitochondrial genome. However, the level of methylation of most genes coding proteins and the A + T region were high. But, the majority of genes encoding tRNAs were hypomethylated, and both the rRNA genes also showed methylation of low or median frequency. Especially, the level of methylation of the intergenic regions is the lowest. Those features indicated that the methylation of DNA may play an important role in gene expressing regulation in the mitochondrial genome of immature testis in E. sinensis.

Epigenetic Research in Neuropsychiatric Disorders: the "Tissue Issue"

PURPOSE OF REVIEW

Evidence has linked neuropsychiatric disorders with epigenetic marks as either a biomarker of disease, biomarker of exposure, or mechanism of disease processes. Neuropsychiatric epidemiologic studies using either target brain tissue or surrogate blood tissue each have methodological challenges and distinct advantages.

RECENT FINDINGS

Brain tissue studies are challenged by small sample sizes of cases and controls, incomplete phenotyping, post-mortem timing, and cellular heterogeneity, but the use of a primary disease relevant tissue is critical. Blood-based studies have access to much larger sample sizes and more replication opportunities, as well as the potential for longitudinal measurements, both prior to onset and during the course of treatments. Yet, blood studies also are challenged by cell-type heterogeneity, and many question the validity of using peripheral tissues as a brain biomarker. Emerging evidence suggests that these limitations to blood-based epigenetic studies are surmountable, but confirmation in target tissue remains important.

SUMMARY

Epigenetic mechanisms have the potential to help elucidate biology connecting experiential risk factors with neuropsychiatric disease manifestation. Cross-tissue studies as well as advanced epidemiologic methods should be employed to more effectively conduct neuropsychiatric epigenetic research.

How to make DNA methylome wide association studies more powerful

Genome-wide association studies had a troublesome adolescence, while researchers increased statistical power, in part by increasing subject numbers. Interrogating the interaction of genetic and environmental influences raised new challenges of statistical power, which were not easily bested by the addition of subjects. Screening the DNA methylome offers an attractive alternative as methylation can be thought of as a proxy for the combined influences of genetics and environment. There are statistical challenges unique to DNA methylome data and also multiple features, which can be exploited to increase power. We anticipate the development of DNA methylome association study designs and new analytical methods, together with integration of data from other molecular species and other studies, which will boost statistical power and tackle causality. In this way, the molecular trajectories that underlie disease development will be uncovered.