Allele-specific methylation in the human genome: implications for genetic studies of complex disease

Epigenetics of prenatal stress in humans: the current research landscape

Fetal exposure to prenatal stress can have significant consequences on short- and long-term health. Epigenetic mechanisms, especially DNA methylation (DNAm), are a possible process how these adverse environmental events could be biologically embedded. We evaluated candidate gene as well as epigenome-wide association studies associating prenatal stress and DNAm changes in peripheral tissues; however, most of these findings lack robust replication. Prenatal stress-associated epigenetic changes have also been linked to child health including internalizing problems, neurobehavioral outcomes and stress reactivity. Future studies should focus on refined measurement and definition of prenatal stress and its timing, ideally also incorporating genomic as well as longitudinal information. This will provide further opportunities to enhance our understanding of the biological embedding of prenatal stress exposure.

Transition of allele-specific DNA hydroxymethylation at regulatory loci is associated with phenotypic variation in monozygotic twins discordant for psychiatric disorders

BACKGROUND

Major psychiatric disorders such as schizophrenia (SCZ) and bipolar disorder (BPD) are complex genetic mental illnesses. Their non-Mendelian features, such as those observed in monozygotic twins discordant for SCZ or BPD, are likely complicated by environmental modifiers of genetic effects. 5-Hydroxymethylcytosine (5hmC) is an important epigenetic mark in gene regulation, and whether it is linked to genetic variants that contribute to non-Mendelian features remains largely unexplored.

METHODS

We combined the 5hmC-selective chemical labeling method (5hmC-seq) and whole-genome sequencing (WGS) analysis of peripheral blood DNA obtained from monozygotic (MZ) twins discordant for SCZ or BPD to identify allelic imbalances in hydroxymethylome maps, and examined association of allele-specific hydroxymethylation (AShM) transition with disease susceptibility based on Bayes factors (BF) derived from the Bayesian generalized additive linear mixed model. We then performed multi-omics integrative analysis to determine the molecular pathogenic basis of those AShM sites. We finally employed luciferase reporter, CRISPR/Cas9 technology, electrophoretic mobility shift assay (EMSA), chromatin immunoprecipitation (ChIP), PCR, FM4-64 imaging analysis, and RNA sequencing to validate the function of interested AShM sites in the human neuroblastoma SK-N-SH cells and human embryonic kidney 293T (HEK293T) cells.

RESULTS

We identified thousands of genetic variants associated with AShM imbalances that exhibited phenotypic variation-associated AShM changes at regulatory loci. These AShM marks showed plausible associations with SCZ or BPD based on their effects on interactions among transcription factors (TFs), DNA methylation levels, or other epigenomic marks and thus contributed to dysregulated gene expression, which ultimately increased disease susceptibility. We then validated that competitive binding of POU3F2 on the alternative allele at the AShM site rs4558409 (G/T) in PLLP-enhanced PLLP expression, while the hydroxymethylated alternative allele, which alleviated the POU3F2 binding activity at the rs4558409 site, might be associated with the downregulated PLLP expression observed in BPD or SCZ. Moreover, disruption of rs4558409 promoted neural development and vesicle trafficking.

CONCLUSION

Our study provides a powerful strategy for prioritizing regulatory risk variants and contributes to our understanding of the interplay between genetic and epigenetic factors in mediating SCZ or BPD susceptibility.

Genome-Environment Interactions and Psychiatric Disorders

Environmental factors are known to interact with the genome by altering epigenetic mechanisms regulating gene expression and contributing to the pathogenesis of psychiatric disorders. This article is a narrative review of how the major environmental factors contribute to the pathogenesis of common psychiatric disorders such as schizophrenia, bipolar disorder, major depressive disorder, and anxiety disorder this way. The cited articles were published between 1 January 2000 and 31 December 2022 and were obtained from PubMed and Google Scholar. The search terms used were as follows: gene or genetic; genome; environment; mental or psychiatric disorder; epigenetic; and interaction. The following environmental factors were found to act epigenetically on the genome to influence the pathogenesis of psychiatric disorders: social determinants of mental health, maternal prenatal psychological stress, poverty, migration, urban dwelling, pregnancy and birth complications, alcohol and substance abuse, microbiota, and prenatal and postnatal infections. The article also discusses the ways by which factors such as drugs, psychotherapy, electroconvulsive therapy, and physical exercise act epigenetically to alleviate the symptoms of psychiatric disorders in affected patients. These data will be useful information for clinical psychiatrists and those researching the pathogenesis and treatment of psychiatric disorders.

Long-read sequencing of diagnosis and post-therapy medulloblastoma reveals complex rearrangement patterns and epigenetic signatures

Cancer genomes harbor a broad spectrum of structural variants (SVs) driving tumorigenesis, a relevant subset of which escape discovery using short-read sequencing. We employed Oxford Nanopore Technologies (ONT) long-read sequencing in a paired diagnostic and post-therapy medulloblastoma to unravel the haplotype-resolved somatic genetic and epigenetic landscape. We assembled complex rearrangements, including a 1.55-Mbp chromothripsis event, and we uncover a complex SV pattern termed templated insertion (TI) thread, characterized by short (mostly <1 kb) insertions showing prevalent self-concatenation into highly amplified structures of up to 50 kbp in size. TI threads occur in 3% of cancers, with a prevalence up to 74% in liposarcoma, and frequent colocalization with chromothripsis. We also perform long-read-based methylome profiling and discover allele-specific methylation (ASM) effects, complex rearrangements exhibiting differential methylation, and differential promoter methylation in cancer-driver genes. Our study shows the advantage of long-read sequencing in the discovery and characterization of complex somatic rearrangements.

The Role of Cluster C19MC in Pre-Eclampsia Development

Pre-eclampsia is a placenta-related complication occurring in 2-10% of all pregnancies. miRNAs are a group of non-coding RNAs regulating gene expression. There is evidence that C19MC miRNAs are involved in the development of the placenta. Deregulation of chromosome 19 microRNA cluster (C19MC) miRNAs expression leads to impaired cell differentiation, abnormal trophoblast invasion and pathological angiogenesis, which can lead to the development of pre-eclampsia. Information was obtained through a review of articles available in PubMed Medline. Articles on the role of the C19MC miRNA in the development of pre-eclampsia published in 2009-2022 were analyzed. This review article summarizes the current data on the role of the C19MC miRNA in the development of pre-eclampsia. They indicate a significant increase in the expression of most C19MC miRNAs in placental tissue and a high level of circulating fractions in serum and plasma, both in the first and/or third trimester in women with PE. Only for miR-525-5p, low levels of plasma expression were noted in the first trimester, and in the placenta in the third trimester. The search for molecular factors indicating the development of pre-eclampsia before the onset of clinical symptoms seems to be a promising diagnostic route. Identifying women at risk of developing pre-eclampsia at the pre-symptomatic stage would avoid serious complications in both mothers and fetuses. We believe that miRNAs belonging to cluster C19MC could be promising biomarkers of pre-eclampsia development.

Identification of genetic variants influencing methylation in brain with pleiotropic effects on psychiatric disorders

Psychiatric disorders affect 29% of the global population at least once in the lifespan, and genetic studies have proved a shared genetic basis among them, although the underlying molecular mechanisms remain largely unknown. DNA methylation plays an important role in complex disorders and, remarkably, enrichment of common genetic variants influencing allele-specific methylation (ASM) has been reported among variants associated with specific psychiatric disorders. In the present study we assessed the contribution of ASM to a set of eight psychiatric disorders by combining genetic, epigenetic and expression data. We interrogated a list of 3896 ASM tagSNPs in the brain in the summary statistics of a cross-disorder GWAS meta-analysis of eight psychiatric disorders from the Psychiatric Genomics Consortium, including more than 162,000 cases and 276,000 controls. We identified 80 SNPs with pleiotropic effects on psychiatric disorders that show an opposite directional effect on methylation and gene expression. These SNPs converge on eight candidate genes: ZSCAN29, ZSCAN31, BTN3A2, DDAH2, HAPLN4, ARTN, FAM109B and NAGA. ZSCAN29 shows the broadest pleiotropic effects, showing associations with five out of eight psychiatric disorders considered, followed by ZSCAN31 and BTN3A2, associated with three disorders. All these genes overlap with CNVs related to cognitive phenotypes and psychiatric traits, they are expressed in the brain, and seven of them have previously been associated with specific psychiatric disorders, supporting our results. To sum up, our integrative functional genomics analysis identified eight psychiatric disease risk genes that impact a broad list of disorders and highlight an etiologic role of SNPs that influence DNA methylation and gene expression in the brain.

Early adversity as the prototype gene × environment interaction in mental disorders?

Childhood adversity (CA) as a significant stressor has consistently been associated with the development of mental disorders. The interaction between CA and genetic variants has been proposed to play a substantial role in disease etiology. In this review, we focus on the gene by environment (GxE) paradigm, its background and interpretation and stress the necessity of its implementation in psychiatric research. Further, we discuss the findings supporting GxCA interactions, ranging from candidate gene studies to polygenic and genome-wide approaches, their strengths and limitations. To illustrate potential underlying epigenetic mechanisms by which GxE effects are translated, we focus on results from FKBP5 × CA studies and discuss how molecular evidence can supplement previous GxE findings. In conclusion, while GxE studies constitute a valuable line of investigation, more harmonized GxE studies in large, deep-phenotyped, longitudinal cohorts, and across different developmental stages are necessary to further substantiate and understand reported GxE findings.

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.

Genome-wide sequencing-based identification of methylation quantitative trait loci and their role in schizophrenia risk

DNA methylation (DNAm) is an epigenetic regulator of gene expression and a hallmark of gene-environment interaction. Using whole-genome bisulfite sequencing, we have surveyed DNAm in 344 samples of human postmortem brain tissue from neurotypical subjects and individuals with schizophrenia. We identify genetic influence on local methylation levels throughout the genome, both at CpG sites and CpH sites, with 86% of SNPs and 55% of CpGs being part of methylation quantitative trait loci (meQTLs). These associations can further be clustered into regions that are differentially methylated by a given SNP, highlighting the genes and regions with which these loci are epigenetically associated. These findings can be used to better characterize schizophrenia GWAS-identified variants as epigenetic risk variants. Regions differentially methylated by schizophrenia risk-SNPs explain much of the heritability associated with risk loci, despite covering only a fraction of the genomic space. We provide a comprehensive, single base resolution view of association between genetic variation and genomic methylation, and implicate schizophrenia GWAS-associated variants as influencing the epigenetic plasticity of the brain.

The authors provide a comprehensive, single base resolution view of association between genetic variation and DNA methylation in human brain. They also show that heritability attributed to schizophrenia GWAS-associated variants reflects the epigenetic plasticity of the brain.

Diversification of the CRISPR Toolbox: Applications of CRISPR-Cas Systems Beyond Genome Editing

The discovery of CRISPR has revolutionized the field of genome engineering, but the potential of this technology is far from reaching its limits. In this review, we explore the broad range of applications of CRISPR technology to highlight the rapid expansion of the field beyond gene editing alone. It has been demonstrated that CRISPR technology can control gene expression, spatiotemporally image the genome in vivo, and detect specific nucleic acid sequences for diagnostics. In addition, new technologies are under development to improve CRISPR quality controls for gene editing, thereby improving the reliability of these technologies for therapeutics and beyond. These are just some of the many CRISPR tools that have been developed in recent years, and the toolbox continues to diversify.

A forensic case study for body fluid identification using DNA methylation analysis

Recently, a method of identifying body fluids using DNA methylation has been developed (Frumkin et al., 2011). An existing multiplex assay using 9 CpG markers could differentiate 5 body fluids: semen, blood, saliva, menstrual blood, and vaginal fluid. To validate this technique, we evaluated the previously described body fluid identification method by means of single base extension (SBE). DNA methylation was applied to 22 samples in 18 forensic cases; seven of these were semen, three were blood, eight were saliva, three were vaginal fluid, and one was menstrual blood. Total of 18 samples were tested, the DNA methylation profiles were coincident from preliminary tests (acid phosphatase (AP), leucomalachite green (LMG, Sigma Aldrich, St Louis, MO, USA) and SALIgAE®) except one sample which displayed an all-negative result. After applying the DNA methylation method to forensic samples, we determined that it could be very useful for differentiating vaginal secretions from menstrual blood, for which there is no conventional preliminary testing method.

Exploring allele specific methylation in drug dependence susceptibility

Drug dependence is a neuropsychiatric condition that involves genetic, epigenetic and environmental factors. Allele-specific methylation (ASM) is a common and stable epigenetic mechanism that involves genetic variants correlating with differential levels of methylation at CpG sites. We selected 182 single-nucleotide polymorphisms (SNPs) described to influence cis ASM in human brain regions to evaluate their possible contribution to drug dependence susceptibility. We performed a case-control association study in a discovery sample of 578 drug-dependent patients (including 428 cocaine-dependent subjects) and 656 controls from Spain, and then, we followed-up the significant associations in an independent sample of 1119 cases (including 589 cocaine-dependent subjects) and 1092 controls. In the discovery sample, we identified five nominal associations, one of them replicated in the follow-up sample (rs6020251). The pooled analysis revealed an association between drug dependence and rs6020251 but also rs11585570, both overcoming the Bonferroni correction for multiple testing. We performed the same analysis considering only cocaine-dependent patients and obtained similar results. The rs6020251 variant correlates with differential methylation levels of cg17974185 and lies in the first intron of the CTNNBL1 gene, in a genomic region with multiple histone marks related to enhancer and promoter regions in brain. Rs11585570 is an eQTL in brain and blood for the SCP2 and ECHDC2 genes and correlates with differential methylation of cg27535305 and cg13461509, located in the promoter regions of both genes. To conclude, using an approach that combines genetic and epigenetic data, we highlighted the CTNNBL1, SCP2 and ECHDC2 genes as potential contributors to drug dependence susceptibility.

CTCF-Mediated Genome Architecture Regulates the Dosage of Mitotically Stable Mono-allelic Expression of Autosomal Genes

The mechanisms that guide the clonally stable random mono-allelic expression of autosomal genes remain enigmatic. We show that (1) mono-allelically expressed (MAE) genes are assorted and insulated from bi-allelically expressed (BAE) genes through CTCF-mediated chromatin loops; (2) the cell-type-specific dynamics of mono-allelic expression coincides with the gain and loss of chromatin insulator sites; (3) dosage of MAE genes is more sensitive to the loss of chromatin insulation than that of BAE genes; and (4) inactive alleles of MAE genes are significantly more insulated than active alleles and are de-repressed upon CTCF depletion. This alludes to a topology wherein the inactive alleles of MAE genes are insulated from the spatial interference of transcriptional states from the neighboring bi-allelic domains via CTCF-mediated loops. We propose that CTCF functions as a typical insulator on inactive alleles, but facilitates transcription through enhancer-linking on active allele of MAE genes, indicating widespread allele-specific regulatory roles of CTCF.

pWGBSSimla: a profile-based whole-genome bisulfite sequencing data simulator incorporating methylation QTLs, allele-specific methylations and differentially methylated regions

MOTIVATION

DNA methylation plays an important role in regulating gene expression. DNA methylation is commonly analyzed using bisulfite sequencing (BS-seq)-based designs, such as whole-genome bisulfite sequencing (WGBS), reduced representation bisulfite sequencing (RRBS) and oxidative bisulfite sequencing (oxBS-seq). Furthermore, there has been growing interest in investigating the roles that genetic variants play in changing the methylation levels (i.e. methylation quantitative trait loci or meQTLs), how methylation regulates the imprinting of gene expression (i.e. allele-specific methylation or ASM) and the differentially methylated regions (DMRs) among different cell types. However, none of the current simulation tools can generate different BS-seq data types (e.g. WGBS, RRBS and oxBS-seq) while modeling meQTLs, ASM and DMRs.

RESULTS

We developed profile-based whole-genome bisulfite sequencing data simulator (pWGBSSimla), a profile-based bisulfite sequencing data simulator, which simulates WGBS, RRBS and oxBS-seq data for different cell types based on real data. meQTLs and ASM are modeled based on the block structures of the methylation status at CpGs, whereas the simulation of DMRs is based on observations of methylation rates in real data. We demonstrated that pWGBSSimla adequately simulates data and allows performance comparisons among different methylation analysis methods.

AVAILABILITY AND IMPLEMENTATION

pWGBSSimla is available at https://omicssimla.sourceforge.io.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Bumblebee Workers Show Differences in Allele-Specific DNA Methylation and Allele-Specific Expression

Allele-specific expression is when one allele of a gene shows higher levels of expression compared with the other allele, in a diploid organism. Recent work has identified allele-specific expression in a number of Hymenopteran species. However, the molecular mechanism which drives this allelic expression bias remains unknown. In mammals, DNA methylation is often associated with genes which show allele-specific expression. DNA methylation systems have been described in species of Hymenoptera, providing a candidate mechanism. Using previously generated RNA-Seq and whole-genome bisulfite sequencing from reproductive and sterile bumblebee (Bombus terrestris) workers, we have identified genome-wide allele-specific expression and allele-specific DNA methylation. The majority of genes displaying allele-specific expression are common between reproductive and sterile workers and the proportion of allele-specific expression bias generally varies between genetically distinct colonies. We have also identified genome-wide allele-specific DNA methylation patterns in both reproductive and sterile workers, with reproductive workers showing significantly more genes with allele-specific methylation. Finally, there is no significant overlap between genes showing allele-specific expression and allele-specific methylation. These results indicate that cis-acting DNA methylation does not directly drive genome-wide allele-specific expression in this species.

SMADs binding site polymorphisms rs9911630 is associated with susceptibility but not prognosis of gastric cancer: a case control study

Background: Single nucleotide polymorphisms (SNPs) in transcription factor binding sites (TFBS) can change their binding strength, affecting the function of transcription factors (TFs). Small mother against decapentaplegic (SMAD) proteins are known as a family of TFs involved in tumorigenesis. We performed this study to investigate whether SNPs in SMADs binding sites affect the susceptibility or prognosis of gastric cancer (GC).

Methods: Using bioinformatics tools, we focused on the association between rs9911630 polymorphism and GC. We performed this case-control study in 1275 GC patients and 1426 cancer-free subjects using TaqMan allelic discrimination method.

Results: We found that rs9911630 A>G polymorphism was associate to an increased risk of gastric cancer (adjusted OR for additive model = 1.16; 95% CI = 1.03-1.30). Furthermore, we assess whether rs9911630 polymorphism affected the prognosis of GC. However, no significant association was discovered between rs9911630 A>G polymorphism and overall survival time of GC patients (HR for addictive model = 1.01; 95%CI = 0.88-1.15).

Conclusions: Our results suggested that rs9911630 polymorphism in SMADs target site might influence susceptibility but not prognosis of gastric cancer.

DNA methylation near the INS gene is associated with INS genetic variation (rs689) and type 1 diabetes in the Diabetes Autoimmunity Study in the Young

OBJECTIVE

Mechanisms underlying the role of non-human leukocyte antigen (HLA) genetic risk variants in type 1 diabetes (T1D) are poorly understood. We aimed to test the association between methylation and non-HLA genetic risk.

METHODS

We conducted a methylation quantitative trait loci (mQTL) analysis in a nested case-control study from the Dietary Autoimmunity Study in the Young. Controls (n = 83) were frequency-matched to T1D cases (n = 83) based on age, race/ethnicity, and sample availability. We evaluated 13 non-HLA genetic markers known be associated with T1D. Genome-wide methylation profiling was performed on peripheral blood samples collected prior to T1D using the Illumina 450 K (discovery set) and infinium methylation EPIC beadchip (EPIC validation) platforms. Linear regression models, adjusting for age and sex, were used to test to each single nucleotide polymorphism (SNP) -probe combination. Logistic regression models were used to test the association between T1D and methylation levels among probes with a significant mQTL. A meta-analysis was used to combine odds ratios from the two platforms.

RESULTS

We identified 10 SNP-methylation probe pairs (false discovery rate (FDR) adjusted P < .05 and validation P < .05). Probes were associated with the GSDMB, C1QTNF6, IL27, and INS genes. The cg03366382 (OR: 1.9, meta-P = .0495), cg21574853 (OR: 2.5, meta-P = .0232), and cg25336198 (odds ratio: 6.6, meta-P = .0081) probes were significantly associated with T1D. The three probes were located upstream from the INS transcription start site.

CONCLUSIONS

We confirmed an association between DNA methylation and rs689 that has been identified in related studies. Measurements in our study preceded the onset of T1D suggesting methylation may have a role in the relationship between INS variation and T1D development.

The applied implications of epigenetics in anxiety, affective and stress-related disorders - A review and synthesis on psychosocial stress, psychotherapy and prevention

Mental disorders are highly complex and multifactorial in origin, comprising an elaborate interplay of genetic and environmental factors. Epigenetic mechanisms such as DNA modifications (e.g. CpG methylation), histone modifications (e.g. acetylation) and microRNAs function as a translator between genes and the environment. Indeed, environmental influences such as exposure to stress shape epigenetic patterns, and lifetime experiences continue to alter the function of the genome throughout the lifespan. Here, we summarize the recently burgeoning body of research regarding the involvement of aberrant epigenetic signatures in mediating an increased vulnerability to a wide range of mental disorders. We review the current knowledge of epigenetic changes to constitute useful markers predicting the clinical response to psychotherapeutic interventions, and of psychotherapy to alter - and potentially reverse - epigenetic risk patterns. Given first evidence pointing to a transgenerational transmission of epigenetic information, epigenetic alterations arising from successful psychotherapy might be transferred to future generations and thus contribute to the prevention of mental disorders. Findings are integrated into a multi-level framework highlighting challenges pertaining to the mechanisms of action and clinical implications of epigenetic research. Promising future directions regarding the prediction, prevention, and personalized treatment of mental disorders in line with a 'precision medicine' approach are discussed.

Epigenetics Underlying Susceptibility and Resilience Relating to Daily Life Stress, Work Stress, and Socioeconomic Status

Susceptibility and resilience to mental disorders result from a complex choreography of gene-environment interactions with epigenetics at the intersection of external psychological stressors and internal biological systems. Increasing awareness of the growing disease burden influenced by daily life stress ("daily hassles"), work-related stress, and low socioeconomic status (SES) has resulted in a novel interest into their underlying molecular signatures. This review offers a brief outline of psychiatric epigenetics and a comprehensive overview of recent findings exploring the relationship of various occupational stressors and DNA methylation in epigenome-wide association studies (EWAS) and in candidate gene studies including the serotonin transporter (SLC6A4; 5-HTTLPR), melatonin receptor 1A (MTNR1A), brain-derived neurotrophic factor (BDNF), tyrosine hydroxylase (TH), and the protein family of DNA methyltransferases (DNMTs). Conceptual and methodological challenges of epigenetic investigations with a special focus on gene-environment interactions are highlighted and discussed. The findings are integrated into a pathophysiological framework featuring epigenetic plasticity factors and work-related stress as a possible central detrimental component targetable by workplace interventions. Finally, the potential of dynamic epigenetic biomarkers of treatment response to pharmacotherapy or psychotherapy is expanded upon.

Exploring genetic variation that influences brain methylation in attention-deficit/hyperactivity disorder

Attention-deficit/hyperactivity disorder (ADHD) is a neurodevelopmental disorder caused by an interplay of genetic and environmental factors. Epigenetics is crucial to lasting changes in gene expression in the brain. Recent studies suggest a role for DNA methylation in ADHD. We explored the contribution to ADHD of allele-specific methylation (ASM), an epigenetic mechanism that involves SNPs correlating with differential levels of DNA methylation at CpG sites. We selected 3896 tagSNPs reported to influence methylation in human brain regions and performed a case-control association study using the summary statistics from the largest GWAS meta-analysis of ADHD, comprising 20,183 cases and 35,191 controls. We observed that genetic risk variants for ADHD are enriched in ASM SNPs and identified associations with eight tagSNPs that were significant at a 5% false discovery rate (FDR). These SNPs correlated with methylation of CpG sites lying in the promoter regions of six genes. Since methylation may affect gene expression, we inspected these ASM SNPs together with 52 ASM SNPs in high LD with them for eQTLs in brain tissues and observed that the expression of three of those genes was affected by them. ADHD risk alleles correlated with increased expression (and decreased methylation) of ARTN and PIDD1 and with a decreased expression (and increased methylation) of C2orf82. Furthermore, these three genes were predicted to have altered expression in ADHD, and genetic variants in C2orf82 correlated with brain volumes. In summary, we followed a systematic approach to identify risk variants for ADHD that correlated with differential cis-methylation, identifying three novel genes contributing to the disorder.

Methylation of the C19MC microRNA locus in the placenta: association with maternal and chilhood body size

OBJECTIVES

To study DNA methylation at the C19MC locus in the placenta and its association with (1) parental body size, (2) transmission of haplotypes for the C19MC rs55765443 SNP, and (3) offspring's body size and/or body composition at birth and in childhood.

SUBJECTS AND METHODS

Seventy-two pregnant women-infant pairs and 63 fathers were included in the study. Weight and height of mothers, fathers and newborns were registered during pregnancy or at birth (n = 72). Placental DNA methylation at the C19MC imprinting control region (ICR) was quantified by bisulfite pyrosequencing. Genotyping of the SNP was performed using restriction fragment length polymorphisms. The children's body size and composition were reassessed at age 6 years (n = 32).

RESULTS

Lower levels of placental C19MC methylation were associated with increased body size of mother, specifically with higher pregestational and predelivery weights and height of the mother (β from -0.294 to -0.371; R² from 0.04 to 0.10 and all p < 0.019), and with higher weight, height, waist and hip circumferences, and fat mass of the child (β from -0.428 to -0.552; R² from 0.33 to 0.56 and all p < 0.009). Parental transmission of the SNP did not correlate with an altered placental methylation status at the C19MC ICR.

CONCLUSIONS

Increased maternal size is associated with reduced placental C19MC methylation, which, in turn, relate to larger body size of the child.

A method for identifying allele-specific hydroxymethylation

We previously identified sequence-dependent allele-specific methylation (sd-ASM) in adult human peripheral blood leukocytes, in which ASM occurs in cis depending on adjacent polymorphic sequences. A number of groups have identified sd-ASM sites in the human and mouse genomes, illustrating the prevalence of sd-ASM in mammalian genomes. In addition, sd-ASM can lead to sequence-dependent allele-specific expression of neighbouring genes. Imprinted genes also often exhibit parent-of-origin-dependent allele-specific methylation (pd-ASM), which causes parent-of-origin-dependent allele-specific expression. However, whether most of the already known sd-ASM and pd-ASM sites are methylated or hydroxymethylated remains unclear due to technical restrictions. Accordingly, a novel method that enables examination of allelic methylation and hydroxymethylation status and also overcomes the drawbacks of conventional methods is needed. Such a method could also be used to elucidate the mechanisms underlying polymorphism-associated inter-individual differences in disease susceptibility and the mechanism of genomic imprinting. Here, we developed a simple method to determine allelic hydroxymethylation status and identified novel sequence- and parent-of-origin-dependent allele-specific hydroxymethylation sites. Correlation analyses of TF binding sequences and methylation or hydroxymethylation between three mouse strains revealed the involvement of Pax5 in strain-specific methylation and hydroxymethylation in exon 7 of Pdgfrb.

Interrogating Mutant Allele Expression via Customized Reference Genomes to Define Influential Cancer Mutations

Genetic alterations are essential for cancer initiation and progression. However, differentiating mutations that drive the tumor phenotype from mutations that do not affect tumor fitness remains a fundamental challenge in cancer biology. To better understand the impact of a given mutation within cancer, RNA-sequencing data was used to categorize mutations based on their allelic expression. For this purpose, we developed the MAXX (Mutation Allelic Expression Extractor) software, which is highly effective at delineating the allelic expression of both single nucleotide variants and small insertions and deletions. Results from MAXX demonstrated that mutations can be separated into three groups based on their expression of the mutant allele, lack of expression from both alleles, or expression of only the wild-type allele. By taking into consideration the allelic expression patterns of genes that are mutated in PDAC, it was possible to increase the sensitivity of widely used driver mutation detection methods, as well as identify subtypes that have prognostic significance and are associated with sensitivity to select classes of therapeutic agents in cell culture. Thus, differentiating mutations based on their mutant allele expression via MAXX represents a means to parse somatic variants in tumor genomes, helping to elucidate a gene’s respective role in cancer.

Unique DNA Methylation Profiles Are Associated with cis-Variation in Honey Bees

DNA methylation is an important epigenetic modification that mediates diverse processes such as cellular differentiation, phenotypic plasticity, and genomic imprinting. Mounting evidence suggests that local DNA sequence variation can be associated with particular DNA methylation states, indicating that the interplay between genetic and epigenetic factors may contribute synergistically to the phenotypic complexity of organisms. Social insects such as ants, bees, and wasps have extensive phenotypic plasticity manifested in their different castes, and this plasticity has been associated with variation in DNA methylation. Yet, the influence of genetic variation on DNA methylation state remains mostly unknown. Here we examine the importance of sequence-specific methylation at the genome-wide level, using whole-genome bisulfite sequencing of the semen of individual honey bee males. We find that individual males harbor unique DNA methylation patterns in their semen, and that genes that are more variable at the epigenetic level are also more likely to be variable at the genetic level. DNA sequence variation can affect DNA methylation by modifying CG sites directly, but can also be associated with local variation in cis that is not CG-site specific. We show that covariation in sequence polymorphism and DNA methylation state contributes to the individual-specificity of epigenetic marks in social insects, which likely promotes their retention across generations, and their capacity to influence evolutionary adaptation.

Who's your daddy? Behavioral and epigenetic consequences of paternal drug exposure

Substance use disorders (SUDs) reflect genetic and environmental factors. While identifying reliable genetic variants that predispose individuals to developing SUDs has been challenging, epigenetic factors may also contribute to the heritability of SUDs. Familial drug use associates with a wide range of problems in children, including an increased risk for developing a SUD. The implications of maternal drug use on offspring development are a well-studied area; however, paternal drug use prior to conception has received relatively little attention. Paternal exposure to several environmental stimuli (i.e. stress or diet manipulations) results in behavioral and epigenetic changes in offspring. The purpose of this review is to determine the state of the preclinical literature on the behavioral and epigenetic consequences of paternal drug exposure. Drug-sired offspring show several developmental and physiological abnormalities. These offspring also show deficits in cognitive and emotional domains. Examining sensitivity to drugs in offspring is a growing area of research. Drug-sired offspring are resistant to the rewarding and reinforcing properties of drugs. However, greater paternal motivation for the drug, combined with high drug intake, can result in addiction-like behaviors in offspring. Drug-sired offspring also show altered histone modifications and DNA methylation levels of imprinted genes and microRNAs; epigenetic-mediated changes were also noted in genes related to glutamatergic and neurotrophic factor signaling. In some instances, drug use resulted in aberrant epigenetic modifications in sire sperm, and these changes were maintained in the brains of offspring. Thus, paternal drug exposure has long-lasting consequences that include altered drug sensitivity in subsequent generations. We discuss factors (i.e. maternal behaviors) that may moderate these paternal drug-induced effects as well as ideas for future directions.

Examining the Impact of Imputation Errors on Fine-Mapping Using DNA Methylation QTL as a Model Trait

This study highlights dangers in over-interpreting fine-mapping results. Chundru et al. show that genotype imputation accuracy has a large impact on fine-mapping accuracy. They used DNA methylation at CpG-sites with a variant...

Genetic variants disrupting DNA methylation at CpG dinucleotides (CpG-SNP) provide a set of known causal variants to serve as models to test fine-mapping methodology. We use 1716 CpG-SNPs to test three fine-mapping approaches (Bayesian imputation-based association mapping, Bayesian sparse linear mixed model, and the J-test), assessing the impact of imputation errors and the choice of reference panel by using both whole-genome sequence (WGS), and genotype array data on the same individuals (n = 1166). The choice of imputation reference panel had a strong effect on imputation accuracy, with the 1000 Genomes Project Phase 3 (1000G) reference panel (n = 2504 from 26 populations) giving a mean nonreference discordance rate between imputed and sequenced genotypes of 3.2% compared to 1.6% when using the Haplotype Reference Consortium (HRC) reference panel (n = 32,470 Europeans). These imputation errors had an impact on whether the CpG-SNP was included in the 95% credible set, with a difference of ∼23% and ∼7% between the WGS and the 1000G and HRC imputed datasets, respectively. All of the fine-mapping methods failed to reach the expected 95% coverage of the CpG-SNP. This is attributed to secondary cis genetic effects that are unable to be statistically separated from the CpG-SNP, and through a masking mechanism where the effect of the methylation disrupting allele at the CpG-SNP is hidden by the effect of a nearby SNP that has strong linkage disequilibrium with the CpG-SNP. The reduced accuracy in fine-mapping a known causal variant in a low-level biological trait with imputed genetic data has implications for the study of higher-order complex traits and disease.

The interaction between oxytocin receptor gene methylation and maternal behavior on children's early theory of mind abilities

Theory of mind, the ability to represent the mental states of others, is an important social cognitive process, which contributes to the development of social competence. Recent research suggests that interactions between gene and environmental factors, such as oxytocin receptor gene (OXTR) polymorphisms and maternal parenting behavior, may underlie individual differences in children's theory of mind. However, the potential influence of DNA methylation of OXTR remains unclear. The current study investigated the roles of OXTR methylation, maternal behavior, and their statistical interaction on toddlers' early emerging theory of mind abilities. Participants included a community sample of 189 dyads of mothers and their 2- to 3-year-old children, whose salivary DNA was analyzed. Results indicated that more maternal structuring behavior was associated with better performance, on a battery of three theory of mind tasks, while higher OXTR methylation within exon 3 was associated with poorer performance. A significant interaction also emerged, such that OXTR methylation was related to theory of mind among children whose mothers displayed less structuring, when controlling for children's age, sex, ethnicity, number of child-aged siblings, verbal ability, and maternal education. Maternal structuring behavior may buffer the potential negative impact of hypermethylation on OXTR gene expression and function.

Early-life adversity-induced long-term epigenetic programming associated with early onset of chronic physical aggression: Studies in humans and animals

Objectives: To examine whether chronic physical aggression (CPA) in adulthood can be epigenetically programmed early in life due to exposure to early-life adversity. Methods: Literature search of public databases such as PubMed/MEDLINE and Scopus. Results: Children/adolescents susceptible for CPA and exposed to early-life abuse fail to efficiently cope with stress that in turn results in the development of CPA later in life. This phenomenon was observed in humans and animal models of aggression. The susceptibility to aggression is a complex trait that is regulated by the interaction between environmental and genetic factors. Epigenetic mechanisms mediate this interaction. Subjects exposed to stress early in life exhibited long-term epigenetic programming that can influence their behaviour in adulthood. This programming affects expression of many genes not only in the brain but also in other systems such as neuroendocrine and immune. Conclusions: The propensity to adult CPA behaviour in subjects experienced to early-life adversity is mediated by epigenetic programming that involves long-term systemic epigenetic alterations in a whole genome.

Relationship between Alzheimer's disease-associated SNPs within the CLU gene, local DNA methylation and episodic verbal memory in healthy and schizophrenia subjects

Genetic variation may impact on local DNA methylation patterns. Therefore, information about allele-specific DNA methylation (ASM) within disease-related loci has been proposed to be useful for the interpretation of GWAS results. To explore mechanisms that may underlie associations between Alzheimer's disease (AD) and schizophrenia risk CLU gene and verbal memory, one of the most affected cognitive domains in both conditions, we studied DNA methylation in a region between AD-associated SNPs rs9331888 and rs9331896 in 72 healthy individuals and 73 schizophrenia patients. Using single-molecule real-time bisulfite sequencing we assessed the haplotype-dependent ASM in this region. We then investigated whether its methylation could influence episodic verbal memory measured with the Rey Auditory Verbal Learning Test in these two cohorts. The region showed a complex methylation pattern, which was similar in healthy and schizophrenia individuals and unrelated to haplotypes. The pattern predicted memory scores in controls. The results suggest that epigenetic modifications within the CLU locus may play a role in memory variation, independent of ASM.

Allele-specific methylation of imprinted genes in fetal cord blood is influenced by cis-acting genetic variants and parental factors

Aim:

To examine the effects of genetic variation, parental age and BMI on parental allele-specific methylation of imprinted genes in fetal cord blood samples.

Methodology:

We have developed SNP genotyping and deep bisulphite sequencing assays for six imprinted genes to determine parental allele-specific methylation patterns in diploid somatic tissues.

Results:

Multivariate linear regression analyses revealed a negative correlation of paternal age with paternal MEG3 allele methylation in fetal cord blood. Methylation of the maternal PEG3 allele showed a positive correlation with maternal age. Paternal BMI was positively correlated with paternal MEST allele methylation. In addition to parental origin, allele-specific methylation of most imprinted genes was largely dependent on the underlying SNP haplotype.

Conclusion:

Our study supports the idea that parental factors can have an impact, although of small effect size, on the epigenome of the next generation, providing an additional layer of complexity to phenotypic diversity.

Methylomic profiling and replication implicates deregulation of PCSK9 in alcohol use disorder

Alcohol use disorder (AUD) is a common and chronic disorder with substantial effects on personal and public health. The underlying pathophysiology is poorly understood but strong evidence suggests significant roles of both genetic and epigenetic components. Given that alcohol affects many organ systems, we performed a cross-tissue and cross-phenotypic analysis of genome-wide methylomic variation in AUD using samples from 3 discovery, 4 replication, and 2 translational cohorts. We identified a differentially methylated region in the promoter of the proprotein convertase subtilisin/kexin 9 (PCSK9) gene that was associated with disease phenotypes. Biological validation showed that PCSK9 promoter methylation is conserved across tissues and positively correlated with expression. Replication in AUD datasets confirmed PCSK9 hypomethylation and a translational mouse model of AUD showed that alcohol exposure leads to PCSK9 downregulation. PCSK9 is primarily expressed in the liver and regulates low-density lipoprotein cholesterol (LDL-C). Our finding of alcohol-induced epigenetic regulation of PCSK9 represents one of the underlying mechanisms between the well-known effects of alcohol on lipid metabolism and cardiovascular risk, with light alcohol use generally being protective while chronic heavy use has detrimental health outcomes.

FADS1-FADS2 genetic polymorphisms are associated with fatty acid metabolism through changes in DNA methylation and gene expression

Background

Genome-wide association studies (GWASs) have shown that genetic variants are important determinants of free fatty acid levels. The mechanisms underlying the associations between genetic variants and free fatty acid levels are incompletely understood. Here, we aimed to identify genetic markers that could influence diverse fatty acid levels in a Chinese population and uncover the molecular mechanisms in terms of DNA methylation and gene expression.

Results

We identified strong associations between single-nucleotide polymorphisms (SNPs) in the fatty acid desaturase (FADS) region and multiple polyunsaturated fatty acids. Expression quantitative trait locus (eQTL) analysis of rs174570 on FADS1 and FADS2 mRNA levels proved that minor allele of rs174570 was associated with decreased FADS1 and FADS2 expression levels (P < 0.05). Methylation quantitative trait locus (mQTL) analysis of rs174570 on DNA methylation levels in three selected regions of FADS region showed that the methylation levels at four CpG sites in FADS1, one CpG site in intragenic region, and three CpG sites in FADS2 were strongly associated with rs174570 (P < 0.05). Then, we demonstrated that methylation levels at three CpG sites in FADS1 were negatively associated with FADS1 and FADS2 expression, while two CpG sites in FADS2 were positively associated with FADS1 and FADS2 expression. Using mediation analysis, we further show that the observed effect of rs174570 on gene expression was tightly correlated with the effect predicted through association with methylation.

Conclusions

Our findings suggest that genetic variants in the FADS region are major genetic modifiers that can regulate fatty acid metabolism through epigenetic gene regulation.

Electronic supplementary material

The online version of this article (10.1186/s13148-018-0545-5) contains supplementary material, which is available to authorized users.

Genetic Variation in the H19-IGF2 Cluster Might Confer Risk of Developing Impaired Renal Function

The H19-IGF2 imprinted gene region could be implicated in the risk of developing impaired renal function (IRF). Our aim was to determine the association of several common H19-IGF2 variants and IRF in a cohort of elderly healthy individuals. The study involved 675 individuals >65 years of age, 184 with type 2 diabetes mellitus (T2DM), and 105 with IRF (estimated glomerular filtration rate [eGFR] <60). They were genotyped for two common H19 single nucleotide polymorphisms (SNPs) (rs2839698 and rs10732516), one H19-IGF2 intergenic indel (rs201858505), and one indel in the 3'UTR of the IGF2. For the H19 SNPs, we also determined the allele present in the methylated chromosome through genotyping the DNA digested with a methylation-sensitive endonuclease. None of the four H19-IGF2 variants was associated with IRF in our cohort. We found a significantly higher frequency of the 3'UTR IGF2 deletion (D) in the eGFR <60 group (p = 0.01; odds ratio = 1.16, 95% confidence interval = 1.10-2.51). This association was independent of age and T2DM, two strong predictors of IRF. In conclusion, a common indel variant in the 3'UTR of the IGF2 gene was associated with the risk of IRF. This association could be explained by the role of IGF2 in podocyte survival, through regulation of IGF2 expression by differential binding of miRNAs to the indel sequences. Functional studies should be necessary to clarify this issue.

DNA Methylation and Psychiatric Disorders

DNA methylation has been an important area of research in the study of molecular mechanism to psychiatric disorders. Recent evidence has suggested that abnormalities in global methylation, methylation of genes, and pathways could play a role in the etiology of many forms of mental illness. In this article, we review the mechanisms of DNA methylation, including the genetic and environmental factors affecting methylation changes. We report and discuss major findings regarding DNA methylation in psychiatric patients, both within the context of global methylation studies and gene-specific methylation studies. Finally, we discuss issues surrounding data quality improvement, the limitations of current methylation analysis methods, and the possibility of using DNA methylation-based treatment for psychiatric disorders in the future.

Gene-environment Interactions in Late Life: Linking Psychosocial Stress with Brain Aging

Gene-environment interactions (GxE) can have lasting consequences on brain structure and function, potentially contributing to diverse neuropsychiatric phenotypes. This has been extensively demonstrated by studies examining GxE in childhood and early adulthood, whereas much fewer studies have addressed this question in late life. The relative paucity of studies examining GxE in late life may stem from the working hypothesis that brains become less malleable to environmental inputs as life progresses. However, while some components of brain plasticity decline with increasing age, others are retained and may even become more pronounced in old ages. Moreover, the micro- and macro-structural brain changes that accrue as a result of aging-related morbidities are likely to accentuate the susceptibility of neural circuits to environmental stressors as life advances. Supporting this hypothesis, psychosocial stress can increase the risk for late-life neuropsychiatric syndromes, especially when afflicting genetically predisposed individuals. This article reviews evidence showing how gene-stress interactions can impact the aging brain and related phenotypes in late life, and it discusses the potential mechanisms underlying such GxE and their implications for the prevention and treatment of late-life neuropsychiatric syndromes.

Allele specific expression and methylation in the bumblebee, Bombus terrestris

The social hymenoptera are emerging as models for epigenetics. DNA methylation, the addition of a methyl group, is a common epigenetic marker. In mammals and flowering plants methylation affects allele specific expression. There is contradictory evidence for the role of methylation on allele specific expression in social insects. The aim of this paper is to investigate allele specific expression and monoallelic methylation in the bumblebee, Bombus terrestris. We found nineteen genes that were both monoallelically methylated and monoallelically expressed in a single bee. Fourteen of these genes express the hypermethylated allele, while the other five express the hypomethylated allele. We also searched for allele specific expression in twenty-nine published RNA-seq libraries. We found 555 loci with allele-specific expression. We discuss our results with reference to the functional role of methylation in gene expression in insects and in the as yet unquantified role of genetic cis effects in insect allele specific methylation and expression.

Genetic-epigenetic interactions in cis: a major focus in the post-GWAS era

Studies on genetic-epigenetic interactions, including the mapping of methylation quantitative trait loci (mQTLs) and haplotype-dependent allele-specific DNA methylation (hap-ASM), have become a major focus in the post-genome-wide-association-study (GWAS) era. Such maps can nominate regulatory sequence variants that underlie GWAS signals for common diseases, ranging from neuropsychiatric disorders to cancers. Conversely, mQTLs need to be filtered out when searching for non-genetic effects in epigenome-wide association studies (EWAS). Sequence variants in CCCTC-binding factor (CTCF) and transcription factor binding sites have been mechanistically linked to mQTLs and hap-ASM. Identifying these sites can point to disease-associated transcriptional pathways, with implications for targeted treatment and prevention.

Molecular mechanisms underlying noncoding risk variations in psychiatric genetic studies

Recent large-scale genetic approaches such as genome-wide association studies have allowed the identification of common genetic variations that contribute to risk architectures of psychiatric disorders. However, most of these susceptibility variants are located in noncoding genomic regions that usually span multiple genes. As a result, pinpointing the precise variant(s) and biological mechanisms accounting for the risk remains challenging. By reviewing recent progresses in genetics, functional genomics and neurobiology of psychiatric disorders, as well as gene expression analyses of brain tissues, here we propose a roadmap to characterize the roles of noncoding risk loci in the pathogenesis of psychiatric illnesses (that is, identifying the underlying molecular mechanisms explaining the genetic risk conferred by those genomic loci, and recognizing putative functional causative variants). This roadmap involves integration of transcriptomic data, epidemiological and bioinformatic methods, as well as in vitro and in vivo experimental approaches. These tools will promote the translation of genetic discoveries to physiological mechanisms, and ultimately guide the development of preventive, therapeutic and prognostic measures for psychiatric disorders.

epiG: statistical inference and profiling of DNA methylation from whole-genome bisulfite sequencing data

The study of epigenetic heterogeneity at the level of individual cells and in whole populations is the key to understanding cellular differentiation, organismal development, and the evolution of cancer. We develop a statistical method, epiG, to infer and differentiate between different epi-allelic haplotypes, annotated with CpG methylation status and DNA polymorphisms, from whole-genome bisulfite sequencing data, and nucleosome occupancy from NOMe-seq data. We demonstrate the capabilities of the method by inferring allele-specific methylation and nucleosome occupancy in cell lines, and colon and tumor samples, and by benchmarking the method against independent experimental data.

Oxytocin and Anxiety Disorders

In the present chapter, we review the literature focusing on oxytocin (OT)-centered research in anxiety spectrum conditions, comprising separation anxiety disorder, specific phobias, social anxiety disorder (SAD), panic disorder, generalized anxiety disorder, and anxiety-related endophenotypes (e.g., trust behavior, behavioral inhibition, neuroticism, and state/trait anxiety). OT receptor gene (OXTR) polymorphisms have been implicated in gene-environment interactions with attachment style and childhood maltreatment and to influence clinical outcomes, including SAD intensity and limbic responsiveness. Epigenetic OXTR DNA methylation patterns have emerged as a link between categorical, dimensional, neuroendocrinological, and neuroimaging SAD correlates, highlighting them as potential peripheral surrogates of the central oxytocinergic tone. A pathophysiological framework of OT integrating the dynamic nature of epigenetic biomarkers and the summarized genetic and peripheral evidence is proposed. Finally, we emphasize opportunities and challenges of OT as a key network node of social interaction and fear learning in social contexts. In conjunction with multi-level investigations incorporating a dimensional understanding of social affiliation and avoidance in anxiety spectrum disorders, these concepts will help to promote research for diagnostic, state, and treatment response biomarkers of the OT system, advancing towards indicated preventive interventions and personalized treatment approaches.

Heritable DNA Methylation in CD4+ Cells among Complex Families Displays Genetic and Non-Genetic Effects

DNA methylation at CpG sites is both heritable and influenced by environment, but the relative contributions of each to DNA methylation levels are unclear. We conducted a heritability analysis of CpG methylation in human CD4+ cells across 975 individuals from 163 families in the Genetics of Lipid-lowering Drugs and Diet Network (GOLDN). Based on a broad-sense heritability (H2) value threshold of 0.4, we identified 20,575 highly heritable CpGs among the 174,445 most variable autosomal CpGs (SD > 0.02). Tests for associations of heritable CpGs with genotype at 2,145,360 SNPs using 717 of 975 individuals showed that ~74% were cis-meQTLs (< 1 Mb away from the CpG), 6% of CpGs exhibited trans-meQTL associations (>1 Mb away from the CpG or located on a different chromosome), and 20% of CpGs showed no strong significant associations with genotype (based on a p-value threshold of 1e-7). Genes proximal to the genotype independent heritable CpGs were enriched for functional terms related to regulation of T cell activation. These CpGs were also among those that distinguished T cells from other blood cell lineages. Compared to genes proximal to meQTL-associated heritable CpGs, genotype independent heritable CpGs were moderately enriched in the same genomic regions that escape erasure during primordial germ cell development and could carry potential for generational transmission.

Prenatal stress exposure, oxytocin receptor gene (OXTR) methylation, and child autistic traits: The moderating role of OXTR rs53576 genotype

Findings of studies investigating OXTR SNP rs53576 (G-A) variation in social behavior have been inconsistent, possibly because DNA methylation after stress exposure was eliminated from consideration. Our goal was to examine OXTR rs53576 allele-specific sensitivity for neonatal OXTR DNA methylation in relation to (1) a prenatal maternal stress composite, and (2) child autistic traits. Prospective data from fetal life to age 6 years were collected in a total of 743 children participating in the Generation R Study. Prenatal maternal stress exposure was uniquely associated with child autistic traits but was unrelated to OXTR methylation across both OXTR rs53576 G-allele homozygous children and A-allele carriers. For child autistic traits in general and social communication problems in particular, we observed a significant OXTR rs53576 genotype by OXTR methylation interaction in the absence of main effects, suggesting that opposing effects cancelled each other out. Indeed, OXTR methylation levels were positively associated with social problems for OXTR rs53576 G-allele homozygous children but not for A-allele carriers. These results highlight the importance of incorporating epi-allelic information and support the role of OXTR methylation in child autistic traits. Autism Res 2017, 10: 430-438. © 2016 International Society for Autism Research, Wiley Periodicals, Inc.

Allele-Skewed DNA Modification in the Brain: Relevance to a Schizophrenia GWAS

Numerous recent studies have suggested that phenotypic effects of DNA sequence variants can be mediated or modulated by their epigenetic marks, such as allele-skewed DNA modification (ASM). Using Affymetrix SNP microarrays, we performed a comprehensive search of ASM effects in human post-mortem brain and sperm samples (total n = 256) from individuals with major psychosis and control individuals. Depending on the phenotypic category of the brain samples, 1.4%-7.5% of interrogated SNPs exhibited ASM effects. Next, we investigated ASM in the context of genetic studies of schizophrenia and detected that brain ASM SNPs were significantly overrepresented among sub-threshold SNPs from a schizophrenia genome-wide association study (GWAS). Brain ASM SNPs showed a much stronger enrichment in a schizophrenia GWAS than in 17 large GWASs of non-psychiatric diseases and traits, arguing that ASM effects are at least partially tissue specific. Studies of germline and control brain ASM SNPs supported a causal association between ASM and schizophrenia. Finally, significantly higher proportions of ASM SNPs than of non-ASM SNPs were detected at loci exhibiting epigenetic signatures of enhancers and promoters, and they were overrepresented within transcription factor binding regions and DNase I hypersensitive sites. All of these findings collectively indicate that ASM SNPs should be prioritized in follow-up GWASs.

X inactivation and reactivation in X-linked diseases

X chromosome inactivation (XCI) is the phenomenon by which mammals compensate for dosage of X-linked genes in females (XX) versus males (XY). XCI patterns can be random or show extreme skewing, and can modify the mode of inheritance of X-driven phenotypes, which contributes to the variability of human pathologies. Recent findings have shown reversibility of the XCI process, which has opened new avenues in the approaches used for the treatment of X-linked diseases.

Parent-of-origin effects on genome-wide DNA methylation in the Cape honey bee (Apis mellifera capensis) may be confounded by allele-specific methylation

BACKGROUND

Intersexual genomic conflict sometimes leads to unequal expression of paternal and maternal alleles in offspring, resulting in parent-of-origin effects. In honey bees reciprocal crosses can show strong parent-of-origin effects, supporting theoretical predictions that genomic imprinting occurs in this species. Mechanisms behind imprinting in honey bees are unclear but differential DNA methylation in eggs and sperm suggests that DNA methylation could be involved. Nonetheless, because DNA methylation is multifunctional, it is difficult to separate imprinting from other roles of methylation. Here we use a novel approach to investigate parent-of-origin DNA methylation in honey bees. In the subspecies Apis mellifera capensis, reproduction of females occurs either sexually by fertilization of eggs with sperm, or via thelytokous parthenogenesis, producing female embryos derived from two maternal genomes.

RESULTS

We compared genome-wide methylation patterns of sexually-produced, diploid embryos laid by a queen, with parthenogenetically-produced diploid embryos laid by her daughters. Thelytokous embryos inheriting two maternal genomes had fewer hypermethylated genes compared to fertilized embryos, supporting the prediction that fertilized embryos have increased methylation due to inheritance of a paternal genome. However, bisulfite PCR and sequencing of a differentially methylated gene, Stan (GB18207) showed strong allele-specific methylation that was maintained in both fertilized and thelytokous embryos. For this gene, methylation was associated with haplotype, not parent of origin.

CONCLUSIONS

The results of our study are consistent with predictions from the kin theory of genomic imprinting. However, our demonstration of allele-specific methylation based on sequence shows that genome-wide differential methylation studies can potentially confound imprinting and allele-specific methylation. It further suggests that methylation patterns are heritable or that specific sequence motifs are targets for methylation in some genes.