Abstract
Preterm birth is associated with dysconnectivity of structural brain networks and is a leading cause of neurocognitive impairment in childhood. Variation in DNA methylation (DNAm) is associated with early exposure to extrauterine life but there has been little research exploring its relationship with brain development.
Using genome-wide DNA methylation data from saliva of 258 neonates, we investigated the impact of gestational age on the methylome and performed functional analysis to identify enriched gene sets from probes that contributed to differentially methylated probes (DMPs) or regions (DMRs). We tested the hypothesis that variation in DNAm could underpin the association between preterm birth and atypical brain development by linking DMPs with measures of white matter connectivity derived from diffusion MRI metrics: peak width of skeletonised mean diffusivity (PSMD), fractional anisotropy (PSFA) and neurite density index (PSNDI).
Gestational age at birth was associated with widespread differential methylation, with genome-wide significant associations observed for 8,870 CpG probes (p<3.6×10−8) and 1,767 differentially methylated regions. Functional analysis identified 14 enriched gene ontology terms pertaining to cell-cell contacts and cell-extracellular matrix contacts. Principal component analysis of probes with genome-wide significance revealed a first principal component (PC1) that explained 23.5% of variance in DNAm, and this was negatively associated with gestational age at birth. PC1 was associated with PSMD (β=0.349, p=8.37×10−10) and PSNDI (β=0.364, p=4.15×10−5), but not with PSFA (β=−0.035, p=0.510); these relationships mirrored the imaging metrics’ associations with gestational age at birth.
Gestational age at birth has a profound and widely distributed effect on the neonatal saliva methylome. Enriched gene ontology terms related to cell-cell contacts reveal pathways that could mediate the effect of early life environmental exposures on development. Finally, associations between differential DNAm and image markers of white matter tract microstructure suggest that variation in DNAm may provide a link between preterm birth and the dysconnectivity of developing brain networks that characterises atypical brain development in preterm infants.
Competing Interest Statement
Riccardo Marioni has received a speaker fee from Illumina and is an advisor to the Epigenetic Clock Development Foundation. Lee Murphy has carried out paid presentations and consultancy for Illumina. The remaining authors report no competing interests.
Funding Statement
Emily Wheater is supported by the Wellcome Trust Translational Neuroscience PhD fellowship programme at the University of Edinburgh (203769/Z/16/A). This work was supported by Theirworld (www.theirworld.org) and was carried out in the MRC Centre for Reproductive Health, which was funded by MRC Centre Grant (MRC G1002033).
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:
National Research Ethics Service, South East Scotland Research Ethics Committee (11/55/0061, 13/SS/0143 and 16/SS/0154).
All necessary patient/participant consent has been obtained and the appropriate institutional forms have been archived.
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 and uploaded the relevant EQUATOR Network research reporting checklist(s) and other pertinent material as supplementary files, if applicable.
Yes
Data Availability
The atlas with templates can be found at https://git.ecdf.ed.ac.uk/jbrl/ena and the code necessary to calculate histogram-based metrics is at https://git.ecdf.ed.ac.uk/jbrl/psmd. Requests for original image data will be considered through the BRAINS governance process: www.brainsimagebank.ac.uk. DNA methylation data will be deposited in NCBI’s
Gene Expression Omnibus upon publication.
https://git.ecdf.ed.ac.uk/jbrl/ena
Abbreviations
- ADHD
- attention deficit hyperactivity disorder
- FA
- fractional anisotropy
- IQ
- intelligence quotient
- MD
- mean diffusivity
- NODDI
- neurite orientation dispersion and density imaging
- NDI
- neurite density index
- PS
- peak width skeletonised
- PCA
- principal component analysis
- SVA
- surrogate variable analysis.
