Immune cell type and DNA methylation vary with reproductive status in women: possible pathways for costs of reproduction

Background

Consistent with evolutionarily theorized costs of reproduction (CoR), reproductive history in women is associated with life expectancy and susceptibility to certain cancers, autoimmune disorders and metabolic disease. Immunological changes originating during reproduction may help explain some of these relationships.

Methodology

To explore the potential role of the immune system in female CoR, we characterized leukocyte composition and regulatory processes using DNA methylation (DNAm) in a cross-sectional cohort of young (20–22 years old) women differing in reproductive status.

Results

Compared to nulliparity, pregnancy was characterized by differential methylation at 828 sites, 96% of which were hypomethylated and enriched for genes associated with T-cell activation, innate immunity, pre-eclampsia and neoplasia. Breastfeeding was associated with differential methylation at 1107 sites (71% hypermethylated), enriched for genes involved in metabolism, immune self-recognition and neurogenesis. There were no significant differences in DNAm between nulliparous and parous women. However, compared to nullipara, pregnant women had lower proportions of B, CD4T, CD8T and natural killer (NK) cells, and higher proportions of granulocytes and monocytes. Monocyte counts were lower and NK counts higher among breastfeeding women, and remained so among parous women.

Implications

Our findings point to widespread differences in DNAm during pregnancy and lactation. These effects appear largely transient, but may accumulate with gravidity become detectable as women age. Nulliparous and parous women differed in leukocyte composition, consistent with more persistent effects of reproduction on cell type. These findings support transient (leukocyte DNAm) and persistent (cell composition) changes associated with reproduction in women, illuminating potential pathways contributing to CoR.

Lay Summary: Evolutionary theory and epidemiology support costs of reproduction (CoR) to women’s health that may involve changes in immune function. We report differences in immune cell composition and gene regulation during pregnancy and breastfeeding. While many of these differences appear transient, immune cell composition may remain, suggesting mechanisms for female CoR.

Culture-Associated DNA Methylation Changes Impact on Cellular Function of Human Intestinal Organoids

BACKGROUND & AIMS

Human intestinal epithelial organoids (IEOs) are a powerful tool to model major aspects of intestinal development, health, and diseases because patient-derived cultures retain many features found in vivo. A necessary aspect of the organoid model is the requirement to expand cultures in vitro through several rounds of passaging. This is of concern because the passaging of cells has been shown to affect cell morphology, ploidy, and function.

METHODS

Here, we analyzed 173 human IEO lines derived from the small and large bowel and examined the effect of culture duration on DNA methylation (DNAm). Furthermore, we tested the potential impact of DNAm changes on gene expression and cellular function.

RESULTS

Our analyses show a reproducible effect of culture duration on DNAm in a large discovery cohort as well as 2 publicly available validation cohorts generated in different laboratories. Although methylation changes were seen in only approximately 8% of tested cytosine-phosphate-guanine dinucleotides (CpGs) and global cellular function remained stable, a subset of methylation changes correlated with altered gene expression at baseline as well as in response to inflammatory cytokine exposure and withdrawal of Wnt agonists. Importantly, epigenetic changes were found to be enriched in genomic regions associated with colonic cancer and distant to the site of replication, indicating similarities to malignant transformation.

CONCLUSIONS

Our study shows distinct culture-associated epigenetic changes in mucosa-derived human IEOs, some of which appear to impact gene transcriptomic and cellular function. These findings highlight the need for future studies in this area and the importance of considering passage number as a potentially confounding factor.

Transcription and DNA Methylation Patterns of Blood-Derived CD8+ T Cells Are Associated With Age and Inflammatory Bowel Disease But Do Not Predict Prognosis

BACKGROUND & AIMS

Gene expression patterns of CD8⁺ T cells have been reported to correlate with clinical outcomes of adults with inflammatory bowel diseases (IBD). We aimed to validate these findings in independent patient cohorts.

METHODS

We obtained peripheral blood samples from 112 children with a new diagnosis of IBD (71 with Crohn's disease and 41 with ulcerative colitis) and 19 children without IBD (controls) and recorded medical information on disease activity and outcomes. CD8⁺ T cells were isolated from blood samples by magnetic bead sorting at the point of diagnosis and during the course of disease. Genome-wide transcription (n = 192) and DNA methylation (n = 66) profiles were generated using Affymetrix and Illumina arrays, respectively. Publicly available transcriptomes and DNA methylomes of CD8⁺ T cells from 3 adult patient cohorts with and without IBD were included in data analyses.

RESULTS

Previously reported CD8⁺ T-cell prognostic expression and exhaustion signatures were only found in the original adult IBD patient cohort. These signatures could not be detected in either a pediatric or a second adult IBD cohort. In contrast, an association between CD8⁺ T-cell gene expression with age and sex was detected across all 3 cohorts. CD8⁺ gene transcription was clearly associated with IBD in the 2 cohorts that included non-IBD controls. Lastly, DNA methylation profiles of CD8⁺ T cells from children with Crohn's disease correlated with age but not with disease outcome.

CONCLUSIONS

We were unable to validate previously reported findings of an association between CD8⁺ T-cell gene transcription and disease outcome in IBD. Our findings reveal the challenges of developing prognostic biomarkers for patients with IBD and the importance of their validation in large, independent cohorts before clinical application.

Multi-Omic Data Integration Allows Baseline Immune Signatures to Predict Hepatitis B Vaccine Response in a Small Cohort

Background

Vaccination remains one of the most effective means of reducing the burden of infectious diseases globally. Improving our understanding of the molecular basis for effective vaccine response is of paramount importance if we are to ensure the success of future vaccine development efforts.

Methods

We applied cutting edge multi-omics approaches to extensively characterize temporal molecular responses following vaccination with hepatitis B virus (HBV) vaccine. Data were integrated across cellular, epigenomic, transcriptomic, proteomic, and fecal microbiome profiles, and correlated to final HBV antibody titres.

Results

Using both an unsupervised molecular-interaction network integration method (NetworkAnalyst) and a data-driven integration approach (DIABLO), we uncovered baseline molecular patterns and pathways associated with more effective vaccine responses to HBV. Biological associations were unravelled, with signalling pathways such as JAK-STAT and interleukin signalling, Toll-like receptor cascades, interferon signalling, and Th17 cell differentiation emerging as important pre-vaccination modulators of response.

Conclusion

This study provides further evidence that baseline cellular and molecular characteristics of an individual's immune system influence vaccine responses, and highlights the utility of integrating information across many parallel molecular datasets.

Systems Biology Methods Applied to Blood and Tissue for a Comprehensive Analysis of Immune Response to Hepatitis B Vaccine in Adults

Conventional vaccine design has been based on trial-and-error approaches, which have been generally successful. However, there have been some major failures in vaccine development and we still do not have highly effective licensed vaccines for tuberculosis, HIV, respiratory syncytial virus, and other major infections of global significance. Approaches at rational vaccine design have been limited by our understanding of the immune response to vaccination at the molecular level. Tools now exist to undertake in-depth analysis using systems biology approaches, but to be fully realized, studies are required in humans with intensive blood and tissue sampling. Methods that support this intensive sampling need to be developed and validated as feasible. To this end, we describe here a detailed approach that was applied in a study of 15 healthy adults, who were immunized with hepatitis B vaccine. Sampling included ~350 mL of blood, 12 microbiome samples, and lymph node fine needle aspirates obtained over a ~7-month period, enabling comprehensive analysis of the immune response at the molecular level, including single cell and tissue sample analysis. Samples were collected for analysis of immune phenotyping, whole blood and single cell gene expression, proteomics, lipidomics, epigenetics, whole blood response to key immune stimuli, cytokine responses, in vitro T cell responses, antibody repertoire analysis and the microbiome. Data integration was undertaken using different approaches-NetworkAnalyst and DIABLO. Our results demonstrate that such intensive sampling studies are feasible in healthy adults, and data integration tools exist to analyze the vast amount of data generated from a multi-omics systems biology approach. This will provide the basis for a better understanding of vaccine-induced immunity and accelerate future rational vaccine design.

The PedBE clock accurately estimates DNA methylation age in pediatric buccal cells

The development of biological markers of aging has primarily focused on adult samples. Epigenetic clocks are a promising tool for measuring biological age that show impressive accuracy across most tissues and age ranges. In adults, deviations from the DNA methylation (DNAm) age prediction are correlated with several age-related phenotypes, such as mortality and frailty. In children, however, fewer such associations have been made, possibly because DNAm changes are more dynamic in pediatric populations as compared to adults. To address this gap, we aimed to develop a highly accurate, noninvasive, biological measure of age specific to pediatric samples using buccal epithelial cell DNAm. We gathered 1,721 genome-wide DNAm profiles from 11 different cohorts of typically developing individuals aged 0 to 20 y old. Elastic net penalized regression was used to select 94 CpG sites from a training dataset (n = 1,032), with performance assessed in a separate test dataset (n = 689). DNAm at these 94 CpG sites was highly predictive of age in the test cohort (median absolute error = 0.35 y). The Pediatric-Buccal-Epigenetic (PedBE) clock was characterized in additional cohorts, showcasing the accuracy in longitudinal data, the performance in nonbuccal tissues and adult age ranges, and the association with obstetric outcomes. The PedBE tool for measuring biological age in children might help in understanding the environmental and contextual factors that shape the DNA methylome during child development, and how it, in turn, might relate to child health and disease.

The biological embedding of early-life socioeconomic status and family adversity in children's genome-wide DNA methylation

Aim:

To examine variation in child DNA methylation to assess its potential as a pathway for effects of childhood social adversity on health across the life course.

Materials & methods:

In a diverse, prospective community sample of 178 kindergarten children, associations between three types of social experience and DNA methylation within buccal epithelial cells later in childhood were examined.

Results:

Family income, parental education and family psychosocial adversity each associated with increased or decreased DNA methylation (488, 354 and 102 sites, respectively) within a unique set of genomic CpG sites. Gene ontology analyses pointed to genes serving immune and developmental regulation functions.

Conclusion:

Findings provided support for DNA methylation as a biomarker linking early-life social experiences with later life health in humans.

Systematic evaluation of DNA methylation age estimation with common preprocessing methods and the Infinium MethylationEPIC BeadChip array

BACKGROUND

The capacity of technologies measuring DNA methylation (DNAm) is rapidly evolving, as are the options for applicable bioinformatics methods. The most commonly used DNAm microarray, the Illumina Infinium HumanMethylation450 (450K array), has recently been replaced by the Illumina Infinium HumanMethylationEPIC (EPIC array), nearly doubling the number of targeted CpG sites. Given that a subset of 450K CpG sites is absent on the EPIC array and that several tools for both data normalization and analyses were developed on the 450K array, it is important to assess their utility when applied to EPIC array data. One of the most commonly used 450K tools is the pan-tissue epigenetic clock, a multivariate predictor of biological age based on DNAm at 353 CpG sites. Of these CpGs, 19 are missing from the EPIC array, thus raising the question of whether EPIC data can be used to accurately estimate DNAm age. We also investigated a 71-CpG epigenetic age predictor, referred to as the Hannum method, which lacks 6 probes on the EPIC array. To evaluate these epigenetic clocks in EPIC data properly, a prior assessment of the effects of data preprocessing methods on DNAm age is also required.

METHODS

DNAm was quantified, on both the 450K and EPIC platforms, from human primary monocytes derived from 172 individuals. We calculated DNAm age from raw, and three different preprocessed data forms to assess the effects of different processing methods on the DNAm age estimate. Using an additional cohort, we also investigated DNAm age of peripheral blood mononuclear cells, bronchoalveolar lavage, and bronchial brushing samples using the EPIC array.

RESULTS

Using monocyte-derived data from subjects on both the 450K and EPIC, we found that DNAm age was highly correlated across both raw and preprocessing methods (r > 0.91). Thus, the correlation between chronological age and the DNAm age estimate is largely unaffected by platform differences and normalization methods. However, we found that the choice of normalization method and measurement platform can lead to a systematic offset in the age estimate which in turn leads to an increase in the median error. Comparing the 450K and EPIC DNAm age estimates, we observed that the median absolute difference was 1.44-3.10 years across preprocessing methods.

CONCLUSIONS

Here, we have provided evidence that the epigenetic clock is resistant to the lack of 19 CpG sites missing from the EPIC array as well as highlighted the importance of considering the technical variance of the epigenetic when interpreting group differences below the reported error. Furthermore, our study highlights the utility of epigenetic age acceleration measure, the residuals from a linear regression of DNAm age on chronological age, as the resulting values are robust with respect to normalization methods and measurement platforms.

Maternal prenatal anxiety and child COMT genotype predict working memory and symptoms of ADHD

Maternal prenatal anxiety is an important risk factor for altered child neurodevelopment but there is uncertainty concerning the biological mechanisms involved and sources of individual differences in children's responses. We sought to determine the role of functional genetic variation in COMT, which encodes catechol-O-methyltransferase, in the association between maternal prenatal anxiety and child symptoms of ADHD and working memory. We used the prospectively-designed ALSPAC cohort (n = 6,969) for our primary data analyses followed by replication analyses in the PREDO cohort (n = 425). Maternal prenatal anxiety was based on self-report measures; child symptoms of ADHD were collected from 4-15 years of age; working memory was assessed from in-person testing at age 8 years; and genetic variation in COMT at rs4680 was determined in both mothers and children. The association between maternal prenatal anxiety and child attention/hyperactivity symptoms and working memory was moderated by the child's rs4680 genotype, with stronger effects obtained for the val/val (G:G) genotype relative to val/met (A:G) (all p<0.01) and met/met (A:A) groups (all p<0.05). Similar findings were observed in the PREDO cohort where maternal prenatal anxiety interacted with child rs4680 to predict symptoms of ADHD at 3.5 years of age. The findings, from two cohorts, show a robust gene-environment interaction, which may contribute to inter-individual differences in the effects of maternal prenatal anxiety on developmental outcomes from childhood to mid-adolescence.

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.

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.

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.

Adjusting for Cell Type Composition in DNA Methylation Data Using a Regression-Based Approach

Analysis of DNA methylation in a population context has the potential to uncover novel gene and environment interactions as well as markers of health and disease. In order to find such associations it is important to control for factors which may mask or alter DNA methylation signatures. Since tissue of origin and coinciding cell type composition are major contributors to DNA methylation patterns, and can easily confound important findings, it is vital to adjust DNA methylation data for such differences across individuals. Here we describe the use of a regression method to adjust for cell type composition in DNA methylation data. We specifically discuss what information is required to adjust for cell type composition and then provide detailed instructions on how to perform cell type adjustment on high dimensional DNA methylation data. This method has been applied mainly to Illumina 450K data, but can also be adapted to pyrosequencing or genome-wide bisulfite sequencing data.

The effect of the extended (3-year) anesthesia curriculum on anesthesia subspecialty education

The purpose of this study was to determine the effect of lengthening the anesthesia residency by the American Board of Anesthesiology on the education of anesthesia subspecialists. A survey of anesthesia residency programs was conducted from 1987 to 1991. The most frequent subspecialty practice in the clinical anesthesia (CA) 3 year is cardiovascular anesthesia. However, the percentage of CA3 residents spending extended time (greater than or equal to 6 months) in subspecialty education has significantly decreased. For example, the percentage of CA3 residents spending 12 months in subspecialty education has decreased 83%. There appears to be a slight increase in the number of CA4 or PGY5 residents (fellows) electing subspecialty practice. It is concluded that the 3-year curriculum has produced a negative impact on the education of anesthesia subspecialists.