Genome-wide neonatal epigenetic changes associated with maternal exposure to the COVID-19 pandemic
BMC Medical Genomics volume 16, Article number: 268 (2023)
During gestation, stressors to the fetus, including viral exposure or maternal psychological distress, can fundamentally alter the neonatal epigenome, and may be associated with long-term impaired developmental outcomes. The impact of in utero exposure to the COVID-19 pandemic on the newborn epigenome has yet to be described.
This study aimed to determine whether there are unique epigenetic signatures in newborns who experienced otherwise healthy pregnancies that occurred during the COVID-19 pandemic (Project RESCUE). The pre-pandemic control and pandemic cohorts (Project RESCUE) included in this study are part of a prospective observational and longitudinal cohort study that evaluates the impact of elevated prenatal maternal stress during the COVID-19 pandemic on early childhood neurodevelopment.
Using buccal swabs collected at birth, differential DNA methylation analysis was performed using the Infinium MethylationEPIC arrays and linear regression analysis. Pathway analysis and gene ontology enrichment were performed on resultant gene lists.
Widespread differential methylation was found between neonates exposed in utero to the pandemic and pre-pandemic neonates. In contrast, there were no apparent epigenetic differences associated with maternal COVID-19 infection during pregnancy. Differential methylation was observed among genomic sites that underpin important neurological pathways that have been previously reported in the literature to be differentially methylated because of prenatal stress, such as NR3C1.
The present study reveals potential associations between exposure to the COVID-19 pandemic during pregnancy and subsequent changes in the newborn epigenome. While this finding warrants further investigation, it is a point that should be considered in any study assessing newborn DNA methylation studies obtained during this period, even in otherwise healthy pregnancies.
We and others have documented that maternal psychological distress, such as stress, anxiety, and/or depression, during the prenatal period has been strongly correlated with impaired fetal development and long-term neurodevelopmental and cognitive consequences throughout the child’s life [1,2,3,4,5,6,7,8,9]. Large-scale traumatic events, such as natural disasters, war, genocide, famine, and pandemics, have been studied for their acute impact on pregnancy, maternal mental health, and subsequent neonatal outcomes [1, 7, 8]. The Dutch Hunger Winter (1944–45) and Canadian Ice Storm (1998) are exemplar traumatic events that have been studied for the postnatal consequences associated with maternal distress and trauma [10,11,12,13,14,15,16]. These studies have leveraged DNA methylation analysis to uncover mechanisms associated with environmental stress, the impact on maternal mental health, and how differential methylation of the fetal genome may affect developmental outcomes [17,18,19,20]. Numerous studies of maternal perinatal and early life stress suggest that there are critical changes to the entire newborn epigenome that may underlie lifelong developmental and neurological outcomes in offspring [3, 7,8,9,10, 21,22,23]. The present study builds on what is known about epigenetic changes detected in cord blood, by assessing epigenetic associations with exposure to an unprecedented global pandemic utilizing buccal swabs collected from neonates in the first month of life [22, 23]. Notably, Sammallahti et al. have highlighted that cord blood shows inconclusive evidence of the association between DNA methylation changes associated with maternal anxiety . It is thus of the utmost importance to expand the study perinatal maternal distress and other exposures on DNA methylation to multiple tissue types and timepoints, such as newborn buccal swabs, as the epigenome is known to vary considerably as a function of time and cell type. The nexus of perinatal maternal stress and epigenomic research is critical to providing a foundation for therapeutic intervention during pregnancy and beyond to improve both maternal and fetal outcomes.
The sudden onset of the COVID-19 pandemic has provided multiple causes of both acute and chronic stress with recession, isolation, reduced access to healthcare, and fear of, or exposure to, the COVID-19 virus. Since March 2020 when the COVID-19 pandemic was deemed a “natural disaster” in the United States by FEMA, the global health burden of disease has climbed to over 400 million confirmed cases and over 6 million deaths, with new variants of the virus emerging frequently . As a result, the COVID-19 pandemic has been correlated with a trend of increased stress, anxiety, and depression among the global population , with pregnant women being no exception from the trend [26, 27].
To document the impact of maternal exposure during the COVID-19 pandemic on fetal development, we began recruitment of pregnant women into the RESCUE Project (Reducing Elevated Stress from COVID-19 Exposure Project) research study in June 2020. For the scope of this study, we sought to characterize the neonatal epigenetic landscape in healthy pregnancies delivered before and after the onset of COVID-19, using pre-pandemic newborn controls (CTL cohort) as a direct comparison to Project RESCUE newborns (RES cohort) recruited from the greater Washington DC Metropolitan region. Due to the varying policies during the COVID-19 pandemic along with varying periods of related surges, we sought to reduce potential confounding factors by focusing on a small, controlled regional cohort with participants largely exposed to the same environment and recruitment schema before and during the COVID-19 pandemic . This study provided a unique opportunity to uncover potential signatures of distress during pregnancy on infant DNA methylation and whether these may predict long-term resilience or susceptibility to adverse neurodevelopmental and behavioral outcomes.
Material and methods
Demographic and clinical data were extracted from maternal questionnaires and medical records. Women eligible for this study were recruited from the Washington DC-Metropolitan area (Virginia, Washington DC, Maryland) and met the following criteria: to be over the age of 18 years with a singleton pregnancy of 8 weeks gestation or greater and to be diagnosed with or are exhibiting symptoms of stress, anxiety, or depression, or have no other health risk factors outside of being pregnant during the COVID-19 pandemic. Additionally, newborn infants born to mothers who tested COVID-19 positive during pregnancy were enrolled and participate in all post-natal portions of the study. Exclusionary criteria for all cohorts included: fetuses/newborns with known chromosomal syndromic conditions, women unable to enter the MRI scanner for physical or psychological reasons, or women who have health conditions that make their pregnancy high-risk. Maternal exposure to and fetal gestational age (GA) at time of COVID-19 infection was captured from medical records and test results. The Coronavirus Perinatal Experiences- Impact Survey (COPE-IS) and the COPE Impact Update (COPE-IU) were also completed to measure the experiences of new and expectant mothers in the time of the Coronavirus COVID-19 (SARS-CoV-2) pandemic . Infant gestational age (GA), sex, birthweight, head circumference, and length at birth, as well as maternal age were recorded at the time of visit. Race and ethnicity information was self-reported. These metrics were assessed for statistical differences using a paired t-test and p-value ≤ 0.05, as well as z-score, when applicable (Table 1). Maternal mental health assessments included several self-reported surveys; The State-Trait Anxiety Inventory was used to document acute (STAI-S) and chronic (STAI-T) anxiety [30, 31]; Perceived Stress Scale (PSS) to measure stress [32, 33]; and the Impact of Events Scale-Revised  to assess for post-traumatic stress disorder (PTSD) symptoms. Scoring was performed and aggregate data was analyzed for “peak,” or highest possible scoring metric, throughout pregnancy and at birth and that “peak” score was used to label the sample with "low", “moderate”, or “high” STAI-S, STAI-T, and PSS (or “yes” or “no” PTSD). (The specific calculations for all maternal mental health surveys can be found in Supplemental Table S1).
Neonatal buccal DNA extraction
Buccal swabs from the RESCUE cohort were collected for global DNA methylation analysis between June 2020 and February 2021 (i.e., prior to deployment of COVID-19 vaccines to the general population). All CTL swabs were collected prior to December 2019. ORAcollect buccal swabs (OG-175, DNA Genotek) were collected from neonates (average day of life 4.72 weeks ± 2.4), just prior to feeds and stored at 4˚C for batch processing. DNA was extracted with the PrepIT-L2P kit (DNA Genotek) and DNA was quantified using the Qubit Broad Range dsDNA assay kit and Qubit 4 fluorometer (Thermofisher Scientific). DNA was stored at -20˚C until ready for downstream applications.
Genome-wide DNA methylation array
An input of 300 ng of DNA was bisulfite-converted using the DNA Methylation-Lightning kit (Zymo Research). After whole-genome amplification and enzymatic fragmentation, samples were hybridized to BeadChip arrays using the Infinium MethylationEPIC BeadChip kit according to the manufacturer’s protocol (Illumina). Intensity values at the over 850,000 methylation sites on the BeadChips were measured across the genome at single-nucleotide resolution using iScan (Illumina) or Nextseq 550 platform with BeadChip adapter (Illumina). Of note, all laboratory procedures from DNA extraction to DNA methylation array preparation and scanning were completed by one individual, in order to reduce variability between batches.
Differential methylation analysis
We adopted the following pipeline for all further studies (Supplemental Fig. S1A). Quality filtration, functional normalization, and differential methylation analysis were performed on all arrays [35, 36]. Quality filtration was performed using wateRmelon, which removes probes if bead count of less than 3 in ≥ 5% of the samples . Functional normalization of the data was executed using the minfi package with ssnoob background and dye correction . Additional filtration using minfi excluded all probes that had a signal not significantly above background control probes (p-value > 0.05), all probes representing single nucleotide variants present in the general population with a minor allelic frequency > 5% (dropLociWithSnps), and cross-reactive probes (maxprobes) [39, 40].
DNA methylation array pipeline
Because CTL samples were originally processed on the iScan platform, we included several replicates of CTL samples with each batch analysis of RES samples on the Nextseq platform to identify potential platform interference. None of the sample rows on either system failed during the scanning process and batch correction was used during bioinformatic analysis to account for the two distinct scanning platforms (Supplemental Fig. S1). To validate the analytic platform, we included all samples processed on these platforms in the same time-frame as those reported in this study; these include the 44 RES and CTL samples reported here, plus another 58 samples from other unrelated stress cohorts.
As seen in Supplemental Fig. S1, prior to any processing (865,859 probes), hierarchical clustering showed clustering of samples by platform and replicates did not cluster together (asterisks in Fig. S1B). After functional normalization, density plots of raw data (777,755 probes) were greatly improved (compare Supplemental Fig. S1C, E) but duplicates did not cluster together (Supplemental Fig. S1D). After batch correction, density plot is further improved (Supplemental Fig. S1G) and replicates cluster together (Supplemental Fig. S1F). We tested an alternative method for normalization, SWAN , but after SWAN and batch correction, while the replicates clustered (Supplemental Fig. S1H), the density plots appeared less ordered (Supplemental Fig. S1I).
Beta-values were extracted, followed by batch correction and differential methylation analysis through linear regression models via limma. Replicate samples were run on each platform and in each batch to remove potential batch effects computationally. Our group previously reported significance thresholding using a commonly reported standard in the literature: results were considered statistically significant for differentially methylated probes (DMPs) with a false discovery rate (FDR) adjusted p-value < 0.05 and log2 fold-change ≥ 1 or ≤ -1 [40,41,42,43]. The networks and functional analyses were generated with QIAGEN IPA and Gene Ontology Enrichment Analysis GO with Enrichr [44, 45]; both tools used a statistical significance threshold of p ≤ 0.05. For pathway analysis (QIAGEN IPA), z-scores were calculated based off mean methylation differences at the multiple CpG sites (Illumina IDs) associated with annotated genes or pathways between CTL and RES cohorts (log2 fold-change). The z-score was calculated through the IPA software by assessing the overall log2 fold-change of all CpG sites associated with a specific pathway (i.e., predominantly negative fold-change differences are associated with hypomethylation in the RES compared to CTL and a negative z-score, where a majority positive fold-change indicates a positive z-score and hypermethylation in RES compared to control). A z-score of 0 means that the sites associated with a particular pathway included a mix of hyper- and hypomethylated sites (positive and negative fold-changes) that netted an overall score of zero. Graphs were made with ChAMP and ggplot2 in R [46, 47]. Comparisons of gene lists were performed on annotated probes only using Venny 2.1 .
We performed genome-wide DNA methylation analysis in two newborn cohorts, one born before the onset of the COVID-19 pandemic (CTL) and one recruited starting June 2020 (RES). Samples used in this study were age- and sex-matched (Table 1). Maternal age at delivery and height, weight, and length at delivery were not significantly different between the two cohorts (Table 1). Both cohorts were ethnically diverse and had a similar average maternal age; only one participant in this study indicated use of tobacco during pregnancy.
Unsupervised hierarchical clustering (Fig. 1a) of DMPs showed the pre-pandemic CTL cohort (in yellow) clustering away from the samples of the RES cohort (in pink). This was also visualized in principal component analysis (Fig. 1b) with clear separation and clustering of pre-pandemic newborns (circles) and newborns exposed to the COVID-19 pandemic during pregnancy (triangles). Linear regression analysis of CTL (n = 12) and RES cohort (n = 33) resulted in significant (FDR p ≤ 0.05) differential DNA methylation at 519 annotated sites (675 total probes) throughout the genome (Fig. 1c, Supplemental Table S2). Among the differentially methylated CpG sites, hypermethylation of NR3C1, a well-known epigenetic marker of prenatal stress, was observed in the RES cohort. Pathway analysis of these 519 differentially methylated sites suggests that increased methylation (positive z-score) may be associated with synaptogenesis signaling, PFKB4 signaling, Gs alpha signaling, and PPARɑ/RXRɑ activation, while decreased methylation (negative z-score) may be associated with neurodegenerative pathways (Huntington’s Disease), opioid signaling, and IL-15 production (Fig. 1d, Supplemental Table S3). Additional analysis of annotated DMPs using a distinct curated database (Supplemental Figure S3, Supplemental Table S4) confirmed association of DMPs with immune system pathways (IL-5 and IL-13 production) and nervous system developmental processes (synapse structural plasticity, blood–brain-barrier permeability).
Interestingly, while the CTL and RES cohorts showed striking differences in the methylation landscape genome-wide, the infection status of the mother during pregnancy did not seem to result in specific epigenetic signatures in the newborn (hashed vs. solid pink in Fig. 1a; green vs. orange triangles in Fig. 1b). There were no significant DMPs among the RES cohort between COVID-19 positive (n = 15) and negative (n = 17) pregnancies at the FDR p-value ≤ 0.05 and log2 fold-change ≥ 1 or ≤ -1 significance threshold, nor differences found in newborns exposed to asymptomatic, symptomatic, or hospitalized maternal COVID-19 infection (Fig. 1c green, yellow, orange, and red boxes).
Maternal mental health metrics showed that overall, the RES cohort had a higher number of participants who scored with “high” state (STAI-S) and trait (STAI-T) anxiety and “moderate" or "high” stress (PSS) when compared to CTL (Fig. 2; data in Supplemental Table S1). Anxiety during pregnancy is estimated to affect 15–23% of women [49,50,51], which is a trend observable in the CTL cohort (~ 18% for both state and trait anxiety). While not significantly different from CTL, the RES cohort notably showed a trend toward a higher rate of both state and trait anxiety (30% and 21%, respectively) and stress (86% vs. 67% in CTL). Within RES, over a quarter of the women experienced PTSD symptoms at some point during pregnancy, even higher than the available report of 15% in the general population during the COVID-19 pandemic . (The IESR questionnaire had not been used for the pre-pandemic CTL cohort). Using these additional mental health metrics (anxiety, stress, and PTSD), differential methylation analysis within the RES cohort (grouping high scores for STAI-S, STAI-T and high-moderate scores for PSS) did not identify significant methylation differences among these subcategories (Supplemental Figure S2; data in Supplemental Table S1). This can be further observed in the initial RES and CTL analysis, where there is no obvious clustering of DMPs associated with high stress, anxiety, or PTSD (Fig. 1c, fuchsia, lime, and aqua boxes).
Finally, due to the observed clustering among the RES cohort (Fig. 1a, pink) with one sub-cluster on either side of the (yellow) CTL cluster, we explored numerous metrics and found the timing of the pandemic during pregnancy to be linked with this pattern. "Early" pregnancies (gestation starting in late 2019-January 2020, before widespread knowledge of the pandemic, n = 13) and “late” pregnancies (conceived February-June 2020, n = 19) clearly segregated not only away from CTLs but also from each other (Fig. 1a). In the overall DNA methylation analysis between RES and CTL samples, there is distinct clustering of DMPs by timing (Fig. 1c, blue and navy boxes). Differential methylation analysis between RES newborns based on timing of pregnancy and CTLs identified 462 unique DMPs (632 total probes; Supplemental Table S5) associated with “early” 2020 pregnancies and 795 unique DMPs (1119 total probes; Supplemental Table S5) with “late” 2020 pregnancies, suggesting that the temporal aspect of the pandemic during pregnancy may have a unique impact on newborn DNA methylation. Further, we noted there was an overlap of 35.4% (331 annotated probes) in shared DMPs among all the resultant gene list from RES, “early” 2020, and “late” 2020 pandemic pregnancy newborns analysis against CTLs (Supplemental Figure S4). Finally, we performed a linear regression analysis between the early and late pandemic groups and found 324 sites of significant differential methylation, which further emphasizes the differences associated with timing of the pandemic on DNA methylation in these newborn cases (Supplemental Table 5).
In the present genome-wide DNA methylation study, significant differences were observed in otherwise healthy pregnancies during the pandemic when compared to pre-pandemic controls. With more than 500 annotated sites of differential methylation between these two groups, we observed significant global changes to the neonatal epigenomic landscape. Among the annotated sites of differential methylation, we noted methylation within the RES cohort at NR3C1. NR3C1 is a glucocorticoid receptor gene, which has been associated with increased methylation in other studies of perinatal stress to newborns [53, 54], as well as trauma , psychological disorders, including major depression, post-traumatic stress disorder, anxiety, and personality disorders . The findings presented here are consistent with existing literature that highlights the impact of prenatal stress on the newborn epigenome. Future studies will focus on validation of NR3C1, as well as the other differentially methylated sites presented in this study, as biomarkers of perinatal stress and aim to assess any potential associations between neonatal epigenomics with long-term neurodevelopmental outcomes.
Top pathways associated with differential methylation in the RES cohort included IL-15 production. Maternal stress during pregnancy has been strongly linked to increased inflammation and risk for future neuropsychiatric outcomes in offspring . IL-15 is a proinflammatory cytokine that is critical for healthy maternal and fetal outcomes in pregnancy , and upregulation of IL-15 during pregnancy has been correlated with adverse outcomes of pregnancy, such as preeclampsia and gestational diabetes mellitus [59, 60], and with neonatal neurodevelopmental delay . While IL-15 production is a key biological response to viral infection, this pathway was more methylated regardless of viral infection status in the mother. Its potential downregulation in newborns of the pandemic may function as a protective mechanism against negative developmental brain outcomes.
The neurodegenerative disease pathway, Huntington’s disease signaling and synaptogenesis signaling were also top terms significantly associated with differential methylation in pandemic pregnancy newborns. Impaired neurodevelopmental outcomes are a strongly correlated with fetal stress and often associated with reduced infant brain plasticity and dysregulated developmental processes, such as synaptogenesis, neuronal migration, and myelination [62, 63]. Mitochondrial dysfunction, inflammation, and oxidative stress are key players in neurodegeneration and neurodegenerative diseases, such as Huntington’s disease [64, 65], that also have been shown to be triggered by stress during pregnancy and perhaps lead to impaired neurodevelopment in newborns [66,67,68]. Together, this combination of pathways may provide insight into the mechanisms that underpin the impact of fetal stressors during pregnancy and long-term impaired neurodevelopmental outcomes that follow.
Another top hypomethylated pathway in the RES cohort was the PFKB4 pathway. The PFKB4 is both a regulator of glycolysis and a regulator of progression of ectodermal patterning toward specific fates including neural plate and neural crest during embryonic development, independently of glycolysis . It also promotes tumor growth and metastasis , and whether this will have an impact on developmental outcomes remains to be seen. It will be critical to document longitudinal developmental outcomes in both CTL and RES cohort to test these mechanistic hypotheses.
As there were no significant differences in DNA methylation associated with maternal COVID-19 infection in this analysis, it may be inferred that the changes occurring to the neonatal are epigenome are likely a direct result of the COVID-19 pandemic environment, rather than maternal COVID-19 infection status during pregnancy. This is congruent with the literature, which suggests that there is likely only rare instances of vertical intrauterine transmission of the SARS-CoV-2 virus from infected mothers to their fetuses [53, 71].
Finally, after noticing that two distinct clusters emerged among the pandemic RES cohort, we subset the DNA methylation analysis to examine whether the onset ("early" group) and continuation ("late" group) of the COVID-19 pandemic may have resulted in different and significant changes in this cohort. Here, we noted over 462 unique annotated sites of differential methylation between the “early” 2020 pandemic group and CTL and 795 annotated sites between the “late” 2020 pandemic group and CTL (Supplemental Table S5). The “early” group included women who were in their first trimester immediately before to the onset of the COVID-19 pandemic (between September 2019 and January 2019), while the late 2020 group included women who were in their first trimester between February 2020 and May 2020 and thus directly impacted by the COVID-19 pandemic for the entirety of their pregnancy. When looking at all groups combined, we noted 77 annotated sites unique to the “early” pandemic RES cohort and 346 unique to the “late” pandemic RES cohort, with an overlap of 331 annotated sites between early, late, and all RES pandemic pregnancy cohort newborns (Supplemental Figure S4). To further emphasize the differences associated with the timing of pandemic exposure, we conducted an additional differential methylation analysis between the early and late pandemic groups, which resulted in a difference of 324 sites (Supplemental Table 5). Published reports suggest that the timing of stress during pregnancy and length of exposure may play an important role in the ultimate postnatal outcomes of a child and, as such, could play a role the differences in DNA methylation within the RES cohort [72, 73]. For example, a 2018 study focused on “military stress” during pregnancy subset their analysis between “new onset” and “chronic” stress exposure, as well as the trimester that the stress occurred . While the exact nature of these differences has yet to be fully described, the present data suggest a distinct and significant impact that the timing and type of distress induced by the COVID-19 pandemic during pregnancy had on the newborn epigenetic landscape. Additionally, future studies may benefit from considering the concept of chronoepigenetics, or study of epigenetic changes over time, to better understand the impact of insults and stressors during vulnerable periods of development and how they may cause critical changes to the fetal and neonatal epigenome, as well as considering additional factors such regionally diverse datasets that may have been differentially impacted by COVID-19 pandemic policies and guidelines .
A few limitations impact this study: availability of COVID-19 testing, missing or non-compliance with completing self-reported psychosocial surveys, and cohort size. Maternal COVID-19 infection status was determined by positive testing noted in the medical records, as well self-reported testing by participants of the RES study. Participants among the RES cohort without known exposures or asymptomatic COVID-19 infections may not have been accurately captured. Also, availability of testing evolved during enrollment for this study. Maternal mental health measures between RES and CTL did not reach statistical significance, however we were likely underpowered given the sample size of the pre-pandemic cohort (n = 12). The maternal mental health surveys (STAI-S, STAI-T, PSS, and IESR) used may not fully capture the unprecedented or dynamic nature of the type of distress that occurred throughout the COVID-19 pandemic. Restriction in access to hospital settings for healthy research volunteers imposed during the pandemic limited study visits for the CTL cohort and missing data, as well as recruitment bias specific to the Washington DC-Metropolitan area, may skew the mental health data acquired in this study . Of note, while there was no difference in average gestational ages at time of collection between the RES and CTL cohorts, buccal swabs were collected within a 2–4 week period after birth (Table 1). Therefore, the data presented here may not be representative of the newborn genome immediately after birth, nor does it represent all COVID-19 pandemic or otherwise healthy pregnancies . Our group has previously reported that the overall paucity and disparity of publicly available data during the onset and continuation of the COVID-19 pandemic has been a pervasive issue and ultimately hampers medical, research, and public health policy decisions [75, 76]. The present study included a total of 44 participants, which is a small sample size for global DNA methylation array studies. To compensate for this, we employed stringent differential methylation criteria (i.e., FDR p-value < 0.05, effect size cutoff of log2 fold-change cutoff of greater than 1 or less than -1). While this improves the statistical rigor of the subsequent differential methylation analysis, we once again emphasize that a larger, well-powered sample population is preferable for a study of this nature. We also acknowledge that while stringent filtering criteria, statistical thresholds, and batch correction was employed in our bioinformatic analysis, potential batch effects, false positives, and confounding factors (time of collection, early infant environmental exposures, etc.) may still exist and impact the data presented here . The breadth of the changes, however, warranted rapid reporting even with these limitations, as the data presented in this study could be of critical importance from a research and public health perspective.
In this study, we found a striking alteration of the neonatal epigenome coincident with the COVID-19 pandemic, but not maternal COVID-19 infection, a factor that must be considered in future DNA methylation studies comparing pre-pandemic data. While there was not significant correlation between maternal mental health survey data (stress, anxiety, depression), fetuses exposed to the onset and continuation of COVID-19 pandemic throughout 2020 had over 500 sites of significant differential DNA methylation, including NR3C1 and other genes associated with pathways involved in neurodevelopment and immune system response. Of note, there also appears to be a temporal aspect of the pandemic on the neonatal epigenome with divergent signatures appearing among infants exposed to just the onset of the pandemic in early 2020 and those exposed to the continuation of the pandemic throughout mid-late 2020. Overall, the widespread differences seen in DNA methylation between pre-pandemic and pandemic-exposed newborns provides critical insight into future development of these cohorts, which should be carefully documented longitudinally.
Availability of data and materials
The datasets generated and analyzed in the current study are pending approval but will be available in the Gene Expression Omnibus (GEO) repository. Until the public release date, 2023–12-31, please use this temporary link to access the data used in this study (GSE229463).
GEO data access token: azadmymilhcvlcb.
Coronavirus Disease 2019
- RESCUE, RES:
Reducing Elevated Stress from COVID-19 Exposure
State-Trait Anxiety Inventory
Perceived Stress Scale
Impact of Event Scale Revised
Post-traumatic stress disorder
Differentially methylated probe
Wu Y, Lu YC, Jacobs M, Pradhan S, Kapse K, Zhao L, et al. Association of Prenatal Maternal Psychological Distress with Fetal Brain Growth, Metabolism, and Cortical Maturation. JAMA Netw Open. 2020;3(1):e1919940.
Abbasi H, Puusepp-Benazzouz H. Pregnancy stress, anxiety, and depression sequela on neonatal brain development. JAMA Pediatrics. 2020;174:908.
Graignic-Philippe R, Dayan J, Chokron S, Jacquet AY, Tordjman S. Effects of prenatal stress on fetal and child development: A critical literature review. Neurosci Biobehav Rev. 2014;43:137–62.
Kinsella MT, Monk C. Impact of maternal stress, depression and anxiety on fetal neurobehavioral development. Clin Obstet Gynecol. 2009;52:425–40.
Weinstock M. The long-term behavioural consequences of prenatal stress. Neurosci Biobehav Rev. 2008;32:1073–86.
Van Den Bergh BRH, Mulder EJH, Mennes M, Glover V. Antenatal maternal anxiety and stress and the neurobehavioural development of the fetus and child: Links and possible mechanisms. A review. In: Neuroscience and Biobehavioral Reviews. 2005.
Wu Y, Limperopoulos C. Pregnancy stress, anxiety, and depression sequela on neonatal brain development - Reply. JAMA Pediatrics. 2020;174:908–9.
Wu Y, Kapse K, Jacobs M, Niforatos-Andescavage N, Donofrio MT, Krishnan A, et al. Association of Maternal Psychological Distress with in Utero Brain Development in Fetuses with Congenital Heart Disease. JAMA Pediatr. 2020;174(3):908–9.
Wu Y, Espinosa KM, Barnett SD, Kapse A, Quistorff JL, Lopez C, et al. Association of Elevated Maternal Psychological Distress, Altered Fetal Brain, and Offspring Cognitive and Social-Emotional Outcomes at 18 Months. JAMA Netw Open. 2022;5(4):e229244.
Yong Ping E, Laplante DP, Elgbeili G, Jones SL, Brunet A, King S. Disaster-related prenatal maternal stress predicts HPA reactivity and psychopathology in adolescent offspring: Project Ice Storm. Psychoneuroendocrinology. 2020;117:104697.
Dashorst P, Mooren TM, Kleber RJ, de Jong PJ, Huntjens RJC. Intergenerational consequences of the Holocaust on offspring mental health: a systematic review of associated factors and mechanisms. Eur J Psychotraumatol. 2019;10:1654065.
Guerra O, Eboreime E. The impact of economic recessions on depression, anxiety, and trauma-related disorders and illness outcomes—A scoping review. Behav Sci. 2021;11:119.
King S, Dancause K, Turcotte-Tremblay AM, Veru F, Laplante DP. Using natural disasters to study the effects of prenatal maternal stress on child health and development. Birth Defects Res Part C Embryo Today. 2012;96:273–88.
Laplante DP, Brunet A, King S. The effects of maternal stress and illness during pregnancy on infant temperament: Project Ice Storm. Pediatr Res. 2016;79(1):107–13.
Roseboom TJ, Painter RC, Van Abeelen AFM, Veenendaal MVE, De Rooij SR. Hungry in the womb: What are the consequences? Lessons from the Dutch famine. Maturitas. 2011;70:141–5.
Hart N. Famine, maternal nutrition and infant mortality: a re-examination of the dutch hunger winter. Popul Stud (NY). 1993;47(1):27–46.
Cao-Lei L, Massart R, Suderman MJ, Machnes Z, Elgbeili G, Laplante DP, et al. DNA methylation signatures triggered by prenatal maternal stress exposure to a natural disaster: Project ice storm. PLoS One. 2014;9(9):e107653.
Cao-Lei L, Dancause KN, Elgbeili G, Massart R, Szyf M, Liu A, et al. DNA methylation mediates the impact of exposure to prenatal maternal stress on BMI and central adiposity in children at age 13 1/2 years: Project Ice Storm. Epigenetics. 2015;10(8):749–61.
Stenz L, Schechter DS, Serpa SR, Paoloni-Giacobino A. Intergenerational Transmission of DNA Methylation Signatures Associated with Early Life Stress. Curr Genomics. 2018;19(8):665–75.
Heijmans BT, Tobi EW, Stein AD, Putter H, Blauw GJ, Susser ES, et al. Persistent epigenetic differences associated with prenatal exposure to famine in humans. Proc Natl Acad Sci U S A. 2008;105(44):17046–9.
Nomura Y, Rompala G, Pritchett L, Aushev V, Chen J, Hurd YL. Natural disaster stress during pregnancy is linked to reprogramming of the placenta transcriptome in relation to anxiety and stress hormones in young offspring. Mol Psychiatry. 2021;26(11):6520–30.
Kotsakis Ruehlmann A, Sammallahti S, Cortés Hidalgo AP, Bakulski KM, Binder EB, Campbell ML, et al. Epigenome-wide meta-analysis of prenatal maternal stressful life events and newborn DNA methylation. Mol Psychiatry. 2023.
Sammallahti S, Cortes Hidalgo AP, Tuominen S, Malmberg A, Mulder RH, Brunst KJ, et al. Maternal anxiety during pregnancy and newborn epigenome-wide DNA methylation. Mol Psychiatry. 2021;26(6):1832–45.
World Health Organization. Coronavirus disease (COVID-19) pandemic. 2021. Available from: https://www.who.int/emergencies/diseases/novel-coronavirus-2019. [cited 2021 May 10].
Patel DI, Gamez Y, Shah L, Patel J. Decline of Psychological Health following the Designation of COVID-19 as a Pandemic: Descriptive Study. J Med Internet Res. 2021;23(4):e24964.
Shah SMA, Mohammad D, Qureshi MFH, Abbas MZ, Aleem S. Prevalence, Psychological Responses and Associated Correlates of Depression, Anxiety and Stress in a Global Population, During the Coronavirus Disease (COVID-19) Pandemic. Community Ment Health J. 2021;57(1):101–10.
Fiorillo A, Gorwood P. The consequences of the COVID-19 pandemic on mental health and implications for clinical practice. Eur Psychiatry. 2020;63(1):e32.
Alexander M, Unruh L, Koval A, Belanger W. United States response to the COVID-19 pandemic, January-November 2020. Health Econ Policy Law. 2022;17(1):62–75.
Thomason ME, Ph D, Graham A, Ph D, Smyser CD, Rogers CE. COPE : Coronavirus Perinatal Experiences – Impact Update ( COPE-IU ). 2020.
Ugalde A, Krishnasamy M, Schofield P. The relationship between self-efficacy and anxiety and general distress in caregivers of people with advanced cancer. J Palliat Med. 2014;17(8):939–41.
Spielberger CD. State-Trait Anxiety Inventory: Bibliography. 2nd ed. Palo Alto: Consulting Psychologists Press; 1989.
Cohen S, Kamarck T, Mermelstein R. A global measure of perceived stress. J Health Soc Behav. 1983;24(4):385–96.
Cohen S, Williamson G. Perceived stress in a probability sample of the United States. Soc Psychol Health. 1988;13:399–411.
Weiss DS, Marmar CR. The impact of event scale: revised BT - Assessing Psychological Trauma and PTSD. In: Assessing Psychological Trauma and PTSD. 1997.
Yuan V, Price EM, Del Gobbo G, Mostafavi S, Cox B, Binder AM, et al. Accurate ethnicity prediction from placental DNA methylation data. Epigenetics Chromatin. 2019;12:51.
Solomon O, MacIsaac J, Quach H, Tindula G, Kobor MS, Huen K, et al. Comparison of DNA methylation measured by Illumina 450K and EPIC BeadChips in blood of newborns and 14-year-old children. Epigenetics. 2018;13(6):655–64.
Pidsley R, Y Wong CC, Volta M, Lunnon K, Mill J, Schalkwyk LC. A data-driven approach to preprocessing Illumina 450K methylation array data. BMC Genomics. 2013;14:293.
Aryee MJ, Jaffe AE, Corrada-Bravo H, Ladd-Acosta C, Feinberg AP, Hansen KD, et al. Minfi: A flexible and comprehensive Bioconductor package for the analysis of Infinium DNA methylation microarrays. Bioinformatics. 2014;30(10):1363–9.
McCartney DL, Walker RM, Morris SW, McIntosh AM, Porteous DJ, Evans KL. Identification of polymorphic and off-target probe binding sites on the Illumina Infinium MethylationEPIC BeadChip. Genom Data. 2016;9:22–4.
Pidsley R, Zotenko E, Peters TJ, Lawrence MG, Risbridger GP, Molloy P, et al. Critical evaluation of the Illumina MethylationEPIC BeadChip microarray for whole-genome DNA methylation profiling. Genome Biol. 2016;17:208.
Maksimovic J, Gordon L, Oshlack A. SWAN: Subset-quantile Within Array Normalization for Illumina Infinium HumanMethylation450 BeadChips. Genome Biol. 2012;13(6):R44.
Ritchie ME, Phipson B, Wu D, Hu Y, Law CW, Shi W, et al. Limma powers differential expression analyses for RNA-sequencing and microarray studies. Nucleic Acids Res. 2015;43(7):e47.
Kocher K, Bhattacharya S, Matthew S, Bramble D, Okitundu-Luwa D, Ngoyi M, et al. DNA methylation patterns associated with konzo in Sub-Saharan Africa. Clin Epigenetics. 2022;14(1). https://doi.org/10.1186/s13148-022-01372-x.
Chen EY, Tan CM, Kou Y, Duan Q, Wang Z, Meirelles GV, et al. Enrichr: Interactive and collaborative HTML5 gene list enrichment analysis tool. BMC Bioinformatics. 2013;14:128.
Kuleshov MV, Jones MR, Rouillard AD, Fernandez NF, Duan Q, Wang Z, et al. Enrichr: a comprehensive gene set enrichment analysis web server 2016 update. Nucleic Acids Res. 2016;44:W90-7.
Krämer A, Green J, Pollard J, Tugendreich S. Causal analysis approaches in ingenuity pathway analysis. Bioinformatics. 2014;30(4):523–30.
Plotly Technologies Inc. Collaborative data science. Montréal: Plotly Technologies Inc; 2015.
Oliveros JC. VENNY. An interactive tool for comparing lists with Venn Diagrams. http://bioinfogp.cnb.csic.es/tools/venny/index.html. BioinfogpCnbCsicEs/Tools/Venny/IndexHtml. 2007.
Dennis CL, Falah-Hassani K, Shiri R. Prevalence of antenatal and postnatal anxiety: Systematic review and meta-analysis. Br J Psychiatry. 2017;210:315–23.
Heron J, O’Connor TG, Evans J, Golding J, Glover V. The course of anxiety and depression through pregnancy and the postpartum in a community sample. J Affect Disord. 2004;80(1):65–73.
Dunkel Schetter C, Lobel M. Pregnancy and birth outcomes: a multilevel analysis of prenatal maternal stress and birth weight. Handbook of health psychology. 2012.
Zhang L, Pan R, Cai Y, Pan J. The prevalence of post-traumatic stress disorder in the general population during the COVID-19 pandemic: a systematic review and single-arm meta-analysis. Psychiatry Investig. 2021;18(5):426–33.
Schwartz DA. An analysis of 38 pregnant women with COVID-19, their newborn infants, and maternal-fetal transmission of SARS-CoV-2: Maternal coronavirus infections and pregnancy outcomes. Arch Pathol Lab Med. 2020;144(7):799–805.
Palma-Gudiel H, Córdova-Palomera A, Eixarch E, Deuschle M, Fañanás L. Maternal psychosocial stress during pregnancy alters the epigenetic signature of the glucocorticoid receptor gene promoter in their offspring: A meta-analysis. Epigenetics. 2015;10:893–902.
Miller O, Shakespeare-Finch J, Bruenig D, Mehta D. DNA methylation of NR3C1 and FKBP5 is associated with posttraumatic stress disorder, posttraumatic growth, and resilience. Psychol Trauma. 2020;12(7):750–5.
Watkeys OJ, Kremerskothen K, Quidé Y, Fullerton JM, Green MJ. Glucocorticoid receptor gene (NR3C1) DNA methylation in association with trauma, psychopathology, transcript expression, or genotypic variation: a systematic review. Neurosci Biobehav Rev. 2018;95:85–122.
Hantsoo L, Kornfield S, Anguera MC, Epperson CN. Inflammation: a proposed intermediary between maternal stress and offspring neuropsychiatric risk. Biological Psychiatry. 2019;85:97–106.
Gordon SM. Interleukin-15 in outcomes of pregnancy. Int J Mol Sci. 2021;22:11094.
Kalantar F, Rajaei S, Heidari AB, Mansouri R, Rashidi N, Izad MH, et al. Serum levels of tumor necrosis factor-α, interleukin-15 and interleukin-10 in patients with pre-eclampsia in comparison with normotensive pregnant women. Iran J Nurs Midwifery Res. 2013;18(6):463–6.
Li J, Li Y, Zhou X, Wei L, Zhang J, Zhu S, et al. Upregulation of IL-15 in the placenta alters trophoblasts behavior contributing to gestational diabetes mellitus. Cell Biosci. 2021;11(1):33.
Bodnar TS, Raineki C, Wertelecki W, Yevtushok L, Plotka L, Zymak-Zakutnya N, et al. Altered maternal immune networks are associated with adverse child neurodevelopment: Impact of alcohol consumption during pregnancy. Brain Behav Immun. 2018;73:205–15.
Lammertink F, Vinkers CH, Tataranno ML, Benders MJNL. Premature Birth and Developmental Programming: Mechanisms of Resilience and Vulnerability. Front Psychiatry. 2021;11:531571.
Lopatina OL, Panina YA, Malinovskaya NA, Salmina AB. Early life stress and brain plasticity: From molecular alterations to aberrant memory and behavior. Rev Neurosci. 2021;32(2):131–42.
Lin MT, Beal MF. Mitochondrial dysfunction and oxidative stress in neurodegenerative diseases. Nature. 2006;443:787–95.
Desplats P, Gutierrez AM, Antonelli MC, Frasch MG. Microglial memory of early life stress and inflammation: Susceptibility to neurodegeneration in adulthood. Neurosci Biobehav Rev. 2020;117:232–42.
Eick SM, Barrett ES, van ‘t Erve TJ, Nguyen RHN, Bush NR, Milne G, et al. Association between prenatal psychological stress and oxidative stress during pregnancy. Paediatr Perinat Epidemiol. 2018;32(4):318–26.
Dowell J, Elser BA, Schroeder RE, Stevens HE. Cellular stress mechanisms of prenatal maternal stress: Heat shock factors and oxidative stress. Neurosci Lett. 2019;709:134368.
Buss C. Maternal oxidative stress during pregnancy and offspring neurodevelopment. Brain Behav Immun. 2021;93:6–7.
Pegoraro C, Figueiredo AL, Maczkowiak F, Pouponnot C, Eychène A, Monsoro-Burq AH. PFKFB4 controls embryonic patterning via Akt signalling independently of glycolysis. Nat Commun. 2015;6:5953.
Goncalves MD, Cantley LC. A Glycolysis Outsider Steps into the Cancer Spotlight. Cell Metabolism. 2018;28:3–4.
Prochaska E, Jang M, Burd I. COVID-19 in pregnancy: Placental and neonatal involvement. Am J Reprod Immunol. 2020;84(5):e13306.
Zhang W, Li Q, Deyssenroth M, Lambertini L, Finik J, Ham J, et al. Timing of prenatal exposure to trauma and altered placental expressions of hypothalamic-pituitary-adrenal axis genes and genes driving neurodevelopment. J Neuroendocrinol. 2018;30(4):e12581.
Gluck O, Grin L, Mizrachi Y, Leytes S, Namazov A, Anteby E, et al. Impact on pregnancy outcomes of exposure to military stress during the first or second trimester as compared with the third trimester. Int J Gynecol Obstet. 2019;146(3):315–20.
Lin YH, Chen CY, Wu SI. Efficiency and quality of data collection among public mental health surveys conducted during the COVID-19 pandemic: Systematic review. J Med Int Res. 2021;23:e25118.
Kocher K, Delot-Vilain A, Spencer DA, LoTempio J, Délot EC. Paucity and Disparity of Publicly Available Sex-Disaggregated Data for the COVID-19 Epidemic Hamper Evidence-Based Decision-Making. Arch Sex Behav. 2021;50(2):407–26.
LoTempio J, Spencer D, Délot-Vilain A, Vilain E, Délot E. We can do better: Lessons learned from the first COVID-19 pandemic. Sci Dipl. 2020;9(2).
Tsai PC, Bell JT. Power and sample size estimation for epigenome-wide association scans to detect differential DNA methylation. Int J Epidemiol. 2015;44(4):1429–41.
The authors would like to sincerely thank all the participants of this study for their time and generous contributions to this research. Additionally, we want to acknowledge the A. James & Alice B. Clark Foundation for funding the DC Mother-Baby Wellness Network, which facilitates recruitment throughout the great DC area to create new interventions to reduce mental health issues commonly experienced by pregnant, low-income women. Additionally, we want to thank the A. James Clark Distinguished Chair of Molecular Genetics for providing funding for the epigenomics experiments presented here.
This study was funded by National Institutes of Health (NHLBI R01 HL116585-01), Intellectual and Developmental Disabilities Research Center (IDDRC), A. James & Alice B. Clark Foundation, and A. James Clark Distinguished Chair of Molecular Genetics. The Institute for Clinical and Translational Science at the University of California, Irvine supports effort for EV, National Center for Advancing Translational Sciences, grant number 5UL1TR001414. The content is solely the responsibility of the authors and does not necessarily represent the official views of the NIH.
Ethics approval and consent to participate
This study was approved by the Children’s National Institutional Review Board under Protocol 14257 and all methods were carried out in accordance with institutional guidelines and regulations. Informed consent to participate in the study was obtained from all participants and recorded in writing.
Consent for publication
The authors declare no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Supplemental Figure S1. Methodology of DNA methylation analysis using 850k EPIC array.
Supplemental Figure S2. Principal Component Analysis of differentially methylated sites for each maternal psychological status assessment instrument: A) STAIS, B) STAI-T, D) PSS, and D) IESR (for RES only). A-C) triangles represent RES while circles represent CTL; all scores (STAI-S, STAI-T, PSS) are broken down in the following: “low” is green, “moderate” is blue, and “high” is red. D) IES-R scores are dichotomized into “no PTSD” (red) and “PTSD” (teal).There were no sites of differential methylation associated with any maternal mental health metric at a statistical threshold of FDR p < 0.05 and log2 fold-change of 1> or <-1.
Supplemental Figure S3. GO enrichment analysis of annotated DMPs between RES and CTL cohorts.
Supplemental Figure S4. Venn diagram of unique and overlapping differentially methylated genes between early and late 2020 pandemic RES subcohorts.
Supplementary Table S1. Maternal psychological stress surveyed throughout pregnancy in pre-pandemic and pandemic cohorts.
Supplementary Table S2. Differentially methylated probes and annotated genes between pre-pandemic and pandemic newborns.
Supplementary Table S3. Pathway analysis of differentially methylated genes associated with pandemic newborns.
Supplementary Table S4. Gene Ontology Enrichment analysis of differentially methylated genes associated with pandemic newborns.
Sites of differential methylation between early and late 2020 pandemic newborns.
About this article
Cite this article
Kocher, K., Bhattacharya, S., Niforatos-Andescavage, N. et al. Genome-wide neonatal epigenetic changes associated with maternal exposure to the COVID-19 pandemic. BMC Med Genomics 16, 268 (2023). https://doi.org/10.1186/s12920-023-01707-4