Epigenome-wide association study of dietary fatty acid intake

BACKGROUND

Dietary intake of n-3 polyunsaturated fatty acids (PUFA) may have a protective effect on the development of cardiovascular diseases, diabetes, depression and cancer, while a high intake of n-6 PUFA was often reported to be associated with inflammation-related traits. The effect of PUFAs on health outcomes might be mediated by DNA methylation (DNAm). The aim of our study is to identify the impact of PUFA intake on DNAm in the Cooperative Health Research in the Region of Augsburg (KORA) FF4 cohort and the Leiden Longevity Study (LLS).

RESULTS

DNA methylation levels were measured in whole blood from the population-based KORA FF4 study (N = 1354) and LLS (N = 448), using the Illumina MethylationEPIC BeadChip and Illumina HumanMethylation450 array, respectively. We assessed associations between DNAm and intake of eight and four PUFAs in KORA and LLS, respectively. Where possible, results were meta-analyzed. Below the Bonferroni correction threshold (p < 7.17 × 10-8), we identified two differentially methylated positions (DMPs) associated with PUFA intake in the KORA study. The DMP cg19937480, annotated to gene PRDX1, was positively associated with docosahexaenoic acid (DHA) in model 1 (beta: 2.00 × 10-5, 95%CI: 1.28 × 10-5-2.73 × 10-5, P value: 6.98 × 10-8), while cg05041783, annotated to gene MARK2, was positively associated with docosapentaenoic acid (DPA) in our fully adjusted model (beta: 9.80 × 10-5, 95%CI: 6.25 × 10-5-1.33 × 10-4, P value: 6.75 × 10-8). In the meta-analysis, we identified the CpG site (cg15951061), annotated to gene CDCA7L below Bonferroni correction (1.23 × 10-7) associated with eicosapentaenoic acid (EPA) intake in model 1 (beta: 2.00 × 10-5, 95% CI: 1.27 × 10-5-2.73 × 10-5, P value = 5.99 × 10-8) and we confirmed the association of cg19937480 with DHA in both models 1 and 2 (beta: 2.07 × 10-5, 95% CI: 1.31 × 10-5-2.83 × 10-5, P value = 1.00 × 10-7 and beta: 2.19 × 10-5, 95% CI: 1.41 × 10-5-2.97 × 10-5, P value = 5.91 × 10-8 respectively).

CONCLUSIONS

Our study identified three CpG sites associated with PUFA intake. The mechanisms of these sites remain largely unexplored, highlighting the novelty of our findings. Further research is essential to understand the links between CpG site methylation and PUFA outcomes.

Maternal Mediterranean diet in pregnancy and newborn DNA methylation: a meta-analysis in the PACE Consortium

Higher adherence to the Mediterranean diet during pregnancy is related to a lower risk of preterm birth and to better offspring cardiometabolic health. DNA methylation may be an underlying biological mechanism. We evaluated whether maternal adherence to the Mediterranean diet was associated with offspring cord blood DNA methylation.We meta-analysed epigenome-wide association studies (EWAS) of maternal adherence to the Mediterranean diet during pregnancy and offspring cord blood DNA methylation in 2802 mother-child pairs from five cohorts. We calculated the relative Mediterranean diet (rMED) score with range 0-18 and an adjusted rMED excluding alcohol (rMEDp, range 0-16). DNA methylation was measured using Illumina 450K arrays. We used robust linear regression modelling adjusted for child sex, maternal education, age, smoking, body mass index, energy intake, batch, and cell types. We performed several functional analyses and examined the persistence of differential DNA methylation into childhood (4.5-7.8 y).rMEDp was associated with cord blood DNA methylation at cg23757341 (0.064% increase in DNA methylation per 1-point increase in the rMEDp score, SE = 0.011, P = 2.41 × 10-8). This cytosine-phosphate-guanine (CpG) site maps to WNT5B, associated with adipogenesis and glycaemic phenotypes. We did not identify associations with childhood gene expression, nor did we find enriched biological pathways. The association did not persist into childhood.In this meta-analysis, maternal adherence to the Mediterranean diet (excluding alcohol) during pregnancy was associated with cord blood DNA methylation level at cg23757341. Potential mediation of DNA methylation in associations with offspring health requires further study.

Sex steroid hormones and risk of breast cancer: a two-sample Mendelian randomization study

BACKGROUND

Breast cancer (BC) has the highest cancer incidence and mortality in women worldwide. Observational epidemiological studies suggest a positive association between testosterone, estradiol, dehydroepiandrosterone sulphate (DHEAS) and other sex steroid hormones with postmenopausal BC. We used a two-sample Mendelian randomization analysis to investigate this association.

METHODS

Genetic instruments for nine sex steroid hormones and sex hormone-binding globulin (SHBG) were obtained from genome-wide association studies (GWAS) of UK Biobank (total testosterone (TT) N: 230,454, bioavailable testosterone (BT) N: 188,507 and SHBG N: 189,473), The United Kingdom Household Longitudinal Study (DHEAS N: 9722), the LIFE-Adult and LIFE-Heart cohorts (estradiol N: 2607, androstenedione N: 711, aldosterone N: 685, progesterone N: 1259 and 17-hydroxyprogesterone N: 711) and the CORtisol NETwork (CORNET) consortium (cortisol N: 25,314). Outcome GWAS summary statistics were obtained from the Breast Cancer Association Consortium (BCAC) for overall BC risk (N: 122,977 cases and 105,974 controls) and subtype-specific analyses.

RESULTS

We found that a standard deviation (SD) increase in TT, BT and estradiol increased the risk of overall BC (OR 1.14, 95% CI 1.09-1.21, OR 1.19, 95% CI 1.07-1.33 and OR 1.03, 95% CI 1.01-1.06, respectively) and ER + BC (OR 1.19, 95% CI 1.12-1.27, OR 1.25, 95% CI 1.11-1.40 and OR 1.06, 95% CI 1.03-1.09, respectively). An SD increase in DHEAS also increased ER + BC risk (OR 1.09, 95% CI 1.03-1.16). Subtype-specific analyses showed similar associations with ER+ expressing subtypes: luminal A-like BC, luminal B-like BC and luminal B/HER2-negative-like BC.

CONCLUSIONS

TT, BT, DHEAS and estradiol increase the risk of ER+ type BCs similar to observational studies. Understanding the role of sex steroid hormones in BC risk, particularly subtype-specific risks, highlights the potential importance of attempts to modify and/or monitor hormone levels in order to prevent BC.

Salicylic Acid and Risk of Colorectal Cancer: A Two-Sample Mendelian Randomization Study

Salicylic acid (SA) has observationally been shown to decrease colorectal cancer (CRC) risk. Aspirin (acetylsalicylic acid, that rapidly deacetylates to SA) is an effective primary and secondary chemopreventive agent. Through a Mendelian randomization (MR) approach, we aimed to address whether levels of SA affected CRC risk, stratifying by aspirin use. A two-sample MR analysis was performed using GWAS summary statistics of SA (INTERVAL and EPIC-Norfolk, N = 14,149) and CRC (CCFR, CORECT, GECCO and UK Biobank, 55,168 cases and 65,160 controls). The DACHS study (4410 cases and 3441 controls) was used for replication and stratification of aspirin-use. SNPs proxying SA were selected via three methods: (1) functional SNPs that influence the activity of aspirin-metabolising enzymes; (2) pathway SNPs present in enzymes' coding regions; and (3) genome-wide significant SNPs. We found no association between functional SNPs and SA levels. The pathway and genome-wide SNPs showed no association between SA and CRC risk (OR: 1.03, 95% CI: 0.84-1.27 and OR: 1.08, 95% CI: 0.86-1.34, respectively). Results remained unchanged upon aspirin use stratification. We found little evidence to suggest that an SD increase in genetically predicted SA protects against CRC risk in the general population and upon stratification by aspirin use.

A Combined Proteomics and Mendelian Randomization Approach to Investigate the Effects of Aspirin-Targeted Proteins on Colorectal Cancer

BACKGROUND

Evidence for aspirin's chemopreventative properties on colorectal cancer (CRC) is substantial, but its mechanism of action is not well-understood. We combined a proteomic approach with Mendelian randomization (MR) to identify possible new aspirin targets that decrease CRC risk.

METHODS

Human colorectal adenoma cells (RG/C2) were treated with aspirin (24 hours) and a stable isotope labeling with amino acids in cell culture (SILAC) based proteomics approach identified altered protein expression. Protein quantitative trait loci (pQTLs) from INTERVAL (N = 3,301) and expression QTLs (eQTLs) from the eQTLGen Consortium (N = 31,684) were used as genetic proxies for protein and mRNA expression levels. Two-sample MR of mRNA/protein expression on CRC risk was performed using eQTL/pQTL data combined with CRC genetic summary data from the Colon Cancer Family Registry (CCFR), Colorectal Transdisciplinary (CORECT), Genetics and Epidemiology of Colorectal Cancer (GECCO) consortia and UK Biobank (55,168 cases and 65,160 controls).

RESULTS

Altered expression was detected for 125/5886 proteins. Of these, aspirin decreased MCM6, RRM2, and ARFIP2 expression, and MR analysis showed that a standard deviation increase in mRNA/protein expression was associated with increased CRC risk (OR: 1.08, 95% CI, 1.03-1.13; OR: 3.33, 95% CI, 2.46-4.50; and OR: 1.15, 95% CI, 1.02-1.29, respectively).

CONCLUSIONS

MCM6 and RRM2 are involved in DNA repair whereby reduced expression may lead to increased DNA aberrations and ultimately cancer cell death, whereas ARFIP2 is involved in actin cytoskeletal regulation, indicating a possible role in aspirin's reduction of metastasis.

IMPACT

Our approach has shown how laboratory experiments and population-based approaches can combine to identify aspirin-targeted proteins possibly affecting CRC risk.

Discovering protein drug targets using knowledge graph embeddings

MOTIVATION

Computational approaches for predicting drug-target interactions (DTIs) can provide valuable insights into the drug mechanism of action. DTI predictions can help to quickly identify new promising (on-target) or unintended (off-target) effects of drugs. However, existing models face several challenges. Many can only process a limited number of drugs and/or have poor proteome coverage. The current approaches also often suffer from high false positive prediction rates.

RESULTS

We propose a novel computational approach for predicting drug target proteins. The approach is based on formulating the problem as a link prediction in knowledge graphs (robust, machine-readable representations of networked knowledge). We use biomedical knowledge bases to create a knowledge graph of entities connected to both drugs and their potential targets. We propose a specific knowledge graph embedding model, TriModel, to learn vector representations (i.e. embeddings) for all drugs and targets in the created knowledge graph. These representations are consequently used to infer candidate drug target interactions based on their scores computed by the trained TriModel model. We have experimentally evaluated our method using computer simulations and compared it to five existing models. This has shown that our approach outperforms all previous ones in terms of both area under ROC and precision-recall curves in standard benchmark tests.

AVAILABILITY AND IMPLEMENTATION

The data, predictions and models are available at: drugtargets.insight-centre.org.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Predicting The Effects of Chemical-Protein Interactions On Proteins Using Tensor Factorisation

Understanding the different effects of chemical substances on human proteins is fundamental for designing new drugs. It is also important for elucidating the different mechanisms of action of drugs that can cause side-effects. In this context, computational methods for predicting chemical-protein interactions can provide valuable insights on the relation between therapeutic chemical substances and proteins. Their predictions therefore can help in multiple tasks such as drug repurposing, identifying new drug side-effects, etc. Despite their useful predictions, these methods are unable to predict the different implications - such as change in protein expression, abundance, etc, - of chemical - protein interactions. Therefore, In this work, we study the modelling of chemical-protein interactions' effects on proteins activity using computational approaches. We hereby propose using 3D tensors to model chemicals, their target proteins and the effects associated to their interactions. We then use multi-part embedding tensor factorisation to predict the different effects of chemicals on human proteins. We assess the predictive accuracy of our proposed method using a benchmark dataset that we built. We then show by computational experimental evaluation that our approach outperforms other tensor factorisation methods in the task of predicting effects of chemicals on human proteins.

Biological applications of knowledge graph embedding models

Complex biological systems are traditionally modelled as graphs of interconnected biological entities. These graphs, i.e. biological knowledge graphs, are then processed using graph exploratory approaches to perform different types of analytical and predictive tasks. Despite the high predictive accuracy of these approaches, they have limited scalability due to their dependency on time-consuming path exploratory procedures. In recent years, owing to the rapid advances of computational technologies, new approaches for modelling graphs and mining them with high accuracy and scalability have emerged. These approaches, i.e. knowledge graph embedding (KGE) models, operate by learning low-rank vector representations of graph nodes and edges that preserve the graph's inherent structure. These approaches were used to analyse knowledge graphs from different domains where they showed superior performance and accuracy compared to previous graph exploratory approaches. In this work, we study this class of models in the context of biological knowledge graphs and their different applications. We then show how KGE models can be a natural fit for representing complex biological knowledge modelled as graphs. We also discuss their predictive and analytical capabilities in different biology applications. In this regard, we present two example case studies that demonstrate the capabilities of KGE models: prediction of drug-target interactions and polypharmacy side effects. Finally, we analyse different practical considerations for KGEs, and we discuss possible opportunities and challenges related to adopting them for modelling biological systems.

Association Between Telomere Length and Risk of Cancer and Non-Neoplastic Diseases: A Mendelian Randomization Study

IMPORTANCE

The causal direction and magnitude of the association between telomere length and incidence of cancer and non-neoplastic diseases is uncertain owing to the susceptibility of observational studies to confounding and reverse causation.

OBJECTIVE

To conduct a Mendelian randomization study, using germline genetic variants as instrumental variables, to appraise the causal relevance of telomere length for risk of cancer and non-neoplastic diseases.

DATA SOURCES

Genomewide association studies (GWAS) published up to January 15, 2015.

STUDY SELECTION

GWAS of noncommunicable diseases that assayed germline genetic variation and did not select cohort or control participants on the basis of preexisting diseases. Of 163 GWAS of noncommunicable diseases identified, summary data from 103 were available.

DATA EXTRACTION AND SYNTHESIS

Summary association statistics for single nucleotide polymorphisms (SNPs) that are strongly associated with telomere length in the general population.

MAIN OUTCOMES AND MEASURES

Odds ratios (ORs) and 95% confidence intervals (CIs) for disease per standard deviation (SD) higher telomere length due to germline genetic variation.

RESULTS

Summary data were available for 35 cancers and 48 non-neoplastic diseases, corresponding to 420 081 cases (median cases, 2526 per disease) and 1 093 105 controls (median, 6789 per disease). Increased telomere length due to germline genetic variation was generally associated with increased risk for site-specific cancers. The strongest associations (ORs [95% CIs] per 1-SD change in genetically increased telomere length) were observed for glioma, 5.27 (3.15-8.81); serous low-malignant-potential ovarian cancer, 4.35 (2.39-7.94); lung adenocarcinoma, 3.19 (2.40-4.22); neuroblastoma, 2.98 (1.92-4.62); bladder cancer, 2.19 (1.32-3.66); melanoma, 1.87 (1.55-2.26); testicular cancer, 1.76 (1.02-3.04); kidney cancer, 1.55 (1.08-2.23); and endometrial cancer, 1.31 (1.07-1.61). Associations were stronger for rarer cancers and at tissue sites with lower rates of stem cell division. There was generally little evidence of association between genetically increased telomere length and risk of psychiatric, autoimmune, inflammatory, diabetic, and other non-neoplastic diseases, except for coronary heart disease (OR, 0.78 [95% CI, 0.67-0.90]), abdominal aortic aneurysm (OR, 0.63 [95% CI, 0.49-0.81]), celiac disease (OR, 0.42 [95% CI, 0.28-0.61]) and interstitial lung disease (OR, 0.09 [95% CI, 0.05-0.15]).

CONCLUSIONS AND RELEVANCE

It is likely that longer telomeres increase risk for several cancers but reduce risk for some non-neoplastic diseases, including cardiovascular diseases.