We conducted bioannotation and bioinformatic analyses (Supplementary Information section 6). The method FUMA (version 1.3.5e)26 was applied to explore the functional consequences of the 579 SNPs (Supplementary Table 9), which included ANNOVAR categories (i.e., the functional consequence of SNPs on genes), Combined Annotation Dependent Depletion (CADD) scores, RegulomeDB scores, expression quantitative trait loci (eQTLs), and chromatin states. The default external reference data for FUMA is described elsewhere26.
