Association of HK2 and NCK2 with Normal Tension Glaucoma in the Japanese Population

Although family studies and genome-wide association studies have shown that genetic factors play a role in glaucoma, it has been difficult to identify the specific genetic variants involved. We tested 669 single nucleotide polymorphisms (SNPs) from the region of chromosome 2 that includes the GLC1B glaucoma locus for association with primary open-angle glaucoma (POAG) and normal tension glaucoma (NTG) in the Japanese population. We performed a two-stage case-control study. The first cohort consisted of 123 POAG cases, 121 NTG cases and 120 controls: the second cohort consisted of 187 POAG cases, 286 NTG cases, and 271 controls. Out of six SNPs showing significant association with POAG in the first round screening, seven SNPs were tested in the second round. Rs678350 in the HK2 gene coding sequence showed significant allelic (p = 0.0027 in Stage Two, 2.7XE-4 in meta-analysis) association with POAG, and significant allelic (p = 4.7XE-4 in Stage Two, 1.0XE-5 in meta-analysis) association with NTG. Although alleles in the TMEM182 gene did not show significant association with glaucoma in the second round, subjects with the A/A allele in TMEM182 rs869833 showed worse visual field mean deviation (p = 0.01). Even though rs2033008 in the NCK2 gene coding sequence did not show significant association in the first round, it had previously shown association with NTG so it was tested for association with NTG in round 2 (p = 0.0053 in Stage Two). Immunohistochemistry showed that both HK2 and NCK2 are expressed in the retinal ganglion cell layer. Once multi-testing was taken into account, only HK2 showed significant association with POAG and NTG in Stage Two. Our data also support previous reports of NCK2 association with NTG, and raise questions about what role TMEM182 might play in phenotypic variability. Our data suggest that HK2 may play an important role in NTG in the Japanese population.


Introduction
Glaucoma is a complex, heterogeneous disease characterized by a progressive degeneration of the optic nerve fibers, and is the second highest cause of blindness worldwide affecting approximately 70 million people [1]. The most common type of openangle glaucoma, primary open-angle glaucoma (POAG), is associated with elevated intraocular pressure (IOP) [2], and another less-common subgroup of open-angle glaucoma, called low-tension glaucoma (LTG) [3,4] or normal tension glaucoma (NTG) is associated with IOP that does not rise outside of the normal range [5]. The prevalence of NTG is reported to be higher among the Japanese than among Caucasians [6,7]. This is an important medical and public health problem because simple screening programs based on detection of elevated IOP are not effective in a population where NTG is highly prevalent. Thus, an accurate diagnostic test for presymptomatic detection of individuals at risk for glaucoma, especially NTG in Japan, is urgently needed.
Open-angle glaucoma is a genetically heterogeneous disorder attributed to the interaction of multiple genes and environmental factors [8,9]. More than 15 POAG loci have been identified by linkage, and five open-angle glaucoma genes located within those loci have been identified [10,11]. More recently, genome-wide association studies (GWAS) using high-density single nucleotide polymorphism (SNP) arrays have been used to identify genetic risk factors involved in the common, complex forms of open-angle glaucoma that do not show classical Mendelian inheritance patterns. Burdon et al. identified susceptibility loci at TMCO1 and CDKN2B-AS1 that contribute to severe forms of glaucoma [12]. Ramdas et al. used meta-analysis of data from six separate studies to find significant evidence that three common variants of CDKN2B, ATOH7 and SIX1 are associated with POAG [13]. Wiggs et al. found significant evidence that genetic variants in CDKN2B-AS1 and a gene desert on 8q22 are associated with optic nerve damage in glaucoma [14].
Based on a linkage study involving 6 Caucasian families in the UK, the GLC1B locus for adult-onset open-angle glaucoma was identified at chromosome 2cen-q13 [15]. The patients in these families had clinical characteristics of low to moderate IOP, disease onset in their late 40 s, and a good response to medical therapy, and those phenotypes mimic the majority of Japanese NTG cases. Thus, the screening of the gene around GLC1B locus may be useful for diagnosis of POAG and NTG in the general Japanese population.
The purpose of this study was to screen for candidate genes for POAG and NTG on chromosome 2, around the GLC1B (glaucoma 1, open angle, B) locus in unrelated Japanese patients, using high density SNP scanning and case-control association.
Here we report one gene that shows significant association with POAG and NTG, support for a previously reported association with NTG, and a gene for which genotype is predictive of severity of mean deviation on the visual field test.

A Two-stage Case-control Study of SNPs on Chromosome2
To identify a gene associated with glaucoma we did a highdensity scan of the region around GLC1B on chromosome 2 by screening 669 SNPs on chromosome 2 in a two-stage case-control study design (Figure 1). We were especially interested in whether any SNPs that fall within genes in the GLC1B region might be associated with POAG or NTG (Table 1). Among genes from this region we found fourteen SNPs that show significant evidence of association with POAG, and nine SNPs that show significant evidence of association with NTG. Four of the SNP alleles which show significant evidence of association are identical between the POAG and NTG subjects.
We identified six SNPs in Stage One that showed evidence of association with POAG (rs1239066, rs1529385, rs869833, rs960011, rs1027003) and two SNPs (rs678350, rs2033008) that showed association with NTG ( Figure 2 and Table 2). The SNPs rs869833 and rs960011 are located within the TMEM182 gene, which contains 5 exons and 229 amino acids. The amino acid sequence of TMEM182 predicts an evolutionarily-conserved novel transmembrane protein, which consists of four putative membrane-spanning regions indicative of an integral membrane topology. The SNP rs678350 is located within the Hexokinase 2 (HK2) gene which contains 18 exons and 917 amino acids. The SNP rs678350 exists on intron1 of HK2 gene. The HK2 gene produces a protein product localizes on the outer membrane of mitochondria and plays an important role in intracellular glucose metabolism by catalyzing the conversion of glucose to glucose-6phosphate. There are no known genes closely neighboring the SNP rs1239066 and rs1027003. We have checked the SNPs near rs1529385 within LOC129293 and TMSB10 and did not find get any positive polymorphisms.
We selected three genes for second stage mutation screening. TMEM182 and HK2 genes were selected because they contain SNPs that showed significant evidence of association in Stage One in this study. Even though SNPs in the NCK2 gene showed significant evidence of association (p = 0.014), the NCK2 gene was selected for second stage mutation screening based of the previous report that it is associated with NTG [16].
The SNP rs1239066 showed significant evidence of association in both Stage One and Stage Two screenings. In the meta-analysis it showed significant evidence of association (POAG, NTG; P = 0.001, 0.005), as did rs1027003 (POAG; P = 0,010) ( Table 2).

HK2 Variants Detected in this Study
The SNP rs678350 in the HK2 gene coding sequence showed significant allelic (p = 0.0027 in Stage Two, 2.7XE-4 in metaanalysis) association with POAG, and significant allelic (p = 4.7XE-4 in Stage Two, 1.0XE-5 in meta-analysis) association Figure 1. Experimental Study Design. The study used a first stage of SNP screening in one population to identify SNPs and genes to be tested in a second population through SNP association testing and mutation screening of genes containing SNPs associated with glaucoma. Stage Two tested SNPs for association in a second population and did mutation screening in that second population in genes containing SNPs that showed significant association with POAG or NTG (p,0.01) in Stage One of our study or that had been previously reported to show significant evidence of association [16]. doi:10.1371/journal.pone.0054115.g001 with NTG ( Table 2). The rs678350 showed also a significant difference in genotype frequency (p = 0.0046 and 0.0039) in the POAG and NTG groups (Table 3). In the second round, we screened the HK2 coding sequence and intron-exon boundaries for mutations in POAG and NTG patients. After direct sequencing, we found 2 coding SNPs; p.Gln142His (A/T at the third nucleotide; rs2229621) in exon 4 and p.Arg844Lys (G/A at the second nucleotide; rs2229629) in exon 17. The allelic frequency of the p.Gln142His (A/T) variant was significantly higher in the NTG group than in the control group (p = 0.025), but it was not higher in the POAG group than in the control group (p = 0.181). The genotype frequency of the p.Gln142His (A/T) variant (dominant model) was significantly higher in the NTG group than in the control group (p = 0.019) but we did not find evidence that the frequency in the POAG group was different from the control group (p = 0.179). There was no evidence of a significant difference between POAG and NTG for p.Arg844Lys. No other mutation was found. We tested the LD block and found no linkage disequilibrium between SNPs rs678350 and rs2229621 (D9 = 0.08). We tested the correlation between the phenotypes POAG or NTG and the genotypes screened in the second stage screening of HK2, and found no association with any phenotypes including age at diagnosis, maximum IOP under medication, and MD value of the visual field ( Table 4). None of the polymorphisms showed deviation from Hardy-Weinberg equilibrium (P,0.05).

NCK2 Variants Detected in this Study
The SNP rs2033008 in the NCK2 gene showed a significant difference in allelic frequency (p = 0.015 in Stage One, p = 0.0053 in Stage Two, and 2.2XE-4 in meta-analysis) between controls and NTG, but not between control and POAG status (p = 0.147 in Stage One, 0.35 in Stage Two, 0.12 in meta-analysis) ( Table 2 and  Table 5). The odds ratio for association with NTG supports a model in which NCK2 is associated with NTG; OR = 0.69 (0.53-0.89), but the odds ratio for association with POAG was not significant; 0.87 (0.65-1.16) in Stage Two. Only this rs2033008 polymorphism in Stage One showed deviation from Hardy-Weinberg equilibrium (p = 0.030).
We screened the sequence of the NCK2 coding sequence and intron-exon boundaries for mutations in POAG and NTG patients and found 1 synonymous coding base change: Thr14Thr (ACC .ACT) in one NTG subject. Although SNP rs2033008 showed significant association but the Thr14Thr variant showed no statistical difference in allele frequency between NTG and normal subjects (p = 0.33). The Thr14Thr heterozygotes (A/A) in the NTG subjects have the worse Mean Deviation value of the visual field compared with human reference sequence (T/T) (p = 0.05) ( Table 6).
The single NTG subject with the Thr14Thr variant was a 50 year old woman whose father also had NTG. Her age at diagnosis was 40 years old. Her initial IOP was 15 mmHg in each eye. The mean deviation (MD) of the visual field test was 23.87 dB in the right eye, and 21.97 dB in the left eye, reflecting a mild NTG phenotype.

TMEM182 Variants Detected in this Study
We found no mutations in TMEM182 coding sequence and intron-exon boundaries, for mutations in POAG and NTG patients. After finding association for SNPs rs869833 and rs960011 in the POAG and NTG subjects in Stage One, we did not find this association confirmed in Stage Two SNP testing (Table 7). When we tested for correlation between Stage Two case endophenotypes and TMEM182 genotype, we found association with the MD value of the visual field in POAG subjects. The POAG subjects homozygous for the A/A allele of SNP rs869833 have worse Mean Deviation value of the visual field test compared with those who carry the G/G genotype in Stage Two subjects (p = 0.01) ( Table 8).

Immunohistochemistory of the HK2 and NCK2
Representative immunohistochemistory (IHC) photographs with Hk2, Nck2, astrocyte marker (GFAP) and retinal ganglion cell marker (C38) on the retinas of untreated mice were shown (Fig. 3, A, B). Hk2 and Nck2 were strongly immunoreactive in the ganglion cell layer. C38 signals co-localized with Hk2 or Nck2 in the ganglion cell layer, as indicated by arrows (Fig. 3 A, B). Hk2 expression is only located in the ganglion cell layer. Nck is expressed in the ganglion cell layer, inner nuclear layer, and outer plexiform layer, with the highest level of expression in the ganglion cell layer.

HK2 in POAG and NTG subjects
The rs678350 in the HK2 gene coding sequence showed significant allelic (p = 0.043 in Stage One, p = 0.0027 in Stage  Two, 2.7XE-4 in meta-analysis association with POAG, and significant allelic (p = 0.004 in Stage One, p = 4.7XE-4 in Stage Two, 1.0XE-5 in meta-analysis) association with NTG ( Table 2). The rs678350 polymorphism showed a significant case-control difference in genotype frequency (p = 0.0046 and 0.0039) in the POAG and NTG groups (Table 3). However, there was no association of this SNP with glaucoma endophenotypes including age at diagnosis, maximum intra ocular pressure under medication, and MD value of the visual field. So the HK2 gene may contribute to disease susceptibility to POAG and NTG, but may not account for all of the phenotypic variability between individuals whose glaucoma results from variants in this gene.
Our association findings suggest that the HK2 gene that contains this polymorphism might play a role in POAG and NTG in the Japanese population, but it remains to be seen whether rs678350 is actually causative, perhaps through altering transcription or splicing, or whether another allele(s) in this gene or its regulatory region might actually be causing the disease. There remains a possibility that the p.Gln142His (A/T) SNP in HK2 may play a role in disease pathology, but our study can only show association, not causation. Because none of our subjects come from the original families used to map the GLC1B locus, we can only draw conclusions regarding the possible role of this gene in the Japanese population, but this finding raises questions about whether this could be the GLC1B gene. The HK2 gene product plays an important role in intracellular glucose metabolism by catalyzing the conversion of glucose to glucose-6-phosphate. The HK2 gene localizes to the outer membrane of mitochondria. Since reduced glucose-6-phosphate content in muscle has been demonstrated in pre-non-insulindependent diabetes mellitus (pre-NIDDM) and NIDDM subjects, HK2 was investigated as a promising candidate gene for noninsulin-dependent diabetes mellitus (NIDDM; OMIM125853) [17,18]; however, those studies concluded that mutations of the HK2 gene, including a common p.Gln142His polymorphism is not a major etiologic factor for NIDDM in the Finnish [17,18,19,20], British [19], and Danish [20] populations. In brain, mitochondrial-hexokinase activity plays a key antioxidant role protecting against oxidative stress (ROS) [21], and complements the classical antioxidant enzymes that protect against oxidative stress [22]. Hexokinase antagonizes the release of mitochondrial cytochrome C activation of Akt, which is recognized as a potent inhibitor of apoptosis. HK2 is probably associated with an anti-oxidative reaction and inhibition of apoptosis through Bax/Bak-mediated cytochrome c release [23]. Leber's hereditary optic neuropathy (LHON) -associated mitochondrial DNA mutations were found in Japanese patients with POAG [24], so it is reasonable to consider a gene whose product plays a role in mitochondria as a candidate gene for other phenotypes involving optic neuropathy.

NCK2 genes in POAG and NTG subjects
The NCK2 gene, which was previously reported to be associated with NTG [16], encodes a member of the NCK family of adaptor proteins, and the adaptor protein which associates with tyrosinephosphorylated growth factor receptors of their cellular substrates. SH2/SH3 domain-containing adapter proteins, such as the NCK family, play a major role in regulating tyrosine kinase signaling [25]. Previously, microsatellite marker D2S176 within the GLC1B locus showed significant association with NTG in the Japanese population, and D2S176 is located 24 kb from the NCK2 gene [16]. Brain-derived neurotrophic factor (BDNF) binds to and activates the TrkB tyrosine kinase receptor to regulate cell differentiation and survival in the nervous system. BDNF Table 3. Stage Two HK2 SNPs Allele Frequencies in Japanese POAG, NTG and Control Subjects.  stimulation promotes interaction of Nck2 with TrkB in cortical neuron [26]. And BDNF signaling in glia is known to play important roles in neural protection and regeneration, particularly in conversion of Muller glia to photoreceptors [27]. In our study, it is interesting that the NCK2 variant rs2033008 showed a significantly difference from the control population in the NTG group, where the disease pathology seems to be focused on the retinal ganglion cells and the optic nerve, but not in the POAG group, where a substantial disease component localizes to the anterior chamber of the eye (Table 5). Thus although our study falls short of achieving a level of significance needed to identify NCK2 de novo as a glaucoma gene and this polymorphism showed deviation from Hardy-Weinberg equilibrium (p = 0.030 in Stage One), our data do support the prior finding of significant allele frequency differences between NTG cases and normal controls in the Japanese population [16]. It is unclear whether this deviation from Hardy-Weinberg equilibrium in Stage One might represent the absence of some alleles from this population because they are associated with diagnoses specifically excluded from this study, such as ocular hypertension.

Immunohistochemistory of the Hk2 and Nck2
Although KH2 and NCK2 had previously been detected in retina, more precise localization to specific cell types is needed to begin understanding how the gene products might play a role in disease pathology. In the representative IHC photographs, antibody against Hk2 was strongly immunoreactive in the ganglion cell layer (GCL). The Hk2 protein localizes to the outer membrane of mitochondria, and interestingly the Hk2 protein appears in the GCL. Nck2, which interacts with BDNF, is expressed in ganglion cell layer (GCL), inner nuclear layer (INL) and outer plexiform layer (OPL), and most expressed in GCL. This localization makes it highly conceivable that the Hk2 gene products could each play a role in glaucoma, and there is possibility that Nck2 could have relationship with glaucoma.

TMEM182 in POAG and NTG
SNPs rs869833 and rs960011 in the TMEM182 gene showed significant association with POAG and NTG in Stage One that was not confirmed in the second stage. The primary amino acid sequence of TMEM182 predicts an evolutionarily-conserved novel transmembrane protein, which consists of four putative membrane-spanning regions indicative of an integral membrane topology. The TMEM182 protein sequence lacks homologies with previously-defined protein families. However, the proinflammatory cytokine TNFa down-regulates TMEM182 transcript expression in adipocytes [28]. Its transcript expresses in white adipose tissues, heart, muscle, and lower relative levels of TMEM182 transcript are found in kidney, testis, and brain. Identification of the intracellular signaling pathway involved in the TNFa-mediated decrease might be one clue offering insights into association between POAG and TMEM182 function. Failure to confirm the association with POAG in the second stage could be attributable to clinical heterogeneity, but this result still needs to be confirmed in a second population. Nakano et al. demonstrated heterogeneity in the Japanese POAG population when their genome-wide association study of 1,575 Japanese POAG and normal subjects identified significant evidence of association with 6 SNPs on Chromosome 1, 10, 12 [29]. They did not report evidence of association with SNPs on Chromosome 2. Thus these SNPs on Chr2 might be the variants for which our study is not well-powered, or clinical heterogeneity might be complicating our ability to detect the association in our limited sample size. In our data set, HK2 shows the strongest evidence of association with NTG in the Japanese population out of all of the genes that have SNPs represented on the screening panel we used.
Even when a single simple Mendelian locus causes a disease, variants in other genes may contribute to phenotypic variability, and phenotypic complexity along with locus and allele heterogeneity can complicate the problem of identifying the underlying causes of the disease. Our findings raise questions about whether additional genes in this region may be contributing to phenotypic heterogeneity within the NTG and POAG populations. The MD values of the visual field in these studies indicate middle to advanced stages of the disease, with the range of values possibly resulting from a combination of genetic complexity and genetic heterogeneity.
On the other hand, we also have to consider the importance of apparent gene deserts. The SNP rs1239066 shows significant evidence of association in both Stage One and Stage Two screenings, and in the meta-analysis (POAG, NTG; P = 0.001, 0.005). Significant evidence was also found for rs1027003 (POAG; P = 0,010) ( Table 2), but neither one is in the immediate vicinity of a known gene. Wiggs et al. found significant evidence that genetic variants in a gene desert on 8q22 are associated with optic nerve damage in glaucoma [14]. So additional follow up studies will need to explore this gene desert region to determine whether any functional sequences there are playing a role in the disease.
Glaucoma is a complex disease, and it involves genetic variants that confer moderate to low effect sizes (e.g., OR = 1.2-1.5). The OR which was identified in the first stage with the P value cutoff of 0.01 was 1.68. This cutoff value was a bit strict to exclude the false positives.
Further investigations of the structure and function of the HK2, NCK2 and TMEM182 proteins would be helpful in understanding the pathogenesis of POAG and NTG. Our data suggest that HK2 may play an important role in NTG in the Japanese population; although our data suggest that HK2 might be the GLC1B gene, a firm conclusion on the subject awaits screening of members of the families originally used to map the GLC1B locus.

Ethics statement
This study was approved by the Institutional Review Board of Tohoku University, Keio University, Tokyo Metropolitan Police Hospital, Niigata University, Ideta Eye Hospital, and all procedures were conducted in accordance with the Declaration of Helsinki. All participants provided written informed consent after an explanation of the purpose and procedures to publish these case details.

Patient Recruitment and Characteristics
The samples used in the first screening (Stage One) were collected in Keio University hospitals, Tokyo Metropolitan Police Hospital, Niigata University, and Ideta Eye Hospital, and the Table 6. Correlation between the POAG or NTG Endophenotypes and NCK2 SNP rs2033008 Screened in Stage Two.

Endophenotype
Age at diagnosis (y.o.) Maximum IOP* (mmHg) MD value of the visual field (dB) samples used in the second screening (Stage Two) were collected at Tohoku University. Routine ophthalmic examinations were performed on all subjects. Individuals were included as POAG cases if they fulfilled the following inclusion criteria: 1) applanation IOP greater than 22 mm Hg in each eye; 2) spherical equivalent more than 28 diopter; 3) glaucomatous cupping in each eye including cup-to-disc ratio greater than 0.7; 4) visual field loss measured by Goldmann perimetry or Humphrey automated field analyzer (Carl Zeiss Meditec, Dublin, CA) in Stage One. The severity of the visual field defects was scored from 1 to 5 according to previously reported criteria ref. The data obtained by two types of perimetry were combined using a five-point scale: 1, no alterations; 2, early defects; 3, moderate defects; 4, severe defects; and 5, light perception only or no light perception. The first four groups on this severity scale followed Kozaki's classification based on Goldmann perimetry or the classification was based on results of Humphrey automated field analyzer [30,31]. Kozaki's classification is widely used in Japan. In Stage Two, all of the visual field loss were measured by Humphrey automated field analyzer according to Anderson-Pattela classification [30] consistent with the glaucomatous cupping in at least one eye; and 5) open anterior chamber angles; and exclusion of secondary causes (e.g., trauma, uveitis, or steroid-induced glaucoma). The criteria for NTG were the same as for POAG except that NTG subjects showed applanation IOP less than 22 mm Hg in both eyes at each examination. Baseline clinical parameters including age, gender, spherical equivalent base line visual acuity (VA), IOP measured by Goldmann applanation tonometry were recorded at the time of first diagnosis of POAG or NTG in each patient. Mean deviation (MD) values indicative of visual field damage were obtained by the Swedish interactive threshold algorithm (SITA)-standard strategy of the 30-2 program of HFA (Carl Zeiss Meditec, Dublin, Californuia, USA). MD was used on reliable visual field test results (,20% fixation errors, ,33% false-positive results, and ,33% false-negative results). Control subjects had these characteristics: IOP less than 22 mm Hg, normal optic discs, and no family history of glaucoma. To decrease the chance of enrolling individuals with pre-symp-  tomatic glaucoma, we limited this group to individuals older than 60 years.

Two Stage Screening Protocol
The two-stage screening strategy is shown in Figure 1. The first stage screened used 669 SNPs from the GLC1B region, using the Illumina panel #8 (Chromosome 2; 77-117 Mbp) (Illumina, San Diego, CA, USA) carried by Illumina (San Diego, CA, USA), and each SNP was tested for association with POAG or NTG. The first stage used a cohort of 368 unrelated Japanese individuals: 123 POAG cases (63 men and 60 women), 121 NTG cases (61 men and 60 women) and 120 normal subjects (61 men and 59 women). Mean age of POAG cases was 56.9611.5 years. Mean age of NTG cases was 54.0612.3 years. Mean age of controls was 70.3610.2 years. The visual field scores were 2.861.0 in POAG cases and 2.760.9 in NTG cases ( Table 9). Single-nucleotide polymorphisms (SNPs) with a call rate ,90% were excluded from the analysis. All of the polymorphisms showed no deviation from Hardy-Weinberg equilibrium (P.0.05) except rs2033008 in Stage One (p = 0.030).
The second stage screened SNPs that showed significant evidence for association in the first round in this study (p,0.01) and in the NCK2 gene, which previously showed association with NTG [16]. The odds ratio which was identified in the first stage with the P value cutoff of 0  (Table 9).

Sample Preparation
Genomic DNA was extracted from leukocytes of the peripheral blood. It was purified by the Qiagen QIAamp Blood Kit (Qiagen, Valencia, CA, USA).

Mutation Screening
Mutation screening was carried out in genes that contained SNPs that showed significant evidence of association in the first stage (HK2 and TMEM182) plus the previously-reported NCK2 gene [16]. All of the exons of the HK2, NCK2, TMSB10 and TMEM182 genes, and positive SNPs were amplified by a polymerase chain reaction (PCR) using 0.5 mM concentration of primers in an amplification mixture (25 ml) containing 0.2 mM dNTPs and 0.5 U Ex Taq polymerase (Takara Bio, Shiga, Japan) with 30 ng template DNA. Oligonucleotides for amplification and sequencing were selected using Primer3 software (http://frodo.wi. mit.edu/cgi-bin/primer3/primer3_www.cgi/ provided in the public domain by the Massachusetts Institute of Technology, Cambridge, MA). Primers for amplification and sequencing of coding sequence were placed in introns far enough from the intron/exon junctions to allow for visualization of the splice site sequences. The PCR fragments were purified with ExoSAP-IT (USB, Cleveland, Ohio, USA), sequenced by the BigDye TM Terminator v3.1 Cycle Sequencing Kit (Perkin-Elmer, Foster City, CA, USA) on an automated DNA sequencer (ABI PRISM TM 3100 Genetic Analyzer, Perkin-Elmer, Waltham, MA, USA).

Statistical Analysis
The significance of association was determined by contingency table analysis using Fisher's exact test or Chi-square test, depending on cell counts. In estimation of genotype-phenotype correlation, we used Dunnett's test to compare group means of those carrying the mutant variant being tested against the group means of those carrying the normal, reference sequence. Odds ratios (approximating to relative risk) were calculated as a measure of the association between the allele frequency and the phenotype of POAG/NTG, estimated using the SNPAlyze program version 7.0 (Dynacom, Yokohama, Japan). Hardy-Weinberg equilibrium was analyzed using gene frequencies obtained by simple gene counting and the chi-square test with Yates' correction for comparing observed and expected values.

Immunohistochemistory of the HK2 and NCK2
Murine retinas were fixed with 4% PFA at 4uC overnight and then cryoprotected in phosphate buffered saline (PBS) with 20% sucrose. Cryosections (thickness 10 mm) were mounted on the slides and incubated with blocking buffer (10% goat serum, 0.5% gelatin, 3% BSA and 0.2% Tween 20 in PBS). Next, they were incubated with primary antibodies at 4uC overnight. Primary Table 9. Clinical Characteristics of Subjects studied in Stage One and Two Screening.