- Open Access
Birthweight DNA methylation signatures in infant saliva
Clinical Epigenetics volume 13, Article number: 57 (2021)
Low birthweight has been repeatedly associated with long-term adverse health outcomes and many non-communicable diseases. Our aim was to look-up cord blood birthweight-associated CpG sites identified by the PACE Consortium in infant saliva, and to explore saliva-specific DNA methylation signatures of birthweight.
DNA methylation was assessed using Infinium HumanMethylation450K array in 135 saliva samples collected from children of the NINFEA birth cohort at an average age of 10.8 (range 7–17) months. The association analyses between birthweight and DNA methylation variations were carried out using robust linear regression models both in the exploratory EWAS analyses and in the look-up of the PACE findings in infant saliva.
None of the cord blood birthweight-associated CpGs identified by the PACE Consortium was associated with birthweight when analysed in infant saliva. In saliva EWAS analyses, considering a false discovery rate p-values < 0.05, birthweight as continuous variable was associated with DNA methylation in 44 CpG sites; being born small for gestational age (SGA, lower 10th percentile of birthweight for gestational age according to WHO reference charts) was associated with DNA methylation in 44 CpGs, with only one overlapping CpG between the two analyses. Despite no overlap with PACE results at the CpG level, two of the top saliva birthweight CpGs mapped at genes associated with birthweight with the same direction of the effect also in the PACE Consortium (MACROD1 and RPTOR).
Our study provides an indication of the birthweight and SGA epigenetic salivary signatures in children around 10 months of age. DNA methylation signatures in cord blood may not be comparable with saliva DNA methylation signatures at about 10 months of age, suggesting that the birthweight epigenetic marks are likely time and tissue specific.
The existence of a relationship between intrauterine or early life exposures and health during the lifecourse has come to attention in the 1990s  and is nowadays recognized as the developmental origins of health and diseases (DOHAD) hypothesis . The intrauterine life is a critical period for adverse exposures to exert their effect , as fetal organs start developing and are sensitive to environmental stimuli that may cause an indelible imprint on future development and function.
In a hostile uterine environment caused by insults, for example poor nutrition, the fetus, to slow down its growth rate in order to match the nutrient supply, responds by developing adaptations including down-regulation of metabolic or organs function. The adaptive process, however, may cause irreversible changes in the development of some tissues and organs and predispose an individual to a higher risk of diseases not only early in life, but during the entire lifecourse .
Low birthweight may be associated with accumulation of adipose tissue and rapid weight gain during childhood , and the risk of respiratory [4, 5], metabolic [6, 7] and cardiovascular diseases [8, 9], hypertension , and neurobehavioral disorders . Low birthweight has also been associated with an increased overall mortality , while cancer incidence rises with increasing birthweight for most types of cancer [13,14,15,16,17].
The Pregnancy and Childhood Epigenetics (PACE) Consortium conducted so far the largest cord blood epigenome-wide DNA methylation study of birthweight using 8825 neonatal blood samples from 24 birth cohorts  and found that 914 CpGs, located in or near 729 genes, were associated with birthweight treated as a continuous variable. In the same study, methylation variation in 51 CpG sites was associated with high birthweight, as compared to normal birthweight, and 4 CpGs appeared to be associated with low versus normal birthweight. In additional analyses conducted on blood from 7278 children at later ages, only 1.3% of the 914 birthweight-associated differentially methylated CpGs in cord blood remained associated in childhood (2–13 years; n = 2756 children from ten studies), one in adolescence (16–18 years; n = 2906 from six studies), and none in adulthood (30–45 years; n = 1616 from three studies).
Although it seems that cord blood DNA methylation markers of birthweight do not persist at later ages, the associations observed at birth are extensive and it is important to confirm their persistence or variations over time. Also, DNA methylation profiles are tissue specific, and it would be optimal to analyse the DNA profile linked to birthweight in different tissues. Most tissues are not accessible with non-invasive methods, and blood is typically used as a surrogate, with the assumption that, being a universal body fluid, it may capture epigenetic changes of target tissues . However, blood samples collection may be difficult in large population studies, as parents may be less prone to expose their children, and especially infants, to vein puncture for blood donation. Saliva and nasal brushes are easily accessible tissues, but to date they have been much less studied [20, 21]. Replication in saliva of blood epigenetic signatures of exposures and/or outcomes and identification of saliva-specific methylation markers could make saliva an optimal candidate for studies in infants and children. Thanks to its non-invasive collection method, it would allow obtaining repeated samples in short time periods, which is necessary requisite for monitoring changes in DNA methylation over time.
Our aim was to perform a look-up, in infant saliva samples, of birthweight-related DNA methylation variation identified in cord blood, and to conduct, to our knowledge, the first saliva epigenome-wide association study (EWAS) of birthweight, to identify methylation markers that may be specific for saliva collected in infancy.
The data were derived from an epigenome-wide case–control study on wheezing nested within the NINFEA birth cohort . The NINFEA study is an Italian web-based multi-purpose mother–child cohort, aimed at exploring the relationship between early-life exposures and long-term health outcomes . Members of the cohort are children born from approximately 7500 pregnant women who between 2005 and 2016 volunteered to participate in the study, had Internet access, and had enough knowledge of Italian to complete web-based questionnaires. The children are followed up with six questionnaires completed by their mothers 6 and 18 months after delivery, and when they turn 4, 7, 10, and 13 years of age. When children were aged approximately 6 months, mothers were asked to donate their and their children’s saliva samples using a mailed sponge Oragene™ DNA self‐collection kits (OG‐250; DNA Genotek, Inc, Ottawa, Ontario, Canada). Approximately half of the participating mothers donated their and their child saliva samples, which are stored in a biobank at −80 ºC.
The case–control study was designed as EWAS of early childhood wheezing, consisting of 72 cases with at least one episode of wheezing between 6 and 18 months of age, and 72 infants without wheezing matched to cases by sex, age, and seasonality/calendar year of saliva sampling. Cases and controls were selected from a subpopulation of singletons, residents in the City of Turin, Italy, and born to mothers who did not report having asthma active during pregnancy. The baseline NINFEA questionnaire, completed at any time during pregnancy, was used to derive information on maternal and pregnancy factors, while child-related variables were collected at the first follow-up questionnaire completed approximately 6 months after delivery. Although information on children ethnic background was not available in the NINFEA cohort, almost the entire study population has both parents born in Italy, and only few study children have one of the parents born in other European countries. Therefore, the ethnic background of the children included in the study is, if not entirely, largely European. In this study, we used the following variables: maternal age at delivery (years), maternal education (low = no education, primary or secondary school vs. high = university degree or higher), parity (nulliparous vs. at least one previous pregnancy > 22 gestational weeks), maternal pre-pregnancy body mass index (BMI; kg/m2), sex, gestational age at birth (weeks), birthweight (grams), small for gestational age (yes vs. no), age at saliva sampling (continuous in months) (see Table 1 for detail).
The Illumina Infinium ® HumanMethylation450 BeadChip (Illumina,Inc, San Diego, CA, USA) was employed to evaluate DNA methylation status of over 485,000 probes in saliva samples. Details on pre‐processing of samples and data quality control can be found in the Additional File 1 (Methods, DNA methylation measurement, data pre-processing, and quality control). Quality controls and probes filtering led to the exclusion of three samples and 63,218 probes, leading to a total of 141 saliva samples and 421,782 probes for analyses.
The NINFEA study was approved by Ethical Committee of the San Giovanni Battista and CTO/CRF/Maria Adelaide Hospital of Turin, and all participating mothers gave informed consent at enrolment and at saliva donation.
For all the analyses, we pooled together cases and controls, leading to a total of 141 subjects. The analyses were, however, based on 135 subjects, as 6 (4.3%) infants had missing values in at least one of the variables included in multivariable models. Methylation levels were analysed as β-values (ratio of methylated probe intensity to overall intensity, representing 0 to 100% methylation at each probe). Although, for normality assumption, log 2-transformed β-values (M-values) may perform better when DNA methylation is used as an outcome, the interpretation of coefficients may be less intuitive. We, therefore, used β-values for more intuitive biological interpretation and easier comparability with other studies, in particular with the PACE birthweight EWAS. All the analyses were performed using R statistical computing software (version 3.6.0) and RStudio (version 1.2.1335) .
Look-up of PACE findings in infant saliva
Out of the 914 birthweight-associated CpGs identified by the PACE Consortium, 891 (97.5%) were available in the NINFEA study after quality checks and probes filtering. We used the same confounding variables of the PACE analysis, but, differently than in the PACE study, we used birthweight as the exposure and DNA methylation variation as the outcome to meet the biological temporality from birth to infancy. Both birthweight and DNA methylation were modeled as continuous variables. We used robust linear regression models adjusted for maternal age at delivery, maternal education, parity, maternal pre-pregnancy BMI, child sex, gestational age at birth, age at saliva sampling, batch, and case–control status of the original nested case–control study (wheezing between 6 and 18 months of age). Using vcovHC function in the sandwich R package , we calculated heteroscedasticity-consistent standard errors.
Although maternal smoking during pregnancy affects both birthweight  and offspring DNA methylation , we did not adjust for smoking, as it was rather infrequent in our study sample (2% prevalence). In order to account for residual technical variability and for cell-type heterogeneity, we performed surrogate variables analysis (sva)  and estimated 7 surrogate variables that were also included as covariates in the model. Cell composition was additionally estimated using the reference-based projections for saliva proposed by Zheng , and given a high correlation between the epithelial tissue component estimate and the first sva component (rho = 0.99), we only used the 7 estimated surrogate variables in the analyses. p-Values adjusted for multiple comparisons were calculated using the Bonferroni correction and the Benjamini and Hochberg false discovery rate (FDR), while histograms and quantile–quantile (QQ) plots were used to graphically compare the observed distribution of p-values versus the expected uniform distribution under the null hypothesis. Given that the direction of the association was determined by the PACE study, we also calculated one sided p-values for each CpG. It has been shown that CpGs associated with a certain trait tend to be highly correlated and that this correlation affects standard procedures for multiple testing, such as Bonferroni and Benjamini and Hochberg corrections. In the NINFEA saliva samples, the correlation between the 891 CpGs identified in the PACE EWAS was 0.48, which is indeed much higher than the reported saliva genome-wide mean correlation of 0.12 . We, therefore, reported the permutation p-values that take into account the distribution of p-values under the null hypothesis and are suggested as a gold standard for the settings where the underlying correlation between CpGs is high.
Epigenome-wide association analyses
Two epigenome-wide association analyses were conducted, first with birthweight as a continuous exposure variable and, second, with small for gestational age (SGA) as a binary exposure variable. The latter was defined as the lowest 10th percentile of the World Health Organization birthweight for sex and gestational age charts , with the remaining population above the 10th percentile as the reference group. As in our study population only 7 children (5%) were classified as large for gestational age (based on the 90th percentile), we had no power to analyse this trait separately and opted for a more conservative approach by including them in the reference group. In both EWAS analyses, we used robust linear regression models adjusted for sex, age at saliva sampling, gestational age, maternal age, parity, maternal pre-pregnancy BMI, maternal education, batch, estimated surrogate variables, and case–control status of the original nested case–control study. p-Values adjusted for multiple comparisons were calculated using the Benjamini and Hochberg false discovery rate (FDR), and Volcano plots were used to visually present the results. In the EWAS of SGA, as a sensitivity analysis we provisionally excluded children born pre-term (< 37 gestational weeks at birth).
To assess whether the age at saliva sampling could have influenced the findings on the top CpG sites identified in the EWAS analyses, we tested the associations of the age at saliva sampling as a continuous variable (in months) with the methylation levels in the top CpG sites using the robust linear regression models adjusted for sex, batch and cell type composition estimated with the reference-based projections for saliva proposed by Zheng .
Finally, we performed a look-up of the saliva top findings in the PACE summary results, publicly available via Zenodo: https://zenodo.org/record/2222287#.YDOhCtWSmUk.
CpGs annotation and functional analysis
Gene Ontology (GO) and Kyoto Encyclopedia of gene and Genomes (KEGG) enrichment analyses were carried out to identify possible functional pathways in the saliva birthweight-associated CpGs set.
In order to compare previously reported associations of epigenome-wide birthweight-associated CpGs and our own results, we searched for findings reported in the EWAS Catalog (http://www.ewascatalog.org; accessed on 12 February 2021) and EWAS Atlas (https://bigd.big.ac.cn/ewas/tools; accessed on 12 February 2021), looking for both CpG and gene level match.
After this first step based on the overlap between our results and those of other EWAS on the same trait, we further accessed EWAS Atlas to examine whether CpGs identified in our study were previously associated with traits other than birthweight.
We also looked in the GWAS Catalog for traits associated with specific single-nucleotide polymorphisms (SNPS) on genes on which these CpGs mapped.
Table 1 shows the characteristics of the study population. The mean maternal age was 35 years; 72% of the mothers had a high educational level. The mean age at saliva sampling was 10.8 months (median 10.4, range 7–17). The mean birthweight was 3242 g, 5% of the children were born pre-term and 15.6% of children were born small for gestational age, respectively).
Look-up of PACE findings in infant saliva
Of the 891 CpG sites associated with birthweight in the PACE study, 52 (5.8%) had a one sided p-value < 0.05 in our study (Fig. 1, Table 2). However, none of the CpGs survived the Bonferroni correction (p < 5.61 × 10−5) or had a FDR < 0.05 (Table 2, Additional file 2: Table S1). There was a 47% concordance in the direction of the coefficients between the PACE and our study (binomial sign test p-value = 0.11, 95% confidence intervals (CIs) 0.44; 0.51). When we used the permutation test to take into account strong mean pairwise correlation between the DNA methylation values of the 891 CpG sites, the minimum permutation p-value was 0.56.
EWAS for continuous birthweight
Out of the 421,782 probes analysed, 8.9% (N = 37,365) were associated with birthweight with a nominal p-value < 0.05. After correction for multiple testing, 44 CpG sites had a FDR < 0.05 (Table 3 and Fig. 2a). Their coefficient estimates ranged from – 0.31 to 0.57, which corresponds to a methylation increase of 0.57% with a 100 gr increase in birthweight. The largest effect was observed for cg02727104 located 13 kb down SOHLH2.
None of the 44 saliva birthweight-associated CpGs overlapped with the birthweight related CpGs identified in cord blood in the PACE study under the Bonferroni threshold. When considering a nominal p-value < 0.10 in the PACE study (N = 78,489 CpGs), there was an overlap in 6 CpG sites, but only one with the same direction of the effect (cg22896429) (Additional file 3: Table S1). At gene level, there was an overlap in 12 genes identified in PACE under the FDR-corrected p-value threshold (with at least one CpG with the same direction of the effect). The average absolute pairwise Pearson's correlation coefficient between the β-values of the identified 44 CpGs was 0.23, which is higher than the mean pairwise genome-wide correlation coefficient in the NINFEA saliva samples but lower than the mean pairwise correlation coefficient between 891 CpGs identified by the PACE Consortium .
EWAS for small for gestational age
We found 44 CpGs associated with SGA at a FDR of less than 0.05 (Table 4 and Fig. 2b). The largest coefficient showed a 4.1% difference in methylation when comparing small with non-small for gestational age (cg12322146, located 125 kb down RBMS3); the largest negative association was -2.3% (cg03066788, located at PNOC). Only one CpG (cg18072629, located at GATA3) overlapped with the 44 CpGs found to be associated with continuous birthweight in our sample. None of the 44 saliva SGA-associated CpGs overlapped with the 914 CpGs associated with continuous birthweight (at Bonferroni threshold), or with the 4 CpGs associated with low birthweight (< 2500 g) in the PACE study conducted on cord blood samples. When considering a nominal p-value < 0.10 in the PACE study, there was an overlap in 7 CpG sites, but only two with consistent direction of the effect (cg06234201, cg15847996). (Additional file 3: Table S1). At gene level, there was an overlap in 9 genes identified in PACE under the FDR-corrected p-value threshold (with at least one CpG with consistent direction of the effect).
Findings were practically unchanged when we excluded preterm infants from the EWAS analysis (data not shown). We found no association between age at saliva sampling and methylation status for any of the top CpGs associated with continuous birthweight or with SGA in our two EWAS analyses.
Out of 44 CpG sites associated with continuous birthweight in our EWAS, 28 (64%) showed a nominal p-value < 0.05 in the EWAS for SGA; all with a direction of the effect that was consistent with the results for birthweight. Likewise, 26 (59%) out of 44 SGA-related CpGs reveal nominal p-value < 0.05 in the EWAS of continuous birthweight, with a consistent direction of the effect.
CpGs Annotation and functional analysis
None of the two sets of CpGs identified in our study, the 44 birthweight-associated CpGs and the 44 SGA-associated CpGs, showed functional enrichment of GO or KEGG terms.
Also, there was no overlap between the CpGs identified in our EWAS on continuous birthweight and 995 birthweight-related CpGs reported in 4 cord-blood and subcutaneous adipose tissue studies from the EWAS Atlas data [32,33,34], including the PACE study . In EWAS Atlas and EWAS Catalog, DNA methylation in six out of 87 genes was previously associated with birthweight, all in the PACE study at Bonferroni threshold. In addition, some genetic variants in these six genes have been previously associated with obesity , adult BMI [36,37,38,39,40,41,42], BMI-adjusted waist-hip ratio [37, 38, 41, 43, 44], high-density lipoprotein cholesterol measurement [45,46,47], and visceral adipose tissue measurement .
In the EWAS Catalog, there were 34 CpGs associated with birthweight in four studies using DNA methylation in cord blood [49,50,51,52], but none of these overlapped with the 87 CpGs found to be associated with birthweight or SGA in our study. In the EWAS Atlas, DNA methylation in 29 out of these 87 CpGs was associated with different traits (Additional file 1: Tables S1 and S2). DNA methylation variation at eleven of them has been previously associated with insulin resistance  (cg03045325), colorectal cancer  (cg02727104, cg26332310, cg12322146), obesity  (cg03066788), bariatric surgery  (cg20515787, cg02547025), mortality  (cg06234201), gestational diabetes mellitus  (cg00383136), amount of visceral adipose tissue  (cg20388707), and gestational age  (cg00701706).
Some of the birthweight- and SGA-associated CpGs in our study map in genes which variants were associated with the following traits: birthweight (USH2A), BMI (CENPO, E2F3, RPTOR, SNTB2, PNOC, LGR4), body fat distribution (CENPO), waist-hip ratio (SYTL2, ZNF423, FOXA3, LMNB2, COL5A1,LGR4), BMI-adjusted waist hip ratio (ZNF423, AFF3), cardiovascular disease (SIPA1L2, FOXA3, RAB37, INPP5A), subcutaneous or visceral adipose tissue measurement (FOXA3, RPTOR), gestational age (CFAP46, INPP5A), lipoprotein cholesterol measurement (DMTN, SNTB2), total cholesterol measurement (E2F3, DMTN), type I diabetes nephropathy (AFF3), type II diabetes mellitus (RAMP1, SYCE1L, ZNF710). Finally at gene level, out of saliva 87 CpGs associated with birthweight and SGA, DNA methylation at 14 genes was associated with birthweight, 22 with BMI and 32 with obesity in previous studies found in EWAS Atlas and EWAS Catalog. Full results are shown in Additional file 1: Tables S1 and S2.
In this study, we investigated the association between birthweight and methylation patterns in saliva samples taken at around 10 months of age. The cord blood methylation signatures of birthweight, found by a large study of the PACE Consortium, were not confirmed in infant saliva.
We identified 87 infant saliva-specific signatures of birthweight or SGA, of which two overlap with PACE results at gene level with the same direction of effect (MACROD1 and RPTOR). Interestingly, single-nucleotide polymorphisms in these genes have been previously associated with obesity, adult BMI, BMI-adjusted waist-hip ratio, high-density lipoprotein cholesterol and visceral adipose tissue levels. Moreover, expanding the look-up in the full PACE results, 13 CpGs find correspondence with a nominal p-value < 0.10 but only 3 have the same direction of the effect, namely cg22896429, cg06234201, cg15847996.
DNA methylation variation at some of the 87 loci identified in our study has been previously associated with multiple traits, such as insulin resistance, colorectal cancer, obesity, and gestational diabetes mellitus. Moreover, the SGA-associated locus cg26615232 maps within USH2A, which variant has been associated with birthweight in a study  that performed GWAS meta-analyses of fetal genetic variants in 321,223 individuals of European ancestry.
It has been repeatedly shown that DNA methylation at many sites is not temporally stable and that each tissue has its unique epigenetic landscape that likely reflects its specific function and response to environmental exposures [21, 64]. For example, a cross-sectional study on 1019 infants  compared the associations between preterm birth and genome-wide DNA methylation profiles using both cord tissue and cord blood samples. The results highlighted differences between the two tissues in DNA methylation variation associated with preterm birth, with only a minority of overlapping CpGs. In DNA from cord tissue, DNA methylation analysis showed enrichment of differentially methylated regions in genes involved in molecular pathways related to fetal growth and development (i.e. Wnt signaling, bone remodeling, and extracellular matrix organization), while in cord blood immune response pathways (i.e. regulation of T cell differentiation, inositol lipid-mediated signaling, and regulation of RNA stability) were enriched. Therefore, it is reasonable to speculate that saliva and blood, which have different functions, include different cell types and have different embryonic origin, and different mechanisms of exposure and responses to environmental factors do not share the same DNA methylation response to fetal growth, birthweight, and their risk factors. Although our results suggest different cord blood and saliva methylation patterns related to birthweight, we cannot exclude that the relatively small sample size of our study contributed to the lack of overlap with the PACE results, due to a reduced power to detect small effects.
In addition to the tissue specificity of DNA methylation, there are also age-related changes (birth vs. infancy) that could explain differences between our and the PACE Consortium findings.
In most tissues, DNA methylation may also vary substantially with time, especially during periods of life associated with high plasticity and fast development. Consistently, the PACE study found that differential methylation associated with birthweight in neonates persisted only minimally across childhood and disappeared by adulthood .
This result is consistent with another longitudinal study in which birthweight- and gestational age-related DNA methylation changes were investigated in cord blood and peripheral blood at ages 7 and 17 in more than 900 children . Across the majority of CpG sites that showed differential methylation in cord blood, a pattern of fast evolution was observed during early childhood that stops with adolescence, providing evidence for the lack of persistence of early life methylation differences.
This is evident also in two studies that analysed infant saliva samples. A longitudinal study with repeated saliva samples collected at birth and at 1 year of age from 50 preterm and 40 infants born at term showed that DNA methylation at the differentially methylated region (DMR) of IGF2 and FKBP5 at birth was lower in preterm infants compared with term infants, but these differences did not persist at 1 year of age . Also, another study on 214 infant saliva samples (62 collected at 6 weeks of age, 30 collected at 52 weeks of age and 61 collected at both ages) showed an age-dependent variation of DNA methylation, with a clear difference between saliva DNA methylation at 6 and 52 weeks of age .
In our study, we could not distinguish between the tissue- and the time-related differences in DNA methylation of birthweight-associated CpG sites, as we had no repeated saliva samples or blood samples collected simultaneously with saliva. It would be interesting in future studies with repeat samples and also with different tissues, to investigate whether the differences between our and PACE findings are due to the time- and tissue- dependent DNA methylation changes, as reported by previous studies both for birthweight and gestational age. Even though the age at saliva sampling varied between 7 and 17 months in our study, age was not associated with DNA methylation variation in any of the CpG sites associated with continuous birthweight and SGA, suggesting that the observed associations are unlikely influenced by the different age of saliva collection. It should be, however, noted that the age range of our population was rather narrow (7–17 months), so we cannot exclude that age could influence DNA methylation in these CpGs when considered across childhood/adolescence.
We supplemented the analysis of birthweight with the analysis on SGA and found different CpG sites associated with the two exposures, with only 1 overlapping CpG when considering FDR-corrected p-values. Although this is in line with the PACE analyses, where the 4 CpGs associated with low birthweight (< 2500 g) did not overlap with the 914 associated with continuous birthweight, approximately 60% of CpGs identified in our two EWAS overlapped at the nominal p-value of 0.05, with seemingly the same direction.
The main strength of our study was the possibility to analyse DNA methylation in saliva samples, which are easy to collect at any age in childhood and are, therefore, good candidates for future large DNA methylation studies. Moreover, while it is probably unfeasible to collect repeated blood samples over relatively short time periods, and especially in the context of large population studies based on children, it is possible to use saliva to monitor changes in DNA methylation over time. As salivary DNA methylation is poorly studied, especially in newborns and infants, it would be important to understand if, and for which specific traits and exposures, salivary and blood DNA methylation can be used interchangeably, and when they clearly show distinct methylation signatures. Previous DNA methylation studies using saliva samples focused on neurobehavioural conditions [69,70,71], respiratory traits , and cancer research , but to our knowledge the associations between birthweight and saliva DNA methylation has not been studied so far.
Although the small sample size of our study, especially in comparison with the PACE study, may have had an impact on EWAS analyses, after an FDR-correction, we identified some novel CpGs associated with birthweight, which are likely to be saliva-specific birthweight signatures. These findings need replication in independent saliva EWAS, and we cannot exclude that additional saliva DNA methylation variations related to offspring birthweight may be identified by larger studies. We argue that the sample size had less impact on our look-up analysis of PACE findings in saliva, as the number of tests performed in the look-up analysis was much lower than the number of tests of the two EWAS.
Our study population was not selected at random from the entire NINFEA cohort, but on the presence/absence of infant wheezing between 6 and 18 months of age, which was a selection factor for the original case–control study. As SGA is a well-known risk factors for wheezing [73, 74], this selection led to a relatively high prevalence of SGA infants (15.6%; 21% in wheezing cases and 10% in controls). To account for this selection, all the analyses were adjusted also for infant wheezing. Moreover, SGA definition is based on an arbitrary population-based reference cut-off, while, biologically, size for gestational age can be considered as a continuous trait. Therefore, the high prevalence of SGA in our study population allowed us to capture the lowest quintile of the size for gestational age distribution, providing more power for EWAS.
Finally, it should be noted that the prevalence of maternal smoking during pregnancy was rather low in our study population (2% compared with 7.7% of maternal smoking during pregnancy and with 14% of maternal smoking just before pregnancy in the entire cohort) . In addition to random variation, the selection process may have played a role also in this case, although with an unexpected direction. Therefore, our findings are less generalizable to the population of women who smoke during pregnancy, but remain applicable to all non-smoking women and those who stopped smoking before pregnancy, which represent the majority of pregnant women.
In conclusion, our study provides an indication of the birthweight and small for gestational age epigenetic salivary signatures in children around 10 months of age and suggests that DNA methylation signatures of birthweight likely differ between cord blood and infant saliva. Further insights are needed to understand whether these differences are due to biological differences between the two tissues or could be attributed to age-related DNA methylation changes.
Availability of data and materials
Data from the NINFEA cohort underlying the findings reported in this study are available to researchers who meet the criteria for access to confidential data and upon reasonable request. Data availability contact: Prof. Lorenzo Richiardi (email@example.com).
Adequate for gestational age
Small for gestational age
False discovery rate
Epigenome-wide association study
Hales CN, Barker DJP. The thrifty phenotype hypothesis. Br Med Bull. 2001;60:5–20. https://doi.org/10.1093/bmb/60.1.5.
Gluckman PD, Hanson MA, Cooper C, Thornburg KL. Effect of in utero and early-life conditions on adult health and disease. New Engl J Med. 2008;359(1):61. https://doi.org/10.1056/NEJMra0708473.
Simeoni U, Armengaud JB, Siddeek B, Tolsa JF. Perinatal origins of adult disease. Neonatology. 2018;113(4):393–9. https://doi.org/10.1159/000487618.
Popovic M, et al. Infant weight trajectories and early childhood wheezing: the NINFEA birth cohort study. Throax. 2016;71:1091–6. https://doi.org/10.1136/thoraxjnl-2015-208208.
Kwinta P, Pietrzyk JJ. Preterm birth and respiratory disease in later life. Expert Rev Respir Med. 2010;4(5):593–604. https://doi.org/10.1586/ers.10.59.
Dabelea D, et al. Association of intrauterine exposure to maternal diabetes and obesity with type 2 diabetes in youth: the SEARCH case–control study. Diabetes Care. 2008;31(7):1422–6. https://doi.org/10.2337/dc07-2417.
Chen PY, et al. Prenatal growth patterns and birthweight are associated with differential DNA methylation and gene expression of cardiometabolic risk genes in human placentas: a discovery-based approach. Reprod Sci. 2018;25(4):523–39. https://doi.org/10.1177/1933719117716779.
Painter RC, Roseboom TJ, Bleker OP. Prenatal exposure to the Dutch famine and disease in later life: an overview. Reprod Toxicol. 2005;20(3):345–52. https://doi.org/10.1016/j.reprotox.2005.04.005.
Szathmári M, Vásárhelyi B, Reusz G, Tulassay T. Adult cardiovascular risk factors in premature babies. Lancet. 2000;356(9233):939–40. https://doi.org/10.1016/S0140-6736(05)73920-9.
Zhang H, et al. In utero and postnatal exposure to environmental tobacco smoke, blood pressure, and hypertension in children: the Seven Northeastern Cities study. Int J Environ Health Res. 2019. https://doi.org/10.1080/09603123.2019.1612043.
Aarnoudse-Moens CSH, Weisglas-Kuperus N, Van Goudoever JB, Oosterlaan J. Meta-analysis of neurobehavioral outcomes in very preterm and/or very low birth weight children. Pediatrics. 2009;124(2):717–28. https://doi.org/10.1542/peds.2008-2816.
Risnes KR, et al. Birthweight and mortality in adulthood: a systematic review and meta-analysis. Int J Epidemiol. 2011;40(3):647–61. https://doi.org/10.1093/ije/dyq267.
Ahlgren M, Wohlfahrt J, Olsen LW, Sørensen TIA, Melbye M. Birth weight and risk of cancer. Cancer. 2007;110(2):412–9. https://doi.org/10.1002/cncr.22773.
Vatten LJ, et al. Birth weight as a predictor of breast cancer: a case–control study in Norway. Br J Cancer. 2002;86(1):89–91. https://doi.org/10.1038/sj.bjc.6600011.
McCormack VA, Silva IDS, Koupil I, Leon DA, Lithell HO. Birth characteristics and adult cancer incidence: Swedish cohort of over 11,000 men and women. Int J Cancer. 2005;115(4):611–7. https://doi.org/10.1002/ijc.20915.
Xue F, Michels KB. Intrauterine factors and risk of breast cancer: a systematic review and meta-analysis of current evidence. Lancet Oncol. 2007;8(12):1088–100. https://doi.org/10.1016/S1470-2045(07)70377-7.
Paltiel O, et al. Birthweight and childhood cancer: preliminary findings from the international childhood cancer cohort consortium (I4C). Paediatr Perinat Epidemiol. 2015;29(4):335–45. https://doi.org/10.1111/ppe.12193.
Küpers LK, et al. Meta-analysis of epigenome-wide association studies in neonates reveals widespread differential DNA methylation associated with birthweight. Nat Commun. 2019;10(1):1–11. https://doi.org/10.1038/s41467-019-09671-3.
Jin Z, Liu Y. DNA methylation in human diseases Introduction to DNA methylation. 2018. https://doi.org/10.1016/j.gendis.2018.01.002.
Lowe R, et al. Buccals are likely to be a more informative surrogate tissue than blood for epigenome-wide association studies. Epigenetics. 2013. https://doi.org/10.4161/epi.24362.
Lin X, et al. Choice of surrogate tissue influences neonatal EWAS findings. BMC Med. 2017;15(1):211. https://doi.org/10.1186/s12916-017-0970-x.
Popovic M, et al. Differentially methylated DNA regions in early childhood wheezing: an epigenome-wide study using saliva. Pediatr Allergy Immunol. 2019;30(3):305–14. https://doi.org/10.1111/pai.13023.
Richiardi L, Baussano I, Vizzini L, Douwes J, Pearce N, Merletti F. Feasibility of recruiting a birth cohort through the Internet: The experience of the NINFEA cohort. Eur J Epidemiol. 2007;22(12):831–7. https://doi.org/10.1007/s10654-007-9194-2.
R. F. for S. C., R Core Team, R: A Language and Environment for Statistical Computing; 2019.
Zeileis A. Econometric computing with HC and HAC covariance matrix estimators.
Brand JS, et al. Associations of maternal quitting, reducing, and continuing smoking during pregnancy with longitudinal fetal growth: findings from Mendelian randomization and parental negative control studies. PLOS Med. 2019;16(11):e1002972. https://doi.org/10.1371/journal.pmed.1002972.
Joubert BR, et al. DNA methylation in newborns and maternal smoking in pregnancy: genome-wide consortium meta-analysis. Am J Hum Genet. 2016;98(4):680–96. https://doi.org/10.1016/j.ajhg.2016.02.019.
Leek JT, Storey JD. Capturing heterogeneity in gene expression studies by surrogate variable analysis. PLoS Genet. 2007;3(9):1724–35. https://doi.org/10.1371/journal.pgen.0030161.
Zheng SC, et al. A novel cell-type deconvolution algorithm reveals substantial contamination by immune cells in saliva, buccal and cervix. Epigenomics. 2018;10(7):925–40. https://doi.org/10.2217/epi-2018-0037.
Popovic M, Fasanelli F, Fiano V, Biggeri A, Richiardi L. Increased correlation between methylation sites in epigenome-wide replication studies: impact on analysis and results. Epigenomics. 2017;9(12):1489–502. https://doi.org/10.2217/epi-2017-0073.
Kiserud T, et al. The World Health Organization fetal growth charts: a multinational longitudinal study of ultrasound biometric measurements and estimated fetal weight. PLoS Med. 2017;14(1):e1002220. https://doi.org/10.1371/journal.pmed.1002220.
Van Dijk SJ, et al. DNA methylation in blood from neonatal screening cards and the association with BMI and insulin sensitivity in early childhood. Int J Obes. 2018;42(1):28–35. https://doi.org/10.1038/ijo.2017.228.
Gillberg L, et al. Adipose tissue transcriptomics and epigenomics in low birthweight men and controls: role of high-fat overfeeding. Diabetologia. 2016;59(4):799–812. https://doi.org/10.1007/s00125-015-3852-9.
Hannon E, et al. Variable DNA methylation in neonates mediates the association between prenatal smoking and birth weight. Philos Trans R Soc B Biol Sci. 2019. https://doi.org/10.1098/rstb.2018.0120.
Berndt SI, et al. Genome-wide meta-analysis identifies 11 new loci for anthropometric traits and provides insights into genetic architecture. Nat Genet. 2013;45(5):501–12. https://doi.org/10.1038/ng.2606.
Locke AE, et al. Genetic studies of body mass index yield new insights for obesity biology. Nature. 2015;518(7538):197–206. https://doi.org/10.1038/nature14177.
Justice AE, et al. Genome-wide meta-analysis of 241,258 adults accounting for smoking behaviour identifies novel loci for obesity traits. Nat Commun. 2017. https://doi.org/10.1038/ncomms14977.
Winkler TW, et al. The influence of age and sex on genetic associations with adult body size and shape: a large-scale genome-wide interaction study. PLoS Genet. 2015;11(10):1–42. https://doi.org/10.1371/journal.pgen.1005378.
Pulit SL, et al. Meta-analysis of genome-wide association studies for body fat distribution in 694 649 individuals of European ancestry. Hum Mol Genet. 2019;28(1):166–74. https://doi.org/10.1093/hmg/ddy327.
Akiyama M, et al. Genome-wide association study identifies 112 new loci for body mass index in the Japanese population. Nat Genet. 2017;49(10):1458–67. https://doi.org/10.1038/ng.3951.
Zhu Z, et al. Shared genetic and experimental links between obesity-related traits and asthma subtypes in UK Biobank. J Allergy Clin Immunol. 2020;145(2):537–49. https://doi.org/10.1016/j.jaci.2019.09.035.
Hoffmann TJ, Choquet H, Yin J, Banda Y, Kvale MN, Glymour M. A large multiethnic genome-wide association study. Genetics. 2018;210(October):499–515. https://doi.org/10.1534/genetics.118.301479.
Shungin D, et al. New genetic loci link adipose and insulin biology to body fat distribution. Nature. 2015;518(7538):187–96. https://doi.org/10.1038/nature14132.
Lotta LA, et al. Association of genetic variants related to gluteofemoral vs abdominal fat distribution with type 2 diabetes, coronary disease, and cardiovascular risk factors. JAMA J Am Med Assoc. 2018;320(24):2553–63. https://doi.org/10.1001/jama.2018.19329.
De Vries PS, et al. Multiancestry genome-wide association study of lipid levels incorporating gene–alcohol interactions. Am J Epidemiol. 2019;188(6):1033–54. https://doi.org/10.1093/aje/kwz005.
Hoffmann TJ, et al. A large electronic-health-record-based genome-wide study of serum lipids. Nat Genet. 2018;50(3):401–13. https://doi.org/10.1038/s41588-018-0064-5.
Qi G, Chatterjee N. Heritability informed power optimization (HIPO) leads to enhanced detection of genetic associations across multiple traits. PLoS Genet. 2018;14(10):e1007549. https://doi.org/10.1371/journal.pgen.1007549.
Karlsson T, et al. Contribution of genetics to visceral adiposity and its relation to cardiovascular and metabolic disease. Nat Med. 2019;25(9):1390–5. https://doi.org/10.1038/s41591-019-0563-7.
Lin X, et al. Developmental pathways to adiposity begin before birth and are influenced by genotype, prenatal environment and epigenome. BMC Med. 2017;15(1):1–18. https://doi.org/10.1186/s12916-017-0800-1.
Engel SM et al. Original contribution neonatal genome-wide methylation patterns in relation to birth weight in the Norwegian mother and child cohort. https://doi.org/10.1093/aje/kwt433
Agha G, et al. Birth weight-for-gestational age is associated with DNA methylation at birth and in childhood. Clin Epigenetics. 2016;8(1):118. https://doi.org/10.1186/s13148-016-0285-3.
Tsai PC, et al. DNA methylation changes in the IGF1R gene in birth weight discordant adult monozygotic twins. Twin Res Hum Genet. 2015;18(6):635–46. https://doi.org/10.1017/thg.2015.76.
Arpón A, et al. Epigenome-wide association study in peripheral white blood cells involving insulin resistance. Sci Rep. 2019. https://doi.org/10.1038/s41598-019-38980-2.
Zhu L, et al. Genome-wide DNA methylation profiling of primary colorectal laterally spreading tumors identifies disease-specific epimutations on common pathways. Int J Cancer. 2018;143(10):2488–98. https://doi.org/10.1002/ijc.31765.
Kirchner H, et al. Altered DNA methylation of glycolytic and lipogenic genes in liver from obese and type 2 diabetic patients. Mol Metab. 2016;5(3):171–83. https://doi.org/10.1016/j.molmet.2015.12.004.
Fraszczyk E, et al. The effects of bariatric surgery on clinical profile, DNA methylation, and ageing in severely obese patients. Clin Epigenetics. 2020;12(1):14. https://doi.org/10.1186/s13148-019-0790-2.
Svane AM, et al. DNA methylation and all-cause mortality in middle-aged and elderly Danish twins. Genes (Basel). 2018. https://doi.org/10.3390/genes9020078.
Alexander J, et al. Offspring sex impacts DNA methylation and gene expression in placentae from women with diabetes during pregnancy. PLoS ONE. 2018. https://doi.org/10.1371/journal.pone.0190698.
Hohos NM, et al. CD4+ and CD8+ T-cell-specific DNA cytosine methylation differences associated with obesity. Obesity. 2018;26(8):1312–21. https://doi.org/10.1002/oby.22225.
Gross AM, et al. Methylome-wide analysis of chronic HIV infection reveals five-year increase in biological age and epigenetic targeting of HLA. Mol Cell. 2016;62(2):157–68. https://doi.org/10.1016/j.molcel.2016.03.019.
Henneman P, et al. Widespread domain-like perturbations of DNA methylation in whole blood of Down syndrome neonates. PLoS ONE. 2018. https://doi.org/10.1371/journal.pone.0194938.
Wozniak MB, et al. Integrative genome-wide gene expression profiling of clear cell renal cell carcinoma in Czech Republic and in the United States. PLoS ONE. 2013. https://doi.org/10.1371/journal.pone.0057886.
Warrington NM, et al. Maternal and fetal genetic effects on birth weight and their relevance to cardio-metabolic risk factors. Nat Genet. 2019;51(5):804–14. https://doi.org/10.1038/s41588-019-0403-1.
Armstrong DA, Lesseur C, Conradt E, Lester BM, Marsit CJ. Global and gene-specific DNA methylation across multiple tissues in early infancy: Implications for children’s health research. FASEB J. 2014;28(5):2088–97. https://doi.org/10.1096/fj.13-238402.
Wu Y, et al. Analysis of two birth tissues provides new insights into the epigenetic landscape of neonates born preterm. Clin Epigenetics. 2019;11(1):1–12. https://doi.org/10.1186/s13148-018-0599-4.
Simpkin AJ, et al. Longitudinal analysis of DNA methylation associated with birth weight and gestational age. Hum Mol Genet. 2015;24(13):3752–63. https://doi.org/10.1093/hmg/ddv119.
C. Piyasena et al., Dynamic changes in DNA methylation occur during the first year of life in preterm infants. Front Endocrinol Lausanne 2016. doi: https://doi.org/10.3389/fendo.2016.00158.
Wikenius E, Moe V, Smith L, Heiervang ER, Berglund A. DNA methylation changes in infants between 6 and 52 weeks. Sci Rep. 2019;9(1):1–12. https://doi.org/10.1038/s41598-019-54355-z.
Lester BM, Conradt E, LaGasse LL, Tronick EZ, Padbury JF, Marsit CJ. Epigenetic programming by maternal behavior in the human infant. Pediatrics. 2018. https://doi.org/10.1542/peds.2017-1890.
Conradt E, et al. DNA methylation of NR3c1 in infancy: Associations between maternal caregiving and infant sex. Infant Ment Health J. 2019;40(4):513–22. https://doi.org/10.1002/imhj.21789.
Sherwood WB, et al. Duration of breastfeeding is associated with leptin (LEP) DNA methylation profiles and BMI in 10-year-old children. Clin Epigenetics. 2019;11(1):128. https://doi.org/10.1186/s13148-019-0727-9.
Lim Y, et al. Salivary DNA methylation panel to diagnose HPV-positive and HPV-negative head and neck cancers. BMC Cancer. 2016;16(1):749. https://doi.org/10.1186/s12885-016-2785-0.
Mebrahtu TF, Feltbower RG, Parslow RC. Effects of birth weight and growth on childhood wheezing disorders: findings from the Born in Bradford Cohort. BMJ Open. 2015;5(11):e009553. https://doi.org/10.1136/bmjopen-2015-009553.
Xu XF, Li YJ, Sheng YJ, Liu JL, Tang LF, Chen ZM. Effect of low birth weight on childhood asthma: a meta-analysis. BMC Pediatr. 2014. https://doi.org/10.1186/1471-2431-14-275.
The authors are grateful to all the participants of the NINFEA cohort.
The NINFEA study was partially funded by the Compagnia San Paolo Foundation. This research was partially funded by the Italian Ministry for Education, University and Research (Ministero dell’Istruzione, dell’Università e della Ricerca – MIUR) under the programme “Dipartimenti di Eccellenza 2018–2022", and by the European Union’s Horizon2020 research and innovation programme under grant agreement no. 733206, LIFE-CYCLE project.
Ethics approval and consent to participate
NINFEA study was approved by Ethical Committee of the San Giovanni Battista and CTO/CRF/Maria Adelaide Hospital of Turin, and all participating mothers gave informed consent at enrolment and at saliva donation.
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
. DNA methylation measurement, data pre-processing, and quality control, traits associated in EWAS Atlas with 44-saliva birthweight related CpGs, traits associated in EWAS Atlas with saliva44-SGA-related CpGs.
. Results from the look-up of PACE findings.
. Look-up of the 87 saliva CpG sites in the PACE birthweight cord blood EWAS.
About this article
Cite this article
Moccia, C., Popovic, M., Isaevska, E. et al. Birthweight DNA methylation signatures in infant saliva. Clin Epigenet 13, 57 (2021). https://doi.org/10.1186/s13148-021-01053-1
- Saliva DNA methylation
- Birth cohort