The influence of genetic variation on complex diseases is potentially mediated through a range of highly dynamic epigenetic processes exhibiting temporal variation during development and later life. ...Here we present a catalogue of the genetic influences on DNA methylation (methylation quantitative trait loci (mQTL)) at five different life stages in human blood: children at birth, childhood, adolescence and their mothers during pregnancy and middle age.
We show that genetic effects on methylation are highly stable across the life course and that developmental change in the genetic contribution to variation in methylation occurs primarily through increases in environmental or stochastic effects. Though we map a large proportion of the cis-acting genetic variation, a much larger component of genetic effects influencing methylation are acting in trans. However, only 7 % of discovered mQTL are trans-effects, suggesting that the trans component is highly polygenic. Finally, we estimate the contribution of mQTL to variation in complex traits and infer that methylation may have a causal role consistent with an infinitesimal model in which many methylation sites each have a small influence, amounting to a large overall contribution.
DNA methylation contains a significant heritable component that remains consistent across the lifespan. Our results suggest that the genetic component of methylation may have a causal role in complex traits. The database of mQTL presented here provide a rich resource for those interested in investigating the role of methylation in disease.
Maternal smoking during pregnancy has been found to influence newborn DNA methylation in genes involved in fundamental developmental processes. It is pertinent to understand the degree to which the ...offspring methylome is sensitive to the intensity and duration of prenatal smoking. An investigation of the persistence of offspring methylation associated with maternal smoking and the relative roles of the intrauterine and postnatal environment is also warranted. In the Avon Longitudinal Study of Parents and Children, we investigated associations between prenatal exposure to maternal smoking and offspring DNA methylation at multiple time points in approximately 800 mother-offspring pairs. In cord blood, methylation at 15 CpG sites in seven gene regions (AHRR, MYO1G, GFI1, CYP1A1, CNTNAP2, KLF13 and ATP9A) was associated with maternal smoking, and a dose-dependent response was observed in relation to smoking duration and intensity. Longitudinal analysis of blood DNA methylation in serial samples at birth, age 7 and 17 years demonstrated that some CpG sites showed reversibility of methylation (GFI1, KLF13 and ATP9A), whereas others showed persistently perturbed patterns (AHRR, MYO1G, CYP1A1 and CNTNAP2). Of those showing persistence, we explored the effect of postnatal smoke exposure and found that the major contribution to altered methylation was attributed to a critical window of in utero exposure. A comparison of paternal and maternal smoking and offspring methylation showed consistently stronger maternal associations, providing further evidence for causal intrauterine mechanisms. These findings emphasize the sensitivity of the methylome to maternal smoking during early development and the long-term impact of such exposure.
Evidence suggests that in utero exposure to undernutrition and overnutrition might affect adiposity in later life. Epigenetic modification is suggested as a plausible mediating mechanism.
We used ...multivariable linear regression and a negative control design to examine offspring epigenome-wide DNA methylation in relation to maternal and offspring adiposity in 1018 participants.
Compared with neonatal offspring of normal weight mothers, 28 and 1621 CpG sites were differentially methylated in offspring of obese and underweight mothers, respectively false discovert rate (FDR)-corrected P-value < 0.05), with no overlap in the sites that maternal obesity and underweight relate to. A positive association, where higher methylation is associated with a body mass index (BMI) outside the normal range, was seen at 78.6% of the sites associated with obesity and 87.9% of the sites associated with underweight. Associations of maternal obesity with offspring methylation were stronger than associations of paternal obesity, supporting an intrauterine mechanism. There were no consistent associations of gestational weight gain with offspring DNA methylation. In general, sites that were hypermethylated in association with maternal obesity or hypomethylated in association with maternal underweight tended to be positively associated with offspring adiposity, and sites hypomethylated in association with maternal obesity or hypermethylated in association with maternal underweight tended to be inversely associated with offspring adiposity.
Our data suggest that both maternal obesity and, to a larger degree, underweight affect the neonatal epigenome via an intrauterine mechanism, but weight gain during pregnancy has little effect. We found some evidence that associations of maternal underweight with lower offspring adiposity and maternal obesity with greater offspring adiposity may be mediated via increased DNA methylation.
Epigenetic markings acquired in early life may have phenotypic consequences later in development through their role in transcriptional regulation with relevance to the developmental origins of ...diseases including obesity. The goal of this study was to investigate whether DNA methylation levels at birth are associated with body size later in childhood.
A study design involving two birth cohorts was used to conduct transcription profiling followed by DNA methylation analysis in peripheral blood. Gene expression analysis was undertaken in 24 individuals whose biological samples and clinical data were collected at a mean ± standard deviation (SD) age of 12.35 (0.95) years, the upper and lower tertiles of body mass index (BMI) were compared with a mean (SD) BMI difference of 9.86 (2.37) kg/m(2). This generated a panel of differentially expressed genes for DNA methylation analysis which was then undertaken in cord blood DNA in 178 individuals with body composition data prospectively collected at a mean (SD) age of 9.83 (0.23) years. Twenty-nine differentially expressed genes (>1.2-fold and p<10(-4)) were analysed to determine DNA methylation levels at 1-3 sites per gene. Five genes were unmethylated and DNA methylation in the remaining 24 genes was analysed using linear regression with bootstrapping. Methylation in 9 of the 24 (37.5%) genes studied was associated with at least one index of body composition (BMI, fat mass, lean mass, height) at age 9 years, although only one of these associations remained after correction for multiple testing (ALPL with height, p(Corrected) = 0.017).
DNA methylation patterns in cord blood show some association with altered gene expression, body size and composition in childhood. The observed relationship is correlative and despite suggestion of a mechanistic epigenetic link between in utero life and later phenotype, further investigation is required to establish causality.
Allergic disease is very common and carries substantial public-health burdens. We conducted a meta-analysis of genome-wide associations with self-reported cat, dust-mite and pollen allergies in ...53,862 individuals. We used generalized estimating equations to model shared and allergy-specific genetic effects. We identified 16 shared susceptibility loci with association P<5×10(-8), including 8 loci previously associated with asthma, as well as 4p14 near TLR1, TLR6 and TLR10 (rs2101521, P=5.3×10(-21)); 6p21.33 near HLA-C and MICA (rs9266772, P=3.2×10(-12)); 5p13.1 near PTGER4 (rs7720838, P=8.2×10(-11)); 2q33.1 in PLCL1 (rs10497813, P=6.1×10(-10)), 3q28 in LPP (rs9860547, P=1.2×10(-9)); 20q13.2 in NFATC2 (rs6021270, P=6.9×10(-9)), 4q27 in ADAD1 (rs17388568, P=3.9×10(-8)); and 14q21.1 near FOXA1 and TTC6 (rs1998359, P=4.8×10(-8)). We identified one locus with substantial evidence of differences in effects across allergies at 6p21.32 in the class II human leukocyte antigen (HLA) region (rs17533090, P=1.7×10(-12)), which was strongly associated with cat allergy. Our study sheds new light on the shared etiology of immune and autoimmune disease.
Statistical models that use an individual's DNA methylation levels to estimate their age (known as epigenetic clocks) have recently been developed, with 96% correlation found between epigenetic and ...chronological age. We postulate that differences between estimated and actual age age acceleration (AA) can be used as a measure of developmental age in early life.
We obtained DNA methylation measures at three time points (birth, age 7 years and age 17 years) in 1018 children from the Avon Longitudinal Study of Parents and Children (ALSPAC). Using an online calculator, we estimated epigenetic age, and thus AA, for each child at each time point. We then investigated whether AA was prospectively associated with repeated measures of height, weight, body mass index (BMI), bone mineral density, bone mass, fat mass, lean mass and Tanner stage.
Positive AA at birth was associated with higher average fat mass 1321 g per year of AA, 95% confidence interval (CI) 386, 2256 g from birth to adolescence (i.e. from age 0-17 years) and AA at age 7 was associated with higher average height (0.23 cm per year of AA, 95% CI 0.04, 0.41 cm). Conflicting evidence for the role of AA (at birth and in childhood) on changes during development was also found, with higher AA being positively associated with changes in weight, BMI and Tanner stage, but negatively with changes in height and fat mass.
We found evidence that being ahead of one's epigenetic age acceleration is related to developmental characteristics during childhood and adolescence. This demonstrates the potential for using AA as a measure of development in future research.
Humans display structural and functional asymmetries in brain organization, strikingly with respect to language and handedness. The molecular basis of these asymmetries is unknown. We report a ...genome-wide association study meta-analysis for a quantitative measure of relative hand skill in individuals with dyslexia reading disability (RD) (n = 728). The most strongly associated variant, rs7182874 (P = 8.68 × 10(-9)), is located in PCSK6, further supporting an association we previously reported. We also confirmed the specificity of this association in individuals with RD; the same locus was not associated with relative hand skill in a general population cohort (n = 2,666). As PCSK6 is known to regulate NODAL in the development of left/right (LR) asymmetry in mice, we developed a novel approach to GWAS pathway analysis, using gene-set enrichment to test for an over-representation of highly associated variants within the orthologs of genes whose disruption in mice yields LR asymmetry phenotypes. Four out of 15 LR asymmetry phenotypes showed an over-representation (FDR ≤ 5%). We replicated three of these phenotypes; situs inversus, heterotaxia, and double outlet right ventricle, in the general population cohort (FDR ≤ 5%). Our findings lead us to propose that handedness is a polygenic trait controlled in part by the molecular mechanisms that establish LR body asymmetry early in development.
Craniofacial morphology is highly heritable, but little is known about which genetic variants influence normal facial variation in the general population. We aimed to identify genetic variants ...associated with normal facial variation in a population-based cohort of 15-year-olds from the Avon Longitudinal Study of Parents and Children. 3D high-resolution images were obtained with two laser scanners, these were merged and aligned, and 22 landmarks were identified and their x, y, and z coordinates used to generate 54 3D distances reflecting facial features. 14 principal components (PCs) were also generated from the landmark locations. We carried out genome-wide association analyses of these distances and PCs in 2,185 adolescents and attempted to replicate any significant associations in a further 1,622 participants. In the discovery analysis no associations were observed with the PCs, but we identified four associations with the distances, and one of these, the association between rs7559271 in PAX3 and the nasion to midendocanthion distance (n-men), was replicated (p = 4 × 10−7). In a combined analysis, each G allele of rs7559271 was associated with an increase in n-men distance of 0.39 mm (p = 4 × 10−16), explaining 1.3% of the variance. Independent associations were observed in both the z (nasion prominence) and y (nasion height) dimensions (p = 9 × 10−9 and p = 9 × 10−10, respectively), suggesting that the locus primarily influences growth in the yz plane. Rare variants in PAX3 are known to cause Waardenburg syndrome, which involves deafness, pigmentary abnormalities, and facial characteristics including a broad nasal bridge. Our findings show that common variants within this gene also influence normal craniofacial development.
Systematic reviews of randomised controlled trials (RCTs) have suggested that maternal vitamin D (25OHD) and calcium supplementation increase birth weight. However, limitations of many trials were ...highlighted in the reviews. Our aim was to combine genetic and RCT data to estimate causal effects of these two maternal traits on offspring birth weight.
We performed two-sample mendelian randomisation (MR) using genetic instrumental variables associated with 25(OH)D and calcium that had been identified in genome-wide association studies (GWAS; sample 1; N = 122,123 for 25OHD and N = 61,275 for calcium). Associations between these maternal genetic variants and offspring birth weight were calculated in the UK Biobank (UKB) (sample 2; N = 190,406). We used data on mother-child pairs from two United Kingdom birth cohorts (combined N = 5,223) in sensitivity analyses to check whether results were influenced by fetal genotype, which is correlated with the maternal genotype (r ≈ 0.5). Further sensitivity analyses to test the reliability of the results included MR-Egger, weighted-median estimator, 'leave-one-out', and multivariable MR analyses. We triangulated MR results with those from RCTs, in which we used randomisation to supplementation with vitamin D (24 RCTs, combined N = 5,276) and calcium (6 RCTs, combined N = 543) as an instrumental variable to determine the effects of 25(OH)D and calcium on birth weight. In the main MR analysis, there was no strong evidence of an effect of maternal 25(OH)D on birth weight (difference in mean birth weight -0.03 g 95% CI -2.48 to 2.42 g, p = 0.981 per 10% higher maternal 25OHD). The effect estimate was consistent across our MR sensitivity analyses. Instrumental variable analyses applied to RCTs suggested a weak positive causal effect (5.94 g 95% CI 2.15-9.73, p = 0.002 per 10% higher maternal 25OHD), but this result may be exaggerated because of risk of bias in the included RCTs. The main MR analysis for maternal calcium also suggested no strong evidence of an effect on birth weight (-20 g 95% CI -44 to 5 g, p = 0.116 per 1 SD higher maternal calcium level). Some sensitivity analyses suggested that the genetic instrument for calcium was associated with birth weight via exposures that are independent of calcium levels (horizontal pleiotropy). Application of instrumental variable analyses to RCTs suggested that calcium has a substantial effect on birth weight (178 g 95% CI 121-236 g, p = 1.43 × 10-9 per 1 SD higher maternal calcium level) that was not consistent with any of the MR results. However, the RCT instrumental variable estimate may have been exaggerated because of risk of bias in the included RCTs. Other study limitations include the low response rate of UK Biobank, which may bias MR estimates, and the lack of suitable data to test whether the effects of genetic instruments on maternal calcium levels during pregnancy were the same as those outside of pregnancy.
Our results suggest that maternal circulating 25(OH)D does not influence birth weight in otherwise healthy newborns. However, the effect of maternal circulating calcium on birth weight is unclear and requires further exploration with more research including RCT and/or MR analyses with more valid instruments.