latent growth curve models (one for each SNP) were fitted, where the same set of covariates, intervention condition, SNP genotype, as well as their interaction term were included to predict the intercept and linear slope. A modified Bonferroni correction was used to adjust for multiple testing involving 10 SNPs (Nyholt 2004). After accounting for the linkage structure (correlations among SNPs), two sets of SNP coding schemes (additive and dominance), and separate testing for AAs and EAs, the adjusted significance criterion was p = 0.0014 for EAs and p = 0.00135 for AAs, respectively.