Genetics and epigenetics of human aggression

There is substantial variation between humans in aggressive behavior, with its biological etiology and molecular genetic basis mostly unknown. This review chapter offers an overview of genomic and omics studies revealing the genetic contribution to aggression and first insights into associations with epigenetic and other omics (e.g., metabolomics) profiles. We allowed for a broad phenotype definition including studies on "aggression," "aggressive behavior," or "aggression-related traits," "antisocial behavior," "conduct disorder," and "oppositional defiant disorder." Heritability estimates based on family and twin studies in children and adults of this broadly defined phenotype of aggression are around 50%, with relatively small fluctuations around this estimate. Next, we review the genome-wide association studies (GWAS) which search for associations with alleles and also allow for gene-based tests and epigenome-wide association studies (EWAS) which seek to identify associations with differently methylated regions across the genome. Both GWAS and EWAS allow for construction of Polygenic and DNA methylation scores at an individual level. Currently, these predict a small percentage of variance in aggression. We expect that increases in sample size will lead to additional discoveries in GWAS and EWAS, and that multiomics approaches will lead to a more comprehensive understanding of the molecular underpinnings of aggression.

Epigenome-wide association study of incident type 2 diabetes: a meta-analysis of five prospective European cohorts

AIMS/HYPOTHESIS

Type 2 diabetes is a complex metabolic disease with increasing prevalence worldwide. Improving the prediction of incident type 2 diabetes using epigenetic markers could help tailor prevention efforts to those at the highest risk. The aim of this study was to identify predictive methylation markers for incident type 2 diabetes by combining epigenome-wide association study (EWAS) results from five prospective European cohorts.

METHODS

We conducted a meta-analysis of EWASs in blood collected 7-10 years prior to type 2 diabetes diagnosis. DNA methylation was measured with Illumina Infinium Methylation arrays. A total of 1250 cases and 1950 controls from five longitudinal cohorts were included: Doetinchem, ESTHER, KORA1, KORA2 and EPIC-Norfolk. Associations between DNA methylation and incident type 2 diabetes were examined using robust linear regression with adjustment for potential confounders. Inverse-variance fixed-effects meta-analysis of cohort-level individual CpG EWAS estimates was performed using METAL. The methylGSA R package was used for gene set enrichment analysis. Confirmation of genome-wide significant CpG sites was performed in a cohort of Indian Asians (LOLIPOP, UK).

RESULTS

The meta-analysis identified 76 CpG sites that were differentially methylated in individuals with incident type 2 diabetes compared with control individuals (p values <1.1 × 10-7). Sixty-four out of 76 (84.2%) CpG sites were confirmed by directionally consistent effects and p values <0.05 in an independent cohort of Indian Asians. However, on adjustment for baseline BMI only four CpG sites remained genome-wide significant, and addition of the 76 CpG methylation risk score to a prediction model including established predictors of type 2 diabetes (age, sex, BMI and HbA1c) showed no improvement (AUC 0.757 vs 0.753). Gene set enrichment analysis of the full epigenome-wide results clearly showed enrichment of processes linked to insulin signalling, lipid homeostasis and inflammation.

CONCLUSIONS/INTERPRETATION

By combining results from five European cohorts, and thus significantly increasing study sample size, we identified 76 CpG sites associated with incident type 2 diabetes. Replication of 64 CpGs in an independent cohort of Indian Asians suggests that the association between DNA methylation levels and incident type 2 diabetes is robust and independent of ethnicity. Our data also indicate that BMI partly explains the association between DNA methylation and incident type 2 diabetes. Further studies are required to elucidate the underlying biological mechanisms and to determine potential causal roles of the differentially methylated CpG sites in type 2 diabetes development.

Increase the Power of Epigenome-Wide Association Testing Using ICC-Based Hypothesis Weighting

For large-scale hypothesis testing such as epigenome-wide association testing, adaptively focusing power on the more promising hypotheses can lead to a much more powerful multiple testing procedure. In this chapter, we introduce a multiple testing procedure that weights each hypothesis based on the intraclass correlation coefficient (ICC), a measure of "noisiness" of CpG methylation measurement, to increase the power of epigenome-wide association testing. Compared to the traditional multiple testing procedure on a filtered CpG set, the proposed procedure circumvents the difficulty to determine the optimal ICC cutoff value and is overall more powerful. We illustrate the procedure and compare the power to classical multiple testing procedures using an example data.

Data Analysis of DNA Methylation Epigenome-Wide Association Studies (EWAS): A Guide to the Principles of Best Practice

Array-based EWAS have become an increasingly popular technique to identify population epigenetic effects, particularly in humans. With the arrival of nonhuman species arrays, such as the mouse, this is likely to become an even more widely used technology. This chapter provides the less experienced researcher a guide to the analysis of data from the most widely used platform, the Illumina Infinium Methylation assay. This includes an overview of quality filtering, data normalization, analysis options, and techniques to improve the interpretation of results.

Meta-analysis of epigenome-wide association studies of carotid intima-media thickness

Common carotid intima-media thickness (cIMT) is an index of subclinical atherosclerosis that is associated with ischemic stroke and coronary artery disease (CAD). We undertook a cross-sectional epigenome-wide association study (EWAS) of measures of cIMT in 6400 individuals. Mendelian randomization analysis was applied to investigate the potential causal role of DNA methylation in the link between atherosclerotic cardiovascular risk factors and cIMT or clinical cardiovascular disease. The CpG site cg05575921 was associated with cIMT (beta = -0.0264, p value = 3.5 × 10^-8) in the discovery panel and was replicated in replication panel (beta = -0.07, p value = 0.005). This CpG is located at chr5:81649347 in the intron 3 of the aryl hydrocarbon receptor repressor gene (AHRR). Our results indicate that DNA methylation at cg05575921 might be in the pathway between smoking, cIMT and stroke. Moreover, in a region-based analysis, 34 differentially methylated regions (DMRs) were identified of which a DMR upstream of ALOX12 showed the strongest association with cIMT (p value = 1.4 × 10^-13). In conclusion, our study suggests that DNA methylation may play a role in the link between cardiovascular risk factors, cIMT and clinical cardiovascular disease.

RnBeads 2.0: comprehensive analysis of DNA methylation data

DNA methylation is a widely investigated epigenetic mark with important roles in development and disease. High-throughput assays enable genome-scale DNA methylation analysis in large numbers of samples. Here, we describe a new version of our RnBeads software - an R/Bioconductor package that implements start-to-finish analysis workflows for Infinium microarrays and various types of bisulfite sequencing. RnBeads 2.0 (https://rnbeads.org/) provides additional data types and analysis methods, new functionality for interpreting DNA methylation differences, improved usability with a novel graphical user interface, and better use of computational resources. We demonstrate RnBeads 2.0 in four re-runnable use cases focusing on cell differentiation and cancer.

BayesCCE: a Bayesian framework for estimating cell-type composition from DNA methylation without the need for methylation reference

We introduce a Bayesian semi-supervised method for estimating cell counts from DNA methylation by leveraging an easily obtainable prior knowledge on the cell-type composition distribution of the studied tissue. We show mathematically and empirically that alternative methods which attempt to infer cell counts without methylation reference only capture linear combinations of cell counts rather than provide one component per cell type. Our approach allows the construction of components such that each component corresponds to a single cell type, and provides a new opportunity to investigate cell compositions in genomic studies of tissues for which it was not possible before.

Epigenome-wide association study of total serum immunoglobulin E in children: a life course approach

Background

IgE-mediated sensitization may be epigenetically programmed in utero, but early childhood environment may further alter complex traits and disease phenotypes through epigenetic plasticity. However, the epigenomic footprint underpinning IgE-mediated type-I hypersensitivity has not been well-understood, especially under a longitudinal early-childhood life-course framework.

Methods

We used epigenome-wide DNA methylation (IlluminaHumanMethylation450 BeadChip) in cord blood and mid-childhood peripheral blood to investigate pre- and post-natal methylation marks associated with mid-childhood (age 6.7-10.2) total serum IgE levels in 217 mother-child pairs in Project Viva-a prospective longitudinal pre-birth cohort in eastern Massachusetts, USA. We identified methylation sites associated with IgE using covariate-adjusted robust linear regressions.

Results

Nineteen methylation marks in cord blood were associated with IgE in mid-childhood (FDR < 0.05) in genes implicated in cell signaling, growth, and development. Among these, two methylation sites (C7orf50, ZAR1) remained robust after the adjustment for the change in DNA methylation from birth to mid-childhood (FDR < 0.05). An analysis of the change in methylation between cord blood and mid-childhood DNA (Δ-DNAm) identified 395 methylation marks in 272 genes associated with mid-childhood IgE (FDR < 0.05), with multiple sites located within ACOT7 (4 sites), EPX (5 sites), EVL (3 sites), KSR1 (4 sites), ZFPM1 (3 sites), and ZNF862 (3 sites). Several of these methylation loci were previously associated with asthma (ADAM19, EPX, IL4, IL5RA, and PRG2).

Conclusion

This study identified fetally programmed and mid-childhood methylation signals associated with mid-childhood IgE. Epigenetic priming during fetal development and early childhood likely plays an important role in IgE-mediated type-I hypersensitivity.

DNA methylation markers associated with type 2 diabetes, fasting glucose and HbA1c levels: a systematic review and replication in a case-control sample of the Lifelines study

AIMS/HYPOTHESIS

Epigenetic mechanisms may play an important role in the aetiology of type 2 diabetes. Recent epigenome-wide association studies (EWASs) identified several DNA methylation markers associated with type 2 diabetes, fasting glucose and HbA_1c levels. Here we present a systematic review of these studies and attempt to replicate the CpG sites (CpGs) with the most significant associations from these EWASs in a case-control sample of the Lifelines study.

METHODS

We performed a systematic literature search in PubMed and EMBASE for EWASs to test the association between DNA methylation and type 2 diabetes and/or glycaemic traits and reviewed the search results. For replication purposes we selected 100 unique CpGs identified in peripheral blood, pancreas, adipose tissue and liver from 15 EWASs, using study-specific Bonferroni-corrected significance thresholds. Methylation data (Illumina 450K array) in whole blood from 100 type 2 diabetic individuals and 100 control individuals from the Lifelines study were available. Multivariate linear models were used to examine the associations of the specific CpGs with type 2 diabetes and glycaemic traits.

RESULTS

From the 52 CpGs identified in blood and selected for replication, 15 CpGs showed nominally significant associations with type 2 diabetes in the Lifelines sample (p < 0.05). The results for five CpGs (in ABCG1, LOXL2, TXNIP, SLC1A5 and SREBF1) remained significant after a stringent multiple-testing correction (changes in methylation from -3% up to 3.6%, p < 0.0009). All associations were directionally consistent with the original EWAS results. None of the selected CpGs from the tissue-specific EWASs were replicated in our methylation data from whole blood. We were also unable to replicate any of the CpGs associated with HbA_1c levels in the healthy control individuals of our sample, while two CpGs (in ABCG1 and CCDC57) for fasting glucose were replicated at a nominal significance level (p < 0.05).

CONCLUSIONS/INTERPRETATION

A number of differentially methylated CpGs reported to be associated with type 2 diabetes in the EWAS literature were replicated in blood and show promise for clinical use as disease biomarkers. However, more prospective studies are needed to support the robustness of these findings.

A Guide to Illumina BeadChip Data Analysis

The Illumina Infinium BeadChips are a powerful array-based platform for genome-wide DNA methylation profiling at approximately 485,000 (450K) and 850,000 (EPIC) CpG sites across the genome. The platform is used in many large-scale population-based epigenetic studies of complex diseases, environmental exposures, or other experimental conditions. This chapter provides an overview of the key steps in analyzing Illumina BeadChip data. We describe key preprocessing steps including data extraction and quality control as well as normalization strategies. We further present principles and guidelines for conducting association analysis at the individual CpG level as well as more sophisticated pathway-based association tests.

"Gap hunting" to characterize clustered probe signals in Illumina methylation array data

Background

The Illumina 450k array has been widely used in epigenetic association studies. Current quality-control (QC) pipelines typically remove certain sets of probes, such as those containing a SNP or with multiple mapping locations. An additional set of potentially problematic probes are those with DNA methylation distributions characterized by two or more distinct clusters separated by gaps. Data-driven identification of such probes may offer additional insights for downstream analyses.

Results

We developed a procedure, termed “gap hunting,” to identify probes showing clustered distributions. Among 590 peripheral blood samples from the Study to Explore Early Development, we identified 11,007 “gap probes.” The vast majority (9199) are likely attributed to an underlying SNP(s) or other variant in the probe, although SNP-affected probes exist that do not produce a gap signals. Specific factors predict which SNPs lead to gap signals, including type of nucleotide change, probe type, DNA strand, and overall methylation state. These expected effects are demonstrated in paired genotype and 450k data on the same samples. Gap probes can also serve as a surrogate for the local genetic sequence on a haplotype scale and can be used to adjust for population stratification.

Conclusions

The characteristics of gap probes reflect potentially informative biology. QC pipelines may benefit from an efficient data-driven approach that “flags” gap probes, rather than filtering such probes, followed by careful interpretation of downstream association analyses. Our results should translate directly to the recently released Illumina EPIC array given the similar chemistry and content design.

Electronic supplementary material

The online version of this article (doi:10.1186/s13072-016-0107-z) contains supplementary material, which is available to authorized users.