To estimate a rg for each pair of phenotypes, we used linkage disequilibrium score regression (LDSR)27 and local analysis of covariant annotation (LAVA).28 LDSR is a method for estimating a rg at a genome-wide level. LAVA estimates rgs at a “local” level within 2,495 genomic regions. We used the default heritability thresholds for LAVA (p=0.05). The Benjamini-Hochberg correction (q<0.05) was applied.