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

Bi-Gaussian analysis reveals distinct education-related alterations in spherical equivalent and axial length-results from the Gutenberg Health Study

PURPOSE

The aim of this study is to investigate the distribution of spherical equivalent and axial length in the general population and to analyze the influence of education on spherical equivalent with a focus on ocular biometric parameters.

METHODS

The Gutenberg Health Study is a population-based cohort study in Mainz, Germany. Participants underwent comprehensive ophthalmologic examinations as part of the 5-year follow-up examination in 2012-2017 including genotyping. The spherical equivalent and axial length distributions were modeled with gaussian mixture models. Regression analysis (on person-individual level) was performed to analyze associations between biometric parameters and educational factors. Mendelian randomization analysis explored the causal effect between spherical equivalent, axial length, and education. Additionally, effect mediation analysis examined the link between spherical equivalent and education.

RESULTS

A total of 8532 study participants were included (median age: 57 years, 49% female). The distribution of spherical equivalent and axial length follows a bi-Gaussian function, partially explained by the length of education (i.e., < 11 years education vs. 11-20 years). Mendelian randomization indicated an effect of education on refractive error using a genetic risk score of education as an instrument variable (- 0.35 diopters per SD increase in the instrument, 95% CI, - 0.64-0.05, p = 0.02) and an effect of education on axial length (0.63 mm per SD increase in the instrument, 95% CI, 0.22-1.04, p = 0.003). Spherical equivalent, axial length and anterior chamber depth were associated with length of education in regression analyses. Mediation analysis revealed that the association between spherical equivalent and education is mainly driven (70%) by alteration in axial length.

CONCLUSIONS

The distribution of axial length and spherical equivalent is represented by subgroups of the population (bi-Gaussian). This distribution can be partially explained by length of education. The impact of education on spherical equivalent is mainly driven by alteration in axial length.

Association of Risk Variants in the CFH Gene With Elevated Levels of Coagulation and Complement Factors in Idiopathic Multifocal Choroiditis

Importance

Idiopathic multifocal choroiditis (MFC) is poorly understood, thereby hindering optimal treatment and monitoring of patients.

Objective

To identify the genes and pathways associated with idiopathic MFC.

Design, Setting, and Participants

This was a case-control genome-wide association study (GWAS) and protein study of blood plasma samples conducted from March 2006 to February 2022. This was a multicenter study involving 6 Dutch universities. Participants were grouped into 2 cohorts: cohort 1 consisted of Dutch patients with idiopathic MFC and controls, and cohort 2 consisted of patients with MFC and controls. Plasma samples from patients with idiopathic MFC who had not received treatment were subjected to targeted proteomics. Idiopathic MFC was diagnosed according to the Standardization of Uveitis Nomenclature (SUN) Working Group guidelines for punctate inner choroidopathy and multifocal choroiditis with panuveitis. Data were analyzed from July 2021 to October 2022.

Main outcomes and measures

Genetic variants associated with idiopathic MFC and risk variants associated with plasma protein concentrations in patients.

Results

This study included a total of 4437 participants in cohort 1 (170 [3.8%] Dutch patients with idiopathic MFC and 4267 [96.2%] controls; mean [SD] age, 55 [18] years; 2443 female [55%]) and 1344 participants in cohort 2 (52 [3.9%] patients with MFC and 1292 [96.1%] controls; 737 male [55%]). The primary GWAS association mapped to the CFH gene with genome-wide significance (lead variant the A allele of rs7535263; odds ratio [OR], 0.52; 95% CI, 0.41-0.64; P = 9.3 × 10-9). There was no genome-wide significant association with classical human leukocyte antigen (HLA) alleles (lead classical allele, HLA-A*31:01; P = .002). The association with rs7535263 showed consistent direction of effect in an independent cohort of 52 cases and 1292 control samples (combined meta-analysis OR, 0.58; 95% CI, 0.38-0.77; P = 3.0 × 10-8). In proteomic analysis of 87 patients, the risk allele G of rs7535263 in the CFH gene was strongly associated with increased plasma concentrations of factor H-related (FHR) proteins (eg, FHR-2, likelihood ratio test, adjusted P = 1.1 × 10-3) and proteins involved in platelet activation and the complement cascade.

Conclusions and relevance

Results suggest that CFH gene variants increase systemic concentrations of key factors of the complement and coagulation cascades, thereby conferring susceptibility to idiopathic MFC. These findings suggest that the complement and coagulation pathways may be key targets for the treatment of idiopathic MFC.

Ocular biological parameters and prevalence of myopia in vocational high school and general high school in China

Significance

Higher prevalence of myopia is possibly associated with more extended schooling schedules. Therefore, adjustments to high school curricula may aid in reducing the prevalence of myopia among adolescents.

Purpose

To investigate the prevalence of myopia among 15- to 18-year-old adolescents in Tianjin, China, and to evaluate the impact of different educational schedules on the prevalence of myopia among high school students.

Methods

This is a school-based epidemiological study with a cross-sectional design. Ocular biological parameters and noncycloplegic photorefraction were examined using optical biometry devices and photoscreener devices. Each student’s spherical equivalent (SE) and ocular biometry were recorded, and the prevalence of myopia was calculated.

Results

A total of 2,867 participants (1,519 males and 1,348 females) were tested for non-cycloplegic refraction, axial length (AL), central corneal thickness (CCT), anterior chamber depth (ACD) and lens thickness (LT). In this research, the overall prevalence of myopia was 81.6%, with high myopia accounting for 11.8%. Myopia prevalence was substantially higher in general high schools than in vocational high schools, with 86.1 and 70.1%, respectively. There were no significant differences in the prevalence of myopia (p = 0.744) or high myopia (p = 0.851) across the three vocational school years. In the general high school, however, there was an increase of 4.6% (p &lt; 0.05) in myopia prevalence between year 10 and year12.

Conclusion

Comparing vocational and standard high school students, there are considerable disparities in prevalence of myopia, spherical equivalent, and ocular biological parameters. The prevalence of myopia and high myopia increased among standard high school students, but remained relatively consistent among students in vocational schools.

Applications of Genomics and Transcriptomics in Precision Medicine for Myopia Control or Prevention

Myopia is a globally emerging concern accompanied by multiple medical and socio-economic burdens with no well-established causal treatment to control thus far. The study of the genomics and transcriptomics of myopia treatment is crucial to delineate disease pathways and provide valuable insights for the design of precise and effective therapeutics. A strong understanding of altered biochemical pathways and underlying pathogenesis leading to myopia may facilitate early diagnosis and treatment of myopia, ultimately leading to the development of more effective preventive and therapeutic measures. In this review, we summarize current data about the genomics and transcriptomics of myopia in human and animal models. We also discuss the potential applicability of these findings to precision medicine for myopia treatment.

Rare variant analyses across multiethnic cohorts identify novel genes for refractive error

Refractive error, measured here as mean spherical equivalent (SER), is a complex eye condition caused by both genetic and environmental factors. Individuals with strong positive or negative values of SER require spectacles or other approaches for vision correction. Common genetic risk factors have been identified by genome-wide association studies (GWAS), but a great part of the refractive error heritability is still missing. Some of this heritability may be explained by rare variants (minor allele frequency [MAF] ≤ 0.01.). We performed multiple gene-based association tests of mean Spherical Equivalent with rare variants in exome array data from the Consortium for Refractive Error and Myopia (CREAM). The dataset consisted of over 27,000 total subjects from five cohorts of Indo-European and Eastern Asian ethnicity. We identified 129 unique genes associated with refractive error, many of which were replicated in multiple cohorts. Our best novel candidates included the retina expressed PDCD6IP, the circadian rhythm gene PER3, and P4HTM, which affects eye morphology. Future work will include functional studies and validation. Identification of genes contributing to refractive error and future understanding of their function may lead to better treatment and prevention of refractive errors, which themselves are important risk factors for various blinding conditions.

Gene-level association analysis of ordinal traits with functional ordinal logistic regressions

In this paper, we develop functional ordinal logistic regression (FOLR) models to perform gene-based analysis of ordinal traits. In the proposed FOLR models, genetic variant data are viewed as stochastic functions of physical positions and the genetic effects are treated as a function of physical positions. The FOLR models are built upon functional data analysis which can be revised to analyze the ordinal traits and high dimension genetic data. The proposed methods are capable of dealing with dense genotype data which is usually encountered in analyzing the next-generation sequencing data. The methods are flexible and can analyze three types of genetic data: (1) rare variants only, (2) common variants only, and (3) a combination of rare and common variants. Simulation studies show that the likelihood ratio test statistics of the FOLR models control type I errors well and have good power performance. The proposed methods achieve the goals of analyzing ordinal traits directly, reducing high dimensionality of dense genetic variants, being computationally manageable, facilitating model convergence, properly controlling type I errors, and maintaining high power levels. The FOLR models are applied to analyze Age-Related Eye Disease Study data, in which two genes are found to strongly associate with four ordinal traits.

The Role of GJD2(Cx36) in Refractive Error Development

Refractive errors are common eye disorders characterized by a mismatch between the focal power of the eye and its axial length. An increased axial length is a common cause of the refractive error myopia (nearsightedness). The substantial increase in myopia prevalence over the last decades has raised public health concerns because myopia can lead to severe ocular complications later in life. Genomewide association studies (GWAS) have made considerable contributions to the understanding of the genetic architecture of refractive errors. Among the hundreds of genetic variants identified, common variants near the gap junction delta-2 (GJD2) gene have consistently been reported as one of the top hits. GJD2 encodes the connexin 36 (Cx36) protein, which forms gap junction channels and is highly expressed in the neural retina. In this review, we provide current evidence that links GJD2(Cx36) to the development of myopia. We summarize the gap junctional communication in the eye and the specific role of GJD2(Cx36) in retinal processing of visual signals. Finally, we discuss the pathways involving dopamine and gap junction phosphorylation and coupling as potential mechanisms that may explain the role of GJD2(Cx36) in refractive error development.

Phenotypic Consequences of the GJD2 Risk Genotype in Myopia Development

Purpose

To study the relatively high effect of the refractive error gene GJD2 in human myopia, and to assess its relationship with refractive error, ocular biometry and lifestyle in various age groups.

Methods

The population-based Rotterdam Study (RS), high myopia case-control study MYopia STudy, and the birth-cohort study Generation R were included in this study. Spherical equivalent (SER), axial length (AL), axial length/corneal radius (AL/CR), vitreous depth (VD), and anterior chamber depth (ACD) were measured using standard ophthalmologic procedures. Biometric measurements were compared between GJD2 (rs524952) genotype groups; education and environmental risk score (ERS) were calculated to estimate gene-environment interaction effects, using the Synergy index (SI).

Results

RS adults carrying two risk alleles had a lower SER and longer AL, ACD and VD (AA versus TT, 0.23D vs. 0.70D; 23.79 mm vs. 23.52 mm; 2.72 mm vs. 2.65 mm; 16.12 mm vs. 15.87 mm; all P < 0.001). Children carrying two risk alleles had larger AL/CR at ages 6 and 9 years (2.88 vs. 2.87 and 3.00 vs. 2.96; all P < 0.001). Education and ERS both negatively influenced myopia and the biometric outcomes, but gene-environment interactions did not reach statistical significance (SI 1.25 [95% confidence interval {CI}, 0.85-1.85] and 1.17 [95% CI, 0.55-2.50] in adults and children).

Conclusions

The elongation of the eye caused by the GJD2 risk genotype follows a dose-response pattern already visible at the age of 6 years. These early effects are an example of how a common myopia gene may drive myopia.

Myopia in African Americans Is Significantly Linked to Chromosome 7p15.2-14.2

Purpose

The purpose of this study was to perform genetic linkage analysis and association analysis on exome genotyping from highly aggregated African American families with nonpathogenic myopia. African Americans are a particularly understudied population with respect to myopia.

Methods

One hundred six African American families from the Philadelphia area with a family history of myopia were genotyped using an Illumina ExomePlus array and merged with previous microsatellite data. Myopia was initially measured in mean spherical equivalent (MSE) and converted to a binary phenotype where individuals were identified as affected, unaffected, or unknown. Parametric linkage analysis was performed on both individual variants (single-nucleotide polymorphisms [SNPs] and microsatellites) as well as gene-based markers. Family-based association analysis and transmission disequilibrium test (TDT) analysis modified for rare variants was also performed.

Results

Genetic linkage analysis identified 2 genomewide significant variants at 7p15.2 and 7p14.2 (in the intergenic region between MIR148A and NFE2L3 and in the noncoding RNA LOC401324) and 2 genomewide significant genes (CRHR2 and AVL9) both at 7p14.3. No genomewide results were found in the association analyses.

Conclusions

This study identified a significant linkage peak in African American families for myopia at 7p15.2 to 7p14.2, the first potential risk locus for myopia in African Americans. Interesting candidate genes are located in the region, including PDE1C, which is highly expressed in the eyes, and known to be involved in retinal development. Further identification of the causal variants at this linkage peak will help elucidate the genetics of myopia in this understudied population.

A Bibliometric and Citation Network Analysis of Myopia Genetics

BACKGROUND

To aim of the study was describe the growth of publications on genetic myopia and understand the current research landscape through the analysis of citation networks, as well as determining the different research areas and the most cited publications.

METHODS

The Web of Science database was used to perform the publication search, looking for the terms "genetic*" AND "myopia" within the period between 2009 and October 2020. The CitNetExplorer and CiteSpace software were then used to conduct the publication analysis. To obtain the graphics, the VOSviewer software was used.

RESULTS

A total of 721 publications were found with 2999 citations generated within the network. The year 2019 was singled out as a "key year", taking into account the number of publications that emerged in that year and given that in 2019, 200 loci associated with refractive errors and myopia were found, which is considered to be great progress. The most widely cited publication was "Genome-wide meta-analyses of multiancestry cohorts identify multiple new susceptibility loci for refractive error and myopia", an article by Verhoeven et al., which was published in 2013. By using the clustering function, we were able to establish three groups that encompassed the different research areas within this field: heritability rate of myopia and its possible association with environmental factors, retinal syndromes associated with myopia and the genetic factors that control and influence axial growth of the eye.

CONCLUSIONS

The citation network offers a comprehensive and objective analysis of the main papers that address genetic myopia.

Genome-Wide Association Study in Asians Identifies Novel Loci for High Myopia and Highlights a Nervous System Role in Its Pathogenesis

PURPOSE

To identify novel susceptibility loci for high myopia.

DESIGN

Genome-wide association study (GWAS) followed by replication and meta-analysis.

PARTICIPANTS

A total of 14 096 samples from East and Southeast Asian populations (2549 patients with high myopia and 11 547 healthy controls).

METHODS

We performed a GWAS in 3269 Japanese individuals (1668 with high myopia and 1601 control participants), followed by replication analysis in a total of 10 827 additional samples (881 with high myopia and 9946 control participants) from Japan, Singapore, and Taiwan. To confirm the biological role of the identified loci in the pathogenesis of high myopia, we performed functional annotation and Gene Ontology (GO) analyses.

MAIN OUTCOME MEASURES

We evaluated the association of single nucleotide polymorphisms with high myopia and GO terms enriched among genes identified in the current study.

RESULTS

We identified 9 loci with genome-wide significance (P < 5.0 × 10^-8). Three loci were previously reported myopia-related loci (ZC3H11B on 1q41, GJD2 on 15q14, and RASGRF1 on 15q25.1), and the other 6 were novel (HIVEP3 on 1p34.2, NFASC/CNTN2 on 1q32.1, CNTN4/CNTN6 on 3p26.3, FRMD4B on 3p14.1, LINC02418 on 12q24.33, and AKAP13 on 15q25.3). The GO analysis revealed a significant role of the nervous system related to synaptic signaling, neuronal development, and Ras/Rho signaling in the pathogenesis of high myopia.

CONCLUSIONS

The current study identified 6 novel loci associated with high myopia and demonstrated an important role of the nervous system in the disease pathogenesis. Our findings give new insight into the genetic factors underlying myopia, including high myopia, by connecting previous findings and allowing for a clarified interpretation of the cause and pathophysiologic features of myopia at the molecular level.

Mutation screening of 17 candidate genes in a cohort of 67 probands with early-onset high myopia

PURPOSE

To detect variants in 17 known potentially causative genes for non-syndromic myopia in 67 Tujia Chinese patients with early-onset high myopia (eo-HM).

METHODS

DNA from 67 unrelated patients with early onset (<7 years old) high myopia (refraction error ≤ -6.00D or axial length > 26 mm) were subjected to whole-exome sequencing (WES). Variants in 17 candidate genes were analysed by multistep bioinformatics analysis. Subsequently, Sanger sequencing was used to verify identified candidate mutations and to assess available family members for co-segregation with myopia.

RESULTS

A multistep systematic analysis of variants in 17 potentially causative genes for eo-HM revealed four novel pathogenic mutations and three potential pathogenic mutations in 4 of 17 genes in 7 of 67 (10.4%) probands. The pathogenic group included one missense mutation (c.100G > C, p.Asp34His) and one splice donor mutation (c.989 + 1G >A) in ARR3, one missense mutation (c.995C > A, p.Thr332Lys) in NDUFAF7 and one novel frameshift mutation (c.726dupA, p.Arg243fs*140) in SLC39A5. The potential pathogenic group included two missense mutations (c.3266A > G, p.Tyr1089Cys; c.913G > A, p.Glu305Lys) in ZNF644 and one missense mutation (c.960T > A, p.His320Gln) in NDUFAF7. Sequence changes were confirmed by Sanger sequencing; all had an allele frequency <0.01 in the 1000G, EVS, ExAC and gnomAD databases. Additionally, both the pathogenic and potentially pathogenic mutations were predicted to be damaging by SIFT, Polyphen-2, PROVEAN, MutationTaster2, CADD and REVEL except the p.Tyr1089Cys and p.Glu305Lys changes were predicted to be neutral by PROVEAN.

CONCLUSION

Our research provides more evidence to support the hypothesis that mutations in ARR3, SLC39A5 and NDUFAF7 are disease-causing genes for eo-HM and broadens the eo-HM mutation spectrum among different ethnic groups. It also deepens understanding of the contributions of ARR3, SLC39A5, and NDUFAF7 to eo-HM.

Myopia genetics in genome-wide association and post-genome-wide association study era

Genome-wide association studies (GWAS) of myopia and refractive error have generated exciting results and identified novel risk-associated loci. However, the interpretation of the findings of GWAS of complex diseases is not straightforward and has remained challenging. This review provides a brief summary of the main focus on the advantages and limitations of GWAS of myopia, with potential strategies that may contribute to further insight into the genetics of myopia in the post-GWAS or omics era.

Frequency and distribution of corneal astigmatism and keratometry features in adult life: Methodology and findings of the UK Biobank study

Purpose

To describe corneal astigmatism in the UK Biobank population and to look for associations with other biometric variables and socio-demographic factors.

Methods

This analysis included a subsample of 107,452 participants of the UK Biobank study who underwent an enhanced ophthalmic examination including autorefractor keratometry (Tomey RC 5000, Tomey Corp., Nagoya, Japan). Participants were recruited from across the United Kingdom between 2006 and 2010, and all were between 40 to 69 years. After quality control and applying relevant exclusions, data on corneal astigmatism on 83,751 participants were included for analysis. Potential associations were tested through univariable regression and significant parameters carried forward for multivariable analysis.

Results

In univariable analysis, the characteristics significantly associated with higher corneal astigmatism (P<0.001), by order of magnitude were, female gender, white ethnicity, lighter skin colour, use of UV protection, lower alcohol intake, lower corneal-compensated intraocular pressure (ccIOP), older age at completion of education, younger age, higher Townsend deprivation index, lower height and lower systolic blood pressure. After inclusion in the multivariable analysis, gender, skin colour, alcohol intake, age at completion of full-time education, ccIOP, age and Townsend deprivation score remained significant (all P<0.001). Increased corneal astigmatism was also found to be significantly associated with amblyopia or strabismus.

Conclusions

This analysis confirms previous associations with astigmatism such as younger age and female gender, and identified novel risk factors including lighter skin colour, lower alcohol intake, later age having completed full time education later, lower ccIOP and higher Townsend deprivation index. Further research is needed to investigate these novel associations.

Genetic variants linked to myopic macular degeneration in persons with high myopia: CREAM Consortium

Purpose

To evaluate the roles of known myopia-associated genetic variants for development of myopic macular degeneration (MMD) in individuals with high myopia (HM), using case-control studies from the Consortium of Refractive Error and Myopia (CREAM).

Methods

A candidate gene approach tested 50 myopia-associated loci for association with HM and MMD, using meta-analyses of case-control studies comprising subjects of European and Asian ancestry aged 30 to 80 years from 10 studies. Fifty loci with the strongest associations with myopia were chosen from a previous published GWAS study. Highly myopic (spherical equivalent [SE] ≤ -5.0 diopters [D]) cases with MMD (N = 348), and two sets of controls were enrolled: (1) the first set included 16,275 emmetropes (SE ≤ -0.5 D); and (2) second set included 898 highly myopic subjects (SE ≤ -5.0 D) without MMD. MMD was classified based on the International photographic classification for pathologic myopia (META-PM).

Results

In the first analysis, comprising highly myopic cases with MMD (N = 348) versus emmetropic controls without MMD (N = 16,275), two SNPs were significantly associated with high myopia in adults with HM and MMD: (1) rs10824518 (P = 6.20E-07) in KCNMA1, which is highly expressed in human retinal and scleral tissues; and (2) rs524952 (P = 2.32E-16) near GJD2. In the second analysis, comprising highly myopic cases with MMD (N = 348) versus highly myopic controls without MMD (N = 898), none of the SNPs studied reached Bonferroni-corrected significance.

Conclusions

Of the 50 myopia-associated loci, we did not find any variant specifically associated with MMD, but the KCNMA1 and GJD2 loci were significantly associated with HM in highly myopic subjects with MMD, compared to emmetropes.

Decreased expression of gap junction delta-2 (GJD2) messenger RNA and connexin 36 protein in form-deprivation myopia of guinea pigs

BACKGROUND

More than ten genome-wide association studies have identified the significant association between the gap junction delta-2 (GJD2) gene and myopia. However, no functional studies have been performed to confirm that this gene is correlated with myopia. This study aimed to observe how this gene changed in mRNA and protein level in the form-deprivation myopia (FDM) animal model.

METHODS

Four-week-old guinea pigs were randomly divided into two groups: control and FDM groups (n = 12 for each group). The right eyes of the FDM group were covered with opaque hemispherical plastic lenses for 3 weeks. For all the animals, refractive status, axial length (AL), and corneal radius of curvature were measured at baseline and 3 weeks later by streak retinoscope, A-scan ultrasonography, and keratometer, respectively. Retinal GJD2 mRNA expression and connexin 36 (Cx36) levels in FDM and control groups were measured by quantitative real-time PCR and Western blot analyses, respectively. Those results were compared using independent t test, Mann-Whitney U test, or paired t test. A significance level of P < 0.05 was used.

RESULTS

Three weeks later, the FDM group (form-deprived eyes) showed about a myopic shift of approximately -6.75 (-7.94 to -6.31) D, while the control group remained hyperopic with only a shift of -0.50 (-0.75 to 0.25) D (Z = -3.38, P < 0.01). The AL increased by 0.74 (0.61-0.76) and 0.10 (0.05-0.21) mm in FDM and control groups, respectively (Z = -3.37, P < 0.01). The relative mRNA expression of GJD2 in the FDM group decreased 31.58% more than the control group (t = 11.44, P < 0.01). The relative protein expression of CX36 on the retina was lowered by 37.72% in form-deprivation eyes as compared to the controls (t = 17.74, P < 0.01).

CONCLUSION

Both the mRNA expression of GJD2 and Cx36 protein amount were significantly decreased in the retina of FDM guinea pigs. This indicates that Cx36 is involved in FDM development, providing compensating evidence for the results obtained from genome-wide association studies.

Genome-wide scans of myopia in Pennsylvania Amish families reveal significant linkage to 12q15, 8q21.3 and 5p15.33

Myopia is one of the most common ocular disorders in the world, yet the genetic etiology of the disease remains poorly understood. Specialized founder populations, such as the Pennsylvania Amish, provide the opportunity to utilize exclusive genomic architecture, like unique haplotypes, to better understand the genetic causes of myopia. We perform genetic linkage analysis on Pennsylvania Amish families that have a strong familial history of myopia to map any potential causal variants and genes for the disease. 293 individuals from 25 extended families were genotyped on the Illumina ExomePlus array and merged with previous microsatellite data. We coded myopia affection as a binary phenotype; myopia was defined as having a mean spherical equivalent (MSE) of less than or equal to - 1 D (diopters). Two-point and multipoint parametric linkage analyses were performed under an autosomal dominant model. When allowing for locus heterogeneity, we identified two novel genome-wide significantly linked variants at 12q15 (heterogeneity LOD, HLOD = 3.77) in PTPRB and at 8q21.3 (HLOD = 3.35) in CNGB3. We identified further three genome-wide significant variants within a single family. These three variants were located in exons of SLC6A18 at 5p15.33 (LODs ranged from 3.51 to 3.37). Multipoint analysis confirmed the significant signal at 5p15.33 with six genome-wide significant variants (LODs ranged from 3.6 to 3.3). Further suggestive evidence of linkage was observed in several other regions of the genome. All three novel linked regions contain strong candidate genes, especially CNGB3 on 8q21.3, which has been shown to affect photoreceptors and cause complete color blindness. Whole genome sequencing on these regions is planned to conclusively elucidate the causal variants.

IMI - Myopia Genetics Report

The knowledge on the genetic background of refractive error and myopia has expanded dramatically in the past few years. This white paper aims to provide a concise summary of current genetic findings and defines the direction where development is needed. We performed an extensive literature search and conducted informal discussions with key stakeholders. Specific topics reviewed included common refractive error, any and high myopia, and myopia related to syndromes. To date, almost 200 genetic loci have been identified for refractive error and myopia, and risk variants mostly carry low risk but are highly prevalent in the general population. Several genes for secondary syndromic myopia overlap with those for common myopia. Polygenic risk scores show overrepresentation of high myopia in the higher deciles of risk. Annotated genes have a wide variety of functions, and all retinal layers appear to be sites of expression. The current genetic findings offer a world of new molecules involved in myopiagenesis. As the missing heritability is still large, further genetic advances are needed. This Committee recommends expanding large-scale, in-depth genetic studies using complementary big data analytics, consideration of gene-environment effects by thorough measurement of environmental exposures, and focus on subgroups with extreme phenotypes and high familial occurrence. Functional characterization of associated variants is simultaneously needed to bridge the knowledge gap between sequence variance and consequence for eye growth.

Association analysis of exome variants and refraction, axial length, and corneal curvature in a European-American population

Refractive errors, myopia, and hyperopia are common visual disorders greatly affecting older individuals. Refraction is determined by genetic factors but only a small percentage of its variation has been explained. We performed a genetic association analysis with three ocular phenotypes: spherical equivalent (a continous measure of refraction), axial length, and corneal curvature in 1,871 European-Americans from the Beaver Dam Eye Study. Individuals were genotyped on the Illumina exome array and imputed to the Haplotype Reference Consortium reference panel. After increasing the number of analyzed variants in targeted protein-coding regions 10-fold via imputation, we confirmed associations for two previously known loci with corneal curvature (chr4q12, rs2114039; g.55092626T > C, β = -0.03 (95% confidence interval [CI]): -0.06, -0.01, P value = 0.01) and spherical equivalent (chr15q14, rs634990; g.35006073T > C, β = -0.27, 95% CI: -0.45, -0.09, P value = 3.79 × 10^-3 ). Despite increased single nucleotide polymorphism (SNP) density, we did not detect any novel significant variants after correction for multiple comparisons. In summary, we confirmed two previous loci associated with corneal curvature and spherical equivalent in a European-American population highlighting the potential biological role of those regions in these traits.

Vitamin D and its pathway genes in myopia: systematic review and meta-analysis

OBJECTIVE

To conduct a systematic review and meta-analysis of the association of blood vitamin D (25-hydroxyvitamin D, 25(OH)D) concentration and vitamin D pathway genes with myopia.

METHODS

We searched the MEDLINE and EMBASE databases for studies published up to 29 January 2018. Cross-sectional or cohort studies which evaluated the blood 25(OH)D concentration, blood 25(OH)D3 concentration or vitamin D pathway genes, in relation to risk of myopia or refractive errors were included. Standard mean difference (SMD) of blood 25(OH)D concentrations between the myopia and non-myopia groups was calculated. The associations of blood 25(OH)D concentrations and polymorphisms in vitamin D pathway genes with myopia using summary ORs were evaluated.

RESULTS

We summarised seven studies involving 25 008 individuals in the meta-analysis. The myopia group had lower 25(OH)D concentration than the non-myopia group (SMD=-0.27 nmol/L, p=0.001). In the full analysis, the risk of myopia was inversely associated with blood 25(OH)D concentration after adjusting for sunlight exposure or time spent outdoors (adjusted odds ratio (AOR)=0.92 per 10 nmol/L, p<0.0001). However, the association was not statistically significant for the <18 years subgroup (AOR=0.91 per 10 nmol/L, p=0.13) and was significant only for 25(OH)D3 (likely to be mainly sunlight derived), but not total 25(OH)D (AOR=0.93 per 10 nmol/L, p=0.00007; AOR=0.91 per 10 nmol/L, p=0.15). We analysed four single nucleotide polymorphisms in the VDR gene from two studies; there was no significant association with myopia.

CONCLUSIONS

Lower 25(OH)D is associated with increased risk of myopia; the lack of a genetic association suggests that 25(OH)D level may be acting as a proxy for time outdoors.

Genome-wide association meta-analysis highlights light-induced signaling as a driver for refractive error

Refractive errors, including myopia, are the most frequent eye disorders worldwide and an increasingly common cause of blindness. This genome-wide association meta-analysis in 160,420 participants and replication in 95,505 participants increased the number of established independent signals from 37 to 161 and showed high genetic correlation between Europeans and Asians (>0.78). Expression experiments and comprehensive in silico analyses identified retinal cell physiology and light processing as prominent mechanisms, and also identified functional contributions to refractive-error development in all cell types of the neurosensory retina, retinal pigment epithelium, vascular endothelium and extracellular matrix. Newly identified genes implicate novel mechanisms such as rod-and-cone bipolar synaptic neurotransmission, anterior-segment morphology and angiogenesis. Thirty-one loci resided in or near regions transcribing small RNAs, thus suggesting a role for post-transcriptional regulation. Our results support the notion that refractive errors are caused by a light-dependent retina-to-sclera signaling cascade and delineate potential pathobiological molecular drivers.

Genetic Association Study Between the COL11A1 and COL18A1 Genes and High Myopia in a Han Chinese Population

PURPOSE

To evaluate the association between high myopia (HM) and single nucleotide polymorphisms (SNPs) in collagen, type XI, alpha 1 (COL11A1) and collagen, type XVIII, alpha 1 (COL18A1) genes in a Han Chinese population.

MATERIALS AND METHODS

A total of 869 patients with HM and 804 controls were recruited for this study. The genotyping of five SNPs in COL11A1 and COL18A1 was performed using the SNaPshot method. The genotyping data were analyzed using the χ² test, and the linkage disequilibrium block structure was calculated and examined by Haploview software.

RESULTS

No statistically significant differences (p > 0.05) were identified between HM cases and controls after a Bonferroni correction for multiple tests in the allele frequencies of COL11A1 and COL18A1 SNPs. However, the G allele of rs2236475 showed a susceptible effect for HM (p = 0.016, corrected p = 0.08, odds ratio [OR] = 1.26). Moreover, the carriers of rs2236475GG genotype displayed an increased risk of HM compared with the rs2236475AA and rs2236475AG+AA genotypes (p = 0.008, OR = 1.79, confidence interval [95% CI] = 1.18-2.64, uncorrected; p = 0.012, OR = 1.74, 95% CI = 1.12-2.57, corrected, respectively).

CONCLUSIONS

Our results suggested that common polymorphisms in these two candidate genes were unlikely to play major roles in the genetic susceptibility to HM. Nevertheless, to avoid filtering real myopia genes, the role of COL11A1 and COL18A1 in the pathogenesis of myopia requires more refinement in both animal models and human genetic epidemiological studies.

Myopia in Chinese families shows linkage to 10q26.13

Purpose

To determine genetic linkage between myopia and Han Chinese patients with a family history of the disease.

Methods

One hundred seventy-six Han Chinese patients from 34 extended families were given eye examinations, and mean spherical equivalent (MSE) in diopters (D) was calculated by adding the spherical component of the refraction to one-half the cylindrical component and taking the average of both eyes. The MSE was converted to a binary phenotype, where all patients with an MSE of -1.00 D or less were coded as affected. Unaffected individuals had an MSE greater than 0.00 D (ages 21 years and up), +1.50 (ages 11-20), or +2.00 D (ages 6-10 years). Individuals between the given upper threshold and -1.00 were coded as unknown. Patients were genotyped on an exome chip. Three types of linkage analyses were performed: single-variant two-point, multipoint, and collapsed haplotype pattern (CHP) variant two-point.

Results

The CHP variant two-point results identified a significant peak (heterogeneity logarithm of the odds [HLOD] = 3.73) at 10q26.13 in TACC2. The single-variant two-point and multipoint analyses showed highly suggestive linkage to the same region. The single-variant two-point results identified 25 suggestive variants at HTRA1, also at 10q26.13.

Conclusions

We report a significant genetic linkage between myopia and Han Chinese patients at 10q26.13. 10q26.13 contains several good candidate genes, such as TACC2 and the known age-related macular degeneration gene HTRA1. Targeted sequencing of the region is planned to identify the causal variant(s).

Insight into the molecular genetics of myopia

Myopia is the most common cause of visual impairment worldwide. Genetic and environmental factors contribute to the development of myopia. Studies on the molecular genetics of myopia are well established and have implicated the important role of genetic factors. With linkage analysis, association studies, sequencing analysis, and experimental myopia studies, many of the loci and genes associated with myopia have been identified. Thus far, there has been no systemic review of the loci and genes related to non-syndromic and syndromic myopia based on the different approaches. Such a systemic review of the molecular genetics of myopia will provide clues to identify additional plausible genes for myopia and help us to understand the molecular mechanisms underlying myopia. This paper reviews recent genetic studies on myopia, summarizes all possible reported genes and loci related to myopia, and suggests implications for future studies on the molecular genetics of myopia.

Development of Refractive Errors-What Can We Learn From Inherited Retinal Dystrophies?

PURPOSE

It is unknown which retinal cells are involved in the retina-to-sclera signaling cascade causing myopia. As inherited retinal dystrophies (IRD) are characterized by dysfunction of a single retinal cell type and have a high risk of refractive errors, a study investigating the affected cell type, causal gene, and refractive error in IRDs may provide insight herein.

DESIGN

Case-control study.

METHODS

Study Population: Total of 302 patients with IRD from 2 ophthalmogenetic centers in the Netherlands. Reference Population: Population-based Rotterdam Study-III and Erasmus Rucphen Family Study (N = 5550). Distributions and mean spherical equivalent (SE) were calculated for main affected cell type and causal gene; and risks of myopia and hyperopia were evaluated using logistic regression.

RESULTS

Bipolar cell-related dystrophies were associated with the highest risk of SE high myopia 239.7; odds ratio (OR) mild hyperopia 263.2, both P < .0001; SE -6.86 diopters (D) (standard deviation [SD] 6.38), followed by cone-dominated dystrophies (OR high myopia 19.5, P < .0001; OR high hyperopia 10.7, P = .033; SE -3.10 D [SD 4.49]); rod dominated dystrophies (OR high myopia 10.1, P < .0001; OR high hyperopia 9.7, P = .001; SE -2.27 D [SD 4.65]), and retinal pigment epithelium (RPE)-related dystrophies (OR low myopia 2.7; P = .001; OR high hyperopia 5.8; P = .025; SE -0.10 D [SD 3.09]). Mutations in RPGR (SE -7.63 D [SD 3.31]) and CACNA1F (SE -5.33 D [SD 3.10]) coincided with the highest degree of myopia and in CABP4 (SE 4.81 D [SD 0.35]) with the highest degree of hyperopia.

CONCLUSIONS

Refractive errors, in particular myopia, are common in IRD. The bipolar synapse and the inner and outer segments of the photoreceptor may serve as critical sites for myopia development.

Caucasian Families Exhibit Significant Linkage of Myopia to Chromosome 11p

Purpose

Myopia is a common visual disorder caused by eye overgrowth, resulting in blurry vision. It affects one in four Americans, and its prevalence is increasing. The genetic mechanisms that underpin myopia are not completely understood. Here, we use genotype data and linkage analyses to identify high-risk genetic loci that are significantly linked to myopia.

Methods

Individuals from 56 Caucasian families with a history of myopia were genotyped on an exome-based array, and the single nucleotide polymorphism (SNP) data were merged with microsatellite genotype data. Refractive error measures on the samples were converted into binary phenotypes consisting of affected, unaffected, or unknown myopia status. Parametric linkage analyses assuming an autosomal dominant model with 90% penetrance and 10% phenocopy rate were performed.

Results

Single variant two-point analyses yielded three significantly linked SNPs at 11p14.1 and 11p11.2; a further 45 SNPs at 11p were found to be suggestive. No other chromosome had any significant SNPs or more than seven suggestive linkages. Two of the significant SNPs were located in BBOX1-AS1 and one in the intergenic region between ORA47 and TRIM49B. Collapsed haplotype pattern two-point analysis and multipoint analyses also yielded multiple suggestively linked genes at 11p. Multipoint analysis also identified suggestive evidence of linkage on 20q13.

Conclusions

We identified three genome-wide significant linked variants on 11p for myopia in Caucasians. Although the novel specific signals still need to be replicated, 11p is a promising region that has been identified by other linkage studies with a number of potentially interesting candidate genes. We hope that the identification of these regions on 11p as potential causal regions for myopia will lead to more focus on these regions and maybe possible replication of our specific linkage peaks in other studies. We further plan targeted sequencing on 11p for our most highly linked families to more clearly understand the source of the linkage in this region.

Epistatic Gene-Based Interaction Analyses for Glaucoma in eMERGE and NEIGHBOR Consortium

Primary open angle glaucoma (POAG) is a complex disease and is one of the major leading causes of blindness worldwide. Genome-wide association studies have successfully identified several common variants associated with glaucoma; however, most of these variants only explain a small proportion of the genetic risk. Apart from the standard approach to identify main effects of variants across the genome, it is believed that gene-gene interactions can help elucidate part of the missing heritability by allowing for the test of interactions between genetic variants to mimic the complex nature of biology. To explain the etiology of glaucoma, we first performed a genome-wide association study (GWAS) on glaucoma case-control samples obtained from electronic medical records (EMR) to establish the utility of EMR data in detecting non-spurious and relevant associations; this analysis was aimed at confirming already known associations with glaucoma and validating the EMR derived glaucoma phenotype. Our findings from GWAS suggest consistent evidence of several known associations in POAG. We then performed an interaction analysis for variants found to be marginally associated with glaucoma (SNPs with main effect p-value <0.01) and observed interesting findings in the electronic MEdical Records and GEnomics Network (eMERGE) network dataset. Genes from the top epistatic interactions from eMERGE data (Likelihood Ratio Test i.e. LRT p-value <1e-05) were then tested for replication in the NEIGHBOR consortium dataset. To replicate our findings, we performed a gene-based SNP-SNP interaction analysis in NEIGHBOR and observed significant gene-gene interactions (p-value <0.001) among the top 17 gene-gene models identified in the discovery phase. Variants from gene-gene interaction analysis that we found to be associated with POAG explain 3.5% of additional genetic variance in eMERGE dataset above what is explained by the SNPs in genes that are replicated from previous GWAS studies (which was only 2.1% variance explained in eMERGE dataset); in the NEIGHBOR dataset, adding replicated SNPs from gene-gene interaction analysis explain 3.4% of total variance whereas GWAS SNPs alone explain only 2.8% of variance. Exploring gene-gene interactions may provide additional insights into many complex traits when explored in properly designed and powered association studies.

The complex nature of primary-open angle glaucoma (POAG) has left researchers exploring the genetic architecture and searching for the missing heritability using a number of different study designs. Over the past decade, many studies have been conducted to explain the etiology of POAG; however, a high proportion of estimated heritability still remains unexplained. GWA studies for POAG have identified significant associations but these associations have only explained a small proportion of the genetic risk (odds ratios range between 1–3). In this paper, we sought to confirm the primary genome-wide significant associations that have been discovered so far for glaucoma in phenotypes developed from EMR data in an effort to show that EMR data can be a powerful resource for finding genetic variants influencing POAG susceptibility. Next, we tested for statistical interactions, which can be presented as an important tool in an attempt to explain POAG heritability. We used a reduced list of variants filtered by marginal main effect analysis to look for epistatic interactions. We present our results from replication of gene-based interaction analyses performed in eMERGE and the NEIGHBOR consortium data. Using expression data and annotations from various publicly available databases, the most significant genes that replicated in our analyses show expression in the eye and trabecular meshwork. Analysis for estimation of genetic variance explained by significant associations from previous GWAS and replicated variants from gene-based interactions suggest that these explain 5.6% of variance in eMERGE dataset and also explain 3.4% variance in NEIGHBOR dataset.

Variation in PTCHD2, CRISP3, NAP1L4, FSCB, and AP3B2 associated with spherical equivalent

PURPOSE

Ocular refraction is measured in spherical equivalent as the power of the external lens required to focus images on the retina. Myopia (nearsightedness) and hyperopia (farsightedness) are the most common refractive errors, and the leading causes of visual impairment and blindness in the world. The goal of this study is to identify rare and low-frequency variants that influence spherical equivalent.

METHODS

We conducted variant-level and gene-level quantitative trait association analyses for mean spherical equivalent, using data from 1,560 individuals in the Beaver Dam Eye Study. Genotyping was conducted using the Illumina exome array. We analyzed 34,976 single nucleotide variants and 11,571 autosomal genes across the genome, using single-variant tests as well as gene-based tests.

RESULTS

Spherical equivalent was significantly associated with five genes in gene-based analysis: PTCHD2 at 1p36.22 (p = 3.6 × 10(-7)), CRISP3 at 6p12.3 (p = 4.3 × 10(-6)), NAP1L4 at 11p15.5 (p = 3.6 × 10(-6)), FSCB at 14q21.2 (p = 1.5 × 10(-7)), and AP3B2 at 15q25.2 (p = 1.6 × 10(-7)). The variant-based tests identified evidence suggestive of association with two novel variants in linkage disequilibrium (pairwise r(2) = 0.80) in the TCTE1 gene region at 6p21.1 (rs2297336, minor allele frequency (MAF) = 14.1%, β = -0.62 p = 3.7 × 10(-6); rs324146, MAF = 16.9%, β = -0.55, p = 1.4 × 10(-5)). In addition to these novel findings, we successfully replicated a previously reported association with rs634990 near GJD2 at 15q14 (MAF = 47%, β = -0.29, p=1.8 × 10(-3)). We also found evidence of association with spherical equivalent on 2q37.1 in PRSS56 at rs1550094 (MAF = 31%, β = -0.33, p = 1.7 × 10(-3)), a region previously associated with myopia.

CONCLUSIONS

We identified several novel candidate genes that may play a role in the control of spherical equivalent. However, further studies are needed to replicate these findings. In addition, our results contribute to the increasing evidence that variation in the GJD2 and PRSS56 genes influence the development of refractive errors. Identifying that variation in these genes is associated with spherical equivalent may provide further insight into the etiology of myopia and consequent vision loss.

Polymorphism in the RASGRF1 gene with high myopia: A meta-analysis

PURPOSE

To evaluate the association between polymorphisms of RASGRF1 rs8027411 and high myopia in Chinese and Japanese populations.

METHODS

All eligible studies investigating the association between the RASGRF1 gene and high myopia listed in PubMed, EMBASE, the Cochrane Library, Web of Science, and the China Biologic Medical Database were retrieved. The effects were assessed with the pooled odds ratio (OR) and 95% confidence interval (CI). Heterogeneity between studies was evaluated with the Q-statistic test. Publication bias was tested with Begg's and Egger's linear regression tests. Subgroup analysis and sensitivity analysis were performed to identify the sources of heterogeneity.

RESULTS

In the present meta-analysis, 2,529 individuals with high myopia and 3,127 controls from four studies were included and divided into seven groups. The results indicated that RASGRF1 rs8027411 was significantly associated with high myopia in Chinese and Japanese populations. Carriers of the rs8027411 G allele had a lower risk of high myopia compared to carriers with the T allele (G versus T, OR=0.83, 95% CI=0.77-0.89; p<0.001). Low but not significant heterogeneity was found in a recessive model. No heterogeneity was found in other genetic models. The subgroup analysis indicated that the protective effect of rs8027411 variants was more prominent in Chinese populations (G versus T, OR=0.80 in Chinese and OR=0.86 in Japanese; GG versus TT, OR=0.65 in Chinese and OR=0.77 in Japanese; GT versus TT, OR=0.76 in Chinese and OR=0.81 in Japanese; (GG+GT) versus TT, OR=0.73 in Chinese and OR=0.80 in Japanese; and GG versus (GT+TT), OR=0.77 in Chinese and OR=0.87 in Japanese). Sensitivity analysis indicated that the study results were stable in allelic, homozygote, heterozygote, and dominant models but were not stable in the recessive model. No evidence of publication bias was found.

CONCLUSIONS

Carriers of the rs8027411 G allele in the RASGRF1 gene may be at a lower risk of high myopia in Chinese and Japanese populations. The RASGRF1 gene may play a role in the development of high myopia, especially in Asians. Additional studies are required to validate these results.

Distribution of astigmatism as a function of age in an Australian population

PURPOSE

Astigmatism is a common cause of refractive error and is known to vary in prevalence with age. Although the search for genes associated with spherical refractive errors (especially myopia) has met with limited success, current efforts to identify genetic variants implicated in astigmatism development have been less rewarding. We aimed to assess the association between astigmatism and age to identify appropriate age cut-offs for maximizing power in genetic studies of astigmatism.

METHODS

We performed a cross-sectional analysis of right eye astigmatism data from four Australian-based eye studies comprising 3841 participants aged 5-90 years. Measurements were performed under cycloplegia using an autorefractor, and individuals with a history of cataract, refractive surgery or corneal pathology were excluded from the analysis. In addition to the magnitude and type (against-the-rule, with-the-rule, and oblique) of astigmatism, we calculated the vector components (J0 , J45 ) and evaluated the association of these outcome measures with age.

RESULTS

The magnitude of refractive astigmatism (RA) remained relatively stable [mean ± SD (-0.44 D ± 0.50)] until individuals reached the age of 50, thereafter increasing in average magnitude by approximately 1.00 D for those subjects aged 90. In contrast, corneal astigmatism (CA) remained relatively stable from childhood until the age of 80 (-0.76 D ± 0.61). The prevalence of clinically significant RA (≥1.00 D) increased with age and was highest in those aged >70 years [55.1% (47.2-62.7%)]. Age was significantly associated with RA in adults [odds ratio (OR) = 1.04 per 1 year, p < 0.001]. A weaker relationship was observed between CA and age (OR = 1.007 per 1 year, p = 0.02).

CONCLUSIONS

We have confirmed the previously documented association between RA and age. Our results indicate that most of the observed change occurs after the age of 50, providing a recommended cut-off for participants in genetic studies of this refractive condition.

The progression of myopia from its onset at age 8-12 to adulthood and the influence of heredity and external factors on myopic progression. A 23-year follow-up study

PURPOSE

To examine myopic progression and factors connected with myopic progression.

METHODS

Myopic schoolchildren, with no previous spectacles, 119 boys and 121 girls, were recruited during 1983-1984 to a randomized 3-year clinical trial of bifocal treatment of myopia with a subsequent 20-year follow-up. Participants' mean age at Baseline was 10.9, ranging from 8.7 to 12.8 years. An ophthalmological examination was carried out annually for 3 years and twice thereafter at ca. 10-year intervals. Additional refraction values were received from prescriptions issued by different ophthalmologists and opticians. Altogether, 1915 refraction values were available. Reading distance and accommodation were measured at each control visit. Data on parents' myopia, daily time spent on reading and close work, outdoor activities and watching television were gathered with a structured questionnaire.

RESULTS

Using bifocals (+1.75 add) or reading without glasses or accommodation stimulus during the 3-year period in childhood did not correlate with adulthood refraction. Short reading distance in childhood predicted higher adulthood myopia among females. The factors predicting faster myopic progression were parents' myopia and less time spent on sports and outdoor activities at childhood. Time spent on reading and close work in childhood was related to myopic progression during the first 3 years but did not predict adulthood myopia. Myopia throughout follow-up was higher among those who watched television <3 hr daily than those who spent more time watching television. Mean myopic progression 8 years after age 20-24 was -0.45 D ± 0.71 (SD), and in 45% of cases, progression was ≥0.5 D.

CONCLUSIONS

In nearly half of the cases, myopia beginning at school continued to progress into adulthood. Higher adulthood myopia was mainly related to parents' myopia and less time spent on sports and outdoor activities in childhood.

Genome-wide meta-analysis of myopia and hyperopia provides evidence for replication of 11 loci

Refractive error (RE) is a complex, multifactorial disorder characterized by a mismatch between the optical power of the eye and its axial length that causes object images to be focused off the retina. The two major subtypes of RE are myopia (nearsightedness) and hyperopia (farsightedness), which represent opposite ends of the distribution of the quantitative measure of spherical refraction. We performed a fixed effects meta-analysis of genome-wide association results of myopia and hyperopia from 9 studies of European-derived populations: AREDS, KORA, FES, OGP-Talana, MESA, RSI, RSII, RSIII and ERF. One genome-wide significant region was observed for myopia, corresponding to a previously identified myopia locus on 8q12 (p = 1.25×10⁻⁸), which has been reported by Kiefer et al. as significantly associated with myopia age at onset and Verhoeven et al. as significantly associated to mean spherical-equivalent (MSE) refractive error. We observed two genome-wide significant associations with hyperopia. These regions overlapped with loci on 15q14 (minimum p value = 9.11×10⁻¹¹) and 8q12 (minimum p value 1.82×10⁻¹¹) previously reported for MSE and myopia age at onset. We also used an intermarker linkage- disequilibrium-based method for calculating the effective number of tests in targeted regional replication analyses. We analyzed myopia (which represents the closest phenotype in our data to the one used by Kiefer et al.) and showed replication of 10 additional loci associated with myopia previously reported by Kiefer et al. This is the first replication of these loci using myopia as the trait under analysis. “Replication-level” association was also seen between hyperopia and 12 of Kiefer et al.'s published loci. For the loci that show evidence of association to both myopia and hyperopia, the estimated effect of the risk alleles were in opposite directions for the two traits. This suggests that these loci are important contributors to variation of refractive error across the distribution.

Severe anisometropic myopia in identical twins

High anisometropic myopia is a rare condition in twins. Genetic factors have been implicated in its development and there may be an association with vision-threatening complications. A pair of 11-year-old twins presented with poor distance vision in both eyes. Detailed ocular examination was performed including slit lamp examination, dilated funduscopy, cycloplegic refraction, keratometry and axial length measurement. The objective refraction was -6.50DS -2.00DC × 180 (right eye), -1.00DS (left eye) for the first twin; -13.75DS -2.25DC × 180 (right eye), -0.50DS -0.75DC × 04 (left eye) for the second twin. This case suggests an underlying genetic defect in the development of myopia.

Genome-wide association studies: applications and insights gained in Ophthalmology

Genome-wide association studies (GWAS) use high-throughput genotyping technologies to genotype thousands of single-nucleotide polymorphisms (SNPs) and relate them to the development of clinical and quantitative traits. Their use has been highly successful in the field of ophthalmology, and since the advent of GWAS in 2005, many genes not previously suspected of having a role in disease have been identified and the findings replicated. We conducted an extensive literature review and describe the concept, design, advantages, and limitations of GWAS and provide a detailed description of the applications and discoveries of GWAS in the field of eye disease to date. There have been many novel findings revealing previously unknown biological insights in a diverse range of common ocular conditions. GWAS have been a highly successful modality for investigating the pathogenesis of a wide variety of ophthalmic conditions. The insights gained into the pathogenesis of disease provide not only a better understanding of underlying disease mechanism but also offer a rationale for targeted treatment and preventative strategies. Expansive international collaboration and standardised phenotyping will permit the continued success of this investigative technique.

Meta-analysis of the association between lumican gene polymorphisms and susceptibility to high Myopia

Backgrounds

Many studies have evaluated the association between lumican (LUM) gene polymorphisms and high myopia. However, the results remain controversial. This meta-analysis aims to comprehensively evaluate the relationship between two common LUM polymorphisms (rs3759223 and rs3759222) and the risk of high myopia.

Methods

A comprehensive literature search for studies published up until September of 2013 was performed. Data were extracted independently by two investigators, and the weighted Odds Ratios (ORs) and 95% Confidence Intervals (CIs) for the associations were obtained by using a random-effects model.

Results

Eight studies (1425cases and 1271 controls) were identified for the analysis of the association between rs3759223 polymorphism and high myopia. The results indicated that rs3759223 polymorphism was associated with high myopia under a recessive model (OR = 1.71, 95%CI 1.04–2.81). Further subgroup analysis indicated that this polymorphism was associated with high myopia among Chinese people in the additive model (OR = 1.17, 95%CI 1.06–1.29) and a recessive model (OR = 1.75, 95%CI 1.00–3.06) with control group coming from hospital based population. Four studies (1024 cases and 1163 controls) were identified for the analysis of the association between rs3759222 polymorphism and high myopia. The results indicated that rs3759222 polymorphism was not associated with high myopia in all genetic models, even the subgroup analysis couldn't provide relative proof to assure the outcome.

Conclusion

This meta-analysis suggests that LUM polymorphisms are associated with the risk of high myopia. However, well-designed studies with larger sample sizes and more ethnic groups are required to further validate this association.

Myopia and level of education: results from the Gutenberg Health Study

PURPOSE

To analyze the association between myopia and educational level in an adult European cohort.

DESIGN

Population-based cross-sectional study.

PARTICIPANTS

A cohort of the Gutenberg Health Study, including 4658 eligible enrollees between 35 and 74 years of age.

METHODS

We applied a standardized protocol entailing a comprehensive questionnaire; thorough ophthalmic, general, cardiovascular, and psychological examinations; and laboratory tests, including genetic analyses. We documented achievement levels in school education and post-school professional education. The spherical equivalent (SE) was determined by noncycloplegic autorefractometry. We fitted mixed linear models including age, gender, and 45 myopia-associated single nucleotide polymorphisms (SNP) as covariates.

MAIN OUTCOME MEASURES

Prevalence and magnitude of myopia in association with years spent in school and level of post-school professional education.

RESULTS

Individuals who graduated from school after 13 years were more myopic (median, -0.5 diopters [D]; first quartile [Q1]/third quartile [Q3], -2.1/0.3 D) than those who graduated after 10 years (median, -0.2 D; Q1/Q3, -1.3/0.8 D), than those who graduated after 9 years (median, 0.3 D; Q1/Q3, -0.6/1.4 D), and than those who never finished secondary school (median, 0.2 D; Q1/Q3, -0.5/1.8 D; P<0.001, respectively). The same holds true for persons with a university degree (median, -0.6 D; Q1/Q3, -2.3/0.3 D) versus those who finished secondary vocational school (median, 0 D; Q1/Q3, -1.1/0.8 D) or primary vocational school (median, 0 D; Q1/Q3, -0.9/1.1 D) versus persons without any post-school professional qualification (median, 0.6 D; Q1/Q3, -0.4/1.7 D; P<0.001, respectively). Of persons who graduated from school after 13 years, 50.9% were myopic (SE, ≤-0.5 D) versus 41.6%, 27.1%, and 26.9% after 10 years, in those who graduated after 9 years, and in those who never graduated from secondary school, respectively (P<0.001). In university graduates, the proportion of myopic persons was higher (53%) than that of those who graduated from secondary (34.8%) or primary (34.7%) vocational schools and than in those without any professional training (23.9%; P<0.001, respectively). In multivariate analyses: higher school and professional levels of education were associated with a more myopic SE independent of gender. There was a small effect of age and SNPs.

CONCLUSIONS

Higher levels of school and post-school professional education are associated with a more myopic refraction. Participants with higher educational achievements more often were myopic than individuals with less education.

Genome-wide association studies of refractive error and myopia, lessons learned, and implications for the future

The investigation of the genetic basis of refractive error and myopia entered a new stage with the introduction of genome-wide association studies (GWAS). Multiple GWAS on many ethnic groups have been published over the years, providing new insight into the genetic architecture and pathophysiology of refractive error. This is a review of the GWAS published to date, the main lessons learned, and future possible directions of genetic studies of myopia and refractive error.

Education influences the association between genetic variants and refractive error: a meta-analysis of five Singapore studies

Refractive error is a complex ocular trait governed by both genetic and environmental factors and possibly their interplay. Thus far, data on the interaction between genetic variants and environmental risk factors for refractive errors are largely lacking. By using findings from recent genome-wide association studies, we investigated whether the main environmental factor, education, modifies the effect of 40 single nucleotide polymorphisms on refractive error among 8461 adults from five studies including ethnic Chinese, Malay and Indian residents of Singapore. Three genetic loci SHISA6-DNAH9, GJD2 and ZMAT4-SFRP1 exhibited a strong association with myopic refractive error in individuals with higher secondary or university education (SHISA6-DNAH9: rs2969180 A allele, β = -0.33 D, P = 3.6 × 10(-6); GJD2: rs524952 A allele, β = -0.31 D, P = 1.68 × 10(-5); ZMAT4-SFRP1: rs2137277 A allele, β = -0.47 D, P = 1.68 × 10(-4)), whereas the association at these loci was non-significant or of borderline significance in those with lower secondary education or below (P for interaction: 3.82 × 10(-3)-4.78 × 10(-4)). The evidence for interaction was strengthened when combining the genetic effects of these three loci (P for interaction = 4.40 × 10(-8)), and significant interactions with education were also observed for axial length and myopia. Our study shows that low level of education may attenuate the effect of risk alleles on myopia. These findings further underline the role of gene-environment interactions in the pathophysiology of myopia.

[Ophthalmological aspects of the Gutenberg Health Study (GHS): an interdisciplinary prospective population-based cohort study]

The Gutenberg Health Study (GHS) is a prospective, single center, population-based cohort study. It is an interdisciplinary project for investigation of ocular, cardiovascular, psychosomatic and immune diseases in the population of the City of Mainz and the district of Mainz-Bingen in central Germany. The main goals of the ophthalmological branch of GHS are determination of the prevalence and incidence of major ocular risk factors and diseases, exploring the genetic determinants and assessing complex interdisciplinary associations. The study cohort includes 15,010 participants aged between 35 and 74 years at the time of inclusion. After completing the baseline investigations the 5-year follow-up of the whole study cohort started in April 2012. The GHS is the most extensive data set of major ophthalmological conditions and risk factors in Germany and will help researchers in understanding complex medical associations.

Association study of fibroblast growth factor 10 (FGF10) polymorphisms with susceptibility to extreme myopia in a Japanese population

PURPOSE

The fibroblast growth factor 10 (FGF10) gene polymorphism rs339501 was previously reported to be associated with high myopia in a Chinese population. In the present study, we investigated whether FGF10 polymorphisms are associated with extreme myopia in a Japanese population as well.

METHODS

A total of 433 Japanese patients with extreme myopia (≤ -10.00 diopters) and 542 Japanese healthy controls (+1.50 to -1.50 diopters) were recruited. We genotyped seven tagging single-nucleotide polymorphisms (SNPs), including rs339501, in FGF10. We also performed an imputation analysis to evaluate the potential association of ungenotyped FGF10 SNPs, and 34 SNPs were imputed.

RESULTS

It was found that rs339501 and rs12517396 exhibited the strongest association with extreme myopia (p=3.9 × 10⁻⁴, corrected p [Pc]=0.0030). A significant association was also observed for rs10462070 (p=6.5 × 10⁻⁴, Pc=0.0059). These three SNPs were in strong linkage disequilibrium (D' ≥0.99, r² ≥0.96). However, the frequency of the A allele of rs339501 was increased in cases compared to controls, which differs from the increased frequency of the G allele in cases in the previous Chinese population.

CONCLUSIONS

Three FGF10 SNPs in complete linkage disequilibrium--rs339501, rs12517396, and rs10462070--were associated with extreme myopia in the Japanese population, and the risk allele of rs339501 differed from the previous Chinese population. Therefore, these three SNPs may not be an important risk factor for susceptibility to extreme myopia. Further studies are needed to elucidate the possible contribution of the FGF10 region in the development of extreme myopia.

Regional replication of association with refractive error on 15q14 and 15q25 in the Age-Related Eye Disease Study cohort

PURPOSE

Refractive error is a complex trait with multiple genetic and environmental risk factors, and is the most common cause of preventable blindness worldwide. The common nature of the trait suggests the presence of many genetic factors that individually may have modest effects. To achieve an adequate sample size to detect these common variants, large, international collaborations have formed. These consortia typically use meta-analysis to combine multiple studies from many different populations. This approach is robust to differences between populations; however, it does not compensate for the different haplotypes in each genetic background evidenced by different alleles in linkage disequilibrium with the causative variant. We used the Age-Related Eye Disease Study (AREDS) cohort to replicate published significant associations at two loci on chromosome 15 from two genome-wide association studies (GWASs). The single nucleotide polymorphisms (SNPs) that exhibited association on chromosome 15 in the original studies did not show evidence of association with refractive error in the AREDS cohort. This paper seeks to determine whether the non-replication in this AREDS sample may be due to the limited number of SNPs chosen for replication.

METHODS

We selected all SNPs genotyped on the Illumina Omni2.5v1_B array or custom TaqMan assays or imputed from the GWAS data, in the region surrounding the SNPs from the Consortium for Refractive Error and Myopia study. We analyzed the SNPs for association with refractive error using standard regression methods in PLINK. The effective number of tests was calculated using the Genetic Type I Error Calculator.

RESULTS

Although use of the same SNPs used in the Consortium for Refractive Error and Myopia study did not show any evidence of association with refractive error in this AREDS sample, other SNPs within the candidate regions demonstrated an association with refractive error. Significant evidence of association was found using the hyperopia categorical trait, with the most significant SNPs rs1357179 on 15q14 (p=1.69×10⁻³) and rs7164400 on 15q25 (p=8.39×10⁻⁴), which passed the replication thresholds.

CONCLUSIONS

This study adds to the growing body of evidence that attempting to replicate the most significant SNPs found in one population may not be significant in another population due to differences in the linkage disequilibrium structure and/or allele frequency. This suggests that replication studies should include less significant SNPs in an associated region rather than only a few selected SNPs chosen by a significance threshold.

Association study of IGF1 polymorphisms with susceptibility to high myopia in a Japanese population

PURPOSE

Polymorphisms in the insulin-like growth factor 1 (IGF1) gene were previously associated with high or extreme myopia in Caucasian and Chinese populations. In the present study, we investigated whether IGF1 polymorphisms are associated with high myopia in a Japanese population.

METHODS

A total of 446 Japanese patients with high myopia (≤-9.00 diopters) and 481 Japanese healthy controls (+1.50 diopters to -1.50 diopters) were recruited. We genotyped seven tagging single-nucleotide polymorphisms (SNPs) in IGF1 and assessed allelic and haplotypic diversity in cases and controls.

RESULTS

There were no statistically significant differences in the allele frequencies of IGF1 SNPs and genotypes between cases and controls (P>0.05). However, the A allele of rs5742629 and the G allele of rs12423791 were associated with a moderately increased risk of high myopia (odds ratio [OR] =1.20 and OR =1.21, respectively) with borderline statistical significance (P=0.0502, corrected P (Pc) =0.21 and P=0.064, Pc=0.29, respectively). The haplotype consisting of the A allele of rs5742629 and the G allele of rs12423791 was marginally associated with the risk of high myopia (P=0.041; OR =1.21); this association was not significant after correction (Pc=0.19).

CONCLUSION

We found that the IGF1 SNPs are not significantly associated with high myopia in our Japanese population. Our results are in contrast to a previous study in which extreme myopia cases had significantly higher frequencies of the G allele of rs5742629 and the C allele of rs12423791 than controls. Therefore, the IGF1 SNPs may not be important factors for susceptibility to high myopia in all populations. Further genetic studies are needed to elucidate the possible contributions of the IGF1 region to the development of high myopia.

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.

Meta-analysis of genome-wide association studies in five cohorts reveals common variants in RBFOX1, a regulator of tissue-specific splicing, associated with refractive error

Visual refractive errors (REs) are complex genetic traits with a largely unknown etiology. To date, genome-wide association studies (GWASs) of moderate size have identified several novel risk markers for RE, measured here as mean spherical equivalent (MSE). We performed a GWAS using a total of 7280 samples from five cohorts: the Age-Related Eye Disease Study (AREDS); the KORA study ('Cooperative Health Research in the Region of Augsburg'); the Framingham Eye Study (FES); the Ogliastra Genetic Park-Talana (OGP-Talana) Study and the Multiethnic Study of Atherosclerosis (MESA). Genotyping was performed on Illumina and Affymetrix platforms with additional markers imputed to the HapMap II reference panel. We identified a new genome-wide significant locus on chromosome 16 (rs10500355, P = 3.9 × 10(-9)) in a combined discovery and replication set (26 953 samples). This single nucleotide polymorphism (SNP) is located within the RBFOX1 gene which is a neuron-specific splicing factor regulating a wide range of alternative splicing events implicated in neuronal development and maturation, including transcription factors, other splicing factors and synaptic proteins.