PhenoScanner (36) was used to identify variants that are associated with other traits. All proxy SNVs in high LD (r2 ≥ 0.8) with the lead SNVs at our five unreported novel loci were investigated in the PhenoScanner 1000 Genomes reference dataset. Results were filtered to those reaching a genome-wide significance P-value ≤ 5 × 10−8.