Genetic insights into resting heart rate and its role in cardiovascular disease

Resting heart rate is associated with cardiovascular diseases and mortality in observational and Mendelian randomization studies. The aims of this study are to extend the number of resting heart rate associated genetic variants and to obtain further insights in resting heart rate biology and its clinical consequences. A genome-wide meta-analysis of 100 studies in up to 835,465 individuals reveals 493 independent genetic variants in 352 loci, including 68 genetic variants outside previously identified resting heart rate associated loci. We prioritize 670 genes and in silico annotations point to their enrichment in cardiomyocytes and provide insights in their ECG signature. Two-sample Mendelian randomization analyses indicate that higher genetically predicted resting heart rate increases risk of dilated cardiomyopathy, but decreases risk of developing atrial fibrillation, ischemic stroke, and cardio-embolic stroke. We do not find evidence for a linear or non-linear genetic association between resting heart rate and all-cause mortality in contrast to our previous Mendelian randomization study. Systematic alteration of key differences between the current and previous Mendelian randomization study indicates that the most likely cause of the discrepancy between these studies arises from false positive findings in previous one-sample MR analyses caused by weak-instrument bias at lower P-value thresholds. The results extend our understanding of resting heart rate biology and give additional insights in its role in cardiovascular disease development.

Understanding the genetic complexity of puberty timing across the allele frequency spectrum

Pubertal timing varies considerably and has been associated with a range of health outcomes in later life. To elucidate the underlying biological mechanisms, we performed multi-ancestry genetic analyses in ~800,000 women, identifying 1,080 independent signals associated with age at menarche. Collectively these loci explained 11% of the trait variance in an independent sample, with women at the top and bottom 1% of polygenic risk exhibiting a ~11 and ~14-fold higher risk of delayed and precocious pubertal development, respectively. These common variant analyses were supported by exome sequence analysis of ~220,000 women, identifying several genes, including rare loss of function variants in ZNF483 which abolished the impact of polygenic risk. Next, we implicated 660 genes in pubertal development using a combination of in silico variant-to-gene mapping approaches and integration with dynamic gene expression data from mouse embryonic GnRH neurons. This included an uncharacterized G-protein coupled receptor GPR83, which we demonstrate amplifies signaling of MC3R, a key sensor of nutritional status. Finally, we identified several genes, including ovary-expressed genes involved in DNA damage response that co-localize with signals associated with menopause timing, leading us to hypothesize that the ovarian reserve might signal centrally to trigger puberty. Collectively these findings extend our understanding of the biological complexity of puberty timing and highlight body size dependent and independent mechanisms that potentially link reproductive timing to later life disease.

Genome-wide analysis identifies genetic effects on reproductive success and ongoing natural selection at the FADS locus

Identifying genetic determinants of reproductive success may highlight mechanisms underlying fertility and identify alleles under present-day selection. Using data in 785,604 individuals of European ancestry, we identified 43 genomic loci associated with either number of children ever born (NEB) or childlessness. These loci span diverse aspects of reproductive biology, including puberty timing, age at first birth, sex hormone regulation, endometriosis and age at menopause. Missense variants in ARHGAP27 were associated with higher NEB but shorter reproductive lifespan, suggesting a trade-off at this locus between reproductive ageing and intensity. Other genes implicated by coding variants include PIK3IP1, ZFP82 and LRP4, and our results suggest a new role for the melanocortin 1 receptor (MC1R) in reproductive biology. As NEB is one component of evolutionary fitness, our identified associations indicate loci under present-day natural selection. Integration with data from historical selection scans highlighted an allele in the FADS1/2 gene locus that has been under selection for thousands of years and remains so today. Collectively, our findings demonstrate that a broad range of biological mechanisms contribute to reproductive success.

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.

Genetic insights into biological mechanisms governing human ovarian ageing

Reproductive longevity is essential for fertility and influences healthy ageing in women1,2, but insights into its underlying biological mechanisms and treatments to preserve it are limited. Here we identify 290 genetic determinants of ovarian ageing, assessed using normal variation in age at natural menopause (ANM) in about 200,000 women of European ancestry. These common alleles were associated with clinical extremes of ANM; women in the top 1% of genetic susceptibility have an equivalent risk of premature ovarian insufficiency to those carrying monogenic FMR1 premutations3. The identified loci implicate a broad range of DNA damage response (DDR) processes and include loss-of-function variants in key DDR-associated genes. Integration with experimental models demonstrates that these DDR processes act across the life-course to shape the ovarian reserve and its rate of depletion. Furthermore, we demonstrate that experimental manipulation of DDR pathways highlighted by human genetics increases fertility and extends reproductive life in mice. Causal inference analyses using the identified genetic variants indicate that extending reproductive life in women improves bone health and reduces risk of type 2 diabetes, but increases the risk of hormone-sensitive cancers. These findings provide insight into the mechanisms that govern ovarian ageing, when they act, and how they might be targeted by therapeutic approaches to extend fertility and prevent disease.

Genome-wide association study in almost 195,000 individuals identifies 50 previously unidentified genetic loci for eye color

Human eye color is highly heritable, but its genetic architecture is not yet fully understood. We report the results of the largest genome-wide association study for eye color to date, involving up to 192,986 European participants from 10 populations. We identify 124 independent associations arising from 61 discrete genomic regions, including 50 previously unidentified. We find evidence for genes involved in melanin pigmentation, but we also find associations with genes involved in iris morphology and structure. Further analyses in 1636 Asian participants from two populations suggest that iris pigmentation variation in Asians is genetically similar to Europeans, albeit with smaller effect sizes. Our findings collectively explain 53.2% (95% confidence interval, 45.4 to 61.0%) of eye color variation using common single-nucleotide polymorphisms. Overall, our study outcomes demonstrate that the genetic complexity of human eye color considerably exceeds previous knowledge and expectations, highlighting eye color as a genetically highly complex human trait.

GWAS on longitudinal growth traits reveals different genetic factors influencing infant, child, and adult BMI

Early childhood growth patterns are associated with adult health, yet the genetic factors and the developmental stages involved are not fully understood. Here, we combine genome-wide association studies with modeling of longitudinal growth traits to study the genetics of infant and child growth, followed by functional, pathway, genetic correlation, risk score, and colocalization analyses to determine how developmental timings, molecular pathways, and genetic determinants of these traits overlap with those of adult health. We found a robust overlap between the genetics of child and adult body mass index (BMI), with variants associated with adult BMI acting as early as 4 to 6 years old. However, we demonstrated a completely distinct genetic makeup for peak BMI during infancy, influenced by variation at the LEPR/LEPROT locus. These findings suggest that different genetic factors control infant and child BMI. In light of the obesity epidemic, these findings are important to inform the timing and targets of prevention strategies.

Publisher Correction: Genome-wide association meta-analysis of individuals of European ancestry identifies new loci explaining a substantial fraction of hair color variation and heritability

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

Association of maternal circulating 25(OH)D and calcium with birth weight: A mendelian randomisation analysis

BACKGROUND

Systematic reviews of randomised controlled trials (RCTs) have suggested that maternal vitamin D (25[OH]D) 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.

METHODS AND FINDINGS

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 25[OH]D 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 25[OH]D). 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 25[OH]D), 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.

CONCLUSIONS

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.

Maternal and fetal genetic effects on birth weight and their relevance to cardio-metabolic risk factors

Birth weight variation is influenced by fetal and maternal genetic and non-genetic factors, and has been reproducibly associated with future cardio-metabolic health outcomes. In expanded genome-wide association analyses of own birth weight (n = 321,223) and offspring birth weight (n = 230,069 mothers), we identified 190 independent association signals (129 of which are novel). We used structural equation modeling to decompose the contributions of direct fetal and indirect maternal genetic effects, then applied Mendelian randomization to illuminate causal pathways. For example, both indirect maternal and direct fetal genetic effects drive the observational relationship between lower birth weight and higher later blood pressure: maternal blood pressure-raising alleles reduce offspring birth weight, but only direct fetal effects of these alleles, once inherited, increase later offspring blood pressure. Using maternal birth weight-lowering genotypes to proxy for an adverse intrauterine environment provided no evidence that it causally raises offspring blood pressure, indicating that the inverse birth weight-blood pressure association is attributable to genetic effects, and not to intrauterine programming.

Coronary artery disease, genetic risk and the metabolome in young individuals

Background: Genome-wide association studies have identified genetic variants associated with coronary artery disease (CAD) in adults - the leading cause of death worldwide. It often occurs later in life, but variants may impact CAD-relevant phenotypes early and throughout the life-course. Cohorts with longitudinal and genetic data on thousands of individuals are letting us explore the antecedents of this adult disease. Methods: 148 metabolites, with a focus on the lipidome, measured using nuclear magnetic resonance ( 1H-NMR) spectroscopy, and genotype data were available from 5,907 individuals at ages 7, 15, and 17 years from the Avon Longitudinal Study of Parents and Children (ALSPAC) cohort. Linear regression was used to assess the association between the metabolites and an adult-derived genetic risk score (GRS) of CAD comprising 146 variants. Individual variant-metabolite associations were also examined. Results: The CAD-GRS associated with 118 of 148 metabolites (false discovery rate [FDR] < 0.05), the strongest associations being with low-density lipoprotein (LDL) and atherogenic non-LDL subgroups. Nine of 146 variants in the GRS associated with one or more metabolites (FDR < 0.05). Seven of these are within lipid loci: rs11591147 PCSK9, rs12149545 HERPUD1-CETP, rs17091891 LPL, rs515135 APOB, rs602633 CELSR2-PSRC1, rs651821 APOA5, rs7412 APOE-APOC1. All associated with metabolites in the LDL or atherogenic non-LDL subgroups or both including aggregate cholesterol measures. The other two variants identified were rs112635299 SERPINA1 and rs2519093 ABO. Conclusions: Genetic variants that influence CAD risk in adults are associated with large perturbations in metabolite levels in individuals as young as seven. The variants identified are mostly within lipid-related loci and the metabolites they associated with are primarily linked to lipoproteins. Along with further research, this knowledge could allow for preventative measures, such as increased monitoring of at-risk individuals and perhaps treatment earlier in life, to be taken years before any symptoms of the disease arise.

A genome-wide association study of mitochondrial DNA copy number in two population-based cohorts

BACKGROUND

Mitochondrial DNA copy number (mtDNA CN) exhibits interindividual and intercellular variation, but few genome-wide association studies (GWAS) of directly assayed mtDNA CN exist. We undertook a GWAS of qPCR-assayed mtDNA CN in the Avon Longitudinal Study of Parents and Children (ALSPAC) and the UK Blood Service (UKBS) cohort. After validating and harmonising data, 5461 ALSPAC mothers (16-43 years at mtDNA CN assay) and 1338 UKBS females (17-69 years) were included in a meta-analysis. Sensitivity analyses restricted to females with white cell-extracted DNA and adjusted for estimated or assayed cell proportions. Associations were also explored in ALSPAC children and UKBS males.

RESULTS

A neutrophil-associated locus approached genome-wide significance (rs709591 [MED24], β (change in SD units of mtDNA CN per allele) [SE] - 0.084 [0.016], p = 1.54e-07) in the main meta-analysis of adult females. This association was concordant in magnitude and direction in UKBS males and ALSPAC neonates. SNPs in and around ABHD8 were associated with mtDNA CN in ALSPAC neonates (rs10424198, β [SE] 0.262 [0.034], p = 1.40e-14), but not other study groups. In a meta-analysis of unrelated individuals (N = 11,253), we replicated a published association in TFAM (β [SE] 0.046 [0.017], p = 0.006), with an effect size much smaller than that observed in the replication analysis of a previous in silico GWAS.

CONCLUSIONS

In a hypothesis-generating GWAS, we confirm an association between TFAM and mtDNA CN and present putative loci requiring replication in much larger samples. We discuss the limitations of our work, in terms of measurement error and cellular heterogeneity, and highlight the need for larger studies to better understand nuclear genomic control of mtDNA copy number.

Genome-wide association study of offspring birth weight in 86 577 women identifies five novel loci and highlights maternal genetic effects that are independent of fetal genetics

Genome-wide association studies of birth weight have focused on fetal genetics, whereas relatively little is known about the role of maternal genetic variation. We aimed to identify maternal genetic variants associated with birth weight that could highlight potentially relevant maternal determinants of fetal growth. We meta-analysed data on up to 8.7 million SNPs in up to 86 577 women of European descent from the Early Growth Genetics (EGG) Consortium and the UK Biobank. We used structural equation modelling (SEM) and analyses of mother–child pairs to quantify the separate maternal and fetal genetic effects. Maternal SNPs at 10 loci (MTNR1B, HMGA2, SH2B3, KCNAB1, L3MBTL3, GCK, EBF1, TCF7L2, ACTL9, CYP3A7) were associated with offspring birth weight at P < 5 × 10⁻⁸. In SEM analyses, at least 7 of the 10 associations were consistent with effects of the maternal genotype acting via the intrauterine environment, rather than via effects of shared alleles with the fetus. Variants, or correlated proxies, at many of the loci had been previously associated with adult traits, including fasting glucose (MTNR1B, GCK and TCF7L2) and sex hormone levels (CYP3A7), and one (EBF1) with gestational duration. The identified associations indicate that genetic effects on maternal glucose, cytochrome P450 activity and gestational duration, and potentially on maternal blood pressure and immune function, are relevant for fetal growth. Further characterization of these associations in mechanistic and causal analyses will enhance understanding of the potentially modifiable maternal determinants of fetal growth, with the goal of reducing the morbidity and mortality associated with low and high birth weights.

DNA methylation derived systemic inflammation indices are associated with head and neck cancer development and survival

OBJECTIVES

Head and neck squamous cell carcinoma (HNSCC) is often associated with chronic systemic inflammation (SI). In the present study, we assessed if DNA methylation-derived SI (mdSI) indices: Neutrophil-to-Lymphocyte ratio (mdNLR) and Lymphocyte-to-Monocyte ratio (mdLMR) are associated with the presence of HNSCC and overall survival (OS).

MATERIALS AND METHODS

We used two peripheral blood DNA methylation datasets: an HNSCC case-control dataset (n = 183) and an HNSCC survival dataset (n = 407) to estimate mdSI indices. We then performed multivariate regressions to test the association between mdSI indices, HNSCC development and OS.

RESULTS

Multivariate logistic regression revealed that elevated mdNLR was associated with increased odds of being an HNSCC case (OR = 3.25, 95% CI = 2.14-5.34, P = 4 × 10-7) while the converse was observed for mdLMR (OR = 0.88, 95% CI = 0.81-0.90, P = 2 × 10-3). In the HNSCC survival dataset, HPV16-E6 seropositive HNSCC cases had an elevated mdLMR (P = 9 × 10-5) and a lower mdNLR (P = 0.003) compared to seronegative patients. Multivariate Cox regression in the HNSCC survival dataset revealed that lower mdLMR (HR = 1.96, 95% CI = 1.30-2.95, P = 0.0013) but not lower mdNLR (HR = 0.68, 95% CI = 0.46-1.00, P = 0.0501) was associated with increased risk of death.

CONCLUSION

Our results indicate that mdSI estimated by DNA methylation data is associated with the presence of HNSCC and overall survival. The mdSI indices may be used as a valuable research tool to reliably estimate SI in the absence of cell-based estimates. Rigorous validation of our findings in large prospective studies is warranted in the future.

Coronary artery disease, genetic risk and the metabolome in young individuals

Background: Genome-wide association studies have identified genetic variants associated with coronary artery disease (CAD) in adults - the leading cause of death worldwide. It often occurs later in life, but variants may impact CAD-relevant phenotypes early and throughout the life-course. Cohorts with longitudinal and genetic data on thousands of individuals are letting us explore the antecedents of this adult disease. Methods: 149 metabolites, with a focus on the lipidome, measured using nuclear magnetic resonance ( 1H-NMR) spectroscopy, and genotype data were available from 5,905 individuals at ages 7, 15, and 17 years from the Avon Longitudinal Study of Parents and Children (ALSPAC) cohort. Linear regression was used to assess the association between the metabolites and an adult-derived genetic risk score (GRS) of CAD comprising 146 variants. Individual variant-metabolite associations were also examined. Results: The CAD-GRS associated with 118 of 149 metabolites (false discovery rate [FDR] < 0.05), the strongest associations being with low-density lipoprotein (LDL) and atherogenic non-LDL subgroups. Nine of 146 variants in the GRS associated with one or more metabolites (FDR < 0.05). Seven of these are within lipid loci: rs11591147 PCSK9, rs12149545 HERPUD1-CETP, rs17091891 LPL, rs515135 APOB, rs602633 CELSR2-PSRC1, rs651821 APOA5, rs7412 APOE-APOC1. All associated with metabolites in the LDL or atherogenic non-LDL subgroups or both including aggregate cholesterol measures. The other two variants identified were rs112635299 SERPINA1 and rs2519093 ABO. Conclusions: Genetic variants that influence CAD risk in adults are associated with large perturbations in metabolite levels in individuals as young as seven. The variants identified are mostly within lipid-related loci and the metabolites they associated with are primarily linked to lipoproteins. This knowledge could allow for preventative measures, such as increased monitoring of at-risk individuals and perhaps treatment earlier in life, to be taken years before any symptoms of the disease arise.

Genome-wide survey of parent-of-origin effects on DNA methylation identifies candidate imprinted loci in humans

Genomic imprinting is an epigenetic mechanism leading to parent-of-origin silencing of alleles. So far, the precise number of imprinted regions in humans is uncertain. In this study, we leveraged genome-wide DNA methylation in whole blood measured longitudinally at three time points (birth, childhood and adolescence) and genome-wide association studies (GWAS) data in 740 mother-child duos from the Avon Longitudinal Study of parents and children to identify candidate imprinted loci. We reasoned that cis-meQTLs at genomic regions that were imprinted would show strong evidence of parent-of-origin associations with DNA methylation, enabling the detection of imprinted regions. Using this approach, we identified genome-wide significant cis-meQTLs that exhibited parent-of-origin effects (POEs) at 82 loci, 34 novel and 48 regions previously implicated in imprinting (3.7-10<P < 10-300). Using an independent dataset from the Brisbane Systems Genetic Study, we replicated 76 out of the 82 identified loci. POEs were remarkably consistent across time points and were so strong at some loci that methylation levels enabled good discrimination of parental transmissions at these and surrounding genomic regions. The implication is that parental allelic transmissions could be modelled at many imprinted (and linked) loci in GWAS of unrelated individuals given a combination of genetic and methylation data. Novel regions showing parent of origin effects on methylation will require replication using a different technology and further functional experiments to confirm that such effects arise through a genomic imprinting mechanism.

Longitudinal serological measures of common infection in the Avon Longitudinal Study of Parents and Children cohort

Antibodies against pathogens provide information on exposure to infectious agents and are meaningful measures of past and present infection. Antibodies were measured in the plasma of children that are the offspring in a population-based birth cohort, the Avon Longitudinal Study of Parents and Children (ALSPAC). Plasma was collected during clinics at age 5, 7, 11 and 15 years. The antigens examined include: fungal ( Saccharomyces cerevisiae); protozoan ( Toxoplasma gondii and surface antigen 1 of  T. gondii); herpes viruses (cytomegalovirus, Epstein-Barr virus, herpes simplex virus type 1); common colds (influenza virus subtypes H1N1 and H3N2); other antigens (measles); animal (feline herpes virus, Theiler's virus); bacteria ( Helicobacter pylori); dietary antigens (bovine casein alpha protein, bovine casein beta protein). Alongside the depth of data available within the ALSPAC cohort, this longitudinal resource will enable the investigation of the association between infections and a wide variety of outcomes.

Longitudinal serological measures of common infection in the Avon Longitudinal Study of Parents and Children cohort

Antibodies against pathogens provide information on exposure to infectious agents and are meaningful measures of past and present infection. Antibodies were measured in the plasma of children that are the offspring in a population-based birth cohort, the Avon Longitudinal Study of Parents and Children (ALSPAC). Plasma was collected during clinics at age 5, 7, 11 and 15 years. The antigens examined include: fungal ( Saccharomyces cerevisiae); protozoan ( Toxoplasma gondii and surface antigen 1 of  T. gondii); herpes viruses (cytomegalovirus, Epstein-Barr virus, herpes simplex virus type 1); common colds (influenza virus subtypes H1N1 and H3N2); other antigens (measles); animal (feline herpes virus, Theiler’s virus); bacteria ( Helicobacter pylori); dietary antigens (bovine casein alpha protein, bovine casein beta protein). Alongside the depth of data available within the ALSPAC cohort, this longitudinal resource will enable the investigation of the association between infections and a wide variety of outcomes.

Genome-wide association meta-analysis of individuals of European ancestry identifies new loci explaining a substantial fraction of hair color variation and heritability

Hair color is one of the most recognizable visual traits in European populations and is under strong genetic control. Here we report the results of a genome-wide association study meta-analysis of almost 300,000 participants of European descent. We identified 123 autosomal and one X-chromosome loci significantly associated with hair color; all but 13 are novel. Collectively, single-nucleotide polymorphisms associated with hair color within these loci explain 34.6% of red hair, 24.8% of blond hair, and 26.1% of black hair heritability in the study populations. These results confirm the polygenic nature of complex phenotypes and improve our understanding of melanin pigment metabolism in humans.

Better governance, better access: practising responsible data sharing in the METADAC governance infrastructure

BACKGROUND

Genomic and biosocial research data about individuals is rapidly proliferating, bringing the potential for novel opportunities for data integration and use. The scale, pace and novelty of these applications raise a number of urgent sociotechnical, ethical and legal questions, including optimal methods of data storage, management and access. Although the open science movement advocates unfettered access to research data, many of the UK's longitudinal cohort studies operate systems of managed data access, in which access is governed by legal and ethical agreements between stewards of research datasets and researchers wishing to make use of them. Amongst other things, these agreements aim to respect the reasonable expectations of the research participants who provided data and samples, as expressed in the consent process. Arguably, responsible data management and governance of data and sample use are foundational to the consent process in longitudinal studies and are an important source of trustworthiness in the eyes of those who contribute data to genomic and biosocial research.

METHODS

This paper presents an ethnographic case study exploring the foundational principles of a governance infrastructure for Managing Ethico-social, Technical and Administrative issues in Data ACcess (METADAC), which are operationalised through a committee known as the METADAC Access Committee. METADAC governs access to phenotype, genotype and 'omic' data and samples from five UK longitudinal studies.

FINDINGS

Using the example of METADAC, we argue that three key structural features are foundational for practising responsible data sharing: independence and transparency; interdisciplinarity; and participant-centric decision-making. We observe that the international research community is proactively working towards optimising the use of research data, integrating/linking these data with routine data generated by health and social care services and other administrative data services to improve the analysis, interpretation and utility of these data. The governance of these new complex data assemblages will require a range of expertise from across a number of domains and disciplines, including that of study participants. Human-mediated decision-making bodies will be central to ensuring achievable, reasoned and responsible decisions about the use of these data; the METADAC model described in this paper provides an example of how this could be realised.

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.

Developmental Changes Within the Genetic Architecture of Social Communication Behavior: A Multivariate Study of Genetic Variance in Unrelated Individuals

BACKGROUND

Recent analyses of trait-disorder overlap suggest that psychiatric dimensions may relate to distinct sets of genes that exert maximum influence during different periods of development. This includes analyses of social communication difficulties that share, depending on their developmental stage, stronger genetic links with either autism spectrum disorder or schizophrenia. We developed a multivariate analysis framework in unrelated individuals to model directly the developmental profile of genetic influences contributing to complex traits, such as social communication difficulties, during an approximately 10-year period spanning childhood and adolescence.

METHODS

Longitudinally assessed quantitative social communication problems (N ≤ 5551) were studied in participants from a United Kingdom birth cohort (Avon Longitudinal Study of Parents and Children; age range, 8-17 years). Using standardized measures, genetic architectures were investigated with novel multivariate genetic-relationship-matrix structural equation models incorporating whole-genome genotyping information. Analogous to twin research, genetic-relationship-matrix structural equation models included Cholesky decomposition, common pathway, and independent pathway models.

RESULTS

A two-factor Cholesky decomposition model described the data best. One genetic factor was common to Social Communication Disorder Checklist measures across development; the other accounted for independent variation at 11 years and later, consistent with distinct developmental profiles in trait-disorder overlap. Importantly, genetic factors operating at 8 years explained only approximately 50% of genetic variation at 17 years.

CONCLUSIONS

Using latent factor models, we identified developmental changes in the genetic architecture of social communication difficulties that enhance the understanding of autism spectrum disorder- and schizophrenia-related dimensions. More generally, genetic-relationship-matrix structural equation models present a framework for modeling shared genetic etiologies between phenotypes and can provide prior information with respect to patterns and continuity of trait-disorder overlap.

Maternal and fetal genetic contribution to gestational weight gain

BACKGROUND

Clinical recommendations to limit gestational weight gain (GWG) imply high GWG is causally related to adverse outcomes in mother or offspring, but GWG is the sum of several inter-related complex phenotypes (maternal fat deposition and vascular expansion, placenta, amniotic fluid and fetal growth). Understanding the genetic contribution to GWG could help clarify the potential effect of its different components on maternal and offspring health. Here we explore the genetic contribution to total, early and late GWG.

PARTICIPANTS AND METHODS

A genome-wide association study was used to identify maternal and fetal variants contributing to GWG in up to 10 543 mothers and 16 317 offspring of European origin, with replication in 10 660 mothers and 7561 offspring. Additional analyses determined the proportion of variability in GWG from maternal and fetal common genetic variants and the overlap of established genome-wide significant variants for phenotypes relevant to GWG (for example, maternal body mass index (BMI) and glucose, birth weight).

RESULTS

Approximately 20% of the variability in GWG was tagged by common maternal genetic variants, and the fetal genome made a surprisingly minor contribution to explain variation in GWG. Variants near the pregnancy-specific beta-1 glycoprotein 5 (PSG5) gene reached genome-wide significance (P=1.71 × 10-8) for total GWG in the offspring genome, but did not replicate. Some established variants associated with increased BMI, fasting glucose and type 2 diabetes were associated with lower early, and higher later GWG. Maternal variants related to higher systolic blood pressure were related to lower late GWG. Established maternal and fetal birth weight variants were largely unrelated to GWG.

CONCLUSIONS

We found a modest contribution of maternal common variants to GWG and some overlap of maternal BMI, glucose and type 2 diabetes variants with GWG. These findings suggest that associations between GWG and later offspring/maternal outcomes may be due to the relationship of maternal BMI and diabetes with GWG.

Maternal thyroid function and child educational attainment: prospective cohort study

OBJECTIVE

To determine if first trimester maternal thyroid dysfunction is a critical determinant of child scholastic performance and overall educational attainment.

DESIGN

Prospective cohort study.

SETTING

Avon Longitudinal Study of Parents and Children cohort in the UK.

PARTICIPANTS

4615 mother-child pairs with an available first trimester sample (median 10 weeks gestation, interquartile range 8-12).

EXPOSURES

Free thyroxine, thyroid stimulating hormone, and thyroid peroxidase antibodies assessed as continuous measures and the seven clinical categories of maternal thyroid function.

MAIN OUTCOME MEASURES

Five age-specific national curriculum assessments in 3580 children at entry stage assessment at 54 months, increasing up to 4461 children at their final school assessment at age 15.

RESULTS

No strong evidence of clinically meaningful associations of first trimester free thyroxine and thyroid stimulating hormone levels with entry stage assessment score or Standard Assessment Test scores at any of the key stages was found. Associations of maternal free thyroxine or thyroid stimulating hormone with the total number of General Certificates of Secondary Education (GCSEs) passed (range 0-16) were all close to the null: free thyroxine, rate ratio per pmol/L 1.00 (95% confidence interval 1.00 to 1.01); and thyroid stimulating hormone, rate ratio 0.98 (0.94 to 1.02). No important relationship was observed when more detailed capped scores of GCSEs allowing for both the number and grade of pass or when language, mathematics, and science performance were examined individually or when all educational assessments undertaken by an individual from school entry to leaving were considered. 200 (4.3%) mothers were newly identified as having hypothyroidism or subclinical hypothyroidism and 97 (2.1%) subclinical hyperthyroidism or hyperthyroidism. Children of mothers with thyroid dysfunction attained an equivalent number of GCSEs and equivalent grades as children of mothers with euthyroidism.

CONCLUSIONS

Maternal thyroid dysfunction in early pregnancy does not have a clinically important association with impaired child performance at school or educational achievement.

Age at menarche and lung function: a Mendelian randomization study

A trend towards earlier menarche in women has been associated with childhood factors (e.g. obesity) and hypothesised environmental exposures (e.g. endocrine disruptors present in household products). Observational evidence has shown detrimental effects of early menarche on various health outcomes including adult lung function, but these might represent spurious associations due to confounding. To address this we used Mendelian randomization where genetic variants are used as proxies for age at menarche, since genetic associations are not affected by classical confounding. We estimated the effects of age at menarche on forced vital capacity (FVC), a proxy for restrictive lung impairment, and ratio of forced expiratory volume in one second to FVC (FEV1/FVC), a measure of airway obstruction, in both adulthood and adolescence. We derived SNP-age at menarche association estimates for 122 variants from a published genome-wide meta-analysis (N = 182,416), with SNP-lung function estimates obtained by meta-analysing three studies of adult women (N = 46,944) and two of adolescent girls (N = 3025). We investigated the impact of departures from the assumption of no pleiotropy through sensitivity analyses. In adult women, in line with previous evidence, we found an effect on restrictive lung impairment with a 24.8 mL increase in FVC per year increase in age at menarche (95% CI 1.8-47.9; p = 0.035); evidence was stronger after excluding potential pleiotropic variants (43.6 mL; 17.2-69.9; p = 0.001). In adolescent girls we found an opposite effect (-56.5 mL; -108.3 to -4.7; p = 0.033), suggesting that the detrimental effect in adulthood may be preceded by a short-term post-pubertal benefit. Our secondary analyses showing results in the same direction in men and boys, in whom age at menarche SNPs have also shown association with sexual development, suggest a role for pubertal timing in general rather than menarche specifically. We found no effect on airway obstruction (FEV1/FVC).

Genomic analyses identify hundreds of variants associated with age at menarche and support a role for puberty timing in cancer risk

The timing of puberty is a highly polygenic childhood trait that is epidemiologically associated with various adult diseases. Using 1000 Genomes Project-imputed genotype data in up to ∼370,000 women, we identify 389 independent signals (P < 5 × 10-8) for age at menarche, a milestone in female pubertal development. In Icelandic data, these signals explain ∼7.4% of the population variance in age at menarche, corresponding to ∼25% of the estimated heritability. We implicate ∼250 genes via coding variation or associated expression, demonstrating significant enrichment in neural tissues. Rare variants near the imprinted genes MKRN3 and DLK1 were identified, exhibiting large effects when paternally inherited. Mendelian randomization analyses suggest causal inverse associations, independent of body mass index (BMI), between puberty timing and risks for breast and endometrial cancers in women and prostate cancer in men. In aggregate, our findings highlight the complexity of the genetic regulation of puberty timing and support causal links with cancer susceptibility.

Whole-Genome Sequencing Coupled to Imputation Discovers Genetic Signals for Anthropometric Traits

Deep sequence-based imputation can enhance the discovery power of genome-wide association studies by assessing previously unexplored variation across the common- and low-frequency spectra. We applied a hybrid whole-genome sequencing (WGS) and deep imputation approach to examine the broader allelic architecture of 12 anthropometric traits associated with height, body mass, and fat distribution in up to 267,616 individuals. We report 106 genome-wide significant signals that have not been previously identified, including 9 low-frequency variants pointing to functional candidates. Of the 106 signals, 6 are in genomic regions that have not been implicated with related traits before, 28 are independent signals at previously reported regions, and 72 represent previously reported signals for a different anthropometric trait. 71% of signals reside within genes and fine mapping resolves 23 signals to one or two likely causal variants. We confirm genetic overlap between human monogenic and polygenic anthropometric traits and find signal enrichment in cis expression QTLs in relevant tissues. Our results highlight the potential of WGS strategies to enhance biologically relevant discoveries across the frequency spectrum.

Role of DNA Methylation in Type 2 Diabetes Etiology: Using Genotype as a Causal Anchor

Several studies have investigated the relationship between genetic variation and DNA methylation with respect to type 2 diabetes, but it is unknown if DNA methylation is a mediator in the disease pathway or if it is altered in response to disease state. This study uses genotypic information as a causal anchor to help decipher the likely role of DNA methylation measured in peripheral blood in the etiology of type 2 diabetes. Illumina HumanMethylation450 BeadChip data were generated on 1,018 young individuals from the Avon Longitudinal Study of Parents and Children (ALSPAC) cohort. In stage 1, 118 unique associations between published type 2 diabetes single nucleotide polymorphisms (SNPs) and genome-wide methylation (methylation quantitative trait loci [mQTLs]) were identified. In stage 2, a further 226 mQTLs were identified between 202 additional independent non-type 2 diabetes SNPs and CpGs identified in stage 1. Where possible, associations were replicated in independent cohorts of similar age. We discovered that around half of known type 2 diabetes SNPs are associated with variation in DNA methylation and postulated that methylation could either be on a causal pathway to future disease or could be a noncausal biomarker. For one locus (KCNQ1), we were able to provide further evidence that methylation is likely to be on the causal pathway to disease in later life.

Shared genetic influences between dimensional ASD and ADHD symptoms during child and adolescent development

BACKGROUND

Shared genetic influences between attention-deficit/hyperactivity disorder (ADHD) symptoms and autism spectrum disorder (ASD) symptoms have been reported. Cross-trait genetic relationships are, however, subject to dynamic changes during development. We investigated the continuity of genetic overlap between ASD and ADHD symptoms in a general population sample during childhood and adolescence. We also studied uni- and cross-dimensional trait-disorder links with respect to genetic ADHD and ASD risk.

METHODS

Social-communication difficulties (N ≤ 5551, Social and Communication Disorders Checklist, SCDC) and combined hyperactive-impulsive/inattentive ADHD symptoms (N ≤ 5678, Strengths and Difficulties Questionnaire, SDQ-ADHD) were repeatedly measured in a UK birth cohort (ALSPAC, age 7 to 17 years). Genome-wide summary statistics on clinical ASD (5305 cases; 5305 pseudo-controls) and ADHD (4163 cases; 12,040 controls/pseudo-controls) were available from the Psychiatric Genomics Consortium. Genetic trait variances and genetic overlap between phenotypes were estimated using genome-wide data.

RESULTS

In the general population, genetic influences for SCDC and SDQ-ADHD scores were shared throughout development. Genetic correlations across traits reached a similar strength and magnitude (cross-trait rg ≤ 1, pmin  = 3 × 10-4) as those between repeated measures of the same trait (within-trait rg ≤ 0.94, pmin  = 7 × 10-4). Shared genetic influences between traits, especially during later adolescence, may implicate variants in K-RAS signalling upregulated genes (p-meta = 6.4 × 10-4). Uni-dimensionally, each population-based trait mapped to the expected behavioural continuum: risk-increasing alleles for clinical ADHD were persistently associated with SDQ-ADHD scores throughout development (marginal regression R2 = 0.084%). An age-specific genetic overlap between clinical ASD and social-communication difficulties during childhood was also shown, as per previous reports. Cross-dimensionally, however, neither SCDC nor SDQ-ADHD scores were linked to genetic risk for disorder.

CONCLUSIONS

In the general population, genetic aetiologies between social-communication difficulties and ADHD symptoms are shared throughout child and adolescent development and may implicate similar biological pathways that co-vary during development. Within both the ASD and the ADHD dimension, population-based traits are also linked to clinical disorder, although much larger clinical discovery samples are required to reliably detect cross-dimensional trait-disorder relationships.

TRPA1 gene polymorphisms and childhood asthma

BACKGROUND

Animal data have suggested that the transient receptor potential ankyrin-1 (TRPA1) ion channel plays a key role in promoting airway inflammation in asthma and may mediate effects of paracetamol on asthma, yet confirmatory human data are lacking. To study associations of TRPA1 gene variants with childhood asthma and total IgE concentration, and interactions between TRPA1 and prenatal paracetamol exposure on these outcomes.

METHODS

We analysed associations between 31 TRPA1 single nucleotide polymorphisms (SNPs) and current doctor-diagnosed asthma and total IgE concentration at 7.5 years in the Avon Longitudinal Study of Parents and Children (ALSPAC) birth cohort. We sought to confirm the most significant associations with comparable outcomes in the Prevention and Incidence of Asthma and Mite Allergy (PIAMA) and Generation R birth cohorts. In ALSPAC, we explored interactions with prenatal paracetamol exposure.

RESULTS

In ALSPAC, there was strong evidence for association between six SNPs and asthma: rs959974 and rs1384001 (per-allele odds ratio for both: 1.30 (95% CI: 1.15-1.47), p = 0.00001), rs7010969 (OR 1.28 (1.13-1.46), p = 0.00004), rs3735945 (OR 1.30 (1.09-1.55), p = 0.003), rs920829 (OR 1.30 (1.09-1.54), p = 0.004) and rs4738202 (OR 1.22 (1.07-1.39), p = 0.004). In a meta-analysis across the three cohorts, the pooled effect estimates confirmed that all six SNPs were significantly associated with asthma. In ALSPAC, TRPA1 associations with asthma were not modified by prenatal paracetamol, although associations with IgE concentration were.

CONCLUSION

This study suggests that TRPA1 may play a role in the development of childhood asthma. (249 words).

Assessment of rates of recanting and hair testing as a biological measure of drug use in a general population sample of young people

AIMS

We investigate the extent of and factors associated with denial of previously reported cannabis and other illicit drug use, and assess the potential of hair testing for measuring substance use in general population samples.

DESIGN

Birth cohort study.

SETTING

United Kingdom, 1991-present.

PARTICIPANTS

A total of 3643 participants who provided hair and self-report measures of cannabis and other illicit drug use in the Avon Longitudinal Study of Parents and Children (ALSPAC) at age 18 years.

MEASUREMENTS

Denial of ever use of cannabis and other illicit drugs at age 18 following previously reported use. Positive hair drug tests for cannabis and other illicit drugs, and expected numbers of false positives and false negatives based on expected sensitivity and specificity.

FINDINGS

Cannabis and other illicit drug use was reported by 1223 and 393 individuals, respectively, before age 18 years. Of these 176 (14.4%) and 99 (25.2%), respectively, denied use at age 18. Denial of cannabis use decreased with the reporting of other substances and antisocial behaviour. Cannabis and other illicit drug use at age 18 was reported by 547 (22.5%) and 203 (8.4%) individuals, respectively. Of these, 111 (20.3%) and 13 (6.4%) were hair-positive for cannabis and other illicit drugs, respectively. Based on hair testing for cannabis use we expect 0 [95% confidence interval (CI) = 0-169] false positives and 394 (95% CI = 323-449) false negatives compared to observed 362 potential false positives and 436 potential false negatives based on self-report. In hair-positive individuals, reporting the use of other substances and antisocial behaviour decreased the odds of a negative self-report.

CONCLUSIONS

Hair analysis provides an unreliable marker of substance use in general population samples. People who report more frequent substance use before age 18 are less likely to later deny previous substance use at age 18 than people who report occasional use.

Programming of Adiposity in Childhood and Adolescence: Associations With Birth Weight and Cord Blood Adipokines

CONTEXT

Exposure to maternal adiposity during pregnancy is associated with higher offspring birth weight and greater adiposity through childhood and adult life. As birth weight reflects the summation of lean and fat mass, the extent to which fat mass at birth tracks into later life is unknown.

OBJECTIVE

To determine whether fat mass at birth is associated with child and adolescent adiposity.

DESIGN, SETTING, AND PARTICIPANTS

UK birth cohort with markers of neonatal fat mass; cord blood leptin, adiponectin, and birth weight and adiposity outcomes at age 9 (n = 2775) and 17 years (n = 2138).

MAIN OUTCOMES

Offspring body mass index (BMI), waist circumference, dual-energy X-ray absorptiometry-determined fat mass, and obesity at age 9 and 17 years.

RESULTS

Higher cord blood leptin was associated with higher z scores of fat mass [difference in mean per 10 pg/mL: 0.03 standard deviation (SD); 95% confidence interval (CI), 0.00 to 0.06], waist circumference (0.04 SD; 95% CI, 0.00 to 0.07), and BMI (0.04 SD; 95% CI, 0.00 to 0.08) at age 9. However, by age 17 the adjusted results were attenuated to the null. Cord blood adiponectin was not associated with measures of adiposity at age 9. At age 17, cord blood adiponectin was positively associated with fat mass (0.02 SD per 10 μg/mL; 95% CI, 0.02 to 0.03) and waist circumference (0.04 SD per 10 μg/mL; 95% CI, 0.03 to 0.05). Birth weight was positively associated with waist circumference (0.03 SD per 100 g; 95% CI, 0.02 to 0.04) and BMI (0.02 SD per 100 g; 95% CI, 0.00 to 0.03), but not fat mass or odds of obesity. Cord blood leptin and adiponectin were not associated with obesity at either age.

CONCLUSIONS

Increased cord blood leptin and adiponectin, known surrogates of fetal fat mass, were weakly associated with increased fat mass in late childhood and adolescence, respectively.

Genome-wide associations for birth weight and correlations with adult disease

Birth weight (BW) has been shown to be influenced by both fetal and maternal factors and in observational studies is reproducibly associated with future risk of adult metabolic diseases including type 2 diabetes (T2D) and cardiovascular disease. These life-course associations have often been attributed to the impact of an adverse early life environment. Here, we performed a multi-ancestry genome-wide association study (GWAS) meta-analysis of BW in 153,781 individuals, identifying 60 loci where fetal genotype was associated with BW (P < 5 × 10-8). Overall, approximately 15% of variance in BW was captured by assays of fetal genetic variation. Using genetic association alone, we found strong inverse genetic correlations between BW and systolic blood pressure (Rg = -0.22, P = 5.5 × 10-13), T2D (Rg = -0.27, P = 1.1 × 10-6) and coronary artery disease (Rg = -0.30, P = 6.5 × 10-9). In addition, using large -cohort datasets, we demonstrated that genetic factors were the major contributor to the negative covariance between BW and future cardiometabolic risk. Pathway analyses indicated that the protein products of genes within BW-associated regions were enriched for diverse processes including insulin signalling, glucose homeostasis, glycogen biosynthesis and chromatin remodelling. There was also enrichment of associations with BW in known imprinted regions (P = 1.9 × 10-4). We demonstrate that life-course associations between early growth phenotypes and adult cardiometabolic disease are in part the result of shared genetic effects and identify some of the pathways through which these causal genetic effects are mediated.

A Genome-Wide Association Meta-Analysis of Attention-Deficit/Hyperactivity Disorder Symptoms in Population-Based Pediatric Cohorts

OBJECTIVE

The aims of this study were to elucidate the influence of common genetic variants on childhood attention-deficit/hyperactivity disorder (ADHD) symptoms, to identify genetic variants that explain its high heritability, and to investigate the genetic overlap of ADHD symptom scores with ADHD diagnosis.

METHOD

Within the EArly Genetics and Lifecourse Epidemiology (EAGLE) consortium, genome-wide single nucleotide polymorphisms (SNPs) and ADHD symptom scores were available for 17,666 children (<13 years of age) from nine population-based cohorts. SNP-based heritability was estimated in data from the three largest cohorts. Meta-analysis based on genome-wide association (GWA) analyses with SNPs was followed by gene-based association tests, and the overlap in results with a meta-analysis in the Psychiatric Genomics Consortium (PGC) case-control ADHD study was investigated.

RESULTS

SNP-based heritability ranged from 5% to 34%, indicating that variation in common genetic variants influences ADHD symptom scores. The meta-analysis did not detect genome-wide significant SNPs, but three genes, lying close to each other with SNPs in high linkage disequilibrium (LD), showed a gene-wide significant association (p values between 1.46 × 10(-6) and 2.66 × 10(-6)). One gene, WASL, is involved in neuronal development. Both SNP- and gene-based analyses indicated overlap with the PGC meta-analysis results with the genetic correlation estimated at 0.96.

CONCLUSION

The SNP-based heritability for ADHD symptom scores indicates a polygenic architecture, and genes involved in neurite outgrowth are possibly involved. Continuous and dichotomous measures of ADHD appear to assess a genetically common phenotype. A next step is to combine data from population-based and case-control cohorts in genetic association studies to increase sample size and to improve statistical power for identifying genetic variants.

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.

A genome-wide approach to children's aggressive behavior: The EAGLE consortium

Individual differences in aggressive behavior emerge in early childhood and predict persisting behavioral problems and disorders. Studies of antisocial and severe aggression in adulthood indicate substantial underlying biology. However, little attention has been given to genome-wide approaches of aggressive behavior in children. We analyzed data from nine population-based studies and assessed aggressive behavior using well-validated parent-reported questionnaires. This is the largest sample exploring children's aggressive behavior to date (N = 18,988), with measures in two developmental stages (N = 15,668 early childhood and N = 16,311 middle childhood/early adolescence). First, we estimated the additive genetic variance of children's aggressive behavior based on genome-wide SNP information, using genome-wide complex trait analysis (GCTA). Second, genetic associations within each study were assessed using a quasi-Poisson regression approach, capturing the highly right-skewed distribution of aggressive behavior. Third, we performed meta-analyses of genome-wide associations for both the total age-mixed sample and the two developmental stages. Finally, we performed a gene-based test using the summary statistics of the total sample. GCTA quantified variance tagged by common SNPs (10-54%). The meta-analysis of the total sample identified one region in chromosome 2 (2p12) at near genome-wide significance (top SNP rs11126630, P = 5.30 × 10(-8) ). The separate meta-analyses of the two developmental stages revealed suggestive evidence of association at the same locus. The gene-based analysis indicated association of variation within AVPR1A with aggressive behavior. We conclude that common variants at 2p12 show suggestive evidence for association with childhood aggression. Replication of these initial findings is needed, and further studies should clarify its biological meaning. © 2015 Wiley Periodicals, Inc.

Meta-analysis of 375,000 individuals identifies 38 susceptibility loci for migraine

Migraine is a debilitating neurological disorder affecting around one in seven people worldwide, but its molecular mechanisms remain poorly understood. There is some debate about whether migraine is a disease of vascular dysfunction or a result of neuronal dysfunction with secondary vascular changes. Genome-wide association (GWA) studies have thus far identified 13 independent loci associated with migraine. To identify new susceptibility loci, we carried out a genetic study of migraine on 59,674 affected subjects and 316,078 controls from 22 GWA studies. We identified 44 independent single-nucleotide polymorphisms (SNPs) significantly associated with migraine risk (P < 5 × 10(-8)) that mapped to 38 distinct genomic loci, including 28 loci not previously reported and a locus that to our knowledge is the first to be identified on chromosome X. In subsequent computational analyses, the identified loci showed enrichment for genes expressed in vascular and smooth muscle tissues, consistent with a predominant theory of migraine that highlights vascular etiologies.

Genome-wide association study identifies 74 loci associated with educational attainment

Educational attainment is strongly influenced by social and other environmental factors, but genetic factors are estimated to account for at least 20% of the variation across individuals. Here we report the results of a genome-wide association study (GWAS) for educational attainment that extends our earlier discovery sample of 101,069 individuals to 293,723 individuals, and a replication study in an independent sample of 111,349 individuals from the UK Biobank. We identify 74 genome-wide significant loci associated with the number of years of schooling completed. Single-nucleotide polymorphisms associated with educational attainment are disproportionately found in genomic regions regulating gene expression in the fetal brain. Candidate genes are preferentially expressed in neural tissue, especially during the prenatal period, and enriched for biological pathways involved in neural development. Our findings demonstrate that, even for a behavioural phenotype that is mostly environmentally determined, a well-powered GWAS identifies replicable associated genetic variants that suggest biologically relevant pathways. Because educational attainment is measured in large numbers of individuals, it will continue to be useful as a proxy phenotype in efforts to characterize the genetic influences of related phenotypes, including cognition and neuropsychiatric diseases.

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.

Electronic supplementary material

The online version of this article (doi:10.1186/s13059-016-0926-z) contains supplementary material, which is available to authorized users.

Genetic Evidence for Causal Relationships Between Maternal Obesity-Related Traits and Birth Weight

IMPORTANCE

Neonates born to overweight or obese women are larger and at higher risk of birth complications. Many maternal obesity-related traits are observationally associated with birth weight, but the causal nature of these associations is uncertain.

OBJECTIVE

To test for genetic evidence of causal associations of maternal body mass index (BMI) and related traits with birth weight.

DESIGN, SETTING, AND PARTICIPANTS

Mendelian randomization to test whether maternal BMI and obesity-related traits are potentially causally related to offspring birth weight. Data from 30,487 women in 18 studies were analyzed. Participants were of European ancestry from population- or community-based studies in Europe, North America, or Australia and were part of the Early Growth Genetics Consortium. Live, term, singleton offspring born between 1929 and 2013 were included.

EXPOSURES

Genetic scores for BMI, fasting glucose level, type 2 diabetes, systolic blood pressure (SBP), triglyceride level, high-density lipoprotein cholesterol (HDL-C) level, vitamin D status, and adiponectin level.

MAIN OUTCOME AND MEASURE

Offspring birth weight from 18 studies.

RESULTS

Among the 30,487 newborns the mean birth weight in the various cohorts ranged from 3325 g to 3679 g. The maternal genetic score for BMI was associated with a 2-g (95% CI, 0 to 3 g) higher offspring birth weight per maternal BMI-raising allele (P = .008). The maternal genetic scores for fasting glucose and SBP were also associated with birth weight with effect sizes of 8 g (95% CI, 6 to 10 g) per glucose-raising allele (P = 7 × 10(-14)) and -4 g (95% CI, -6 to -2 g) per SBP-raising allele (P = 1×10(-5)), respectively. A 1-SD ( ≈ 4 points) genetically higher maternal BMI was associated with a 55-g higher offspring birth weight (95% CI, 17 to 93 g). A 1-SD ( ≈ 7.2 mg/dL) genetically higher maternal fasting glucose concentration was associated with 114-g higher offspring birth weight (95% CI, 80 to 147 g). However, a 1-SD ( ≈ 10 mm Hg) genetically higher maternal SBP was associated with a 208-g lower offspring birth weight (95% CI, -394 to -21 g). For BMI and fasting glucose, genetic associations were consistent with the observational associations, but for systolic blood pressure, the genetic and observational associations were in opposite directions.

CONCLUSIONS AND RELEVANCE

In this mendelian randomization study, genetically elevated maternal BMI and blood glucose levels were potentially causally associated with higher offspring birth weight, whereas genetically elevated maternal SBP was potentially causally related to lower birth weight. If replicated, these findings may have implications for counseling and managing pregnancies to avoid adverse weight-related birth outcomes.

Genome-wide association analysis identifies three new susceptibility loci for childhood body mass index

A large number of genetic loci are associated with adult body mass index. However, the genetics of childhood body mass index are largely unknown. We performed a meta-analysis of genome-wide association studies of childhood body mass index, using sex- and age-adjusted standard deviation scores. We included 35 668 children from 20 studies in the discovery phase and 11 873 children from 13 studies in the replication phase. In total, 15 loci reached genome-wide significance (P-value < 5 × 10(-8)) in the joint discovery and replication analysis, of which 12 are previously identified loci in or close to ADCY3, GNPDA2, TMEM18, SEC16B, FAIM2, FTO, TFAP2B, TNNI3K, MC4R, GPR61, LMX1B and OLFM4 associated with adult body mass index or childhood obesity. We identified three novel loci: rs13253111 near ELP3, rs8092503 near RAB27B and rs13387838 near ADAM23. Per additional risk allele, body mass index increased 0.04 Standard Deviation Score (SDS) [Standard Error (SE) 0.007], 0.05 SDS (SE 0.008) and 0.14 SDS (SE 0.025), for rs13253111, rs8092503 and rs13387838, respectively. A genetic risk score combining all 15 SNPs showed that each additional average risk allele was associated with a 0.073 SDS (SE 0.011, P-value = 3.12 × 10(-10)) increase in childhood body mass index in a population of 1955 children. This risk score explained 2% of the variance in childhood body mass index. This study highlights the shared genetic background between childhood and adult body mass index and adds three novel loci. These loci likely represent age-related differences in strength of the associations with body mass index.

Genome-wide association study identifies common and low-frequency variants at the AMH gene locus that strongly predict serum AMH levels in males

Anti-Müllerian hormone (AMH) is an essential messenger of sexual differentiation in the foetus and is an emerging biomarker of postnatal reproductive function in females. Due to a paucity of adequately sized studies, the genetic determinants of circulating AMH levels are poorly characterized. In samples from 2815 adolescents aged 15 from the ALSPAC study, we performed the first genome-wide association study of serum AMH levels across a set of ∼9 m '1000 Genomes Reference Panel' imputed genetic variants. Genetic variants at the AMH protein-coding gene showed considerable allelic heterogeneity, with both common variants [rs4807216 (P(Male) = 2 × 10(-49), Beta: ∼0.9 SDs per allele), rs8112524 (P(Male) = 3 × 10(-8), Beta: ∼0.25)] and low-frequency variants [rs2385821 (P(Male) = 6 × 10(-31), Beta: ∼1.2, frequency 3.6%)] independently associated with apparently large effect sizes in males, but not females. For all three SNPs, we highlight mechanistic links to AMH gene function and demonstrate highly significant sex interactions (P(Het) 0.0003-6.3 × 10(-12)), culminating in contrasting estimates of trait variance explained (24.5% in males versus 0.8% in females). Using these SNPs as a genetic proxy for AMH levels, we found no evidence in additional datasets to support a biological role for AMH in complex traits and diseases in men.

The regulatory element READ1 epistatically influences reading and language, with both deleterious and protective alleles

Background

Reading disability (RD) and language impairment (LI) are heritable learning disabilities that obstruct acquisition and use of written and spoken language, respectively. We previously reported that two risk haplotypes, each in strong linkage disequilibrium (LD) with an allele of READ1, a polymorphic compound short tandem repeat within intron 2 of risk gene DCDC2, are associated with RD and LI. Additionally, we showed a non-additive genetic interaction between READ1 and KIAHap, a previously reported risk haplotype in risk gene KIAA0319, and that READ1 binds the transcriptional regulator ETV6.

Objective

To examine the hypothesis that READ1 is a transcriptional regulator of KIAA0319.

Methods

We characterised associations between READ1 alleles and RD and LI in a large European cohort, and also assessed interactions between READ1 and KIAHap and their effect on performance on measures of reading, language and IQ. We also used family-based data to characterise the genetic interaction, and chromatin conformation capture (3C) to investigate the possibility of a physical interaction between READ1 and KIAHap.

Results and conclusions

READ1 and KIAHap show interdependence—READ1 risk alleles synergise with KIAHap, whereas READ1 protective alleles act epistatically to negate the effects of KIAHap. The family data suggest that these variants interact in trans genetically, while the 3C results show that a region of DCDC2 containing READ1 interacts physically with the region upstream of KIAA0319. These data support a model in which READ1 regulates KIAA0319 expression through KIAHap and in which the additive effects of READ1 and KIAHap alleles are responsible for the trans genetic interaction.

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.

Lymphoblastoid cell lines reveal associations of adult DNA methylation with childhood and current adversity that are distinct from whole blood associations

BACKGROUND

Some cohort studies bank lymphoblastoid cell lines (LCLs) as a renewable source of participant DNA. However, although LCL DNA has proved valuable for genetic studies, its utility in epigenetic epidemiology research is unknown.

METHODS

To assess whether LCL DNA can be used for life-course environmental epigenomics, we carried out a pilot methylomic study (using the Illumina Infinium Human Methylation 450 BeadChip) of nil-passage, Epstein-Barr virus (EBV)-transformed LCLs (n = 42) and 28 matched whole-blood (WB) samples. These were from adult male participants of the British 1958 birth cohort, selected for extremes of social economic position (SEP) in childhood and adulthood, with additional information available on childhood abuse and prenatal tobacco exposure.

RESULTS

We identified a small number of weak associations between these exposures and methylation levels of both individual CpG sites and genomic regions in WB and LCLs. However, only one of the regional, and none of the individual CpG site associations were common to both sample types. The lack of overlap between the associations detected in LCL compared with those found in WB could either be due to the EBV-transformation process, or to the fact that, unlike WB, LCLs are essentially a single (CD19+) cell type. We provide evidence that the latter is the more potent explanation, by showing that CpG sites known to be differentially methylated between different types of blood cell have significantly lower correlations (R = 0.11) than average (R = 0.2) between WB and LCLs in our datasets, whereas sites known to be affected by EBV-transformation have significantly higher correlations (R = 0.3).

CONCLUSIONS

This small pilot study suggests that the DNA methylation profile of LCLs is more closely related to that of B cells than WB and, additionally, that LCLs may nevertheless be useful for life-course environmental epigenomics.

Genome-wide association study identifies SNPs in the MHC class II loci that are associated with self-reported history of whooping cough

Whooping cough is currently seeing resurgence in countries despite high vaccine coverage. There is considerable variation in subject-specific response to infection and vaccine efficacy, but little is known about the role of human genetics. We carried out a case–control genome-wide association study of adult or parent-reported history of whooping cough in two cohorts from the UK: the ALSPAC cohort and the 1958 British Birth Cohort (815/758 cases and 6341/4308 controls, respectively). We also imputed HLA alleles using dense SNP data in the MHC region and carried out gene-based and gene-set tests of association and estimated the amount of additive genetic variation explained by common SNPs. We observed a novel association at SNPs in the MHC class II region in both cohorts [lead SNP rs9271768 after meta-analysis, odds ratio [95% confidence intervals (CIs)] 1.47 (1.35, 1.6), P-value 1.21E − 18]. Multiple strong associations were also observed at alleles at the HLA class II loci. The majority of these associations were explained by the lead SNP rs9271768. Gene-based and gene-set tests and estimates of explainable common genetic variation could not establish the presence of additional associations in our sample. Genetic variation at the MHC class II region plays a role in susceptibility to whooping cough. These findings provide additional perspective on mechanisms of whooping cough infection and vaccine efficacy.