Collapsed methylation quantitative trait loci analysis for low frequency and rare variants.

BACKGROUND: Single variant approaches have been successful in identifying DNA methylation quantitative trait loci (mQTL), although as with complex traits they lack the statistical power to identify the effects from rare genetic variants. We have undertaken extensive analyses to identify regions of low frequency and rare variants that are associated with DNA methylation levels. METHODS: We used repeated measurements of DNA methylation from five different life stages in human blood, taken from the Avon Longitudinal Study of Parents and Children (ALSPAC) cohort. Variants were collapsed across CpG islands and their flanking regions to identify variants collectively associated with methylation, where no single variant was individually responsible for the observed signal. All analyses were undertaken using the sequence kernel association test. RESULTS: For loci where no individual variant mQTL was observed based on a single variant analysis, we identified 95 unique regions where the combined effect of low frequency variants (MAF ≤ 5%) provided strong evidence of association with methylation. For loci where there was previous evidence of an individual variant mQTL, a further 3 regions provided evidence of association between multiple low frequency variants and methylation levels. Effects were observed consistently across 5 different time points in the lifecourse and evidence of replication in the TwinsUK and Exeter cohorts was also identified. CONCLUSION: We have demonstrated the potential of this novel approach to mQTL analysis by analysing the combined effect of multiple low frequency or rare variants. Future studies should benefit from applying this approach as a complementary follow up to single variant analyses.

Prenatal and early life influences on epigenetic age in children: a study of mother-offspring pairs from two cohort studies.

DNA methylation-based biomarkers of aging are highly correlated with actual age. Departures of methylation-estimated age from actual age can be used to define epigenetic measures of child development or age acceleration (AA) in adults. Very little is known about genetic or environmental determinants of these epigenetic measures of aging. We obtained DNA methylation profiles using Infinium HumanMethylation450 BeadChips across five time-points in 1018 mother-child pairs from the Avon Longitudinal Study of Parents and Children. Using the Horvath age estimation method, we calculated epigenetic age for these samples. AA was defined as the residuals from regressing epigenetic age on actual age. AA was tested for associations with cross-sectional clinical variables in children. We identified associations between AA and sex, birth weight, birth by caesarean section and several maternal characteristics in pregnancy, namely smoking, weight, BMI, selenium and cholesterol level. Offspring of non-drinkers had higher AA on average but this difference appeared to resolve during childhood. The associations between sex, birth weight and AA found in ARIES were replicated in an independent cohort (GOYA). In children, epigenetic AA measures are associated with several clinically relevant variables, and early life exposures appear to be associated with changes in AA during adolescence. Further research into epigenetic aging, including the use of causal inference methods, is required to better our understanding of aging.

Identification of loci where DNA methylation potentially mediates genetic risk of type 1 diabetes.

The risk of Type 1 Diabetes (T1D) comprises both genetic and environmental components. We investigated whether genetic susceptibility to T1D could be mediated by changes in DNA methylation, an epigenetic mechanism that potentially plays a role in autoimmune diabetes. From enrichment analysis, we found that there was a common genetic influence for both DNA methylation and T1D across the genome, implying that methylation could be either on the causal pathway to T1D or a non-causal biomarker of T1D genetic risk. Using data from a general population comprising blood samples taken at birth (n = 844), childhood (n = 846) and adolescence (n = 907), we then evaluated the associations between 64 top GWAS single nucleotide polymorphisms (SNPs) and DNA methylation levels at 55 non-HLA loci. We identified 95 proximal SNP-cytosine phosphate guanine (CpG) pairs (cis) and 1 distal SNP-CpG association (trans) consistently at birth, childhood, and adolescence. Combining genetic co-localization and Mendelian Randomization analysis, we provided evidence that at 5 loci, ITGB3BP, AFF3, PTPN2, CTSH and CTLA4, DNA methylation is potentially mediating the genetic risk of T1D mainly by influencing local gene expression.

Longitudinal analysis of DNA methylation associated with birth weight and gestational age.

Gestational age (GA) and birth weight have been implicated in the determination of long-term health. It has been hypothesized that changes in DNA methylation may mediate these long-term effects. We obtained DNA methylation profiles from cord blood and peripheral blood at ages 7 and 17 in the same children from the Avon Longitudinal Study of Parents and Children. Repeated-measures data were used to investigate changes in birth-related methylation during childhood and adolescence. Ten developmental phenotypes (e.g. height) were analysed to identify possible mediation of health effects by DNA methylation. In cord blood, methylation at 224 CpG sites was found to be associated with GA and 23 CpG sites with birth weight. Methylation changed in the majority of these sites over time, but neither birth characteristic was strongly associated with methylation at age 7 or 17 (using a conservative correction for multiple testing of P < 1.03 × 10(-7)), suggesting resolution of differential methylation by early childhood. Associations were observed between birth weight-associated CpG sites and phenotypic characteristics in childhood. One strong association involved birth weight, methylation of a CpG site proximal to the NFIX locus and bone mineral density at age 17. Analysis of serial methylation from birth to adolescence provided evidence for a lack of persistence of methylation differences beyond early childhood. Sites associated with birth weight were linked to developmental genes and have methylation levels which are associated with developmental phenotypes. Replication and interrogation of causal relationships are needed to substantiate whether methylation differences at birth influence the association between birth weight and development.

Prenatal exposure to maternal smoking and offspring DNA methylation across the lifecourse: findings from the Avon Longitudinal Study of Parents and Children (ALSPAC).

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.

Sex-stratified genome-wide association studies including 270,000 individuals show sexual dimorphism in genetic loci for anthropometric traits.

Given the anthropometric differences between men and women and previous evidence of sex-difference in genetic effects, we conducted a genome-wide search for sexually dimorphic associations with height, weight, body mass index, waist circumference, hip circumference, and waist-to-hip-ratio (133,723 individuals) and took forward 348 SNPs into follow-up (additional 137,052 individuals) in a total of 94 studies. Seven loci displayed significant sex-difference (FDR<5%), including four previously established (near GRB14/COBLL1, LYPLAL1/SLC30A10, VEGFA, ADAMTS9) and three novel anthropometric trait loci (near MAP3K1, HSD17B4, PPARG), all of which were genome-wide significant in women (P<5×10(-8)), but not in men. Sex-differences were apparent only for waist phenotypes, not for height, weight, BMI, or hip circumference. Moreover, we found no evidence for genetic effects with opposite directions in men versus women. The PPARG locus is of specific interest due to its role in diabetes genetics and therapy. Our results demonstrate the value of sex-specific GWAS to unravel the sexually dimorphic genetic underpinning of complex traits.

Genome-wide joint meta-analysis of SNP and SNP-by-smoking interaction identifies novel loci for pulmonary function.

Genome-wide association studies have identified numerous genetic loci for spirometic measures of pulmonary function, forced expiratory volume in one second (FEV(1)), and its ratio to forced vital capacity (FEV(1)/FVC). Given that cigarette smoking adversely affects pulmonary function, we conducted genome-wide joint meta-analyses (JMA) of single nucleotide polymorphism (SNP) and SNP-by-smoking (ever-smoking or pack-years) associations on FEV(1) and FEV(1)/FVC across 19 studies (total N = 50,047). We identified three novel loci not previously associated with pulmonary function. SNPs in or near DNER (smallest P(JMA = )5.00×10(-11)), HLA-DQB1 and HLA-DQA2 (smallest P(JMA = )4.35×10(-9)), and KCNJ2 and SOX9 (smallest P(JMA = )1.28×10(-8)) were associated with FEV(1)/FVC or FEV(1) in meta-analysis models including SNP main effects, smoking main effects, and SNP-by-smoking (ever-smoking or pack-years) interaction. The HLA region has been widely implicated for autoimmune and lung phenotypes, unlike the other novel loci, which have not been widely implicated. We evaluated DNER, KCNJ2, and SOX9 and found them to be expressed in human lung tissue. DNER and SOX9 further showed evidence of differential expression in human airway epithelium in smokers compared to non-smokers. Our findings demonstrated that joint testing of SNP and SNP-by-environment interaction identified novel loci associated with complex traits that are missed when considering only the genetic main effects.

Association between variants of PRDM1 and NDP52 and Crohn's disease, based on exome sequencing and functional studies.

BACKGROUND & AIMS: Genome-wide association studies (GWAS) have identified 140 Crohn's disease (CD) susceptibility loci. For most loci, the variants that cause disease are not known and the genes affected by these variants have not been identified. We aimed to identify variants that cause CD through detailed sequencing, genetic association, expression, and functional studies. METHODS: We sequenced whole exomes of 42 unrelated subjects with CD and 5 healthy subjects (controls) and then filtered single nucleotide variants by incorporating association results from meta-analyses of CD GWAS and in silico mutation effect prediction algorithms. We then genotyped 9348 subjects with CD, 2868 subjects with ulcerative colitis, and 14,567 control subjects and associated variants analyzed in functional studies using materials from subjects and controls and in vitro model systems. RESULTS: We identified rare missense mutations in PR domain-containing 1 (PRDM1) and associated these with CD. These mutations increased proliferation of T cells and secretion of cytokines on activation and increased expression of the adhesion molecule L-selectin. A common CD risk allele, identified in GWAS, correlated with reduced expression of PRDM1 in ileal biopsy specimens and peripheral blood mononuclear cells (combined P = 1.6 × 10(-8)). We identified an association between CD and a common missense variant, Val248Ala, in nuclear domain 10 protein 52 (NDP52) (P = 4.83 × 10(-9)). We found that this variant impairs the regulatory functions of NDP52 to inhibit nuclear factor κB activation of genes that regulate inflammation and affect the stability of proteins in Toll-like receptor pathways. CONCLUSIONS: We have extended the results of GWAS and provide evidence that variants in PRDM1 and NDP52 determine susceptibility to CD. PRDM1 maps adjacent to a CD interval identified in GWAS and encodes a transcription factor expressed by T and B cells. NDP52 is an adaptor protein that functions in selective autophagy of intracellular bacteria and signaling molecules, supporting the role of autophagy in the pathogenesis of CD.

Genome-wide association studies identify CHRNA5/3 and HTR4 in the development of airflow obstruction.

RATIONALE: Genome-wide association studies (GWAS) have identified loci influencing lung function, but fewer genes influencing chronic obstructive pulmonary disease (COPD) are known. OBJECTIVES: Perform meta-analyses of GWAS for airflow obstruction, a key pathophysiologic characteristic of COPD assessed by spirometry, in population-based cohorts examining all participants, ever smokers, never smokers, asthma-free participants, and more severe cases. METHODS: Fifteen cohorts were studied for discovery (3,368 affected; 29,507 unaffected), and a population-based family study and a meta-analysis of case-control studies were used for replication and regional follow-up (3,837 cases; 4,479 control subjects). Airflow obstruction was defined as FEV(1) and its ratio to FVC (FEV(1)/FVC) both less than their respective lower limits of normal as determined by published reference equations. MEASUREMENTS AND MAIN RESULTS: The discovery meta-analyses identified one region on chromosome 15q25.1 meeting genome-wide significance in ever smokers that includes AGPHD1, IREB2, and CHRNA5/CHRNA3 genes. The region was also modestly associated among never smokers. Gene expression studies confirmed the presence of CHRNA5/3 in lung, airway smooth muscle, and bronchial epithelial cells. A single-nucleotide polymorphism in HTR4, a gene previously related to FEV(1)/FVC, achieved genome-wide statistical significance in combined meta-analysis. Top single-nucleotide polymorphisms in ADAM19, RARB, PPAP2B, and ADAMTS19 were nominally replicated in the COPD meta-analysis. CONCLUSIONS: These results suggest an important role for the CHRNA5/3 region as a genetic risk factor for airflow obstruction that may be independent of smoking and implicate the HTR4 gene in the etiology of airflow obstruction.

Genetic predictors of fibrin D-dimer levels in healthy adults.

BACKGROUND: Fibrin fragment D-dimer, one of several peptides produced when crosslinked fibrin is degraded by plasmin, is the most widely used clinical marker of activated blood coagulation. To identity genetic loci influencing D-dimer levels, we performed the first large-scale, genome-wide association search. METHODS AND RESULTS: A genome-wide investigation of the genomic correlates of plasma D-dimer levels was conducted among 21 052 European-ancestry adults. Plasma levels of D-dimer were measured independently in each of 13 cohorts. Each study analyzed the association between ≈2.6 million genotyped and imputed variants across the 22 autosomal chromosomes and natural-log­transformed D-dimer levels using linear regression in additive genetic models adjusted for age and sex. Among all variants, 74 exceeded the genome-wide significance threshold and marked 3 regions. At 1p22, rs12029080 (P=6.4×10(-52)) was 46.0 kb upstream from F3, coagulation factor III (tissue factor). At 1q24, rs6687813 (P=2.4×10(-14)) was 79.7 kb downstream of F5, coagulation factor V. At 4q32, rs13109457 (P=2.9×10(-18)) was located between 2 fibrinogen genes: 10.4 kb downstream from FGG and 3.0 kb upstream from FGA. Variants were associated with a 0.099-, 0.096-, and 0.061-unit difference, respectively, in natural-log­transformed D-dimer and together accounted for 1.8% of the total variance. When adjusted for nonsynonymous substitutions in F5 and FGA loci known to be associated with D-dimer levels, there was no evidence of an additional association at either locus. CONCLUSIONS: Three genes were associated with fibrin D-dimer levels. Of these 3, the F3 association was the strongest, and has not been previously reported.

Molecular and population analysis of natural selection on the human haptoglobin duplication.

Haptoglobin binds free haemoglobin that prevents oxidative damage produced by haemolysis. There is a copy number variant (CNV) in the haptoglobin gene (HP) consisting of two alleles, Hp1 (no duplication), and Hp2 (1.7kb duplication involving two exons). The spread of the Hp2 allele is believed to have taken place under selective pressures conferred by malaria resistance. However, molecular evidence is lacking and Hp did not emerge in genomewide SNPs surveys for evidence of selection. In Europe, there is geographical constancy of Hp2 frequency, indicative of absence of clinal pressures and that modern day European alleles represent a "snapshot" of their out-of-Africa migrations. In this work we test for signatures of natural selection acting on the Hp CNV in a sample from the UK population (Avon Longitudinal Study of Parents and Children, ALSPAC). We present here heterozygosity decay, pairwise F(ST) values observed between ALSPAC and 301 populations from all five populated continents, extended haplotype homozygosity analyses involving the CNV and 80 SNPs surrounding the CNV ~500kb in each direction, and linkage disequilibrium and pairwise haplotypic analyses involving 160 SNPs on chromosome 16q22.1. Taken together, our results represent the first molecular analysis of natural selection in the Hp CNV genetic region.

Meta-analysis of genome-wide scans for human adult stature identifies novel Loci and associations with measures of skeletal frame size.

Recent genome-wide (GW) scans have identified several independent loci affecting human stature, but their contribution through the different skeletal components of height is still poorly understood. We carried out a genome-wide scan in 12,611 participants, followed by replication in an additional 7,187 individuals, and identified 17 genomic regions with GW-significant association with height. Of these, two are entirely novel (rs11809207 in CATSPER4, combined P-value = 6.1x10(-8) and rs910316 in TMED10, P-value = 1.4x10(-7)) and two had previously been described with weak statistical support (rs10472828 in NPR3, P-value = 3x10(-7) and rs849141 in JAZF1, P-value = 3.2x10(-11)). One locus (rs1182188 at GNA12) identifies the first height eQTL. We also assessed the contribution of height loci to the upper- (trunk) and lower-body (hip axis and femur) skeletal components of height. We find evidence for several loci associated with trunk length (including rs6570507 in GPR126, P-value = 4x10(-5) and rs6817306 in LCORL, P-value = 4x10(-4)), hip axis length (including rs6830062 at LCORL, P-value = 4.8x10(-4) and rs4911494 at UQCC, P-value = 1.9x10(-4)), and femur length (including rs710841 at PRKG2, P-value = 2.4x10(-5) and rs10946808 at HIST1H1D, P-value = 6.4x10(-6)). Finally, we used conditional analyses to explore a possible differential contribution of the height loci to these different skeletal size measurements. In addition to validating four novel loci controlling adult stature, our study represents the first effort to assess the contribution of genetic loci to three skeletal components of height. Further statistical tests in larger numbers of individuals will be required to verify if the height loci affect height preferentially through these subcomponents of height.

A genome-wide meta-analysis of genetic variants associated with allergic rhinitis and grass sensitization and their interaction with birth order.

BACKGROUND: Hay fever or seasonal allergic rhinitis (AR) is a chronic disorder associated with IgE sensitization to grass. The underlying genetic variants have not been studied comprehensively. There is overwhelming evidence that those who have older siblings have less AR, although the mechanism for this remains unclear. OBJECTIVE: We sought to identify common genetic variant associations with prevalent AR and grass sensitization using existing genome-wide association study (GWAS) data and to determine whether genetic variants modify the protective effect of older siblings. METHOD: Approximately 2.2 million genotyped or imputed single nucleotide polymorphisms were investigated in 4 large European adult cohorts for AR (3,933 self-reported cases vs 8,965 control subjects) and grass sensitization (2,315 cases vs 10,032 control subjects). RESULTS: Three loci reached genome-wide significance for either phenotype. The HLA variant rs7775228, which cis-regulates HLA-DRB4, was strongly associated with grass sensitization and weakly with AR (P(grass) = 1.6 × 10(-9); P(AR) = 8.0 × 10(-3)). Variants in a locus near chromosome 11 open reading frame 30 (C11orf30) and leucine-rich repeat containing 32 (LRRC32), which was previously associated with atopic dermatitis and eczema, were also strongly associated with both phenotypes (rs2155219; P(grass) = 9.4 × 10(-9); P(AR) = 3.8 × 10(-8)). The third genome-wide significant variant was rs17513503 (P(grass) = 1.2 × 10(-8); PAR = 7.4 × 10(-7)) which was located near transmembrane protein 232 (TMEM232) and solute carrier family 25, member 46 (SLC25A46). Twelve further loci with suggestive associations were also identified. Using a candidate gene approach, where we considered variants within 164 genes previously thought to be important, we found variants in 3 further genes that may be of interest: thymic stromal lymphopoietin (TSLP), Toll-like receptor 6 (TLR6) and nucleotide-binding oligomerization domain containing 1 (NOD1/CARD4). We found no evidence for variants that modified the effect of birth order on either phenotype. CONCLUSIONS: This relatively large meta-analysis of GWASs identified few loci associated with AR and grass sensitization. No birth order interaction was identified in the current analyses.

Novel associations of multiple genetic loci with plasma levels of factor VII, factor VIII, and von Willebrand factor: The CHARGE (Cohorts for Heart and Aging Research in Genome Epidemiology) Consortium.

BACKGROUND: Plasma levels of coagulation factors VII (FVII), VIII (FVIII), and von Willebrand factor (vWF) influence risk of hemorrhage and thrombosis. We conducted genome-wide association studies to identify new loci associated with plasma levels. METHODS AND RESULTS: The setting of the study included 5 community-based studies for discovery comprising 23 608 European-ancestry participants: Atherosclerosis Risk In Communities Study, Cardiovascular Health Study, British 1958 Birth Cohort, Framingham Heart Study, and Rotterdam Study. All subjects had genome-wide single-nucleotide polymorphism (SNP) scans and at least 1 phenotype measured: FVII activity/antigen, FVIII activity, and vWF antigen. Each study used its genotype data to impute to HapMap SNPs and independently conducted association analyses of hemostasis measures using an additive genetic model. Study findings were combined by meta-analysis. Replication was conducted in 7604 participants not in the discovery cohort. For FVII, 305 SNPs exceeded the genome-wide significance threshold of 5.0x10(-8) and comprised 5 loci on 5 chromosomes: 2p23 (smallest P value 6.2x10(-24)), 4q25 (3.6x10(-12)), 11q12 (2.0x10(-10)), 13q34 (9.0x10(-259)), and 20q11.2 (5.7x10(-37)). Loci were within or near genes, including 4 new candidate genes and F7 (13q34). For vWF, 400 SNPs exceeded the threshold and marked 8 loci on 6 chromosomes: 6q24 (1.2x10(-22)), 8p21 (1.3x10(-16)), 9q34 (<5.0x10(-324)), 12p13 (1.7x10(-32)), 12q23 (7.3x10(-10)), 12q24.3 (3.8x10(-11)), 14q32 (2.3x10(-10)), and 19p13.2 (1.3x10(-9)). All loci were within genes, including 6 new candidate genes, as well as ABO (9q34) and VWF (12p13). For FVIII, 5 loci were identified and overlapped vWF findings. Nine of the 10 new findings were replicated. CONCLUSIONS: New genetic associations were discovered outside previously known biological pathways and may point to novel prevention and treatment targets of hemostasis disorders.

Genomic and phenotypic insights from an atlas of genetic effects on DNA methylation.

Characterizing genetic influences on DNA methylation (DNAm) provides an opportunity to understand mechanisms underpinning gene regulation and disease. In the present study, we describe results of DNAm quantitative trait locus (mQTL) analyses on 32,851 participants, identifying genetic variants associated with DNAm at 420,509 DNAm sites in blood. We present a database of >270,000 independent mQTLs, of which 8.5% comprise long-range (trans) associations. Identified mQTL associations explain 15-17% of the additive genetic variance of DNAm. We show that the genetic architecture of DNAm levels is highly polygenic. Using shared genetic control between distal DNAm sites, we constructed networks, identifying 405 discrete genomic communities enriched for genomic annotations and complex traits. Shared genetic variants are associated with both DNAm levels and complex diseases, but only in a minority of cases do these associations reflect causal relationships from DNAm to trait or vice versa, indicating a more complex genotype-phenotype map than previously anticipated.

Beta2-adrenoceptor polymorphisms and asthma from childhood to middle age in the British 1958 birth cohort: a genetic association study.

BACKGROUND: Functionally relevant polymorphisms of the beta2-adrenoceptor gene (ADRB2) are common in white populations, but their contribution to the burden of airways disease in the population is uncertain. We aimed to relate the long-term prevalence of asthma or wheeze to functional coding region polymorphisms in the ADRB2 gene. METHODS: The British 1958 birth cohort consisted of all people born in Britain during a week in 1958. Asthma, wheezy bronchitis, and wheezing were ascertained by interview at ages 7, 11, 16, 23, 33, and 42 years, and lung function tests at 35 and 45 years. DNA samples from 8018 participants in the 45-year follow-up were genotyped for three coding variants in the ADRB2 gene. We extend the follow-up of this nationwide cohort by a further 10 years and relate asthma prevalence, prognosis, and lung function to functional coding region polymorphisms in the ADRB2 gene in the cohort members who contributed DNA samples. We also compared and combined our findings with those reaching significance in two previous meta-analyses. FINDINGS: Half the cohort (4105 of 8018) had some history of wheezing illness by age 42 years. Neither lifetime prevalence nor age at onset were related to ADRB2 coding variants. However, the common polymorphisms Arg16Gly (rs1042713, Arg 16 allele frequency 36.3%) and Gln27Glu (rs1042714, Glu 27 allele frequency 44.6%) were significantly associated with persistence of asthmatic symptoms from childhood to middle age. Among homozygotes for the Arg16-Gln27 haplotype at these loci, 19.3% (41 of 212) childhood wheezers had five or more wheezing episodes in the past year at age 42, compared with 11.9% (71 of 599) with no copy of this haplotype. However, only 3% of all frequent adult wheezing was statistically attributable to this haplotype. The less common Thr164Ile polymorphism (rs1800888, Ile allele frequency 1.5%) was not a major predictor of either frequency or prognosis of asthma. Our data do not support the findings of previous meta-analyses when considered in isolation or when combined with their contributory studies. INTERPRETATION: ADRB2 polymorphisms might predict a small component of the long-term prognosis in childhood asthma, but are not important determinants of asthma incidence or prevalence in the British population.

The epigenetic clock and physical development during childhood and adolescence: longitudinal analysis from a UK birth cohort.

Background: 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. Methods: 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. Results: 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. Conclusions: 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.

DNA Methylation and BMI: Investigating Identified Methylation Sites at HIF3A in a Causal Framework.

Multiple differentially methylated sites and regions associated with adiposity have now been identified in large-scale cross-sectional studies. We tested for replication of associations between previously identified CpG sites at HIF3A and adiposity in ∼1,000 mother-offspring pairs from the Avon Longitudinal Study of Parents and Children (ALSPAC). Availability of methylation and adiposity measures at multiple time points, as well as genetic data, allowed us to assess the temporal associations between adiposity and methylation and to make inferences regarding causality and directionality. Overall, our results were discordant with those expected if HIF3A methylation has a causal effect on BMI and provided more evidence for causality in the reverse direction (i.e., an effect of BMI on HIF3A methylation). These results are based on robust evidence from longitudinal analyses and were also partially supported by Mendelian randomization analysis, although this latter analysis was underpowered to detect a causal effect of BMI on HIF3A methylation. Our results also highlight an apparent long-lasting intergenerational influence of maternal BMI on offspring methylation at this locus, which may confound associations between own adiposity and HIF3A methylation. Further work is required to replicate and uncover the mechanisms underlying the direct and intergenerational effect of adiposity on DNA methylation.

Systematic identification of genetic influences on methylation across the human life course.

BACKGROUND: 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. RESULTS: 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. CONCLUSIONS: 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.