PT - JOURNAL ARTICLE AU - Aguilar-Lacasaña, Sofía AU - Cosin-Tomas, Marta AU - Raimbault, Bruno AU - Gómez, Laura AU - Sánchez, Olga AU - Zanini, Maria Julia AU - Capdevila, Rosalia Pascal AU - Foraster, Maria AU - Gascon, Mireia AU - Rivas, Ioar AU - Llurba, Elisa AU - Gómez-Roig, Maria Dolores AU - Sunyer, Jordi AU - Vrijheid, Martine AU - Bustamante, Mariona AU - Dadvand, Payam TI - Epigenome-wide association study of pregnancy exposure to green space and placental DNA methylation AID - 10.1101/2024.07.21.24310780 DP - 2024 Jan 01 TA - medRxiv PG - 2024.07.21.24310780 4099 - http://medrxiv.org/content/early/2024/07/22/2024.07.21.24310780.short 4100 - http://medrxiv.org/content/early/2024/07/22/2024.07.21.24310780.full AB - 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 StatementThe authors have declared no competing interest.Funding StatementNutriPROGRAM (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 DeclarationsI confirm all relevant ethical guidelines have been followed, and any necessary IRB and/or ethics committee approvals have been obtained.YesThe 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.YesI 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).YesI have followed all appropriate research reporting guidelines, such as any relevant EQUATOR Network research reporting checklist(s) and other pertinent material, if applicable.YesGenome 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.AbbreviationsADCY4Adenylate Cyclase 4 geneAIPAryl Hydrocarbon Receptor Interacting ProteinBiSCBarcelona Life Study CohortBMIQBeta-mixture quantileBNBonferroniBpbase pairsC3orf55Solute Carrier Family 66 Member 1 Like, PseudogeneCpGCytosine-phosphate-guanineCOPDchronic obstructive pulmonary diseaseDMPDifferentially methylated positionDMRDifferentially methylated regionDNAmDNA methylationIQRInterquartile rangeeQTMExpression quantitative methylated regionsESCAPEEuropean Study of Cohorts for Air Pollution EffectsEnhEnhancerEWASEpigenome-wide association studygDMRGermline differently methylated regionGLRBGlycine Receptor Beta geneGPSGlobal Positioning SystemGSGreen spaceGOGene OntologyHSPA1BHeat Shock Protein Family A (Hsp70) Member 1B geneHTR2A5-Hydroxytryptamine Receptor 2A geneHUGE-FHuman Genome facilityICGCCartographic and Geology Institute of CataloniaKEGGKyoto Encyclopedia of Genes and GenomesLOXLysile Oxidase geneLURLand Use Regression ModelsmQTLmethylation Quantitative Trait LocimtDNAcmitochondria DNA contentNR2E1Nuclear Receptor Subfamily 2 Group E Member 1 2OSBPL9Oxysterol Binding Protein Like 9 genePCAPrincipal Component AnalysisPDE8APhosphodiesterase 8A genePM2.5Particulate matter with aerodynamic diameter 2.5 micromPMDPartially methylated domainQCQuality controlRRed valuesMAFminor allele frequencyNIRnear infrared valuesNDVINormalized Difference Vegetation IndexSKA1Spindle And Kinetochore Associated Complex Subunit 1 geneSLC25A10Solute carrier family 25 member 10 geneSNPssingle nucleotide polymorphismsSESsocioeconomic statusTSStranscription start siteTssAactive transcription start siteTssAFlinkflanking active transcription start siteUCSCUniversity of California, Santa CruzUTRuntranslated regionZNF22Zinc Finger Protein 22ZNF22-AS1Zinc Finger Protein 22 Antisense RNA 1ZNF562Zinc Finger Protein 562 gene.