Next Article in Journal
Angio-Long Noncoding RNA MALAT1 (rs3200401) and MIAT (rs1061540) Gene Variants in Ovarian Cancer
Next Article in Special Issue
Opposite and Differently Altered Postmortem Changes in H3 and H3K9me3 Patterns in the Rat Frontal Cortex and Hippocampus
Previous Article in Journal / Special Issue
Epigenetics of Genes Preferentially Expressed in Dissimilar Cell Populations: Myoblasts and Cerebellum
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Epigenomes
Volume 8
Issue 1
10.3390/epigenomes8010003
epigenomes-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Sex-Specific Associations between Prenatal Exposure to Di(2-ethylhexyl) Phthalate, Epigenetic Age Acceleration, and Susceptibility to Early Childhood Upper Respiratory Infections

by
Sarah M. Merrill
1,2,3,
Nicole Letourneau
4,5,6,7,8,
Gerald F. Giesbrecht
6,7,9,10,
Karlie Edwards
2,3,
Julia L. MacIsaac
2,3,
Jonathan W. Martin
11,
Amy M. MacDonald
12,
David W. Kinniburgh
12,13,
Michael S. Kobor
2,3,14,
Deborah Dewey
6,7,8,10,
Gillian England-Mason
6,7,* and
The APrON Study Team
15,16,†
1
Department of Psychiatry and Human Behavior, The Warren Alpert Medical School at Brown University, Providence, RI 02903, USA
2
Department of Medical Genetics, British Columbia Children’s Hospital Research Institute, University of British Columbia, Vancouver, BC V6T 1Z4, Canada
3
Centre for Molecular Medicine and Therapeutics, Vancouver, BC V6H 0B3, Canada
4
Faculty of Nursing, University of Calgary, Calgary, AB T2N 1N4, Canada
5
Department of Psychiatry, Cumming School of Medicine, University of Calgary, Calgary, AB T2N 1N4, Canada
6
Department of Pediatrics, Cumming School of Medicine, University of Calgary, Calgary, AB T2N 1N4, Canada
7
Owerko Centre, Alberta Children’s Hospital Research Institute, University of Calgary, Calgary, AB T2N 1N4, Canada
8
Hotchkiss Brain Institute, Calgary, AB T2N 4N1, Canada
9
Department of Psychology, Faculty of Arts, University of Calgary, Calgary, AB T2N 1N4, Canada
10
Department of Community Health Sciences, Cumming School of Medicine, University of Calgary, Calgary, AB T2N 1N4, Canada
11
Science for Life Laboratory, Department of Environmental Science, Stockholm University, 106 91 Stockholm, Sweden
12
Alberta Centre for Toxicology, University of Calgary, Calgary, AB T2N 1N4, Canada
13
Department of Laboratory Medicine and Pathology, University of Alberta, Edmonton, AB T6G 2R3, Canada
14
Program in Child and Brain Development, Canadian Institute for Advanced Research, Toronto, ON M5G 1M1, Canada
15
University of Calgary, Calgary, AB T2N 1N4, Canada
16
University of Alberta, Edmonton, AB T6G 2R3, Canada
 
add Show full affiliation list
 
remove Hide full affiliation list
*
Author to whom correspondence should be addressed.
†
Membership of the Group is provided in the Acknowledgments.
Epigenomes 2024, 8(1), 3; https://doi.org/10.3390/epigenomes8010003
Submission received: 4 December 2023 / Revised: 21 January 2024 / Accepted: 23 January 2024 / Published: 26 January 2024
(This article belongs to the Collection Feature Papers in Epigenomes)
Browse Figures
Versions Notes

Abstract

:
Di(2-ethylhexyl) phthalate (DEHP) is a common plasticizer that can affect immune system development and susceptibility to infection. Aging processes (measured as epigenetic age acceleration (EAA)) may mediate the immune-related effects of prenatal exposure to DEHP. This study’s objective was to examine associations between prenatal DEHP exposure, EAA at three months of age, and the number of upper respiratory infections (URIs) from 12 to 18 months of age using a sample of 69 maternal–child pairs from a Canadian pregnancy cohort. Blood DNA methylation data were generated using the Infinium HumanMethylation450 BeadChip; EAA was estimated using Horvath’s pan-tissue clock. Robust regressions examined overall and sex-specific associations. Higher prenatal DEHP exposure (B = 6.52, 95% CI = 1.22, 11.81) and increased EAA (B = 2.98, 95% CI = 1.64, 4.32) independently predicted more URIs. In sex-specific analyses, some similar effects were noted for boys, and EAA mediated the association between prenatal DEHP exposure and URIs. In girls, higher prenatal DEHP exposure was associated with decreased EAA, and no mediation was noted. Higher prenatal DEHP exposure may be associated with increased susceptibility to early childhood URIs, particularly in boys, and aging biomarkers such as EAA may be a biological mechanism. Larger cohort studies examining the potential developmental immunotoxicity of phthalates are needed.

Graphical Abstract

1. Introduction

An upper respiratory infection (URI), also known as the ‘common cold’, is an infection of the nose, sinuses, or throat that can occur year round [1]. URIs are a heterogenous group of infections caused by viruses belonging to several different families, with the most common viral pathogens being rhinoviruses, influenzas A and B, parainfluenza virus, respiratory syncytial virus (RSV), and coronaviruses [2]. Before the age of two, infants and toddlers generally get six to eight colds each year, but some children may experience 10 or more each year [3]. Although there is variation in the symptoms (e.g., chills, nasal congestion, runny nose, coughing, sore throat) and duration of URIs (e.g., one to three weeks), frequent and persistent URIs in young children can indicate an increased susceptibility to infection related to atopic sensitization and subsequent development of atopic diseases such as asthma [4,5]. The increasing prevalence of both acute respiratory illnesses [6] and atopic diseases in children, particularly in those under five years of age [7,8], has prompted researchers to explore environmental factors that may be contributing to this trend [9,10]. It has been recommended that research is needed to investigate whether prenatal exposure to environmental chemicals may be linked to the frequency of early childhood URIs [11,12].
Phthalates are one class of environmental chemicals that may affect susceptibility to URIs and the later development of atopic diseases in children [11,13,14]. They are endocrine-disrupting chemicals (EDCs), and early exposure to phthalates has the potential to alter health and developmental trajectories [15]. One of the most prevalent phthalates is di(2-ethylhexyl) phthalate (DEHP). Exposure to DEHP occurs mainly through ingestion of contaminated food and drinking water, as well as dermal contact with polyvinyl chloride (PVC) plastic products such as vinyl flooring and children’s toys [16,17]. DEHP is commonly incorporated into PVC products in amounts from 40 to 65% and was once thought to have “no known adverse or toxic effects” [18], but increasing epidemiological evidence reports that DEHP is associated with negative health effects, including respiratory system complications [19]. Prenatal exposure to DEHP is of particular concern, as phthalates can transfer to the developing fetus across the placenta, are detectable in placental tissue, amniotic fluid, umbilical cord blood, and meconium, and could have adverse developmental effects on the fetus [20,21,22,23,24]. Based on findings from animal models, prenatal exposure to DEHP has been hypothesized to enhance susceptibility to infections by weakening innate immunity [25,26] and, thus, may contribute to the developmental origins of atopic diseases [27].
The potential biological mechanisms through which prenatal exposure to phthalates may contribute to immunity and infection are not well understood. Some studies in mice have reported that DEHP exposure can increase levels of immunological biomarkers (e.g., immunoglobulin G1 (IgG1), immunoglobulin E (IgE), interleukin 4 (IL4)) and induce airway inflammation [28,29]. However, no consistent pattern of results has emerged from murine research investigating the effects of phthalates on immune system function; an earlier review noted that studies had reported immunopotentiation, immunosuppression, and even null effects [13]. In a recent review of the potential immunotoxicity of phthalates in model organisms, higher DEHP exposure was reported to inhibit immune system function and interfere with immune signaling pathways, with some research indicating that these changes were mediated by DEHP’s effects on gene expression and regulation [30]. Further, DEHP has also been purported to affect other innate immune pathways (e.g., macrophage response) in animal models of disease [31].
Several studies have examined prenatal exposure to DEHP and early immune system function in humans. In studies investigating prenatal exposure to phthalates and immune system indices in umbilical cord blood, prenatal DEHP was not associated with elevated immune biomarkers (i.e., IgE, thymic stromal lymphopoietin (TSL), interleukin-33 (IL-33)) [32,33] or altered immune cell profiles (e.g., cytokines, chemokines) [34]. Further, little is known about the potential effects of prenatal DEHP exposure on childhood susceptibility to respiratory illnesses. A study in the INMA (Infancia y Medio Ambiente; Environment and Childhood) birth cohort from Spain reported that maternal average exposure to DEHP (median concentration = 101.7 μg/g creatinine) during the first and third trimesters of pregnancy was associated with higher odds of children developing wheeze, respiratory tract infections, bronchitis, and asthma from birth until seven years of age; although, the association with respiratory infections was attenuated in adjusted models [11]. In another study in the Odense Child Cohort from Denmark, prenatal DEHP exposure was not significantly associated with risk of the URI symptom of rhinitis in children at age five; although, a trend between higher exposure and greater risk was noted [35]. Given the documented increased susceptibility to respiratory tract infections in children with atopic diseases such as asthma [36,37] and research indicating that DEHP exposure likely plays a role in innate immune response to respiratory infections [26,38], more research is needed that examines the associations between prenatal exposure to DEHP and the frequency of early childhood URIs, as well as potential biological mediators of these associations.
DNA methylation is a plausible biomolecular mediator of the associations between prenatal DEHP exposure and health outcomes, such as susceptibility to URIs. DNA methylation is a stable and mitotically heritable epigenetic marker that can regulate gene expression, and as such, it is ideal for studying the potential biological embedding of early life chemical exposures with later life health [39]. There has been a call for investigations of the associations between early phthalate exposure and epigenetic mechanisms [40,41], particularly using novel tools such as epigenetic clocks (i.e., clocks that use algorithms to calculate the age of tissues based on DNA methylation) [42]. Epigenetic clocks are accurate predictors of DNA methylation-based age, also known as biological age [43]. These epigenetic biomarkers can provide insights into how environmental chemical exposures may affect biological aging and be examined as potential biological mediators. While most epigenetic clocks were developed using adult populations, several perform well in pediatric populations. For example, the Horvath pan-tissue clock [44] accurately estimates DNA methylation age across the lifespan in most human tissues. Epigenetic age acceleration (EAA) as estimated by the Horvath clock has been robustly associated with both early life exposures and later life morbidity and mortality [45], including respiratory illness outcomes [46]. EAA has also been consistently associated with, and even hailed as an emerging biomarker of asthma and atopic diseases in children [47,48] and may provide insight into the biological embedding of early life chemical exposures on immune system development. For example, one study found that accelerated epigenetic age in children was associated both with an increase in allergy and serum IgE [47], an antibody involved in allergic response mechanisms, which has also previously been found to increase with higher prenatal exposure to phthalates [32].
Due to the endocrine-disrupting nature of phthalates, DEHP could have sex-specific effects on epigenetic age. This is supported by a recent study, which reported that higher prenatal DEHP was associated with decreased epigenetic age acceleration (as calculated using the Horvath pan-tissue clock) in males, but not females, at seven years of age [42]. However, to our knowledge, no prior research has examined prenatal DEHP exposure and epigenetic aging rates in infants, or whether EAA may mediate associations between prenatal DEHP exposure and early childhood immune health outcomes, such as susceptibility to infections. Based on the previous literature, it is hypothesized that higher prenatal exposure to DEHP would be associated with a greater number of URIs in early childhood, biological aging would mediate this association, and associations would be sex-specific. To help address the current knowledge gaps, the present study examined the overall and sex-specific associations between prenatal exposure to DEHP, DNA methylation-based aging (i.e., EAA), and early childhood frequency of URIs in a sample of 69 maternal-child pairs from a Canadian pregnancy cohort.

2. Results

2.1. Sample Characteristics

Maternal and child characteristics for the overall sample (Table 1) and sex-stratified groups are presented (Tables S1 and S2). The participant flowchart is also provided (Figure S1). In the overall sample, mothers were mainly university-educated (73.91%), born in Canada (89.86), and had a median household income of >CAD 70,000 (CAD70k) (81.16%); all mothers were married/cohobating (100.00%). For just over half of the women, this was their first pregnancy (52.17% = primiparous). Children (49.28% female) were mainly white (92.75%), born between 35 and 42 weeks of gestation (M = 39.49, SD = 1.40 weeks), and weighed between 2260 and 4904 g (M = 3479.42, SD = 505.91 g). Sample characteristics did not differ substantially between the groups stratified by child sex (n = 34 girls and n = 35 boys).

2.2. Phthalate Exposure

Detection rates, minimums, maximums, geometric means (GMs), and percentiles for the individual DEHP metabolites (μg/g creatinine), mono(2-ethylhexyl) phthalate (MEHP), mono(2-ethyl-5-hydroxy-hexyl) phthalate (MEHHP), mono(2-ethyl-5-oxohexyl) phthalate (MEOHP), and mono(2-ethyl-5-carboxypentyl) phthalate (MECPP), and the molar sum of these four metabolites (μmol/g creatinine) are reported (to three significant figures) for the overall sample population (Table 2) and sex-stratified groups (Table S3). All metabolites were above the limit of detection (LOD) in 100% of maternal second trimester urine samples. Of the metabolites, MECPP had the highest GM (16.6 μg/g creatinine in the overall sample; 16.7 and 16.6 μg/g creatinine in girls and boys, respectively), and the molar sum of DEHP (herein referred to as DEHP)had a GM of 0.0814 μmol/g creatinine in the overall sample and similar levels in each sub-group (0.0824 and 0.0805 μmol/g creatinine in girls and boys, respectively).

2.3. Estimated Cell Type Proportions

Estimated cell type proportions were within the expected ranges for infants who were approximately three-months of age [49]. Preparatory analyses showed that across all 12 estimated cell type proportions, there were no significant associations with either prenatal DEHP exposure (Figure S2) or childhood URIs (Figure S3) within the overall sample. There were also no significant associations in sex-stratified analyses. In girls, cell type proportions were not associated with prenatal DEHP (Figure S4) or childhood URIs (Figure S5). In boys, cell type proportions were also not associated with prenatal DEHP (Figure S6) or childhood URIs (Figure S7).

2.4. Epigenetic Age Acceleration

Descriptive analysis of EAA showed that it was normally distributed and centered around zero (Figure S8). As the age range for the sample collection was narrow (i.e., M = 12.5 weeks, SD = 1.0 week), EAA (i.e., derived from the residuals of a regression of chronological age on predicted biological age) and epigenetic age difference (i.e., derived from the difference of chronological age minus predicted biological age) were identical (Figure S9).

2.5. Overall Associations between Prenatal DEHP Exposure, Epigenetic Age, and Early Childhood URIs

In the main analyses, adjusted robust regression models were used to determine the associations between prenatal DEHP, EAA, and children’s number of early childhood URIs (i.e., colds). In the overall sample, higher prenatal exposure to DEHP was associated with a greater number of early childhood colds (B = 6.52, 95% CI = 1.22, 11.81, adjusted p = 0.07). Prenatal DEHP exposure was not associated with EAA (B = −0.24, 95% CI = −1.24, 0.76, adjusted p = 0.63), but increased EAA was associated with a greater number of early childhood colds in the overall sample (B = 2.98, 95% CI = 1.64, 4.32, adjusted p = 0.001). No evidence of mediation was found in the overall sample (indirect effect = −1.40, 95% CI = −3.94, 1.65, adjusted p = 0.37) (Figure 1).

2.6. Sex-Specific Associations between Prenatal DEHP Exposure, Epigenetic Age, and Early Childhood URIs

In the main analyses, adjusted robust regressions were also used to determine the associations between prenatal DEHP, EAA, and children’s number of colds in sex-stratified models. In girls, prenatal exposure to DEHP was unrelated to the number of early childhood colds (B = −4.67, 95% CI = −13.33, 3.99, adjusted p = 0.28). Prenatal DEHP exposure was associated with decreased EAA (B = −2.24, 95% CI = −3.43, −1.05, adjusted p = 0.001), but EAA was not associated with the number of early childhood colds in girls (B = 1.42, 95% CI = −1.38, 4.22, adjusted p = 0.41). No evidence of mediation was found in girls (indirect effect = −3.32, 95% CI = −8.76, 2.22, adjusted p = 0.22) (Figure 2).
In boys, higher prenatal exposure to DEHP was associated with a greater number of early childhood colds (B = 9.66, 95% CI = 2.39, 16.94, adjusted p = 0.05). Prenatal DEHP exposure was not associated with EAA (B = 0.88, 95% CI = −0.55, 2.31, adjusted p = 0.22), but increased EAA was associated with a greater number of early childhood colds in boys (B = 3.61, 95% CI = 1.64, 5.58, adjusted p = 0.001). We found evidence of mediation (indirect effect = 3.26, 95% CI = 1.80, 4.51, adjusted p = 0.09); specifically, EAA mediated the association between prenatal DEHP and the number of early childhood colds in boys (Figure 3).

2.7. Post Hoc and Sensitivity Analyses

Post hoc analyses examined the associations between prenatal DEHP, other biological aging predictions, and children’s number of colds. The direction and magnitude of effects were similar in models that examined associations with intrinsic epigenetic age acceleration (IEAA), which was EAA corrected for estimated cell type proportions [44] (Table 3). In the analyses, which examined associations with EAA estimated using a different epigenetic clock tool that is tissue-appropriate, the Horvath skin and blood clock, which is less accurate in pediatric samples than the Horvath pan-tissue clock [50], a mediating effect was noted in boys (indirect effect = 3.15, 95% CI = 0.04, 9.52) and in girls (indirect effect = −3.90, 95% CI = −10.56 − 0.23) (Table 3). Estimates were largely unchanged in the sensitivity analysis that adjusted for prenatal exposure to bisphenol A (BPA) (Table S4).

3. Discussion

In this prospective study using data from the Alberta Pregnancy Outcomes and Nutrition (APrON) cohort, we observed that higher maternal second trimester urinary DEHP concentrations and increased epigenetic age acceleration (EAA) were independently associated with a greater number of early childhood upper respiratory infections (URIs). Further, our analyses stratified by child sex revealed several unique associations. In girls, higher prenatal exposure to DEHP was associated with decreased EAA (i.e., slower biological aging). While, in boys, we found some similar effects as those noted in the overall models (i.e., positive total effect and positive path b), but the magnitude of these effects was larger in boys. Also in boys, we found an indirect effect; EAA mediated the association between prenatal exposure to DEHP and early childhood URIs. To our knowledge, no prior research has examined the association between prenatal exposure to DEHP and epigenetic aging rates in infants, nor whether EAA may mediate associations between prenatal DEHP and immune health outcomes. The present analyses were exploratory and limited by the small sample size. Although these findings require replication, this study provides novel evidence suggesting that prenatal exposure to DEHP may increase susceptibility to the ‘common cold’ in young children and that EAA may be a potential biological mechanism, which transmits risk for viral infections.
The common cold is the most frequent infection experienced in childhood and poses an enormous economic burden in terms of healthcare utilization, pediatric consultations, and absences from childcare, school, and the parental workplace [51,52]. We found that higher prenatal exposure to DEHP was associated with a greater number of early childhood colds, particularly in boys (B’s 6.52 and 9.66 in the overall sample and boys, respectively). This finding suggests that in utero exposure to DEHP may exert immunosuppressive effects in children, with potentially more adverse effects in boys. Immune-related effects have been documented following phthalate exposure in animal models; although results are mixed, predominantly immunosuppressive effects are reported [13]. Similarly, phthalate exposure is acknowledged to affect human respiratory health and is thought to play a role in susceptibility to respiratory infections, but the overall evidence is weak and inconsistent [53]. To our knowledge, no prior work has linked prenatal exposure to phthalates to increased incidence of the common cold. A nationally representative cross-sectional study from the United States reported that exposure to other endocrine-disrupting chemicals (EDCs), per- and polyfluoroalkyl substances (PFAS), was associated with greater odds of developing colds in children and adolescents, with larger effects in boys [54]. It is possible that prenatal exposure to endocrine disruptors could alter immune system development and susceptibility to infectious diseases in young children by disrupting the activity of sex steroid hormones [30,55]. Toxicology research indicates that DEHP and its metabolites disrupt androgen signaling pathways and can exert anti-androgenic effects [56,57]. In animal models, prenatal DEHP exposure has been shown to reduce androgen (e.g., testosterone) activity in pregnant dams and male offspring [58,59]. In humans, higher exposure to DEHP has been associated with lower serum testosterone concentrations in pregnant women, regardless of fetal sex [60], as well as in cord blood [61]. Generally, androgens such as testosterone have been shown to demonstrate immunosuppressive effects [55], but the effects of in utero exposure to DEHP on androgen levels in early development and the implications for children’s immune system function are still unknown. Given the significant societal and economic burden associated with the common cold, more research into the underlying mechanisms that transmit risk for respiratory illnesses following early exposure to environmental chemicals is needed. As a male bias is often observed in infectious disease incidence [62], further research examining the role of sex in immune function is strongly recommended.
The possible immune-related dysfunctions that may confer susceptibility to infection following prenatal exposure to DEHP, potentially in a sex-specific manner, are not well understood. In model organisms, DEHP has been shown to contribute to the decline of immune functions through increasing reactive oxygen species (ROS), interfering with the expression of immune genes (e.g., TSLP, TSLPR, IL-7R) and affecting immune cells (e.g., T cells, BK cells) and cytokines [30]. The mechanisms by which prenatal exposure to phthalates interferes with immune cells and immune system function in young girls and boys remains unclear. The Hokkaido Study on Environment and Children’s Health from Japan found that maternal second trimester concentrations of one DEHP metabolite, mono(2-ethylhexyl) phthalate (MEHP), were associated with lower cord blood immunoglobin E (IgE) levels in overall models, but not in sex-stratified analyses [63]. Conversely, both the Maternal-Infant Research on Environmental Chemicals (MIREC) study from Canada and Barwon Infant Study (BIS) from Australia did not find associations between prenatal exposure to DEHP and umbilical cord blood immune system biomarkers in overall or sex-stratified analyses [32] or models adjusted for child sex [34], respectively. Notably, we also did not find any associations with DEHP exposure and altered proportions of immune cell types. However, the MIREC study found that maternal urinary mono(3-carboxypropyl) phthalate (MCPP) concentrations during pregnancy were associated with increased levels of IgE and interleukin-33 (IL-33)/ thymic stromal lymphopoietin (TSLP) in their overall, but not their sex-stratified, analyses [32]. Similarly, other phthalates (e.g., di-n-butyl phthalate (DnBP), molar sum of dibutyl phthalates (DBPs), and benzyl butyl phthalate (BBzP)) have been linked to altered immune cell (e.g., cytokines, chemokines) profiles in cord blood [33,34]. These altered immune profiles are likely related to infection susceptibility [33] and are thought to explain the role of prenatal exposure to phthalates in the etiology of childhood atopy [32,63]. Although the current examination of prenatal DEHP exposure and childhood colds (i.e., a measure related to immune health) is not directly comparable to these other examinations of immunological profiles, this evidence cumulatively suggests that prenatal phthalate exposure may have immunotoxic effects.
One of the proposed modes of action of DEHP during pregnancy that contributes to adverse health effects is through activation of the peroxisome proliferator activated receptor (PPAR) gamma (γ) signaling pathway [64]. This is a potential pathway through which prenatal exposure to DEHP may affect children’s immune system development. For example, the DEHP metabolite MEHP has been shown to inhibit and induce apoptosis of immune cells (i.e., B cells, which produce antibodies to protect from infection) in bone marrow, which is thought to occur through activation of PPARγ [65]. Similarly, a recent study indicated that activation of PPAR alpha (α) and PPARγ by DEHP and MEHP suppresses macrophage (i.e., white blood cells whose primary function is to phagocytize foreign particles such as viruses) activity in the liver [66]. A recent review acknowledged that phthalates likely cause dysregulation of the immune system through diverse mechanisms, including PPARγ, reduced biosynthesis of testosterone, and aryl hydrocarbon receptor (AHR) dysregulation [67]. Thus, given the mixed evidence of cell type and immune marker differences in blood, and the relevance of DEHP exposure for mast cell maturation and macrophage activity, it is possible that other biological samples, rather than blood, may be more appropriate for direct examination of phthalate immunotoxicity in future studies [31,67]. It is tempting to hypothesize that EAA, a cumulative biomarker, was able to detect perturbation in immune function. Although our preparatory analyses did not show that the blood cell proportions accounted for this association, it is possible that EAA may be able to detect differences in immune cell maturation or function due to DEHP exposure. More research is needed to broadly examine the potential immunotoxicity of prenatal phthalates to better understand if they impair innate immune responses [25,26] in a sexually dimorphic manner [68] and whether they contribute to the developmental origins of respiratory and atopic diseases [27].
A central focus of this study was the investigation of the influence of an environmental factor, prenatal exposure to phthalates, on EAA, as it related to a relevant clinical phenotype, early childhood colds. We found thought-provoking associations between epigenetic aging rates in infants and the number of colds from 12 to 18 months of age in the overall sample and in boys; we found that increased epigenetic age acceleration was associated with a greater number of early childhood colds (B’s 2.98 and 3.61 in the overall sample and boys, respectively). In girls, no evidence of an association between EAA and early childhood URIs was found (i.e., adjusted p > 0.10). Previously, DNA methylation alterations have been linked to susceptibility to URIs [69,70], and our findings may elucidate epigenetic aging as a potential DNA methylation-based indicator of infection risk in young children, but additional research is needed to replicate and clarify these sex-specific findings. Epigenetic age deviations in children are likely associated with developmental trajectories, diseases, and certain environmental conditions that may accelerate or decelerate biological development in early life and childhood [71]. Notably, many epigenetic clocks were developed for use with adult populations and have been found to not accurately predict the age of individuals under the age of 20 [43]. The Horvath pan-tissue clock is recommended for pediatric blood samples [50], and positive EAA estimated using this clock has been previously associated with greater odds of atopic sensitization in children from Project Viva [47]. At this time, more work is required to better understand whether pediatric epigenetic age estimates are a driving force or the result of specific pediatric phenotypes [71]. After it has been clarified whether epigenetic age estimates are a driving force or resultant from phenotypes how increased epigenetic age acceleration relates to susceptibility to infection and the experience of early life immune events can be further interpreted.
The developmental origins of health and disease (DOHaD) conceptual model suggests that epigenetic modifications, such as DNA methylation alterations, are a mechanism through which early environmental exposures transmit risk for later health and disease [72]. There has been a call for investigations of early phthalate exposure on DNA methylation profiles [40,41], particularly in infants and young children. We found that higher prenatal exposure to DEHP was associated with slower biological aging, or DNA methylation-predicted age that was ‘younger’ than chronological age in infant girls. This finding is intriguing, as accelerated epigenetic age is typically associated with greater morbidity and mortality in adult populations [73,74], but what ‘faster’ or ‘slower’ biological aging means in pediatric populations is uncertain. In children, it is possible that any deviation from concordance between biological and chronological age—regardless of faster or slower—could be indicative of environmental experiences that may influence future health and development. Potentially, decreased EAA in young children may reflect developmental deceleration or reflect a developmentally younger pattern of DNA methylation [75]. The current understanding of the biological implications of EAA in pediatric populations remains inconclusive and may very well be situationally dependent. The present findings provide additional support for the contention that an epigenetic biomarker of prenatal exposure to phthalates may be evident in infants as young as three months of age [76], and indicates that epigenetic markers of biological aging, EAA, in young infants are sensitive to prenatal exposure to DEHP. This may be influenced by sex differences in aging, as men have been found to display positive EAA compared to women, across epigenetic clocks (e.g., Horvath pan-tissue, Horvath skin and blood, Hannum, GrimAge, PhenoAge, and Pace of Aging) [77,78]. Additional investigations examining associations between prenatal phthalate exposures, as well as other environmental chemicals (e.g., metals, air pollutants, pesticides), and epigenetic clocks in very young pediatric samples may provide novel insights into the early molecular programming of disease susceptibility and contribute to the continued evolution of the DOHaD paradigm [79].
Given the potential sensitivity of the developing immune system to early phthalate exposure, researchers have called for critical policy reforms and stricter regulations on phthalates such as DEHP to protect healthy development in children. The limited evidence on potential biological mechanisms through which phthalates may adversely affect immune system development may be one of the reasons that these chemicals are still considered safe exposures for pregnant women and children [80,81]. This was the first study to examine infant epigenetic age as a potential biomolecular mediator of the immunosuppressive effects of prenatal exposure to DEHP, and future well-powered investigations are needed to replicate and validate the present findings. Based on results from animal models, DNA methylation alterations are hypothesized to underlie the associations between exposure to phthalates and adverse health outcomes in human populations [40,82], including susceptibility to viral infections [83]. We found support for these hypothesized associations in our mediation analyses. Specifically, we found that EAA mediated the association between prenatal DEHP exposure and the number of early childhood colds in boys (indirect effect = 3.26); this finding should be cautiously interpreted. In our supplementary analyses, which examined other tissue-appropriate epigenetic age indices, we also found an indication that EAA (estimated from the Horvath skin and blood clock) mediated this association in boys and girls (indirect effects 3.15 and −3.90 for boys and girls, respectively). These findings suggest that biological aging may be a factor underlying the association between prenatal DEHP exposure and children’s development of upper respiratory illnesses and supports the examination of epigenetic age as a potential mediator of child health outcomes in future investigations. A recent epigenome-wide investigation indicated that prenatal exposure to phthalates may program the activity and signaling pathways of immune cells and sex steroid hormones in infants [76], both of which have important implications for sex differences in the capacity to cope with infections and the later development of inflammatory conditions and autoimmune disorders [84,85]. Additional work is needed to understand how prenatal DEHP exposure contributes to developmental immunotoxicity [80,81] and the biological (e.g., immunological, epigenetic) mechanisms that transmit risk in young girls and boys.
Our statistical confidence in these preliminary findings is strengthened by our analytical approach, which included the use of robust models and correction for multiple comparisons. However, the small sample size, especially given our adjustment for relevant covariates and the smaller sub-samples used for the sex-stratified analyses, is a notable limitation. Our main analyses examined DNA methylation-based age, as predicted by the Horvath pan-tissue clock [44], which is the most widely-accepted clock for estimating epigenetic aging in all age groups, including pediatric populations [71]. We also found similar effect estimates in our post hoc analyses, which used other epigenetic age predictions, and our sensitivity analysis, which adjusted for prenatal exposure to BPA. Our outcome measure, the number of early childhood colds, was obtained from maternal reports and is subject to recall bias. It is possible that unmeasured confounders biased effect estimates. Phthalate monoesters are detectable in human breastmilk, and their concentrations have been found to correlate with important predictors of later immune system outcomes such as hormone concentrations in infants [86]. It is unclear how maternal DEHP exposure may influence breastmilk composition, but it could impact the endocrine control of lactation, thereby disturbing the immune components in breast milk and impacting infant immune system development and response to infections [87]. It is also possible that maternal exposure to DEHP may influence maternal susceptibility to infections and as a result, impact the frequency of URIs in infants. This study did not have data available on maternal breast milk composition or maternal susceptibility to URIs, and these factors should be examined in future research. This was an exploratory examination and additional well-powered studies (e.g., larger cohort studies) are needed to clarify the associations between prenatal phthalate exposures, epigenetic biomarkers, and children’s susceptibility to infection.

4. Materials and Methods

4.1. Participants and Procedure

Participants included a sample of maternal–child pairs (n = 69) recruited between 2009 and 2012 from the APrON study [88,89]. Inclusion criteria for the present study were as follows: (i) a maternal spot urine sample provided during the second trimester of pregnancy, (ii) mothers did not report smoking, consuming alcohol, or receiving steroids during pregnancy, (iii) a venous blood sample obtained from infants at three-months of age, and (iv) data available on the number of colds children had from 12- to 18-months of age (see Figure S1 for participant flowchart). The research protocol was approved by the Conjoint Health Research Ethics Board at the University of Calgary (Ethics ID: REB14-1702). Written, informed consent was obtained from families prior to the collection of biospecimens and the completion of questionnaires.

4.2. Urinary DEHP Assessment

Maternal spot urine samples were collected during the second trimester of pregnancy (M gestational age of 17.0 ± 2.1 weeks). The methods describing the urinary sample collection protocol, quality control experiments, and quantification of phthalate metabolites at the Alberta Centre for Toxicology have been previously described [90,91]. Four monoester phthalate metabolites of DEHP were quantified using liquid chromatography-tandem mass spectrometry (QTRAP 5500, AB Sciex, Concord, ON, Canada). The analytes included MEHP, MEHHP, MEOHP, and MECPP. The limit of detection (LOD) for all metabolites was 0.10 μg/L, and for statistical modeling, all values below the LOD were assigned the value of the LOD/√2 [92]. The molar sum of DEHP metabolites was also calculated [91,93]. Aliquots (1 mL) of the same urine samples were also analyzed quantitatively for creatinine at the Clinical Trials Laboratory, Alberta Health Services (Edmonton, AB, Canada). To account for urinary dilution, creatinine-adjusted concentrations of DEHP (μmol/g creatinine) were used in the analyses.

4.3. Children’s URIs

A sub-sample of APrON mothers and infants participated in an 18-month clinic visit [94]. During this visit, mothers completed a brief infant health and development questionnaire. One item from this questionnaire was used to determine children’s URIs for the present study: “How many colds has your child had in the last 6 months?” Mothers reported on the number of colds their child had in the 6-month (1/2 year) period from 12-months to 18-months of age (mean number of colds = 2.5, SD = 1.5, range: 1–7).

4.4. Epigenetic Age Calculation

4.4.1. Infant Blood Sample Collection

Venous blood samples were collected when infants were three-months of age (M = 2.9, SD = 0.2) using established protocol [49,76]. Venipuncture blood draws were performed by a certified pediatric phlebotomist at the Alberta Children’s Hospital. Whole blood samples were collected and stored for less than 6 h at −80 °C prior to processing. To separate whole blood into plasma, buffy coat, and erythrocytes, samples were centrifuged at 3000 RCF for 15 min. This was completed to separate and aliquot the buffy coat preparation, which consists of peripheral blood mononuclear cells (PBMCs) and some granulocytes, though a much lower proportion of granulocytes than in whole blood preparations. The buffy coat samples were stored at −80 °C and shipped on dry ice to the Department of Medical Genetics at the University of British Columbia for DNA extraction and quantification.

4.4.2. DNA Methylation Assay and Data Processing

DNA methylation data were extracted and processed using established protocol [49,76]. In brief, a sample of 750 ng of genomic DNA was used for bisulfite conversion using the EZ DNA Methylation Kit (Zymo Research, Irvine, CA, USA). Following the manufacturer’s protocol, quantification of DNA methylation status of purified bisulfite converted DNA was performed using the Infinium HumanMethylation450 (450k) Bead Chip Assay (Illumina, San Diego, CA, USA). This platform quantified DNA methylation status at >450,000 CpG sites; this covers 99% of available reference sequence genes [95]. A sample of 160 ng of bisulfite converted DNA was whole genome amplified, fragmented, and hybridized to the BeadChip arrays. Via the 450k platform, each CpG site is targeted by a probe to distinguish “methylated” and “unmethylated” intensity through different dye colors (i.e., green and red). Methylation levels, also known as beta (β) values, were calculated. This was completed in Illumina GenomeStudio (San Diego, CA, USA) by dividing the methylation probe signal intensity by the sum of the methylated and unmethylated probe signal intensities). The β values range from 0 (i.e., completely unmethylated site) to 1 (i.e., fully methylated site). The β values were extracted and imported into R statistical software for processing and analysis. Further detail on the 450k arrays for the APrON cohort has been previously published [76]. The raw data were pre-processed through the same pipeline as previously used [49,76]; this study used the pre-processed DNA methylation data to calculate epigenetic age.

4.4.3. Estimated Cell Type Proportions

Defining cell type identity is a major biological function of DNA methylation, and as such, proportions of immune cells making up blood can be defined using DNA methylation profiles. Therefore, cell type proportions in buffy coat blood were estimated using constraint projection [96] on the IDOL (Identifying Optimal Libraries) extended cell type reference data set and version 1.8.0 of the R package FlowSorted.BloodExtended.EPIC [97,98,99] to derive 12 cell type proportions: neutrophils, eosinophils, basophils, naïve and memory B cells, natural killer (NK) cells, naïve and memory CD4-T cells, T regulatory cells, naïve and memory CD8-T cells, and monocytes. Our samples were from young infant blood. Thus, based on the approach previously used [49,100], available cord blood references for deconvolution were not used, as these references account for nucleated red blood cells (NRBCs; immature blood cells), a cell type not found in the blood of infants at three-months of age [101].

4.4.4. Epigenetic Clock

In the main analysis, epigenetic age, or biological age, as determined by the pattern of β values in the methylome, was calculated using the Horvath pan-tissue clock, and in the post hoc analysis using the Horvath skin and blood clock, using the online age calculator: https://dnamage.genetics.ucla.edu/ (accessed on 22 February 2020) [44]. These measures are appropriately applied to venous blood and have been validated and employed frequently in research [102], as well as in similarly aged children [102,103]. A recent review concluded that the Horvath pan-tissue clock, a widely used age estimator based on DNA methylation at 353 CpG sites developed using DNA methylation data from over 8000 samples, was the most accurate in pediatric blood samples [50]. A normalized, background- and color-corrected β file was provided for the calculator. Epigenetic age acceleration (EAA) was defined as the residuals from regressing the calculated epigenetic age on the infants’ actual chronological age. This provides a measure of whether infants are aging faster (increased EAA) or slower (decreased EAA) biologically than their chronological age; this was the main biological aging estimate used in the present analyses.
Other biological aging predictions were also calculated and compared with EAA. Epigenetic age difference was calculated as the estimated epigenetic age subtracted from the reported chronological age. The intrinsic epigenetic age acceleration (IEAA) measure was outputted by the calculator as residuals from the regression of biological age on chronological age, accounting for reference-estimated cell type proportions [44]. Lastly, in a post hoc analysis, EAA was calculated using the Horvath skin and blood clock to attempt to replicate our findings in this same sample with a tool employing a unique set of 391 CpGs [104].

4.5. Covariates

Based on variables previously identified to be associated with prenatal phthalate exposure, epigenetic modifications, and/or children’s colds [11,76,105], this study considered several covariates. Potential covariates included the following maternal characteristics: age, education, household income, self-reported race, parity, history of asthma, and pre-pregnancy body mass index (BMI). Potential covariates also included the following child characteristics: gestational age at birth, birthweight, month of birth, and childcare arrangements at 12-months of age. Dichotomous variables were coded as follows: household income (i.e., <$70,000 CAD versus >$70,000 CAD), maternal education (i.e., undergraduate university degree or higher versus trade school/high school diploma/lower), maternal self-reported race (i.e., White versus not White), maternal parity (i.e., primiparous versus multiparous), maternal history of asthma (i.e., yes versus no), and childcare arrangements (i.e., childcare utilized versus not utilized). All mothers (n = 69) reported breastfeeding at three-months postpartum; this was not examined as a covariate.

4.6. Statistical Analyses

All analyses were performed in R (version 4.3.1). To describe study sample characteristics, means (i.e., geometric means for phthalate metabolites, otherwise arithmetic means are reported) and standard deviations (SD) for continuous variables and proportions for categorical variables were used.
Preparatory analyses were conducted to describe and characterize the estimated cell type proportions and EAA. The associations between estimated cell type proportions and DEHP and URIs were examined using linear regressions with Bonferroni-corrected p-value significance threshold of p < 0.004.
The main analyses involved examining the associations between prenatal exposure to DEHP, EAA, and URIs. This was completed utilizing robust multivariable regressions. These regressions were completed using the R package MASS and with Huber M-estimation [106,107]. Overall models (n = 69) and models stratified by child sex (n = 34 girls; n = 35 boys) were examined. Covariates were included in final robust models if they were associated with the outcome at p < 0.20 [108]. Overall models were adjusted for three covariates (i.e., maternal self-reported race, maternal parity, and birth month), female-stratified models were adjusted for two covariates (i.e., maternal age and self-reported race), and male-stratified models were adjusted for four covariates (i.e., maternal education, maternal parity, maternal history of asthma, and child month of birth). Robust mediation tests were performed (using ZYmediate from the R package WRS2) using 10,000 bootstrapped samples and the method proposed by Zu and Yuan [109]. In regression approaches, mediating effects (i.e., indirect effects) can be examined even in the absence of other effects (e.g., direct effects, relation of X to M) [110]. Mediation analysis results are reported using the widely applied model notation [111]: relation of X to Y (path c; total effect), relation of X to M (path a), relation of M to Y adjusted for X (path b), relation of X to Y adjusted for M (path c’; direct effect), mediating effect (path ab; indirect effect). All tests were 2-sided. We corrected for multiple comparisons in the robust regression analyses using the Benjamini–Hochberg procedure to control the false discovery rate (FDR; i.e., proportion of false positives) [112]. Given the limited evidence on the effect of prenatal exposure to phthalates on early childhood EAA and immune health outcomes, we adopted an exploratory approach and report adjusted p-values (q-values) from 0.05 to 0.1. The adjusted p-value (q-value) thresholds of 0.05 and 0.10 yield FDRs of 5% and 10%, respectively. Thus, we can have the most statistical confidence in results at adjusted p < 0.05. Simulation analysis indicated that a sample size of 59 would be adequately powered (0.80) to detect a mediation effect, when using a percentile bootstrap method, and considering medium-to-large path a and b effects [113]. A larger sample size (≥78) would be required when considering smaller path a and b effects [113], and thus, these analyses are considered exploratory.
We did not examine EAA predictions from other epigenetic clocks in the main analyses, due to their poor performance in children [50]. However, post hoc analyses explored associations with (i) intrinsic epigenetic age acceleration (IEAA), which is corrected for cell types [44], and (ii) EAA estimated using Horvath’s skin and blood clock [104]. Sensitivity analysis also examined associations when prenatal exposure to another non-persistent endocrine disruptor, bisphenol A (BPA), was included as a covariate.

5. Conclusions

This exploratory study found that higher prenatal DEHP exposure and increased EAA were independently associated with a greater number of early childhood URIs, particularly in boys. Further, we found that EAA mediated the association between prenatal DEHP exposure and the number of URIs in boys. Lastly, we found that prenatal exposure to DEHP was associated with a decreased EAA in girls, possibly indicating developmental deceleration. As early childhood constitutes a window of susceptibility for the maturation and priming of the immune system, it is critical to examine the potential developmental immunotoxicity of phthalates and the sex-specific effects of these EDCs.

Supplementary Materials

The following supporting information can be downloaded at https://www.mdpi.com/article/10.3390/epigenomes8010003/s1. Figure S1. Flowchart of participants included in the present study; Figure S2. Estimated cell type proportions in three-month-old infant venous blood were not associated with prenatal DEHP exposure across the entire cohort (n = 69). Each graph represents the association between each of the 12 estimated cell type proportions (on the y-axis) and DEHP exposure (on the x-axis). On each graph is the Pearson’s correlation (R) and p-value of the correlation between estimated cell type proportion and DEHP exposure; Figure S3. Estimated cell type proportions in three-month-old infant venous blood were not associated with childhood upper respiratory infections (URIs) across the entire cohort (n = 69). Each graph represents the association between each of the 12 estimated cell type proportions (on the y-axis) and childhood URIs (on the x-axis). On each graph is the Pearson’s correlation (R) and p-value of the correlation between estimated cell type proportion and childhood URIs; Figure S4. Estimated cell type proportions in three-month-old infant girls’ venous blood were not associated with prenatal DEHP exposure (n = 34). Each graph represents the association between each of the 12 estimated cell type proportions (on the y-axis) and DEHP exposure (on the x-axis). On each graph is the Pearson’s correlation (R) and p-value of the correlation between estimated cell type proportion and DEHP exposure; Figure S5. Estimated cell type proportions in three-month-old infant girls’ venous blood were not associated with childhood upper respiratory infections (URIs) (n = 34). Each graph represents the association between each of the 12 estimated cell type proportions (on the y-axis) and childhood URIs (on the x-axis). On each graph is the Pearson’s correlation (R) and p-value of the correlation between estimated cell type proportion and childhood URIs; Figure S6. Estimated cell type proportions in three-month-old infant boys’ venous blood were not associated with prenatal DEHP exposure (n = 35). Each graph represents the association between each of the 12 estimated cell type proportions (on the y-axis) and DEHP exposure (on the x-axis). On each graph is the Pearson’s correlation (R) and p-value of the correlation between estimated cell type proportion and DEHP exposure; Figure S7. Estimated cell type proportions in three-month-old infant boys’ venous blood were not associated with childhood upper respiratory infections (URIs) (n = 35). Each graph represents the association between each of the 12 estimated cell type proportions (on the y-axis) and childhood URIs (on the x-axis). On each graph is the Pearson’s correlation (R) and p-value of the correlation between estimated cell type proportion and childhood URIs; Figure S8. Histogram with kernel density plot of epigenetic age acceleration (EAA) values estimated from the Horvath pan-tissue clock in these samples. The dashed line indicates the mean, which is centered at approximately zero, and the data are approximately normally distributed; Figure S9. Estimated epigenetic age acceleration (EAA) and epigenetic age difference from the Horvath pan-tissue clock were the same in this sample due to the narrow age range during blood sample collection. The Horvath pan-tissue EAA was derived from residuals of a regression of chronological age on predicted biological age is on the x-axis. While Horvath pan-tissue epigenetic age difference was derived from the difference of chronological age minus predicted biological age. As the correlation is so high, these measures of EAA are synonymous in this sample; Table S1. Maternal and child characteristics for the sub-sample stratified by child sex (n = 34 girls); Table S2. Maternal and child characteristics for the sub-sample stratified by child sex (n = 35 boys); Table S3. Creatinine-adjusted prenatal phthalate metabolite concentrations (μg/g creatinine) and DEHP (μmol/g creatinine) in maternal second trimester urine for the sex-stratified groups; Table S4. Sensitivity analysis adjusting for prenatal exposure to bisphenol A (BPA).

Author Contributions

S.M.M. conceptualized and designed the study, carried out the data analyses, and critically reviewed and revised the manuscript. N.L., G.F.G., J.W.M., M.S.K. and D.D. obtained study funding, designed the data collection instruments, coordinated and supervised collection of data, and critically reviewed and revised the manuscript. K.E., J.L.M., A.M.M. and D.W.K. carried out the data analyses and critically reviewed and revised the manuscript. G.E.-M. conceptualized and designed the study, carried out the data analyses, drafted the initial manuscript, and critically reviewed and revised the manuscript. The APrON Study Team obtained study funding and designed data collection instruments. All authors have read and agreed to the published version of the manuscript.

Funding

The APrON cohort was established by an interdisciplinary team grant provided by Alberta Innovates Health Solutions (formally the Alberta Heritage Foundation for Medical Research). The collection and analysis of data presented in this manuscript was enabled through additional funding provided by the Canadian Institutes of Health Research (CIHR) (MOP-123535), Alberta Children’s Hospital Foundation, Palix Foundation, Allergy, Genes and Environment Network (AllerGen NCE Inc., Calgary, AB, USA), and U.S. National Institutes of Health (Exploration/Development Grant 1R21ES021295-01R21). G. England-Mason received salary support through a Fellowship provided by CIHR (HTA-472411) and Postgraduate Fellowship in Health Innovation provided by Alberta Innovates, the Ministry of Economic Development, Trade and Tourism, and the Government of Alberta.

Institutional Review Board Statement

The study was conducted in accordance with the Declaration of Helsinki/The Study protocol was approved by the research ethics board at the University of Calgary (Ethics ID: REB14-1702) on 15 January 2009.

Informed Consent Statement

Women gave written, informed consent at the time of recruitment and prior to the collection of biospecimens and completion of questionnaires.

Data Availability Statement

Data are contained within the article and Supplementary Materials.

Acknowledgments

We are extremely grateful to all the families who took part in this study and the APrON team (http://www.apronstudy.ca/), investigators, research assistants, graduate and undergraduate students, volunteers, clerical staff, and managers. We acknowledge the significant contributions of the APrON Study Team whose individual members are B.J. Kaplan, D. Dewey, C.J. Field, R.C. Bell, F.P. Bernier, M. Cantell, L.M. Casey, M. Eliasziw, A. Farmer, L. Gagnon, G.F. Giesbrecht, L. Goonewardene, D. Johnston, L. Kooistra, N. Letourneau, D.P. Manca, J. W. Martin, L.J. McCargar, M. O’Beirne, V.J. Pop, A.J. Deane, and N. Singha, and the APrON Management Team: N. Letourneau (current PI), R.C. Bell, D. Dewey, C.J. Field, L. Forbes, G. Giesbrecht, C. Lebel, B. Leung, C. McMorris, and K. Ross. The extended cell type reference used for this research was enabled, in part, by the use of the Software packaged developed by at Dartmouth College, which software to the licensing terms made available by Dartmouth Technology Transfer and which software is provided “AS IS” with no warranties whatsoever.

Conflicts of Interest

The authors do not have any conflicts of interest to declare. The funding sources did not have any role in the design of the study; in the data collection, analyses, or interpretation; in the writing of the manuscript; or in the decision to publish the results.

References

  1. Worrall, G. Common Cold. Can. Fam. Physician 2011, 57, 1289–1290. [Google Scholar]
  2. Pappas, D.E. The Common Cold. Princ. Pract. Pediatr. Infect. Dis. 2018, 199–202.e1. [Google Scholar] [CrossRef]
  3. Canadian Paediatric Society Infectious Diseases and Immunization Committee Colds in Children. Paediatr. Child. Health 2005, 10, 493–495.
  4. Di Cicco, M.; D’Elios, S.; Peroni, D.G.; Comberiati, P. The Role of Atopy in Asthma Development and Persistence. Curr. Opin. Allergy Clin. Immunol. 2020, 20, 131–137. [Google Scholar] [CrossRef] [PubMed]
  5. Jartti, T.; Bønnelykke, K.; Elenius, V.; Feleszko, W. Role of Viruses in Asthma. Semin. Immunopathol. 2020, 42, 61–74. [Google Scholar] [CrossRef]
  6. Ma, K.C.; Winn, A.; Moline, H.L.; Scobie, H.M.; Midgley, C.M.; Kirking, H.L.; Adjemian, J.; Hartnett, K.P.; Johns, D.; Jones, J.M.; et al. Increase in Acute Respiratory Illnesses Among Children and Adolescents Associated with Rhinoviruses and Enteroviruses, Including Enterovirus D68—United States, July–September 2022. MMWR Morb. Mortal. Wkly. Rep. 2022, 71, 1265–1270. [Google Scholar] [CrossRef]
  7. Mohn, C.H.; Blix, H.S.; Halvorsen, J.A.; Nafstad, P.; Valberg, M.; Lagerløv, P. Incidence Trends of Atopic Dermatitis in Infancy and Early Childhood in a Nationwide Prescription Registry Study in Norway. JAMA Netw. Open 2018, 1, e184145. [Google Scholar] [CrossRef] [PubMed]
  8. Yang, E.J.; Sekhon, S.; Sanchez, I.M.; Beck, K.M.; Bhutani, T. Recent Developments in Atopic Dermatitis. Pediatrics 2018, 142, e20181102. [Google Scholar] [CrossRef] [PubMed]
  9. Garner, R.; Kohen, D. Changes in the Prevalence of Asthma among Canadian Children. Health Rep. 2008, 19, 45–50. [Google Scholar]
  10. Gershon, A.S.; Guan, J.; Wang, C.; To, T. Trends in Asthma Prevalence and Incidence in Ontario, Canada, 1996–2005: A Population Study. Am. J. Epidemiol. 2010, 172, 728–736. [Google Scholar] [CrossRef]
  11. Gascon, M.; Casas, M.; Morales, E.; Valvi, D.; Ballesteros-Gómez, A.; Luque, N.; Rubio, S.; Monfort, N.; Ventura, R.; Martínez, D.; et al. Prenatal Exposure to Bisphenol A and Phthalates and Childhood Respiratory Tract Infections and Allergy. J. Allergy Clin. Immunol. 2015, 135, 370–378. [Google Scholar] [CrossRef]
  12. Sotir, M.; Yeatts, K.; Shy, C. Presence of Asthma Risk Factors and Environmental Exposures Related to Upper Respiratory Infection-Triggered Wheezing in Middle School-Age Children. Environ. Health Perspect. 2003, 111, 657–662. [Google Scholar] [CrossRef] [PubMed]
  13. Kimber, I.; Dearman, R.J. An Assessment of the Ability of Phthalates to Influence Immune and Allergic Responses. Toxicology 2010, 271, 73–82. [Google Scholar] [CrossRef] [PubMed]
  14. Robinson, L.; Miller, R. The Impact of Bisphenol A and Phthalates on Allergy, Asthma, and Immune Function: A Review of Latest Findings. Curr. Environ. Health Rep. 2015, 2, 379–387. [Google Scholar] [CrossRef]
  15. Kahn, L.G.; Philippat, C.; Nakayama, S.F.; Slama, R.; Trasande, L. Endocrine-Disrupting Chemicals: Implications for Human Health. Lancet Diabetes Endocrinol. 2020, 8, 703–718. [Google Scholar] [CrossRef]
  16. Koch, H.M.; Lorber, M.; Christensen, K.L.Y.; Pälmke, C.; Koslitz, S.; Brüning, T. Identifying Sources of Phthalate Exposure with Human Biomonitoring: Results of a 48h Fasting Study with Urine Collection and Personal Activity Patterns. Int. J. Hydrogen Environ. Health 2013, 216, 672–681. [Google Scholar] [CrossRef]
  17. Pacyga, D.C.; Sathyanarayana, S.; Strakovsky, R.S. Dietary Predictors of Phthalate and Bisphenol Exposures in Pregnant Women. Adv. Nutr. 2019, 10, 803–815. [Google Scholar] [CrossRef] [PubMed]
  18. McKeen, L.W. 3—Plastics Used in Medical Devices. In Handbook of Polymer Applications in Medicine and Medical Devices; Modjarrad, K., Ebnesajjad, S., Eds.; Plastics Design Library; William Andrew Publishing: Oxford, UK, 2014; pp. 21–53. ISBN 978-0-323-22805-3. [Google Scholar]
  19. Zarean, M.; Keikha, M.; Poursafa, P.; Khalighinejad, P.; Amin, M.; Kelishadi, R. A Systematic Review on the Adverse Health Effects of Di-2-Ethylhexyl Phthalate. Environ. Sci. Pollut. Res. 2016, 23, 24642–24693. [Google Scholar] [CrossRef]
  20. Guo, J.; Wu, M.; Gao, X.; Chen, J.; Li, S.; Chen, B.; Dong, R. Meconium Exposure to Phthalates, Sex and Thyroid Hormones, Birth Size and Pregnancy Outcomes in 251 Mother–Infant Pairs from Shanghai. Int. J. Environ. Res. Public. Health 2020, 17, 7711. [Google Scholar] [CrossRef]
  21. Katsikantami, I.; Tzatzarakis, M.N.; Alegakis, A.K.; Karzi, V.; Hatzidaki, E.; Stavroulaki, A.; Vakonaki, E.; Xezonaki, P.; Sifakis, S.; Rizos, A.K.; et al. Phthalate Metabolites Concentrations in Amniotic Fluid and Maternal Urine: Cumulative Exposure and Risk Assessment. Toxicol. Rep. 2020, 7, 529–538. [Google Scholar] [CrossRef]
  22. Li, L.-X.; Chen, L.; Meng, X.-Z.; Chen, B.-H.; Chen, S.-Q.; Zhao, Y.; Zhao, L.-F.; Liang, Y.; Zhang, Y.-H. Exposure Levels of Environmental Endocrine Disruptors in Mother-Newborn Pairs in China and Their Placental Transfer Characteristics. PLoS ONE 2013, 8, e62526. [Google Scholar] [CrossRef]
  23. Mose, T.; Mortensen, G.K.; Hedegaard, M.; Knudsen, L.E. Phthalate Monoesters in Perfusate from a Dual Placenta Perfusion System, the Placenta Tissue and Umbilical Cord Blood. Reprod. Toxicol. 2007, 23, 83–91. [Google Scholar] [CrossRef]
  24. Mose, T.; Knudsen, L.E.; Hedegaard, M.; Mortensen, G.K. Transplacental Transfer of Monomethyl Phthalate and Mono(2-Ethylhexyl) Phthalate in a Human Placenta Perfusion System. Int. J. Toxicol. 2007, 26, 221–229. [Google Scholar] [CrossRef]
  25. Shin, I.-S.; Lee, M.-Y.; Cho, E.-S.; Choi, E.; Son, H.-Y.; Lee, K.-Y. Effects of Maternal Exposure to Di(2-Ethylhexyl)Phthalate (DEHP) during Pregnancy on Susceptibility to Neonatal Asthma. Toxicol. Appl. Pharmacol. 2014, 274, 402–407. [Google Scholar] [CrossRef]
  26. Walker, C.; Ghazisaeidi, S.; Collet, B.; Boisvert, A.; Culty, M. In Utero Exposure to Low Doses of Genistein and Di-(2-Ethylhexyl) Phthalate (DEHP) Alters Innate Immune Cells in Neonatal and Adult Rat Testes. Andrology 2020, 8, 943–964. [Google Scholar] [CrossRef]
  27. Yen, P.-L.; Yang, C.-R.; Huang, M.-L.; Lin, T.-A.; Liao, V.H.-C. Chronic Exposure to Di(2-Ethylhexyl) Phthalate (DEHP) Weakens Innate Immunity and Leads to Immunosenescence in C. Elegans. Environ. Toxicol. Pharmacol. 2023, 98, 104071. [Google Scholar] [CrossRef] [PubMed]
  28. Guo, J.; Han, B.; Qin, L.; Li, B.; You, H.; Yang, J.; Liu, D.; Wei, C.; Nanberg, E.; Bornehag, C.-G.; et al. Pulmonary Toxicity and Adjuvant Effect of Di-(2-Exylhexyl) Phthalate in Ovalbumin-Immunized BALB/c Mice. PLoS ONE 2012, 7, e39008. [Google Scholar] [CrossRef] [PubMed]
  29. Larsen, S.T.; Hansen, J.S.; Hansen, E.W.; Clausen, P.A.; Nielsen, G.D. Airway Inflammation and Adjuvant Effect after Repeated Airborne Exposures to Di-(2-Ethylhexyl)Phthalate and Ovalbumin in BALB/c Mice. Toxicology 2007, 235, 119–129. [Google Scholar] [CrossRef] [PubMed]
  30. Zhang, Y.; Lyu, L.; Tao, Y.; Ju, H.; Chen, J. Health Risks of Phthalates: A Review of Immunotoxicity. Environ. Pollut. 2022, 313, 120173. [Google Scholar] [CrossRef] [PubMed]
  31. Schmidt, S. Targeting the Macrophage: Immune Cells May Be the Key to Phthalate-Induced Liver Toxicity. Environ. Health Perspect. 2022, 130, 034003. [Google Scholar] [CrossRef] [PubMed]
  32. Ashley-Martin, J.; Dodds, L.; Levy, A.R.; Platt, R.W.; Marshall, J.S.; Arbuckle, T.E. Prenatal Exposure to Phthalates, Bisphenol A and Perfluoroalkyl Substances and Cord Blood Levels of IgE, TSLP and IL-33. Environ. Res. 2015, 140, 360–368. [Google Scholar] [CrossRef]
  33. Tsai, C.-K.; Cheng, H.-H.; Hsu, T.-Y.; Wang, J.-Y.; Hung, C.-H.; Tsai, C.-C.; Lai, Y.-J.; Lin, Y.-J.; Huang, H.-C.; Chan, J.Y.H.; et al. Prenatal Exposure to Di-Ethyl Phthalate (DEP) Is Related to Increasing Neonatal IgE Levels and the Altering of the Immune Polarization of Helper-T Cells. Int. J. Environ. Res. Public. Health 2021, 18, 6364. [Google Scholar] [CrossRef]
  34. Eisner, A.; Gao, Y.; Collier, F.; Drummond, K.; Thomson, S.; Burgner, D.; Vuillermin, P.; Tang, M.L.; Mueller, J.; Symeonides, C.; et al. Cord Blood Immune Profile: Associations with Higher Prenatal Plastic Chemical Levels. Environ. Pollut. 2022, 315, 120332. [Google Scholar] [CrossRef]
  35. Jøhnk, C.; Høst, A.; Husby, S.; Schoeters, G.; Timmermann, C.A.G.; Kyhl, H.B.; Beck, I.H.; Andersson, A.-M.; Frederiksen, H.; Jensen, T.K. Maternal Phthalate Exposure and Asthma, Rhinitis and Eczema in 552 Children Aged 5 Years; a Prospective Cohort Study. Environ. Health 2020, 19, 32. [Google Scholar] [CrossRef] [PubMed]
  36. James, K.M.; Peebles, R.S.; Hartert, T.V. Response to Infections in Patients with Asthma and Atopic Disease: An Epiphenomenon or Reflection of Host Susceptibility? J. Allergy Clin. Immunol. 2012, 130, 343–351. [Google Scholar] [CrossRef]
  37. Kansen, H.M.; Lebbink, M.A.; Mul, J.; van Erp, F.C.; van Engelen, M.; de Vries, E.; Prevaes, S.M.P.J.; Le, T.M.; van der Ent, C.K.; Verhagen, L.M. Risk Factors for Atopic Diseases and Recurrent Respiratory Tract Infections in Children. Pediatr. Pulmonol. 2020, 55, 3168–3179. [Google Scholar] [CrossRef]
  38. Winans, B.; Humble, M.C.; Lawrence, B.P. Environmental Toxicants and the Developing Immune System: A Missing Link in the Global Battle against Infectious Disease? Reprod. Toxicol. 2011, 31, 327–336. [Google Scholar] [CrossRef]
  39. Aristizabal, M.J.; Anreiter, I.; Halldorsdottir, T.; Odgers, C.L.; McDade, T.W.; Goldenberg, A.; Mostafavi, S.; Kobor, M.S.; Binder, E.B.; Sokolowski, M.B.; et al. Biological Embedding of Experience: A Primer on Epigenetics. Proc. Natl. Acad. Sci. USA 2020, 117, 23261–23269. [Google Scholar] [CrossRef]
  40. Benjamin, S.; Masai, E.; Kamimura, N.; Takahashi, K.; Anderson, R.C.; Faisal, P.A. Phthalates Impact Human Health: Epidemiological Evidences and Plausible Mechanism of Action. J. Hazard. Mater. 2017, 340, 360–383. [Google Scholar] [CrossRef] [PubMed]
  41. Dutta, S.; Haggerty, D.K.; Rappolee, D.A.; Ruden, D.M. Phthalate Exposure and Long-Term Epigenomic Consequences: A Review. Front. Genet. 2020, 11, 405. [Google Scholar] [CrossRef] [PubMed]
  42. Khodasevich, D.; Holland, N.; Hubbard, A.; Harley, K.; Deardorff, J.; Eskenazi, B.; Cardenas, A. Associations between Prenatal Phthalate Exposure and Childhood Epigenetic Age Acceleration. Environ. Res. 2023, 231, 116067. [Google Scholar] [CrossRef] [PubMed]
  43. Simpson, D.J.; Chandra, T. Epigenetic Age Prediction. Aging Cell 2021, 20, e13452. [Google Scholar] [CrossRef]
  44. Horvath, S. DNA Methylation Age of Human Tissues and Cell Types. Genome Biol. 2013, 14, R115. [Google Scholar] [CrossRef]
  45. Binder, A.M.; Horvath, S. Epigenetic Clocks. In Epigenetic Epidemiology; Michels, K.B., Ed.; Springer International Publishing: Cham, Switzerland, 2022; pp. 261–276. ISBN 978-3-030-94475-9. [Google Scholar]
  46. Bejaoui, Y.; Amanullah, F.H.; Saad, M.; Taleb, S.; Bradic, M.; Megarbane, A.; Hssain, A.A.; Khalil, C.A.; Hajj, N.E. Epigenetic Age Acceleration in Surviving versus Deceased COVID-19 Patients with Acute Respiratory Distress Syndrome Following Hospitalization. Clin. Epigenetics 2023, 15, 186. [Google Scholar] [CrossRef] [PubMed]
  47. Peng, C.; Cardenas, A.; Rifas-Shiman, S.L.; Hivert, M.-F.; Gold, D.R.; Platts-Mills, T.A.; Lin, X.; Oken, E.; Avila, L.; Celedón, J.C.; et al. Epigenetic Age Acceleration Is Associated with Allergy and Asthma in Children in Project Viva. J. Allergy Clin. Immunol. 2019, 143, 2263–2270.e14. [Google Scholar] [CrossRef]
  48. Vasileva, D.; Greenwood, C.M.T.; Daley, D. A Review of the Epigenetic Clock: Emerging Biomarkers for Asthma and Allergic Disease. Genes 2023, 14, 1724. [Google Scholar] [CrossRef] [PubMed]
  49. Merrill, S.M.; Gladish, N.; Fu, M.P.; Moore, S.R.; Konwar, C.; Giesbrecht, G.F.; MacIssac, J.L.; Kobor, M.S.; Letourneau, N.L. Associations of Peripheral Blood DNA Methylation and Estimated Monocyte Proportion Differences during Infancy with Toddler Attachment Style. Attach. Hum. Dev. 2023, 25, 132–161. [Google Scholar] [CrossRef]
  50. Fang, F.; Zhou, L.; Perng, W.; Marsit, C.J.; Knight, A.K.; Cardenas, A.; Aung, M.T.; Hivert, M.-F.; Aris, I.M.; Goodrich, J.M.; et al. Evaluation of Pediatric Epigenetic Clocks across Multiple Tissues. Clin. Epigenetics 2023, 15, 142. [Google Scholar] [CrossRef]
  51. Heikkinen, T.; Järvinen, A. The Common Cold. Lancet 2003, 361, 51–59. [Google Scholar] [CrossRef]
  52. Hellgren, J.; Cervin, A.; Nordling, S.; Bergman, A.; Cardell, L.O. Allergic Rhinitis and the Common Cold—High Cost to Society. Allergy 2010, 65, 776–783. [Google Scholar] [CrossRef]
  53. Yu, Y.; Wang, J.Q. Phthalate Exposure and Lung Disease: The Epidemiological Evidences, Plausible Mechanism and Advocacy of Interventions. Rev. Environ. Health 2022. [Google Scholar] [CrossRef]
  54. Zhang, Y.; Mustieles, V.; Sun, Y.; Oulhote, Y.; Wang, Y.-X.; Messerlian, C. Association between Serum Per- and Polyfluoroalkyl Substances Concentrations and Common Cold among Children and Adolescents in the United States. Environ. Int. 2022, 164, 107239. [Google Scholar] [CrossRef]
  55. Taneja, V. Sex Hormones Determine Immune Response. Front. Immunol. 2018, 9, 1931. [Google Scholar] [CrossRef] [PubMed]
  56. Beg, M.A.; Sheikh, I.A. Endocrine Disruption: Structural Interactions of Androgen Receptor against Di(2-Ethylhexyl) Phthalate and Its Metabolites. Toxics 2020, 8, 115. [Google Scholar] [CrossRef]
  57. Borch, J.; Metzdorff, S.B.; Vinggaard, A.M.; Brokken, L.; Dalgaard, M. Mechanisms Underlying the Anti-Androgenic Effects of Diethylhexyl Phthalate in Fetal Rat Testis. Toxicology 2006, 223, 144–155. [Google Scholar] [CrossRef]
  58. Foster, P.M.D. Mode of Action: Impaired Fetal Leydig Cell Function—Effects on Male Reproductive Development Produced by Certain Phthalate Esters. Crit. Rev. Toxicol. 2005, 35, 713–719. [Google Scholar] [CrossRef] [PubMed]
  59. Parks, L.G.; Ostby, J.S.; Lambright, C.R.; Abbott, B.D.; Klinefelter, G.R.; Barlow, N.J.; Gray, L.E. The Plasticizer Diethylhexyl Phthalate Induces Malformations by Decreasing Fetal Testosterone Synthesis during Sexual Differentiation in the Male Rat. Toxicol. Sci. 2000, 58, 339–349. [Google Scholar] [CrossRef]
  60. Sathyanarayana, S.; Barrett, E.; Butts, S.; Wang, C.; Swan, S.H. Phthalate Exposure and Reproductive Hormone Concentrations in Pregnancy. Reproduction 2014, 147, 401–409. [Google Scholar] [CrossRef] [PubMed]
  61. Lin, L.-C.; Wang, S.-L.; Chang, Y.-C.; Huang, P.-C.; Cheng, J.-T.; Su, P.-H.; Liao, P.-C. Associations between Maternal Phthalate Exposure and Cord Sex Hormones in Human Infants. Chemosphere 2011, 83, 1192–1199. [Google Scholar] [CrossRef]
  62. Guerra-Silveira, F.; Abad-Franch, F. Sex Bias in Infectious Disease Epidemiology: Patterns and Processes. PLoS ONE 2013, 8, e62390. [Google Scholar] [CrossRef]
  63. Ait Bamai, Y.; Miyashita, C.; Araki, A.; Nakajima, T.; Sasaki, S.; Kishi, R. Effects of Prenatal Di(2-Ethylhexyl) Phthalate Exposure on Childhood Allergies and Infectious Diseases: The Hokkaido Study on Environment and Children’s Health. Sci. Total Environ. 2018, 618, 1408–1415. [Google Scholar] [CrossRef]
  64. Huang, Y.; Garcia, J.M.; Shu, W.; Rong, H.; Zhang, L.; Wang, Y.; Tan, Y.; Lin, H.; Zeng, H.; Chen, J. Peroxisome Proliferator Activated Receptor Gamma in Human Placenta May Mediate the Adverse Effects of Phthalates Exposure in Pregnancy. Reprod. Toxicol. 2018, 75, 121–126. [Google Scholar] [CrossRef]
  65. Schlezinger, J.J.; Howard, G.J.; Hurst, C.H.; Emberley, J.K.; Waxman, D.J.; Webster, T.; Sherr, D.H. Environmental and Endogenous Peroxisome Proliferator-Activated Receptor γ Agonists Induce Bone Marrow B Cell Growth Arrest and Apoptosis: Interactions between Mono(2-Ethylhexyl)Phthalate, 9-Cis-Retinoic Acid, and 15-Deoxy-Δ12,14-Prostaglandin J21. J. Immunol. 2004, 173, 3165–3177. [Google Scholar] [CrossRef]
  66. Xu, M.; Li, Y.; Wang, X.; Zhang, Q.; Wang, L.; Zhang, X.; Cui, W.; Han, X.; Ma, N.; Li, H.; et al. Role of Hepatocyte- and Macrophage-Specific PPARγ in Hepatotoxicity Induced by Diethylhexyl Phthalate in Mice. Environ. Health Perspect. 2022, 130, 017005. [Google Scholar] [CrossRef]
  67. Palacios-Arreola, M.I.; Morales-Montor, J.; Cazares-Martinez, C.J.; Gomez-Arroyo, S.; Nava-Castro, K.E. Environmental Pollutants: An Immunoendocrine Perspective on Phthalates. Front. Biosci.-Landmark 2020, 26, 401–430. [Google Scholar] [CrossRef]
  68. Shepherd, R.; Cheung, A.S.; Pang, K.; Saffery, R.; Novakovic, B. Sexual Dimorphism in Innate Immunity: The Role of Sex Hormones and Epigenetics. Front. Immunol. 2020, 11, 604000. [Google Scholar] [CrossRef]
  69. AbdelHamid, S.G.; Refaat, A.A.; Benjamin, A.M.; Elmawardy, L.A.; Elgendy, L.A.; Manolly, M.M.; Elmaksoud, N.A.; Sherif, N.; Hamdy, N.M. Deciphering Epigenetic(s) Role in Modulating Susceptibility to and Severity of COVID-19 Infection and/or Outcome: A Systematic Rapid Review. Environ. Sci. Pollut. Res. 2021, 28, 54209–54221. [Google Scholar] [CrossRef] [PubMed]
  70. Crimi, E.; Benincasa, G.; Figueroa-Marrero, N.; Galdiero, M.; Napoli, C. Epigenetic Susceptibility to Severe Respiratory Viral Infections and Its Therapeutic Implications: A Narrative Review. Br. J. Anaesth. 2020, 125, 1002–1017. [Google Scholar] [CrossRef] [PubMed]
  71. Wang, J.; Zhou, W.-H. Epigenetic Clocks in the Pediatric Population: When and Why They Tick? Chin. Med. J. 2021, 134, 2901–2910. [Google Scholar] [CrossRef] [PubMed]
  72. Mansell, T.; Saffery, R. The End of the Beginning: Epigenetic Variation in Utero as a Mediator of Later Human Health and Disease. Epigenomics 2017, 9, 217–222. [Google Scholar] [CrossRef] [PubMed]
  73. Faul, J.D.; Kim, J.K.; Levine, M.E.; Thyagarajan, B.; Weir, D.R.; Crimmins, E.M. Epigenetic-Based Age Acceleration in a Representative Sample of Older Americans: Associations with Aging-Related Morbidity and Mortality. Proc. Natl. Acad. Sci. USA 2023, 120, e2215840120. [Google Scholar] [CrossRef] [PubMed]
  74. Yusupov, N.; Dieckmann, L.; Erhart, M.; Sauer, S.; Rex-Haffner, M.; Kopf-Beck, J.; Brückl, T.M.; Czamara, D.; Binder, E.B. Transdiagnostic Evaluation of Epigenetic Age Acceleration and Burden of Psychiatric Disorders. Neuropsychopharmacology 2023, 48, 1409–1417. [Google Scholar] [CrossRef]
  75. Roubinov, D.; Meaney, M.J.; Boyce, W.T. Change of Pace: How Developmental Tempo Varies to Accommodate Failed Provision of Early Needs. Neurosci. Biobehav. Rev. 2021, 131, 120–134. [Google Scholar] [CrossRef]
  76. England-Mason, G.; Merrill, S.M.; Gladish, N.; Moore, S.R.; Giesbrecht, G.F.; Letourneau, N.; MacIsaac, J.L.; MacDonald, A.M.; Kinniburgh, D.W.; Ponsonby, A.-L.; et al. Prenatal Exposure to Phthalates and Peripheral Blood and Buccal Epithelial DNA Methylation in Infants: An Epigenome-Wide Association Study. Environ. Int. 2022, 163, 107183. [Google Scholar] [CrossRef]
  77. Engelbrecht, H.-R.; Merrill, S.M.; Gladish, N.; MacIsaac, J.L.; Lin, D.T.S.; Ecker, S.; Chrysohoou, C.A.; Pes, G.M.; Kobor, M.S.; Rehkopf, D.H. Sex Differences in Epigenetic Age in Mediterranean High Longevity Regions. Front. Aging 2022, 3, 1007098. [Google Scholar] [CrossRef] [PubMed]
  78. Kankaanpää, A.; Tolvanen, A.; Saikkonen, P.; Heikkinen, A.; Laakkonen, E.K.; Kaprio, J.; Ollikainen, M.; Sillanpää, E. Do Epigenetic Clocks Provide Explanations for Sex Differences in Life Span? A Cross-Sectional Twin Study. J. Gerontol. A Biol. Sci. Med. Sci. 2022, 77, 1898–1906. [Google Scholar] [CrossRef]
  79. Haugen, A.C.; Schug, T.T.; Collman, G.; Heindel, J.J. Evolution of DOHaD: The Impact of Environmental Health Sciences. J. Dev. Orig. Health Dis. 2015, 6, 55–64. [Google Scholar] [CrossRef] [PubMed]
  80. Tran, N.Q.V.; Miyake, K. Neurodevelopmental Disorders and Environmental Toxicants: Epigenetics as an Underlying Mechanism. Int. J. Genom. 2017, 2017, 7526592. [Google Scholar] [CrossRef]
  81. Bulka, C.M.; Enggasser, A.E.; Fry, R.C. Epigenetics at the Intersection of COVID-19 Risk and Environmental Chemical Exposures. Curr. Environ. Health Rep. 2022, 9, 477–489. [Google Scholar] [CrossRef]
  82. Laffont, S.; Guéry, J.-C. Deconstructing the Sex Bias in Allergy and Autoimmunity: From Sex Hormones and Beyond. Adv. Immunol. 2019, 142, 35–64. [Google Scholar] [CrossRef]
  83. Zazara, D.E.; Arck, P.C. Developmental Origin and Sex-Specific Risk for Infections and Immune Diseases Later in Life. Semin. Immunopathol. 2019, 41, 137–151. [Google Scholar] [CrossRef]
  84. Dietert, R.R. Developmental Immunotoxicity, Perinatal Programming, and Noncommunicable Diseases: Focus on Human Studies. Adv. Med. 2014, 2014, 867805. [Google Scholar] [CrossRef]
  85. Lucaccioni, L.; Trevisani, V.; Passini, E.; Righi, B.; Plessi, C.; Predieri, B.; Iughetti, L. Perinatal Exposure to Phthalates: From Endocrine to Neurodevelopment Effects. Int. J. Mol. Sci. 2021, 22, 4063. [Google Scholar] [CrossRef] [PubMed]
  86. Main, K.M.; Mortensen, G.K.; Kaleva, M.M.; Boisen, K.A.; Damgaard, I.N.; Chellakooty, M.; Schmidt, I.M.; Suomi, A.-M.; Virtanen, H.E.; Petersen, J.H.; et al. Human Breast Milk Contamination with Phthalates and Alterations of Endogenous Reproductive Hormones in Infants Three Months of Age. Environ. Health Perspect. 2006, 114, 270–276. [Google Scholar] [CrossRef] [PubMed]
  87. Bermejo-Haro, M.Y.; Camacho-Pacheco, R.T.; Brito-Pérez, Y.; Mancilla-Herrera, I. The Hormonal Physiology of Immune Components in Breast Milk and Their Impact on the Infant Immune Response. Mol. Cell. Endocrinol. 2023, 572, 111956. [Google Scholar] [CrossRef] [PubMed]
  88. Kaplan, B.J.; Giesbrecht, G.F.; Leung, B.M.Y.; Field, C.J.; Dewey, D.; Bell, R.C.; Manca, D.P.; O’Beirne, M.; Johnston, D.W.; Pop, V.J.; et al. The Alberta Pregnancy Outcomes and Nutrition (APrON) Cohort Study: Rationale and Methods. Matern. Child. Nutr. 2014, 10, 44–60. [Google Scholar] [CrossRef] [PubMed]
  89. Letourneau, N.; Aghajafari, F.; Bell, R.C.; Deane, A.J.; Dewey, D.; Field, C.; Giesbrecht, G.; Kaplan, B.; Leung, B.; Ntanda, H. The Alberta Pregnancy Outcomes and Nutrition (APrON) Longitudinal Study: Cohort Profile and Key Findings from the First Three Years. BMJ Open 2022, 12, e047503. [Google Scholar] [CrossRef] [PubMed]
  90. England-Mason, G.; Grohs, M.N.; Reynolds, J.E.; MacDonald, A.; Kinniburgh, D.; Liu, J.; Martin, J.W.; Lebel, C.; Dewey, D.; APrON Study Team. White Matter Microstructure Mediates the Association between Prenatal Exposure to Phthalates and Behavior Problems in Preschool Children. Environ. Res. 2020, 182, 109093. [Google Scholar] [CrossRef] [PubMed]
  91. England-Mason, G.; Martin, J.W.; MacDonald, A.; Kinniburgh, D.; Giesbrecht, G.F.; Letourneau, N.; Dewey, D. Similar Names, Different Results: Consistency of the Associations between Prenatal Exposure to Phthalates and Parent-Ratings of Behavior Problems in Preschool Children. Environ. Int. 2020, 142, 105892. [Google Scholar] [CrossRef]
  92. Hornung, R.W.; Reed, L.D. Estimation of Average Concentration in the Presence of Nondetectable Values. Appl. Occup. Environ. Hyg. 1990, 5, 46–51. [Google Scholar] [CrossRef]
  93. Messerlian, C.; Wylie, B.J.; Minguez-Alarcon, L.; Williams, P.L.; Ford, J.B.; Souter, I.C.; Calafat, A.M.; Hauser, R. Urinary Concentrations of Phthalate Metabolites in Relation to Pregnancy Loss among Women Conceiving with Medically Assisted Reproduction. Epidemiology 2016, 27, 879–888. [Google Scholar] [CrossRef] [PubMed]
  94. Letourneau, N.L.; Kozyrskyj, A.L.; Cosic, N.; Ntanda, H.N.; Anis, L.; Hart, M.J.; Campbell, T.S.; Giesbrecht, G.F. Maternal Sensitivity and Social Support Protect against Childhood Atopic Dermatitis. Allergy Asthma Clin. Immunol. 2017, 13, 26. [Google Scholar] [CrossRef] [PubMed]
  95. Bibikova, M.; Barnes, B.; Tsan, C.; Ho, V.; Klotzle, B.; Le, J.M.; Delano, D.; Zhang, L.; Schroth, G.P.; Gunderson, K.L.; et al. High Density DNA Methylation Array with Single CpG Site Resolution. Genomics 2011, 98, 288–295. [Google Scholar] [CrossRef] [PubMed]
  96. Houseman, E.A.; Accomando, W.P.; Koestler, D.C.; Christensen, B.C.; Marsit, C.J.; Nelson, H.H.; Wiencke, J.K.; Kelsey, K.T. DNA Methylation Arrays as Surrogate Measures of Cell Mixture Distribution. BMC Bioinform. 2012, 13, 86. [Google Scholar] [CrossRef] [PubMed]
  97. Gervin, K.; Salas, L.A.; Bakulski, K.M.; Van Zelm, M.C.; Koestler, D.C.; Wiencke, J.K.; Duijts, L.; Moll, H.A.; Kelsey, K.T.; Kobor, M.S.; et al. Systematic Evaluation and Validation of Reference and Library Selection Methods for Deconvolution of Cord Blood DNA Methylation Data. Clin. Epigenetics 2019, 11, 125. [Google Scholar] [CrossRef]
  98. Koestler, D.C.; Jones, M.J.; Usset, J.; Christensen, B.C.; Butler, R.A.; Kobor, M.S.; Wiencke, J.K.; Kelsey, K.T. Improving Cell Mixture Deconvolution by Identifying Optimal DNA Methylation Libraries (IDOL). BMC Bioinform. 2016, 17, 120. [Google Scholar] [CrossRef]
  99. Salas, L.A.; Zhang, Z.; Koestler, D.C.; Butler, R.A.; Hansen, H.M.; Molinaro, A.M.; Wiencke, J.K.; Kelsey, K.T.; Christensen, B.C. Enhanced Cell Deconvolution of Peripheral Blood Using DNA Methylation for High-Resolution Immune Profiling. Nat. Commun. 2022, 13, 761. [Google Scholar] [CrossRef]
  100. Merrill, S.M.; Moore, S.R.; Gladish, N.; Giesbrecht, G.F.; Dewey, D.; Konwar, C.; MacIssac, J.L.; Kobor, M.S.; Letourneau, N.L. Paternal Adverse Childhood Experiences: Associations with Infant DNA Methylation. Dev. Psychobiol. 2021, 63, e22174. [Google Scholar] [CrossRef]
  101. Hermansen, M.C. Nucleated Red Blood Cells in the Fetus and Newborn. Arch. Dis. Child. Fetal Neonatal Ed. 2001, 84, F211–F215. [Google Scholar] [CrossRef]
  102. Bell, C.G.; Lowe, R.; Adams, P.D.; Baccarelli, A.A.; Beck, S.; Bell, J.T.; Christensen, B.C.; Gladyshev, V.N.; Heijmans, B.T.; Horvath, S.; et al. DNA Methylation Aging Clocks: Challenges and Recommendations. Genome Biol. 2019, 20, 249. [Google Scholar] [CrossRef]
  103. Moore, S.R.; McEwen, L.M.; Quirt, J.; Morin, A.; Mah, S.M.; Barr, R.G.; Boyce, W.T.; Kobor, M.S. Epigenetic Correlates of Neonatal Contact in Humans. Dev. Psychopathol. 2017, 29, 1517–1538. [Google Scholar] [CrossRef]
  104. Horvath, S.; Oshima, J.; Martin, G.M.; Lu, A.T.; Quach, A.; Cohen, H.; Felton, S.; Matsuyama, M.; Lowe, D.; Kabacik, S.; et al. Epigenetic Clock for Skin and Blood Cells Applied to Hutchinson Gilford Progeria Syndrome and Ex Vivo Studies. Aging 2018, 10, 1758–1775. [Google Scholar] [CrossRef]
  105. Impinen, A.; Nygaard, U.C.; Lødrup Carlsen, K.C.; Mowinckel, P.; Carlsen, K.H.; Haug, L.S.; Granum, B. Prenatal Exposure to Perfluoralkyl Substances (PFASs) Associated with Respiratory Tract Infections but Not Allergy- and Asthma-Related Health Outcomes in Childhood. Environ. Res. 2018, 160, 518–523. [Google Scholar] [CrossRef]
  106. Lourenço, V.M.; Pires, A.M.; Kirst, M. Robust Linear Regression Methods in Association Studies. Bioinformatics 2011, 27, 815–821. [Google Scholar] [CrossRef] [PubMed]
  107. Venables; Ripley, B.D. Modern Applied Statistics with S, 4th ed.; Springer: Berlin/Heidelberg, Germany, 2002. [Google Scholar]
  108. Greenland, S.; Daniel, R.; Pearce, N. Outcome Modelling Strategies in Epidemiology: Traditional Methods and Basic Alternatives. Int. J. Epidemiol. 2016, 45, 565–575. [Google Scholar] [CrossRef] [PubMed]
  109. Zu, J.; Yuan, K.-H. Local Influence and Robust Procedures for Mediation Analysis. Multivar. Behav. Res. 2010, 45, 1–44. [Google Scholar] [CrossRef] [PubMed]
  110. Hayes, A.F. Introduction to Mediation, Moderation, and Conditional Process Analysis: A Regression-Based Approach; Guilford Publications: New York, NY, USA, 2022; ISBN 978-1-4625-4903-0. [Google Scholar]
  111. MacKinnon, D.P.; Fairchild, A.J.; Fritz, M.S. Mediation Analysis. Annu. Rev. Psychol. 2007, 58, 593. [Google Scholar] [CrossRef] [PubMed]
  112. Benjamini, Y.; Hochberg, Y. Controlling the False Discovery Rate: A Practical and Powerful Approach to Multiple Testing. J. R. Stat. Soc. Ser. B 1995, 57, 289–300. [Google Scholar] [CrossRef]
  113. Fritz, M.S.; Mackinnon, D.P. Required Sample Size to Detect the Mediated Effect. Psychol. Sci. 2007, 18, 233–239. [Google Scholar] [CrossRef] [PubMed]
Epigenomes 08 00003 g001
Figure 1. Mediation model in which infant EAA was examined as a mediator of prenatal exposure to DEHP and the number of early childhood colds in the overall sample (n = 69). This model depicts the relation of prenatal DEHP exposure to early childhood colds (path c; total effect), the relation of prenatal DEHP exposure to EAA (path a), the relation of EAA to early childhood colds adjusted for prenatal DEHP exposure (path b), and the relation of prenatal DEHP exposure to early childhood colds adjusted for EAA (path c’; direct effect). Unstandardized regression coefficients and 95% CIs are reported. ** adjusted p < 0.05. * adjusted p < 0.10.
Figure 1. Mediation model in which infant EAA was examined as a mediator of prenatal exposure to DEHP and the number of early childhood colds in the overall sample (n = 69). This model depicts the relation of prenatal DEHP exposure to early childhood colds (path c; total effect), the relation of prenatal DEHP exposure to EAA (path a), the relation of EAA to early childhood colds adjusted for prenatal DEHP exposure (path b), and the relation of prenatal DEHP exposure to early childhood colds adjusted for EAA (path c’; direct effect). Unstandardized regression coefficients and 95% CIs are reported. ** adjusted p < 0.05. * adjusted p < 0.10.
Epigenomes 08 00003 g001
Epigenomes 08 00003 g002
Figure 2. Mediation model in which infant EAA was examined as a mediator of prenatal exposure to DEHP and the number of early childhood colds in girls only (n = 34). For this sub-sample, this model depicts the relation of prenatal DEHP exposure to early childhood colds (path c; total effect), the relation of prenatal DEHP exposure to EAA (path a), the relation of EAA to early childhood colds adjusted for prenatal DEHP exposure (path b), and the relation of prenatal DEHP exposure to early childhood colds adjusted for EAA (path c’; direct effect). Unstandardized regression coefficients and 95% CIs are reported. ** adjusted p < 0.05.
Figure 2. Mediation model in which infant EAA was examined as a mediator of prenatal exposure to DEHP and the number of early childhood colds in girls only (n = 34). For this sub-sample, this model depicts the relation of prenatal DEHP exposure to early childhood colds (path c; total effect), the relation of prenatal DEHP exposure to EAA (path a), the relation of EAA to early childhood colds adjusted for prenatal DEHP exposure (path b), and the relation of prenatal DEHP exposure to early childhood colds adjusted for EAA (path c’; direct effect). Unstandardized regression coefficients and 95% CIs are reported. ** adjusted p < 0.05.
Epigenomes 08 00003 g002
Epigenomes 08 00003 g003
Figure 3. Mediation model in which infant EAA was examined as a mediator of prenatal exposure to DEHP and number of early childhood URIs colds in boys only (n = 35). For this sub-sample, this model depicts the relation of prenatal DEHP exposure to early childhood colds (path c; total effect), the relation of prenatal DEHP exposure to EAA (path a), the relation of EAA to early childhood colds adjusted for prenatal DEHP exposure (path b), and the relation of prenatal DEHP exposure to early childhood colds adjusted for EAA (path c’; direct effect). Unstandardized regression coefficients and 95% CIs are reported. ** adjusted p < 0.05. * adjusted p < 0.10.
Figure 3. Mediation model in which infant EAA was examined as a mediator of prenatal exposure to DEHP and number of early childhood URIs colds in boys only (n = 35). For this sub-sample, this model depicts the relation of prenatal DEHP exposure to early childhood colds (path c; total effect), the relation of prenatal DEHP exposure to EAA (path a), the relation of EAA to early childhood colds adjusted for prenatal DEHP exposure (path b), and the relation of prenatal DEHP exposure to early childhood colds adjusted for EAA (path c’; direct effect). Unstandardized regression coefficients and 95% CIs are reported. ** adjusted p < 0.05. * adjusted p < 0.10.
Epigenomes 08 00003 g003
Table 1. Maternal and child characteristics for the overall sample (n = 69).
Table 1. Maternal and child characteristics for the overall sample (n = 69).
n (%)Mean (SD)
Maternal Characteristics
Age (years)-31.24 (4.01)
White64 (92.75%)-
Married/Cohabiting69 (100.00%)-
Household Income > CAD 70k 156 (81.16%)-
Child Characteristics
Birthweight (g)-3479.42 (505.91)
Gestational age at birth (weeks)-39.49 (1.40)
Sex (Female)34 (49.28%)-
Age at blood draw (weeks)-12.52 (0.96)
Number of colds-2.54 (1.51)
1 CAD (Canadian dollars).
Table 2. Creatinine-adjusted phthalate metabolite concentrations (μg/g creatinine) and DEHP (μmol/g creatinine) in maternal second trimester urines (n = 69).
Table 2. Creatinine-adjusted phthalate metabolite concentrations (μg/g creatinine) and DEHP (μmol/g creatinine) in maternal second trimester urines (n = 69).
Metabolite % > LODMinimumMaximumGM25th Percentile50th Percentile75th Percentile
MEHP 1100%0.71231.03.522.263.124.95
MEHHP 1100%2.3233.510.67.059.8116.4
MEOHP 1100%2.5825.79.136.538.8313.0
MECPP 1100%6.4162.416.611.916.023.3
DEHP 2-0.02220.2850.08140.05340.07550.126
1 LOD = 0.10 μg/L.; 2 Molar sum of individual DEHP metabolites.
Table 3. Post hoc analyses which considered intrinsic epigenetic age acceleration (IEAA) and EAA estimated using Horvath’s skin and blood clock.
Table 3. Post hoc analyses which considered intrinsic epigenetic age acceleration (IEAA) and EAA estimated using Horvath’s skin and blood clock.
IEAAEAA (Skin and Blood Clock)
B (95% CI)B (95% CI)
Overall Sample (n = 69)
Total Effect (path c) 1--
Path a−0.34 (−1.30, 0.61)−0.03 (−0.58, 0.51)
Path b2.98 † (1.67, 4.28)4.15 † (1.38, 6.92)
Direct effect (path c’)7.81 † (2.54, 13.07)7.11 † (1.63, 12.58)
Indirect effect (path ab)−1.79 (−5.41, 1.26)−0.07 (−2.34, 2.35)
Girls (n = 34)
Total Effect (path c) 1--
Path a−2.40 † (−3.54, −1.27)−0.80 † (−1.59, −0.02)
Path b1.04 (−1.88, 3.96)3.47 (−0.79, 7.73)
Direct effect (path c’)−3.47 (−15.03, 8.09)−2.24 (−11.86, 7.39)
Indirect effect (path ab)−2.11 (−8.56, 5.87)−3.90 † (−10.56, −0.23)
Boys (n = 35)
Total Effect (path c) 1--
Path a0.85 (−0.54, 2.25)0.50 (−0.16, 1.17)
Path b3.70 † (1.74, 5.66)4.55 † (0.28, 8.81)
Direct effect (path c’)7.29 (−0.12, 14.70)8.38 (−0.11, 16.87)
Indirect effect (path ab)3.27 (−0.78, 8.58)3.15 † (0.04, 9.52)
† Unadjusted p < 0.05. 1 Path not examined again as it is already reported in the results (Section 2.5 and Section 2.6).
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Merrill, S.M.; Letourneau, N.; Giesbrecht, G.F.; Edwards, K.; MacIsaac, J.L.; Martin, J.W.; MacDonald, A.M.; Kinniburgh, D.W.; Kobor, M.S.; Dewey, D.; et al. Sex-Specific Associations between Prenatal Exposure to Di(2-ethylhexyl) Phthalate, Epigenetic Age Acceleration, and Susceptibility to Early Childhood Upper Respiratory Infections. Epigenomes 2024, 8, 3. https://doi.org/10.3390/epigenomes8010003

AMA Style

Merrill SM, Letourneau N, Giesbrecht GF, Edwards K, MacIsaac JL, Martin JW, MacDonald AM, Kinniburgh DW, Kobor MS, Dewey D, et al. Sex-Specific Associations between Prenatal Exposure to Di(2-ethylhexyl) Phthalate, Epigenetic Age Acceleration, and Susceptibility to Early Childhood Upper Respiratory Infections. Epigenomes. 2024; 8(1):3. https://doi.org/10.3390/epigenomes8010003

Chicago/Turabian Style

Merrill, Sarah M., Nicole Letourneau, Gerald F. Giesbrecht, Karlie Edwards, Julia L. MacIsaac, Jonathan W. Martin, Amy M. MacDonald, David W. Kinniburgh, Michael S. Kobor, Deborah Dewey, and et al. 2024. "Sex-Specific Associations between Prenatal Exposure to Di(2-ethylhexyl) Phthalate, Epigenetic Age Acceleration, and Susceptibility to Early Childhood Upper Respiratory Infections" Epigenomes 8, no. 1: 3. https://doi.org/10.3390/epigenomes8010003

APA Style

Merrill, S. M., Letourneau, N., Giesbrecht, G. F., Edwards, K., MacIsaac, J. L., Martin, J. W., MacDonald, A. M., Kinniburgh, D. W., Kobor, M. S., Dewey, D., England-Mason, G., & The APrON Study Team. (2024). Sex-Specific Associations between Prenatal Exposure to Di(2-ethylhexyl) Phthalate, Epigenetic Age Acceleration, and Susceptibility to Early Childhood Upper Respiratory Infections. Epigenomes, 8(1), 3. https://doi.org/10.3390/epigenomes8010003

Article Metrics

Back to TopTop