Abstract

Primary angle-closure glaucoma (PACG) is a common form of glaucoma in the Far East. Its defining feature is iridocorneal angle closure. In addition to PACG, indications of angle closure are included in the diagnostic criteria of related conditions primary angle-closure suspect (PACS) and primary angle closure (PAC). To the best of our knowledge, a causative gene for iridocorneal angle closure in humans has not been identified. This study aimed to identify the genetic cause of iridocorneal angle closure in a pedigree with at least 10 individuals diagnosed with PACS, PAC or PACG. Results of linkage analysis, segregation analysis of 44 novel variations, whole exome sequencing of 10 individuals, screenings of controls and bioinformatics predictions identified a mutation in COL18A1 that encodes collagen type XVIII as the most likely cause of angle closure in the pedigree. The role of COL18A1 in the etiology of Knobloch syndrome (KS) that is consistently accompanied by optic anomalies, available functional data on the encoded protein and the recognized role of collagens and the extracellular matrix in glaucoma pathogenesis supported the proposed role of the COL18A1 mutation in the pedigree. Subsequent identification of other COL18A1 mutations in PACS affected individuals of two unrelated families further supported that COL18A1 may affect angle closure. These PACS individuals were parents and grandparents of KS-affected children. In conclusion, a gene that affects angle closure in humans, a critical feature of PACG, has been identified. The findings also reinforce the importance of collagens in eye features and functions.

Introduction

Glaucoma is a heterogeneous group of optic neuropathies characterized by a specific pattern of optic nerve degeneration and visual field loss that is usually accompanied by increased intraocular pressure (IOP) (1). It is the major cause of irreversible blindness worldwide (2). Primary glaucoma is classified into three major forms based on the anatomy of the anterior chamber drainage angle of the eye and age of onset: primary congenital glaucoma (PCG), primary open-angle glaucoma (POAG) and primary angle-closure glaucoma (PACG) (1). In glaucoma patients with increased IOP, the increase is thought to be due mainly to impaired drainage of aqueous humor from the anterior chamber (3). In PCG, the most severe and the least-common form of glaucoma, an anatomical defect (trabeculodysgenesis) in the trabecular meshwork (TM) is thought to affect poor drainage (4). As the angle between the iris root and the posterior TM remains open in POAG patients, increased IOP in these individuals is thought to be due to TM anomalies that hinder the flow of aqueous humor. In contrast, eyes with PACG reveal a closed-anterior chamber angle due to contact between the peripheral iris and the posterior TM. This closure obstructs outflow of the aqueous humor. With regards to genetics, PCG is usually a monogenetic Mendelian disease for which three causative genes, cytochrome p4501B1 (CYP1B1) (OMIM 601771), latent transforming growth factor beta binding protein 2 (LTBP2) (OMIM 602091) and TEK receptor tyrosine kinase (TEK) (OMIM 600221), have been identified (58). POAG and PACG are generally considered complex multifactorial disorders (9). Several POAG-causing genes, including myocilin (MYOC) (OMIM 601652), optineurin (OPTN) (OMIM 602432), WD repeat containing protein 36 (WDR36) (OMIM 609669), neurotrophin-4 (NTF4) (OMIM 162662) and TANK-binding kinase 1 (TBK1) (OMIM 604834) have been identified (1014), but mutations in these account for disease in <10% of patients. Contrary to PCG and POAG, a definitive PACG causing gene has not been identified in humans. However, high prevalence of this condition in Asians, observations on familial clustering and results of sib-pair studies suggest that genetic factors are involved in the etiology of PACG (1518). Association studies have implicated possible roles for several genes including collagen type XI alpha 1 (COL11A1) (OMIM 120280), pleckstrin homology domain containing, family A member 7 (PLEKHA7) (OMIM 612686), ATP binding cassette subfamily C member 5 (ABCC5) (OMIM 605251), protein-L-isoaspartate (D-aspartate) O-methyltransferase domain containing (PCMTD1)–C2H2C-type zinc finger (ST18) (OMIM 617155), ependymin related 1 (EPDR1), choline O-acetyltransferase (CHAT) (OMIM 118490), GLIS family zinc finger 3 (GLIS3) (OMIM 610192), fermitin family member 2 (FERMT2) (OMIM 607746) and dolichyl-phosphate mannosyltransferase subunit 2, regulatory (DPM2) (OMIM 603564)– family with sequence similarity 102 member A (FAM102A) (OMIM 610891) encoding genes (1921). Results of candidate gene screenings have also suggested possible roles for some genes including genes that code for matrix metalloproteinase-9 (MMP-9) (OMIM 120361), hepatocyte growth factor (HGF) (OMIM 142409), methylene tetrahydrofolate reductase (MTHFR) (OMIM 607093), nitric oxide synthase 3 (eNOS) (OMIM 163729) and heat shock protein family A member 1A (HSP70-1A) (OMIM 140550) (2226). However, these findings are not generally considered definitive, and putative roles of genes suggested by some studies were not confirmed in independent studies (2729). Taken together, a PACG-causing gene or a strongly associated locus has not yet been identified (17).

Failure to find a strong PACG-associated locus likely reflects the complexity of the disease, meaning that multiple genes that affect various anatomical and functional features that associate with PACG contribute to its etiology. Factors known to be associated with PACG include hyperoptic refractive error, shallow anterior chamber, thick crystalline lens, short axial length, small corneal diameter and narrow iridocorneal angle (3032). Clearly, some of these also associate with each other. Ultimately, the defining feature of PACG in individuals with glaucomatous optic nerve damage is angle closure. In addition to PACG, indications of angle closure are included in the diagnostic criteria of conditions known as primary angle-closure suspect (PACS) and primary angle closure (PAC) (see Methods section) (33). Angle closure evidences in these conditions in the absence of glaucomatous optic neuropathy.

The focus of this study is genetic analysis of a pedigree-labeled ANG200 in which angle closure is prevalent. The individuals were diagnosed with PACS, PAC or PACG. We report the presence of a mutation in COL18A1 in all the individuals diagnosed with these conditions. We found supporting evidence for the role of COL18A1 in pathogenesis of angle closure by observation of COL18A1 mutations in PACS-affected individuals of two additional unrelated families.

Results

Family and clinical data

The ANG200 pedigree with at least 1 members who presented with varying degrees of iridocorneal angle closure is shown in Figure 1. Their phenotypic data are presented in Table 1. Certain features of pedigree ANG200 are to be noted. First, nearly all individuals in generation III had been or were affected with cataract and/or a degree of angle closure. The high incidence of cataract in these individuals may be related to the fact that they were elderly, aged 75 years or older. Cataract was not observed in individuals of generation IV whose ages range from 41 to 69 years. The angle-closure status of many individuals in the pedigree was unknown to them until eye examinations were performed after they were recruited into the study. Age of onset of the condition is unknown in all affected individuals. Rather unexpectedly, many individuals in generations III and IV, including ANG200 IV-5–IV-9, have remained unwed.

ANG200 pedigree. Black-filled shapes: individuals with various degrees of ridocorneal angle closure and diagnosed with PACS, PAC or PACG (see Table 1); grey-filled shape: glaucoma diagnosis based on presence of glaucomatous optic nerve damage, but angle closure not assessable due to prior intracapsular cataract extraction eyes resulting from surgery; hatched shapes: reported by family members to have been blind in one or both eyes at time of death; open shapes of generations I–IV: reported by family members not to have had ocular problems; open shapes of generation V: normal eye exam. Arrow, proband; *, individuals whose whole genome SNP data were included in the linkage analysis, ° individuals for whom exome sequencing was performed. Clinical diagnoses and age at time of examination (years) are indicated; NEE, normal eye examination. Genotypes of individuals ascertained by direct sequencing are presented: N, normal COL18A1 allele; M, mutated COL18A1 allele.
Figure 1

ANG200 pedigree. Black-filled shapes: individuals with various degrees of ridocorneal angle closure and diagnosed with PACS, PAC or PACG (see Table 1); grey-filled shape: glaucoma diagnosis based on presence of glaucomatous optic nerve damage, but angle closure not assessable due to prior intracapsular cataract extraction eyes resulting from surgery; hatched shapes: reported by family members to have been blind in one or both eyes at time of death; open shapes of generations I–IV: reported by family members not to have had ocular problems; open shapes of generation V: normal eye exam. Arrow, proband; *, individuals whose whole genome SNP data were included in the linkage analysis, ° individuals for whom exome sequencing was performed. Clinical diagnoses and age at time of examination (years) are indicated; NEE, normal eye examination. Genotypes of individuals ascertained by direct sequencing are presented: N, normal COL18A1 allele; M, mutated COL18A1 allele.

Table 1

Phenotypic and genotypic data on pedigree ANG200 individuals with mutation in COL18A1

Individual IDIII-3III-4III-5III-6IV-2IV-3IV-5IV-6IV-7IV-8IV-9V-1V-2V-3&
COL18A1 genotypeNM?M#NMNMNMNMNMNMMMNMNMNMNMNN
Clinical examination
SexFMMFMFFMFMFFMF
Age at examination (yrs)82777569665649534841373540
Angle1
circumference OD/OS*360°/VLE270°/ 270°360°/ 360°360°/ 360°270°/ 270°270°/ 270°360°/ 360°360°/ 360°360°/ 360°270°/ 270°360°/360°360°/360°360°/360°
GradeOD/OS*0-1/VLE1/ 10-1/ 0-10-1/ 0-11 /11 /10-1/ 0-10-1/ 0-10-1/ 0-11 /13/33/33/3
Cornea1OD/OSEdema/EdemaClear/ClearClear/ ClearClear/ ClearClear/ ClearClear/ ClearClear/ ClearClear/ ClearClear/ ClearClear/ ClearClear/ ClearClear/ClearClear/ClearClear/Clear
Anterior chamber1OD/OSDeep/DeepShallow & clear/ Shallow & clearShallow & clear/ Shallow & clearShallow & clear/ Shallow & clearShallow & clear/ Shallow & clearShallow & clear/ Shallow & clearShallow & clear/ Shallow & clearShallow & clear/ Shallow & clearShallow & clear/ Shallow & clearShallow & clear/ Shallow & clearShallow & clear/ Shallow clearDeep/DeepDeep/DeepDeep/Deep
Lens1OD/OSAphakic/ AphakicSevere cataract/
VLE
Severe cataract/
moderate cataract
Severe cataract (OU)Clear/ ClearClear/ ClearClear/ ClearClear/ ClearClear/ ClearClear/ ClearClear/ ClearClear/ClearClear/ClearClear/Clear
Vertical C/D ratio1OD/OS0.80/ 0.800.40/VLE0.25/ 0.250.45/ 0.450.55/ 0.600.40/ 0.400.50/ 0.500.50/ 0.500.30/ 0.300.40/ 0.450.60/ 0.600.25/0.250.50/0.600.30/0.30
IOP2(mmHg)OD/OS144/14425/VLE26/2812/1218/1818/1916/1617/1729/3023/2314/1712/1211/13
RefractionOD/OS+13.75/ no reflex-1.25 -1.25x65/ VLE+2.25 -0.75x135/ +3.50 -0.25x60Plano -2.75/ +1.50+1.50 -0.75x90/ +1.75-0.50x85+0.75 -1.00x80/ +0.25 -1.00x180Plano -.25/ +0.25 -0.50x50Plano -0.50x3/ plano -0.25x180+1.25 -0.50x85/ +1.25 -0.25x85-0.75 -0.50x62/ -0.75 -1.00x95Plano/ plano -0.75x170Plano/PlanoPlano/Plano-3.50/ -3.50
Individual IDIII-3III-4III-5III-6IV-2IV-3IV-5IV-6IV-7IV-8IV-9V-1V-2V-3&
COL18A1 genotypeNM?M#NMNMNMNMNMNMMMNMNMNMNMNN
Clinical examination
SexFMMFMFFMFMFFMF
Age at examination (yrs)82777569665649534841373540
Angle1
circumference OD/OS*360°/VLE270°/ 270°360°/ 360°360°/ 360°270°/ 270°270°/ 270°360°/ 360°360°/ 360°360°/ 360°270°/ 270°360°/360°360°/360°360°/360°
GradeOD/OS*0-1/VLE1/ 10-1/ 0-10-1/ 0-11 /11 /10-1/ 0-10-1/ 0-10-1/ 0-11 /13/33/33/3
Cornea1OD/OSEdema/EdemaClear/ClearClear/ ClearClear/ ClearClear/ ClearClear/ ClearClear/ ClearClear/ ClearClear/ ClearClear/ ClearClear/ ClearClear/ClearClear/ClearClear/Clear
Anterior chamber1OD/OSDeep/DeepShallow & clear/ Shallow & clearShallow & clear/ Shallow & clearShallow & clear/ Shallow & clearShallow & clear/ Shallow & clearShallow & clear/ Shallow & clearShallow & clear/ Shallow & clearShallow & clear/ Shallow & clearShallow & clear/ Shallow & clearShallow & clear/ Shallow & clearShallow & clear/ Shallow clearDeep/DeepDeep/DeepDeep/Deep
Lens1OD/OSAphakic/ AphakicSevere cataract/
VLE
Severe cataract/
moderate cataract
Severe cataract (OU)Clear/ ClearClear/ ClearClear/ ClearClear/ ClearClear/ ClearClear/ ClearClear/ ClearClear/ClearClear/ClearClear/Clear
Vertical C/D ratio1OD/OS0.80/ 0.800.40/VLE0.25/ 0.250.45/ 0.450.55/ 0.600.40/ 0.400.50/ 0.500.50/ 0.500.30/ 0.300.40/ 0.450.60/ 0.600.25/0.250.50/0.600.30/0.30
IOP2(mmHg)OD/OS144/14425/VLE26/2812/1218/1818/1916/1617/1729/3023/2314/1712/1211/13
RefractionOD/OS+13.75/ no reflex-1.25 -1.25x65/ VLE+2.25 -0.75x135/ +3.50 -0.25x60Plano -2.75/ +1.50+1.50 -0.75x90/ +1.75-0.50x85+0.75 -1.00x80/ +0.25 -1.00x180Plano -.25/ +0.25 -0.50x50Plano -0.50x3/ plano -0.25x180+1.25 -0.50x85/ +1.25 -0.25x85-0.75 -0.50x62/ -0.75 -1.00x95Plano/ plano -0.75x170Plano/PlanoPlano/Plano-3.50/ -3.50

(Continued)

Table 1

Phenotypic and genotypic data on pedigree ANG200 individuals with mutation in COL18A1

Individual IDIII-3III-4III-5III-6IV-2IV-3IV-5IV-6IV-7IV-8IV-9V-1V-2V-3&
COL18A1 genotypeNM?M#NMNMNMNMNMNMMMNMNMNMNMNN
Clinical examination
SexFMMFMFFMFMFFMF
Age at examination (yrs)82777569665649534841373540
Angle1
circumference OD/OS*360°/VLE270°/ 270°360°/ 360°360°/ 360°270°/ 270°270°/ 270°360°/ 360°360°/ 360°360°/ 360°270°/ 270°360°/360°360°/360°360°/360°
GradeOD/OS*0-1/VLE1/ 10-1/ 0-10-1/ 0-11 /11 /10-1/ 0-10-1/ 0-10-1/ 0-11 /13/33/33/3
Cornea1OD/OSEdema/EdemaClear/ClearClear/ ClearClear/ ClearClear/ ClearClear/ ClearClear/ ClearClear/ ClearClear/ ClearClear/ ClearClear/ ClearClear/ClearClear/ClearClear/Clear
Anterior chamber1OD/OSDeep/DeepShallow & clear/ Shallow & clearShallow & clear/ Shallow & clearShallow & clear/ Shallow & clearShallow & clear/ Shallow & clearShallow & clear/ Shallow & clearShallow & clear/ Shallow & clearShallow & clear/ Shallow & clearShallow & clear/ Shallow & clearShallow & clear/ Shallow & clearShallow & clear/ Shallow clearDeep/DeepDeep/DeepDeep/Deep
Lens1OD/OSAphakic/ AphakicSevere cataract/
VLE
Severe cataract/
moderate cataract
Severe cataract (OU)Clear/ ClearClear/ ClearClear/ ClearClear/ ClearClear/ ClearClear/ ClearClear/ ClearClear/ClearClear/ClearClear/Clear
Vertical C/D ratio1OD/OS0.80/ 0.800.40/VLE0.25/ 0.250.45/ 0.450.55/ 0.600.40/ 0.400.50/ 0.500.50/ 0.500.30/ 0.300.40/ 0.450.60/ 0.600.25/0.250.50/0.600.30/0.30
IOP2(mmHg)OD/OS144/14425/VLE26/2812/1218/1818/1916/1617/1729/3023/2314/1712/1211/13
RefractionOD/OS+13.75/ no reflex-1.25 -1.25x65/ VLE+2.25 -0.75x135/ +3.50 -0.25x60Plano -2.75/ +1.50+1.50 -0.75x90/ +1.75-0.50x85+0.75 -1.00x80/ +0.25 -1.00x180Plano -.25/ +0.25 -0.50x50Plano -0.50x3/ plano -0.25x180+1.25 -0.50x85/ +1.25 -0.25x85-0.75 -0.50x62/ -0.75 -1.00x95Plano/ plano -0.75x170Plano/PlanoPlano/Plano-3.50/ -3.50
Individual IDIII-3III-4III-5III-6IV-2IV-3IV-5IV-6IV-7IV-8IV-9V-1V-2V-3&
COL18A1 genotypeNM?M#NMNMNMNMNMNMMMNMNMNMNMNN
Clinical examination
SexFMMFMFFMFMFFMF
Age at examination (yrs)82777569665649534841373540
Angle1
circumference OD/OS*360°/VLE270°/ 270°360°/ 360°360°/ 360°270°/ 270°270°/ 270°360°/ 360°360°/ 360°360°/ 360°270°/ 270°360°/360°360°/360°360°/360°
GradeOD/OS*0-1/VLE1/ 10-1/ 0-10-1/ 0-11 /11 /10-1/ 0-10-1/ 0-10-1/ 0-11 /13/33/33/3
Cornea1OD/OSEdema/EdemaClear/ClearClear/ ClearClear/ ClearClear/ ClearClear/ ClearClear/ ClearClear/ ClearClear/ ClearClear/ ClearClear/ ClearClear/ClearClear/ClearClear/Clear
Anterior chamber1OD/OSDeep/DeepShallow & clear/ Shallow & clearShallow & clear/ Shallow & clearShallow & clear/ Shallow & clearShallow & clear/ Shallow & clearShallow & clear/ Shallow & clearShallow & clear/ Shallow & clearShallow & clear/ Shallow & clearShallow & clear/ Shallow & clearShallow & clear/ Shallow & clearShallow & clear/ Shallow clearDeep/DeepDeep/DeepDeep/Deep
Lens1OD/OSAphakic/ AphakicSevere cataract/
VLE
Severe cataract/
moderate cataract
Severe cataract (OU)Clear/ ClearClear/ ClearClear/ ClearClear/ ClearClear/ ClearClear/ ClearClear/ ClearClear/ClearClear/ClearClear/Clear
Vertical C/D ratio1OD/OS0.80/ 0.800.40/VLE0.25/ 0.250.45/ 0.450.55/ 0.600.40/ 0.400.50/ 0.500.50/ 0.500.30/ 0.300.40/ 0.450.60/ 0.600.25/0.250.50/0.600.30/0.30
IOP2(mmHg)OD/OS144/14425/VLE26/2812/1218/1818/1916/1617/1729/3023/2314/1712/1211/13
RefractionOD/OS+13.75/ no reflex-1.25 -1.25x65/ VLE+2.25 -0.75x135/ +3.50 -0.25x60Plano -2.75/ +1.50+1.50 -0.75x90/ +1.75-0.50x85+0.75 -1.00x80/ +0.25 -1.00x180Plano -.25/ +0.25 -0.50x50Plano -0.50x3/ plano -0.25x180+1.25 -0.50x85/ +1.25 -0.25x85-0.75 -0.50x62/ -0.75 -1.00x95Plano/ plano -0.75x170Plano/PlanoPlano/Plano-3.50/ -3.50

(Continued)

Table 1

(Continued)

Individual IDIII-3III-4III-5III-6IV-2IV-3IV-5IV-6IV-7IV-8IV-9V-1V-2V-3&
COL18A1 genotypeNM?M#NMNMNMNMNMNMMMNMNMNMNMNN
Biometrics
CCT3(μm)OD/OS554/ 564578/ 578530/ 527520/ 516534/ 531561/ 562559/ 557557/ 553550/ 555
Axial length3(mm)OD/OS22.69/ 22.9121.48/ 21.6023.24/ 23.0422.87/ 22.9122.82/ 22.7721.85/ 21.9023.12/ 23.2123.32/ 23.3522.79/ 22.6324.00/ 24.1024.11/ 24.07
ACD3(mm) OD/OS2.85/ 2.772.43/ 2.402.19/ 2.022.30/ 2.352.53/ 2.602.55/ 2.532.70/ 2.682.40/ 2.452.41/ 2.422.46/ 2.492.61/ 2.63
Lens thickness3(mm) OD/OS4.60/ 5.014.44/ 4.674.46/ 4.564.48/ 4.414.47/ 4.473.99/ 4.024.30/ 4.244.11/ 4.113.91/ 3.914.47/ 4.43
OCT
PP-OCTcOD/OSnormal/ normalminimal RNFL defect (OU)normal/ normal
Globe thicknessc(μm)OD/OS98/ 99
Visual fieldearly superior arcuate defect (OU)suspicious early defects (OU)normal/ normalnormal/ normalnormal/ Normal
TreatmentCataract surgery (OU)Cataract surgery (OD)Cataract surgery (OU)Cataract surgery (OU)Laser PILaser PILaser PILaser PI + medicationLaser PILaser PI
CommentsDiabetic retinopathy (OU)PAS observed (OU)Epiretinal membrane (OU)
DiagnosisGlaucoma (OU)fPACGPAC (OU)PACG (OU)PACS (OU)PACS (OU)PACS (OU)PACS (OU)PACS (OU)PAC (OU)PAC (OU)NEENEENEE/Low myopia (OU)
Individual IDIII-3III-4III-5III-6IV-2IV-3IV-5IV-6IV-7IV-8IV-9V-1V-2V-3&
COL18A1 genotypeNM?M#NMNMNMNMNMNMMMNMNMNMNMNN
Biometrics
CCT3(μm)OD/OS554/ 564578/ 578530/ 527520/ 516534/ 531561/ 562559/ 557557/ 553550/ 555
Axial length3(mm)OD/OS22.69/ 22.9121.48/ 21.6023.24/ 23.0422.87/ 22.9122.82/ 22.7721.85/ 21.9023.12/ 23.2123.32/ 23.3522.79/ 22.6324.00/ 24.1024.11/ 24.07
ACD3(mm) OD/OS2.85/ 2.772.43/ 2.402.19/ 2.022.30/ 2.352.53/ 2.602.55/ 2.532.70/ 2.682.40/ 2.452.41/ 2.422.46/ 2.492.61/ 2.63
Lens thickness3(mm) OD/OS4.60/ 5.014.44/ 4.674.46/ 4.564.48/ 4.414.47/ 4.473.99/ 4.024.30/ 4.244.11/ 4.113.91/ 3.914.47/ 4.43
OCT
PP-OCTcOD/OSnormal/ normalminimal RNFL defect (OU)normal/ normal
Globe thicknessc(μm)OD/OS98/ 99
Visual fieldearly superior arcuate defect (OU)suspicious early defects (OU)normal/ normalnormal/ normalnormal/ Normal
TreatmentCataract surgery (OU)Cataract surgery (OD)Cataract surgery (OU)Cataract surgery (OU)Laser PILaser PILaser PILaser PI + medicationLaser PILaser PI
CommentsDiabetic retinopathy (OU)PAS observed (OU)Epiretinal membrane (OU)
DiagnosisGlaucoma (OU)fPACGPAC (OU)PACG (OU)PACS (OU)PACS (OU)PACS (OU)PACS (OU)PACS (OU)PAC (OU)PAC (OU)NEENEENEE/Low myopia (OU)

M, mutated COL18A1 allele; N, normal COL18A1 allele; OD, oculus dexter (right eye), OS, oculus sinister (left eye); OU, oculus uterque (both eyes); C/D ratio, cup to disc ratio; IOP, intraocular pressure; CCT, central corneal thickness; ACD, anterior chamber depth; OCT, optical coherence tomography; PAS; peripheral anterior synechiae; PP-OCT, peripapillary optical thickness of retinal nerve fiber layer; VLE, virtually lost eye (due to severe glaucoma); laser PI, laser peripheral iridotomy; NEE, normal eye exam. 1, by slit lamp biomicroscopy; 2, by Goldmann applanation tonometry; 3, by laser interferometry (LenStar); 4, with medication; #, genotype inferred; *, glaucomatous optic nerve damage is evident, but angle closure not assessable due to prior intracapsular cataract extraction eyes resulting from surgery; &, example of data on individual with normal genotype. Blanks indicate missing data.

Table 1

(Continued)

Individual IDIII-3III-4III-5III-6IV-2IV-3IV-5IV-6IV-7IV-8IV-9V-1V-2V-3&
COL18A1 genotypeNM?M#NMNMNMNMNMNMMMNMNMNMNMNN
Biometrics
CCT3(μm)OD/OS554/ 564578/ 578530/ 527520/ 516534/ 531561/ 562559/ 557557/ 553550/ 555
Axial length3(mm)OD/OS22.69/ 22.9121.48/ 21.6023.24/ 23.0422.87/ 22.9122.82/ 22.7721.85/ 21.9023.12/ 23.2123.32/ 23.3522.79/ 22.6324.00/ 24.1024.11/ 24.07
ACD3(mm) OD/OS2.85/ 2.772.43/ 2.402.19/ 2.022.30/ 2.352.53/ 2.602.55/ 2.532.70/ 2.682.40/ 2.452.41/ 2.422.46/ 2.492.61/ 2.63
Lens thickness3(mm) OD/OS4.60/ 5.014.44/ 4.674.46/ 4.564.48/ 4.414.47/ 4.473.99/ 4.024.30/ 4.244.11/ 4.113.91/ 3.914.47/ 4.43
OCT
PP-OCTcOD/OSnormal/ normalminimal RNFL defect (OU)normal/ normal
Globe thicknessc(μm)OD/OS98/ 99
Visual fieldearly superior arcuate defect (OU)suspicious early defects (OU)normal/ normalnormal/ normalnormal/ Normal
TreatmentCataract surgery (OU)Cataract surgery (OD)Cataract surgery (OU)Cataract surgery (OU)Laser PILaser PILaser PILaser PI + medicationLaser PILaser PI
CommentsDiabetic retinopathy (OU)PAS observed (OU)Epiretinal membrane (OU)
DiagnosisGlaucoma (OU)fPACGPAC (OU)PACG (OU)PACS (OU)PACS (OU)PACS (OU)PACS (OU)PACS (OU)PAC (OU)PAC (OU)NEENEENEE/Low myopia (OU)
Individual IDIII-3III-4III-5III-6IV-2IV-3IV-5IV-6IV-7IV-8IV-9V-1V-2V-3&
COL18A1 genotypeNM?M#NMNMNMNMNMNMMMNMNMNMNMNN
Biometrics
CCT3(μm)OD/OS554/ 564578/ 578530/ 527520/ 516534/ 531561/ 562559/ 557557/ 553550/ 555
Axial length3(mm)OD/OS22.69/ 22.9121.48/ 21.6023.24/ 23.0422.87/ 22.9122.82/ 22.7721.85/ 21.9023.12/ 23.2123.32/ 23.3522.79/ 22.6324.00/ 24.1024.11/ 24.07
ACD3(mm) OD/OS2.85/ 2.772.43/ 2.402.19/ 2.022.30/ 2.352.53/ 2.602.55/ 2.532.70/ 2.682.40/ 2.452.41/ 2.422.46/ 2.492.61/ 2.63
Lens thickness3(mm) OD/OS4.60/ 5.014.44/ 4.674.46/ 4.564.48/ 4.414.47/ 4.473.99/ 4.024.30/ 4.244.11/ 4.113.91/ 3.914.47/ 4.43
OCT
PP-OCTcOD/OSnormal/ normalminimal RNFL defect (OU)normal/ normal
Globe thicknessc(μm)OD/OS98/ 99
Visual fieldearly superior arcuate defect (OU)suspicious early defects (OU)normal/ normalnormal/ normalnormal/ Normal
TreatmentCataract surgery (OU)Cataract surgery (OD)Cataract surgery (OU)Cataract surgery (OU)Laser PILaser PILaser PILaser PI + medicationLaser PILaser PI
CommentsDiabetic retinopathy (OU)PAS observed (OU)Epiretinal membrane (OU)
DiagnosisGlaucoma (OU)fPACGPAC (OU)PACG (OU)PACS (OU)PACS (OU)PACS (OU)PACS (OU)PACS (OU)PAC (OU)PAC (OU)NEENEENEE/Low myopia (OU)

M, mutated COL18A1 allele; N, normal COL18A1 allele; OD, oculus dexter (right eye), OS, oculus sinister (left eye); OU, oculus uterque (both eyes); C/D ratio, cup to disc ratio; IOP, intraocular pressure; CCT, central corneal thickness; ACD, anterior chamber depth; OCT, optical coherence tomography; PAS; peripheral anterior synechiae; PP-OCT, peripapillary optical thickness of retinal nerve fiber layer; VLE, virtually lost eye (due to severe glaucoma); laser PI, laser peripheral iridotomy; NEE, normal eye exam. 1, by slit lamp biomicroscopy; 2, by Goldmann applanation tonometry; 3, by laser interferometry (LenStar); 4, with medication; #, genotype inferred; *, glaucomatous optic nerve damage is evident, but angle closure not assessable due to prior intracapsular cataract extraction eyes resulting from surgery; &, example of data on individual with normal genotype. Blanks indicate missing data.

Among 10 pedigree members in whom some degree of angle closure was definitively observed, 5 were diagnosed with PACS, 3 with PAC and 2 with PACG. An additional member, ANG200 III-3, was diagnosed with glaucoma based on presence of glaucomatous optic nerve damage, but angle closure was not assessable due to prior intracapsular cataract extraction surgery. Images pertaining to the proband (ANG200 IV-9) consistent with PAC diagnosis in this individual are presented in Figures 2 and 3. The slit-lamp photographs show narrow angle configuration and a shallow anterior chamber in both of her eyes. The TM is not visible due to an anteriorly positioned iris root (Fig. 2). Fundus images of her eyes show normal optic nerve heads (Fig. 3A), and optical coherence tomography demonstrates normal thickness of the peripapillary retinal nerve fiber layer in both eyes (Fig. 3B). Visual field data evidence a full and normal field, without definite signs of glaucomatous damage in either eye (Fig. 3C). Slit-lamp photographs of another patient (ANG200 IV-8) that show narrow angle configuration are presented in Supplementary Figure S1. Unfortunately, images from the PACG-diagnosed individuals are not available.

Slit-lamp photographs of proband IV-9 affected with PAC. (A) Presence of a shallow anterior chamber in both eyes is evident. (B) The angle configuration is narrow as observed using Van Herick's method. (C) The gonioscopic appearance is consistent with angle closure disease; the TM is not visible due to iris root apposition. Also, the bulging lens–iris diaphragm is to be noted. Images to the left and right, respectively, pertain to the right and left eye.
Figure 2

Slit-lamp photographs of proband IV-9 affected with PAC. (A) Presence of a shallow anterior chamber in both eyes is evident. (B) The angle configuration is narrow as observed using Van Herick's method. (C) The gonioscopic appearance is consistent with angle closure disease; the TM is not visible due to iris root apposition. Also, the bulging lens–iris diaphragm is to be noted. Images to the left and right, respectively, pertain to the right and left eye.

Data that evidence absence of end-organ glaucomatous damage in eyes of proband IV-9. (A) Fundus photographs of the right (image at left) and left (image at right) eyes reveal normal and healthy-appearing optic nerve heads in both eyes. (B) Peripapillary optical coherence tomography images depicting the retinal nerve fiber layer of the right (OD) and left eyes (OS) demonstrates normal and stable thickness values from 2014 to 2015 in both eyes. (C) Visual field examination of the right and left eyes (left and right images, respectively) shows a full and normal field with no sign of glaucomatous damage.
Figure 3

Data that evidence absence of end-organ glaucomatous damage in eyes of proband IV-9. (A) Fundus photographs of the right (image at left) and left (image at right) eyes reveal normal and healthy-appearing optic nerve heads in both eyes. (B) Peripapillary optical coherence tomography images depicting the retinal nerve fiber layer of the right (OD) and left eyes (OS) demonstrates normal and stable thickness values from 2014 to 2015 in both eyes. (C) Visual field examination of the right and left eyes (left and right images, respectively) shows a full and normal field with no sign of glaucomatous damage.

Attention to the paraclinical measurements presented in Table 1, including those pertaining to angle closure, IOP, cup/disc ratio, refraction, central corneal thickness, anterior chamber depth, lens thickness and axial length, evidences that only measurements pertaining to angle closure itself (angle circumference/angle grade) are outside the normal range in all 10 individuals diagnosed with PACS, PAC or PACG. IOP is within the normal range among the PACS individuals and elevated in the PAC and PACG patients. There are deviations from a normal cup/disc ratio in all three groups of patients, but extent of deviations do not strictly correlate with disease severity (Table 1). The three individuals of generation V had normal eye examinations.

Genetic analysis of ANG200 pedigree

Whole genome linkage analysis with single nucleotide polymorphism (SNP) markers was performed to identify one or more candidate loci that segregated with phenotypic presentations in ANG200. A monogenic inheritance pattern of angle closure in ANG200 was not indisputably evident. Given that the proband and her four siblings with some degree of angle closure and diagnosed with PACS or PAC were offsprings of a consanguineous marriage, the possibility of autosomal recessive inheritance was considered at least for this branch of the pedigree. However, autozygosity mapping did not reveal a homozygous region >1 Mb common to the five siblings and a positive logarithm of odds (LOD) score was not obtained anywhere in the genome under a recessive model. Also, a positive LOD score under a recessive model was not attained when all individuals genotyped were considered, further reducing the likelihood of recessive inheritance. Quality control reports on the SNP genotypings were good (call rates, >0.99), suggesting homozygous regions were not undetected because of technical issues. Contrary to analysis under a recessive model, a 2.1 Mb locus on chromosome 21q22.3 with a LOD score of 0.941 was identified upon analysis of the SNP genotype data of the branch of the pedigree that included the proband, her siblings and her parents under an autosomal dominant model. When data on all individuals genotyped was included in the analysis under a dominant model, the LOD score increased to 1.672 (Fig. 4A). The locus contained 255 marker SNPs and was bound by proximal and distal markers, rs403603 (21:44554443) and rs7278087 (21:46678912), respectively. The haplotype defined by the nucleotides at the 255 SNP positions and shared by all the affected individuals included in the analysis is presented in Supplementary Table S1-A. Marker rs7278087 was the right-most marker on chromosome 21. The region from the position 21:44554443 to the end of chromosome contained 70 annotated protein-encoding genes and 166 open-reading frames (http://www.ncbi.nlm.nih.gov/genome/gdv/).

Disease-associated locus and disease segregating sequence variations in pedigree ANG200. (A) LOD plot of chromosome 21 under autosomal dominant model. (B) Chromatograms of sequence variation in COL18A1 (left) and DIP2A (right). Top: homozygous for wild-type allele; middle: heterozygous; bottom: homozygous for variant allele.
Figure 4

Disease-associated locus and disease segregating sequence variations in pedigree ANG200. (A) LOD plot of chromosome 21 under autosomal dominant model. (B) Chromatograms of sequence variation in COL18A1 (left) and DIP2A (right). Top: homozygous for wild-type allele; middle: heterozygous; bottom: homozygous for variant allele.

Given the large number of genes and the nominal LOD score of the locus, all the exons in the genomes of the three siblings ANG200 IV-5, IV-8 and IV-9 were sequenced. After filtering the exome data, no homozygous region ≥1 Mb common to the three siblings was found, consistent with results of linkage analysis that had ruled out recessive inheritance. Many heterozygous sequence variations that cause amino acid or splicing changes distributed throughout the genome with minimal-allele frequencies of ≤0.03 common to the three exomed individuals were identified. A total of 44 of these heterozygous variations were novel at time of analysis, and all of these were genotyped by Sanger sequencing in the 8 individuals who had been included in the linkage analysis and also in ANG200 IV-2, ANG200 IV-3 who had newly been recruited (Tables S2 and S3). Two variations, c.550G>A in COL18A1 and c.2489C>T in DIP2A, that encode, collagen type XVIII, alpha-1 protein and Disco interacting protein 2 homolog A, were present in all the individuals with angle closure and diagnosed with PACG, PAC or PACS, respectively (Fig. 4B). The variations were also present in individual III-3 who had been diagnosed with glaucoma, but for whom definitive information on the closure status of the iridocorneal angle in her eyes prior to surgery was not available (Tables 1 and S3). The nucleotide variations in COL18A1 and DIP2A cause p.Glu184Lys and p.Ala830Val in the coded proteins, respectively. Both genes are positioned on chromosome 21, within the locus that had been identified by linkage analysis. The variations in the two genes are separated by 1090929 bp, suggesting close linkage. Both variations were screened in ethnically matched healthy elderly control individuals; the c.550G>A change in COL18A1 was not observed in 400 controls whereas the c.2489C>T change in DIP2A was observed in 9 of 100 controls. DIP2A was an interesting candidate gene because a reporter gene study had suggested high expression of Dip2a in many mouse tissues, especially in nervous system tissues including the retina (34). But the DIP2A variation in ANG200 appears to be a relatively common variation in the Iranian population that is not associated with disease status.

Glycine at position corresponding to amino acid 184 in human collagen type XVIII that was changed by the mutation in pedigree ANG200 is conserved among mammals (Table S4). This amino acid is positioned within a non-collagenous domain specific to collagen type XVIII and absent in paralogous human proteins (Fig. 5) (see Discussion). The in silico bioinformatics tools SIFT (http://sift.jcvi.org/) and PolyPhen-2 (http://genetics.bwh.harvard.edu/pph2/), predicted the amino acid change p.Glu184Lys as damaging and probably damaging, respectively. Combined Annotation Dependent Depletion (CADD) http://cadd.gs.washington.edu/info), which is a tool for scoring the deleteriousness of single nucleotide variants as well as insertion/deletion variants in the human genome based on over 63 distinct annotations, calculated a C-score of 18 for the mutation and marked it to be among the 5.3% most deleterious variations (35).

Schematic representation of protein structure of long isoform of collagen type XVIII. Positions of the disease causing mutations including p.Glu184Lys in pedigree ANG200, p.Arg505ValfsX34 in KS100 and p.Leu612TrpfsX23 in KS101 are shown. Col, collagenous domain; NC, non-collagenous domain.
Figure 5

Schematic representation of protein structure of long isoform of collagen type XVIII. Positions of the disease causing mutations including p.Glu184Lys in pedigree ANG200, p.Arg505ValfsX34 in KS100 and p.Leu612TrpfsX23 in KS101 are shown. Col, collagenous domain; NC, non-collagenous domain.

Although the findings strongly implicated the variation in COL18A1 as the causative mutation of angle closure in pedigree ANG200, there remained due motivation for concern. Specifically, it could be argued that screening of only novel sequence variations as candidate causative mutations for a non-rare presentation like angle closure may be inappropriate. The fact that extent of angle closure in affected ANG200 individuals was variable also begged care in analysis and interpretation of the genetic data (Table 1). With these considerations and given the importance of angle closure in the etiology of PACG, we sequenced the exomes of all available ANG200 members who had not already been sequenced to unequivocally and definitively identify the causative mutation. All the exome data was filtered to include sequence variations in the genome with a minimum allele frequency of ≤0.03 and present in all nine living individuals with angle closure and diagnosed with PACG, PAC or PACS. Although the prevalence of angle closure in Iran is unknown, the rather conservative minor allele frequency (MAF) cutoff of ≤0.03 was used to avoid losing data pertaining to a relatively frequent potential causative allele. The analysis confirmed the earlier findings and showed that c.550G>A in COL18A1 and c.2489C>T in DIP2A were the only sequence variations that satisfied the selection criteria. The variation in DIP2A was ruled out as a candidate causative mutation because of its high frequency in control individuals. The haplotype defined by the nucleotides at SNP positions within the exome sequencing data and shared by all the affected individuals included in the analysis is presented in Supplementary Table S1-B. PACG-diagnosed individual ANG200 III-4 (who had not been exomed) was inferred to carry a COL18A1-mutated allele because one of his children was homozygous for the mutation. It was concluded that COL18A1 indeed likely affects angle closure. The putative disease-associated genotype was present in two generation V individuals (ANG200 V-1 and V-2) with a normal eye examination who had not been exomed. We suggest that absence of angle closure in these individuals is because they are relatively young (35 and 37 years old) and we expect that they will later manifest the condition (see Discussion). ANG200 V-3, who was a few years older and diagnosed normal, had a wild-type COL18A1 genotype.

PACS diagnosis in carriers of COL18A1 mutations in KS pedigrees KS100 and KS101

To gain further support for the proposal that COL18A1 affects angle closure, attempts were made to identify mutations in the gene in other unrelated angle closure-affected individuals. As COL18A1 has more than 40 exons and as mutations in the gene were not expected to be common causes of angle closure-related phenotypes, mutation screening in randomly chosen affected individuals was not considered. However, it was known that COL18A1 had earlier been identified as the causative gene of KS (OMIM 267750), which is a very rare autosomal recessive disorder first described in 1971 (36). Although KS is associated with considerable clinical variability, severe optical abnormalities including high myopia that often leads to bilateral blindness and occipital encephalocele are its most common features (37,38). Assuming high penetrance of mutated COL18A1 alleles, it was expected that parents and individuals of earlier generations of KS patients may manifest with angle closure. A look out for KS-affected individuals ultimately resulted in identification of one patient in each of two consanguineous families, KS100 1ne KS101 (Fig. 6A). The patients had been referred to the retina service of Labafi Nejad Hospital because of poor vision. Cyclo-refraction revealed severe high myopia and retinal examination showed macular hypoplasia (posterior staphyloma) and generalized chorio-retinal atrophy in the eyes of both children. A preliminary KS diagnosis was made for the children. Sequencing of COL18A1 resulted in identification of novel homozygous deleterious mutations c.1513delC (p.Arg505ValfsX34) and c.1834delC (p.Leu612TrpfsX23) in the patients, thus confirming the clinical-based diagnosis of KS. Results of clinical examinations and genetic analysis of members of the pedigrees are presented in Supplementary Table S5 as well as Figures 6B and 7. In KS100, the father, maternal grandmother and paternal grandmother of the KS-affected probands were diagnosed with PACS, and these individuals were carriers of the c.1513delC mutation in COL18A1. In KS101, the parents and paternal grandfather of the proband were recruited. The parents had normal eye examination; only the grandparent who was 65 years old and 22 and 30 years older than his son and daughter-in-law, was diagnosed with PACS, respectively. All three individuals were heterozygous carriers of the c.1834delC mutation. Diagnosis of PACS in individuals who are carriers of different mutations in COL18A1 strongly supports the conjecture that mutations in this gene can promote angle closure. Furthermore, absence of angle closure indications in most relatively young carriers of the mutations (individuals in mid-40s or younger) suggests that the phenotypic effects of mutations in COL18A1 with respect to angle closure are relatively late onset.

Knobloch syndrome families KS100 and KS101. (A) Pedigrees of the affected probands. Black-filled shapes represent the KS affected probands. *, individuals in family who underwent ophthalmic examination and were screened for mutation in respective proband. Clinical diagnoses, age at time of examination (years) and genotypes are indicated. NEE, normal eye examination; N, normal COL18A1 allele; M, mutated COL18A1 allele. (B) Chromatograms of sequence variation in COL18A1 in KS100 (left) and KS101 (right). Top: homozygous for wild-type allele; middle: heterozygous; bottom: homozygous for variant allele.
Figure 6

Knobloch syndrome families KS100 and KS101. (A) Pedigrees of the affected probands. Black-filled shapes represent the KS affected probands. *, individuals in family who underwent ophthalmic examination and were screened for mutation in respective proband. Clinical diagnoses, age at time of examination (years) and genotypes are indicated. NEE, normal eye examination; N, normal COL18A1 allele; M, mutated COL18A1 allele. (B) Chromatograms of sequence variation in COL18A1 in KS100 (left) and KS101 (right). Top: homozygous for wild-type allele; middle: heterozygous; bottom: homozygous for variant allele.

Slit-lamp photographs of PACS affected individual KS100 IV2, heterozygous carrier of COL18A1 mutation p.Arg505ValfsX34. (A) Images show shallow anterior chamber in both eyes, and narrow angle configuration is evident. (B) Gonioscopic appearance is consistent with angle closure.
Figure 7

Slit-lamp photographs of PACS affected individual KS100 IV2, heterozygous carrier of COL18A1 mutation p.Arg505ValfsX34. (A) Images show shallow anterior chamber in both eyes, and narrow angle configuration is evident. (B) Gonioscopic appearance is consistent with angle closure.

Discussion

We have identified mutation c.550G>A in COL18A1 that causes p.Glu184Lys in the encoded protein, collagen type XVIII alpha-1 protein, as likely cause of iridocorneal angle closure in pedigree ANG200 with multiple members diagnosed with PACS, PAC or PACG. Inheritance pattern of angle closure caused by this mutation was autosomal dominant. Given the depth of the genetic investigation that includes whole genome linkage analysis and exome sequencing of 10 individuals, and that all novel variations and earlier reported variations with a MAF of ≤0.03 were considered, it is highly likely that the identified mutation could play a role in PACG manifestation in the pedigree. Identification of different COL18A1 mutations in PACS-affected individuals of two additional unrelated families provides additional strong support for an association between mutations in this gene and angle closure. These individuals were parents or grandparents of KS patients. Abnormal clinical presentations in parents of KS patients have not previously been reported, even in cases where genetic analysis confirmed their carrier status (37,38). Potential explanations for this include the possibility that non-severe presentations of PACS and PAC were overlooked in the framework of the more severe disease of offsprings or that they were not evident because the parents were young. Our findings are consistent with the latter and suggest that COL18A1 mutations may usually affect angle closure in the fourth decade of life or later. Four heterozygous carriers of COL18A1 mutations (ANG200 V-1, ANG200 V-2, KS101 III-1 and KS101 III-2) without indications of angle closure were identified, and their age at time of examination ranged from 35 to 43 years.

Additional to presence of different COL18A1 mutations in individuals with angle closure from three unrelated families, considerable supportive evidence exists for the conjecture of roles for COL18A1 in eye structure and/or functions. The first is that ocular anomaly is a defining feature of KS, a disease caused by mutations in COL18A1. Interestingly, glaucoma along with a plethora of other presentations has twice been reported in KS patients (3941). Supportive evidence also derives from several studies in mice. Retinal pigment epithelium degeneration and severe defects in the iris were reported in Col18a1 null mice; the defects were attributed to disturbed proteostasis or to altered properties of the basal membrane (42). One of the mouse studies specifically concluded that lack of type XVIII collagen results in anterior ocular defects (43). Unfortunately, phenotype data on heterozygous carriers of null COL18A1 mutations was not available. Also, COL18A1 has been shown to be expressed in human eye tissues including the cornea and the ciliary body (44,45). The ciliary body is positioned in the anterior chamber of the eye, and is the source of aqueous humor. Equally relevant to the present study, COL18A1 expression in the human TM has been reported in several studies (4649). Finally, roles for collagens and for the extracellular matrix in the pathogenesis of glaucoma have previously been proposed and extensively discussed (see below).

Although our finding of a putative role of COL18A1 in affecting iridocorneal angle closure is interesting in itself, it is of particular importance because of the association between angle closure and PACG. PACG is an important public health entity. It is estimated that 15.7 million individuals in the world are affected with PACG. It is projected that 21 million will be affected by 2020, and that PACG by that time will cause bilateral blindness in 5.3 million people (5052). Most PACG patients are from Asia, particularly China, Mongolia, Singapore and India (5356). As implicated in the name of the disease, angle closure is a prominent manifestation of PACG, and the presence of an occludable angle is in fact frequently used in epidemiological studies to assess PACG (15). It is thus tempting to suggest that COL18A1 may in some cases contribute to the etiology of PACG. Consistent with this, two individuals with the COL18A1 mutation in ANG200 were diagnosed with PACG. Seven others were diagnosed with PACS or PAC, and PACS and PAC affected individuals are predisposed to later manifest PACG. The potential and the extent of COL18A1 contribution to PACG disease burden awaits large scale screenings in PACG affected individuals. The fact that the in silico tool CADD assessed the c.550G>A mutation in COL18A1 to be only among the 5.3% most deleterious variations among known human genome variations, suggests that it is not highly deleterious. In fact, some bioinformatics tools including MutationTaster (http://www.mutationtaster.org/) and Mutation Assessor (http://mutationassessor.org/r3/) did not designate the COL18A1 mutation in ANG200 as damaging. Consistent with the mutation not being highly deleterious, the clinical presentations of the only ANG200 member (ANG200 IV-7) with a homozygous mutant genotype were not more severe than in the heterozygous individuals, and a KS diagnosis in this patient is certainly not called for. Finally, it is to be noted that the four eldest members of ANG200 diagnosed with glaucoma were also affected with cataract, and it cannot be excluded that the mutation in COL18A1 contributed to the latter condition.

The murine and human COL18A1 genes were first identified in 1994 (57,58). Type XVIII collagen that is encoded by the gene is a non-fibrillar proteoglycan collagen that forms homotrimers (59). Type XVIII collagen is expressed in many tissues including ocular tissues (44,6063). This protein has received much attention because endostatin, which inhibits angiogenesis and affects various endothelial functions, is derived from its C-terminus by proteolytic cleavage (6467). Three isoforms of type XVIII collagen that contain 1339, 1519 and 1754 amino acids have been described (59). The protein structure of all isoforms is defined by the presence of an N-terminal non-collagenous domain (NC1), 10 collagenous repeats (Col1-10) interrupted by 9 non-collagenous repeats (NC2-10) and a C-terminal non-collagenous domain (NC11) (Fig. 5). The NC1 and NC11 domains are specific to type XVIII collagen. The isoforms of type XVIII collagen differ only in their N-terminal NC-domain, and are derived from COL18A1 by use of either of two promoters and by alternative splicing (59). NC1 in its entirety includes a signal peptide, a coiled coil motif and a cysteine-rich frizzled domain homologous to the extracellular ligand binding domain of frizzled proteins involved in Wnt signaling (Fig. 5) (59). The cysteine-rich frizzled domain is present only in the longest type XVIII collagen isoform, and fragments containing this domain have been identified in the medium of cultured cells (59). It has been suggested that this fragment may affect Wnt signaling pathways (59). The p.Glu184 that is mutated in pedigree ANG200 exists in the two longest XVIII collagen isoforms and its position with respect to the longest isoform is between the coiled coil motif and the frizzled-like domain (Fig. 5). The effects of the mutation on collagen XVIII functions, including possible effects on roles of the frizzled-like domain are unknown. Interestingly, membrane frizzled-related protein encoded by MFRP that has also been implicated in the etiology of PACG, also contains a frizzled-like domain (68). Additionally, it has been proposed that Wnt signaling is involved in the regulation of eyeball size, an endophenotype of PACG (6972). The effects of the KS families’ COL18A1 mutations on protein structure are not focused as they cause early premature transcriptional termination (Fig. 5).

Finally, it is to be noted that in addition to type XVIII collagen, evidence for roles of various collagens in the pathogenesis of various forms of glaucoma is abundant in the literature (7376). This is perhaps best interpreted in light of the role of the extracellular matrix allotted to the etiology of glaucoma (8,73,77). Eight loci were identified in one or both of two recent large genome-wide association (GWA) studies on PACG patients, and COL11A1 that encodes one of the alpha chains of type XI collagen was one of the genes identified in both of the studies (19,21). An SNP in COL1A1 was associated with increased risk of myopia in Japanese and Chinese individuals (78,79). A GWA study of PACG in a dog breed identified COL1A2 as a susceptibility locus (80). In addition to collagens, genes with roles in the extracellular matrix and also associated with glaucoma include MMP-9 (8,22,77,8186). It has been suggested that inter-individual differences in tolerance to IOP as reflected in glaucoma diagnosis with normal tension in some and ocular hypertension without glaucoma in others reflect variations in biomechanical properties of the extracellular matrix of relevant ocular tissues (73). The finding of COL18A1 mutations as cause of angle closure in ANG200 and two other families reinforces the importance of collagens in eye features and functions.

Materials and Methods

This research was performed in accordance with the Declaration of Helsinki and with approval of the ethics board of the University of Tehran and the Ophthalmic Research Center at Shahid Beheshti University of Medical Sciences. Participants or guardians consented to participate in the study.

Clinical evaluations of members of ANG200 pedigree

Complete ophthalmologic examinations that included slit-lamp biomicroscopy, measurement of IOP, gonioscopy, fundus examination and standard achromatic perimetry were performed by glaucoma specialists (S.Y. & M.J.) on available individuals of the pedigree when possible. IOP measurements were obtained using a Goldmann applanation tonometer (AT900; Haag-Streit AG, Koeniz, Switzerland), and visual fields were assessed using the Humphrey Zeiss 750 Visual Field Analyzer (Carl Zeiss Meditec, Dublin, CA). Biometric measurements were obtained by laser interferometry (LenStar LS 900; Haag Streit, Switzerland). These measurements included central corneal thickness, axial length, anterior chamber depth and lens thickness. No mydriatic or cycloplegic agents were used. Indentation gonioscopy was used for grade classification of angle closure: grade 0, no angle structures visible; grade 1, only Schwalbe’s line visible; grade 2: TM also visible; grade 3: scleral spur also visible; grade 4, ciliary band also visible. Diagnosis was based on criteria recommended by the International Society of Geographical and Epidemiological Ophthalmology (http://iceh.lshtm.ac.uk/isgeo/). PACS was diagnosed when the posterior TM was not visible in at least 180 degrees of the angle circumference on gonioscopy in the absence of synechiae, normal IOP and a healthy optic nerve head. PAC was diagnosed when PACS criteria mentioned above were fulfilled together with structural or functional evidence of drainage angle damage (i.e. presence of peripheral anterior synechiae and raised IOP, respectively). PACG was diagnosed in the presence of PAC criteria and end-organ glaucomatous damage, that is, glaucomatous optic nerve head changes and/or visual field defects.

Genetic analysis of ANG200 pedigree

Genome-wide linkage analysis.

Genome-wide SNP genotyping of eight individuals of the ANG200 pedigree was carried out using HumanCytoSNP-12v2-1_C BeadChips (Illumina, San Diego, CA, USA) (Fig. 1, Table S2). ANG200 III-4, ANG200 IV-2, ANG200 IV-3 and individuals of generation V were not genotyped because they had died prior to start of linkage analysis or because they were identified and recruited after the analysis had been performed. Autozygosity mapping pertaining to data of individuals ANG200 IV-5–IV-9 was performed to identify common homozygous regions with a minimal physical length of 1 Mb using the Homozygous Detector Tool within the GenomeStudio_Genotyping _Module (Illumina). Merlin was used to attain LOD scores under assumption of both autosomal recessive and autosomal dominant inheritance and disease-allele frequency of 0.001 (87).

Exome sequencing.

Exome sequencing was initially performed using the DNA of three members of pedigree ANG200 (ANG200 IV-5, IV-8 and IV-9), and subsequently on the DNA of seven additional members (Fig. 1, Table S2). Exons were enriched using the TruSeq Exome Enrichment kit and sequenced on an Illumina HiSeq 2000 system (Illumina). Sequence alignment and variant calling were performed against human reference genome UCSC NCBI37/hg19 using CASAVA software (Illumina). Additionally, NextBio (http://www.nextbio.com/b/nextbio.nb), ENSEMBL (http://asia.ensembl.org/info/docs/tools/index.html) and Enlis Genomics (http://www.enlis.com/) software programs were used for variant detection. Subsequently, a first file of variations was prepared by removing SNPs with a MAF of >0.03 in the dbSNP database (http://www.ncbi.nlm.nih.gov/), the 1000 Genomes databases (www.1000genomes.org), the NHLBI Exome Sequencing Project (http://evs.gs.washington.edu/EVS/) or the Exome Aggregation Consortium database (http://exac.broadinstitute.org/) and SNPs observed in the exomes of 22 unrelated Iranians affected with non-ocular diseases. Variations in this file that did not affect amino acid change or splicing were also filtered out. The same protocol was used for copy number variations. A second file consisting of only novel variations was then prepared by removing from the first file all variations previously reported in databases described above.

Screening of 44 novel variations in the ANG200 pedigree.

All novel variations common to the first three individuals exomed (44 variations) were screened by Sanger sequencing in all available ANG200 pedigree members. This was done prior to exome sequencing of seven additional individuals.

Screening of c.550G>A in COL18A1 and c.2489C>T in DIP2A in control individuals.

COL18A1 variation c.550G>A and c.2489C>T variation in DIP2A were screened, respectively, in 400 and 100 unrelated control Iranian individuals older than 60 years of age and without self-reported familial history of ocular diseases by allele-specific polymerase chain reaction (PCR) protocols. Older individuals were used as controls because of consideration of possible late onset of manifestations. COL18A1 reference sequences used were NG_011903.1, NM_130444.2 and NP_569711.2, and DIP2A reference sequences used were NG_015996.1, NM_015151.3 and NP_055966.2.

Clinical evaluations and genetic analysis of KS pedigrees KS100 and KS101

KS diagnosis in two young unrelated children was based on presence of severe eye defects including high myopia and severe macular atrophy. Neither patient had other KS affected family members. KS diagnosis in the children was confirmed by sequencing of exons of COL18A1, which is the known causative gene of KS (BGI Clinical Laboratories, Hong Kong). Upon identification of KS causing mutations, available members of the families (KS100 and KS101) of the two patients underwent ophthalmic examinations as described above. Subsequently, the KS causing mutations were screened in the family members by Sanger sequencing.

Acknowledgements

We acknowledge the Iran National Science Foundation and the Ophthalmic Research Center of Shahid Beheshti University of Medical Sciences for funding this research, and the Ophthalmic Research Center of Shahid Beheshti University of Medical Sciences for the fellowship given to F.S.. We also express gratitude to ANG200, KS100 and KS101 pedigrees members for their cooperation in this study.

Conflict of Interest statement. None declared.

Funding

Iran National Science Foundation; Ophthalmic Research Center of Shahid Beheshti University of Medical Sciences (93227 and 90203).

References

1

Ray
,
K.
,
Mukhopadhyay
,
A.
and
Acharya
,
M.
(
2003
)
Recent advances in molecular genetics of glaucoma
.
Mol. Cell. Biochem.
,
253
,
223
231
.

2

Thylefors
,
B.
,
Negrel
,
A.D.
,
Pararajasegaram
,
R.
and
Dadzie
,
K.Y.
(
1995
)
Global data on blindness
.
Bull. World. Health. Organ.
,
73
,
115
121
.

3

Tamm
,
E.R.
,
Braunger
,
B.M.
and
Fuchshofer
,
R.
(
2015
)
Intraocular pressure and the mechanisms involved in resistance of the aqueous humor flow in the trabecular meshwork outflow pathways
.
Prog. Mol. Biol. Transl. Sci.
,
134
,
301
314
.

4

deLuise
,
V.P.
and
Anderson
,
D.R.
(
1983
)
Primary infantile glaucoma (congenital glaucoma)
.
Surv. Ophthalmol.
,
28
,
1
19
.

5

Souma
,
T.
,
Tompson
,
S.W.
,
Thomson
,
B.R.
,
Siggs
,
O.M.
,
Kizhatil
,
K.
,
Yamaguchi
,
S.
,
Feng
,
L.
,
Limviphuvadh
,
V.
,
Whisenhunt
,
K.N.
,
Maurer-Stroh
,
S.
et al. (
2016
)
Angiopoietin receptor TEK mutations underlie primary congenital glaucoma with variable expressivity
.
J. Clin. Invest.
,
126
,
2575
2587
.

6

Stoilov
,
I.
,
Akarsu
,
A.N.
and
Sarfarazi
,
M.
(
1997
)
Identification of three different truncating mutations in cytochrome P4501B1 (CYP1B1) as the principal cause of primary congenital glaucoma (Buphthalmos) in families linked to the GLC3A locus on chromosome 2p21
.
Hum. Mol. Genet.
,
6
,
641
647
.

7

Ali
,
M.
,
McKibbin
,
M.
,
Booth
,
A.
,
Parry
,
D.A.
,
Jain
,
P.
,
Riazuddin
,
S.A.
,
Hejtmancik
,
J.F.
,
Khan
,
S.N.
,
Firasat
,
S.
,
Shires
,
M.
et al. (
2009
)
Null mutations in LTBP2 cause primary congenital glaucoma
.
Am. J. Hum. Genet.
,
84
,
664
671
.

8

Narooie-Nejad
,
M.
,
Paylakhi
,
S.H.
,
Shojaee
,
S.
,
Fazlali
,
Z.
,
Rezaei Kanavi
,
M.
,
Nilforushan
,
N.
,
Yazdani
,
S.
,
Babrzadeh
,
F.
,
Suri
,
F.
,
Ronaghi
,
M.
et al. (
2009
)
Loss of function mutations in the gene encoding latent transforming growth factor beta binding protein 2, LTBP2, cause primary congenital glaucoma
.
Hum. Mol. Genet.
,
18
,
3969
3977
.

9

Sakurada
,
Y.
and
Mabuchi
,
F.
(
2015
)
Advances in glaucoma genetics
.
Prog. Brain. Res.
,
220
,
107
126
.

10

Stone
,
E.M.
,
Fingert
,
J.H.
,
Alward
,
W.L.
,
Nguyen
,
T.D.
,
Polansky
,
J.R.
,
Sunden
,
S.L.
,
Nishimura
,
D.
,
Clark
,
A.F.
,
Nystuen
,
A.
,
Nichols
,
B.E.
et al. (
1997
)
Identification of a gene that causes primary open angle glaucoma
.
Science
,
275
,
668
670
.

11

Monemi
,
S.
,
Spaeth
,
G.
,
DaSilva
,
A.
,
Popinchalk
,
S.
,
Ilitchev
,
E.
,
Liebmann
,
J.
,
Ritch
,
R.
,
Heon
,
E.
,
Crick
,
R.P.
,
Child
,
A.
et al. (
2005
)
Identification of a novel adult-onset primary open-angle glaucoma (POAG) gene on 5q22.1
.
Hum. Mol. Genet.
,
14
,
725
733
.

12

Rezaie
,
T.
,
Child
,
A.
,
Hitchings
,
R.
,
Brice
,
G.
,
Miller
,
L.
,
Coca-Prados
,
M.
,
Heon
,
E.
,
Krupin
,
T.
,
Ritch
,
R.
,
Kreutzer
,
D.
et al. (
2002
)
Adult-onset primary open-angle glaucoma caused by mutations in optineurin
.
Science
,
295
,
1077
1079
.

13

Pasutto
,
F.
,
Matsumoto
,
T.
,
Mardin
,
C.Y.
,
Sticht
,
H.
,
Brandstatter
,
J.H.
,
Michels-Rautenstrauss
,
K.
,
Weisschuh
,
N.
,
Gramer
,
E.
,
Ramdas
,
W.D.
,
van
Koolwijk
,
L.M.
et al. (
2009
)
Heterozygous NTF4 mutations impairing neurotrophin-4 signaling in patients with primary open-angle glaucoma
.
Am. J. Hum. Genet.
,
85
,
447
456
.

14

Fingert
,
J.H.
,
Robin
,
A.L.
,
Stone
,
J.L.
,
Roos
,
B.R.
,
Davis
,
L.K.
,
Scheetz
,
T.E.
,
Bennett
,
S.R.
,
Wassink
,
T.H.
,
Kwon
,
Y.H.
,
Alward
,
W.L.
et al. (
2011
)
Copy number variations on chromosome 12q14 in patients with normal tension glaucoma
.
Hum. Mol. Genet.
,
20
,
2482
2494
.

15

Ahram
,
D.F.
,
Alward
,
W.L.
and
Kuehn
,
M.H.
(
2015
)
The genetic mechanisms of primary angle closure glaucoma
.
Eye (Lond)
,
29
,
1251
1259
.

16

Amerasinghe
,
N.
,
Zhang
,
J.
,
Thalamuthu
,
A.
,
He
,
M.
,
Vithana
,
E.N.
,
Viswanathan
,
A.
,
Wong
,
T.Y.
,
Foster
,
P.J.
and
Aung
,
T.
(
2011
)
The heritability and sibling risk of angle closure in Asians
.
Ophthalmology
,
118
,
480
485
.

17

Shastry
,
B.S.
(
2013
)
Genetic susceptibility to primary angle closure glaucoma (PACG)
.
Discov. Med.
,
15
,
17
22
.

18

Yazdani
,
S.
,
Akbarian
,
S.
,
Pakravan
,
M.
and
Afrouzifar
,
M.
(
2015
)
Prevalence of angle closure in siblings of patients with primary angle-closure glaucoma
.
J. Glaucoma
,
24
,
149
153
.

19

Vithana
,
E.N.
,
Khor
,
C.C.
,
Qiao
,
C.
,
Nongpiur
,
M.E.
,
George
,
R.
,
Chen
,
L.J.
,
Do
,
T.
,
Abu-Amero
,
K.
,
Huang
,
C.K.
,
Low
,
S.
et al. (
2012
) Genome-wide association analyses identify three new susceptibility loci for primary angle closure glaucoma.
Nat. Genet.
,
44
,
1142
1146
.

20

Nongpiur
,
M.E.
,
Khor
,
C.C.
,
Jia
,
H.
,
Cornes
,
B.K.
,
Chen
,
L.J.
,
Qiao
,
C.
,
Nair
,
K.S.
,
Cheng
,
C.Y.
,
Xu
,
L.
,
George
,
R.
et al. (
2014
)
ABCC5, a gene that influences the anterior chamber depth, is associated with primary angle closure glaucoma
.
PLoS. Genet.
,
10
,
e1004089
.

21

Khor
,
C.C.
,
Do
,
T.
,
Jia
,
H.
,
Nakano
,
M.
,
George
,
R.
,
Abu-Amero
,
K.
,
Duvesh
,
R.
,
Chen
,
L.J.
,
Li
,
Z.
,
Nongpiur
,
M.E.
et al. (
2016
) Genome-wide association study identifies five new susceptibility loci for primary angle closure glaucoma.
Nat. Genet.
,
48
,
556
562
.

22

Wang
,
I.J.
,
Chiang
,
T.H.
,
Shih
,
Y.F.
,
Lu
,
S.C.
,
Lin
,
L.L.
,
Shieh
,
J.W.
,
Wang
,
T.H.
,
Samples
,
J.R.
and
Hung
,
P.T.
(
2006
)
The association of single nucleotide polymorphisms in the MMP-9 genes with susceptibility to acute primary angle closure glaucoma in Taiwanese patients
.
Mol. Vis.
,
12
,
1223
1232
.

23

Awadalla
,
M.S.
,
Thapa
,
S.S.
,
Burdon
,
K.P.
,
Hewitt
,
A.W.
and
Craig
,
J.E.
(
2011
)
The association of hepatocyte growth factor (HGF) gene with primary angle closure glaucoma in the Nepalese population
.
Mol. Vis.
,
17
,
2248
2254
.

24

Awadalla
,
M.S.
,
Thapa
,
S.S.
,
Hewitt
,
A.W.
,
Craig
,
J.E.
and
Burdon
,
K.P.
(
2013
)
Association of eNOS polymorphisms with primary angle-closure glaucoma
.
Invest. Ophthalmol. Vis. Sci.
,
54
,
2108
2114
.

25

Ayub
,
H.
,
Khan
,
M.I.
,
Micheal
,
S.
,
Akhtar
,
F.
,
Ajmal
,
M.
,
Shafique
,
S.
,
Ali
,
S.H.
, den
Hollander
,
A.I.
,
Ahmed
,
A.
and
Qamar
,
R.
(
2010
)
Association of eNOS and HSP70 gene polymorphisms with glaucoma in Pakistani cohorts
.
Mol. Vis.
,
16
,
18
25
.

26

Shi
,
H.
,
Zhang
,
J.
,
Zhu
,
R.
,
Hu
,
N.
,
Lu
,
H.
,
Yang
,
M.
,
Qin
,
B.
,
Shi
,
J.
and
Guan
,
H.
(
2016
)
Primary angle closure and sequence variants within microRNA binding sites of genes involved in eye development
.
PLoS. One.
,
11
,
e0166055
.

27

Aung
,
T.
,
Yong
,
V.H.
,
Lim
,
M.C.
,
Venkataraman
,
D.
,
Toh
,
J.Y.
,
Chew
,
P.T.
and
Vithana
,
E.N.
(
2008
)
Lack of association between the rs2664538 polymorphism in the MMP-9 gene and primary angle closure glaucoma in Singaporean subjects
.
J. Glaucoma
,
17
,
257
258
.

28

Gupta
,
S.
,
Bhaskar
,
P.K.
,
Bhardwaj
,
R.
,
Chandra
,
A.
,
Chaudhry
,
V.N.
,
Chaudhry
,
P.
,
Ali
,
A.
,
Mukherjee
,
A.
and
Mutsuddi
,
M.
(
2014
)
MTHFR C677T predisposes to POAG but not to PACG in a North Indian population: a case control study
.
PLoS. One.
,
9
,
e103063
.

29

Nilforoushan
,
N.
,
Aghapour
,
S.
,
Raoofian
,
R.
,
Saee Rad
,
S.
,
Greene
,
W.K.
,
Fakhraie
,
G.
and
Heidari
,
M.
(
2012
)
Lack of association between the C677T single nucleotide polymorphism of the MTHFR gene and glaucoma in Iranian patients
.
Acta. Med. Iran.
,
50
,
208
212
.

30

Congdon
,
N.G.
,
Youlin
,
Q.
,
Quigley
,
H.
,
Hung
,
P.T.
,
Wang
,
T.H.
,
Ho
,
T.C.
and
Tielsch
,
J.M.
(
1997
)
Biometry and primary angle-closure glaucoma among Chinese, white, and black populations
.
Ophthalmology
,
104
,
1489
1495
.

31

Foster
,
P.J.
,
Alsbirk
,
P.H.
,
Baasanhu
,
J.
,
Munkhbayar
,
D.
,
Uranchimeg
,
D.
and
Johnson
,
G.J.
(
1997
)
Anterior chamber depth in Mongolians: variation with age, sex, and method of measurement
.
Am. J. Ophthalmol
,
124
,
53
60
.

32

Salmon
,
J.F.
(
1999
)
Predisposing factors for chronic angle-closure glaucoma
.
Prog. Retin. Eye. Res.
,
18
,
121
132
.

33

Foster
,
P.J.
,
Buhrmann
,
R.
,
Quigley
,
H.A.
and
Johnson
,
G.J.
(
2002
)
The definition and classification of glaucoma in prevalence surveys
.
Br. J. Ophthalmol.
,
86
,
238
242
.

34

Zhang
,
L.
,
Mabwi
,
H.A.
,
Palange
,
N.J.
,
Jia
,
R.
,
Ma
,
J.
,
Bah
,
F.B.
,
Sah
,
R.K.
,
Li
,
D.
,
Wang
,
D.
,
Bah
,
F.B.
et al. (
2015
)
Expression patterns and potential biological roles of Dip2a
.
PLoS. One.
,
10
,
e0143284
.

35

Kircher
,
M.
,
Witten
,
D.M.
,
Jain
,
P.
,
O'Roak
,
B.J.
,
Cooper
,
G.M.
and
Shendure
,
J.
(
2014
)
A general framework for estimating the relative pathogenicity of human genetic variants
.
Nat. Genet.
,
46
,
310
315
.

36

Knobloch
,
W.
and
Layer
,
J.
(
1971
)
Retinal detachment and encephalocele
.
J. Pediatr. Ophthalmol. Strabismus.
,
8
,
181
184
.

37

Passos-Bueno
,
M.R.
,
Suzuki
,
O.T.
,
Armelin-Correa
,
L.M.
,
Sertie
,
A.L.
,
Errera
,
F.I.
,
Bagatini
,
K.
,
Kok
,
F.
and
Leite
,
K.R.
(
2006
)
Mutations in collagen 18A1 and their relevance to the human phenotype
.
An. Acad. Bras. Cienc.
,
78
,
123
131
.

38

Suzuki
,
O.T.
,
Sertie
,
A.L.
,
Der Kaloustian
,
V.M.
,
Kok
,
F.
,
Carpenter
,
M.
,
Murray
,
J.
,
Czeizel
,
A.E.
,
Kliemann
,
S.E.
,
Rosemberg
,
S.
,
Monteiro
,
M.
et al. (
2002
)
Molecular analysis of collagen XVIII reveals novel mutations, presence of a third isoform, and possible genetic heterogeneity in Knobloch syndrome
.
Am. J. Hum. Genet.
,
71
,
1320
1329
.

39

Joyce
,
S.
,
Tee
,
L.
,
Abid
,
A.
,
Khaliq
,
S.
,
Mehdi
,
S.Q.
and
Maher
,
E.R.
(
2010
)
Locus heterogeneity and Knobloch syndrome
.
Am. J. Med. Genet. A.
,
152A
,
2880
2881
.

40

Khaliq
,
S.
,
Abid
,
A.
,
White
,
D.R.
,
Johnson
,
C.A.
,
Ismail
,
M.
,
Khan
,
A.
,
Ayub
,
Q.
,
Sultana
,
S.
,
Maher
,
E.R.
and
Mehdi
,
S.Q.
(
2007
)
Mapping of a novel type III variant of Knobloch syndrome (KNO3) to chromosome 17q11.2
.
Am. J. Med. Genet.
A.,
143A
,
2768
2774
.

41

Paisa'n-Ruiz
,
C.
,
Scopes
,
G.
,
Lee
,
P.
and
Houlden
,
H.
(
2009
)
Homozygosity mapping through whole genome analysis identifies a COL18A1 mutation in an Indian family presenting with an autosomal recessive neurological disorder
.
Am. J. Med. Genet. Part B.
,
150B
,
993
997
.

42

Marneros
,
A.G.
and
Olsen
,
B.R.
(
2003
)
Age-dependent iris abnormalities in collagen XVIII/endostatin deficient mice with similarities to human pigment dispersion syndrome
.
Invest. Ophthalmol. Vis. Sci.
,
44
,
2367
2372
.

43

Ylikarppa
,
R.
,
Eklund
,
L.
,
Sormunen
,
R.
,
Kontiola
,
A.I.
,
Utriainen
,
A.
,
Maatta
,
M.
,
Fukai
,
N.
,
Olsen
,
B.R.
and
Pihlajaniemi
,
T.
(
2003
)
Lack of type XVIII collagen results in anterior ocular defects
.
FASEB. J.
,
17
,
2257
2259
.

44

Maatta
,
M.
,
Heljasvaara
,
R.
,
Sormunen
,
R.
,
Pihlajaniemi
,
T.
,
Autio-Harmainen
,
H.
and
Tervo
,
T.
(
2006
)
Differential expression of collagen types XVIII/endostatin and XV in normal, keratoconus, and scarred human corneas
.
Cornea
,
25
,
341
349
.

45

Goel
,
R.
,
Murthy
,
K.R.
,
Srikanth
,
S.M.
,
Pinto
,
S.M.
,
Bhattacharjee
,
M.
,
Kelkar
,
D.S.
,
Madugundu
,
A.K.
,
Dey
,
G.
,
Mohan
,
S.S.
,
Krishna
,
V.
et al. (
2013
)
Characterizing the normal proteome of human ciliary body
.
Clin. Proteomics.
,
10
,
9
.

46

Zhao
,
X.
,
Pearson
,
K.E.
,
Stephan
,
D.A.
and
Russell
,
P.
(
2003
)
Effects of prostaglandin analogues on human ciliary muscle and trabecular meshwork cells
.
Invest. Ophthalmol. Vis. Sci.
,
44
,
1945
1952
.

47

Zhao
,
X.
,
Ramsey
,
K.E.
,
Stephan
,
D.A.
and
Russell
,
P.
(
2004
)
Gene and protein expression changes in human trabecular meshwork cells treated with transforming growth factor-beta
.
Invest. Ophthalmol. Vis. Sci.
,
45
,
4023
4034
.

48

Liton
,
P.B.
,
Luna
,
C.
,
Challa
,
P.
,
Epstein
,
D.L.
and
Gonzalez
,
P.
(
2006
)
Genome-wide expression profile of human trabecular meshwork cultured cells, nonglaucomatous and primary open angle glaucoma tissue
.
Mol. Vis.
,
12
,
774
790
.

49

Paylakhi
,
S.H.
,
Fan
,
J.B.
,
Mehrabian
,
M.
,
Sadeghizadeh
,
M.
,
Yazdani
,
S.
,
Katanforoush
,
A.
,
Kanavi
,
M.R.
,
Ronaghi
,
M.
and
Elahi
,
E.
(
2011
)
Effect of PITX2 knockdown on transcriptome of primary human trabecular meshwork cell cultures
.
Mol. Vis.
,
17
,
1209
1221
.

50

Bonomi
,
L.
(
2002
)
Epidemiology of angle-closure glaucoma
.
Acta Ophthalmol. Scand. Suppl.
,
236
,
11
13
.

51

Quigley
,
H.A.
and
Broman
,
A.T.
(
2006
)
The number of people with glaucoma worldwide in 2010 and 2020
.
Br. J. Ophthalmol.
,
90
,
262
267
.

52

Cedrone
,
C.
,
Mancino
,
R.
,
Cerulli
,
A.
,
Cesareo
,
M.
and
Nucci
,
C.
(
2008
)
Epidemiology of primary glaucoma: prevalence, incidence, and blinding effects
.
Prog. Brain. Res.
,
173
,
3
14
.

53

Dandona
,
L.
,
Dandona
,
R.
,
Mandal
,
P.
,
Srinivas
,
M.
,
John
,
R.K.
,
McCarty
,
C.A.
and
Rao
,
G.N.
(
2000
)
Angle-closure glaucoma in an urban population in southern India. The Andhra Pradesh eye disease study
. Ophthalmology,
107
,
1710
1716
.

54

Foster
,
P.J.
,
Baasanhu
,
J.
,
Alsbirk
,
P.H.
,
Munkhbayar
,
D.
,
Uranchimeg
,
D.
and
Johnson
,
G.J.
(
1996
)
Glaucoma in Mongolia. A population-based survey in Hovsgol province, northern Mongolia
.
Arch. Ophthalmol.
,
114
,
1235
1241
.

55

Foster
,
P.J.
and
Johnson
,
G.J.
(
2001
)
Glaucoma in China: how big is the problem?
Br. J. Ophthalmol.
,
85
,
1277
1282
.

56

Foster
,
P.J.
,
Oen
,
F.T.
,
Machin
,
D.
,
Ng
,
T.P.
,
Devereux
,
J.G.
,
Johnson
,
G.J.
,
Khaw
,
P.T.
and
Seah
,
S.K.
(
2000
)
The prevalence of glaucoma in Chinese residents of Singapore: a cross-sectional population survey of the Tanjong Pagar district
.
Arch. Ophthalmol.
,
118
,
1105
1111
.

57

Oh
,
S.P.
,
Kamagata
,
Y.
,
Muragaki
,
Y.
,
Timmons
,
S.
,
Ooshima
,
A.
and
Olsen
,
B.R.
(
1994
)
Isolation and sequencing of cDNAs for proteins with multiple domains of Gly-Xaa-Yaa repeats identify a distinct family of collagenous proteins
.
Proc. Natl. Acad. Sci. U S A
,
91
,
4229
4233
.

58

Oh
,
S.P.
,
Warman
,
M.L.
,
Seldin
,
M.F.
,
Cheng
,
S.D.
,
Knoll
,
J.H.
,
Timmons
,
S.
and
Olsen
,
B.R.
(
1994
)
Cloning of cDNA and genomic DNA encoding human type XVIII collagen and localization of the alpha 1(XVIII) collagen gene to mouse chromosome 10 and human chromosome 21
.
Genomics
,
19
,
494
499
.

59

Elamaa
,
H.
,
Snellman
,
A.
,
Rehn
,
M.
,
Autio-Harmainen
,
H.
and
Pihlajaniemi
,
T.
(
2003
)
Characterization of the human type XVIII collagen gene and proteolytic processing and tissue location of the variant containing a frizzled motif
.
Matrix. Biol.
,
22
,
427
442
.

60

Nicolae
,
C.
and
Olsen
,
B.R.
(
2010
)
Unexpected matrix diseases and novel therapeutic strategies
.
Cell. Tissue. Res.
,
339
,
155
165
.

61

Rasi
,
K.
,
Piuhola
,
J.
,
Czabanka
,
M.
,
Sormunen
,
R.
,
Ilves
,
M.
,
Leskinen
,
H.
,
Rysa
,
J.
,
Kerkela
,
R.
,
Janmey
,
P.
,
Heljasvaara
,
R.
et al. (
2010
)
Collagen XV is necessary for modeling of the extracellular matrix and its deficiency predisposes to cardiomyopathy
.
Circ. Res.
,
107
,
1241
1252
.

62

Iozzo
,
R.V.
(
2005
)
Basement membrane proteoglycans: from cellar to ceiling
.
Nat. Rev. Mol. Cell. Biol.
,
6
,
646
656
.

63

Li
,
D.
,
Clark
,
C.C.
and
Myers
,
J.C.
(
2000
)
Basement membrane zone type XV collagen is a disulfide-bonded chondroitin sulfate proteoglycan in human tissues and cultured cells
.
J. Biol. Chem.
,
275
,
22339
22347
.

64

Dhanabal
,
M.
,
Ramchandran
,
R.
,
Waterman
,
M.J.
,
Lu
,
H.
,
Knebelmann
,
B.
,
Segal
,
M.
and
Sukhatme
,
V.P.
(
1999
)
Endostatin induces endothelial cell apoptosis
.
J. Biol. Chem.
,
274
,
11721
11726
.

65

Dixelius
,
J.
,
Larsson
,
H.
,
Sasaki
,
T.
,
Holmqvist
,
K.
,
Lu
,
L.
,
Engstrom
,
A.
,
Timpl
,
R.
,
Welsh
,
M.
and
Claesson-Welsh
,
L.
(
2000
)
Endostatin-induced tyrosine kinase signaling through the Shb adaptor protein regulates endothelial cell apoptosis
.
Blood
,
95
,
3403
3411
.

66

Schmidt
,
A.
,
Wenzel
,
D.
,
Ferring
,
I.
,
Kazemi
,
S.
,
Sasaki
,
T.
,
Hescheler
,
J.
,
Timpl
,
R.
,
Addicks
,
K.
,
Fleischmann
,
B.K.
and
Bloch
,
W.
(
2004
)
Influence of endostatin on embryonic vasculo- and angiogenesis
.
Dev. Dyn.
,
230
,
468
480
.

67

O'Reilly
,
M.S.
,
Boehm
,
T.
,
Shing
,
Y.
,
Fukai
,
N.
,
Vasios
,
G.
,
Lane
,
W.S.
,
Flynn
,
E.
,
Birkhead
,
J.R.
,
Olsen
,
B.R.
and
Folkman
,
J.
(
1997
)
Endostatin: an endogenous inhibitor of angiogenesis and tumor growth
.
Cell
,
88
,
277
285
.

68

Katoh
,
M.
(
2001
)
Molecular cloning and characterization of MFRP, a novel gene encoding a membrane-type Frizzled-related protein
.
Biochem. Biophys. Res. Commun.
,
282
,
116
123
.

69

Tomaselli
,
S.
,
Megiorni
,
F.
,
Lin
,
L.
,
Mazzilli
,
M.C.
,
Gerrelli
,
D.
,
Majore
,
S.
,
Grammatico
,
P.
and
Achermann
,
J.C.
(
2011
)
Human RSPO1/R-spondin1 is expressed during early ovary development and augments beta-catenin signaling
.
PLoS. One.
,
6
,
e16366
.

70

Hao
,
H.X.
,
Xie
,
Y.
,
Zhang
,
Y.
,
Charlat
,
O.
,
Oster
,
E.
,
Avello
,
M.
,
Lei
,
H.
,
Mickanin
,
C.
,
Liu
,
D.
,
Ruffner
,
H.
et al. (
2012
)
ZNRF3 promotes Wnt receptor turnover in an R-spondin-sensitive manner
.
Nature
,
485
,
195
200
.

71

Fuhrmann
,
S.
(
2008
)
Wnt signaling in eye organogenesis
.
Organogenesis
,
4
,
60
67
.

72

Nakagawa
,
S.
(
2004
)
Role of Wnt signaling in eye organogenesis
.
Seikagaku
,
76
,
1573
1576
.

73

Huang
,
W.
,
Fan
,
Q.
,
Wang
,
W.
,
Zhou
,
M.
,
Laties
,
A.M.
and
Zhang
,
X.
(
2013
)
Collagen: a potential factor involved in the pathogenesis of glaucoma
.
Med. Sci. Monit. Basic. Res.
,
19
,
237
240
.

74

Vithana
,
E.N.
,
Aung
,
T.
,
Khor
,
C.C.
,
Cornes
,
B.K.
,
Tay
,
W.T.
,
Sim
,
X.
,
Lavanya
,
R.
,
Wu
,
R.
,
Zheng
,
Y.
,
Hibberd
,
M.L.
et al. (
2011
)
Collagen-related genes influence the glaucoma risk factor, central corneal thickness
.
Hum. Mol. Genet.
,
20
,
649
658
.

75

Aihara
,
M.
,
Lindsey
,
J.D.
and
Weinreb
,
R.N.
(
2003
)
Ocular hypertension in mice with a targeted type I collagen mutation
.
Invest. Ophthalmol. Vis. Sci.
,
44
,
1581
1585
.

76

Sawaguchi
,
S.
,
Yue
,
B.Y.
,
Fukuchi
,
T.
,
Abe
,
H.
,
Suda
,
K.
,
Kaiya
,
T.
and
Iwata
,
K.
(
1999
)
Collagen fibrillar network in the optic nerve head of normal monkey eyes and monkey eyes with laser-induced glaucoma--a scanning electron microscopic study
.
Curr. Eye. Res.
,
18
,
143
149
.

77

Suri
,
F.
,
Yazdani
,
S.
and
Elahi
,
E.
(
2015
)
Glaucoma in iran and contributions of studies in iran to the understanding of the etiology of glaucoma
.
J. Ophthalmic.Vis. Res.
,
10
,
68
76
.

78

Inamori
,
Y.
,
Ota
,
M.
,
Inoko
,
H.
,
Okada
,
E.
,
Nishizaki
,
R.
,
Shiota
,
T.
,
Mok
,
J.
,
Oka
,
A.
,
Ohno
,
S.
and
Mizuki
,
N.
(
2007
)
The COL1A1 gene and high myopia susceptibility in Japanese
.
Hum. Genet.
,
122
,
151
157
.

79

Zhang
,
D.
,
Shi
,
Y.
,
Gong
,
B.
,
He
,
F.
,
Lu
,
F.
,
Lin
,
H.
,
Wu
,
Z.
,
Cheng
,
J.
,
Chen
,
B.
,
Liao
,
S.
et al. (
2011
)
An association study of the COL1A1 gene and high myopia in a Han Chinese population
.
Mol. Vis.
,
17
,
3379
3383
.

80

Ahram
,
D.F.
,
Cook
,
A.C.
,
Kecova
,
H.
,
Grozdanic
,
S.D.
and
Kuehn
,
M.H.
(
2014
)
Identification of genetic loci associated with primary angle-closure glaucoma in the basset hound
.
Mol. Vis.
,
20
,
497
510
.

81

Cong
,
Y.
,
Guo
,
X.
,
Liu
,
X.
,
Cao
,
D.
,
Jia
,
X.
,
Xiao
,
X.
,
Li
,
S.
,
Fang
,
S.
and
Zhang
,
Q.
(
2009
)
Association of the single nucleotide polymorphisms in the extracellular matrix metalloprotease-9 gene with PACG in southern China
.
Mol. Vis.
,
15
,
1412
1417
.

82

Micheal
,
S.
,
Yousaf
,
S.
,
Khan
,
M.I.
,
Akhtar
,
F.
,
Islam
,
F.
,
Khan
,
W.A.
,
den
Hollander
,
A.I.
,
Qamar
,
R.
and
Ahmed
,
A.
(
2013
)
Polymorphisms in matrix metalloproteinases MMP1 and MMP9 are associated with primary open-angle and angle closure glaucoma in a Pakistani population
.
Mol. Vis.
,
19
,
441
447
.

83

Awadalla
,
M.S.
,
Burdon
,
K.P.
,
Kuot
,
A.
,
Hewitt
,
A.W.
and
Craig
,
J.E.
(
2011
)
Matrix metalloproteinase-9 genetic variation and primary angle closure glaucoma in a Caucasian population
.
Mol. Vis.
,
17
,
1420
1424
.

84

Micheal
,
S.
,
Qamar
,
R.
,
Akhtar
,
F.
,
Khan
,
M.I.
,
Khan
,
W.A.
and
Ahmed
,
A.
(
2009
)
MTHFR gene C677T and A1298C polymorphisms and homocysteine levels in primary open angle and primary closed angle glaucoma
.
Mol. Vis.
,
15
,
2268
2278
.

85

Safari
,
I.
,
Suri
,
F.
,
Haji-Seyed-Javadi
,
R.
,
Yazdani
,
S.
and
Elahi
,
E.
(
2016
)
The p.Gly61Glu Mutation in CYP1B1 Affects the Extracellular Matrix in Glaucoma Patients
.
Ophthalmic. Res.
,
56
,
98
103
.

86

Chua
,
J.
,
Seet
,
L.F.
,
Jiang
,
Y.
,
Su
,
R.
,
Htoon
,
H.M.
,
Charlton
,
A.
,
Aung
,
T.
and
Wong
,
T.T.
(
2008
)
Increased SPARC expression in primary angle closure glaucoma iris
.
Mol. Vis.
,
14
,
1886
1892
.

87

Abecasis
,
G.R.
,
Cherny
,
S.S.
,
Cookson
,
W.O.
and
Cardon
,
L.R.
(
2002
)
Merlin--rapid analysis of dense genetic maps using sparse gene flow trees
.
Nat. Genet.
,
30
,
97
101
.

This article is published and distributed under the terms of the Oxford University Press, Standard Journals Publication Model (https://academic.oup.com/journals/pages/open_access/funder_policies/chorus/standard_publication_model)

Supplementary data