ABSTRACT
Inflammation is associated with a range of neuropsychiatric symptoms; however, the nature of the causal relationship is unclear. We used complementary non-genetic, genetic risk score (GRS), and Mendelian randomization (MR) analyses to examine whether inflammatory markers are associated with affect, depressive and anxiety disorders, and cognition. We tested in ≈ 55,098 (59% female) individuals from the Dutch Lifelines cohort the concurrent/prospective associations of C-reactive protein (CRP) with: depressive and anxiety disorders; positive/negative affect; and attention, psychomotor speed, episodic memory, and executive functioning. Additionally, we examined the association between inflammatory GRSs (CRP, interleukin-6 [IL-6], IL-6 receptor [IL-6R and soluble IL-6R (sIL-6R)], glycoprotein acetyls [GlycA]) on these same outcomes (N max =57,946), followed by MR analysis examining evidence of causality of CRP on outcomes (N max =23,268). In non-genetic analyses, higher CRP was associated with a depressive disorder, lower positive/higher negative affect, and worse executive function, attention, and psychomotor speed after adjusting for potential confounders. In genetic analyses, CRP GRS was associated with any anxiety disorder (ß = 0.002, p = 0.037) whereas GlycA GRS was associated with major depressive disorder (ß = 0.001, p = 0.036). Both CRP GRS (ß = 0.006, p = 0.035) and GlycA GRS (ß = 0.006, p = 0.049) were associated with greater negative affect. Inflammatory GRSs were not associated with cognition, except sIL-6R GRS which was associated with poorer memory (ß=-0.009, p = 0.018). There was weak evidence for a CRP-anxiety association using MR (ß = 0.12; p = 0.054). Genetic and non-genetic analyses provide consistent evidence for an association between CRP and negative affect. These results suggest that dysregulated immune physiology may impact a broad range of trans-diagnostic affective symptoms.
ABSTRACT
BACKGROUND: Individuals who have experienced a stroke, or transient ischemic attack, face a heightened risk of future cardiovascular events. Identification of genetic and molecular risk factors for subsequent cardiovascular outcomes may identify effective therapeutic targets to improve prognosis after an incident stroke. METHODS: We performed genome-wide association studies for subsequent major adverse cardiovascular events (MACE; ncases=51 929; ncontrols=39 980) and subsequent arterial ischemic stroke (AIS; ncases=45 120; ncontrols=46 789) after the first incident stroke within the Million Veteran Program and UK Biobank. We then used genetic variants associated with proteins (protein quantitative trait loci) to determine the effect of 1463 plasma protein abundances on subsequent MACE using Mendelian randomization. RESULTS: Two variants were significantly associated with subsequent cardiovascular events: rs76472767 near gene RNF220 (odds ratio, 0.75 [95% CI, 0.64-0.85]; P=3.69×10-8) with subsequent AIS and rs13294166 near gene LINC01492 (odds ratio, 1.52 [95% CI, 1.37-1.67]; P=3.77×10-8) with subsequent MACE. Using Mendelian randomization, we identified 2 proteins with an effect on subsequent MACE after a stroke: CCL27 ([C-C motif chemokine 27], effect odds ratio, 0.77 [95% CI, 0.66-0.88]; adjusted P=0.05) and TNFRSF14 ([tumor necrosis factor receptor superfamily member 14], effect odds ratio, 1.42 [95% CI, 1.24-1.60]; adjusted P=0.006). These proteins are not associated with incident AIS and are implicated to have a role in inflammation. CONCLUSIONS: We found evidence that 2 proteins with little effect on incident stroke appear to influence subsequent MACE after incident AIS. These associations suggest that inflammation is a contributing factor to subsequent MACE outcomes after incident AIS and highlights potential novel targets.
Subject(s)
Biological Specimen Banks , Genome-Wide Association Study , Mendelian Randomization Analysis , Stroke , Veterans , Humans , Male , Stroke/genetics , Stroke/epidemiology , Female , United Kingdom/epidemiology , Middle Aged , Aged , Disease Progression , Polymorphism, Single Nucleotide/genetics , Ischemic Stroke/genetics , Ischemic Stroke/epidemiology , Risk Factors , Quantitative Trait Loci , UK BiobankABSTRACT
Observational studies are rarely representative of their target population because there are known and unknown factors that affect an individual's choice to participate (the selection mechanism). Selection can cause bias in a given analysis if the outcome is related to selection (conditional on the other variables in the model). Detecting and adjusting for selection bias in practice typically requires access to data on nonselected individuals. Here, we propose methods to detect selection bias in genetic studies by comparing correlations among genetic variants in the selected sample to those expected under no selection. We examine the use of four hypothesis tests to identify induced associations between genetic variants in the selected sample. We evaluate these approaches in Monte Carlo simulations. Finally, we use these approaches in an applied example using data from the UK Biobank (UKBB). The proposed tests suggested an association between alcohol consumption and selection into UKBB. Hence, UKBB analyses with alcohol consumption as the exposure or outcome may be biased by this selection.
ABSTRACT
Epigenome-wide association studies (EWAS) aim to identify differentially methylated loci associated with complex traits and disorders. EWAS of cigarette smoking shows some of the most widespread DNA methylation (DNAm) associations in blood. However, traditional EWAS cannot differentiate between causation and confounding, leading to ambiguity in etiological interpretations. Here, we apply an integrated approach combining Mendelian Randomization and twin-based Direction-of-Causation analyses (MR-DoC) to examine causality underlying smoking-associated blood DNAm changes in the Netherlands Twin Register (N=2577). Evidence across models suggests that current smoking's causal effects on DNAm likely drive many of the previous EWAS findings, implicating functional pathways relevant to several adverse health outcomes of smoking, including hemopoiesis, cell- and neuro-development, and immune regulation. Additionally, we find evidence of potential reverse causal influences at some DNAm sites, with 17 of these sites enriched for gene regulatory functional elements in the brain. The top three sites with evidence of DNAm's effects on smoking annotate to genes involved in G protein-coupled receptor signaling (GNG7, RGS3) and innate immune response (SLC15A4), elucidating potential biological risk factors for smoking. This study highlights the utility of integrating genotypic and DNAm measures in twin cohorts to clarify the causal relationships between health behaviors and blood DNAm.
ABSTRACT
Background: Low-grade systemic inflammation is implicated in the pathogenesis of various neuropsychiatric conditions affecting mood and cognition. While much of the evidence concerns depression, large-scale population studies of anxiety, affect, and cognitive function are scarce. Importantly, causality remains unclear. We used complementary non-genetic, genetic risk score (GRS), and Mendelian randomization (MR) analyses to examine whether inflammatory markers are associated with affect, depressive and anxiety disorders, and cognitive performance in the Lifelines Cohort; and whether associations are likely to be causal. Methods: Using data from up to 55,098 (59% female) individuals from the Dutch Lifelines cohort, we tested the cross-sectional and longitudinal associations of C-reactive protein (CRP) with (i) depressive and anxiety disorders; (ii) positive and negative affect scores, and (iii) five cognitive measures assessing attention, psychomotor speed, episodic memory, and executive functioning (figural fluency and working memory). Additionally, we examined the association between inflammatory marker GRSs (CRP, interleukin-6 [IL-6], IL-6 receptor [IL-6R and soluble IL-6R (sIL-6R)], glycoprotein acetyls [GlycA]) on these same outcomes (Nmax=57,946), followed by MR analysis examining evidence of causality of CRP on outcomes (Nmax=23,268). In genetic analyses, all GRSs and outcomes were z-transformed. Results: In non-genetic analyses, higher CRP was associated with diagnosis of any depressive disorder, lower positive and higher negative affect scores, and worse performance on tests of figural fluency, attention, and psychomotor speed after adjusting for potential confounders, although the magnitude of these associations was small. In genetic analyses, CRPGRS was associated with any anxiety disorder (ß=0.002, p=0.037, N=57,047) whereas GlycAGRS was associated with major depressive disorder (ß=0.001, p=0.036; N=57,047). Both CRPGRS (ß=0.006, p=0.035, N=57,946) and GlycAGRS (ß=0.006, p=0.049; N=57,946) were associated with higher negative affect score. Inflammatory marker GRSs were not associated with cognitive performance, except sIL-6RGRS which was associated with poorer memory performance (ß=-0.009, p=0.018, N=36,783). Further examination of the CRP-anxiety association using MR provided some weak evidence of causality (ß=0.12; p=0.054). Conclusions: Genetic and non-genetic analyses provide consistent evidence for an association between CRP and negative affect. Genetic analyses suggest that IL-6 signaling could be relevant for memory, and that the association between CRP and anxiety disorders could be causal. These results suggest that dysregulated immune physiology may impact a broad range of trans-diagnostic affective symptoms. However, given the small effect sizes and multiple tests conducted, future studies are required to investigate whether effects are moderated by sub-groups and whether these findings replicate in other cohorts.
ABSTRACT
Individual sensitivity to environmental exposures may be genetically influenced. This genotype-by-environment interplay implies differences in phenotypic variance across genotypes. However, environmental sensitivity genetic variants have proven challenging to detect. GWAS of monozygotic twin differences is a family-based variance analysis method, which is more robust to systemic biases that impact population-based methods. We combined data from up to 21,792 monozygotic twins (10,896 pairs) from 11 studies to conduct the largest GWAS meta-analysis of monozygotic phenotypic differences in children and adolescents/adults for seven psychiatric and neurodevelopmental phenotypes: attention deficit hyperactivity disorder (ADHD) symptoms, autistic traits, anxiety and depression symptoms, psychotic-like experiences, neuroticism, and wellbeing. The SNP-heritability of variance in these phenotypes were estimated (h2: 0% to 18%), but were imprecise. We identified a total of 13 genome-wide significant associations (SNP, gene, and gene-set), including genes related to stress-reactivity for depression, growth factor-related genes for autistic traits and catecholamine uptake-related genes for psychotic-like experiences. Monozygotic twins are an important new source of evidence about the genetics of environmental sensitivity.
ABSTRACT
Population differences in cardiometabolic disease remain unexplained. Misleading assumptions over genetic explanations are partly due to terminology used to distinguish populations, specifically ancestry, race, and ethnicity. These terms differentially implicate environmental and biological causal pathways, which should inform their use. Genetic variation alone accounts for a limited fraction of population differences in cardiometabolic disease. Research effort should focus on societally driven, lifelong environmental determinants of population differences in disease. Rather than pursuing population stratifiers to personalize medicine, we advocate removing socioeconomic barriers to receipt of and adherence to healthcare interventions, which will have markedly greater impact on improving cardiometabolic outcomes. This requires multidisciplinary collaboration and public and policymaker engagement to address inequalities driven by society rather than biology per se.
Subject(s)
Cardiovascular Diseases , Ethnicity , Racial Groups , Humans , Cardiovascular Diseases/epidemiology , Metabolic Diseases/epidemiology , Metabolic Diseases/genetics , Genetic Predisposition to Disease , Socioeconomic Factors , Healthcare Disparities/ethnologyABSTRACT
BACKGROUND: Alzheimer's disease (AD)-related neuropathological changes can occur decades before clinical symptoms. We aimed to investigate whether neurodevelopment and/or neurodegeneration affects the risk of AD, through reducing structural brain reserve and/or increasing brain atrophy, respectively. METHODS: We used bidirectional two-sample Mendelian randomisation to estimate the effects between genetic liability to AD and global and regional cortical thickness, estimated total intracranial volume, volume of subcortical structures and total white matter in 37 680 participants aged 8-81 years across 5 independent cohorts (Adolescent Brain Cognitive Development, Generation R, IMAGEN, Avon Longitudinal Study of Parents and Children and UK Biobank). We also examined the effects of global and regional cortical thickness and subcortical volumes from the Enhancing NeuroImaging Genetics through Meta-Analysis (ENIGMA) Consortium on AD risk in up to 37 741 participants. RESULTS: Our findings show that AD risk alleles have an age-dependent effect on a range of cortical and subcortical brain measures that starts in mid-life, in non-clinical populations. Evidence for such effects across childhood and young adulthood is weak. Some of the identified structures are not typically implicated in AD, such as those in the striatum (eg, thalamus), with consistent effects from childhood to late adulthood. There was little evidence to suggest brain morphology alters AD risk. CONCLUSIONS: Genetic liability to AD is likely to affect risk of AD primarily through mechanisms affecting indicators of brain morphology in later life, rather than structural brain reserve. Future studies with repeated measures are required for a better understanding and certainty of the mechanisms at play.
ABSTRACT
Mendelian Randomisation (MR) estimates causal effects between risk factors and complex outcomes using genetic instruments. Pleiotropy, heritable confounders, and heterogeneous causal effects violate MR assumptions and can lead to biases. To alleviate these, we propose an approach employing a Phenome-Wide association Clustering of the MR instruments (PWC-MR) and apply this method to revisit the surprisingly large apparent causal effect of body mass index (BMI) on educational attainment (EDU): [Formula: see text] = -0.19 [-0.22, -0.16]. First, we cluster 324 BMI-associated genetic instruments based on their association with 407 traits in the UK Biobank, which yields six distinct groups. Subsequent cluster-specific MR reveals heterogeneous causal effect estimates on EDU. A cluster enriched for socio-economic indicators yields the largest BMI-on-EDU causal effect estimate ([Formula: see text] = -0.49 [-0.56, -0.42]) whereas a cluster enriched for body-mass specific traits provides a more likely estimate ([Formula: see text] = -0.09 [-0.13, -0.05]). Follow-up analyses confirms these findings: within-sibling MR ([Formula: see text] = -0.05 [-0.09, -0.01]); MR for childhood BMI on EDU ([Formula: see text] = -0.03 [-0.06, -0.002]); step-wise multivariable MR ([Formula: see text] = -0.05 [-0.07, -0.02]) where socio-economic indicators are jointly modelled. Here we show how the in-depth examination of the BMI-EDU causal relationship demonstrates the utility of our PWC-MR approach in revealing distinct pleiotropic pathways and confounder mechanisms.
Subject(s)
Genome-Wide Association Study , Mendelian Randomization Analysis , Humans , Child , Genome-Wide Association Study/methods , Mendelian Randomization Analysis/methods , Obesity/epidemiology , Obesity/genetics , Risk Factors , Educational Status , Polymorphism, Single NucleotideABSTRACT
Background: Individuals who have experienced a stroke, or transient ischemic attack, face a heightened risk of future cardiovascular events. Identification of genetic and molecular risk factors for subsequent cardiovascular outcomes may identify effective therapeutic targets to improve prognosis after an incident stroke. Methods: We performed genome-wide association studies (GWAS) for subsequent major adverse cardiovascular events (MACE) (Ncases=51,929, Ncntrl=39,980) and subsequent arterial ischemic stroke (AIS) Ncases=45,120, Ncntrl=46,789) after first incident stroke within the Million Veteran Program and UK Biobank. We then used genetic variants associated with proteins (pQTLs) to determine the effect of 1,463 plasma protein abundances on subsequent MACE using Mendelian randomization (MR). Results: Two variants were significantly associated with subsequent cardiovascular events: rs76472767 (OR=0.75, 95% CI = 0.64-0.85, p= 3.69×10-08) with subsequent AIS and rs13294166 (OR=1.52, 95% CI = 1.37-1.67, p=3.77×10-08) with subsequent MACE. Using MR, we identified 2 proteins with an effect on subsequent MACE after a stroke: CCL27 (effect OR= 0.77, 95% CI = 0.66-0.88, adj. p=0.05), and TNFRSF14 (effect OR=1.42, 95% CI = 1.24-1.60, adj. p=0.006). These proteins are not associated with incident AIS and are implicated to have a role in inflammation. Conclusions: We found evidence that two proteins with little effect on incident stroke appear to influence subsequent MACE after incident AIS. These associations suggest that inflammation is a contributing factor to subsequent MACE outcomes after incident AIS and highlights potential novel targets.
ABSTRACT
BACKGROUND: Although it is known that variation in the aldehyde dehydrogenase 2 (ALDH2) gene family influences the East Asian alcohol flushing response, knowledge about other genetic variants that affect flushing symptoms is limited. METHODS: We performed a genome-wide association study meta-analysis and heritability analysis of alcohol flushing in 15,105 males of East Asian ancestry (Koreans and Chinese) to identify genetic associations with alcohol flushing. We also evaluated whether self-reported flushing can be used as an instrumental variable for alcohol intake. RESULTS: We identified variants in the region of ALDH2 strongly associated with alcohol flushing, replicating previous studies conducted in East Asian populations. Additionally, we identified variants in the alcohol dehydrogenase 1B (ADH1B) gene region associated with alcohol flushing. Several novel variants were identified after adjustment for the lead variants (ALDH2-rs671 and ADH1B-rs1229984), which need to be confirmed in larger studies. The estimated SNP-heritability on the liability scale was 13% (S.E. = 4%) for flushing, but the heritability estimate decreased to 6% (S.E. = 4%) when the effects of the lead variants were controlled for. Genetic instrumentation of higher alcohol intake using these variants recapitulated known associations of alcohol intake with hypertension. Using self-reported alcohol flushing as an instrument gave a similar association pattern of higher alcohol intake and cardiovascular disease-related traits (e.g. stroke). CONCLUSION: This study confirms that ALDH2-rs671 and ADH1B-rs1229984 are associated with alcohol flushing in East Asian populations. Our findings also suggest that self-reported alcohol flushing can be used as an instrumental variable in future studies of alcohol consumption.
Subject(s)
Alcohol Drinking , East Asian People , Flushing , Humans , Male , Alcohol Dehydrogenase/genetics , Alcohol Drinking/genetics , Aldehyde Dehydrogenase, Mitochondrial/genetics , East Asian People/genetics , Genome-Wide Association Study , Polymorphism, Single Nucleotide , Flushing/chemically inducedABSTRACT
Genome-wide association studies (GWAS) are commonly used to identify genomic variants that are associated with complex traits, and estimate the magnitude of this association for each variant. However, it has been widely observed that the association estimates of variants tend to be lower in a replication study than in the study that discovered those associations. A phenomenon known as Winner's Curse is responsible for this upward bias present in association estimates of significant variants in the discovery study. We review existing Winner's Curse correction methods which require only GWAS summary statistics in order to make adjustments. In addition, we propose modifications to improve existing methods and propose a novel approach which uses the parametric bootstrap. We evaluate and compare methods, first using a wide variety of simulated data sets and then, using real data sets for three different traits. The metric, estimated mean squared error (MSE) over significant SNPs, was primarily used for method assessment. Our results indicate that widely used conditional likelihood based methods tend to perform poorly. The other considered methods behave much more similarly, with our proposed bootstrap method demonstrating very competitive performance. To complement this review, we have developed an R package, 'winnerscurse' which can be used to implement these various Winner's Curse adjustment methods to GWAS summary statistics.
Subject(s)
Genome-Wide Association Study , Polymorphism, Single Nucleotide , Genome-Wide Association Study/methods , Likelihood Functions , Genetic Association Studies , Bias , Phenotype , Polymorphism, Single Nucleotide/geneticsABSTRACT
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).
Subject(s)
Colorectal Neoplasms , Esophageal Neoplasms , Esophageal Squamous Cell Carcinoma , Fatty Acids, Omega-3 , Inflammatory Bowel Diseases , Humans , Fatty Acid Desaturases/genetics , Fatty Acid Desaturases/metabolism , Fatty Acids, Unsaturated/metabolism , Polymorphism, Single NucleotideABSTRACT
Participant overlap can induce overfitting bias into Mendelian randomization (MR) and polygenic risk score (PRS) studies. Here, we evaluated a block jackknife resampling framework for genome-wide association studies (GWAS) and PRS construction to mitigate overfitting bias in MR analyses and implemented this study design in a causal inference setting using data from the UK Biobank. We simulated PRS and MR under three scenarios: (1) using weighted SNP estimates from an external GWAS, (2) using weighted SNP estimates from an overlapping GWAS sample and (3) using a block jackknife resampling framework. Based on a P-value threshold to derive genetic instruments for MR studies (P < 5 × 10-8) and a 10% variance in the exposure explained by all SNPs, block-jackknifing PRS did not suffer from overfitting bias (mean R2 = 0.034) compared with the externally weighted PRS (mean R2 = 0.040). In contrast, genetic instruments derived from overlapping samples explained a higher variance (mean R2 = 0.048) compared with the externally derived score. Overfitting became considerably more severe when using a more liberal P-value threshold to construct PRS (e.g. P < 0.05, overlapping sample PRS mean R2 = 0.103, externally weighted PRS mean R2 = 0.086), whereas estimates using jackknife score remained robust to overfitting (mean R2 = 0.084). Using block jackknife resampling MR in an applied analysis, we examined the effects of body mass index on circulating biomarkers which provided comparable estimates to an externally weighted instrument, whereas the overfitted scores typically provided narrower confidence intervals. Furthermore, we extended this framework into sex-stratified, multivariate and bidirectional settings to investigate the effect of childhood body size on adult testosterone levels.
Subject(s)
Genome-Wide Association Study , Mendelian Randomization Analysis , Adult , Humans , Risk Factors , Body Mass Index , Polymorphism, Single Nucleotide/geneticsABSTRACT
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).
Subject(s)
Genome-Wide Association Study , Neoplasms , Humans , Mendelian Randomization Analysis , Reproducibility of Results , Fatty Acids , Quality Control , Neoplasms/epidemiology , Neoplasms/geneticsABSTRACT
Identifying genomic regions pertinent to complex traits is a common goal of genome-wide and epigenome-wide association studies (GWAS and EWAS). GWAS identify causal genetic variants, directly or via linkage disequilibrium, and EWAS identify variation in DNA methylation associated with a trait. While GWAS in principle will only detect variants due to causal genes, EWAS can also identify genes via confounding, or reverse causation. We systematically compare GWAS (N > 50,000) and EWAS (N > 4500) results of 15 complex traits. We evaluate if the genes or gene ontology terms flagged by GWAS and EWAS overlap, and find substantial overlap for diastolic blood pressure, (gene overlap P = 5.2 × 10-6; term overlap P = 0.001). We superimpose our empirical findings against simulated models of varying genetic and epigenetic architectures and observe that in most cases GWAS and EWAS are likely capturing distinct genesets. Our results indicate that GWAS and EWAS are capturing different aspects of the biology of complex traits.
Subject(s)
Epigenome , Genome-Wide Association Study , Genome-Wide Association Study/methods , Epigenesis, Genetic , Multifactorial Inheritance , DNA MethylationABSTRACT
Proteome-wide Mendelian randomization (MR) shows value in prioritizing drug targets in Europeans but with limited evidence in other ancestries. Here, we present a multi-ancestry proteome-wide MR analysis based on cross-population data from the Global Biobank Meta-analysis Initiative (GBMI). We estimated the putative causal effects of 1,545 proteins on eight diseases in African (32,658) and European (1,219,993) ancestries and identified 45 and 7 protein-disease pairs with MR and genetic colocalization evidence in the two ancestries, respectively. A multi-ancestry MR comparison identified two protein-disease pairs with MR evidence in both ancestries and seven pairs with specific effects in the two ancestries separately. Integrating these MR signals with clinical trial evidence, we prioritized 16 pairs for investigation in future drug trials. Our results highlight the value of proteome-wide MR in informing the generalizability of drug targets for disease prevention across ancestries and illustrate the value of meta-analysis of biobanks in drug development.
ABSTRACT
Despite early interest, the evidence linking fatty acids to cardiovascular diseases (CVDs) remains controversial. We used Mendelian randomization to explore the involvement of polyunsaturated (PUFA) and monounsaturated (MUFA) fatty acids biosynthesis in the etiology of several CVD endpoints in up to 1 153 768 European (maximum 123 668 cases) and 212 453 East Asian (maximum 29 319 cases) ancestry individuals. As instruments, we selected single nucleotide polymorphisms mapping to genes with well-known roles in PUFA (i.e. FADS1/2 and ELOVL2) and MUFA (i.e. SCD) biosynthesis. Our findings suggest that higher PUFA biosynthesis rate (proxied by rs174576 near FADS1/2) is related to higher odds of multiple CVDs, particularly ischemic stroke, peripheral artery disease and venous thromboembolism, whereas higher MUFA biosynthesis rate (proxied by rs603424 near SCD) is related to lower odds of coronary artery disease among Europeans. Results were unclear for East Asians as most effect estimates were imprecise. By triangulating multiple approaches (i.e. uni-/multi-variable Mendelian randomization, a phenome-wide scan, genetic colocalization and within-sibling analyses), our results are compatible with higher low-density lipoprotein (LDL) cholesterol (and possibly glucose) being a downstream effect of higher PUFA biosynthesis rate. Our findings indicate that PUFA and MUFA biosynthesis are involved in the etiology of CVDs and suggest LDL cholesterol as a potential mediating trait between PUFA biosynthesis and CVDs risk.
Subject(s)
Cardiovascular Diseases , Humans , Cardiovascular Diseases/genetics , Mendelian Randomization Analysis , Fatty Acids/genetics , Asian People/genetics , Polymorphism, Single Nucleotide/geneticsABSTRACT
Estimating effects of parental and sibling genotypes (indirect genetic effects) can provide insight into how the family environment influences phenotypic variation. There is growing molecular genetic evidence for effects of parental phenotypes on their offspring (e.g. parental educational attainment), but the extent to which siblings affect each other is currently unclear. Here we used data from samples of unrelated individuals, without (singletons) and with biological full-siblings (non-singletons), to investigate and estimate sibling effects. Indirect genetic effects of siblings increase (or decrease) the covariance between genetic variation and a phenotype. It follows that differences in genetic association estimates between singletons and non-singletons could indicate indirect genetic effects of siblings if there is no heterogeneity in other sources of genetic association between singletons and non-singletons. We used UK Biobank data to estimate polygenic score (PGS) associations for height, BMI and educational attainment in self-reported singletons (N = 50,143) and non-singletons (N = 328,549). The educational attainment PGS association estimate was 12% larger (95% C.I. 3%, 21%) in the non-singleton sample than in the singleton sample, but the height and BMI PGS associations were consistent. Birth order data suggested that the difference in educational attainment PGS associations was driven by individuals with older siblings rather than firstborns. The relationship between number of siblings and educational attainment PGS associations was non-linear; PGS associations were 24% smaller in individuals with 6 or more siblings compared to the rest of the sample (95% C.I. 11%, 38%). We estimate that a 1 SD increase in sibling educational attainment PGS corresponds to a 0.025 year increase in the index individual's years in schooling (95% C.I. 0.013, 0.036). Our results suggest that older siblings may influence the educational attainment of younger siblings, adding to the growing evidence that effects of the environment on phenotypic variation partially reflect social effects of germline genetic variation in relatives.