Genome-wide association analyses of risk tolerance and risky behaviors in over 1 million individuals identify hundreds of loci and shared genetic influences.

Manhattan plots.In all panels, the x-axis is chromosomal position; the y-axis is the GWAS P value on a −log10 scale (based on a two-tailed z-test); each lead SNP is marked by a red “×”; each conditional association is marked by a red “o”; and each SNP that is both a lead SNP and a conditional association is marked by a red “⊗”. a, Manhattan plots for the discovery GWAS of general risk tolerance (n = 939,908). b, Local Manhattan plots of a long-range LD region on chromosome 3 and a candidate inversion on chromosome 18 that contain lead SNPs for all seven of our GWAS. The gray background marks the locations of long-range LD or candidate inversion regions. c, Local Manhattan plots of the areas around the 15 most commonly tested candidate genes in the prior literature on the genetics of risk tolerance. Each local plot shows all SNPs within 500 kb of the gene’s borders that are in weak LD (r2 > 0.1) with a SNP in the gene. The 15 plots are concatenated and shown together in the panel, divided by the black vertical lines. The 15 genes are not particularly strongly associated with general risk tolerance or the risky behaviors, as can be seen by comparing the results within each row across panels b and c (the three rows correspond to the GWAS of general risk tolerance, adventurousness (n = 557,923), and the first PC of the four risky behaviors (n = 315,894)).

Genetic correlations with general risk tolerance.The genetic correlations were estimated using bivariate LD Score (LDSC) regression9. Error bars show 95% confidence intervals. For the supplementary phenotypes and the additional risky behaviors, green bars represent significant estimates with the expected signs, where higher risk tolerance is associated with riskier behavior. For the other phenotypes, blue bars represent significant estimates. Light green and light blue bars represent genetic correlations that are statistically significant at the 5% level, and dark green and dark blue bars represent correlations that are statistically significant after Bonferroni correction for 35 tests (the total number of phenotypes tested). Grey bars represent correlations that are not statistically significant at the 5% level. The two dotted vertical lines indicate genetic correlations of −0.5 and 0.5, respectively. All significance tests are two-sided.

Results from selected biological analyses.a, DEPICT gene-set enrichment diagram. We identified 93 reconstituted gene sets that are significantly enriched (FDR < 0.01) for genes overlapping DEPICT-defined loci associated with general risk tolerance; using the Affinity Propagation method43, these were grouped into the 13 clusters displayed in the graph. Each cluster was named after its exemplary gene set, as chosen by the Affinity Propagation tool, and each cluster’s color represents the permutation P value of its most significant gene set. The “synapse part” cluster includes the gene set “glutamate receptor activity,” and several members of the “GABAA receptor activation” cluster are defined by gamma-aminobutyric acid signaling. Overlap between the named representatives of two clusters is represented by an edge. Edge width represents the Pearson correlation ρ between the two respective vectors of gene membership scores (ρ < 0.3, no edge; 0.3 ≤ ρ < 0.5, thin edge; 0.5 ≤ ρ < 0.7, intermediate edge; ρ ≥ 0.7, thick edge). b, Results of DEPICT tissue enrichment analysis using GTEx data. The panel shows whether the genes overlapping DEPICT-defined loci associated with general risk tolerance are significantly overexpressed (relative to genes in random sets of loci matched by gene density) in various tissues. Tissues are grouped by organ or tissue type. The orange bars correspond to tissues with significant overexpression (FDR < 0.01). The y-axis is the significance on a −log10 scale. See Supplementary Note for additional details.