We used Functional Mapping and Annotation of Genome-Wide Association Studies (FUMA) SNP2GENE web tool (Watanabe et al. 2017), to perform annotation of the male- and female-specific genome-wide associations. We used the default settings with minor modifications: the minimum p-value of the lead SNP was set at 1.0E-5, the r2 threshold to define the LD structure of the lead SNPs greater or equal to 0.6, the maximum p-value cutoff at 0.5, MAF>0.01, 250 kb as the maximum distance between LD blocks to merge into a locus, 1000 Genomes Project Phase 3 European population as the reference panel, and variants from the reference panel were included for identification of functional variants in LD with the lead SNP. For the gene-based analysis implemented via generalized-gene set analysis of GWAS data (MAGMA) (de Leeuw et al. 2015) in FUMA, SNPs are mapped to genes if they fall within the gene start and end sites. In this analysis the mean of the chi2-statistic for the SNPs in a gene is calculated and the p-value is obtained from a known approximation of the sampling distribution (Brown
