Statistical analyses were performed using PLINK (Purcell et al. 2007) and R software. A logistic regression model was fit for the binary phenotype, SI-1, adjusting for age, sex, and geographic area under additive, dominant, and recessive models. The cumulative logit model was fit to the ordinal phenotype SI-2 and ND (Bender and Grouven 1997).