The total overlap across all association analyses comprised the genes WNT10A, HOXC13, DSP, LPAR6, ALX4, EDAR, CDH3, HR, and SPINK5. Notably, two of the top associated single variants (albeit not genome-wide significant), i.e., 2:218882368:C:A (Ptwo-as-control = 4.1 × 10−5) and 21:44499878:C:T (Ptwo-as-control = 9.0 × 10−6), which are located in WNT10A and TSPEAR respectively, were reported to be pathogenic for ectodermal dysplasia in previous studies28,29. This evidence concerns the gene HR and ectodermal dysplasia syndrome.