Abstract
Maternal exposure to green space during pregnancy has been associated with lower risk for adverse birth outcomes; however, the underlying biological mechanisms remain largely unknown. Epigenetic changes, such as DNA methylation (DNAm), might be one of the molecular mechanisms contributing to this association. The placenta is the key organ for foetal growth and development, but studies of prenatal green space exposure and placental DNAm are scarce and have not been conducted on a genome-wide scale. Here, we aimed to investigate the association between exposure to green space during pregnancy and epigenome-wide placental DNAm.
This study was conducted in 550 mother-child pairs from the Barcelona Life Study Cohort (BiSC), a population-based birth cohort in Spain (2018-2021). We comprehensively assessed green space exposure during pregnancy as (i) residential surrounding greenness (satellite-based Normalized Difference Vegetation Index (NDVI) in buffers of 100m, 300m and 500m), (ii) residential distance to the nearest major green space (meters), iii) use of green space (hours/week), and iv) visual access to greenery through home window (≥half of the view). Placental DNAm was measured with the EPIC array. To identify differentially methylated positions (DMPs), we fitted robust linear regression models adjusted for relevant covariates, while differentially methylated regions (DMRs) were identified using the dmrff method.
After Bonferroni (BN) correction, cg14852540, annotated to SLC25A10 gene, was inversely associated with residential surrounding greenness within 500m buffer. Additionally, ten unique DMRs were associated with surrounding greenness within 300m and 500m, distance or visual access to green space. No associations were found for surrounding greenness within 100m or use of green space. Genes annotated to these loci are involved in transcriptional regulation, metabolism and mitochondrial respiration.
Overall, we identified associations between prenatal green space exposure and DNAm levels in placenta. Further research is needed to validate these results and understand the underlying biological pathways.
Competing Interest Statement
The authors have declared no competing interest.
Funding Statement
NutriPROGRAM (AC18/00006) ALMA project (PI20/00190) ATHLETE project (874583) AirNB project (785994) Health Effects Institute (HEI)(R-82811201) ENTENTE project (PI20/01116) FI-AGAUR Predoctoral contract (2023 FI-2 00797) Beatriu de Pinos Postdoctoral Contract awarded by Generalitat de Catalunya-AGAUR and European Commission-Horizon 2020 (2019 BP 00107). Miguel Servet fellowship (Grant CP19/00183) We acknowledge support from the grant CEX2023-0001290-S funded by MCIN/AEI/ 10.13039/501100011033 and support from the Generalitat de Catalunya through the CERCA Program.
Author Declarations
I confirm all relevant ethical guidelines have been followed, and any necessary IRB and/or ethics committee approvals have been obtained.
Yes
The details of the IRB/oversight body that provided approval or exemption for the research described are given below:
Clinical Research Ethics Committee of the Parc de Salut Mar project (2018/8050/I) Medical Research Committee of the Fundacio de Gestio Sanitaria del Hospital de la Santa Creu i Sant Pau de Barcelona (EC/18/206/5272) Ethics Committee of the Fundacio Sant Joan de Deu (PIC-27-18).
I confirm that all necessary patient/participant consent has been obtained and the appropriate institutional forms have been archived, and that any patient/participant/sample identifiers included were not known to anyone (e.g., hospital staff, patients or participants themselves) outside the research group so cannot be used to identify individuals.
Yes
I understand that all clinical trials and any other prospective interventional studies must be registered with an ICMJE-approved registry, such as ClinicalTrials.gov. I confirm that any such study reported in the manuscript has been registered and the trial registration ID is provided (note: if posting a prospective study registered retrospectively, please provide a statement in the trial ID field explaining why the study was not registered in advance).
Yes
I have followed all appropriate research reporting guidelines, such as any relevant EQUATOR Network research reporting checklist(s) and other pertinent material, if applicable.
Yes
Footnotes
↵* Joint last authors
Data Availability
Genome wide DNA methylation summarized results can be found at Zenodo repository (THIS WILL BE DONE UPON ACCEPTANCE OF THE MANUSCRIPT). Individual cohort level data may be available by application to the relevant institution after obtaining required approvals.
Abbreviations
- ADCY4
- Adenylate Cyclase 4 gene
- AIP
- Aryl Hydrocarbon Receptor Interacting Protein
- BiSC
- Barcelona Life Study Cohort
- BMIQ
- Beta-mixture quantile
- BN
- Bonferroni
- Bp
- base pairs
- C3orf55
- Solute Carrier Family 66 Member 1 Like, Pseudogene
- CpG
- Cytosine-phosphate-guanine
- COPD
- chronic obstructive pulmonary disease
- DMP
- Differentially methylated position
- DMR
- Differentially methylated region
- DNAm
- DNA methylation
- IQR
- Interquartile range
- eQTM
- Expression quantitative methylated regions
- ESCAPE
- European Study of Cohorts for Air Pollution Effects
- Enh
- Enhancer
- EWAS
- Epigenome-wide association study
- gDMR
- Germline differently methylated region
- GLRB
- Glycine Receptor Beta gene
- GPS
- Global Positioning System
- GS
- Green space
- GO
- Gene Ontology
- HSPA1B
- Heat Shock Protein Family A (Hsp70) Member 1B gene
- HTR2A
- 5-Hydroxytryptamine Receptor 2A gene
- HUGE-F
- Human Genome facility
- ICGC
- Cartographic and Geology Institute of Catalonia
- KEGG
- Kyoto Encyclopedia of Genes and Genomes
- LOX
- Lysile Oxidase gene
- LUR
- Land Use Regression Models
- mQTL
- methylation Quantitative Trait Loci
- mtDNAc
- mitochondria DNA content
- NR2E1
- Nuclear Receptor Subfamily 2 Group E Member 1 2
- OSBPL9
- Oxysterol Binding Protein Like 9 gene
- PCA
- Principal Component Analysis
- PDE8A
- Phosphodiesterase 8A gene
- PM2.5
- Particulate matter with aerodynamic diameter 2.5 microm
- PMD
- Partially methylated domain
- QC
- Quality control
- R
- Red values
- MAF
- minor allele frequency
- NIR
- near infrared values
- NDVI
- Normalized Difference Vegetation Index
- SKA1
- Spindle And Kinetochore Associated Complex Subunit 1 gene
- SLC25A10
- Solute carrier family 25 member 10 gene
- SNPs
- single nucleotide polymorphisms
- SES
- socioeconomic status
- TSS
- transcription start site
- TssA
- active transcription start site
- TssAFlink
- flanking active transcription start site
- UCSC
- University of California, Santa Cruz
- UTR
- untranslated region
- ZNF22
- Zinc Finger Protein 22
- ZNF22-AS1
- Zinc Finger Protein 22 Antisense RNA 1
- ZNF562
- Zinc Finger Protein 562 gene.