Genetics of circulating inflammatory proteins identifies drivers of immune-mediated disease risk and therapeutic targets.

Circulating proteins have important functions in inflammation and a broad range of diseases. To identify genetic influences on inflammation-related proteins, we conducted a genome-wide protein quantitative trait locus (pQTL) study of 91 plasma proteins measured using the Olink Target platform in 14,824 participants. We identified 180 pQTLs (59 cis, 121 trans). Integration of pQTL data with eQTL and disease genome-wide association studies provided insight into pathogenesis, implicating lymphotoxin-α in multiple sclerosis. Using Mendelian randomization (MR) to assess causality in disease etiology, we identified both shared and distinct effects of specific proteins across immune-mediated diseases, including directionally discordant effects of CD40 on risk of rheumatoid arthritis versus multiple sclerosis and inflammatory bowel disease. MR implicated CXCL5 in the etiology of ulcerative colitis (UC) and we show elevated gut CXCL5 transcript expression in patients with UC. These results identify targets of existing drugs and provide a powerful resource to facilitate future drug target prioritization.

An atlas of genetic scores to predict multi-omic traits.

The use of omic modalities to dissect the molecular underpinnings of common diseases and traits is becoming increasingly common. But multi-omic traits can be genetically predicted, which enables highly cost-effective and powerful analyses for studies that do not have multi-omics1. Here we examine a large cohort (the INTERVAL study2; n = 50,000 participants) with extensive multi-omic data for plasma proteomics (SomaScan, n = 3,175; Olink, n = 4,822), plasma metabolomics (Metabolon HD4, n = 8,153), serum metabolomics (Nightingale, n = 37,359) and whole-blood Illumina RNA sequencing (n = 4,136), and use machine learning to train genetic scores for 17,227 molecular traits, including 10,521 that reach Bonferroni-adjusted significance. We evaluate the performance of genetic scores through external validation across cohorts of individuals of European, Asian and African American ancestries. In addition, we show the utility of these multi-omic genetic scores by quantifying the genetic control of biological pathways and by generating a synthetic multi-omic dataset of the UK Biobank3 to identify disease associations using a phenome-wide scan. We highlight a series of biological insights with regard to genetic mechanisms in metabolism and canonical pathway associations with disease; for example, JAK-STAT signalling and coronary atherosclerosis. Finally, we develop a portal ( https://www.omicspred.org/ ) to facilitate public access to all genetic scores and validation results, as well as to serve as a platform for future extensions and enhancements of multi-omic genetic scores.

The risk of venous thromboembolism in oral contraceptive users: the role of genetic factors-a prospective cohort study of 240,000 women in the UK Biobank.

BACKGROUND: More than 150 million women worldwide use oral contraceptives. Women with inherited thrombophilia and carriers of certain thrombophilia gene variants, such as factor V Leiden and the prothrombin, are at an increased risk for venous thromboembolism, especially when combined with oral contraceptive use. Venous thromboembolism is a complex disorder involving many genetic risk factors, and recently, polygenic risk scores have been proposed to capture a significant proportion of the genetic risk of venous thromboembolism. OBJECTIVE: The aim of this study was to estimate the risk for developing venous thromboembolism when initiating oral contraceptive use (first 2 years) and during continued use among women with a high genetic liability. STUDY DESIGN: We used a prospective study design in which 244,420 participants from the UK Biobank were followed from birth. The effect of oral contraceptive use during the first 2 years and in the remaining years of oral contraceptive use on the risk of developing venous thromboembolism was estimated using a Cox regression with a time-dependent exposure variable. Women were stratified according to their polygenic risk scores and whether they were carriers of factor V Leiden and/or prothrombin variants. RESULTS: When genetic risk was not considered, an increased risk for venous thromboembolism was observed during the first 2 years of oral contraceptive use (hazard ratio, 3.09; 95% confidence interval, 3.00-3.20) but not during continued use (hazard ratio, 0.92; 95% confidence interval, 0.80-1.05). However, when genetic risk was considered, women in the highest polygenic risk score category had a more pronounced risk of developing a venous thromboembolism during the first 2 years of oral contraceptive use (hazard ratio, 6.35; 95% confidence interval, 4.98-8.09), and a high risk was also observed among factor V Leiden (hazard ratio, 5.73; 95% confidence interval, 5.31-6.17) and prothrombin variant carriers (hazard ratio, 5.23; 95% confidence interval, 4.67 - 5.87). A high polygenic risk score in combination with being a factor V Leiden and prothrombin variant carrier conferred the highest risk for developing a venous thromboembolism during the first 2 years of oral contraceptive use (hazard ratio, 14.8; 95% confidence interval, 9.28-23.6). Women with a high genetic liability also had an increased risk during continued use but it was less pronounced, and the highest risk was conferred to carriers of both factor V Leiden and the prothrombin variant (hazard ratio, 4.93; 95% confidence interval, 3.16-7.7). CONCLUSION: Evaluating polygenic risk can identify additional venous thromboembolism risk that is not captured in the commonly investigated genes for inherited thrombophilia. Our results indicate that oral contraceptive use is associated with an increased risk for developing a venous thromboembolism, particularly among women with a high genetic predisposition, and that oral contraceptive use dramatically increases the risk thereof short after initiation of use, which decreases with continued use. This suggests that the polygenic risk score could be used to identify women who are at high risk for developing a venous thromboembolism and advise them on alternative methods of contraception.

Decreased neutrophil function in newly diagnosed multiple myeloma patients is restored with lenalidomide therapy.

OBJECTIVES: Bacterial infections are common and a major cause of morbidity and mortality in multiple myeloma (MM). We have investigated the function of polymorphonuclear leukocyte (PMN), the immune system's first line of defense against bacteria, in peripheral blood (PB) and bone marrow (BM) samples from patients with newly diagnosed MM (NDMM), smoldering MM (SMM), monoclonal gammopathy of undetermined significance (MGUS) and healthy controls. METHODS: Phagocytosis and oxidative burst in PMN cells from patients and healthy donors were investigated using PhagoTest and PhagoBurst assay. RESULTS: PMN from NDMM, SMM, and MGUS patients had reduced phagocytosis and oxidative burst ability compared with healthy controls. The dysfunction was most prominent in BM samples from MM, SMM, and MGUS patients. Importantly the reduced phagocytosis in MM patients was restored in patients on lenalidomide therapy. Consistently the ability of Escherichia coli stimulated oxidative burst in BM was reduced for the MM, SMM, and MGUS cohort in contrast to the healthy controls and the patients on lenalidomide treatment. CONCLUSION: Our results show that MM patients have neutrophil dysfunction that could contribute to susceptibility for bacterial infections and that lenalidomide therapy was associated with restored PMN function.

Genetic Landscape of the ACE2 Coronavirus Receptor.

BACKGROUND: SARS-CoV-2, the causal agent of COVID-19, enters human cells using the ACE2 (angiotensin-converting enzyme 2) protein as a receptor. ACE2 is thus key to the infection and treatment of the coronavirus. ACE2 is highly expressed in the heart and respiratory and gastrointestinal tracts, playing important regulatory roles in the cardiovascular and other biological systems. However, the genetic basis of the ACE2 protein levels is not well understood. METHODS: We have conducted the largest genome-wide association meta-analysis of plasma ACE2 levels in >28 000 individuals of the SCALLOP Consortium (Systematic and Combined Analysis of Olink Proteins). We summarize the cross-sectional epidemiological correlates of circulating ACE2. Using the summary statistics-based high-definition likelihood method, we estimate relevant genetic correlations with cardiometabolic phenotypes, COVID-19, and other human complex traits and diseases. We perform causal inference of soluble ACE2 on vascular disease outcomes and COVID-19 severity using mendelian randomization. We also perform in silico functional analysis by integrating with other types of omics data. RESULTS: We identified 10 loci, including 8 novel, capturing 30% of the heritability of the protein. We detected that plasma ACE2 was genetically correlated with vascular diseases, severe COVID-19, and a wide range of human complex diseases and medications. An X-chromosome cis-protein quantitative trait loci-based mendelian randomization analysis suggested a causal effect of elevated ACE2 levels on COVID-19 severity (odds ratio, 1.63 [95% CI, 1.10-2.42]; P=0.01), hospitalization (odds ratio, 1.52 [95% CI, 1.05-2.21]; P=0.03), and infection (odds ratio, 1.60 [95% CI, 1.08-2.37]; P=0.02). Tissue- and cell type-specific transcriptomic and epigenomic analysis revealed that the ACE2 regulatory variants were enriched for DNA methylation sites in blood immune cells. CONCLUSIONS: Human plasma ACE2 shares a genetic basis with cardiovascular disease, COVID-19, and other related diseases. The genetic architecture of the ACE2 protein is mapped, providing a useful resource for further biological and clinical studies on this coronavirus receptor.

Precision medicine in complex diseases-Molecular subgrouping for improved prediction and treatment stratification.

Complex diseases are caused by a combination of genetic, lifestyle, and environmental factors and comprise common noncommunicable diseases, including allergies, cardiovascular disease, and psychiatric and metabolic disorders. More than 25% of Europeans suffer from a complex disease, and together these diseases account for 70% of all deaths. The use of genomic, molecular, or imaging data to develop accurate diagnostic tools for treatment recommendations and preventive strategies, and for disease prognosis and prediction, is an important step toward precision medicine. However, for complex diseases, precision medicine is associated with several challenges. There is a significant heterogeneity between patients of a specific disease-both with regards to symptoms and underlying causal mechanisms-and the number of underlying genetic and nongenetic risk factors is often high. Here, we summarize precision medicine approaches for complex diseases and highlight the current breakthroughs as well as the challenges. We conclude that genomic-based precision medicine has been used mainly for patients with highly penetrant monogenic disease forms, such as cardiomyopathies. However, for most complex diseases-including psychiatric disorders and allergies-available polygenic risk scores are more probabilistic than deterministic and have not yet been validated for clinical utility. However, subclassifying patients of a specific disease into discrete homogenous subtypes based on molecular or phenotypic data is a promising strategy for improving diagnosis, prediction, treatment, prevention, and prognosis. The availability of high-throughput molecular technologies, together with large collections of health data and novel data-driven approaches, offers promise toward improved individual health through precision medicine.

Effects of oral contraceptives and menopausal hormone therapy on the risk of rheumatoid arthritis: a prospective cohort study.

OBJECTIVES: Oral contraceptives (OC) and menopausal hormone therapy (MHT) contain exogenous sex hormones and are used by millions of women around the world. However, their effect on development of rheumatoid arthritis (RA) is still debated and the current literature suggests that they may exert opposite effects on the risk of RA. The present study aimed to estimate the effects of exogenous hormones on development of RA, both during the reproductive lifespan and later in life. METHODS: The association between OC and RA, as well as between MHT and late-onset RA (LORA), was investigated using time-dependent Cox regression modelling in white British women from the UK Biobank (N = 236 602 and N = 102 466, respectively) and replicated in women from all ethnic groups. RESULTS: OC use was associated with a decreased risk of RA in ever-users (hazard ratio [HR]=0.89; 95% CI = 0.82-0.96), as well as in current (HR = 0.81; 0.73-0.91) and former users (HR = 0.92; 0.84 -1.00), compared with never-users. In contrast, MHT use was associated with an increased risk of LORA in ever-users (HR = 1.16; 1.06-1.26) as well as in former users (HR = 1.13; 1.03-1.24) compared with never-users. CONCLUSION: OC use appears to protect against RA, while MHT may increase the risk of LORA. This study provides new insights into the possible inverse effect of exposure to different exogenous sex hormones on the risk of RA.

Oral Contraceptives, Hormone Replacement Therapy, and Stroke Risk.

BACKGROUND: Millions of women worldwide use exogenous hormones as oral contraceptives or hormone replacement therapy. Still, time-dependent and long-term consequences of exogenous hormones on stroke risk remains unclear. METHODS: We examined the association between self-reported oral contraceptive and hormone replacement therapy use and stroke risk in 257 194 women from the UK Biobank, born between 1939 and 1970. Outcomes included any type of stroke, ischemic stroke, intracerebral hemorrhage, and subarachnoid hemorrhage. Exposures were analyzed as time-varying variables in Cox regression models. RESULTS: During first year of oral contraceptive use, an increased event rate of any stroke was observed (hazard ratio [HR], 2.49 [95% CI, 1.44-4.30]), while the hazards were found to be comparable during remaining years of use (HR, 1.00 [95% CI, 0.86-1.14]), compared with nonusers. Similarly, first year of hormone replacement therapy use was associated with higher hazard rates of any stroke (HR, 2.12 [95% CI, 1.66-2.70]), as well as cause-specific stroke, including ischemic stroke (HR, 1.93 [95% CI, 1.05-3.57]) and subarachnoid hemorrhage (HR, 2.17 [95% CI, 1.25-3.78]), which remained increased for any stroke during remaining years of use (HR, 1.18 [95% CI, 1.05-1.31]), and after discontinuation (HR, 1.16 [95% CI, 1.02-1.32]). CONCLUSIONS: Oral contraceptive use and hormone replacement therapy were associated with an increased risk of stroke, especially during the first year of use, possibly due to immediate changes in hemostatic balance. This study provides new insights on the effects of hormone exposure on stroke risk and provide evidence of not only an overall risk but also a pronounced effects seen in the beginning of treatment.

Immune cells lacking Y chromosome show dysregulation of autosomal gene expression.

Epidemiological investigations show that mosaic loss of chromosome Y (LOY) in leukocytes is associated with earlier mortality and morbidity from many diseases in men. LOY is the most common acquired mutation and is associated with aberrant clonal expansion of cells, yet it remains unclear whether this mosaicism exerts a direct physiological effect. We studied DNA and RNA from leukocytes in sorted- and single-cells in vivo and in vitro. DNA analyses of sorted cells showed that men diagnosed with Alzheimer's disease was primarily affected with LOY in NK cells whereas prostate cancer patients more frequently displayed LOY in CD4 + T cells and granulocytes. Moreover, bulk and single-cell RNA sequencing in leukocytes allowed scoring of LOY from mRNA data and confirmed considerable variation in the rate of LOY across individuals and cell types. LOY-associated transcriptional effect (LATE) was observed in ~ 500 autosomal genes showing dysregulation in leukocytes with LOY. The fraction of LATE genes within specific cell types was substantially larger than the fraction of LATE genes shared between different subsets of leukocytes, suggesting that LOY might have pleiotropic effects. LATE genes are involved in immune functions but also encode proteins with roles in other diverse biological processes. Our findings highlight a surprisingly broad role for chromosome Y, challenging the view of it as a "genetic wasteland", and support the hypothesis that altered immune function in leukocytes could be a mechanism linking LOY to increased risk for disease.

Association between Brachyspira and irritable bowel syndrome with diarrhoea.

OBJECTIVE: The incidence of IBS increases following enteric infections, suggesting a causative role for microbial imbalance. However, analyses of faecal microbiota have not demonstrated consistent alterations. Here, we used metaproteomics to investigate potential associations between mucus-resident microbiota and IBS symptoms. DESIGN: Mucus samples were prospectively collected from sigmoid colon biopsies from patients with IBS and healthy volunteers, and their microbial protein composition analysed by mass spectrometry. Observations were verified by immunofluorescence, electron microscopy and real-time PCR, further confirmed in a second cohort, and correlated with comprehensive profiling of clinical characteristics and mucosal immune responses. RESULTS: Metaproteomic analysis of colon mucus samples identified peptides from potentially pathogenic Brachyspira species in a subset of patients with IBS. Using multiple diagnostic methods, mucosal Brachyspira colonisation was detected in a total of 19/62 (31%) patients with IBS from two prospective cohorts, versus 0/31 healthy volunteers (p<0.001). The prevalence of Brachyspira colonisation in IBS with diarrhoea (IBS-D) was 40% in both cohorts (p=0.02 and p=0.006 vs controls). Brachyspira attachment to the colonocyte apical membrane was observed in 20% of patients with IBS and associated with accelerated oro-anal transit, mild mucosal inflammation, mast cell activation and alterations of molecular pathways linked to bacterial uptake and ion-fluid homeostasis. Metronidazole treatment paradoxically promoted Brachyspira relocation into goblet cell secretory granules-possibly representing a novel bacterial strategy to evade antibiotics. CONCLUSION: Mucosal Brachyspira colonisation was significantly more common in IBS and associated with distinctive clinical, histological and molecular characteristics. Our observations suggest a role for Brachyspira in the pathogenesis of IBS, particularly IBS-D.

Genome-wide association analysis of 350 000 Caucasians from the UK Biobank identifies novel loci for asthma, hay fever and eczema.

Even though heritability estimates suggest that the risk of asthma, hay fever and eczema is largely due to genetic factors, previous studies have not explained a large part of the genetics behind these diseases. In this genome-wide association study, we include 346 545 Caucasians from the UK Biobank to identify novel loci for asthma, hay fever and eczema and replicate novel loci in three independent cohorts. We further investigate if associated lead single nucleotide polymorphisms (SNPs) have a significantly larger effect for one disease compared to the other diseases, to highlight possible disease-specific effects. We identified 141 loci, of which 41 are novel, to be associated (P ≤ 3 × 10-8) with asthma, hay fever or eczema, analyzed separately or as disease phenotypes that includes the presence of different combinations of these diseases. The largest number of loci was associated with the combined phenotype (asthma/hay fever/eczema). However, as many as 20 loci had a significantly larger effect on hay fever/eczema only compared to their effects on asthma, while 26 loci exhibited larger effects on asthma compared with their effects on hay fever/eczema. At four of the novel loci, TNFRSF8, MYRF, TSPAN8, and BHMG1, the lead SNPs were in Linkage Disequilibrium (LD) (>0.8) with potentially casual missense variants. Our study shows that a large amount of the genetic contribution is shared between the diseases. Nonetheless, a number of SNPs have a significantly larger effect on one of the phenotypes, suggesting that part of the genetic contribution is more phenotype specific.

A combined genome-wide association and molecular study of age-related hearing loss in H. sapiens.

BACKGROUND: Sensorineural hearing loss is one of the most common sensory deficiencies. However, the molecular contribution to age-related hearing loss is not fully elucidated. METHODS: We performed genome-wide association studies (GWAS) for hearing loss-related traits in the UK Biobank (N = 362,396) and selected a high confidence set of ten hearing-associated gene products for staining in human cochlear samples: EYA4, LMX1A, PTK2/FAK, UBE3B, MMP2, SYNJ2, GRM5, TRIOBP, LMO-7, and NOX4. RESULTS: All proteins were found to be expressed in human cochlear structures. Our findings illustrate cochlear structures that mediate mechano-electric transduction of auditory stimuli, neuronal conductance, and neuronal plasticity to be involved in age-related hearing loss. CONCLUSIONS: Our results suggest common genetic variation to influence structural resilience to damage as well as cochlear recovery after trauma, which protect against accumulated damage to cochlear structures and the development of hearing loss over time.

Characterization of the human ABO genotypes and their association to common inflammatory and cardiovascular diseases in the UK Biobank.

The ABO gene contains three major alleles that encodes different antigens; A, B, and O, which determine an individual's blood group. Previous studies have primarily focused on identifying associations between ABO blood groups and diseases risk. Here, we sought to test for association between ABO genotypes (OO, OA, AA; OB, BB, and AB) and a large set of common inflammatory and cardiovascular diseases in UK Biobank as well as disease-related protein biomarkers in NSPHS. We first tested for association by conducting a likelihood ratio test, testing whether ABO contributed significantly to the risk for 24 diseases, and 438 plasma proteins. For phenotypes with FDR < 0.05, we tested for pair-wise differences between genetically determined ABO genotypes using logistic or linear regression. Our study confirmed previous findings of a strong association between ABO and cardiovascular disease, identified associations for both type 1 and type 2 diabetes, and provide additional evidence of significant differences between heterozygous and homozygous allele carriers for pulmonary embolism, deep vein thrombosis, but also for von Willebrand factor levels. Furthermore, the results indicated an additive effect between genotypes, even between the two most common A subgroups, A1 and A2. Additionally, we found that ABO contributed significantly to 39 plasma proteins, of which 23 have never been linked to the ABO locus before. These results show the need of incorporating ABO genotype information in the consultation and management of patients at risk, rather than classifying patients into blood groups.

Geographical differences in tonsillar carriage rates of Fusobacterium necrophorum - A cross-sectional study in Sweden and Zambia.

While Fusobacterium necrophorum historically has been considered normal tonsillar flora, recent studies from Europe and the US have suggested that carriage occur transiently in adolescence and young adulthood. However, no studies originating from Africa exist. In this cross-sectional study of tonsillar carriage of F. necrophorum, we aimed to investigate geographical differences in tonsillar carriage rates of F. necrophorum in healthy participants aged 15-25 years in Sweden and Zambia and further investigate the age distribution of tonsillar carriage in Zambia. Specimens were obtained by tonsillar swabs and analyzed with real-time PCR for F. necrophorum. In participants aged 15-25 years, tonsillar carriage was more common in Sweden 21/100 (21%) than in Zambia 6/192 (3%), p < 0.001. In Zambian participants aged above 25 years tonsillar carriage was rare 1/76 (1%). In conclusion, the high rate of tonsillar carriage in participants aged 15-25 years in Sweden has implications on the interpretation of tonsillar findings in patients with pharyngotonsillitis. Interestingly, a geographical difference was found with tonsillar carriage rarely identified in Zambia.

Genetic variants influencing phenotypic variance heterogeneity.

Most genetic studies identify genetic variants associated with disease risk or with the mean value of a quantitative trait. More rarely, genetic variants associated with variance heterogeneity are considered. In this study, we have identified such variance single-nucleotide polymorphisms (vSNPs) and examined if these represent biological gene × gene or gene × environment interactions or statistical artifacts caused by multiple linked genetic variants influencing the same phenotype. We have performed a genome-wide study, to identify vSNPs associated with variance heterogeneity in DNA methylation levels. Genotype data from over 10 million single-nucleotide polymorphisms (SNPs), and DNA methylation levels at over 430 000 CpG sites, were analyzed in 729 individuals. We identified vSNPs for 7195 CpG sites (P < 9.4 × 10-11). This is a relatively low number compared to 52 335 CpG sites for which SNPs were associated with mean DNA methylation levels. We further showed that variance heterogeneity between genotypes mainly represents additional, often rare, SNPs in linkage disequilibrium (LD) with the respective vSNP and for some vSNPs, multiple low frequency variants co-segregating with one of the vSNP alleles. Therefore, our results suggest that variance heterogeneity of DNA methylation mainly represents phenotypic effects by multiple SNPs, rather than biological interactions. Such effects may also be important for interpreting variance heterogeneity of more complex clinical phenotypes.

NCF1-339 polymorphism is associated with altered formation of neutrophil extracellular traps, high serum interferon activity and antiphospholipid syndrome in systemic lupus erythematosus.

​OBJECTIVES: A single nucleotide polymorphism in the NCF1 gene (NCF1-339, rs201802880), encoding NADPH oxidase type II subunit NCF1/p47phox, reducing production of reactive oxygen species (ROS) is strongly associated with the development of systemic lupus erythematosus (SLE). This study aimed at characterising NCF1-339 effects on neutrophil extracellular trap (NET) formation, type I interferon activity and antibody profile in patients with SLE. ​METHODS: Neutrophil NET-release pathways (n=31), serum interferon (n=141) and finally antibody profiles (n=305) were investigated in SLE subjects from Lund, genotyped for NCF1-339. Then, 1087 SLE subjects from the rheumatology departments of four Swedish SLE centres, genotyped for NCF1-339, were clinically characterised to validate these findings. ​RESULTS: Compared with patients with normal-ROS NCF1-339 genotypes, neutrophils from patients with SLE with low-ROS NCF1-339 genotypes displayed impaired NET formation (p<0.01) and increased dependence on mitochondrial ROS (p<0.05). Low-ROS patients also had increased frequency of high serum interferon activity (80% vs 21.4%, p<0.05) and positivity for anti-ß2 glycoprotein I (p<0.01) and anticardiolipin antibodies (p<0.05) but were not associated with other antibodies. We confirmed an over-representation of having any antiphospholipid antibody, OR 1.40 (95% CI 1.01 to 1.95), anti-ß2 glycoprotein I, OR 1.82 (95% CI 1.02 to 3.24) and the antiphospholipid syndrome (APS), OR 1.74 (95% CI 1.19 to 2.55) in all four cohorts (n=1087). ​CONCLUSIONS: The NCF1-339 SNP mediated decreased NADPH oxidase function, is associated with high interferon activity and impaired formation of NETs in SLE, allowing dependence on mitochondrial ROS. Unexpectedly, we revealed a striking connection between the ROS deficient NCF1-339 genotypes and the presence of phospholipid antibodies and APS.

Gene-environment interaction study for BMI reveals interactions between genetic factors and physical activity, alcohol consumption and socioeconomic status.

Previous genome-wide association studies (GWAS) have identified hundreds of genetic loci to be associated with body mass index (BMI) and risk of obesity. Genetic effects can differ between individuals depending on lifestyle or environmental factors due to gene-environment interactions. In this study, we examine gene-environment interactions in 362,496 unrelated participants with Caucasian ancestry from the UK Biobank resource. A total of 94 BMI-associated SNPs, selected from a previous GWAS on BMI, were used to construct weighted genetic scores for BMI (GSBMI). Linear regression modeling was used to estimate the effect of gene-environment interactions on BMI for 131 lifestyle factors related to: dietary habits, smoking and alcohol consumption, physical activity, socioeconomic status, mental health, sleeping patterns, as well as female-specific factors such as menopause and childbirth. In total, 15 lifestyle factors were observed to interact with GSBMI, of which alcohol intake frequency, usual walking pace, and Townsend deprivation index, a measure of socioeconomic status, were all highly significant (p = 1.45*10-29, p = 3.83*10-26, p = 4.66*10-11, respectively). Interestingly, the frequency of alcohol consumption, rather than the total weekly amount resulted in a significant interaction. The FTO locus was the strongest single locus interacting with any of the lifestyle factors. However, 13 significant interactions were also observed after omitting the FTO locus from the genetic score. Our analyses indicate that many lifestyle factors modify the genetic effects on BMI with some groups of individuals having more than double the effect of the genetic score. However, the underlying causal mechanisms of gene-environmental interactions are difficult to deduce from cross-sectional data alone and controlled experiments are required to fully characterise the causal factors.

The relative contribution of DNA methylation and genetic variants on protein biomarkers for human diseases.

Associations between epigenetic alterations and disease status have been identified for many diseases. However, there is no strong evidence that epigenetic alterations are directly causal for disease pathogenesis. In this study, we combined SNP and DNA methylation data with measurements of protein biomarkers for cancer, inflammation or cardiovascular disease, to investigate the relative contribution of genetic and epigenetic variation on biomarker levels. A total of 121 protein biomarkers were measured and analyzed in relation to DNA methylation at 470,000 genomic positions and to over 10 million SNPs. We performed epigenome-wide association study (EWAS) and genome-wide association study (GWAS) analyses, and integrated biomarker, DNA methylation and SNP data using between 698 and 1033 samples depending on data availability for the different analyses. We identified 124 and 45 loci (Bonferroni adjusted P < 0.05) with effect sizes up to 0.22 standard units' change per 1% change in DNA methylation levels and up to four standard units' change per copy of the effective allele in the EWAS and GWAS respectively. Most GWAS loci were cis-regulatory whereas most EWAS loci were located in trans. Eleven EWAS loci were associated with multiple biomarkers, including one in NLRC5 associated with CXCL11, CXCL9, IL-12, and IL-18 levels. All EWAS signals that overlapped with a GWAS locus were driven by underlying genetic variants and three EWAS signals were confounded by smoking. While some cis-regulatory SNPs for biomarkers appeared to have an effect also on DNA methylation levels, cis-regulatory SNPs for DNA methylation were not observed to affect biomarker levels. We present associations between protein biomarker and DNA methylation levels at numerous loci in the genome. The associations are likely to reflect the underlying pattern of genetic variants, specific environmental exposures, or represent secondary effects to the pathogenesis of disease.

Mapping of 79 loci for 83 plasma protein biomarkers in cardiovascular disease.

Recent advances in highly multiplexed immunoassays have allowed systematic large-scale measurement of hundreds of plasma proteins in large cohort studies. In combination with genotyping, such studies offer the prospect to 1) identify mechanisms involved with regulation of protein expression in plasma, and 2) determine whether the plasma proteins are likely to be causally implicated in disease. We report here the results of genome-wide association (GWA) studies of 83 proteins considered relevant to cardiovascular disease (CVD), measured in 3,394 individuals with multiple CVD risk factors. We identified 79 genome-wide significant (p<5e-8) association signals, 55 of which replicated at P<0.0007 in separate validation studies (n = 2,639 individuals). Using automated text mining, manual curation, and network-based methods incorporating information on expression quantitative trait loci (eQTL), we propose plausible causal mechanisms for 25 trans-acting loci, including a potential post-translational regulation of stem cell factor by matrix metalloproteinase 9 and receptor-ligand pairs such as RANK-RANK ligand. Using public GWA study data, we further evaluate all 79 loci for their causal effect on coronary artery disease, and highlight several potentially causal associations. Overall, a majority of the plasma proteins studied showed evidence of regulation at the genetic level. Our results enable future studies of the causal architecture of human disease, which in turn should aid discovery of new drug targets.