DNA methylation signatures of chronic alcohol dependence in purified CD3+ T-cells of patients undergoing alcohol treatment

Several studies have shown an association of alcohol dependence with DNA methylation (DNAm), suggesting that environmentally-induced changes on epigenomic variation may play an important role in alcohol dependence. In the present study, we analysed genome-wide DNAm profiles of purified CD3⁺ T-cells from pre- and post-treatment alcohol dependent patients, as well as closely matched healthy controls. We identified 59 differentially methylated CpG sites comparing patients prior to treatment with healthy controls and were able to confirm 8 of those sites in additional analyses for differentially methylated regions. Comparing patients before and after a 3-week alcohol treatment program we revealed another unique set of 48 differentially methylated CpG sites. Additionally, we found that the mean global DNAm was significantly lower in patients prior to treatment compared to controls, but reverted back to levels similar to controls after treatment. We validated top-ranked hits derived from the epigenome-wide analysis by pyrosequencing and further replicated two of them in an independent cohort and confirmed differential DNAm of HECW2 and SRPK3 in whole blood. This study is the first to show widespread DNAm variation in a disease-relevant blood cell type and implicates HECW2 and SRPK3 DNAm as promising blood-based candidates to follow up in future studies.

Maternal blood contamination of collected cord blood can be identified using DNA methylation at three CpGs

Background

Cord blood is a commonly used tissue in environmental, genetic, and epigenetic population studies due to its ready availability and potential to inform on a sensitive period of human development. However, the introduction of maternal blood during labor or cross-contamination during sample collection may complicate downstream analyses. After discovering maternal contamination of cord blood in a cohort study of 150 neonates using Illumina 450K DNA methylation (DNAm) data, we used a combination of linear regression and random forest machine learning to create a DNAm-based screening method. We identified a panel of DNAm sites that could discriminate between contaminated and non-contaminated samples, then designed pyrosequencing assays to pre-screen DNA prior to being assayed on an array.

Results

Maternal contamination of cord blood was initially identified by unusual X chromosome DNA methylation patterns in 17 males. We utilized our DNAm panel to detect contaminated male samples and a proportional amount of female samples in the same cohort. We validated our DNAm screening method on an additional 189 sample cohort using both pyrosequencing and DNAm arrays, as well as 9 publically available cord blood 450K data sets. The rate of contamination varied from 0 to 10% within these studies, likely related to collection specific methods.

Conclusions

Maternal blood can contaminate cord blood during sample collection at appreciable levels across multiple studies. We have identified a panel of markers that can be used to identify this contamination, either post hoc after DNAm arrays have been completed, or in advance using a targeted technique like pyrosequencing.

Electronic supplementary material

The online version of this article (doi:10.1186/s13148-017-0370-2) contains supplementary material, which is available to authorized users.

Genome-Wide Interaction Analysis of Air Pollution Exposure and Childhood Asthma with Functional Follow-up

RATIONALE

The evidence supporting an association between traffic-related air pollution exposure and incident childhood asthma is inconsistent and may depend on genetic factors.

OBJECTIVES

To identify gene-environment interaction effects on childhood asthma using genome-wide single-nucleotide polymorphism (SNP) data and air pollution exposure. Identified loci were further analyzed at epigenetic and transcriptomic levels.

METHODS

We used land use regression models to estimate individual air pollution exposure (represented by outdoor NO₂ levels) at the birth address and performed a genome-wide interaction study for doctors' diagnoses of asthma up to 8 years in three European birth cohorts (n = 1,534) with look-up for interaction in two separate North American cohorts, CHS (Children's Health Study) and CAPPS/SAGE (Canadian Asthma Primary Prevention Study/Study of Asthma, Genetics and Environment) (n = 1,602 and 186 subjects, respectively). We assessed expression quantitative trait locus effects in human lung specimens and blood, as well as associations among air pollution exposure, methylation, and transcriptomic patterns.

MEASUREMENTS AND MAIN RESULTS

In the European cohorts, 186 SNPs had an interaction P < 1 × 10^-4 and a look-up evaluation of these disclosed 8 SNPs in 4 loci, with an interaction P < 0.05 in the large CHS study, but not in CAPPS/SAGE. Three SNPs within adenylate cyclase 2 (ADCY2) showed the same direction of the interaction effect and were found to influence ADCY2 gene expression in peripheral blood (P = 4.50 × 10^-4). One other SNP with P < 0.05 for interaction in CHS, rs686237, strongly influenced UDP-Gal:betaGlcNAc β-1,4-galactosyltransferase, polypeptide 5 (B4GALT5) expression in lung tissue (P = 1.18 × 10^-17). Air pollution exposure was associated with differential discs, large homolog 2 (DLG2) methylation and expression.

CONCLUSIONS

Our results indicated that gene-environment interactions are important for asthma development and provided supportive evidence for interaction with air pollution for ADCY2, B4GALT5, and DLG2.

The relation between DNA methylation patterns and serum cytokine levels in community-dwelling adults: a preliminary study

Background

The levels of circulating cytokines fluctuate with age, acute illness, and chronic disease, and are predictive of mortality; this is also true for patterns of DNA (CpG) methylation. Given that immune cells are particularly sensitive to changes in the concentration of cytokines in their microenvironment, we hypothesized that serum levels of TNF, IL-6, IL-8 and IL-10 would correlate with genome-wide alterations in the DNA methylation levels of blood leukocytes. To test this, we evaluated community-dwelling adults (n = 14; 48–78 years old) recruited to a pilot study for the Canadian Longitudinal Study on Aging (CLSA), examining DNA methylation patterns in peripheral blood mononuclear cells using the Illumina HumanMethylation 450 K BeadChip.

Results

We show that, apart from age, serum IL-10 levels exhibited the most substantial association to DNA methylation patterns, followed by TNF, IL-6 and IL-8. Furthermore, while the levels of these cytokines were higher in elderly adults, no associations with epigenetic accelerated aging, derived using the epigenetic clock, were observed.

Conclusions

As a preliminary study with a small sample size, the conclusions drawn from this work must be viewed with caution; however, our observations are encouraging and certainly warrant more suitably powered studies of this relationship.

Electronic supplementary material

The online version of this article (doi:10.1186/s12863-017-0525-3) contains supplementary material, which is available to authorized users.

Associations between maternal risk factors of adverse pregnancy and birth outcomes and the offspring epigenetic clock of gestational age at birth

Background

A recent study has shown that it is possible to accurately estimate gestational age (GA) at birth from the DNA methylation (DNAm) of fetal umbilical cord blood/newborn blood spots. This DNAm GA predictor may provide additional information relevant to developmental stage. In 814 mother-neonate pairs, we evaluated the associations between DNAm GA and a number of maternal and offspring characteristics. These characteristics reflect prenatal environmental adversity and are expected to influence newborn developmental stage.

Results

DNAm GA acceleration (GAA; i.e., older DNAm GA than chronological GA) of the offspring at birth was associated with maternal age of over 40 years at delivery, pre-eclampsia and fetal demise in a previous pregnancy, maternal pre-eclampsia and treatment with antenatal betamethasone in the index pregnancy, lower neonatal birth size, lower 1-min Apgar score, and female sex. DNAm GA deceleration (GAD; i.e., younger DNAm GA than chronological GA) of the offspring at birth was associated with insulin-treated gestational diabetes mellitus (GDM) in a previous pregnancy and Sjögren’s syndrome. These findings were more accentuated when the DNAm GA calculation was based on the raw difference between DNAm GA and GA than on the residual from the linear regression of DNAm GA on GA.

Conclusions

Our findings show that variations in the DNAm GA of the offspring at birth are associated with a number of maternal and offspring characteristics known to reflect exposure to prenatal environmental adversity. Future studies should be aimed at determining if this biological variation is predictive of developmental adversity.

Electronic supplementary material

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

Differential DNA methylation and lymphocyte proportions in a Costa Rican high longevity region

Background

The Nicoya Peninsula in Costa Rica has one of the highest old-age life expectancies in the world, but the underlying biological mechanisms of this longevity are not well understood. As DNA methylation is hypothesized to be a component of biological aging, we focused on this malleable epigenetic mark to determine its association with current residence in Nicoya versus elsewhere in Costa Rica. Examining a population’s unique DNA methylation pattern allows us to differentiate hallmarks of longevity from individual stochastic variation. These differences may be characteristic of a combination of social, biological, and environmental contexts.

Methods

In a cross-sectional subsample of the Costa Rican Longevity and Healthy Aging Study, we compared whole blood DNA methylation profiles of residents from Nicoya (n = 48) and non-Nicoya (other Costa Rican regions, n = 47) using the Infinium HumanMethylation450 microarray.

Results

We observed a number of differences that may be markers of delayed aging, such as bioinformatically derived differential CD8+ T cell proportions. Additionally, both site- and region-specific analyses revealed DNA methylation patterns unique to Nicoyans. We also observed lower overall variability in DNA methylation in the Nicoyan population, another hallmark of younger biological age.

Conclusions

Nicoyans represent an interesting group of individuals who may possess unique immune cell proportions as well as distinct differences in their epigenome, at the level of DNA methylation.

Electronic supplementary material

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

Developmental pathways to adiposity begin before birth and are influenced by genotype, prenatal environment and epigenome

BACKGROUND

Obesity is an escalating health problem worldwide, and hence the causes underlying its development are of primary importance to public health. There is growing evidence that suboptimal intrauterine environment can perturb the metabolic programing of the growing fetus, thereby increasing the risk of developing obesity in later life. However, the link between early exposures in the womb, genetic susceptibility, and perturbed epigenome on metabolic health is not well understood. In this study, we shed more light on this aspect by performing a comprehensive analysis on the effects of variation in prenatal environment, neonatal methylome, and genotype on birth weight and adiposity in early childhood.

METHODS

In a prospective mother-offspring cohort (N = 987), we interrogated the effects of 30 variables that influence the prenatal environment, umbilical cord DNA methylation, and genotype on offspring weight and adiposity, over the period from birth to 48 months. This is an interim analysis on an ongoing cohort study.

RESULTS

Eleven of 30 prenatal environments, including maternal adiposity, smoking, blood glucose and plasma unsaturated fatty acid levels, were associated with birth weight. Polygenic risk scores derived from genetic association studies on adult adiposity were also associated with birth weight and child adiposity, indicating an overlap between the genetic pathways influencing metabolic health in early and later life. Neonatal methylation markers from seven gene loci (ANK3, CDKN2B, CACNA1G, IGDCC4, P4HA3, ZNF423 and MIRLET7BHG) were significantly associated with birth weight, with a subset of these in genes previously implicated in metabolic pathways in humans and in animal models. Methylation levels at three of seven birth weight-linked loci showed significant association with prenatal environment, but none were affected by polygenic risk score. Six of these birth weight-linked loci continued to show a longitudinal association with offspring size and/or adiposity in early childhood.

CONCLUSIONS

This study provides further evidence that developmental pathways to adiposity begin before birth and are influenced by environmental, genetic and epigenetic factors. These pathways can have a lasting effect on offspring size, adiposity and future metabolic outcomes, and offer new opportunities for risk stratification and prevention of obesity.

CLINICAL TRIAL REGISTRATION

This birth cohort is a prospective observational study, designed to study the developmental origins of health and disease, and was retrospectively registered on 1 July 2010 under the identifier NCT01174875 .

Epigenetics studies of fetal alcohol spectrum disorder: where are we now?

Adverse in utero events can alter the development and function of numerous physiological systems, giving rise to lasting neurodevelopmental deficits. In particular, data have shown that prenatal alcohol exposure can reprogram neurobiological systems, altering developmental trajectories and resulting in increased vulnerability to adverse neurobiological, behavioral and health outcomes. Increasing evidence suggests that epigenetic mechanisms are potential mediators for the reprogramming of neurobiological systems, as they may provide a link between the genome, environmental conditions and neurodevelopmental outcomes. This review outlines the current state of epigenetic research in fetal alcohol spectrum disorder, highlighting the role of epigenetic mechanisms in the reprogramming of neurobiological systems by alcohol and as potential diagnostic tools for fetal alcohol spectrum disorder. We also present an assessment of the current limitations in studies of prenatal alcohol exposure, and highlight the future steps needed in the field.

An empirically driven data reduction method on the human 450K methylation array to remove tissue specific non-variable CpGs

Background

Population based epigenetic association studies of disease and exposures are becoming more common with the availability of economical genome-wide technologies for interrogation of the methylome, such as the Illumina 450K Human Methylation Array (450K). Often, the expected small number of differentially methylated cytosine-guanine pairs (CpGs) in studies of the human methylome presents a statistical challenge, as the large number of CpGs measured on the 450K necessitates careful multiple test correction. While the 450K is a highly useful tool for population epigenetic studies, many of the CpGs tested are not variable and thus of limited information content in the context of the study and tissue. CpGs with observed lack of variability in the tissue under study could be removed to reduce the data dimensionality, limit the severity of multiple test correction and allow for improved detection of differential DNA methylation.

Methods

Here, we performed a meta-analysis of 450K data from three commonly studied human tissues, namely blood (605 samples), buccal epithelial cells (121 samples) and placenta (157 samples). We developed lists of CpGs that are non-variable in each tissue.

Results

These lists are surprisingly large (blood 114,204 CpGs, buccal epithelial cells 120,009 CpGs and placenta 101,367 CpGs) and thus will be valuable filters for epigenetic association studies, considerably reducing the dimensionality of the 450K and subsequently the multiple testing correction severity.

Conclusions

We propose this empirically derived method for data reduction to allow for more power in detecting differential DNA methylation associated with exposures in studies on the human methylome.

Electronic supplementary material

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

Differential DNA methylation in peripheral blood mononuclear cells in adolescents exposed to significant early but not later childhood adversity

Internationally adopted adolescents who are adopted as young children from conditions of poverty and deprivation have poorer physical and mental health outcomes than do adolescents conceived, born, and raised in the United States by families similar to those who adopt internationally. Using a sample of Russian and Eastern European adoptees to control for Caucasian race and US birth, and nonadopted offspring of well-educated and well-resourced parents to control for postadoption conditions, we hypothesized that the important differences in environments, conception to adoption, might be reflected in epigenetic patterns between groups, specifically in DNA methylation. Thus, we conducted an epigenome-wide association study to compare DNA methylation profiles at approximately 416,000 individual CpG loci from peripheral blood mononuclear cells of 50 adopted youth and 33 nonadopted youth. Adopted youth averaged 22 months at adoption, and both groups averaged 15 years at testing; thus, roughly 80% of their lives were lived in similar circumstances. Although concurrent physical health did not differ, cell-type composition predicted using the DNA methylation data revealed a striking difference in the white blood cell-type composition of the adopted and nonadopted youth. After correcting for cell type and removing invariant probes, 30 CpG sites in 19 genes were more methylated in the adopted group. We also used an exploratory functional analysis that revealed that 223 gene ontology terms, clustered in neural and developmental categories, were significantly enriched between groups.

An epigenetic clock analysis of race/ethnicity, sex, and coronary heart disease

Background

Epigenetic biomarkers of aging (the “epigenetic clock”) have the potential to address puzzling findings surrounding mortality rates and incidence of cardio-metabolic disease such as: (1) women consistently exhibiting lower mortality than men despite having higher levels of morbidity; (2) racial/ethnic groups having different mortality rates even after adjusting for socioeconomic differences; (3) the black/white mortality cross-over effect in late adulthood; and (4) Hispanics in the United States having a longer life expectancy than Caucasians despite having a higher burden of traditional cardio-metabolic risk factors.

Results

We analyzed blood, saliva, and brain samples from seven different racial/ethnic groups. We assessed the intrinsic epigenetic age acceleration of blood (independent of blood cell counts) and the extrinsic epigenetic aging rates of blood (dependent on blood cell counts and tracks the age of the immune system). In blood, Hispanics and Tsimane Amerindians have lower intrinsic but higher extrinsic epigenetic aging rates than Caucasians. African-Americans have lower extrinsic epigenetic aging rates than Caucasians and Hispanics but no differences were found for the intrinsic measure. Men have higher epigenetic aging rates than women in blood, saliva, and brain tissue.

Conclusions

Epigenetic aging rates are significantly associated with sex, race/ethnicity, and to a lesser extent with CHD risk factors, but not with incident CHD outcomes. These results may help elucidate lower than expected mortality rates observed in Hispanics, older African-Americans, and women.

Electronic supplementary material

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

DNA methylation signature of human fetal alcohol spectrum disorder

Background

Prenatal alcohol exposure is the leading preventable cause of behavioral and cognitive deficits, which may affect between 2 and 5 % of children in North America. While the underlying mechanisms of alcohol’s effects on development remain relatively unknown, emerging evidence implicates epigenetic mechanisms in mediating the range of symptoms observed in children with fetal alcohol spectrum disorder (FASD). Thus, we investigated the effects of prenatal alcohol exposure on genome-wide DNA methylation in the NeuroDevNet FASD cohort, the largest cohort of human FASD samples to date.

Methods

Genome-wide DNA methylation patterns of buccal epithelial cells (BECs) were analyzed using the Illumina HumanMethylation450 array in a Canadian cohort of 206 children (110 FASD and 96 controls). Genotyping was performed in parallel using the Infinium HumanOmni2.5-Quad v1.0 BeadChip.

Results

After correcting for the effects of genetic background, we found 658 significantly differentially methylated sites between FASD cases and controls, with 41 displaying differences in percent methylation change >5 %. Furthermore, 101 differentially methylated regions containing two or more CpGs were also identified, overlapping with 95 different genes. The majority of differentially methylated genes were highly expressed at the level of mRNA in brain samples from the Allen Brain Atlas, and independent DNA methylation data from cortical brain samples showed high correlations with BEC DNA methylation patterns. Finally, overrepresentation analysis of genes with up-methylated CpGs revealed a significant enrichment for neurodevelopmental processes and diseases, such as anxiety, epilepsy, and autism spectrum disorders.

Conclusions

These findings suggested that prenatal alcohol exposure is associated with distinct DNA methylation patterns in children and adolescents, raising the possibility of an epigenetic biomarker of FASD.

Electronic supplementary material

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

Inhalation of diesel exhaust and allergen alters human bronchial epithelium DNA methylation

BACKGROUND

Allergic disease affects 30% to 40% of the world's population, and its development is determined by the interplay between environmental and inherited factors. Air pollution, primarily consisting of diesel exhaust emissions, has increased at a similar rate to allergic disease. Exposure to diesel exhaust may play a role in the development and progression of allergic disease, in particular allergic respiratory disease. One potential mechanism underlying the connection between air pollution and increased allergic disease incidence is DNA methylation, an epigenetic process with the capacity to integrate gene-environment interactions.

OBJECTIVE

We sought to investigate the effect of allergen and diesel exhaust exposure on bronchial epithelial DNA methylation.

METHODS

We performed a randomized crossover-controlled exposure study to allergen and diesel exhaust in humans, and measured single-site (CpG) resolution global DNA methylation in bronchial epithelial cells.

RESULTS

Exposure to allergen alone, diesel exhaust alone, or allergen and diesel exhaust together (coexposure) led to significant changes in 7 CpG sites at 48 hours. However, when the same lung was exposed to allergen and diesel exhaust but separated by approximately 4 weeks, significant changes in more than 500 sites were observed. Furthermore, sites of differential methylation differed depending on which exposure was experienced first. Functional analysis of differentially methylated CpG sites found genes involved in transcription factor activity, protein metabolism, cell adhesion, and vascular development, among others.

CONCLUSIONS

These findings suggest that specific exposures can prime the lung for changes in DNA methylation induced by a subsequent insult.

A single flexible RNAPII-CTD integrates many different transcriptional programs

The RNAPII-CTD functions as a binding platform for coordinating the recruitment of transcription associated factors. Altering CTD function results in gene expression defects, although mounting evidence suggests that these effects likely vary among species and loci. Here we highlight emerging evidence of species- and loci-specific functions for the RNAPII-CTD.

Improving cell mixture deconvolution by identifying optimal DNA methylation libraries (IDOL)

Background

Confounding due to cellular heterogeneity represents one of the foremost challenges currently facing Epigenome-Wide Association Studies (EWAS). Statistical methods leveraging the tissue-specificity of DNA methylation for deconvoluting the cellular mixture of heterogenous biospecimens offer a promising solution, however the performance of such methods depends entirely on the library of methylation markers being used for deconvolution. Here, we introduce a novel algorithm for Identifying Optimal Libraries (IDOL) that dynamically scans a candidate set of cell-specific methylation markers to find libraries that optimize the accuracy of cell fraction estimates obtained from cell mixture deconvolution.

Results

Application of IDOL to training set consisting of samples with both whole-blood DNA methylation data (Illumina HumanMethylation450 BeadArray (HM450)) and flow cytometry measurements of cell composition revealed an optimized library comprised of 300 CpG sites. When compared existing libraries, the library identified by IDOL demonstrated significantly better overall discrimination of the entire immune cell landscape (p = 0.038), and resulted in improved discrimination of 14 out of the 15 pairs of leukocyte subtypes. Estimates of cell composition across the samples in the training set using the IDOL library were highly correlated with their respective flow cytometry measurements, with all cell-specific R²>0.99 and root mean square errors (RMSEs) ranging from [0.97 % to 1.33 %] across leukocyte subtypes. Independent validation of the optimized IDOL library using two additional HM450 data sets showed similarly strong prediction performance, with all cell-specific R²>0.90 and RMSE<4.00 %. In simulation studies, adjustments for cell composition using the IDOL library resulted in uniformly lower false positive rates compared to competing libraries, while also demonstrating an improved capacity to explain epigenome-wide variation in DNA methylation within two large publicly available HM450 data sets.

Conclusions

Despite consisting of half as many CpGs compared to existing libraries for whole blood mixture deconvolution, the optimized IDOL library identified herein resulted in outstanding prediction performance across all considered data sets and demonstrated potential to improve the operating characteristics of EWAS involving adjustments for cell distribution. In addition to providing the EWAS community with an optimized library for whole blood mixture deconvolution, our work establishes a systematic and generalizable framework for the assembly of libraries that improve the accuracy of cell mixture deconvolution.

Electronic supplementary material

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

DNA methylation profiling in human Huntington's disease brain

Despite extensive progress in Huntington's disease (HD) research, very little is known about the association of epigenetic variation and HD pathogenesis in human brain tissues. Moreover, its contribution to the tissue-specific transcriptional regulation of the huntingtin gene (HTT), in which HTT expression levels are highest in brain and testes, is currently unknown. To investigate the role of DNA methylation in HD pathogenesis and tissue-specific expression of HTT, we utilized the Illumina HumanMethylation450K BeadChip array to measure DNA methylation in a cohort of age-matched HD and control human cortex and liver tissues. In cortex samples, we found minimal evidence of HD-associated DNA methylation at probed sites after correction for cell heterogeneity but did observe an association with the age of disease onset. In contrast, comparison of matched cortex and liver samples revealed tissue-specific DNA methylation of the HTT gene region at 38 sites (FDR < 0.05). Importantly, we identified a novel differentially methylated binding site in the HTT proximal promoter for the transcription factor CTCF. This CTCF site displayed increased occupancy in cortex, where HTT expression is higher, compared with the liver. Additionally, CTCF silencing reduced the activity of an HTT promoter-reporter construct, suggesting that CTCF plays a role in regulating HTT promoter function. Overall, although we were unable to detect HD-associated DNA methylation alterations at queried sites, we found that DNA methylation may be correlated to the age of disease onset in cortex tissues. Moreover, our data suggest that DNA methylation may, in part, contribute to tissue-specific HTT transcription through differential CTCF occupancy.

Comparison of Methyl-capture Sequencing vs. Infinium 450K methylation array for methylome analysis in clinical samples

Interindividual variability in the epigenome has gained tremendous attention for its potential in pathophysiological investigation, disease diagnosis, and evaluation of clinical intervention. DNA methylation is the most studied epigenetic mark in epigenome-wide association studies (EWAS) as it can be detected from limited starting material. Infinium 450K methylation array is the most popular platform for high-throughput profiling of this mark in clinical samples, as it is cost-effective and requires small amounts of DNA. However, this method suffers from low genome coverage and errors introduced by probe cross-hybridization. Whole-genome bisulfite sequencing can overcome these limitations but elevates the costs tremendously. Methyl-Capture Sequencing (MC Seq) is an attractive intermediate solution to increase the methylome coverage in large sample sets. Here we first demonstrate that MC Seq can be employed using DNA amounts comparable to the amounts used for Infinium 450K. Second, to provide guidance when choosing between the 2 platforms for EWAS, we evaluate and compare MC Seq and Infinium 450K in terms of coverage, technical variation, and concordance of methylation calls in clinical samples. Last, since the focus in EWAS is to study interindividual variation, we demonstrate the utility of MC Seq in studying interindividual variation in subjects from different ethnicities.

Natural variation in maternal care and cross-tissue patterns of oxytocin receptor gene methylation in rats

This article is part of a Special Issue "Parental Care". Since the first report of maternal care effects on DNA methylation in rats, epigenetic modifications of the genome in response to life experience have become the subject of intense focus across many disciplines. Oxytocin receptor expression varies in response to early experience, and both oxytocin signaling and methylation status of the oxytocin receptor gene (Oxtr) in blood have been related to disordered social behavior. It is unknown whether Oxtr DNA methylation varies in response to early life experience, and whether currently employed peripheral measures of Oxtr methylation reflect variation in the brain. We examined the effects of early life rearing experience via natural variation in maternal licking and grooming during the first week of life on behavior, physiology, gene expression, and epigenetic regulation of Oxtr across blood and brain tissues (mononucleocytes, hippocampus, striatum, and hypothalamus). Rats reared by "high" licking-grooming (HL) and "low" licking-grooming (LL) rat dams exhibited differences across study outcomes: LL offspring were more active in behavioral arenas, exhibited lower body mass in adulthood, and showed reduced corticosterone responsivity to a stressor. Oxtr DNA methylation was significantly lower at multiple CpGs in the blood of LL versus HL males, but no differences were found in the brain. Across groups, Oxtr transcript levels in the hypothalamus were associated with reduced corticosterone secretion in response to stress, congruent with the role of oxytocin signaling in this region. Methylation of specific CpGs at a high or low level was consistent across tissues, especially within the brain. However, individual variation in DNA methylation relative to these global patterns was not consistent across tissues. These results suggest that blood Oxtr DNA methylation may reflect early experience of maternal care, and that Oxtr methylation across tissues is highly concordant for specific CpGs, but that inferences across tissues are not supported for individual variation in Oxtr methylation.

DNA methylation and healthy human aging

The process of aging results in a host of changes at the cellular and molecular levels, which include senescence, telomere shortening, and changes in gene expression. Epigenetic patterns also change over the lifespan, suggesting that epigenetic changes may constitute an important component of the aging process. The epigenetic mark that has been most highly studied is DNA methylation, the presence of methyl groups at CpG dinucleotides. These dinucleotides are often located near gene promoters and associate with gene expression levels. Early studies indicated that global levels of DNA methylation increase over the first few years of life and then decrease beginning in late adulthood. Recently, with the advent of microarray and next‐generation sequencing technologies, increases in variability of DNA methylation with age have been observed, and a number of site‐specific patterns have been identified. It has also been shown that certain CpG sites are highly associated with age, to the extent that prediction models using a small number of these sites can accurately predict the chronological age of the donor. Together, these observations point to the existence of two phenomena that both contribute to age‐related DNA methylation changes: epigenetic drift and the epigenetic clock. In this review, we focus on healthy human aging throughout the lifetime and discuss the dynamics of DNA methylation as well as how interactions between the genome, environment, and the epigenome influence aging rates. We also discuss the impact of determining ‘epigenetic age’ for human health and outline some important caveats to existing and future studies.

The epigenomic landscape of African rainforest hunter-gatherers and farmers

The genetic history of African populations is increasingly well documented, yet their patterns of epigenomic variation remain uncharacterized. Moreover, the relative impacts of DNA sequence variation and temporal changes in lifestyle and habitat on the human epigenome remain unknown. Here we generate genome-wide genotype and DNA methylation profiles for 362 rainforest hunter-gatherers and sedentary farmers. We find that the current habitat and historical lifestyle of a population have similarly critical impacts on the methylome, but the biological functions affected strongly differ. Specifically, methylation variation associated with recent changes in habitat mostly concerns immune and cellular functions, whereas that associated with historical lifestyle affects developmental processes. Furthermore, methylation variation—particularly that correlated with historical lifestyle—shows strong associations with nearby genetic variants that, moreover, are enriched in signals of natural selection. Our work provides new insight into the genetic and environmental factors affecting the epigenomic landscape of human populations over time.

Genetic and environmental factors affect genome-wide patterns of epigenetic variation. Here, the authors show that while current habitat and historical lifestyle impact the methylome of rainforest hunter-gatherers and sedentary farmers, the biological functions affected and the degree of genetic control differ.

Prenatal alcohol exposure alters steady-state and activated gene expression in the adult rat brain

BACKGROUND

Prenatal alcohol exposure (PAE) is associated with alterations in numerous physiological systems, including the stress and immune systems. We have previously shown that PAE increases the course and severity of arthritis in an adjuvant-induced arthritis (AA) model. While the molecular mechanisms underlying these effects are not fully known, changes in neural gene expression are emerging as important factors in the etiology of PAE effects. As the prefrontal cortex (PFC) and hippocampus (HPC) play key roles in neuroimmune function, PAE-induced alterations to their transcriptome may underlie abnormal steady-state functions and responses to immune challenge. This study examined brains from adult PAE and control females from our recent AA study to determine whether PAE causes long-term alterations in gene expression and whether these mediate the altered severity and course of arthritis in PAE females.

METHODS

Adult females from PAE, pair-fed (PF), and ad libitum-fed control (C) groups were injected with either saline or complete Freund's adjuvant. Animals were terminated at the peak of inflammation or during resolution (Days 16 and 39 postinjection, respectively); cohorts of saline-injected PAE, PF, and C females were terminated in parallel. Gene expression was analyzed in the PFC and HPC using whole-genome mRNA expression microarrays.

RESULTS

Significant changes in gene expression in both the PFC and HPC were found in PAE compared to controls in response to ethanol exposure alone (saline-injected females), including genes involved in neurodevelopment, apoptosis, and energy metabolism. Moreover, in response to inflammation (adjuvant-injected females), PAE animals showed unique expression patterns, while failing to exhibit the activation of genes and regulators involved in the immune response observed in control and pair-fed animals.

CONCLUSIONS

These results support the hypothesis that PAE affects neuroimmune function at the level of gene expression, demonstrating long-term effects of PAE on the central nervous system response under steady-state conditions and following an inflammatory insult.

funtooNorm: an R package for normalization of DNA methylation data when there are multiple cell or tissue types

MOTIVATION

DNA methylation patterns are well known to vary substantially across cell types or tissues. Hence, existing normalization methods may not be optimal if they do not take this into account. We therefore present a new R package for normalization of data from the Illumina Infinium Human Methylation450 BeadChip (Illumina 450 K) built on the concepts in the recently published funNorm method, and introducing cell-type or tissue-type flexibility.

RESULTS

funtooNorm is relevant for data sets containing samples from two or more cell or tissue types. A visual display of cross-validated errors informs the choice of the optimal number of components in the normalization. Benefits of cell (tissue)-specific normalization are demonstrated in three data sets. Improvement can be substantial; it is strikingly better on chromosome X, where methylation patterns have unique inter-tissue variability.

AVAILABILITY AND IMPLEMENTATION

An R package is available at https://github.com/GreenwoodLab/funtooNorm, and has been submitted to Bioconductor at http://bioconductor.org.

The RNAPII-CTD Maintains Genome Integrity through Inhibition of Retrotransposon Gene Expression and Transposition

RNA polymerase II (RNAPII) contains a unique C-terminal domain that is composed of heptapeptide repeats and which plays important regulatory roles during gene expression. RNAPII is responsible for the transcription of most protein-coding genes, a subset of non-coding genes, and retrotransposons. Retrotransposon transcription is the first step in their multiplication cycle, given that the RNA intermediate is required for the synthesis of cDNA, the material that is ultimately incorporated into a new genomic location. Retrotransposition can have grave consequences to genome integrity, as integration events can change the gene expression landscape or lead to alteration or loss of genetic information. Given that RNAPII transcribes retrotransposons, we sought to investigate if the RNAPII-CTD played a role in the regulation of retrotransposon gene expression. Importantly, we found that the RNAPII-CTD functioned to maintaining genome integrity through inhibition of retrotransposon gene expression, as reducing CTD length significantly increased expression and transposition rates of Ty1 elements. Mechanistically, the increased Ty1 mRNA levels in the rpb1-CTD11 mutant were partly due to Cdk8-dependent alterations to the RNAPII-CTD phosphorylation status. In addition, Cdk8 alone contributed to Ty1 gene expression regulation by altering the occupancy of the gene-specific transcription factor Ste12. Loss of STE12 and TEC1 suppressed growth phenotypes of the RNAPII-CTD truncation mutant. Collectively, our results implicate Ste12 and Tec1 as general and important contributors to the Cdk8, RNAPII-CTD regulatory circuitry as it relates to the maintenance of genome integrity.