Maternal epigenetic clocks measured during pregnancy do not predict gestational age at delivery or offspring birth outcomes: a replication study in metropolitan Cebu, Philippines

Adverse birth outcomes, such as early gestational age and low birth weight, can have lasting effects on morbidity and mortality, with impacts that persist into adulthood. Identifying the maternal factors that contribute to adverse birth outcomes in the next generation is thus a priority. Epigenetic clocks, which have emerged as powerful tools for quantifying biological aging and various dimensions of physiological dysregulation, hold promise for clarifying relationships between maternal biology and infant health, including the maternal factors or states that predict birth outcomes. Nevertheless, studies exploring the relationship between maternal epigenetic age and birth outcomes remain few. Here, we attempt to replicate a series of analyses previously reported in a US-based sample, using a larger similarly aged sample (n = 296) of participants of a long-running study in the Philippines. New pregnancies were identified prospectively, dried blood spot samples were collected during the third trimester, and information was obtained on gestational age at delivery and offspring weight after birth. Genome-wide DNA methylation was assessed with the Infinium EPIC array. Using a suite of 15 epigenetic clocks, we only found one significant relationship: advanced age on the epigenetic clock trained on leptin predicted a significantly earlier gestational age at delivery (β = - 0.15, p = 0.009). Of the other 29 relationships tested predicting gestational age and offspring birth weight, none were statistically significant. In this sample of Filipino women, epigenetic clocks capturing multiple dimensions of biology and health do not predict birth outcomes in offspring.

Alpha-synuclein overexpression induces epigenomic dysregulation of glutamate signaling and locomotor pathways

Parkinson’s disease (PD) is a neurological disorder with complex interindividual etiology that is becoming increasingly prevalent worldwide. Elevated alpha-synuclein levels can increase risk of PD and may influence epigenetic regulation of PD pathways. Here, we report genome-wide DNA methylation and hydroxymethylation alterations associated with overexpression of two PD-linked alpha-synuclein variants (wild-type and A30P) in LUHMES cells differentiated to dopaminergic neurons. Alpha-synuclein altered DNA methylation at thousands of CpGs and DNA hydroxymethylation at hundreds of CpGs in both genotypes, primarily in locomotor behavior and glutamate signaling pathway genes. In some cases, epigenetic changes were associated with transcription. SMITE network analysis incorporating H3K4me1 ChIP-seq to score DNA methylation and hydroxymethylation changes across promoters, enhancers, and gene bodies confirmed epigenetic and transcriptional deregulation of glutamate signaling modules in both genotypes. Our results identify distinct and shared impacts of alpha-synuclein variants on the epigenome, and associate alpha-synuclein with the epigenetic etiology of PD.

Prenatal exposure to phthalates and peripheral blood and buccal epithelial DNA methylation in infants: An epigenome-wide association study

BACKGROUND

Prenatal exposure to phthalates has been associated with adverse health and neurodevelopmental outcomes. DNA methylation (DNAm) alterations may be a mechanism underlying these effects, but prior investigations of prenatal exposure to phthalates and neonatal DNAm profiles are limited to placental tissue and umbilical cord blood.

OBJECTIVE

Conduct an epigenome-wide association study (EWAS) of the associations between prenatal exposure to phthalates and DNAm in two accessible infant tissues, venous buffy coat blood and buccal epithelial cells (BECs).

METHODS

Participants included 152 maternal-infant pairs from the Alberta Pregnancy Outcomes and Nutrition (APrON) study. Maternal second trimester urine samples were analyzed for nine phthalate metabolites. Blood (n = 74) or BECs (n = 78) were collected from 3-month-old infants and profiled for DNAm using the Infinium HumanMethylation450 (450K) BeadChip. Robust linear regressions were used to investigate the associations between high (HMWPs) and low molecular weight phthalates (LMWPs) and change in methylation levels at variable Cytosine-phosphate-Guanine (CpG) sites in infant tissues, as well as the sensitivity of associations to potential confounders.

RESULTS

One candidate CpG in gene RNF39 reported by a previous study examining prenatal exposure to phthalates and cord blood DNAm was replicated. The EWAS identified 12 high-confidence CpGs in blood and another 12 in BECs associated with HMWPs and/or LMWPs. Prenatal exposure to bisphenol A (BPA) associated with two of the CpGs associated with HMWPs in BECs.

DISCUSSION

Prenatal exposure to phthalates was associated with DNAm variation at CpGs annotated to genes associated with endocrine hormone activity (i.e., SLCO4A1, TPO), immune pathways and DNA damage (i.e., RASGEF1B, KAZN, HLA-A, MYO18A, DIP2C, C1or109), and neurodevelopment (i.e., AMPH, NOTCH3, DNAJC5). Future studies that characterize the stability of these associations in larger samples, multiple cohorts, across tissues, and investigate the potential associations between these biomarkers and relevant health and neurodevelopmental outcomes are needed.

Integrative Multi-Omics database (iMOMdb) of Asian Pregnant Women

Asians are underrepresented across many omics databases, thereby limiting the potential of precision medicine in nearly 60% of the global population. As such, there is a pressing need for multi-omics derived quantitative trait loci (QTLs) to fill the knowledge gap of complex traits in populations of Asian ancestry. Here, we provide the first blood-based multi-omics analysis of Asian pregnant women, constituting high-resolution genotyping (N = 1079), DNA methylation (N = 915) and transcriptome profiling (N = 238). Integrative omics analysis identified 219 154 CpGs associated with cis-DNA methylation QTLs (meQTLs) and 3703 RNAs associated with cis-RNA expression QTLs (eQTLs). Ethnicity was the largest contributor of inter-individual variation across all omics datasets, with 2561 genes identified as hotspots of this variation; 395 of these hotspot genes also contained both ethnicity-specific eQTLs and meQTLs. Gene set enrichment analysis of these ethnicity QTL hotspots showed pathways involved in lipid metabolism, adaptive immune system and carbohydrate metabolism. Pathway validation by profiling the lipidome (~480 lipids) of antenatal plasma (N = 752) and placenta (N = 1042) in the same cohort showed significant lipid differences among Chinese, Malay and Indian women, validating ethnicity-QTL gene effects across different tissue types. To develop deeper insights into the complex traits and benefit future precision medicine research in Asian pregnant women, we developed iMOMdb, an open-access database.

HIV, pathology and epigenetic age acceleration in different human tissues

Epigenetic clocks based on patterns of DNA methylation have great importance in understanding aging and disease; however, there are basic questions to be resolved in their application. It remains unknown whether epigenetic age acceleration (EAA) within an individual shows strong correlation between different primary tissue sites, the extent to which tissue pathology and clinical illness correlate with EAA in the target organ, and if EAA variability across tissues differs according to sex. Considering the outsized role of age-related illness in Human Immunodeficiency Virus-1 (HIV), these questions were pursued in a sample enriched for tissue from HIV-infected individuals. We used a custom methylation array to generate DNA methylation data from 661 samples representing 11 human tissues (adipose, blood, bone marrow, heart, kidney, liver, lung, lymph node, muscle, spleen and pituitary gland) from 133 clinically characterized, deceased individuals, including 75 infected with HIV. We developed a multimorbidity index based on the clinical disease history. Epigenetic age was moderately correlated across tissues. Blood had the greatest number and degree of correlation, most notably with spleen and bone marrow. However, blood did not correlate with epigenetic age of liver. EAA in liver was weakly correlated with EAA in kidney, adipose, lung and bone marrow. Clinically, hypertension was associated with EAA in several tissues, consistent with the multiorgan impacts of this illness. HIV infection was associated with positive age acceleration in kidney and spleen. Male sex was associated with increased epigenetic acceleration in several tissues. Preliminary evidence indicates that amyotrophic lateral sclerosis is associated with positive EAA in muscle tissue. Finally, greater multimorbidity was associated with greater EAA across all tissues. Blood alone will often fail to detect EAA in other tissues. While hypertension is associated with increased EAA in several tissues, many pathologies are associated with organ-specific age acceleration.

Epigenome-wide association study and epigenetic age acceleration associated with cigarette smoking among Costa Rican adults

Smoking-associated DNA methylation (DNAm) signatures are reproducible among studies of mostly European descent, with mixed evidence if smoking accelerates epigenetic aging and its relationship to longevity. We evaluated smoking-associated DNAm signatures in the Costa Rican Study on Longevity and Healthy Aging (CRELES), including participants from the high longevity region of Nicoya. We measured genome-wide DNAm in leukocytes, tested Epigenetic Age Acceleration (EAA) from five clocks and estimates of telomere length (DNAmTL), and examined effect modification by the high longevity region. 489 participants had a mean (SD) age of 79.4 (10.8) years, and 18% were from Nicoya. Overall, 7.6% reported currently smoking, 35% were former smokers, and 57.4% never smoked. 46 CpGs and five regions (e.g. AHRR, SCARNA6/SNORD39, SNORA20, and F2RL3) were differentially methylated for current smokers. Former smokers had increased Horvath’s EAA (1.69-years; 95% CI 0.72, 2.67), Hannum’s EAA (0.77-years; 95% CI 0.01, 1.52), GrimAge (2.34-years; 95% CI1.66, 3.02), extrinsic EAA (1.27-years; 95% CI 0.34, 2.21), intrinsic EAA (1.03-years; 95% CI 0.12, 1.94) and shorter DNAmTL (− 0.04-kb; 95% CI − 0.08, − 0.01) relative to non-smokers. There was no evidence of effect modification among residents of Nicoya. Our findings recapitulate previously reported and novel smoking-associated DNAm changes in a Latino cohort.

Dichotomy in the Impact of Elevated Maternal Glucose Levels on Neonatal Epigenome

CONTEXT

Antenatal hyperglycemia is associated with increased risk of future adverse health outcomes in both mother and child. Variations in offspring's epigenome can reflect the impact and response to in utero glycemic exposure, and may have different consequences for the child.

OBJECTIVE

We examined possible differences in associations of basal glucose status and glucose handling during pregnancy with both clinical covariates and offspring cord tissue DNA methylation.

RESEARCH DESIGN AND METHODS

This study included 830 mother-offspring dyads from the Growing Up in Singapore Towards Healthy Outcomes cohort. The fetal epigenome of umbilical cord tissue was profiled using Illumina HumanMethylation450 arrays. Associations of maternal mid-pregnancy fasting (fasting plasma glucose [FPG]) and 2-hour plasma glucose (2hPG) after a 75-g oral glucose challenge with both maternal clinical phenotypes and offspring epigenome at delivery were investigated separately.

RESULTS

Maternal age, prepregnancy body mass index, and blood pressure measures were associated with both FPG and 2hPG, whereas Chinese ethnicity (P = 1.9 × 10-4), maternal height (P = 1.1 × 10-4), pregnancy weight gain (P = 2.2 × 10-3), prepregnancy alcohol consumption (P = 4.6 × 10-4), and tobacco exposure (P = 1.9 × 10-3) showed significantly opposite associations between the 2 glucose measures. Most importantly, we observed a dichotomy in the effects of these glycemic indices on the offspring epigenome. Offspring born to mothers with elevated 2hPG showed global hypomethylation. CpGs most associated with the 2 measures also reflected differences in gene ontologies and had different associations with offspring birthweight.

CONCLUSIONS

Our findings suggest that 2 traditionally used glycemic indices for diagnosing gestational diabetes may reflect distinctive pathophysiologies in pregnancy, and have differential impacts on the offspring's DNA methylome.

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.

Maternal Prenatal Anxiety and the Fetal Origins of Epigenetic Aging

BACKGROUND

The fetal origins of mental health is a well-established framework that currently lacks a robust index of the biological embedding of prenatal adversity. The Pediatric-Buccal-Epigenetic (PedBE) clock is a novel epigenetic tool that associates with aspects of the prenatal environment, but additional validation in longitudinal datasets is required. Likewise, the relationship between prenatal maternal mental health and the PedBE clock has not been described.

METHODS

Longitudinal cohorts from the Netherlands (Basal Influences on Baby Development [BIBO] n = 165) and Singapore (Growing Up in Singapore Towards Healthy Outcomes [GUSTO] n = 340) provided data on prenatal maternal anxiety and longitudinal assessments of buccal cell-derived genome-wide DNA methylation assessed at 6 and 10 years of age in BIBO, and at 3, 9, and 48 months of age in GUSTO. Measures of epigenetic age acceleration were calculated using the PedBE clock and benchmarked against an established multi-tissue epigenetic predictor.

RESULTS

Prenatal maternal anxiety predicted child PedBE epigenetic age acceleration in both cohorts, with effects largely restricted to males and not females. These results were independent of obstetric, socioeconomic, and genetic risk factors, with a larger effect size for prenatal anxiety than depression. PedBE age acceleration predicted increased externalizing symptoms in males from mid- to late childhood in the BIBO cohort only.

CONCLUSIONS

These findings point to the fetal origins of epigenetic age acceleration and reveal an increased sensitivity in males. Convergent evidence underscores the societal importance of providing timely and effective mental health support to pregnant individuals, which may have lasting consequences for both mother and child.

Infant DNA methylation: an early indicator of intergenerational trauma?

Exposure to adverse childhood experiences (ACEs) increases risk for mental and physical health problems. Intergenerationally, mothers' ACEs predict children's health problems including neurodevelopmental and behavioural problems and poorer physical health. Theories of intergenerational trauma suggest that ACEs experienced in one generation negatively affect the health and well-being of future generations, with DNA methylation (DNAm) being one of several potential biological explanations. To begin exploring this hypothesis, we tested whether infant DNA methylation associated with intergenerational trauma. Secondary analysis employed data from the Alberta Pregnancy Outcomes and Nutrition (APrON) study. Subsample data were collected from mothers during pregnancy and postpartum on measures of distress, stress and ACEs and from infants at 3 months of age on DNAm from blood (n = 92) and buccal epithelial cells (BECs; n = 124; primarily nonoverlapping individuals between tissues). Blood and BECs were examined in separate analyses. Preliminary associations identified in blood and BECs suggest that infant DNAm patterns may relate to maternal ACEs. For the majority of ACE-related DNAm sites, neither maternal perinatal distress, nor maternal cortisol awakening response (CAR; a measure of hypothalamic-pituitary-adrenocortical axis function), substantially reduced associations between maternal ACEs and infant DNAm. However, accounting for maternal perinatal distress and cortisol substantially changed the effect of ACEs in a greater proportion of blood DNAm sites than BEC DNAm sites in the top ACEs-associated correlated methylated regions (CMRs), as well as across all CMRs and all remaining CpGs (that did not fall into CMRs). Possible DNAm patterns in infants, thus, might capture a signature of maternal intergenerational trauma, and this effect appears to be more dependent on maternal perinatal distress and CAR in blood relative to BECs.

Peroxisome proliferator-activated receptor gamma gene variants modify human airway and systemic responses to indoor dibutyl phthalate exposure

BACKGROUND

Single nucleotide polymorphisms (SNPs) of peroxisome proliferator-activated receptor gamma (PPAR-γ; gene: PPARG) and oxidative stress genes are associated with asthma risk. However, whether such variants modulate responses to dibutyl phthalate (DBP), a common plasticizer associated with increased asthma development, remains unknown. The purpose of this study is to investigate how SNPs in PPARG and oxidative stress genes, as represented by two separate genetic risk scores, modify the impact of DBP exposure on lung function and the airway and systemic response after an inhaled allergen challenge.

METHODS

We conducted a double-blinded human crossover study with sixteen allergen-sensitized participants exposed for three hours to DBP and control air on distinct occasions separated by a 4-week washout. Each exposure was followed by an allergen inhalation challenge; subsequently, lung function was measured, and blood and bronchoalveolar lavage (BAL) were collected and analyzed for cell counts and allergen-specific immunoglobulin E (IgE). Genetic risk scores for PPAR-γ (P-GRS; weighted sum of PPARG SNPs rs10865710, rs709158, and rs3856806) and oxidative stress (OS-GRS; unweighted sum of 16 SNPs across multiple genes) were developed, and their ability to modify DBP effects were assessed using linear mixed-effects models.

RESULTS

P-GRS and OS-GRS modified DBP effects on allergen-specific IgE in blood at 20 h (interaction effect [95% CI]: 1.43 [1.13 to 1.80], p = 0.005) and 3 h (0.99 [0.98 to 1], p = 0.03), respectively. P-GRS also modified DBP effects on Th2 cells in blood at 3 h (- 25.2 [- 47.7 to - 2.70], p = 0.03) and 20 h (- 39.1 [- 57.9 to - 20.3], p = 0.0005), and Th2 cells in BAL at 24 h (- 4.99 [- 8.97 to - 1.01], p = 0.02). An increasing P-GRS associated with reduced DBP effect on Th2 cells. Neither GRS significantly modified DBP effects on lung function parameters.

CONCLUSIONS

PPAR-γ variants modulated several airway and systemic immune responses to the ubiquitous chemical plasticizer DBP. Our results suggest that PPAR-γ variants may play a greater role than those in oxidative stress-related genes in airway allergic responses to DBP.

TRIAL REGISTRATION

This study reports results from The Phthalate-Allergen Immune Response Study that was registered on ClinicalTrials.gov with identification NCT02688478.

Neural correlates of attentional control in social anxiety disorder: the impact of early-life adversity and DNA methylation

BACKGROUND

Social anxiety disorder is characterized by intense fear and avoidance of social interactions and scrutiny by others. Although alterations in attentional control seem to play a central role in the psychopathology of social anxiety disorder, the neural underpinnings in prefrontal brain regions have not yet been fully clarified.

METHODS

The present study used functional MRI in participants (age 18-50 yr) with social anxiety disorder (n = 42, 31 female) and without (n = 58, 33 female). It investigated the interrelation of the effects of social anxiety disorder and early-life adversity (a main environmental risk factor of social anxiety disorder) on brain activity during an attentional control task. We applied DNA methylation analysis to determine whether epigenetic modulation in the gene encoding the glucocorticoid receptor, NR3C1, might play a mediating role in this process.

RESULTS

We identified 2 brain regions in the left and medial prefrontal cortex that exhibited an interaction effect of social anxiety disorder and early-life adversity. In participants with low levels of early-life adversity, neural activity in response to disorder-related stimuli was increased in association with social anxiety disorder. In participants with high levels of early-life adversity, neural activity was increased only in participants without social anxiety disorder. NR3C1 DNA methylation partly mediated the effect of social anxiety disorder on brain activity as a function of early-life adversity.

LIMITATIONS

The absence of behavioural correlates associated with social anxiety disorder limited functional interpretation of the results.

CONCLUSION

These findings demonstrate that the neurobiological processes that underlie social anxiety disorder might be fundamentally different depending on experiences of early-life adversity. Long-lasting effects of early-life adversity might be encoded in NR3C1 DNA methylation and entail alterations in social anxiety disorder-related activity patterns in the neural network of attentional control.

Risk-focused differences in molecular processes implicated in SARS-CoV-2 infection: corollaries in DNA methylation and gene expression

BACKGROUND

Understanding the molecular basis of susceptibility factors to the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infection is a global health imperative. It is well-established that males are more likely to acquire SARS-CoV-2 infection and exhibit more severe outcomes. Similarly, exposure to air pollutants and pre-existing respiratory chronic conditions, such as asthma and chronic obstructive respiratory disease (COPD) confer an increased risk to coronavirus disease 2019 (COVID-19).

METHODS

We investigated molecular patterns associated with risk factors in 398 candidate genes relevant to COVID-19 biology. To accomplish this, we downloaded DNA methylation and gene expression data sets from publicly available repositories (GEO and GTEx Portal) and utilized data from an empirical controlled human exposure study conducted by our team.

RESULTS

First, we observed sex-biased DNA methylation patterns in autosomal immune genes, such as NLRP2, TLE1, GPX1, and ARRB2 (FDR < 0.05, magnitude of DNA methylation difference Δβ > 0.05). Second, our analysis on the X-linked genes identified sex associated DNA methylation profiles in genes, such as ACE2, CA5B, and HS6ST2 (FDR < 0.05, Δβ > 0.05). These associations were observed across multiple respiratory tissues (lung, nasal epithelia, airway epithelia, and bronchoalveolar lavage) and in whole blood. Some of these genes, such as NLRP2 and CA5B, also exhibited sex-biased gene expression patterns. In addition, we found differential DNA methylation patterns by COVID-19 status for genes, such as NLRP2 and ACE2 in an exploratory analysis of an empirical data set reporting on human COVID-9 infections. Third, we identified modest DNA methylation changes in CpGs associated with PRIM2 and TATDN1 (FDR < 0.1, Δβ > 0.05) in response to particle-depleted diesel exhaust in bronchoalveolar lavage. Finally, we captured a DNA methylation signature associated with COPD diagnosis in a gene involved in nicotine dependence (COMT) (FDR < 0.1, Δβ > 0.05).

CONCLUSION

Our findings on sex differences might be of clinical relevance given that they revealed molecular associations of sex-biased differences in COVID-19. Specifically, our results hinted at a potentially exaggerated immune response in males linked to autosomal genes, such as NLRP2. In contrast, our findings at X-linked loci such as ACE2 suggested a potentially distinct DNA methylation pattern in females that may interact with its mRNA expression and inactivation status. We also found tissue-specific DNA methylation differences in response to particulate exposure potentially capturing a nitrogen dioxide (NO2) effect-a contributor to COVID-19 susceptibility. While we identified a molecular signature associated with COPD, all COPD-affected individuals were smokers, which may either reflect an association with the disease, smoking, or may highlight a compounded effect of these two risk factors in COVID-19. Overall, our findings point towards a molecular basis of variation in susceptibility factors that may partly explain disparities in the risk for SARS-CoV-2 infection.

Machine Learning-Based Single Cell and Integrative Analysis Reveals That Baseline mDC Predisposition Correlates With Hepatitis B Vaccine Antibody Response

Vaccination to prevent infectious disease is one of the most successful public health interventions ever developed. And yet, variability in individual vaccine effectiveness suggests that a better mechanistic understanding of vaccine-induced immune responses could improve vaccine design and efficacy. We have previously shown that protective antibody levels could be elicited in a subset of recipients with only a single dose of the hepatitis B virus (HBV) vaccine and that a wide range of antibody levels were elicited after three doses. The immune mechanisms responsible for this vaccine response variability is unclear. Using single cell RNA sequencing of sorted innate immune cell subsets, we identified two distinct myeloid dendritic cell subsets (NDRG1-expressing mDC2 and CDKN1C-expressing mDC4), the ratio of which at baseline (pre-vaccination) correlated with the immune response to a single dose of HBV vaccine. Our results suggest that the participants in our vaccine study were in one of two different dendritic cell dispositional states at baseline – an NDRG2-mDC2 state in which the vaccine elicited an antibody response after a single immunization or a CDKN1C-mDC4 state in which the vaccine required two or three doses for induction of antibody responses. To explore this correlation further, genes expressed in these mDC subsets were used for feature selection prior to the construction of predictive models using supervised canonical correlation machine learning. The resulting models showed an improved correlation with serum antibody titers in response to full vaccination. Taken together, these results suggest that the propensity of circulating dendritic cells toward either activation or suppression, their “dispositional endotype” at pre-vaccination baseline, could dictate response to vaccination.

Prenatal Adversity Alters the Epigenetic Profile of the Prefrontal Cortex: Sexually Dimorphic Effects of Prenatal Alcohol Exposure and Food-Related Stress

Prenatal adversity or stress can have long-term consequences on developmental trajectories and health outcomes. Although the biological mechanisms underlying these effects are poorly understood, epigenetic modifications, such as DNA methylation, have the potential to link early-life environments to alterations in physiological systems, with long-term functional implications. We investigated the consequences of two prenatal insults, prenatal alcohol exposure (PAE) and food-related stress, on DNA methylation profiles of the rat brain during early development. As these insults can have sex-specific effects on biological outcomes, we analyzed epigenome-wide DNA methylation patterns in prefrontal cortex, a key brain region involved in cognition, executive function, and behavior, of both males and females. We found sex-dependent and sex-concordant influences of these insults on epigenetic patterns. These alterations occurred in genes and pathways related to brain development and immune function, suggesting that PAE and food-related stress may reprogram neurobiological/physiological systems partly through central epigenetic changes, and may do so in a sex-dependent manner. Such epigenetic changes may reflect the sex-specific effects of prenatal insults on long-term functional and health outcomes and have important implications for understanding possible mechanisms underlying fetal alcohol spectrum disorder and other neurodevelopmental disorders.

Prefrontal cortex VAMP1 gene network moderates the effect of the early environment on cognitive flexibility in children

During development, genetic and environmental factors interact to modify specific phenotypes. Both in humans and in animal models, early adversities influence cognitive flexibility, an important brain function related to behavioral adaptation to variations in the environment. Abnormalities in cognitive functions are related to changes in synaptic connectivity in the prefrontal cortex (PFC), and altered levels of synaptic proteins. We investigated if individual variations in the expression of a network of genes co-expressed with the synaptic protein VAMP1 in the prefrontal cortex moderate the effect of early environmental quality on the performance of children in cognitive flexibility tasks. Genes overexpressed in early childhood and co-expressed with the VAMP1 gene in the PFC were selected for study. SNPs from these genes (post-clumping) were compiled in an expression-based polygenic score (PFC-ePRS-VAMP1). We evaluated cognitive performance of the 4 years-old children in two cohorts using similar cognitive flexibility tasks. In the first cohort (MAVAN) we utilized two CANTAB tasks: (a) the Intra-/Extra-dimensional Set Shift (IED) task, and (b) the Spatial Working Memory (SWM) task. In the second cohort, GUSTO, we used the Dimensional Change Card Sort (DCCS) task. The results show that in 4 years-old children, the PFC-ePRS-VAMP1 network moderates responsiveness to the effects of early adversities on the performance in attentional flexibility tests. The same result was observed for a spatial working memory task. Compared to attentional flexibility, reversal learning showed opposite effects of the environment, as moderated by the ePRS. A parallel ICA analysis was performed to identify relationships between whole-brain voxel based gray matter density and SNPs that comprise the PFC-ePRS-VAMP1. The early environment predicts differences in gray matter content in regions such as prefrontal and temporal cortices, significantly associated with a genetic component related to Wnt signaling pathways. Our data suggest that a network of genes co-expressed with VAMP1 in the PFC moderates the influence of early environment on cognitive function in children.

Autosomal sex-associated co-methylated regions predict biological sex from DNA methylation

Sex is a modulator of health that has been historically overlooked in biomedical research. Recognizing this knowledge gap, funding agencies now mandate the inclusion of sex as a biological variable with the goal of stimulating efforts to illuminate the molecular underpinnings of sex biases in health and disease. DNA methylation (DNAm) is a strong molecular candidate for mediating such sex biases; however, a robust and well characterized annotation of sex differences in DNAm is yet to emerge. Beginning with a large (n = 3795) dataset of DNAm profiles from normative adult whole blood samples, we identified, validated and characterized autosomal sex-associated co-methylated genomic regions (sCMRs). Strikingly, sCMRs showed consistent sex differences in DNAm over the life course and a subset were also consistent across cell, tissue and cancer types. sCMRs included sites with known sex differences in DNAm and links to health conditions with sex biased effects. The robustness of sCMRs enabled the generation of an autosomal DNAm-based predictor of sex with 96% accuracy. Testing this tool on blood DNAm profiles from individuals with sex chromosome aneuploidies (Klinefelter [47,XXY], Turner [45,X] and 47,XXX syndrome) revealed an intimate relationship between sex chromosomes and sex-biased autosomal DNAm.

Paternal adverse childhood experiences: Associations with infant DNA methylation

Adverse childhood experiences (ACEs), or cumulative childhood stress exposures, such as abuse, neglect, and household dysfunction, predict later health problems in both the exposed individuals and their offspring. One potential explanation suggests exposure to early adversity predicts epigenetic modification, especially DNA methylation (DNAm), linked to later health. Stress experienced preconception by mothers may associate with DNAm in the next generation. We hypothesized that fathers' exposure to ACEs also associates with their offspring DNAm, which, to our knowledge, has not been previously explored. An epigenome-wide association study (EWAS) of blood DNAm (n = 45) from 3-month-old infants was regressed onto fathers' retrospective ACEs at multiple Cytosine-phosphate-Guanosine (CpG) sites to discover associations. This accounted for infants' sex, age, ethnicity, cell type proportion, and genetic variability. Higher ACE scores associated with methylation values at eight CpGs. Post-hoc analysis found no contribution of paternal education, income, marital status, and parental postpartum depression, but did with paternal smoking and BMI along with infant sleep latency. These same CpGs also contributed to the association between paternal ACEs and offspring attention problems at 3 years. Collectively, these findings suggested there were biological associations with paternal early life adversity and offspring DNAm in infancy, potentially affecting offspring later childhood outcomes.

Optimized CRISPR-mediated gene knockin reveals FOXP3-independent maintenance of human Treg identity

Regulatory T cell (Treg) therapy is a promising curative approach for a variety of immune-mediated conditions. CRISPR-based genome editing allows precise insertion of transgenes through homology-directed repair, but its use in human Tregs has been limited. We report an optimized protocol for CRISPR-mediated gene knockin in human Tregs with high-yield expansion. To establish a benchmark of human Treg dysfunction, we target the master transcription factor FOXP3 in naive and memory Tregs. Although FOXP3-ablated Tregs upregulate cytokine expression, effects on suppressive capacity in vitro manifest slowly and primarily in memory Tregs. Moreover, FOXP3-ablated Tregs retain their characteristic protein, transcriptional, and DNA methylation profile. Instead, FOXP3 maintains DNA methylation at regions enriched for AP-1 binding sites. Thus, although FOXP3 is important for human Treg development, it has a limited role in maintaining mature Treg identity. Optimized gene knockin with human Tregs will enable mechanistic studies and the development of tailored, next-generation Treg cell therapies.

DNA methylation differences associated with social anxiety disorder and early life adversity

Social anxiety disorder (SAD) is a psychiatric disorder characterized by extensive fear in social situations. Multiple genetic and environmental factors are known to contribute to its pathogenesis. One of the main environmental risk factors is early life adversity (ELA). Evidence is emerging that epigenetic mechanisms such as DNA methylation might play an important role in the biological mechanisms underlying SAD and ELA. To investigate the relationship between ELA, DNA methylation, and SAD, we performed an epigenome-wide association study for SAD and ELA examining DNA from whole blood of a cohort of 143 individuals using DNA methylation arrays. We identified two differentially methylated regions (DMRs) associated with SAD located within the genes SLC43A2 and TNXB. As this was the first epigenome-wide association study for SAD, it is worth noting that both genes have previously been associated with panic disorder. Further, we identified two DMRs associated with ELA within the SLC17A3 promoter region and the SIAH3 gene and several DMRs that were associated with the interaction of SAD and ELA. Of these, the regions within C2CD2L and MRPL28 showed the largest difference in DNA methylation. Lastly, we found that two DMRs were associated with both the severity of social anxiety and ELA, however, neither of them was found to mediate the contribution of ELA to SAD later in life. Future studies are needed to replicate our findings in independent cohorts and to investigate the biological pathways underlying these effects.

Internalizing symptoms associate with the pace of epigenetic aging in childhood

Childhood psychiatric symptoms may be associated with advanced biological aging. This study examined whether epigenetic age acceleration (EAA) associates with internalizing and externalizing symptoms that were prospectively collected across childhood in a longitudinal cohort study. At age 6 buccal epithelial cells from 148 children (69 girls) were collected to survey genome-wide DNA methylation. EAA was estimated using the Horvath clock. Internalizing symptoms at ages 2.5 and 4 years significantly predicted higher EAA at age 6, which in turn was significantly associated with internalizing symptoms at ages 6-10 years. Similar trends for externalizing symptoms did not reach statistical significance. These findings indicate advanced biological aging in relation to child mental health and may help better identify those at risk for lasting impairments associated with internalizing disorders.

Human T-bet Governs Innate and Innate-like Adaptive IFN-γ Immunity against Mycobacteria

Inborn errors of human interferon gamma (IFN-γ) immunity underlie mycobacterial disease. We report a patient with mycobacterial disease due to inherited deficiency of the transcription factor T-bet. The patient has extremely low counts of circulating Mycobacterium-reactive natural killer (NK), invariant NKT (iNKT), mucosal-associated invariant T (MAIT), and Vδ2+ γδ T lymphocytes, and of Mycobacterium-non reactive classic TH1 lymphocytes, with the residual populations of these cells also producing abnormally small amounts of IFN-γ. Other lymphocyte subsets develop normally but produce low levels of IFN-γ, with the exception of CD8+ αβ T and non-classic CD4+ αβ TH1∗ lymphocytes, which produce IFN-γ normally in response to mycobacterial antigens. Human T-bet deficiency thus underlies mycobacterial disease by preventing the development of innate (NK) and innate-like adaptive lymphocytes (iNKT, MAIT, and Vδ2+ γδ T cells) and IFN-γ production by them, with mycobacterium-specific, IFN-γ-producing, purely adaptive CD8+ αβ T, and CD4+ αβ TH1∗ cells unable to compensate for this deficit.

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.

Occurrence of Accelerated Epigenetic Aging and Methylation Disruptions in Human Immunodeficiency Virus Infection Before Antiretroviral Therapy

BACKGROUND

Whether accelerated aging develops over the course of chronic human immunodeficiency virus (HIV) infection or can be observed before significant immunosuppression on is unknown. We studied DNA methylation in blood to estimate cellular aging in persons living with HIV (PLWH) before the initiation of antiretroviral therapy (ART).

METHODS

A total of 378 ART-naive PLWH who had CD4 T-cell counts >500/µL and were enrolled in the Strategic Timing of Antiretroviral Therapy trial (Pulmonary Substudy) were compared with 34 HIV-negative controls. DNA methylation was performed using the Illumina MethylationEPIC BeadChip. Differentially methylated positions (DMPs) and differentially methylated regions (DMRs) in PLWH compared with controls were identified using a robust linear model. Methylation age was calculated using a previously described epigenetic clock.

RESULTS

There were a total of 56 639 DMPs and 6103 DMRs at a false discovery rate of <0.1. The top 5 DMPs corresponded to genes NLRC5, VRK2, B2M, and GPR6 and were highly enriched for cancer-related pathways. PLWH had significantly higher methylation age than HIV-negative controls (P = .001), with black race, low CD4 and high CD8 T-cell counts, and duration of HIV being risk factors for age acceleration.

CONCLUSIONS

PLWH before the initiation of ART and with preserved immune status show evidence of advanced methylation aging.