Study population and design
This study was conducted in the framework of IPANEMA (Impact of Particulate Matter on Mothers and Babies in Antwerp), a prospective pregnancy cohort study of the Antwerp University Hospital (UZA) in collaboration with the Flemish Institute for Technological Research (VITO) and the University of Antwerp (UA). Pregnant women were recruited between April 2015 and January 2018 at the UZA prenatal clinic by a midwife or obstetrician at a gestational age of 12 to 14 weeks. The inclusion criteria were: a singleton pregnancy; the ability to fill out extensive Dutch questionnaires; delivery planned in the Antwerp University Hospital. All participating mothers gave written informed consent. The study protocol was approved by the ethical committee of the University of Antwerp (14/40/411) and registered under number NCT02592005 at clinicaltrials.gov. Health-related information on mothers and babies was extracted from the hospital records and questionnaires that participants completed at enrolment, during pregnancy and after delivery. These questionnaires provided detailed information on participants’ socio-demographic and lifestyle characteristics. A detailed protocol of the IPANEMA study can be found elsewhere [23].
Hair sample collection and cortisol measurement
At the end of second trimester of pregnancy in-hospital consultation and shortly after delivery, a strand of hair of at least 2mm thick was bound together with a cotton thread and cut close to the scalp from the posterior vertex region of the head. This area of the scalp exhibits the lowest intra-individual variability in HCC [24]. Hair samples were stored in paper envelopes at room temperature until analysis. When protected from ultraviolet light, cortisol concentrations in hair samples remain stable at room temperature for several years [25]. Cortisol concentrations were determined from the 3 cm of hair closest to the scalp. Based on an average hair growth of 1 cm per month, this length represents cortisol secretion in a 3-month period, a trimester, prior to sampling [26]. There is a wide consensus that the first 5–6 cm of hair nearest to a person's scalp can reliably reflect HPA activity [27]. Analysis was performed at the Institute of Public Health, Department of Environmental Medicine of the University of Southern Denmark (SDU), using liquid chromatography and tandem mass spectrometry (LC-MS/MS) as described by Chen et al. [28], after minor modifications. Hair samples were washed with methanol and dried at room temperature. The 3 cm of hair closest to the scalp was cut into segments of 2-3 mm. A typical amount of hair weighed 20-30 mg. Aliquots of 100 mL 20 ng/mL isotope labeled cortisol (cortisol-D4) were added as internal standard, together with 0.9 mL methanol. Samples were incubated in the dark at 25°C while whirl mixed at 2000 revolutions per minute for 5 days and subsequently centrifuged at 3000g for 5 minutes. 20 μL of the supernatant was injected onto a High-Performance Liquid Chromatography (HPLC) column. HPLC was performed using an Accella 1250 pump (Thermo Scientific, San Jose, CA) and a PAL autosampler (CTC analytics, Zwingen, Switzerland). The analytical column was a Kinetex C18 column, 100 x 4.6 mm (2.6 μm) equipped with a 2 x 4 mm C18 SecurityGuard column (Phenomenex, Torrance, CA). Isocratic elution was performed with a mobile phase system consisting of methanol and 0.1 M formic acid (80:20) at a flow rate of 400 μL/min for 6 min. After the peaks were eluted, a wash procedure was performed before the next samples was injected onto the column. The triple quadrupole mass spectrometer utilized was a TSQ Vantage (Thermo Scientific, San Jose, CA). The calibration curve and calculation of the sample concentration were based on the area ratio of the analyte/isotope labeled internal standard. Quality control samples were included in each series of samples. The limit of quantification (LOQ) for cortisol was 1.0 pg/mg hair. The intra-day repeatability coefficient of variation was 8.7% and the inter-day reproducibility coefficient of variation was 9.5%
Residential exposure assessment
Assessment of all residential exposure variables was based on the participants’ geocoded home address. Geographical Information System (GIS) analyses were carried out using ESRI ArcGIS software version 10.4 (Environmental Systems Research Institute, Redlands, California, USA). The residential degree of urbanization was assessed according to the Eurostat definition that classifies local administrative units as cities, towns, suburbs or rural areas based on a combination of geographical contiguity and population density, applied to 1 km² population grid cells [29]. We assessed residential exposure to fine particulate matter (PM2.5), nitrogen dioxide (NO2) and black carbon (BC), primary constituents of traffic-related air pollution. Residential exposure to PM2.5, NO2 and BC was modelled using a spatial temporal interpolation method. In Flanders, atmospheric pollutants are continuously measured by a network of automatic monitoring stations by the Flemish Environment Agency. The Belgian Interregional Environment Agency (IRCEL, Intergewestelijke Cel voor het Leefmilieu) uses these measurements together with information on land cover to interpolate the air pollutant concentrations on a 4x4 km² resolution [30]. These background results are combined with a bi-gaussian dispersion model based on emissions from point sources and line sources, the Immission Frequency Distribution Model (IFDM). The combined RIO-IFDM model chain produces daily averaged pollutant concentrations in Belgium on a high resolution receptor grid [31]. We calculated mean air pollutant concentrations at the residential address over a 3-month period before sampling, similar to the period of cortisol accumulation in the hair samples, and over a 1-year period before sampling. Residential proximity to major roads is often used as a surrogate measure of long-term exposure to traffic-related air pollution [32]. We calculated the straight-line distance from each residence to the nearest major road. Major roads included international motorways (E-roads) and the network of large national and local roads of Belgium (N-roads).
Residential access to a neighborhood greenspace was based on the 2016 version of the land-use map of Flanders, which maps land cover types, i.e. natural vegetated land cover and urban greenery, in 10x10 m² raster cells [33]. Green cells were clustered to assess the area and public accessibility of greenspace in the maternal residential surroundings. Access to a small neighborhood greenspace (NHGS) was defined as access to at least 0.2 hectares (ha) of greenspace within a travel distance of 400 meters (m) from residence, access to a large greenspace was defined as access to at least 10 ha of greenspace within a travel distance of 800 m from residence. In the large greenspace typology, small water bodies are included when surrounded by > 50% greenspace, agricultural land is included when surrounded by > 30% greenspace. More technical background information on the green typology indicators can be found elsewhere [34].
Potential covariates of hair cortisol concentrations
Possible covariates of HCC were identified based on available data within the IPANEMA cohort and on existing literature [35]. Tested covariates included maternal age, parity, maternal socioeconomic status (SES) defined as the highest educational attainment of the mother and categorized as low/intermediate/high, pre-existing chronic diseases (diabetes, asthma, cardiovascular disease), pre-pregnancy body mass index (BMI), gestational week at sampling, season of sampling, smoking and alcohol consumption before pregnancy, systemic use of glucocorticoids and daily hair washing. We assessed maternal ethnic background as European/non-European country of birth since hair growth rate may be influenced by ethnicity [36]. We additionally tested variables that may have a link with both residential environment and biological stress, i.e. neighborhood SES and residential exposure to noise. A systematic review in the World Health Organization (WHO) European Region showed that lower neighborhood SES is usually linked with higher levels of air pollutants [37]. Independent from higher levels of exposure, deprived mothers may have a higher vulnerability, leading to more pronounced adverse health effects of a given environmental exposure [38]. The Area Deprivation Index (ADI) is a yearly calculated indicator for neighborhood SES on a sub-municipality level in Flanders (Statistics Flanders, n.d.). Deprivation is recorded by the Flemish Child and Family Government Agency (www.kindengezin.be). Selection criteria for deprivation are the family’s monthly income, the parents’ educational attainment, the children’s development, the parents’ employment situation, housing and health. If a family fulfils three or more criteria, it is considered to be underprivileged (OECD, 2000). The index of year X (%), i.e. 2017, considers all children born in year X, X-1 and X-2 that live in deprived households in a given area in Flanders, divided by the total number of children born in the area during the same period. The ADI of the participants’ neighborhood was subsequently categorized into tertiles representing low, intermediate and high area deprivation across the range of ADI among all participants.
Residential proximity to major roads may also lead to elevated noise levels [41]. Residential noise exposure levels were assessed using the Flemish strategic noise map of 2016, which includes major road infrastructure as defined in the EU-guideline 2002/49/EG [42]. The strategic noise map expresses noise levels in Lden, the average sound level over a 24 hour period with a penalty of 5 dB added for evening hours and a penalty of 10 dB added for nighttime hours [43]. The WHO guideline for average noise exposure produced by road traffic is set at 53 decibels (dB) Lden, road traffic noise above this level has been associated with adverse health effects, including adverse birth outcomes [44]. Noise exposure was therefore evaluated binary as exposure to a noise level ≥ 53 dB Lden.
Statistical analysis
Statistical analysis was performed using SPSS Statistics (version 25; IBM, Armonk, NY, USA) and R (version 2018; R Foundation for Statistical Computing, Vienna, Austria). Descriptive statistics provide an overview of study population characteristics, residential exposure characteristics and geometric mean HCC concentration with 95% confidence interval. Air pollution variables were logarithmically transformed (ln-scale) because of skewed distributions, distance to major roads was logarithmically transformed to reflect the non-linear distance decay of traffic-related exposure to air pollutants [45].
Spearman rank correlations between residential exposures variables were assessed, since correlations of 0.9 or higher between exposure variables indicate strongly connected exposures that cannot be disentangled [46]. The outcome variable HCC was logarithmically transformed to obtain a normal distribution. For HCC below the LOQ of 1 pg/mg hair, a random imputation from a log-normal probability distribution was performed where the mean was allowed to depend on observed values for hair cortisone concentrations that were measured simultaneously with cortisol, since both glucocorticoids were highly correlated (p < 0.01, Pearson's r = 0.711 for 2nd trimester cortisol and cortisone, p < 0.01, Pearson's r = 0.758 for 3rd trimester cortisol and cortisone). Linear regression models were used to analyze associations between 3-month mean air pollutant concentrations (PM2.5, NO2, BC) , distance to major roads and access to greenspace as a predictor and 2nd and 3rd trimester Hair Cortisol Concentrations as an outcome. Given the limited number of study participants, we decided not to adjust for a set of a priori selected covariates. The final regression models were only adjusted for significant covariates (p < 0.05). All assumptions of linear regression were checked. To quantify the association, the estimated change in HCC (β) with 95% confidence interval (95% CI) is presented for an increase in exposure from the 25th to the 75th percentile.
Effect modification by access the neighborhood greenspace was assessed by adding the interaction term of exposure to air pollution or distance to major roads and access to greenspace into the regression model. The level of significance for estimates was set at p < 0.05.
We conducted several sensitivity analyses to evaluate our results. We tested 1-year mean air pollutant concentrations in relation to 2nd and 3rd trimester HCC to confirm the robustness of 3-month mean results. We additionally adjusted our 2nd trimester models for frequency of hair washing, a significant determinant of 2nd trimester HCC, independent of biological stress.