Specifications of the ACMG/AMP variant curation guidelines for myocilin: Recommendations from the clingen glaucoma expert panel.

The standardization of variant curation criteria is essential for accurate interpretation of genetic results and clinical care of patients. The variant curation guidelines developed by the American College of Medical Genetics and Genomics (ACMG) and the Association for Molecular Pathology (AMP) in 2015 are widely used but are not gene specific. To address this issue, the Clinical Genome Resource (ClinGen) Variant Curation Expert Panels (VCEP) have been tasked with developing gene-specific variant curation guidelines. The Glaucoma VCEP was created to develop rule specifications for genes associated with primary glaucoma, including myocilin (MYOC), the most common cause of Mendelian glaucoma. Of the 28 ACMG/AMP criteria, the Glaucoma VCEP adapted 15 rules to MYOC and determined 13 rules not applicable. Key specifications included determining minor allele frequency thresholds, developing an approach to counting probands and segregations, and reviewing functional assays. The rules were piloted on 81 variants and led to a change in classification in 40% of those that were classified in ClinVar, with functional evidence influencing the classification of 18 variants. The standardized variant curation guidelines for MYOC provide a framework for the consistent application of the rules between laboratories, to improve MYOC genetic testing in the management of glaucoma.

Trio-based exome sequencing arrests de novo mutations in early-onset high myopia.

The etiology of the highly myopic condition has been unclear for decades. We investigated the genetic contributions to early-onset high myopia (EOHM), which is defined as having a refraction of less than or equal to -6 diopters before the age of 6, when children are less likely to be exposed to high educational pressures. Trios (two nonmyopic parents and one child) were examined to uncover pathogenic mutations using whole-exome sequencing. We identified parent-transmitted biallelic mutations or de novo mutations in as-yet-unknown or reported genes in 16 probands. Interestingly, an increased rate of de novo mutations was identified in the EOHM patients. Among the newly identified candidate genes, a BSG mutation was identified in one EOHM proband. Expanded screening of 1,040 patients found an additional four mutations in the same gene. Then, we generated Bsg mutant mice to further elucidate the functional impact of this gene and observed typical myopic phenotypes, including an elongated axial length. Using a trio-based exonic screening study in EOHM, we deciphered a prominent role for de novo mutations in EOHM patients without myopic parents. The discovery of a disease gene, BSG, provides insights into myopic development and its etiology, which expands our current understanding of high myopia and might be useful for future treatment and prevention.

Meta-analysis of Genome-Wide Association Studies Identifies Novel Loci Associated With Optic Disc Morphology.

Primary open-angle glaucoma is the most common optic neuropathy and an important cause of irreversible blindness worldwide. The optic nerve head or optic disc is divided in two parts: a central cup (without nerve fibers) surrounded by the neuroretinal rim (containing axons of the retinal ganglion cells). The International Glaucoma Genetics Consortium conducted a meta-analysis of genome-wide association studies consisting of 17,248 individuals of European ancestry and 6,841 individuals of Asian ancestry. The outcomes of the genome-wide association studies were disc area and cup area. These specific measurements describe optic nerve morphology in another way than the vertical cup-disc ratio, which is a clinically used measurement, and may shed light on new glaucoma mechanisms. We identified 10 new loci associated with disc area (CDC42BPA, F5, DIRC3, RARB, ABI3BP, DCAF4L2, ELP4, TMTC2, NR2F2, and HORMAD2) and another 10 new loci associated with cup area (DHRS3, TRIB2, EFEMP1, FLNB, FAM101, DDHD1, ASB7, KPNB1, BCAS3, and TRIOBP). The new genes participate in a number of pathways and future work is likely to identify more functions related to the pathogenesis of glaucoma.

A genome-wide association study of intra-ocular pressure suggests a novel association in the gene FAM125B in the TwinsUK cohort.

Glaucoma is a major cause of blindness in the world. To date, common genetic variants associated with glaucoma only explain a small proportion of its heritability. We performed a genome-wide association study of intra-ocular pressure (IOP), an underlying endophenotype for glaucoma. The discovery phase of the study was carried out in the TwinsUK cohort (N = 2774) analyzing association between IOP and single nucleotide polymorphisms (SNPs) imputed to HapMap2. The results were validated in 12 independent replication cohorts of European ancestry (combined N = 22 789) that were a part of the International Glaucoma Genetics Consortium. Expression quantitative trait locus (eQTL) analyses of the significantly associated SNPs were performed using data from the Multiple Tissue Human Expression Resource (MuTHER) Study. In the TwinsUK cohort, IOP was significantly associated with a number of SNPs at 9q33.3 (P = 3.48 × 10(-8) for rs2286885, the most significantly associated SNP at this locus), within the genomic sequence of the FAM125B gene. Independent replication in a composite panel of 12 cohorts revealed consistent direction of effect and significant association (P = 0.003, for fixed-effect meta-analysis). Suggestive evidence for an eQTL effect of rs2286885 was observed for one of the probes targeting the coding region of the FAM125B gene. This gene codes for a component of a membrane complex involved in vesicular trafficking process, a function similar to that of the Caveolin genes (CAV1 and CAV2) which have previously been associated with primary open-angle glaucoma. This study suggests a novel association between SNPs in FAM125B and IOP in the TwinsUK cohort, though further studies to elucidate the functional role of this gene in glaucoma are necessary.

A meta-analysis of genome-wide association studies for adiponectin levels in East Asians identifies a novel locus near WDR11-FGFR2.

Blood levels of adiponectin, an adipocyte-secreted protein correlated with metabolic and cardiovascular risks, are highly heritable. Genome-wide association (GWA) studies for adiponectin levels have identified 14 loci harboring variants associated with blood levels of adiponectin. To identify novel adiponectin-associated loci, particularly those of importance in East Asians, we conducted a meta-analysis of GWA studies for adiponectin in 7827 individuals, followed by two stages of replications in 4298 and 5954 additional individuals. We identified a novel adiponectin-associated locus on chromosome 10 near WDR11-FGFR2 (P = 3.0 × 10(-14)) and provided suggestive evidence for a locus on chromosome 12 near OR8S1-LALBA (P = 1.2 × 10(-7)). Of the adiponectin-associated loci previously described, we confirmed the association at CDH13 (P = 6.8 × 10(-165)), ADIPOQ (P = 1.8 × 10(-22)), PEPD (P = 3.6 × 10(-12)), CMIP (P = 2.1 × 10(-10)), ZNF664 (P = 2.3 × 10(-7)) and GPR109A (P = 7.4 × 10(-6)). Conditional analysis at ADIPOQ revealed a second signal with suggestive evidence of association only after conditioning on the lead SNP (Pinitial = 0.020; Pconditional = 7.0 × 10(-7)). We further confirmed the independence of two pairs of closely located loci (<2 Mb) on chromosome 16 at CMIP and CDH13, and on chromosome 12 at GPR109A and ZNF664. In addition, the newly identified signal near WDR11-FGFR2 exhibited evidence of association with triglycerides (P = 3.3 × 10(-4)), high density lipoprotein cholesterol (HDL-C, P = 4.9 × 10(-4)) and body mass index (BMI)-adjusted waist-hip ratio (P = 9.8 × 10(-3)). These findings improve our knowledge of the genetic basis of adiponectin variation, demonstrate the shared allelic architecture for adiponectin with lipids and central obesity and motivate further studies of underlying mechanisms.

Mutations in SCO2 are associated with autosomal-dominant high-grade myopia.

Myopia, or near-sightedness, is an ocular refractive error of unfocused image quality in front of the retinal plane. Individuals with high-grade myopia (dioptric power greater than -6.00) are predisposed to ocular morbidities such as glaucoma, retinal detachment, and myopic maculopathy. Nonsyndromic, high-grade myopia is highly heritable, and to date multiple gene loci have been reported. We performed exome sequencing in 4 individuals from an 11-member family of European descent from the United States. Affected individuals had a mean dioptric spherical equivalent of -22.00 sphere. A premature stop codon mutation c.157C>T (p.Gln53*) cosegregating with disease was discovered within SCO2 that maps to chromosome 22q13.33. Subsequent analyses identified three additional mutations in three highly myopic unrelated individuals (c.341G>A, c.418G>A, and c.776C>T). To determine differential gene expression in a developmental mouse model, we induced myopia by applying a -15.00D lens over one eye. Messenger RNA levels of SCO2 were significantly downregulated in myopic mouse retinae. Immunohistochemistry in mouse eyes confirmed SCO2 protein localization in retina, retinal pigment epithelium, and sclera. SCO2 encodes for a copper homeostasis protein influential in mitochondrial cytochrome c oxidase activity. Copper deficiencies have been linked with photoreceptor loss and myopia with increased scleral wall elasticity. Retinal thinning has been reported with an SC02 variant. Human mutation identification with support from an induced myopic animal provides biological insights of myopic development.

Nine loci for ocular axial length identified through genome-wide association studies, including shared loci with refractive error.

Refractive errors are common eye disorders of public health importance worldwide. Ocular axial length (AL) is the major determinant of refraction and thus of myopia and hyperopia. We conducted a meta-analysis of genome-wide association studies for AL, combining 12,531 Europeans and 8,216 Asians. We identified eight genome-wide significant loci for AL (RSPO1, C3orf26, LAMA2, GJD2, ZNRF3, CD55, MIP, and ALPPL2) and confirmed one previously reported AL locus (ZC3H11B). Of the nine loci, five (LAMA2, GJD2, CD55, ALPPL2, and ZC3H11B) were associated with refraction in 18 independent cohorts (n = 23,591). Differential gene expression was observed for these loci in minus-lens-induced myopia mouse experiments and human ocular tissues. Two of the AL genes, RSPO1 and ZNRF3, are involved in Wnt signaling, a pathway playing a major role in the regulation of eyeball size. This study provides evidence of shared genes between AL and refraction, but importantly also suggests that these traits may have unique pathways.

Genetic variants on chromosome 1q41 influence ocular axial length and high myopia.

As one of the leading causes of visual impairment and blindness, myopia poses a significant public health burden in Asia. The primary determinant of myopia is an elongated ocular axial length (AL). Here we report a meta-analysis of three genome-wide association studies on AL conducted in 1,860 Chinese adults, 929 Chinese children, and 2,155 Malay adults. We identified a genetic locus on chromosome 1q41 harboring the zinc-finger 11B pseudogene ZC3H11B showing genome-wide significant association with AL variation (rs4373767, ß = -0.16 mm per minor allele, P(meta) =2.69 × 10(-10)). The minor C allele of rs4373767 was also observed to significantly associate with decreased susceptibility to high myopia (per-allele odds ratio (OR) =0.75, 95% CI: 0.68-0.84, P(meta) =4.38 × 10(-7)) in 1,118 highly myopic cases and 5,433 controls. ZC3H11B and two neighboring genes SLC30A10 and LYPLAL1 were expressed in the human neural retina, retinal pigment epithelium, and sclera. In an experimental myopia mouse model, we observed significant alterations to gene and protein expression in the retina and sclera of the unilateral induced myopic eyes for the murine genes ZC3H11A, SLC30A10, and LYPLAL1. This supports the likely role of genetic variants at chromosome 1q41 in influencing AL variation and high myopia.

Common genetic determinants of intraocular pressure and primary open-angle glaucoma.

Intraocular pressure (IOP) is a highly heritable risk factor for primary open-angle glaucoma and is the only target for current glaucoma therapy. The genetic factors which determine IOP are largely unknown. We performed a genome-wide association study for IOP in 11,972 participants from 4 independent population-based studies in The Netherlands. We replicated our findings in 7,482 participants from 4 additional cohorts from the UK, Australia, Canada, and the Wellcome Trust Case-Control Consortium 2/Blue Mountains Eye Study. IOP was significantly associated with rs11656696, located in GAS7 at 17p13.1 (p=1.4×10(-8)), and with rs7555523, located in TMCO1 at 1q24.1 (p=1.6×10(-8)). In a meta-analysis of 4 case-control studies (total N = 1,432 glaucoma cases), both variants also showed evidence for association with glaucoma (p=2.4×10(-2) for rs11656696 and p=9.1×10(-4) for rs7555523). GAS7 and TMCO1 are highly expressed in the ciliary body and trabecular meshwork as well as in the lamina cribrosa, optic nerve, and retina. Both genes functionally interact with known glaucoma disease genes. These data suggest that we have identified two clinically relevant genes involved in IOP regulation.

Whole exome sequencing identifies a mutation for a novel form of corneal intraepithelial dyskeratosis.

BACKGROUND: Corneal intraepithelial dyskeratosis is an extremely rare condition. The classical form, affecting Native American Haliwa-Saponi tribe members, is called hereditary benign intraepithelial dyskeratosis (HBID). Herein, we present a new form of corneal intraepithelial dyskeratosis for which we identified the causative gene by using deep sequencing technology. METHODS AND RESULTS: A seven member Caucasian French family with two corneal intraepithelial dyskeratosis affected individuals (6-year-old proband and his mother) was ascertained. The proband presented with bilateral complete corneal opacification and dyskeratosis. Palmoplantar hyperkeratosis and laryngeal dyskeratosis were associated with the phenotype. Histopathology studies of cornea and vocal cord biopsies showed dyskeratotic keratinisation. Quantitative PCR ruled out 4q35 duplication, classically described in HBID cases. Next generation sequencing with mean coverage of 50× using the Illumina Hi Seq and whole exome capture processing was performed. Sequence reads were aligned, and screened for single nucleotide variants and insertion/deletion calls. In-house pipeline filtering analyses and comparisons with available databases were performed. A novel missense mutation M77T was discovered for the gene NLRP1 which maps to chromosome 17p13.2. This was a de novo mutation in the proband's mother, following segregation in the family, and not found in 738 control DNA samples. NLRP1 expression was determined in adult corneal epithelium. The amino acid change was found to destabilise significantly the protein structure. CONCLUSIONS: We describe a new corneal intraepithelial dyskeratosis and how we identified its causative gene. The NLRP1 gene product is implicated in inflammation, autoimmune disorders, and caspase mediated apoptosis. NLRP1 polymorphisms are associated with various diseases.

Genome-wide meta-analysis of five Asian cohorts identifies PDGFRA as a susceptibility locus for corneal astigmatism.

Corneal astigmatism refers to refractive abnormalities and irregularities in the curvature of the cornea, and this interferes with light being accurately focused at a single point in the eye. This ametropic condition is highly prevalent, influences visual acuity, and is a highly heritable trait. There is currently a paucity of research in the genetic etiology of corneal astigmatism. Here we report the results from five genome-wide association studies of corneal astigmatism across three Asian populations, with an initial discovery set of 4,254 Chinese and Malay individuals consisting of 2,249 cases and 2,005 controls. Replication was obtained from three surveys comprising of 2,139 Indians, an additional 929 Chinese children, and an independent 397 Chinese family trios. Variants in PDGFRA on chromosome 4q12 (lead SNP: rs7677751, allelic odds ratio = 1.26 (95% CI: 1.16-1.36), P(meta) = 7.87×10(-9)) were identified to be significantly associated with corneal astigmatism, exhibiting consistent effect sizes across all five cohorts. This highlights the potential role of variants in PDGFRA in the genetic etiology of corneal astigmatism across diverse Asian populations.

Mutation in collagen II alpha 1 isoforms delineates Stickler and Wagner syndrome phenotypes.

PURPOSE: Stickler syndrome is an arthro-ophthalmopathy with phenotypic overlap with Wagner syndrome. The common Stickler syndrome type I is inherited as an autosomal dominant trait, with causal mutations in collagen type II alpha 1 (COL2A1). Wagner syndrome is associated with mutations in versican (VCAN), which encodes for a chondroitin sulfate proteoglycan. A three-generation Caucasian family variably diagnosed with either syndrome was screened for sequence variants in the COL2A1 and VCAN genes. METHODS: Genomic DNA samples derived from saliva were collected from all family members (six affected and four unaffected individuals). Complete sequencing of COL2A1 and VCAN was performed on two affected individuals. Direct sequencing of remaining family members was conducted if the discovered variants followed segregation. RESULTS: A base-pair substitution (c.258C>A) in exon 2 of COL2A1 cosegregated with familial disease status. This known mutation occurs in a highly conserved site that causes a premature stop codon (p.C86X). The mutation was not seen in 1,142 ethnically matched control DNA samples. CONCLUSIONS: Premature stop codons in COL2A1 exon 2 lead to a Stickler syndrome type I ocular-only phenotype with few or no systemic manifestations. Mutation screening of COL2A1 exon 2 in families with autosomal dominant vitreoretinopathy is important for accurate clinical diagnosis.

Interrogation of the platelet-derived growth factor receptor alpha locus and corneal astigmatism in Australians of Northern European ancestry: results of a genome-wide association study.

PURPOSE: Corneal astigmatism is a common eye disorder characterized by irregularities in corneal curvature. Recently, the rs7677751 single nucleotide polymorphism (SNP) at the platelet-derived growth factor receptor alpha (PDGFRA) locus was found to be associated with corneal astigmatism in people of Asian ancestry. In the present study, we sought to replicate this finding and identify other genetic markers of corneal astigmatism in an Australian population of Northern European ancestry. METHODS: Data from two cohorts were included in this study. The first cohort consisted of 1,013 individuals who were part of the Western Australian Pregnancy Cohort (Raine) Study: 20-year follow-up Eye Study. The second cohort comprised 1,788 individuals of 857 twin families who were recruited through the Twins Eye Study in Tasmania and the Brisbane Adolescent Twin Study. Corneal astigmatism was calculated as the absolute difference between the keratometry readings in two meridians, and genotype data were extracted from genome-wide arrays. Initially, each cohort was analyzed separately, before being combined for meta- and subsequent genome-wide pathway analysis. RESULTS: Following meta-analysis, SNP rs7677751 at the PDGFRA locus had a combined p=0.32. No variant was found to be statistically significantly associated with corneal astigmatism at the genome-wide level (p<5.0×10(-8)). The SNP with strongest association was rs1164064 (p=1.86×10(-6)) on chromosome 3q13. Gene-based pathway analysis identified a significant association between the Gene Ontology "segmentation" (GO:0035282) pathway, corrected p=0.009. CONCLUSIONS: Our data suggest that the PDGFRA locus does not transfer a major risk of corneal astigmatism in people of Northern European ancestry. Better-powered studies are required to validate the novel putative findings of our study.

Genetic contributions to myopic refractive error: Insights from human studies and supporting evidence from animal models.

Genetic studies of both population-based and recruited affected patient cohorts have identified a number of genomic regions and candidate genes that may contribute to myopic development. Scientists have developed animal models of myopia, as collection of affected tissues from patents is impractical. Recent advances in whole exome sequencing technology show promise for further elucidation of disease causing variants as in the recent identification of rare variants within ZNF644 segregating with pathological myopia. We present a review of the current research trends and findings on genetic contributions to myopic refraction including candidate loci for myopic development and their genomic convergence with expression studies of animal models inducing myopic development.

Whole genome expression profiling of normal human fetal and adult ocular tissues.

To study growth and development of ocular tissues, gene expression patterns in normal human fetal versus adult eyes were compared. Human retina/retinal pigment epithelium, choroid, sclera, optic nerve* and cornea* tissues were dissected from fetal (24 week gestational age) (N = 9; *N = 6), and adult (N = 6) normal donor eyes. The Illumina(®) whole genome expression microarray platform was used to assess differential expression. Statistical significance for all comparisons was determined using the Benjamin and Hochberg False Discovery Rate (FDR, 5%). Significant gene expression fold changes > 1.5 were found in adult versus fetal retina/RPE (N = 1185), choroid (N = 6446), sclera (N = 1349), and cornea (N = 3872), but not optic nerve. Genes showing differential expression were assessed using Ingenuity Pathway Analysis (IPA) for enriched functions and canonical pathways. In all tissues, development, cell death/growth, cancer functions, and signaling canonical pathways were enriched. There was also a general trend of down-regulation of collagen genes in adult tissues.

Common genetic variants near the Brittle Cornea Syndrome locus ZNF469 influence the blinding disease risk factor central corneal thickness.

Central corneal thickness (CCT), one of the most highly heritable human traits (h(2) typically>0.9), is important for the diagnosis of glaucoma and a potential risk factor for glaucoma susceptibility. We conducted genome-wide association studies in five cohorts from Australia and the United Kingdom (total N = 5058). Three cohorts were based on individually genotyped twin collections, with the remaining two cohorts genotyped on pooled samples from singletons with extreme trait values. The pooled sample findings were validated by individual genotyping the pooled samples together with additional samples also within extreme quantiles. We describe methods for efficient combined analysis of the results from these different study designs. We have identified and replicated quantitative trait loci on chromosomes 13 and 16 for association with CCT. The locus on chromosome 13 (nearest gene FOXO1) had an overall meta-analysis p-value for all the individually genotyped samples of 4.6x10(-10). The locus on chromosome 16 was associated with CCT with p = 8.95x10(-11). The nearest gene to the associated chromosome 16 SNPs was ZNF469, a locus recently implicated in Brittle Cornea Syndrome (BCS), a very rare disorder characterized by abnormal thin corneas. Our findings suggest that in addition to rare variants in ZNF469 underlying CCT variation in BCS patients, more common variants near this gene may contribute to CCT variation in the general population.

Four novel Loci (19q13, 6q24, 12q24, and 5q14) influence the microcirculation in vivo.

There is increasing evidence that the microcirculation plays an important role in the pathogenesis of cardiovascular diseases. Changes in retinal vascular caliber reflect early microvascular disease and predict incident cardiovascular events. We performed a genome-wide association study to identify genetic variants associated with retinal vascular caliber. We analyzed data from four population-based discovery cohorts with 15,358 unrelated Caucasian individuals, who are members of the Cohort for Heart and Aging Research in Genomic Epidemiology (CHARGE) consortium, and replicated findings in four independent Caucasian cohorts (n  =  6,652). All participants had retinal photography and retinal arteriolar and venular caliber measured from computer software. In the discovery cohorts, 179 single nucleotide polymorphisms (SNP) spread across five loci were significantly associated (p<5.0×10(-8)) with retinal venular caliber, but none showed association with arteriolar caliber. Collectively, these five loci explain 1.0%-3.2% of the variation in retinal venular caliber. Four out of these five loci were confirmed in independent replication samples. In the combined analyses, the top SNPs at each locus were: rs2287921 (19q13; p  =  1.61×10(-25), within the RASIP1 locus), rs225717 (6q24; p = 1.25×10(-16), adjacent to the VTA1 and NMBR loci), rs10774625 (12q24; p  =  2.15×10(-13), in the region of ATXN2,SH2B3 and PTPN11 loci), and rs17421627 (5q14; p = 7.32×10(-16), adjacent to the MEF2C locus). In two independent samples, locus 12q24 was also associated with coronary heart disease and hypertension. Our population-based genome-wide association study demonstrates four novel loci associated with retinal venular caliber, an endophenotype of the microcirculation associated with clinical cardiovascular disease. These data provide further insights into the contribution and biological mechanisms of microcirculatory changes that underlie cardiovascular disease.

Large-scale multitrait genome-wide association analyses identify hundreds of glaucoma risk loci.

Glaucoma, a leading cause of irreversible blindness, is a highly heritable human disease. Previous genome-wide association studies have identified over 100 loci for the most common form, primary open-angle glaucoma. Two key glaucoma-associated traits also show high heritability: intraocular pressure and optic nerve head excavation damage quantified as the vertical cup-to-disc ratio. Here, since much of glaucoma heritability remains unexplained, we conducted a large-scale multitrait genome-wide association study in participants of European ancestry combining primary open-angle glaucoma and its two associated traits (total sample size over 600,000) to substantially improve genetic discovery power (263 loci). We further increased our power by then employing a multiancestry approach, which increased the number of independent risk loci to 312, with the vast majority replicating in a large independent cohort from 23andMe, Inc. (total sample size over 2.8 million; 296 loci replicated at P < 0.05, 240 after Bonferroni correction). Leveraging multiomics datasets, we identified many potential druggable genes, including neuro-protection targets likely to act via the optic nerve, a key advance for glaucoma because all existing drugs only target intraocular pressure. We further used Mendelian randomization and genetic correlation-based approaches to identify novel links to other complex traits, including immune-related diseases such as multiple sclerosis and systemic lupus erythematosus.

Genome-wide association identifies ATOH7 as a major gene determining human optic disc size.

Optic nerve assessment is important for many blinding diseases, with cup-to-disc ratio (CDR) assessments commonly used in both diagnosis and progression monitoring of glaucoma patients. Optic disc, cup, rim area and CDR measurements all show substantial variation between human populations and high heritability estimates within populations. To identify loci underlying these quantitative traits, we performed a genome-wide association study in two Australian twin cohorts and identified rs3858145, P=6.2x10(-10), near the ATOH7 gene as associated with the mean disc area. ATOH7 is known from studies in model organisms to play a key role in retinal ganglion cell formation. The association with rs3858145 was replicated in a cohort of UK twins, with a meta-analysis of the combined data yielding P=3.4x10(-10). Imputation further increased the evidence for association for several SNPs in and around ATOH7 (P=1.3x10(-10) to 4.3x10(-11), top SNP rs1900004). The meta-analysis also provided suggestive evidence for association for the cup area at rs690037, P=1.5x10(-7), in the gene RFTN1. Direct sequencing of ATOH7 in 12 patients with optic nerve hypoplasia, one of the leading causes of blindness in children, revealed two novel non-synonymous mutations (Arg65Gly, Ala47Thr) which were not found in 90 unrelated controls (combined Fisher's exact P=0.0136). Furthermore, the Arg65Gly variant was found to have very low frequency (0.00066) in an additional set of 672 controls.

Large scale international replication and meta-analysis study confirms association of the 15q14 locus with myopia. The CREAM consortium.

Myopia is a complex genetic disorder and a common cause of visual impairment among working age adults. Genome-wide association studies have identified susceptibility loci on chromosomes 15q14 and 15q25 in Caucasian populations of European ancestry. Here, we present a confirmation and meta-analysis study in which we assessed whether these two loci are also associated with myopia in other populations. The study population comprised 31 cohorts from the Consortium of Refractive Error and Myopia (CREAM) representing 4 different continents with 55,177 individuals; 42,845 Caucasians and 12,332 Asians. We performed a meta-analysis of 14 single nucleotide polymorphisms (SNPs) on 15q14 and 5 SNPs on 15q25 using linear regression analysis with spherical equivalent as a quantitative outcome, adjusted for age and sex. We calculated the odds ratio (OR) of myopia versus hyperopia for carriers of the top-SNP alleles using a fixed effects meta-analysis. At locus 15q14, all SNPs were significantly replicated, with the lowest P value 3.87 × 10(-12) for SNP rs634990 in Caucasians, and 9.65 × 10(-4) for rs8032019 in Asians. The overall meta-analysis provided P value 9.20 × 10(-23) for the top SNP rs634990. The risk of myopia versus hyperopia was OR 1.88 (95 % CI 1.64, 2.16, P < 0.001) for homozygous carriers of the risk allele at the top SNP rs634990, and OR 1.33 (95 % CI 1.19, 1.49, P < 0.001) for heterozygous carriers. SNPs at locus 15q25 did not replicate significantly (P value 5.81 × 10(-2) for top SNP rs939661). We conclude that common variants at chromosome 15q14 influence susceptibility for myopia in Caucasian and Asian populations world-wide.