We first performed a GWAS in 833 POAG patients and 686 control subjects (hereafter referred to as ''Present GWAS''; Figure 1A ) who were selected and divided into case and control groups based on our strict diagnosis criteria. After genotyping, 653,519 SNPs passed the quality controls and were used for the association study ( Figure 1A) . In a quantile-quantile (Q-Q) plot ( Figure S1A ), the genomic inflation factor (l) showed 1.021, suggesting that the population substructure should not have any substantial effects on the association analysis. Under these conditions, we obtained 8 genome-wide significant SNPs (Table  S1 ) that passed the Bonferroni correction threshold (0.05/ 653,519 = 7.65610 28 ), of which 5 SNPs were located in CDKN2B-AS1 on chromosome 9p21.3 ( Figure 2 ). A few SNPs that showed modest to strong association were distributed across a single linkage disequilibrium (LD) block located in the locus ( Figure S2 ). The remaining genome-wide significant SNPs identified from different chromosomes turned out to be genotyping errors due to the poor 2D clusters ( Figure S3B-D) .