Epigenetic aging & embodying injustice: US My Body My Story and Multi-Ethnic Atherosclerosis Study

Importance

Epigenetic accelerated aging is associated with exposure to social and economic adversity and may increase risk of premature morbidity and mortality. However, no studies have included measures of structural racism and few have compared estimates within or across the 1st and 2nd generation of epigenetic clocks (the latter additionally trained on phenotypic data).

Objective

To determine if accelerated epigenetic aging is associated with exposures to diverse measures of racialized, economic, and environmental injustice measured at different levels and time periods.

Design

Cross-sectional My Body My Story Study (MBMS; US, 2008-2010) and Exam 5 Multi-Ethnic Atherosclerosis Study (MESA; US, 2010-2012). MBMS DNA extraction: 2021; linkage of structural measures to MBMS and MESA: 2022.

Setting

MBMS recruited a random sample of US-born Black non-Hispanic (BNH) and white non-Hispanic (WNH) participants from 4 community health centers in Boston, MA. The MESA Exam 5 epigenetic component included 975 randomly selected US-born BNH, WNH, and Hispanic participants from four field sites: Baltimore, MD; Forsyth County, NC; New York City, NY; St. Paul, MN.

Participants

US-born persons (MBMS: 224 BNH, 69 WNH; MESA: 229 BNH, 555 WNH, 191 Hispanic).

Main outcome and measures

10 epigenetic clocks (six 1st generation; four 2nd generation), computed using DNA methylation data (DNAm) from blood spots (MBMS; N = 293) and purified monocytes (MESA; N = 975).

Results

Among Black non-Hispanic MBMS participants, epigenetic age acceleration was associated with being born in a Jim Crow state by 0.14 standard deviations (95% confidence interval [CI] 0.00, 0.27) and with birth state conservatism (0.06, 95% CI 0.00, 0.05), pooling across all clocks, as was low parental education for both Black non-Hispanic and white non-Hispanic MBMS participants (respectively: 0.24, 95% CI 0.08, 0.39, and 0.27, 95% CI 0.03, 0.51. Adult impoverishment was positively associated with the pooled 2nd generation clocks among the MESA participants (Black non-Hispanic: 0.06, 95% CI 0.01, 0.12; white non-Hispanic: 0.05, 95% CI 0.01, 0.08; Hispanic: 0.07, 95% CI 0.01, 0.14).

Conclusions and Relevance

Epigenetic accelerated aging may be one of the biological mechanisms linking exposure to racialized and economic injustice to well-documented inequities in premature morbidity and mortality.

An epigenome-wide analysis of DNA methylation, racialized and economic inequities, and air pollution

Importance

DNA methylation (DNAm) provides a plausible mechanism by which adverse exposures become embodied and contribute to health inequities, due to its role in genome regulation and responsiveness to social and biophysical exposures tied to societal context. However, scant epigenome-wide association studies (EWAS) have included structural and lifecourse measures of exposure, especially in relation to structural discrimination.

Objective

Our study tests the hypothesis that DNAm is a mechanism by which racial discrimination, economic adversity, and air pollution become biologically embodied.

Design

A series of cross-sectional EWAS, conducted in My Body My Story (MBMS, biological specimens collected 2008-2010, DNAm assayed in 2021); and the Multi Ethnic Study of Atherosclerosis (MESA; biological specimens collected 2010-2012, DNAm assayed in 2012-2013); using new georeferenced social exposure data for both studies (generated in 2022).

Setting

MBMS was recruited from four community health centers in Boston; MESA was recruited from four field sites in: Baltimore, MD; Forsyth County, NC; New York City, NY; and St. Paul, MN.

Participants

Two population-based samples of US-born Black non-Hispanic (Black NH), white non-Hispanic (white NH), and Hispanic individuals (MBMS; n=224 Black NH and 69 white NH) and (MESA; n=229 Black NH, n=555 white NH and n=191 Hispanic).

Exposures

Eight social exposures encompassing racial discrimination, economic adversity, and air pollution.

Main outcome

Genome-wide changes in DNAm, as measured using the Illumina EPIC BeadChip (MBMS; using frozen blood spots) and Illumina 450k BeadChip (MESA; using purified monocytes). Our hypothesis was formulated after data collection.

Results

We observed the strongest associations with traffic-related air pollution (measured via black carbon and nitrogen oxides exposure), with evidence from both studies suggesting that air pollution exposure may induce epigenetic changes related to inflammatory processes. We also found suggestive associations of DNAm variation with measures of structural racial discrimination (e.g., for Black NH participants, born in a Jim Crow state; adult exposure to racialized economic residential segregation) situated in genes with plausible links to effects on health.

Conclusions and Relevance

Overall, this work suggests that DNAm is a biological mechanism through which structural racism and air pollution become embodied and may lead to health inequities.

Design and quality control of large-scale two-sample Mendelian randomization studies

Background

Mendelian randomization (MR) studies are susceptible to metadata errors (e.g. incorrect specification of the effect allele column) and other analytical issues that can introduce substantial bias into analyses. We developed a quality control (QC) pipeline for the Fatty Acids in Cancer Mendelian Randomization Collaboration (FAMRC) that can be used to identify and correct for such errors.

Methods

We collated summary association statistics from fatty acid and cancer genome-wide association studies (GWAS) and subjected the collated data to a comprehensive QC pipeline. We identified metadata errors through comparison of study-specific statistics to external reference data sets (the National Human Genome Research Institute-European Bioinformatics Institute GWAS catalogue and 1000 genome super populations) and other analytical issues through comparison of reported to expected genetic effect sizes. Comparisons were based on three sets of genetic variants: (i) GWAS hits for fatty acids, (ii) GWAS hits for cancer and (iii) a 1000 genomes reference set.

Results

We collated summary data from 6 fatty acid and 54 cancer GWAS. Metadata errors and analytical issues with the potential to introduce substantial bias were identified in seven studies (11.6%). After resolving metadata errors and analytical issues, we created a data set of 219 842 genetic associations with 90 cancer types, generated in analyses of 566 665 cancer cases and 1 622 374 controls.

Conclusions

In this large MR collaboration, 11.6% of included studies were affected by a substantial metadata error or analytical issue. By increasing the integrity of collated summary data prior to their analysis, our protocol can be used to increase the reliability of downstream MR analyses. Our pipeline is available to other researchers via the CheckSumStats package (https://github.com/MRCIEU/CheckSumStats).

Genetic impacts on DNA methylation help elucidate regulatory genomic processes

BACKGROUND

Pinpointing genetic impacts on DNA methylation can improve our understanding of pathways that underlie gene regulation and disease risk.

RESULTS

We report heritability and methylation quantitative trait locus (meQTL) analysis at 724,499 CpGs profiled with the Illumina Infinium MethylationEPIC array in 2358 blood samples from three UK cohorts. Methylation levels at 34.2% of CpGs are affected by SNPs, and 98% of effects are cis-acting or within 1 Mbp of the tested CpG. Our results are consistent with meQTL analyses based on the former Illumina Infinium HumanMethylation450 array. Both SNPs and CpGs with meQTLs are overrepresented in enhancers, which have improved coverage on this platform compared to previous approaches. Co-localisation analyses across genetic effects on DNA methylation and 56 human traits identify 1520 co-localisations across 1325 unique CpGs and 34 phenotypes, including in disease-relevant genes, such as USP1 and DOCK7 (total cholesterol levels), and ICOSLG (inflammatory bowel disease). Enrichment analysis of meQTLs and integration with expression QTLs give insights into mechanisms underlying cis-meQTLs (e.g. through disruption of transcription factor binding sites for CTCF and SMC3) and trans-meQTLs (e.g. through regulating the expression of ACD and SENP7 which can modulate DNA methylation at distal sites).

CONCLUSIONS

Our findings improve the characterisation of the mechanisms underlying DNA methylation variability and are informative for prioritisation of GWAS variants for functional follow-ups. The MeQTL EPIC Database and viewer are available online at https://epicmeqtl.kcl.ac.uk .

Epigenetic clocks and research implications of the lack of data on whom they have been developed: a review of reported and missing sociodemographic characteristics

Epigenetic clocks are increasingly being used as a tool to assess the impact of a wide variety of phenotypes and exposures on healthy ageing, with a recent focus on social determinants of health. However, little attention has been paid to the sociodemographic characteristics of participants on whom these clocks have been based. Participant characteristics are important because sociodemographic and socioeconomic factors are known to be associated with both DNA methylation variation and healthy ageing. It is also well known that machine learning algorithms have the potential to exacerbate health inequities through the use of unrepresentative samples - prediction models may underperform in social groups that were poorly represented in the training data used to construct the model. To address this gap in the literature, we conducted a review of the sociodemographic characteristics of the participants whose data were used to construct 13 commonly used epigenetic clocks. We found that although some of the epigenetic clocks were created utilizing data provided by individuals from different ages, sexes/genders, and racialized groups, sociodemographic characteristics are generally poorly reported. Reported information is limited by inadequate conceptualization of the social dimensions and exposure implications of gender and racialized inequality, and socioeconomic data are infrequently reported. It is important for future work to ensure clear reporting of tangible data on the sociodemographic and socioeconomic characteristics of all the participants in the study to ensure that other researchers can make informed judgements about the appropriateness of the model for their study population.

Use of Correct and Incorrect Methods of Accounting for Age in Studies of Epigenetic Accelerated Aging: Implications and Recommendations for Best Practices

Motivated by our conduct of a literature review on social exposures and accelerated aging as measured by a growing number of epigenetic "clocks" (which estimate age via DNA methylation (DNAm) patterns), we report on 3 different approaches in the epidemiologic literature-1 incorrect and 2 correct-on the treatment of age in these and other studies using other common exposures (i.e., body mass index and alcohol consumption). Among the 50 empirical articles reviewed, the majority (n = 29; 58%) used the incorrect method of analyzing accelerated aging detrended for age as the outcome and did not control for age as a covariate. By contrast, only 42% used correct methods, which are either to analyze accelerated aging detrended for age as the outcome and control for age as a covariate (n = 16; 32%) or to analyze raw DNAm age as the outcome and control for age as a covariate (n = 5; 10%). In accord with prior demonstrations of bias introduced by use of the incorrect approach, we provide simulation analyses and additional empirical analyses to illustrate how the incorrect method can lead to bias towards the null, and we discuss implications for extant research and recommendations for best practices.

The association between genetically elevated polyunsaturated fatty acids and risk of cancer

BACKGROUND

The causal relevance of polyunsaturated fatty acids (PUFAs) for risk of site-specific cancers remains uncertain.

METHODS

Using a Mendelian randomization (MR) framework, we assessed the causal relevance of PUFAs for risk of cancer in European and East Asian ancestry individuals. We defined the primary exposure as PUFA desaturase activity, proxied by rs174546 at the FADS locus. Secondary exposures were defined as omega 3 and omega 6 PUFAs that could be proxied by genetic polymorphisms outside the FADS region. Our study used summary genetic data on 10 PUFAs and 67 cancers, corresponding to 562,871 cases and 1,619,465 controls, collected by the Fatty Acids in Cancer Mendelian Randomization Collaboration. We estimated odds ratios (ORs) for cancer per standard deviation increase in genetically proxied PUFA exposures.

FINDINGS

Genetically elevated PUFA desaturase activity was associated (P < 0.0007) with higher risk (OR [95% confidence interval]) of colorectal cancer (1.09 [1.07-1.11]), esophageal squamous cell carcinoma (1.16 [1.06-1.26]), lung cancer (1.06 [1.03-1.08]) and basal cell carcinoma (1.05 [1.02-1.07]). There was little evidence for associations with reproductive cancers (OR = 1.00 [95% CI: 0.99-1.01]; Pheterogeneity = 0.25), urinary system cancers (1.03 [0.99-1.06], Pheterogeneity = 0.51), nervous system cancers (0.99 [0.95-1.03], Pheterogeneity = 0.92) or blood cancers (1.01 [0.98-1.04], Pheterogeneity = 0.09). Findings for colorectal cancer and esophageal squamous cell carcinoma remained compatible with causality in sensitivity analyses for violations of assumptions. Secondary MR analyses highlighted higher omega 6 PUFAs (arachidonic acid, gamma-linolenic acid and dihomo-gamma-linolenic acid) as potential mediators. PUFA biosynthesis is known to interact with aspirin, which increases risk of bleeding and inflammatory bowel disease. In a phenome-wide MR study of non-neoplastic diseases, we found that genetic lowering of PUFA desaturase activity, mimicking a hypothetical intervention to reduce cancer risk, was associated (P < 0.0006) with increased risk of inflammatory bowel disease but not bleeding.

INTERPRETATION

The PUFA biosynthesis pathway may be an intervention target for prevention of colorectal cancer and esophageal squamous cell carcinoma but with potential for increased risk of inflammatory bowel disease.

FUNDING

Cancer Resesrch UK (C52724/A20138, C18281/A19169). UK Medical Research Council (MR/P014054/1). National Institute for Health Research (NIHR202411). UK Medical Research Council (MC_UU_00011/1, MC_UU_00011/3, MC_UU_00011/6, and MC_UU_00011/4). National Cancer Institute (R00 CA215360). National Institutes of Health (U01 CA164973, R01 CA60987, R01 CA72520, U01 CA74806, R01 CA55874, U01 CA164973 and U01 CA164973).

Role of DNA methylation in the relationship between glioma risk factors and glioma incidence: a two-step Mendelian randomization study

Genetic evidence suggests glioma risk is altered by leukocyte telomere length, allergic disease (asthma, hay fever or eczema), alcohol consumption, childhood obesity, low-density lipoprotein cholesterol (LDLc) and triglyceride levels. DNA methylation (DNAm) variation influences many of these glioma-related traits and is an established feature of glioma. Yet the causal relationship between DNAm variation with both glioma incidence and glioma risk factors is unknown. We applied a two-step Mendelian randomization (MR) approach and several sensitivity analyses (including colocalization and Steiger filtering) to assess the association of DNAm with glioma risk factors and glioma incidence. We used data from a recently published catalogue of germline genetic variants robustly associated with DNAm variation in blood (32,851 participants) and data from a genome-wide association study of glioma risk (12,488 cases and 18,169 controls, sub-divided into 6191 glioblastoma cases and 6305 non-glioblastoma cases). MR evidence indicated that DNAm at 3 CpG sites (cg01561092, cg05926943, cg01584448) in one genomic region (HEATR3) had a putative association with glioma and glioblastoma risk (False discovery rate [FDR] < 0.05). Steiger filtering provided evidence against reverse causation. Colocalization presented evidence against genetic confounding and suggested that differential DNAm at the 3 CpG sites and glioma were driven by the same genetic variant. MR provided little evidence to suggest that DNAm acts as a mediator on the causal pathway between risk factors previously examined and glioma onset. To our knowledge, this is the first study to use MR to appraise the causal link of DNAm with glioma risk factors and glioma onset. Subsequent analyses are required to improve the robustness of our results and rule out horizontal pleiotropy.

Abstract 6504: Germline genetically predicted body mass index is associated with endometrial cancer somatic transcriptomic, immune, and mutational signatures in The Cancer Genome Atlas

High body mass index (BMI) is a causal risk factor for endometrial cancer but the molecular mechanisms underlying this association remain elusive. Here we sought to characterize the tumor genomic landscape of endometrial cancers that have developed on a germline genetic background of predisposition to elevated BMI. We built a polygenic score (PGS) for adult BMI in women using effect size estimates and allele information on 242 independent (r2&lt;0.05) variants associated with BMI at genome-wide significance (P&lt;5x10-9) in 379,501 women of European ancestry. We performed sample and germline (blood) genotype quality control and imputation into the 1000 Genomes reference panel on data from 354 endometrial cancer cases of genetically inferred European ancestry from The Cancer Genome Atlas (TCGA). We assigned each woman in this TCGA cohort her genetically predicted life-course BMI based on the BMI PGS and found this to be modestly correlated with BMI at the point of diagnosis (r2=0.23; P=2.1x10-5). Multivariable linear (default) and quasi-Poisson (for zero-inflated counts) regression models were used to test for associations between the BMI germline PGS and endometrial cancer tumor genomic, transcriptomic, proteomic, and immune traits in TCGA. All analyses were adjusted for age, stage, microsatellite status and 10 genetic principal components. Mutational signature models were also adjusted for signature accuracy. We ranked 18,458 genes based on the association between their tumor expression and the BMI PGS and performed gene set enrichment analysis to identify associations between the BMI PGS and upregulation of genes in the IL6-JAK-STAT3 signaling (false discovery rate (FDR)=8.50x10-7), inflammatory response (FDR=7.03x10-6), interferon gamma response (FDR=5.49x10-5) and glycolysis (FDR=3.28x10-4) pathways. High BMI PGS had an inverse association with endometrial tumor EGFR (FDR=0.07) protein levels of the 131 tumor proteins profiled by reverse phase protein array. Endometrial tumors that had developed on a germline background predictive of high BMI were also associated with increased infiltration of activated mast cells (FDR=9.55x10-3) in our evaluation of 22 tumor immune cell infiltrates quantified by the CIBERSORT algorithm, as well as the mitotic and aging clock-like single base substitution (SBS) signatures 1 (FDR=0.01) and 5 (FDR=0.04). The two SBS signature associations and the activated mast cell association with the BMI PGS were substantially more pronounced in the subgroup of endometrial cancers with microsatellite instability. Thus, we combined germline and somatic data using a novel approach to identify endometrial cancer tumor molecular features associated with genetically predicted higher BMI, providing precision multi-omic portraits of endometrial cancers that develop on a background of adiposity.

Citation Format: George Richenberg, Victoria Gray, Carina Owen, Tom Gaunt, Caroline Relton, Emma Vincent, Siddhartha Kar. Germline genetically predicted body mass index is associated with endometrial cancer somatic transcriptomic, immune, and mutational signatures in The Cancer Genome Atlas [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2023; Part 1 (Regular and Invited Abstracts); 2023 Apr 14-19; Orlando, FL. Philadelphia (PA): AACR; Cancer Res 2023;83(7_Suppl):Abstract nr 6504.

Analysis of Systemic Epigenetic Alterations in Inflammatory Bowel Disease: Defining Geographical, Genetic and Immune-Inflammatory influences on the Circulating Methylome

BACKGROUND

Epigenetic alterations may provide valuable insights into gene-environment interactions in the pathogenesis of inflammatory bowel disease [IBD].

METHODS

Genome-wide methylation was measured from peripheral blood using the Illumina 450k platform in a case-control study in an inception cohort (295 controls, 154 Crohn's disease [CD], 161 ulcerative colitis [UC], 28 IBD unclassified [IBD-U)] with covariates of age, sex and cell counts, deconvoluted by the Houseman method. Genotyping was performed using Illumina HumanOmniExpressExome-8 BeadChips and gene expression using the Ion AmpliSeq Human Gene Expression Core Panel. Treatment escalation was characterized by the need for biological agents or surgery after initial disease remission.

RESULTS

A total of 137 differentially methylated positions [DMPs] were identified in IBD, including VMP1/MIR21 [p = 9.11 × 10-15] and RPS6KA2 [6.43 × 10-13], with consistency seen across Scandinavia and the UK. Dysregulated loci demonstrate strong genetic influence, notably VMP1 [p = 1.53 × 10-15]. Age acceleration is seen in IBD [coefficient 0.94, p < 2.2 × 10-16]. Several immuno-active genes demonstrated highly significant correlations between methylation and gene expression in IBD, in particular OSM: IBD r = -0.32, p = 3.64 × 10-7 vs non-IBD r = -0.14, p = 0.77]. Multi-omic integration of the methylome, genome and transcriptome also implicated specific pathways that associate with immune activation, response and regulation at disease inception. At follow-up, a signature of three DMPs [TAP1, TESPA1, RPTOR] were associated with treatment escalation to biological agents or surgery (hazard ratio of 5.19 [CI: 2.14-12.56], logrank p = 9.70 × 10-4).

CONCLUSION

These data demonstrate consistent epigenetic alterations at diagnosis in European patients with IBD, providing insights into the pathogenetic importance and translational potential of epigenetic mapping in complex disease.

Chromosome 18 Loss of Heterozygosity in Small Intestinal Neuroendocrine Tumours: Multi-Omic and Tumour Composition Analyses

INTRODUCTION

Small intestinal neuroendocrine tumours (siNETs) are rare neoplasms which present with low mutational burden and can be subtyped based on copy number variation (CNV). Currently, siNETs can be molecularly classified as having chromosome 18 loss of heterozygosity (18LOH), multiple CNVs (MultiCNV), or no CNVs. 18LOH tumours have better progression-free survival when compared to MultiCNV and NoCNV tumours, however, the mechanism underlying this is unknown, and clinical practice does not currently consider CNV status.

METHODS

Here, we use genome-wide tumour DNA methylation (n = 54) and gene expression (n = 20 matched to DNA methylation) to better understand how gene regulation varies by 18LOH status. We then use multiple cell deconvolution methods to analyse how cell composition varies between 18LOH status and determine potential associations with progression-free survival.

RESULTS

We identified 27,464 differentially methylated CpG sites and 12 differentially expressed genes between 18LOH and non-18LOH (MultiCNV + NoCNV) siNETs. Although few differentially expressed genes were identified, these genes were highly enriched with the differentially methylated CpG sites compared to the rest of the genome. We identified differences in tumour microenvironment between 18LOH and non-18LOH tumours, including CD14+ infiltration in a subset of non-18LOH tumours which had the poorest clinical outcomes.

CONCLUSIONS

We identify a small number of genes which appear to be linked to the 18LOH status of siNETs, and find evidence of potential epigenetic dysregulation of these genes. We also find a potential prognostic marker for worse progression-free outcomes in the form of higher CD14 infiltration in non-18LOH siNETs.

The impact of low input DNA on the reliability of DNA methylation as measured by the Illumina Infinium MethylationEPIC BeadChip

DNA methylation (DNAm) is commonly assayed using the Illumina Infinium MethylationEPIC BeadChip, but there is currently little published evidence to define the lower limits of the amount of DNA that can be used whilst preserving data quality. Such evidence is valuable for analyses utilizing precious or limited DNA sources. We used a single pooled sample of DNA in quadruplicate at three dilutions to define replicability and noise, and an independent population dataset of 328 individuals (from a community-based study including US-born non-Hispanic Black and white persons) to assess the impact of total DNA input on the quality of data generated using the Illumina Infinium MethylationEPIC BeadChip. We found that data are less reliable and more noisy as DNA input decreases to 40ng, with clear reductions in data quality; and that low DNA input is associated with a reduction in power to detect EWAS associations, requiring larger sample sizes. We conclude that DNA input as low as 40ng can be used with the Illumina Infinium MethylationEPIC BeadChip, provided quality checks and sensitivity analyses are undertaken.

Alternative COVID-19 mitigation measures in school classrooms: analysis using an agent-based model of SARS-CoV-2 transmission

The SARS-CoV-2 epidemic has impacted children's education, with schools required to implement infection control measures that have led to periods of absence and classroom closures. We developed an agent-based epidemiological model of SARS-CoV-2 transmission in a school classroom that allows us to quantify projected infection patterns within primary school classrooms, and related uncertainties. Our approach is based on a contact model constructed using random networks, informed by structured expert judgement. The effectiveness of mitigation strategies in suppressing infection outbreaks and limiting pupil absence are considered. COVID-19 infections in primary schools in England in autumn 2020 were re-examined and the model was then used to estimate infection levels in autumn 2021, as the Delta variant was emerging and it was thought likely that school transmission would play a major role in an incipient new wave of the epidemic. Our results were in good agreement with available data. These findings indicate that testing-based surveillance is more effective than bubble quarantine, both for reducing transmission and avoiding pupil absence, even accounting for insensitivity of self-administered tests. Bubble quarantine entails large numbers of absences, with only modest impact on classroom infections. However, maintaining reduced contact rates within the classroom can have a major benefit for managing COVID-19 in school settings.

DNA methylation signature of chronic low-grade inflammation and its role in cardio-respiratory diseases

We performed a multi-ethnic Epigenome Wide Association study on 22,774 individuals to describe the DNA methylation signature of chronic low-grade inflammation as measured by C-Reactive protein (CRP). We find 1,511 independent differentially methylated loci associated with CRP. These CpG sites show correlation structures across chromosomes, and are primarily situated in euchromatin, depleted in CpG islands. These genomic loci are predominantly situated in transcription factor binding sites and genomic enhancer regions. Mendelian randomization analysis suggests altered CpG methylation is a consequence of increased blood CRP levels. Mediation analysis reveals obesity and smoking as important underlying driving factors for changed CpG methylation. Finally, we find that an activated CpG signature significantly increases the risk for cardiometabolic diseases and COPD.

Complex trait methylation scores in the prediction of major depressive disorder

BACKGROUND

DNA methylation (DNAm) is associated with time-varying environmental factors that contribute to major depressive disorder (MDD) risk. We sought to test whether DNAm signatures of lifestyle and biochemical factors were associated with MDD to reveal dynamic biomarkers of MDD risk that may be amenable to lifestyle interventions.

METHODS

Here, we calculated methylation scores (MS) at multiple p-value thresholds for lifestyle (BMI, smoking, alcohol consumption, and educational attainment) and biochemical (high-density lipoprotein (HDL) and total cholesterol) factors in Generation Scotland (GS) (N=9,502) and in a replication cohort (ALSPAC_adults, N=565), using CpG sites reported in previous well-powered methylome-wide association studies. We also compared their predictive accuracy for MDD to a MDD MS in an independent GS sub-sample (N=4,432).

FINDINGS

Each trait MS was significantly associated with its corresponding phenotype in GS (β_range=0.089-1.457) and in ALSPAC (β_range=0.078-2.533). Each MS was also significantly associated with MDD before and after adjustment for its corresponding phenotype in GS (β_range=0.053-0.145). After accounting for relevant lifestyle factors, MS for educational attainment (β=0.094) and alcohol consumption (MS_p-value<0.01-0.5; β_range=-0.069-0.083) remained significantly associated with MDD in GS. Smoking (AUC=0.569) and educational attainment (AUC=0.585) MSs could discriminate MDD from controls better than the MDD MS (AUC=0.553) in the independent GS sub-sample. Analyses implicating MDD did not replicate across ALSPAC, although the direction of effect was consistent for all traits when adjusting for the MS corresponding phenotypes.

INTERPRETATION

We showed that lifestyle and biochemical MS were associated with MDD before and after adjustment for their corresponding phenotypes (p_nominal<0.05), but not when smoking, alcohol consumption, and BMI were also included as covariates. MDD results did not replicate in the smaller, female-only independent ALSPAC cohort (N_ALSPAC=565; N_GS=9,502), potentially due to demographic differences or low statistical power, but effect sizes were consistent with the direction reported in GS. DNAm scores for modifiable MDD risk factors may contribute to disease vulnerability and, in some cases, explain additional variance to their observed phenotypes.

FUNDING

Wellcome Trust.

DNA methylation changes associated with prenatal mercury exposure: A meta-analysis of prospective cohort studies from PACE consortium

Mercury (Hg) is a ubiquitous heavy metal that originates from both natural and anthropogenic sources and is transformed in the environment to its most toxicant form, methylmercury (MeHg). Recent studies suggest that MeHg exposure can alter epigenetic modifications during embryogenesis. In this study, we examined associations between prenatal MeHg exposure and levels of cord blood DNA methylation (DNAm) by meta-analysis in up to seven independent studies (n = 1462) as well as persistence of those relationships in blood from 7 to 8 year-old children (n = 794). In cord blood, we found limited evidence of differential DNAm at cg24184221 in MED31 (β = 2.28 × 10-4, p-value = 5.87 × 10-5) in relation to prenatal MeHg exposure. In child blood, we identified differential DNAm at cg15288800 (β = 0.004, p-value = 4.97 × 10-5), also located in MED31. This repeated link to MED31, a gene involved in lipid metabolism and RNA Polymerase II transcription function, may suggest a DNAm perturbation related to MeHg exposure that persists into early childhood. Further, we found evidence for association between prenatal MeHg exposure and child blood DNAm levels at two additional CpGs: cg12204245 (β = 0.002, p-value = 4.81 × 10-7) in GRK1 and cg02212000 (β = -0.001, p-value = 8.13 × 10-7) in GGH. Prenatal MeHg exposure was associated with DNAm modifications that may influence health outcomes, such as cognitive or anthropometric development, in different populations.

Epigenetic Regulation of F2RL3 Associates With Myocardial Infarction and Platelet Function

BACKGROUND

DNA hypomethylation at the F2RL3 (F2R like thrombin or trypsin receptor 3) locus has been associated with both smoking and atherosclerotic cardiovascular disease; whether these smoking-related associations form a pathway to disease is unknown. F2RL3 encodes protease-activated receptor 4, a potent thrombin receptor expressed on platelets. Given the role of thrombin in platelet activation and the role of thrombus formation in myocardial infarction, alterations to this biological pathway could be important for ischemic cardiovascular disease.

METHODS

We conducted multiple independent experiments to assess whether DNA hypomethylation at F2RL3 in response to smoking is associated with risk of myocardial infarction via changes to platelet reactivity. Using cohort data (N=3205), we explored the relationship between smoking, DNA hypomethylation at F2RL3, and myocardial infarction. We compared platelet reactivity in individuals with low versus high DNA methylation at F2RL3 (N=41). We used an in vitro model to explore the biological response of F2RL3 to cigarette smoke extract. Finally, a series of reporter constructs were used to investigate how differential methylation could impact F2RL3 gene expression.

RESULTS

Observationally, DNA methylation at F2RL3 mediated an estimated 34% of the smoking effect on increased risk of myocardial infarction. An association between methylation group (low/high) and platelet reactivity was observed in response to PAR4 (protease-activated receptor 4) stimulation. In cells, cigarette smoke extract exposure was associated with a 4.9% to 9.3% reduction in DNA methylation at F2RL3 and a corresponding 1.7-(95% CI, 1.2-2.4, P=0.04) fold increase in F2RL3 mRNA. Results from reporter assays suggest the exon 2 region of F2RL3 may help control gene expression.

CONCLUSIONS

Smoking-induced epigenetic DNA hypomethylation at F2RL3 appears to increase PAR4 expression with potential downstream consequences for platelet reactivity. Combined evidence here not only identifies F2RL3 DNA methylation as a possible contributory pathway from smoking to cardiovascular disease risk but from any feature potentially influencing F2RL3 regulation in a similar manner.

Epigenetic biomarkers of ageing are predictive of mortality risk in a longitudinal clinical cohort of individuals diagnosed with oropharyngeal cancer

BACKGROUND

Epigenetic clocks are biomarkers of ageing derived from DNA methylation levels at a subset of CpG sites. The difference between age predicted by these clocks and chronological age, termed "epigenetic age acceleration", has been shown to predict age-related disease and mortality. We aimed to assess the prognostic value of epigenetic age acceleration and a DNA methylation-based mortality risk score with all-cause mortality in a prospective clinical cohort of individuals with head and neck cancer: Head and Neck 5000. We investigated two markers of intrinsic epigenetic age acceleration (IEAAHorvath and IEAAHannum), one marker of extrinsic epigenetic age acceleration (EEAA), one optimised to predict physiological dysregulation (AgeAccelPheno), one optimised to predict lifespan (AgeAccelGrim) and a DNA methylation-based predictor of mortality (ZhangScore). Cox regression models were first used to estimate adjusted hazard ratios (HR) and 95% confidence intervals (CI) for associations of epigenetic age acceleration with all-cause mortality in people with oropharyngeal cancer (n = 408; 105 deaths). The added prognostic value of epigenetic markers compared to a clinical model including age, sex, TNM stage and HPV status was then evaluated.

RESULTS

IEAAHannum and AgeAccelGrim were associated with mortality risk after adjustment for clinical and lifestyle factors (HRs per standard deviation [SD] increase in age acceleration = 1.30 [95% CI 1.07, 1.57; p = 0.007] and 1.40 [95% CI 1.06, 1.83; p = 0.016], respectively). There was weak evidence that the addition of AgeAccelGrim to the clinical model improved 3-year mortality prediction (area under the receiver operating characteristic curve: 0.80 vs. 0.77; p value for difference = 0.069).

CONCLUSION

In the setting of a large, clinical cohort of individuals with head and neck cancer, our study demonstrates the potential of epigenetic markers of ageing to enhance survival prediction in people with oropharyngeal cancer, beyond established prognostic factors. Our findings have potential uses in both clinical and non-clinical contexts: to aid treatment planning and improve patient stratification.

Meta-analysis of epigenome-wide association studies in newborns and children show widespread sex differences in blood DNA methylation

BACKGROUND

Among children, sex-specific differences in disease prevalence, age of onset, and susceptibility have been observed in health conditions including asthma, immune response, metabolic health, some pediatric and adult cancers, and psychiatric disorders. Epigenetic modifications such as DNA methylation may play a role in the sexual differences observed in diseases and other physiological traits.

METHODS

We performed a meta-analysis of the association of sex and cord blood DNA methylation at over 450,000 CpG sites in 8438 newborns from 17 cohorts participating in the Pregnancy And Childhood Epigenetics (PACE) Consortium. We also examined associations of child sex with DNA methylation in older children ages 5.5-10 years from 8 cohorts (n = 4268).

RESULTS

In newborn blood, sex was associated at Bonferroni level significance with differences in DNA methylation at 46,979 autosomal CpG sites (p < 1.3 × 10-7) after adjusting for white blood cell proportions and batch. Most of those sites had lower methylation levels in males than in females. Of the differentially methylated CpG sites identified in newborn blood, 68% (31,727) met look-up level significance (p < 1.1 × 10-6) in older children and had methylation differences in the same direction.

CONCLUSIONS

This is a large-scale meta-analysis examining sex differences in DNA methylation in newborns and older children. Expanding upon previous studies, we replicated previous findings and identified additional autosomal sites with sex-specific differences in DNA methylation. Differentially methylated sites were enriched in genes involved in cancer, psychiatric disorders, and cardiovascular phenotypes.

High COVID-19 transmission potential associated with re-opening universities can be mitigated with layered interventions

Controlling COVID-19 transmission in universities poses challenges due to the complex social networks and potential for asymptomatic spread. We developed a stochastic transmission model based on realistic mixing patterns and evaluated alternative mitigation strategies. We predict, for plausible model parameters, that if asymptomatic cases are half as infectious as symptomatic cases, then 15% (98% Prediction Interval: 6-35%) of students could be infected during the first term without additional control measures. First year students are the main drivers of transmission with the highest infection rates, largely due to communal residences. In isolation, reducing face-to-face teaching is the most effective intervention considered, however layering multiple interventions could reduce infection rates by 75%. Fortnightly or more frequent mass testing is required to impact transmission and was not the most effective option considered. Our findings suggest that additional outbreak control measures should be considered for university settings.

Do Mass Spectrometry-Derived Metabolomics Improve the Prediction of Pregnancy-Related Disorders? Findings from a UK Birth Cohort with Independent Validation

Many women who experience gestational diabetes (GDM), gestational hypertension (GHT), pre-eclampsia (PE), have a spontaneous preterm birth (sPTB) or have an offspring born small/large for gestational age (SGA/LGA) do not meet the criteria for high-risk pregnancies based upon certain maternal risk factors. Tools that better predict these outcomes are needed to tailor antenatal care to risk. Recent studies have suggested that metabolomics may improve the prediction of these pregnancy-related disorders. These have largely been based on targeted platforms or focused on a single pregnancy outcome. The aim of this study was to assess the predictive ability of an untargeted platform of over 700 metabolites to predict the above pregnancy-related disorders in two cohorts. We used data collected from women in the Born in Bradford study (BiB; two sub-samples, n = 2000 and n = 1000) and the Pregnancy Outcome Prediction study (POPs; n = 827) to train, test and validate prediction models for GDM, PE, GHT, SGA, LGA and sPTB. We compared the predictive performance of three models: (1) risk factors (maternal age, pregnancy smoking, BMI, ethnicity and parity) (2) mass spectrometry (MS)-derived metabolites (n = 718 quantified metabolites, collected at 26-28 weeks' gestation) and (3) combined risk factors and metabolites. We used BiB for the training and testing of the models and POPs for independent validation. In both cohorts, discrimination for GDM, PE, LGA and SGA improved with the addition of metabolites to the risk factor model. The models' area under the curve (AUC) were similar for both cohorts, with good discrimination for GDM (AUC (95% CI) BiB 0.76 (0.71, 0.81) and POPs 0.76 (0.72, 0.81)) and LGA (BiB 0.86 (0.80, 0.91) and POPs 0.76 (0.60, 0.92)). Discrimination was improved for the combined models (compared to the risk factors models) for PE and SGA, with modest discrimination in both studies (PE-BiB 0.68 (0.58, 0.78) and POPs 0.66 (0.60, 0.71); SGA-BiB 0.68 (0.63, 0.74) and POPs 0.64 (0.59, 0.69)). Prediction for sPTB was poor in BiB and POPs for all models. In BiB, calibration for the combined models was good for GDM, LGA and SGA. Retained predictors include 4-hydroxyglutamate for GDM, LGA and PE and glycerol for GDM and PE. MS-derived metabolomics combined with maternal risk factors improves the prediction of GDM, PE, LGA and SGA, with good discrimination for GDM and LGA. Validation across two very different cohorts supports further investigation on whether the metabolites reflect novel causal paths to GDM and LGA.

BMI trajectory in childhood is associated with asthma incidence at young adulthood mediated by DNA methylation

PURPOSE

Body mass index (BMI) is associated with asthma but associations of BMI temporal patterns with asthma incidence are unclear. Previous studies suggest that DNA methylation (DNAm) is associated with asthma status and variation in DNAm is a consequence of BMI changes. This study assessed the direct and indirect (via DNAm) effects of BMI trajectories in childhood on asthma incidence at young adulthood.

METHODS

Data from the Isle of Wight (IoW) birth cohort were included in the analyses. Group-based trajectory modelling was applied to infer latent BMI trajectories from ages 1 to 10 years. An R package, ttscreening, was applied to identify differentially methylated CpGs at age 10 years associated with BMI trajectories, stratified for sex. Logistic regressions were used to further exclude CpGs with DNAm at age 10 years not associated with asthma incidence at 18 years. CpGs discovered via path analyses that mediated the association of BMI trajectories with asthma incidence in the IoW cohort were further tested in an independent cohort, the Avon Longitudinal Study of Children and Parents (ALSPAC).

RESULTS

Two BMI trajectories (high vs. normal) were identified. Of the 442,474 CpG sites, DNAm at 159 CpGs in males and 212 in females were potentially associated with BMI trajectories. Assessment of their association with asthma incidence identified 9 CpGs in males and 6 CpGs in females. DNAm at 4 of these 15 CpGs showed statistically significant mediation effects (p-value < 0.05). At two of the 4 CpGs (cg23632109 and cg10817500), DNAm completely mediated the association (i.e., only statistically significant indirect effects were identified). In the ALSPAC cohort, at all four CpGs, the same direction of mediating effects were observed as those found in the IoW cohort, although statistically insignificant.

CONCLUSION

The association of BMI trajectory in childhood with asthma incidence at young adulthood is possibly mediated by DNAm.

Early-childhood cytomegalovirus infection and children's neurocognitive development

BACKGROUND

Despite a clear association seen in congenitally infected children, the effect of postnatal cytomegalovirus (CMV) infection during early childhood on cognitive development has not yet been determined.

METHODS

CMV-infection status was obtained based on serological measurements when children were 7 years old. Using population-based longitudinal data, we employed multivariate Poisson regression with a robust variance estimator to characterize the relationship between childhood CMV infection and adverse neurocognitive outcomes in children. Suboptimal neurocognitive outcomes were compared between CMV-positive and CMV-negative children using various cognitive assessments from 8 to 15 years of age. Children were evaluated on the cognitive domains of language, reading, memory and general intelligence, with a suboptimal score being >2 standard deviations lower than the mean score. Approximate Bayes factor (ABF) analysis was used to determine the level of evidence for the observed associations.

RESULTS

With adjustment for potential confounders, we observed that early-childhood CMV infection was associated with suboptimal total intelligence quotient (IQ) at 8 years of age [incidence-rate ratio (IRR) = 2.50, 95% confidence interval (CI) 1.35-4.62, ABF = 0.08], but not with suboptimal total IQ at 15 years of age (IRR = 0.97, 95% CI 0.43-2.19, ABF = 1.68). Suboptimal attentional control at 8 years (IRR = 1.74, 95% CI 1.13-2.68, ABF = 0.18) and reading comprehension at 9 years (IRR = 1.93, 95% CI 1.12-3.33, ABF = 0.24) were also associated with CMV infection. ABF analysis provided strong evidence for the association between CMV infection and total IQ at 8 years, and only anecdotal evidence for attentional control at 8 years and reading comprehension at 9 years. All other cognitive measures assessed were not associated with CMV infection.

CONCLUSION

In this large-scale prospective cohort, we observed some evidence for adverse neurocognitive effects of postnatal CMV infection on general intelligence during early childhood, although not with lasting effect. If confirmed, these results could support the implementation of preventative measures to combat postnatal CMV infection.

Contacts and behaviours of university students during the COVID-19 pandemic at the start of the 2020/2021 academic year

University students have unique living, learning and social arrangements which may have implications for infectious disease transmission. To address this data gap, we created CONQUEST (COroNavirus QUESTionnaire), a longitudinal online survey of contacts, behaviour, and COVID-19 symptoms for University of Bristol (UoB) staff/students. Here, we analyse results from 740 students providing 1261 unique records from the start of the 2020/2021 academic year (14/09/2020-01/11/2020), where COVID-19 outbreaks led to the self-isolation of all students in some halls of residences. Although most students reported lower daily contacts than in pre-COVID-19 studies, there was heterogeneity, with some reporting many (median = 2, mean = 6.1, standard deviation = 15.0; 8% had ≥ 20 contacts). Around 40% of students' contacts were with individuals external to the university, indicating potential for transmission to non-students/staff. Only 61% of those reporting cardinal symptoms in the past week self-isolated, although 99% with a positive COVID-19 test during the 2 weeks before survey completion had self-isolated within the last week. Some students who self-isolated had many contacts (mean = 4.3, standard deviation = 10.6). Our results provide context to the COVID-19 outbreaks seen in universities and are available for modelling future outbreaks and informing policy.

Neonatal DNA methylation and childhood low prosocial behavior: An epigenome-wide association meta-analysis

Low prosocial behavior in childhood has been consistently linked to later psychopathology, with evidence supporting the influence of both genetic and environmental factors on its development. Although neonatal DNA methylation (DNAm) has been found to prospectively associate with a range of psychological traits in childhood, its potential role in prosocial development has yet to be investigated. This study investigated prospective associations between cord blood DNAm at birth and low prosocial behavior within and across four longitudinal birth cohorts from the Pregnancy And Childhood Epigenetics (PACE) Consortium. We examined (a) developmental trajectories of "chronic-low" versus "typical" prosocial behavior across childhood in a case-control design (N = 2,095), and (b) continuous "low prosocial" scores at comparable cross-cohort time-points (N = 2,121). Meta-analyses were performed to examine differentially methylated positions and regions. At the cohort-specific level, three CpGs were found to associate with chronic low prosocial behavior; however, none of these associations was replicated in another cohort. Meta-analysis revealed no epigenome-wide significant CpGs or regions. Overall, we found no evidence for associations between DNAm patterns at birth and low prosocial behavior across childhood. Findings highlight the importance of employing multi-cohort approaches to replicate epigenetic associations and reduce the risk of false positive discoveries.

DNA methylation signatures of aggression and closely related constructs: A meta-analysis of epigenome-wide studies across the lifespan

DNA methylation profiles of aggressive behavior may capture lifetime cumulative effects of genetic, stochastic, and environmental influences associated with aggression. Here, we report the first large meta-analysis of epigenome-wide association studies (EWAS) of aggressive behavior (N = 15,324 participants). In peripheral blood samples of 14,434 participants from 18 cohorts with mean ages ranging from 7 to 68 years, 13 methylation sites were significantly associated with aggression (alpha = 1.2 × 10-7; Bonferroni correction). In cord blood samples of 2425 children from five cohorts with aggression assessed at mean ages ranging from 4 to 7 years, 83% of these sites showed the same direction of association with childhood aggression (r = 0.74, p = 0.006) but no epigenome-wide significant sites were found. Top-sites (48 at a false discovery rate of 5% in the peripheral blood meta-analysis or in a combined meta-analysis of peripheral blood and cord blood) have been associated with chemical exposures, smoking, cognition, metabolic traits, and genetic variation (mQTLs). Three genes whose expression levels were associated with top-sites were previously linked to schizophrenia and general risk tolerance. At six CpGs, DNA methylation variation in blood mirrors variation in the brain. On average 44% (range = 3-82%) of the aggression-methylation association was explained by current and former smoking and BMI. These findings point at loci that are sensitive to chemical exposures with potential implications for neuronal functions. We hope these results to be a starting point for studies leading to applications as peripheral biomarkers and to reveal causal relationships with aggression and related traits.

Epigenome-wide association meta-analysis of DNA methylation with coffee and tea consumption

Coffee and tea are extensively consumed beverages worldwide which have received considerable attention regarding health. Intake of these beverages is consistently linked to, among others, reduced risk of diabetes and liver diseases; however, the mechanisms of action remain elusive. Epigenetics is suggested as a mechanism mediating the effects of dietary and lifestyle factors on disease onset. Here we report the results from epigenome-wide association studies (EWAS) on coffee and tea consumption in 15,789 participants of European and African-American ancestries from 15 cohorts. EWAS meta-analysis of coffee consumption reveals 11 CpGs surpassing the epigenome-wide significance threshold (P-value <1.1×10-7), which annotated to the AHRR, F2RL3, FLJ43663, HDAC4, GFI1 and PHGDH genes. Among them, cg14476101 is significantly associated with expression of the PHGDH and risk of fatty liver disease. Knockdown of PHGDH expression in liver cells shows a correlation with expression levels of genes associated with circulating lipids, suggesting a role of PHGDH in hepatic-lipid metabolism. EWAS meta-analysis on tea consumption reveals no significant association, only two CpGs annotated to CACNA1A and PRDM16 genes show suggestive association (P-value <5.0×10-6). These findings indicate that coffee-associated changes in DNA methylation levels may explain the mechanism of action of coffee consumption in conferring risk of diseases.

An exploration of the genetic epidemiology of non-suicidal self-harm and suicide attempt

BACKGROUND

Empirical evidence supporting the distinction between suicide attempt (SA) and non-suicidal self-harm (NSSH) is lacking. Although NSSH is a risk factor for SA, we do not currently know whether these behaviours lie on a continuum of severity, or whether they are discrete outcomes with different aetiologies. We conducted this exploratory genetic epidemiology study to investigate this issue further.

METHODS

We explored the extent of genetic overlap between NSSH and SA in a large, richly-phenotyped cohort (the Avon Longitudinal Study of Parents and Children; N = 4959), utilising individual-level genetic and phenotypic data to conduct analyses of genome-wide complex traits and polygenic risk scores (PRS).

RESULTS

The single nucleotide polymorphism heritability of NSSH was estimated to be 13% (SE 0.07) and that of SA to be 0% (SE 0.07). Of the traits investigated, NSSH was most strongly correlated with higher IQ (rG = 0.31, SE = 0.22), there was little evidence of high genetic correlation between NSSH and SA (rG = - 0.1, SE = 0.54), likely due to the low heritability estimate for SA. The PRS for depression differentiated between those with NSSH and SA in multinomial regression. The optimal PRS prediction model for SA (Nagelkerke R2 0.022, p < 0.001) included ADHD, depression, income, anorexia and neuroticism and explained more variance than the optimal prediction model for NSSH (Nagelkerke R2 0.010, p < 0.001) which included ADHD, alcohol consumption, autism spectrum conditions, depression, IQ, neuroticism and suicide attempt.

CONCLUSIONS

Our findings suggest that SA does not have a large genetic component, and that although NSSH and SA are not discrete outcomes there appears to be little genetic overlap between the two. The relatively small sample size and resulting low heritability estimate for SA was a limitation of the study. Combined with low heritability estimates, this implies that family or population structures in SA GWASs may contribute to signals detected.

DNA methylation at birth is associated with lung function development until age 26 years

Little is known about whether DNA methylation (DNAm) of cytosine-phosphate-guanine (CpG) sites at birth predicts patterns of lung function development. We used heel prick DNAm from the F1-generation of Isle of Wight birth cohort (IOWBC-F1) for discovery of CpGs associated with lung function trajectories (forced expiratory volume in 1 s, forced vital capacity, their ratio, and forced expiratory flow at 25-75% of forced vital capacity) over the first 26 years, stratified by sex. We replicated the findings in the Avon Longitudinal Study of Parents and Children (ALSPAC) using cord blood DNAm.Epigenome-wide screening was applied to identify CpGs associated with lung function trajectories in 396 boys and 390 girls of IOWBC-F1. Replication in ALSPAC focussed on lung function at ages 8, 15 and 24 years. Statistically significantly replicated CpGs were investigated for consistency in direction of association between cohorts, stability of DNAm over time in IOWBC-F1, relevant biological processes and for association with gene expression (n=161) in IOWBC F2-generation (IOWBC-F2).Differential DNAm of eight CpGs on genes GLUL, MYCN, HLX, LHX1, COBL, COL18A1, STRA6, and WNT11 involved in developmental processes, were significantly associated with lung function in the same direction in IOWBC-F1 and ALSPAC, and showed stable patterns at birth, aged 10 and 18 years between high and low lung function trajectories in IOWBC-F1. CpGs on LHX1 and COL18A1 were linked to gene expression in IOWBC-F2.In two large cohorts, novel DNAm at birth were associated with patterns of lung function in adolescence and early adulthood providing possible targets for preventative interventions against adverse pulmonary function development.

Modelling pooling strategies for SARS-CoV-2 testing in a university setting

Background: Pre-symptomatic and asymptomatic transmission of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) are important elements in the coronavirus disease 2019 (COVID-19) pandemic, and there remains a reliance on testing to manage the spread of the disease. In the UK, many universities opened for blended learning for the 2020-2021 academic year, with a mixture of face to face and online teaching. Methods: In this study we present a simulation framework to evaluate the effectiveness of different mass testing strategies within a university setting, across a range of transmission scenarios. Results: The sensitivity of 5x pooled RT-qPCR tests appears to be higher than testing using the lateral flow device with relatively little loss compared to single RT-qPCR tests, and is improved by pooling by social cluster. The range of strategies that we evaluated give comparable results for estimating prevalence. Conclusions: Pooling tests by known social structures, such as student households can substantially improve the cost effectiveness of RT-qPCR tests. We also note that routine recording of quantitative RT-qPCR results would facilitate future modelling studies.

Pre-adolescence DNA methylation is associated with BMI status change from pre- to post-adolescence

BACKGROUND

Previous studies have shown that DNA methylation (DNAm) is associated with body mass index (BMI). However, it is unknown whether DNAm at pre-adolescence is associated with BMI status transition from pre- to post-adolescence. In the Isle of Wight (IoW) birth cohort, genome-wide DNA methylation in whole blood was measured using Illumina Infinium Human450 and EPIC BeadChip arrays in n = 325 subjects, and pre- to post-adolescence BMI transition was classified into four groups: (1) normal to normal, (2) normal to overweight or obese, (3) overweight or obese to normal, and (4) persistent overweight or obese. We used recursive random forest to screen genome-wide Cytosine-phosphate-Guanine (CpG) sites with DNAm potentially associated with BMI transition for each gender, and the association of BMI status transition with DNAm at an earlier age was assessed via logistic regressions. To evaluate gender specificity, interactions between DNAm and gender were included in the model. Findings in the IoW cohort were further tested in an independent cohort, the Avon Longitudinal Study of Parents and Children (ALSPAC).

RESULTS

In total, 174 candidate CpGs were selected including CpGs from screening and CpGs previously associated correctionally with BMI in children and adults. Of these 174 CpGs, pre-adolescent DNAm of 38 CpGs in the IoW cohort was associated with BMI status transition, including 30 CpGs showing gender-specific associations. Thirteen CpGs showed consistent associations between the IoW cohort and the ALSPAC cohort (11 of which were gender-specific).

CONCLUSION

Pre-adolescence DNAm is associated with the change in BMI status from pre- to post-adolescence and such associations are likely to be gender-specific.

Childhood adversity, pubertal timing and self-harm: a longitudinal cohort study

BACKGROUND

The occurrence of early childhood adversity is strongly linked to later self-harm, but there is poor understanding of how this distal risk factor might influence later behaviours. One possible mechanism is through an earlier onset of puberty in children exposed to adversity, since early puberty is associated with an increased risk of adolescent self-harm. We investigated whether early pubertal timing mediates the association between childhood adversity and later self-harm.

METHODS

Participants were 6698 young people from a UK population-based birth cohort (ALSPAC). We measured exposure to nine types of adversity from 0 to 9 years old, and self-harm when participants were aged 16 and 21 years. Pubertal timing measures were age at peak height velocity (aPHV - males and females) and age at menarche (AAM). We used generalised structural equation modelling for analyses.

RESULTS

For every additional type of adversity; participants had an average 12-14% increased risk of self-harm by 16. Relative risk (RR) estimates were stronger for direct effects when outcomes were self-harm with suicidal intent. There was no evidence that earlier pubertal timing mediated the association between adversity and self-harm [indirect effect RR 1.00, 95% confidence interval (CI) 1.00-1.00 for aPHV and RR 1.00, 95% CI 1.00-1.01 for AAM].

CONCLUSIONS

A cumulative measure of exposure to multiple types of adversity does not confer an increased risk of self-harm via early pubertal timing, however both childhood adversity and early puberty are risk factors for later self-harm. Research identifying mechanisms underlying the link between childhood adversity and later self-harm is needed to inform interventions.

Novel DNA methylation signatures of tobacco smoking with trans-ethnic effects

BACKGROUND

Smoking remains one of the leading preventable causes of death. Smoking leaves a strong signature on the blood methylome as shown in multiple studies using the Infinium HumanMethylation450 BeadChip. Here, we explore novel blood methylation smoking signals on the Illumina MethylationEPIC BeadChip (EPIC) array, which also targets novel CpG-sites in enhancers.

METHOD

A smoking-methylation meta-analysis was carried out using EPIC DNA methylation profiles in 1407 blood samples from four UK population-based cohorts, including the MRC National Survey for Health and Development (NSHD) or 1946 British birth cohort, the National Child Development Study (NCDS) or 1958 birth cohort, the 1970 British Cohort Study (BCS70), and the TwinsUK cohort (TwinsUK). The overall discovery sample included 269 current, 497 former, and 643 never smokers. Replication was pursued in 3425 trans-ethnic samples, including 2325 American Indian individuals participating in the Strong Heart Study (SHS) in 1989-1991 and 1100 African-American participants in the Genetic Epidemiology Network of Arteriopathy Study (GENOA).

RESULTS

Altogether 952 CpG-sites in 500 genes were differentially methylated between smokers and never smokers after Bonferroni correction. There were 526 novel smoking-associated CpG-sites only profiled by the EPIC array, of which 486 (92%) replicated in a meta-analysis of the American Indian and African-American samples. Novel CpG sites mapped both to genes containing previously identified smoking-methylation signals and to 80 novel genes not previously linked to smoking, with the strongest novel signal in SLAMF7. Comparison of former versus never smokers identified that 37 of these sites were persistently differentially methylated after cessation, where 16 represented novel signals only profiled by the EPIC array. We observed a depletion of smoking-associated signals in CpG islands and an enrichment in enhancer regions, consistent with previous results.

CONCLUSION

This study identified novel smoking-associated signals as possible biomarkers of exposure to smoking and may help improve our understanding of smoking-related disease risk.

Sex-specific associations of asthma acquisition with changes in DNA methylation during adolescence

BACKGROUND

Underlying biological mechanisms involved in sex differences in asthma status changes from pre- to post-adolescence are unclear. DNA methylation (DNAm) has been shown to be associated with the risk of asthma.

OBJECTIVE

We hypothesized that asthma acquisition from pre- to post-adolescence was associated with changes in DNAm during this period at asthma-associated cytosine-phosphate-guanine (CpG) sites and such an association was sex-specific.

METHODS

Subjects from the Isle of Wight birth cohort (IOWBC) with DNAm in blood at ages 10 and 18 years (n = 124 females, 151 males) were studied. Using a training-testing approach, epigenome-wide CpGs associated with asthma were identified. Logistic regression was used to examine sex-specific associations of DNAm changes with asthma acquisition between ages 10 and 18 at asthma-associated CpGs. The ALSPAC birth cohort was used for independent replication. To assess functional relevance of identified CpGs, association of DNAm with gene expression in blood was assessed.

RESULTS

We identified 535 CpGs potentially associated with asthma. Significant interaction effects of DNAm changes and sex on asthma acquisition in adolescence were found at 13 of the 535 CpGs in IOWBC (P-values <1.0 × 10^-3 ). In the replication cohort, consistent interaction effects were observed at 10 of the 13 CpGs. At 7 of these 10 CpGs, opposite DNAm changes across adolescence were observed between sexes in both cohorts. In both cohorts, cg20891917, located on IFRD1 linked to asthma, shows strong sex-specific effects on asthma transition (P-values <.01 in both cohorts).

CONCLUSION AND CLINICAL RELEVANCE

Gender reversal in asthma acquisition is associated with opposite changes in DNAm (males vs females) from pre- to post-adolescence at asthma-associated CpGs. These CpGs are potential biomarkers of sex-specific asthma acquisition in adolescence.

Do nuclear magnetic resonance (NMR)-based metabolomics improve the prediction of pregnancy-related disorders? Findings from a UK birth cohort with independent validation

BACKGROUND

Prediction of pregnancy-related disorders is usually done based on established and easily measured risk factors. Recent advances in metabolomics may provide earlier and more accurate prediction of women at risk of pregnancy-related disorders.

METHODS

We used data collected from women in the Born in Bradford (BiB; n = 8212) and UK Pregnancies Better Eating and Activity Trial (UPBEAT; n = 859) studies to create and validate prediction models for pregnancy-related disorders. These were gestational diabetes mellitus (GDM), hypertensive disorders of pregnancy (HDP), small for gestational age (SGA), large for gestational age (LGA) and preterm birth (PTB). We used ten-fold cross-validation and penalised regression to create prediction models. We compared the predictive performance of (1) risk factors (maternal age, pregnancy smoking, body mass index (BMI), ethnicity and parity) to (2) nuclear magnetic resonance-derived metabolites (N = 156 quantified metabolites, collected at 24-28 weeks gestation) and (3) combined risk factors and metabolites. The multi-ethnic BiB cohort was used for training and testing the models, with independent validation conducted in UPBEAT, a multi-ethnic study of obese pregnant women.

RESULTS

Maternal age, pregnancy smoking, BMI, ethnicity and parity were retained in the combined risk factor and metabolite models for all outcomes apart from PTB, which did not include maternal age. In addition, 147, 33, 96, 51 and 14 of the 156 metabolite traits were retained in the combined risk factor and metabolite model for GDM, HDP, SGA, LGA and PTB, respectively. These include cholesterol and triglycerides in very low-density lipoproteins (VLDL) in the models predicting GDM, HDP, SGA and LGA, and monounsaturated fatty acids (MUFA), ratios of MUFA to omega 3 fatty acids and total fatty acids, and a ratio of apolipoprotein B to apolipoprotein A-1 (APOA:APOB1) were retained predictors for GDM and LGA. In BiB, discrimination for GDM, HDP, LGA and SGA was improved in the combined risk factors and metabolites models. Risk factor area under the curve (AUC 95% confidence interval (CI)): GDM (0.69 (0.64, 0.73)), HDP (0.74 (0.70, 0.78)) and LGA (0.71 (0.66, 0.75)), and SGA (0.59 (0.56, 0.63)). Combined risk factor and metabolite models AUC 95% (CI): GDM (0.78 (0.74, 0.81)), HDP (0.76 (0.73, 0.79)) and LGA (0.75 (0.70, 0.79)), and SGA (0.66 (0.63, 0.70)). For GDM, HDP and LGA, but not SGA, calibration was good for a combined risk factor and metabolite model. Prediction of PTB was poor for all models. Independent validation in UPBEAT at 24-28 weeks and 15-18 weeks gestation confirmed similar patterns of results, but AUCs were attenuated.

CONCLUSIONS

Our results suggest a combined risk factor and metabolite model improves prediction of GDM, HDP and LGA, and SGA, when compared to risk factors alone. They also highlight the difficulty of predicting PTB, with all models performing poorly.

Fresh street: the development and feasibility of a place-based, subsidy for fresh fruit and vegetables

Background

Many UK communities experience food insecurity, and consume diets high in energy-dense, nutrient poor, processed foods and low in fruit and vegetables (FV). We explored a novel area-based approach to promote FV consumption and healthy eating in one such community.

Methods

We developed a weekly subsidy scheme for fresh FV with key local stakeholders in an area of socioeconomic deprivation in Northern England. The scheme (Fresh Street) offered five £1 vouchers to every household, regardless of income or household type. Vouchers were redeemable with local suppliers of fresh FV (not supermarkets). The feasibility of the scheme was assessed in four streets using rapid ethnographic assessment and voucher redemption information.

Results

Local councillors and public health teams were supportive of the scheme. Most eligible households joined the scheme (n = 80/97, 83%), and 89.3% (17 849/19 982) of vouchers issued were redeemed. Householders reported that the scheme made them think about what they were eating, and prompted them to buy and eat more FV.

Conclusions

This feasibility study reported high levels of acceptance for a place-based, household-level weekly FV subsidy scheme. Further research is required to evaluate the effectiveness of this approach to creating healthy diets, eating behaviours and food systems.

Population neuroimaging: generation of a comprehensive data resource within the ALSPAC pregnancy and birth cohort

Neuroimaging offers a valuable insight into human brain development by allowing in vivo assessment of structure, connectivity and function. Multimodal neuroimaging data have been obtained as part of three sub-studies within the Avon Longitudinal Study of Parents and Children, a prospective multigenerational pregnancy and birth cohort based in the United Kingdom. Brain imaging data were acquired when offspring were between 18 and 24 years of age, and included acquisition of structural, functional and magnetization transfer magnetic resonance, diffusion tensor, and magnetoencephalography imaging. This resource provides a unique opportunity to combine neuroimaging data with extensive phenotypic and genotypic measures from participants, their mothers, and fathers.

Investigating evidence for a causal association between inflammation and self-harm: A multivariable Mendelian Randomisation study

BACKGROUND

The causal role of inflammatory markers on self-harm and suicidal risk has been studied using observational data, with conflicting results. Confounding and reverse causation can lead to bias, so we appraised question from a genetic perspective to protect against these biases. We measured associations between genetic liability for high levels of inflammatory markers Interleukin-6 (IL-6) and C-reactive protein (CRP) on self-harm, and conducted a secondary analysis restricted to self-harm with suicidal intent.

METHODS

We conducted two sample and multivariable Mendelian randomisation (MR) to assess the effects of IL-6 and CRP on self-harm utilising existing data and conducting new genome wide association studies to instrument IL-6 and CRP, and for the outcome of self-harm.

RESULTS

No single nucleotide polymorphisms (SNPs) reached genome-wide significance for self-harm, however 193 SNPs met suggestive significance levels (p < 5 × 10-6). We found no evidence of an association between our instruments for IL-6 and self-harm in the two-sample MR, however we found an inverse association between instruments for CRP and self-harm, indicating that higher levels of circulating CRP may protect against self-harm (inverse variance weighted OR 0.92, 95%CI 0.84, 1.01, p = 0.08; MR Egger OR 0.86, 95% CI 0.74, 1.00, p = 0.05). The direct effect estimate for IL-6 was slightly smaller in the multivariable MR than in the two sample MR, while the CRP effect estimates were consistent with the two sample MR (OR 0.92, SE 1.05, p = 0.09).

CONCLUSIONS

Our findings are conflicting and indicate that IL-6 and CRP are not robust etiological markers of increased self-harm or suicide risk.

Identifying epigenetic biomarkers of established prognostic factors and survival in a clinical cohort of individuals with oropharyngeal cancer

BACKGROUND

Smoking status, alcohol consumption and HPV infection (acquired through sexual activity) are the predominant risk factors for oropharyngeal cancer and are thought to alter the prognosis of the disease. Here, we conducted single-site and differentially methylated region (DMR) epigenome-wide association studies (EWAS) of these factors, in addition to ∼ 3-year survival, using Illumina Methylation EPIC DNA methylation profiles from whole blood in 409 individuals as part of the Head and Neck 5000 (HN5000) study. Overlapping sites between each factor and survival were then assessed using two-step Mendelian randomization to assess whether methylation at these positions causally affected survival.

RESULTS

Using the MethylationEPIC array in an OPC dataset, we found novel CpG associations with smoking, alcohol consumption and ~ 3-year survival. We found no CpG associations below our multiple testing threshold associated with HPV16 E6 serological response (used as a proxy for HPV infection). CpG site associations below our multiple-testing threshold (PBonferroni < 0.05) for both a prognostic factor and survival were observed at four gene regions: SPEG (smoking), GFI1 (smoking), PPT2 (smoking) and KHDC3L (alcohol consumption). Evidence for a causal effect of DNA methylation on survival was only observed in the SPEG gene region (HR per SD increase in methylation score 1.28, 95% CI 1.14 to 1.43, P 2.12 × 10-05).

CONCLUSIONS

Part of the effect of smoking on survival in those with oropharyngeal cancer may be mediated by methylation at the SPEG gene locus. Replication in data from independent datasets and data from HN5000 with longer follow-up times is needed to confirm these findings.

Testing for causality between systematically identified risk factors and glioma: a Mendelian randomization study

BACKGROUND

Whilst epidemiological studies have provided evidence of associations between certain risk factors and glioma onset, inferring causality has proven challenging. Using Mendelian randomization (MR), we assessed whether associations of 36 reported glioma risk factors showed evidence of a causal relationship.

METHODS

We performed a systematic search of MEDLINE from inception to October 2018 to identify candidate risk factors and conducted a meta-analysis of two glioma genome-wide association studies (5739 cases and 5501 controls) to form our exposure and outcome datasets. MR analyses were performed using genetic variants to proxy for candidate risk factors. We investigated whether risk factors differed by subtype diagnosis (either glioblastoma (n = 3112) or non-glioblastoma (n = 2411)). MR estimates for each risk factor were determined using multiplicative random effects inverse-variance weighting (IVW). Sensitivity analyses investigated potential pleiotropy using MR-Egger regression, the weighted median estimator, and the mode-based estimator. To increase power, trait-specific polygenic risk scores were used to test the association of a genetically predicated increase in each risk factor with glioma onset.

RESULTS

Our systematic search identified 36 risk factors that could be proxied using genetic variants. Using MR, we found evidence that four genetically predicted traits increased risk of glioma, glioblastoma or non-glioblastoma: longer leukocyte telomere length, liability to allergic disease, increased alcohol consumption and liability to childhood extreme obesity (> 3 standard deviations from the mean). Two traits decreased risk of non-glioblastoma cancers: increased low-density lipoprotein cholesterol (LDLc) and triglyceride levels. Our findings were similar across sensitivity analyses that made allowance for pleiotropy (genetic confounding).

CONCLUSIONS

Our comprehensive investigation provides evidence of a causal link between both genetically predicted leukocyte telomere length, allergic disease, alcohol consumption, childhood extreme obesity, and LDLc and triglyceride levels, and glioma. The findings from our study warrant further research to uncover mechanisms that implicate these traits in glioma onset.

Circulating markers of cellular immune activation in prediagnostic blood sample and lung cancer risk in the Lung Cancer Cohort Consortium (LC3)

Cell-mediated immune suppression may play an important role in lung carcinogenesis. We investigated the associations for circulating levels of tryptophan, kynurenine, kynurenine:tryptophan ratio (KTR), quinolinic acid (QA) and neopterin as markers of immune regulation and inflammation with lung cancer risk in 5,364 smoking-matched case-control pairs from 20 prospective cohorts included in the international Lung Cancer Cohort Consortium. All biomarkers were quantified by mass spectrometry-based methods in serum/plasma samples collected on average 6 years before lung cancer diagnosis. Odds ratios (ORs) and 95% confidence intervals (CIs) for lung cancer associated with individual biomarkers were calculated using conditional logistic regression with adjustment for circulating cotinine. Compared to the lowest quintile, the highest quintiles of kynurenine, KTR, QA and neopterin were associated with a 20-30% higher risk, and tryptophan with a 15% lower risk of lung cancer (all ptrend < 0.05). The strongest associations were seen for current smokers, where the adjusted ORs (95% CIs) of lung cancer for the highest quintile of KTR, QA and neopterin were 1.42 (1.15-1.75), 1.42 (1.14-1.76) and 1.45 (1.13-1.86), respectively. A stronger association was also seen for KTR and QA with risk of lung squamous cell carcinoma followed by adenocarcinoma, and for lung cancer diagnosed within the first 2 years after blood draw. This study demonstrated that components of the tryptophan-kynurenine pathway with immunomodulatory effects are associated with risk of lung cancer overall, especially for current smokers. Further research is needed to evaluate the role of these biomarkers in lung carcinogenesis and progression.

Epigenetic prediction of complex traits and mortality in a cohort of individuals with oropharyngeal cancer

BACKGROUND

DNA methylation (DNAm) variation is an established predictor for several traits. In the context of oropharyngeal cancer (OPC), where 5-year survival is ~ 65%, DNA methylation may act as a prognostic biomarker. We examined the accuracy of DNA methylation biomarkers of 4 complex exposure traits (alcohol consumption, body mass index [BMI], educational attainment and smoking status) in predicting all-cause mortality in people with OPC.

RESULTS

DNAm predictors of alcohol consumption, BMI, educational attainment and smoking status were applied to 364 individuals with OPC in the Head and Neck 5000 cohort (HN5000; 19.6% of total OPC cases in the study), followed up for median 3.9 years; inter-quartile range (IQR) 3.3 to 5.2 years (time-to-event-death or censor). The proportion of phenotypic variance explained in each trait was as follows: 16.5% for alcohol consumption, 22.7% for BMI, 0.4% for educational attainment and 51.1% for smoking. We then assessed the relationship between each DNAm predictor and all-cause mortality using Cox proportional-hazard regression analysis. DNAm prediction of smoking was most consistently associated with mortality risk (hazard ratio [HR], 1.38 per standard deviation (SD) increase in smoking DNAm score; 95% confidence interval [CI] 1.04 to 1.83; P 0.025, in a model adjusted for demographic, lifestyle, health and biological variables). Finally, we examined the accuracy of each DNAm predictor of mortality. DNAm predictors explained similar levels of variance in mortality to self-reported phenotypes. Receiver operator characteristic (ROC) curves for the DNAm predictors showed a moderate discrimination of alcohol consumption (area under the curve [AUC] 0.63), BMI (AUC 0.61) and smoking (AUC 0.70) when predicting mortality. The DNAm predictor for education showed poor discrimination (AUC 0.57). Z tests comparing AUCs between self-reported phenotype ROC curves and DNAm score ROC curves did not show evidence for difference between the two (alcohol consumption P 0.41, BMI P 0.62, educational attainment P 0.49, smoking P 0.19).

CONCLUSIONS

In the context of a clinical cohort of individuals with OPC, DNAm predictors for smoking, alcohol consumption, educational attainment and BMI exhibit similar predictive values for all-cause mortality compared to self-reported data. These findings may have translational utility in prognostic model development, particularly where phenotypic data are not available.

Smoking, DNA Methylation, and Lung Function: a Mendelian Randomization Analysis to Investigate Causal Pathways

Whether smoking-associated DNA methylation has a causal effect on lung function has not been thoroughly evaluated. We first investigated the causal effects of 474 smoking-associated CpGs on forced expiratory volume in 1 s (FEV1) in UK Biobank (n = 321,047) by using two-sample Mendelian randomization (MR) and then replicated this investigation in the SpiroMeta Consortium (n = 79,055). Second, we used two-step MR to investigate whether DNA methylation mediates the effect of smoking on FEV1. Lastly, we evaluated the presence of horizontal pleiotropy and assessed whether there is any evidence for shared causal genetic variants between lung function, DNA methylation, and gene expression by using a multiple-trait colocalization ("moloc") framework. We found evidence of a possible causal effect for DNA methylation on FEV1 at 18 CpGs (p < 1.2 × 10-4). Replication analysis supported a causal effect at three CpGs (cg21201401 [LIME1 and ZGPAT], cg19758448 [PGAP3], and cg12616487 [EML3 and AHNAK] [p < 0.0028]). DNA methylation did not clearly mediate the effect of smoking on FEV1, although DNA methylation at some sites might influence lung function via effects on smoking. By using "moloc", we found evidence of shared causal variants between lung function, gene expression, and DNA methylation. These findings highlight potential therapeutic targets for improving lung function and possibly smoking cessation, although larger, tissue-specific datasets are required to confirm these results.

Bayesian network analysis incorporating genetic anchors complements conventional Mendelian randomization approaches for exploratory analysis of causal relationships in complex data

Mendelian randomization (MR) implemented through instrumental variables analysis is an increasingly popular causal inference tool used in genetic epidemiology. But it can have limitations for evaluating simultaneous causal relationships in complex data sets that include, for example, multiple genetic predictors and multiple potential risk factors associated with the same genetic variant. Here we use real and simulated data to investigate Bayesian network analysis (BN) with the incorporation of directed arcs, representing genetic anchors, as an alternative approach. A Bayesian network describes the conditional dependencies/independencies of variables using a graphical model (a directed acyclic graph) with an accompanying joint probability. In real data, we found BN could be used to infer simultaneous causal relationships that confirmed the individual causal relationships suggested by bi-directional MR, while allowing for the existence of potential horizontal pleiotropy (that would violate MR assumptions). In simulated data, BN with two directional anchors (mimicking genetic instruments) had greater power for a fixed type 1 error than bi-directional MR, while BN with a single directional anchor performed better than or as well as bi-directional MR. Both BN and MR could be adversely affected by violations of their underlying assumptions (such as genetic confounding due to unmeasured horizontal pleiotropy). BN with no directional anchor generated inference that was no better than by chance, emphasizing the importance of directional anchors in BN (as in MR). Under highly pleiotropic simulated scenarios, BN outperformed both MR (and its recent extensions) and two recently-proposed alternative approaches: a multi-SNP mediation intersection-union test (SMUT) and a latent causal variable (LCV) test. We conclude that BN incorporating genetic anchors is a useful complementary method to conventional MR for exploring causal relationships in complex data sets such as those generated from modern "omics" technologies.

Immune-mediated genetic pathways resulting in pulmonary function impairment increase lung cancer susceptibility

Impaired lung function is often caused by cigarette smoking, making it challenging to disentangle its role in lung cancer susceptibility. Investigation of the shared genetic basis of these phenotypes in the UK Biobank and International Lung Cancer Consortium (29,266 cases, 56,450 controls) shows that lung cancer is genetically correlated with reduced forced expiratory volume in one second (FEV1: rg = 0.098, p = 2.3 × 10-8) and the ratio of FEV1 to forced vital capacity (FEV1/FVC: rg = 0.137, p = 2.0 × 10-12). Mendelian randomization analyses demonstrate that reduced FEV1 increases squamous cell carcinoma risk (odds ratio (OR) = 1.51, 95% confidence intervals: 1.21-1.88), while reduced FEV1/FVC increases the risk of adenocarcinoma (OR = 1.17, 1.01-1.35) and lung cancer in never smokers (OR = 1.56, 1.05-2.30). These findings support a causal role of pulmonary impairment in lung cancer etiology. Integrative analyses reveal that pulmonary function instruments, including 73 novel variants, influence lung tissue gene expression and implicate immune-related pathways in mediating the observed effects on lung carcinogenesis.

Validation and characterisation of a DNA methylation alcohol biomarker across the life course

BACKGROUND

Recently, an alcohol predictor was developed using DNA methylation at 144 CpG sites (DNAm-Alc) as a biomarker for improved clinical or epidemiologic assessment of alcohol-related ill health. We validate the performance and characterise the drivers of this DNAm-Alc for the first time in independent populations.

RESULTS

In N = 1049 parents from the Avon Longitudinal Study of Parents and Children (ALSPAC) Accessible Resource for Integrated Epigenomic Studies (ARIES) at midlife, we found DNAm-Alc explained 7.6% of the variation in alcohol intake, roughly half of what had been reported previously, and interestingly explained a larger 9.8% of Alcohol Use Disorders Identification Test (AUDIT) score, a scale of alcohol use disorder. Explanatory capacity in participants from the offspring generation of ARIES measured during adolescence was much lower. However, DNAm-Alc explained 14.3% of the variation in replication using the Head and Neck 5000 (HN5000) clinical cohort that had higher average alcohol consumption. To investigate whether this relationship was being driven by genetic and/or earlier environment confounding, we examined how earlier versus concurrent DNAm-Alc measures predicted AUDIT scores. In both ARIES parental and offspring generations, we observed associations between AUDIT and concurrent, but not earlier DNAm-Alc, suggesting independence from genetic and stable environmental contributions.

CONCLUSIONS

The stronger relationship between DNAm-Alcs and AUDIT in parents at midlife compared to adolescents despite similar levels of consumption suggests that DNAm-Alc likely reflects long-term patterns of alcohol abuse. Such biomarkers may have potential applications for biomonitoring and risk prediction, especially in cases where reporting bias is a concern.

DNA Methylation Signatures of Breastfeeding in Buccal Cells Collected in Mid-Childhood

Breastfeeding has long-term benefits for children that may be mediated via the epigenome. This pathway has been hypothesized, but the number of empirical studies in humans is small and mostly done by using peripheral blood as the DNA source. We performed an epigenome-wide association study (EWAS) in buccal cells collected around age nine (mean = 9.5) from 1006 twins recruited by the Netherlands Twin Register (NTR). An age-stratified analysis examined if effects attenuate with age (median split at 10 years; n<10 = 517, mean age = 7.9; n>10 = 489, mean age = 11.2). We performed replication analyses in two independent cohorts from the NTR (buccal cells) and the Avon Longitudinal Study of Parents and Children (ALSPAC) (peripheral blood), and we tested loci previously associated with breastfeeding in epigenetic studies. Genome-wide DNA methylation was assessed with the Illumina Infinium MethylationEPIC BeadChip (Illumina, San Diego, CA, USA) in the NTR and with the HumanMethylation450 Bead Chip in the ALSPAC. The duration of breastfeeding was dichotomized ('never' vs. 'ever'). In the total sample, no robustly associated epigenome-wide significant CpGs were identified (α = 6.34 × 10-8). In the sub-group of children younger than 10 years, four significant CpGs were associated with breastfeeding after adjusting for child and maternal characteristics. In children older than 10 years, methylation differences at these CpGs were smaller and non-significant. The findings did not replicate in the NTR sample (n = 98; mean age = 7.5 years), and no nearby sites were associated with breastfeeding in the ALSPAC study (n = 938; mean age = 7.4). Of the CpG sites previously reported in the literature, three were associated with breastfeeding in children younger than 10 years, thus showing that these CpGs are associated with breastfeeding in buccal and blood cells. Our study is the first to show that breastfeeding is associated with epigenetic variation in buccal cells in children. Further studies are needed to investigate if methylation differences at these loci are caused by breastfeeding or by other unmeasured confounders, as well as what mechanism drives changes in associations with age.

Methods for Dealing With Missing Covariate Data in Epigenome-Wide Association Studies

Multiple imputation (MI) is a well-established method for dealing with missing data. MI is computationally intensive when imputing missing covariates with high-dimensional outcome data (e.g., DNA methylation data in epigenome-wide association studies (EWAS)), because every outcome variable must be included in the imputation model to avoid biasing associations towards the null. Instead, EWAS analyses are reduced to only complete cases, limiting statistical power and potentially causing bias. We used simulations to compare 5 MI methods for high-dimensional data under 2 missingness mechanisms. All imputation methods had increased power over complete-case (C-C) analyses. Imputing missing values separately for each variable was computationally inefficient, but dividing sites at random into evenly sized bins improved efficiency and gave low bias. Methods imputing solely using subsets of sites identified by the C-C analysis suffered from bias towards the null. However, if these subsets were added into random bins of sites, this bias was reduced. The optimal methods were applied to an EWAS with missingness in covariates. All methods identified additional sites over the C-C analysis, and many of these sites had been replicated in other studies. These methods are also applicable to other high-dimensional data sets, including the rapidly expanding area of "-omics" studies.

Serum cholesterol, MTHFR methylation, and symptoms of depression in children

Depression is associated with dietary factors and epigenetics. Serum cholesterol, which is prone to dietary influences, has been linked to symptoms of depression. This relationship may be (in part) due to altered epigenetic regulation of Methylenetetrahydrofolate Reductase (MTHFR). MTHFR codes for the MTHFR enzyme, which has diverse metabolic functions, and has recently been linked individually with diet, serum cholesterol levels and depressive symptoms. In 514 mother-child pairs, we examined prospective relationships between maternal (pregnancy) and child (7 years) serum cholesterol, MTHFR DNA methylation (DNAm; birth, 7 years), and development of depression symptoms from 8-15 years. After adjusting for potential confounding, we had three main findings. First, higher prenatal cholesterol associated (at a small effect size) with higher MTHFR DNAm at birth. Second, there was small effect size continuity for MTHFR DNAm between birth and age 7. Third, higher age 7 MTHFR DNAm associated with higher initial symptoms of depression symptoms at age 8, again at a small effect size. Overall, our findings provide preliminary evidence for a relationship between prenatal cholesterol, MTHFR DNAm, and symptoms of depression in children. (PsycINFO Database Record (c) 2019 APA, all rights reserved).

Pathways between early-life adversity and adolescent self-harm: the mediating role of inflammation in the Avon Longitudinal Study of Parents and Children

BACKGROUND

Adverse childhood experiences (ACEs) such as physical and emotional abuse are strongly associated with self-harm, but mechanisms underlying this relationship are unclear. Inflammation has been linked to both the experience of ACEs and self-harm or suicide in prior research. This is the first study to examine whether inflammatory markers mediate the association between exposure to ACEs and self-harm.

METHODS

Participants were 4,308 young people from the Avon Longitudinal Study of Parents and Children (ALSPAC), a population-based birth cohort in the United Kingdom. A structural equation modelling approach was used to fit a mediation model with the number of ACEs experienced between ages 0 and 9 years old (yo), levels of the inflammatory markers interleukin-6 and C-reactive protein measured at 9.5 yo, and self-harm reported at 16 yo.

RESULTS

The mean number of ACEs young people experienced was 1.41 (SE 0.03). Higher ACE scores were associated with an increased risk of self-harm at 16 yo (direct effect relative risk (RR) per additional ACE 1.11, 95% CI 1.05, 1.18, p < 0.001). We did not find evidence of an indirect effect of ACEs on self-harm via inflammation (RR 1.00, 95% CI 1.00, 1.01, p = 0.38).

CONCLUSIONS

Young people who have been exposed to ACEs are a group at high risk of self-harm. The association between ACEs and self-harm does not appear to be mediated by an inflammatory process in childhood, as indexed by peripheral levels of circulating inflammatory markers measured in childhood. Further research is needed to identify alternative psychological and biological mechanisms underlying this relationship.