PRSet: Pathway-based polygenic risk score analyses and software.

Polygenic risk scores (PRSs) have been among the leading advances in biomedicine in recent years. As a proxy of genetic liability, PRSs are utilised across multiple fields and applications. While numerous statistical and machine learning methods have been developed to optimise their predictive accuracy, these typically distil genetic liability to a single number based on aggregation of an individual's genome-wide risk alleles. This results in a key loss of information about an individual's genetic profile, which could be critical given the functional sub-structure of the genome and the heterogeneity of complex disease. In this manuscript, we introduce a 'pathway polygenic' paradigm of disease risk, in which multiple genetic liabilities underlie complex diseases, rather than a single genome-wide liability. We describe a method and accompanying software, PRSet, for computing and analysing pathway-based PRSs, in which polygenic scores are calculated across genomic pathways for each individual. We evaluate the potential of pathway PRSs in two distinct ways, creating two major sections: (1) In the first section, we benchmark PRSet as a pathway enrichment tool, evaluating its capacity to capture GWAS signal in pathways. We find that for target sample sizes of >10,000 individuals, pathway PRSs have similar power for evaluating pathway enrichment as leading methods MAGMA and LD score regression, with the distinct advantage of providing individual-level estimates of genetic liability for each pathway -opening up a range of pathway-based PRS applications, (2) In the second section, we evaluate the performance of pathway PRSs for disease stratification. We show that using a supervised disease stratification approach, pathway PRSs (computed by PRSet) outperform two standard genome-wide PRSs (computed by C+T and lassosum) for classifying disease subtypes in 20 of 21 scenarios tested. As the definition and functional annotation of pathways becomes increasingly refined, we expect pathway PRSs to offer key insights into the heterogeneity of complex disease and treatment response, to generate biologically tractable therapeutic targets from polygenic signal, and, ultimately, to provide a powerful path to precision medicine.

Genetic sensitivity analysis: Adjusting for genetic confounding in epidemiological associations.

Associations between exposures and outcomes reported in epidemiological studies are typically unadjusted for genetic confounding. We propose a two-stage approach for estimating the degree to which such observed associations can be explained by genetic confounding. First, we assess attenuation of exposure effects in regressions controlling for increasingly powerful polygenic scores. Second, we use structural equation models to estimate genetic confounding using heritability estimates derived from both SNP-based and twin-based studies. We examine associations between maternal education and three developmental outcomes - child educational achievement, Body Mass Index, and Attention Deficit Hyperactivity Disorder. Polygenic scores explain between 14.3% and 23.0% of the original associations, while analyses under SNP- and twin-based heritability scenarios indicate that observed associations could be almost entirely explained by genetic confounding. Thus, caution is needed when interpreting associations from non-genetically informed epidemiology studies. Our approach, akin to a genetically informed sensitivity analysis can be applied widely.

Genome-wide gene-environment analyses of major depressive disorder and reported lifetime traumatic experiences in UK Biobank.

Depression is more frequent among individuals exposed to traumatic events. Both trauma exposure and depression are heritable. However, the relationship between these traits, including the role of genetic risk factors, is complex and poorly understood. When modelling trauma exposure as an environmental influence on depression, both gene-environment correlations and gene-environment interactions have been observed. The UK Biobank concurrently assessed Major Depressive Disorder (MDD) and self-reported lifetime exposure to traumatic events in 126,522 genotyped individuals of European ancestry. We contrasted genetic influences on MDD stratified by reported trauma exposure (final sample size range: 24,094-92,957). The SNP-based heritability of MDD with reported trauma exposure (24%) was greater than MDD without reported trauma exposure (12%). Simulations showed that this is not confounded by the strong, positive genetic correlation observed between MDD and reported trauma exposure. We also observed that the genetic correlation between MDD and waist circumference was only significant in individuals reporting trauma exposure (rg = 0.24, p = 1.8 × 10-7 versus rg = -0.05, p = 0.39 in individuals not reporting trauma exposure, difference p = 2.3 × 10-4). Our results suggest that the genetic contribution to MDD is greater when reported trauma is present, and that a complex relationship exists between reported trauma exposure, body composition, and MDD.

Studying individual risk factors for self-harm in the UK Biobank: A polygenic scoring and Mendelian randomisation study.

BACKGROUND: Identifying causal risk factors for self-harm is essential to inform preventive interventions. Epidemiological studies have identified risk factors associated with self-harm, but these associations can be subject to confounding. By implementing genetically informed methods to better account for confounding, this study aimed to better identify plausible causal risk factors for self-harm. METHODS AND FINDINGS: Using summary statistics from 24 genome-wide association studies (GWASs) comprising 16,067 to 322,154 individuals, polygenic scores (PSs) were generated to index 24 possible individual risk factors for self-harm (i.e., mental health vulnerabilities, substance use, cognitive traits, personality traits, and physical traits) among a subset of UK Biobank participants (N = 125,925, 56.2% female) who completed an online mental health questionnaire in the period from 13 July 2016 to 27 July 2017. In total, 5,520 (4.4%) of these participants reported having self-harmed in their lifetime. In binomial regression models, PSs indexing 6 risk factors (major depressive disorder [MDD], attention deficit/hyperactivity disorder [ADHD], bipolar disorder, schizophrenia, alcohol dependence disorder, and lifetime cannabis use) predicted self-harm, with effect sizes ranging from odds ratio (OR) = 1.05 (95% CI 1.02 to 1.07, q = 0.008) for lifetime cannabis use to OR = 1.20 (95% CI 1.16 to 1.23, q = 1.33 × 10-35) for MDD. No systematic differences emerged between suicidal and non-suicidal self-harm. To further probe causal relationships, two-sample Mendelian randomisation (MR) analyses were conducted, with MDD, ADHD, and schizophrenia emerging as the most plausible causal risk factors for self-harm. The genetic liabilities for MDD and schizophrenia were associated with self-harm independently of diagnosis and medication. Main limitations include the lack of representativeness of the UK Biobank sample, that self-harm was self-reported, and the limited power of some of the included GWASs, potentially leading to possible type II error. CONCLUSIONS: In addition to confirming the role of MDD, we demonstrate that ADHD and schizophrenia likely play a role in the aetiology of self-harm using multivariate genetic designs for causal inference. Among the many individual risk factors we simultaneously considered, our findings suggest that systematic detection and treatment of core psychiatric symptoms, including psychotic and impulsivity symptoms, may be beneficial among people at risk for self-harm.

Genetic comorbidity between major depression and cardio-metabolic traits, stratified by age at onset of major depression.

It is imperative to understand the specific and shared etiologies of major depression and cardio-metabolic disease, as both traits are frequently comorbid and each represents a major burden to society. This study examined whether there is a genetic association between major depression and cardio-metabolic traits and if this association is stratified by age at onset for major depression. Polygenic risk scores analysis and linkage disequilibrium score regression was performed to examine whether differences in shared genetic etiology exist between depression case control status (N cases = 40,940, N controls = 67,532), earlier (N = 15,844), and later onset depression (N = 15,800) with body mass index, coronary artery disease, stroke, and type 2 diabetes in 11 data sets from the Psychiatric Genomics Consortium, Generation Scotland, and UK Biobank. All cardio-metabolic polygenic risk scores were associated with depression status. Significant genetic correlations were found between depression and body mass index, coronary artery disease, and type 2 diabetes. Higher polygenic risk for body mass index, coronary artery disease, and type 2 diabetes was associated with both early and later onset depression, while higher polygenic risk for stroke was associated with later onset depression only. Significant genetic correlations were found between body mass index and later onset depression, and between coronary artery disease and both early and late onset depression. The phenotypic associations between major depression and cardio-metabolic traits may partly reflect their overlapping genetic etiology irrespective of the age depression first presents.

Polygenic scores via penalized regression on summary statistics.

Polygenic scores (PGS) summarize the genetic contribution of a person's genotype to a disease or phenotype. They can be used to group participants into different risk categories for diseases, and are also used as covariates in epidemiological analyses. A number of possible ways of calculating PGS have been proposed, and recently there is much interest in methods that incorporate information available in published summary statistics. As there is no inherent information on linkage disequilibrium (LD) in summary statistics, a pertinent question is how we can use LD information available elsewhere to supplement such analyses. To answer this question, we propose a method for constructing PGS using summary statistics and a reference panel in a penalized regression framework, which we call lassosum. We also propose a general method for choosing the value of the tuning parameter in the absence of validation data. In our simulations, we showed that pseudovalidation often resulted in prediction accuracy that is comparable to using a dataset with validation phenotype and was clearly superior to the conservative option of setting the tuning parameter of lassosum to its lowest value. We also showed that lassosum achieved better prediction accuracy than simple clumping and P-value thresholding in almost all scenarios. It was also substantially faster and more accurate than the recently proposed LDpred.

Amelioration of X-Linked Related Autophagy Failure in Danon Disease With DNA Methylation Inhibitor.

BACKGROUND: Danon disease is an X-linked disorder that leads to fatal cardiomyopathy caused by a deficiency in lysosome-associated membrane protein-2 (LAMP2). In female patients, a later onset and less severe clinical phenotype have been attributed to the random inactivation of the X chromosome carrying the mutant diseased allele. We generated a patient-specific induced pluripotent stem cell (iPSCs)-based model of Danon disease to evaluate the therapeutic potential of Xi-chromosome reactivation using a DNA methylation inhibitor. METHODS: Using whole-exome sequencing, we identified a nonsense mutation (c.520C>T, exon 4) of the LAMP2 gene in a family with Danon disease. We generated iPSC lines from somatic cells derived from the affected mother and her 2 sons, and we then differentiated them into cardiomyocytes (iPSC-CMs) for modeling the histological and functional signatures, including autophagy failure of Danon disease. RESULTS: Our iPSC-CM platform provides evidence that random inactivation of the wild-type and mutant LAMP2 alleles on the X chromosome is responsible for the unusual phenotype in female patients with Danon disease. In vitro, iPSC-CMs from these patients reproduced the histological features and autophagy failure of Danon disease. Administration of the DNA demethylating agent 5-aza-2'-deoxycytidine reactivated the silent LAMP2 allele in iPSCs and iPSC-CMs in female patients with Danon disease and ameliorated their autophagy failure, supporting the application of a patient-specific iPSC platform for disease modeling and drug screening. CONCLUSIONS: Our iPSC-CM platform provides novel mechanistic and therapeutic insights into the contribution of random X chromosome inactivation to disease phenotype in X-linked Danon disease.

Novel pre-mRNA splicing of intronically integrated HBV generates oncogenic chimera in hepatocellular carcinoma.

BACKGROUND & AIMS: Hepatitis B virus (HBV) integration is common in HBV-associated hepatocellular carcinoma (HCC) and may play an important pathogenic role through the production of chimeric HBV-human transcripts. We aimed to screen the transcriptome for HBV integrations in HCCs. METHODS: Transcriptome sequencing was performed on paired HBV-associated HCCs and corresponding non-tumorous liver tissues to identify viral-human chimeric sites. Validation was further performed in an expanded cohort of human HCCs. RESULTS: Here we report the discovery of a novel pre-mRNA splicing mechanism in generating HBV-human chimeric protein. This mechanism was exemplified by the formation of a recurrent HBV-cyclin A2 (CCNA2) chimeric transcript (A2S), as detected in 12.5% (6 of 48) of HCC patients, but in none of the 22 non-HCC HBV-associated cirrhotic liver samples examined. Upon the integration of HBV into the intron of the CCNA2 gene, the mammalian splicing machinery utilized the foreign splice sites at 282nt. and 458nt. of the HBV genome to generate a pseudo-exon, forming an in-frame chimeric fusion with CCNA2. The A2S chimeric protein gained a non-degradable property and promoted cell cycle progression, demonstrating its potential oncogenic functions. CONCLUSIONS: A pre-mRNA splicing mechanism is involved in the formation of HBV-human chimeric proteins. This represents a novel and possibly common mechanism underlying the formation of HBV-human chimeric transcripts from intronically integrated HBV genome with functional impact. LAY SUMMARY: HBV is involved in the mammalian pre-mRNA splicing machinery in the generation of potential tumorigenic HBV-human chimeras. This study also provided insight on the impact of intronic HBV integration with the gain of splice sites in the development of HBV-associated HCC.

Patient-specific induced-pluripotent stem cells-derived cardiomyocytes recapitulate the pathogenic phenotypes of dilated cardiomyopathy due to a novel DES mutation identified by whole exome sequencing.

In this paper, we report a novel heterozygous mutation of A285V codon conversion on exon 4 of the desmin (DES), using whole exome sequencing (WES) in an isolated proband with documented dilated cardiomyopathy (DCM). This mutation is predicted to cause three-dimensional structure changes of DES. Immunohistological and electron microscopy studies demonstrated diffuse abnormal DES aggregations in DCM-induced-pluripotent stem cell (iPSC)-derived cardiomyocytes, and control-iPSC-derived cardiomyocytes transduced with A285V-DES. DCM-iPSC-derived cardiomyocytes also exhibited functional abnormalities in vitro. This is the first demonstration that patient-specific iPSC-derived cardiomyocytes can be used to provide histological and functional confirmation of a suspected genetic basis for DCM identified by WES.

BridgePRS leverages shared genetic effects across ancestries to increase polygenic risk score portability.

Here we present BridgePRS, a novel Bayesian polygenic risk score (PRS) method that leverages shared genetic effects across ancestries to increase PRS portability. We evaluate BridgePRS via simulations and real UK Biobank data across 19 traits in individuals of African, South Asian and East Asian ancestry, using both UK Biobank and Biobank Japan genome-wide association study summary statistics; out-of-cohort validation is performed in the Mount Sinai (New York) BioMe biobank. BridgePRS is compared with the leading alternative, PRS-CSx, and two other PRS methods. Simulations suggest that the performance of BridgePRS relative to PRS-CSx increases as uncertainty increases: with lower trait heritability, higher polygenicity and greater between-population genetic diversity; and when causal variants are not present in the data. In real data, BridgePRS has a 61% larger average R2 than PRS-CSx in out-of-cohort prediction of African ancestry samples in BioMe (P = 6 × 10-5). BridgePRS is a computationally efficient, user-friendly and powerful approach for PRS analyses in non-European ancestries.

BridgePRS : A powerful trans-ancestry Polygenic Risk Score method.

Polygenic Risk Scores (PRS) have huge potential to contribute to biomedical research and to a future of precision medicine, but to date their calculation relies largely on Europeanancestry GWAS data. This global bias makes most PRS substantially less accurate in individuals of non-European ancestry. Here we present BridgePRS , a novel Bayesian PRS method that leverages shared genetic effects across ancestries to increase the accuracy of PRS in non-European populations. The performance of BridgePRS is evaluated in simulated data and real UK Biobank (UKB) data across 19 traits in African, South Asian and East Asian ancestry individuals, using both UKB and Biobank Japan GWAS summary statistics. BridgePRS is compared to the leading alternative, PRS-CSx , and two single-ancestry PRS methods adapted for trans-ancestry prediction. PRS trained in the UK Biobank are then validated out-of-cohort in the independent Mount Sinai (New York) Bio Me Biobank. Simulations reveal that BridgePRS performance, relative to PRS-CSx , increases as uncertainty increases: with lower heritability, higher polygenicity, greater between-population genetic diversity, and when causal variants are not present in the data. Our simulation results are consistent with real data analyses in which BridgePRS has better predictive accuracy in African ancestry samples, especially in out-of-cohort prediction (into Bio Me ), which shows a 60% boost in mean R 2 compared to PRS-CSx ( P = 2 × 10 -6 ). BridgePRS performs the full PRS analysis pipeline, is computationally efficient, and is a powerful method for deriving PRS in diverse and under-represented ancestry populations.

EraSOR: a software tool to eliminate inflation caused by sample overlap in polygenic score analyses.

BACKGROUND: Polygenic risk score (PRS) analyses are now routinely applied across biomedical research. However, as PRS studies grow in size, there is an increased risk of sample overlap between the genome-wide association study (GWAS) from which the PRS is derived and the "target sample," in which PRSs are computed and hypotheses are tested. Despite the wide recognition of the sample overlap problem, its potential impact on the results from PRS studies has not yet been quantified, and no analytical solution has been provided. FINDINGS: Here, we first conduct a comprehensive investigation into the scale of the sample overlap problem, finding that PRS results can be substantially inflated even in the presence of minimal overlap. Next, we introduce a method and software, EraSOR (Erase Sample Overlap and Relatedness), which eliminates the inflation caused by sample overlap (and close relatedness) in almost all settings tested here. CONCLUSIONS: EraSOR could be useful in PRS studies (with target sample >1,000) similar to those investigated here, either (i) to mitigate the potential effects of known or unknown intercohort overlap and close relatedness or (ii) as a sensitivity tool to highlight the possible presence of sample overlap before its direct removal, when possible, or else to provide a lower bound on PRS analysis results after accounting for potential sample overlap.

Disentangling Independent and Mediated Causal Relationships Between Blood Metabolites, Cognitive Factors, and Alzheimer's Disease.

Background: Education and cognition demonstrate consistent inverse associations with Alzheimer's disease (AD). The biological underpinnings, however, remain unclear. Blood metabolites reflect the end point of biological processes and are accessible and malleable. Identifying metabolites with etiological relevance to AD and disentangling how these relate to cognitive factors along the AD causal pathway could, therefore, offer unique insights into underlying causal mechanisms. Methods: Using data from the largest metabolomics genome-wide association study (N ≈ 24,925) and three independent AD cohorts (N = 4725), cross-trait polygenic scores were generated and meta-analyzed. Metabolites genetically associated with AD were taken forward for causal analyses. Bidirectional two-sample Mendelian randomization interrogated univariable causal relationships between 1) metabolites and AD; 2) education and cognition; 3) metabolites, education, and cognition; and 4) education, cognition, and AD. Mediating relationships were computed using multivariable Mendelian randomization. Results: Thirty-four metabolites were genetically associated with AD at p < .05. Of these, glutamine and free cholesterol in extra-large high-density lipoproteins demonstrated a protective causal effect (glutamine: 95% confidence interval [CI], 0.70 to 0.92; free cholesterol in extra-large high-density lipoproteins: 95% CI, 0.75 to 0.92). An AD-protective effect was also observed for education (95% CI, 0.61 to 0.85) and cognition (95% CI, 0.60 to 0.89), with bidirectional mediation evident. Cognition as a mediator of the education-AD relationship was stronger than vice versa, however. No evidence of mediation via any metabolite was found. Conclusions: Glutamine and free cholesterol in extra-large high-density lipoproteins show protective causal effects on AD. Education and cognition also demonstrate protection, though education's effect is almost entirely mediated by cognition. These insights provide key pieces of the AD causal puzzle, important for informing future multimodal work and progressing toward effective intervention strategies.

The Value of Rare Genetic Variation in the Prediction of Common Obesity in European Ancestry Populations.

Polygenic risk scores (PRSs) aggregate the effects of genetic variants across the genome and are used to predict risk of complex diseases, such as obesity. Current PRSs only include common variants (minor allele frequency (MAF) ≥1%), whereas the contribution of rare variants in PRSs to predict disease remains unknown. Here, we examine whether augmenting the standard common variant PRS (PRScommon) with a rare variant PRS (PRSrare) improves prediction of obesity. We used genome-wide genotyped and imputed data on 451,145 European-ancestry participants of the UK Biobank, as well as whole exome sequencing (WES) data on 184,385 participants. We performed single variant analyses (for both common and rare variants) and gene-based analyses (for rare variants) for association with BMI (kg/m2), obesity (BMI ≥ 30 kg/m2), and extreme obesity (BMI ≥ 40 kg/m2). We built PRSscommon and PRSsrare using a range of methods (Clumping+Thresholding [C+T], PRS-CS, lassosum, gene-burden test). We selected the best-performing PRSs and assessed their performance in 36,757 European-ancestry unrelated participants with whole genome sequencing (WGS) data from the Trans-Omics for Precision Medicine (TOPMed) program. The best-performing PRScommon explained 10.1% of variation in BMI, and 18.3% and 22.5% of the susceptibility to obesity and extreme obesity, respectively, whereas the best-performing PRSrare explained 1.49%, and 2.97% and 3.68%, respectively. The PRSrare was associated with an increased risk of obesity and extreme obesity (ORobesity = 1.37 per SDPRS, Pobesity = 1.7x10-85; ORextremeobesity = 1.55 per SDPRS, Pextremeobesity = 3.8x10-40), which was attenuated, after adjusting for PRScommon (ORobesity = 1.08 per SDPRS, Pobesity = 9.8x10-6; ORextremeobesity= 1.09 per SDPRS, Pextremeobesity = 0.02). When PRSrare and PRScommon are combined, the increase in explained variance attributed to PRSrare was small (incremental Nagelkerke R2 = 0.24% for obesity and 0.51% for extreme obesity). Consistently, combining PRSrare to PRScommon provided little improvement to the prediction of obesity (PRSrare AUC = 0.591; PRScommon AUC = 0.708; PRScombined AUC = 0.710). In summary, while rare variants show convincing association with BMI, obesity and extreme obesity, the PRSrare provides limited improvement over PRScommon in the prediction of obesity risk, based on these large populations.

Genetic analysis of amyotrophic lateral sclerosis identifies contributing pathways and cell types.

Despite the considerable progress in unraveling the genetic causes of amyotrophic lateral sclerosis (ALS), we do not fully understand the molecular mechanisms underlying the disease. We analyzed genome-wide data involving 78,500 individuals using a polygenic risk score approach to identify the biological pathways and cell types involved in ALS. This data-driven approach identified multiple aspects of the biology underlying the disease that resolved into broader themes, namely, neuron projection morphogenesis, membrane trafficking, and signal transduction mediated by ribonucleotides. We also found that genomic risk in ALS maps consistently to GABAergic interneurons and oligodendrocytes, as confirmed in human single-nucleus RNA-seq data. Using two-sample Mendelian randomization, we nominated six differentially expressed genes (ATG16L2, ACSL5, MAP1LC3A, MAPKAPK3, PLXNB2, and SCFD1) within the significant pathways as relevant to ALS. We conclude that the disparate genetic etiologies of this fatal neurological disease converge on a smaller number of final common pathways and cell types.

Identification of novel locus associated with coronary artery aneurysms and validation of loci for susceptibility to Kawasaki disease.

Kawasaki disease (KD) is a paediatric vasculitis associated with coronary artery aneurysms (CAA). Genetic variants influencing susceptibility to KD have been previously identified, but no risk alleles have been validated that influence CAA formation. We conducted a genome-wide association study (GWAS) for CAA in KD patients of European descent with 200 cases and 276 controls. A second GWAS for susceptibility pooled KD cases with healthy paediatric controls from vaccine trials in the UK (n = 1609). Logistic regression mixed models were used for both GWASs. The susceptibility GWAS was meta-analysed with 400 KD cases and 6101 controls from a previous European GWAS, these results were further meta-analysed with Japanese GWASs at two putative loci. The CAA GWAS identified an intergenic region of chromosome 20q13 with multiple SNVs showing genome-wide significance. The risk allele of the most associated SNV (rs6017006) was present in 13% of cases and 4% of controls; in East Asian 1000 Genomes data, the allele was absent or rare. Susceptibility GWAS with meta-analysis with previously published European data identified two previously associated loci (ITPKC and FCGR2A). Further meta-analysis with Japanese GWAS summary data from the CASP3 and FAM167A genomic regions validated these loci in Europeans showing consistent effects of the top SNVs in both populations. We identified a novel locus for CAA in KD patients of European descent. The results suggest that different genes determine susceptibility to KD and development of CAA and future work should focus on the function of the intergenic region on chromosome 20q13.

Tutorial: a guide to performing polygenic risk score analyses.

A polygenic score (PGS) or polygenic risk score (PRS) is an estimate of an individual's genetic liability to a trait or disease, calculated according to their genotype profile and relevant genome-wide association study (GWAS) data. While present PRSs typically explain only a small fraction of trait variance, their correlation with the single largest contributor to phenotypic variation-genetic liability-has led to the routine application of PRSs across biomedical research. Among a range of applications, PRSs are exploited to assess shared etiology between phenotypes, to evaluate the clinical utility of genetic data for complex disease and as part of experimental studies in which, for example, experiments are performed that compare outcomes (e.g., gene expression and cellular response to treatment) between individuals with low and high PRS values. As GWAS sample sizes increase and PRSs become more powerful, PRSs are set to play a key role in research and stratified medicine. However, despite the importance and growing application of PRSs, there are limited guidelines for performing PRS analyses, which can lead to inconsistency between studies and misinterpretation of results. Here, we provide detailed guidelines for performing and interpreting PRS analyses. We outline standard quality control steps, discuss different methods for the calculation of PRSs, provide an introductory online tutorial, highlight common misconceptions relating to PRS results, offer recommendations for best practice and discuss future challenges.

Cannabis use, depression and self-harm: phenotypic and genetic relationships.

BACKGROUND AND AIMS: The use of cannabis has previously been linked to both depression and self-harm; however, the role of genetics in this relationship is unclear. This study aimed to estimate the phenotypic and genetic associations between cannabis use and depression and self-harm. DESIGN: Cross-sectional data collected through UK Biobank were used to test the phenotypic association between cannabis use, depression and self-harm. UK Biobank genetic data were then combined with consortia genome-wide association study summary statistics to further test the genetic relationships between these traits using LD score regression, polygenic risk scoring and Mendelian randomization methods. SETTING: United Kingdom, with additional international consortia data. PARTICIPANTS: A total of 126 291 British adults aged between 40 and 70 years, recruited into UK Biobank. MEASUREMENTS: Phenotypic outcomes were life-time history of cannabis use (including initial and continued cannabis use), depression (including single-episode and recurrent depression) and self-harm. Genome-wide genetic data were used and assessment centre, batch and the first six principal components were included as key covariates when handling genetic data. FINDINGS: In UK Biobank, cannabis use is associated with an increased likelihood of depression [odds ratio (OR) = 1.64, 95% confidence interval (CI) = 1.59-1.70] and self-harm (OR = 2.85, 95% CI = 2.69-3.01). The strength of this phenotypic association is stronger when more severe trait definitions of cannabis use and depression are considered. Using consortia genome-wide summary statistics, significant genetic correlations are seen between cannabis use and depression [rg = 0.289, standard error (SE) = 0.036]. Polygenic risk scores for cannabis use and depression explain a small but significant proportion of variance in cannabis use, depression and self-harm within a UK Biobank target sample. However, two-sample Mendelian randomization analyses were not significant. CONCLUSIONS: Cannabis use appeared to be both phenotypically and genetically associated with depression and self-harm. Limitations in statistical power mean that conclusions could not be made on the direction of causality between these traits.

Classical Human Leukocyte Antigen Alleles and C4 Haplotypes Are Not Significantly Associated With Depression.

BACKGROUND: The prevalence of depression is higher in individuals with autoimmune diseases, but the mechanisms underlying the observed comorbidities are unknown. Shared genetic etiology is a plausible explanation for the overlap, and in this study we tested whether genetic variation in the major histocompatibility complex (MHC), which is associated with risk for autoimmune diseases, is also associated with risk for depression. METHODS: We fine-mapped the classical MHC (chr6: 29.6-33.1 Mb), imputing 216 human leukocyte antigen (HLA) alleles and 4 complement component 4 (C4) haplotypes in studies from the Psychiatric Genomics Consortium Major Depressive Disorder Working Group and the UK Biobank. The total sample size was 45,149 depression cases and 86,698 controls. We tested for association between depression status and imputed MHC variants, applying both a region-wide significance threshold (3.9 × 10-6) and a candidate threshold (1.6 × 10-4). RESULTS: No HLA alleles or C4 haplotypes were associated with depression at the region-wide threshold. HLA-B*08:01 was associated with modest protection for depression at the candidate threshold for testing in HLA genes in the meta-analysis (odds ratio = 0.98, 95% confidence interval = 0.97-0.99). CONCLUSIONS: We found no evidence that an increased risk for depression was conferred by HLA alleles, which play a major role in the genetic susceptibility to autoimmune diseases, or C4 haplotypes, which are strongly associated with schizophrenia. These results suggest that any HLA or C4 variants associated with depression either are rare or have very modest effect sizes.