We identified three low-frequency variants located at 6p22.2 (rs2298090 in HIST1H1E encoding p.Lys152Arg: OR = 0.30, P = 2.40 × 10(-16)) and 6p21.33 (rs200847762 in FKBPL encoding p.Pro137Leu: OR = 0.11, P = 3.77 × 10(-16); rs11754464 in MSH5: OR = 1.78, P = 3.71 × 10(-7)) associated with NOA risk after Bonferroni correction.