Identification of myopia-associated WNT7B polymorphisms provides insights into the mechanism underlying the development of myopia

Identification of LRRC46 as a novel candidate gene for high myopia

High myopia (HM) is the primary cause of blindness, with the microstructural organization and composition of collagenous fibers in the cornea and sclera playing a crucial role in the biomechanical behavior of these tissues. In a previously reported myopic linkage region, MYP5 (17q21-22), a potential candidate gene, LRRC46 (c.C235T, p.Q79X), was identified in a large Han Chinese pedigree. LRRC46 is expressed in various eye tissues in humans and mice, including the retina, cornea, and sclera. In subsequent cell experiments, the mutation (c.C235T) decreased the expression of LRRC46 protein in human corneal epithelial cells (HCE-T). Further investigation revealed that Lrrc46-/- mice (KO) exhibited a classical myopia phenotype. The thickness of the cornea and sclera in KO mice became thinner and more pronounced with age, the activity of limbal stem cells decreased, and microstructural changes were observed in the fibroblasts of the sclera and cornea. We performed RNA-seq on scleral and corneal tissues of KO and normal control wild-type (WT) mice, which indicated a significant downregulation of the collagen synthesis-related pathway (extracellular matrix, ECM) in KO mice. Subsequent in vitro studies further indicated that LRRC46, a member of the important LRR protein family, primarily affected the formation of collagens. This study suggested that LRRC46 is a novel candidate gene for HM, influencing collagen protein VIII (Col8a1) formation in the eye and gradually altering the biomechanical structure of the cornea and sclera, thereby promoting the occurrence and development of HM.

Age and light damage influence Fzd5 regulation of ocular growth-related genes

Genetic and environmental factors can independently or coordinatively drive ocular axis growth. Mutations in FRIZZLED5 (FZD5) have been associated with microphthalmia, coloboma, and, more recently, high myopia. The molecular mechanism of how Fzd5 participates in ocular growth remains unknown. In this study, we compiled a list of human genes associated with ocular growth abnormalities based on public databases and a literature search. We identified a set of ocular growth-related genes from the list that was altered in the Fzd5 mutant mice by RNAseq analysis at different time points. The Fzd5 regulation of this set of genes appeared to be impacted by age and light damage. Further bioinformatical analysis indicated that these genes are extracellular matrix (ECM)-related; and meanwhile an altered Wnt signaling was detected. Altogether, the data suggest that Fzd5 may regulate ocular growth through regulating ECM remodeling, hinting at a genetic-environmental interaction in gene regulation of ocular axis control.

A multiethnic genome-wide analysis of 19,420 individuals identifies novel loci associated with axial length and shared genetic influences with refractive error and myopia

Introduction: Long axial length (AL) is a risk factor for myopia. Although family studies indicate that AL has an important genetic component with heritability estimates up to 0.94, there have been few reports of AL-associated loci. Methods: Here, we conducted a multiethnic genome-wide association study (GWAS) of AL in 19,420 adults of European, Latino, Asian, and African ancestry from the Genetic Epidemiology Research on Adult Health and Aging (GERA) cohort, with replication in a subset of the Consortium for Refractive Error and Myopia (CREAM) cohorts of European or Asian ancestry. We further examined the effect of the identified loci on the mean spherical equivalent (MSE) within the GERA cohort. We also performed genome-wide genetic correlation analyses to quantify the genetic overlap between AL and MSE or myopia risk in the GERA European ancestry sample. Results: Our multiethnic GWA analysis of AL identified a total of 16 genomic loci, of which 5 are novel. We found that all AL-associated loci were significantly associated with MSE after Bonferroni correction. We also found that AL was genetically correlated with MSE (rg = -0.83; SE, 0.04; p = 1.95 × 10-89) and myopia (rg = 0.80; SE, 0.05; p = 2.84 × 10-55). Finally, we estimated the array heritability for AL in the GERA European ancestry sample using LD score regression, and found an overall heritability estimate of 0.37 (s.e. = 0.04). Discussion: In this large and multiethnic study, we identified novel loci, associated with AL at a genome-wide significance level, increasing substantially our understanding of the etiology of AL variation. Our results also demonstrate an association between AL-associated loci and MSE and a shared genetic basis between AL and myopia risk.

A Genome-Wide Association Study for Susceptibility to Axial Length in Highly Myopic Eyes

High myopia has long been highly prevalent worldwide with a largely yet unexplained genetic contribution. To identify novel susceptibility genes for axial length (AL) in highly myopic eyes, a genome-wide association study (GWAS) was performed using the genomic dataset of 350 deep whole-genome sequencing data from highly myopic patients. Top single nucleotide polymorphisms (SNPs) were functionally annotated. Immunofluorescence staining, quantitative polymerase chain reaction, and western blot were performed using neural retina of form-deprived myopic mice. Enrichment analyses were further performed. We identified the four top SNPs and found that ADAM Metallopeptidase With Thrombospondin Type 1 Motif 16 (ADAMTS16) and Phosphatidylinositol Glycan Anchor Biosynthesis Class Z (PIGZ) had the potential of clinical significance. Animal experiments confirmed that PIGZ expression could be observed and showed higher expression level in form-deprived mice, especially in the ganglion cell layer. The messenger RNA (mRNA) levels of both ADAMTS16 and PIGZ were significantly higher in the neural retina of form-deprived eyes (p = 0.005 and 0.007 respectively), and both proteins showed significantly upregulated expression in the neural retina of deprived eyes (p = 0.004 and 0.042, respectively). Enrichment analysis revealed a significant role of cellular adhesion and signal transduction in AL, and also several AL-related pathways including circadian entrainment and inflammatory mediator regulation of transient receptor potential channels were proposed. In conclusion, the current study identified four novel SNPs associated with AL in highly myopic eyes and confirmed that the expression of ADAMTS16 and PIGZ was significantly upregulated in neural retina of deprived eyes. Enrichment analyses provided novel insight into the etiology of high myopia and opened avenues for future research interest.

Supplementary Information

The online version contains supplementary material available at 10.1007/s43657-022-00082-x.

Shedding light on myopia by studying complete congenital stationary night blindness

Myopia is the most common eye disorder, caused by heterogeneous genetic and environmental factors. Rare progressive and stationary inherited retinal disorders are often associated with high myopia. Genes implicated in myopia encode proteins involved in a variety of biological processes including eye morphogenesis, extracellular matrix organization, visual perception, circadian rhythms, and retinal signaling. Differentially expressed genes (DEGs) identified in animal models mimicking myopia are helpful in suggesting candidate genes implicated in human myopia. Complete congenital stationary night blindness (cCSNB) in humans and animal models represents an ON-bipolar cell signal transmission defect and is also associated with high myopia. Thus, it represents also an interesting model to identify myopia-related genes, as well as disease mechanisms. While the origin of night blindness is molecularly well established, further research is needed to elucidate the mechanisms of myopia development in subjects with cCSNB. Using whole transcriptome analysis on three different mouse models of cCSNB (in Gpr179^-/-, Lrit3^-/- and Grm6^-/-), we identified novel actors of the retinal signaling cascade, which are also novel candidate genes for myopia. Meta-analysis of our transcriptomic data with published transcriptomic databases and genome-wide association studies from myopia cases led us to propose new biological/cellular processes/mechanisms potentially at the origin of myopia in cCSNB subjects. The results provide a foundation to guide the development of pharmacological myopia therapies.

Bi-allelic variants in WNT7B disrupt the development of multiple organs in humans

BACKGROUND

Pulmonary hypoplasia, Diaphragmatic anomalies, Anophthalmia/microphthalmia and Cardiac defects delineate the PDAC syndrome. We aim to identify the cause of PDAC syndrome in patients who do not carry pathogenic variants in RARB and STRA6, which have been previously associated with this disorder.

METHODS

We sequenced the exome of patients with unexplained PDAC syndrome and performed functional validation of candidate variants.

RESULTS

We identified bi-allelic variants in WNT7B in fetuses with PDAC syndrome from two unrelated families. In one family, the fetus was homozygous for the c.292C>T (p.(Arg98*)) variant whereas the fetuses from the other family were compound heterozygous for the variants c.225C>G (p.(Tyr75*)) and c.562G>A (p.(Gly188Ser)). Finally, a molecular autopsy by proxy in a consanguineous couple that lost two babies due to lung hypoplasia revealed that both parents carry the p.(Arg98*) variant. Using a WNT signalling canonical luciferase assay, we demonstrated that the identified variants are deleterious. In addition, we found that wnt7bb mutant zebrafish display a defect of the swimbladder, an air-filled organ that is a structural homolog of the mammalian lung, suggesting that the function of WNT7B has been conserved during evolution for the development of these structures.

CONCLUSION

Our findings indicate that defective WNT7B function underlies a form of lung hypoplasia that is associated with the PDAC syndrome, and provide evidence for involvement of the WNT-β-catenin pathway in human lung, tracheal, ocular, cardiac, and renal development.

Association between central serous chorioretinopathy susceptibility genes and choroidal parameters

PURPOSE

To evaluate the association between central serous chorioretinopathy (CSC) susceptibility genes and choroidal parameters in a large Japanese cohort.

STUDY DESIGN

Retrospective cohort study.

METHODS

Of the 9850 individuals in the Nagahama study whose second visit was between 2013 and 2016, those with optical coherence tomography (OCT) images with enhanced depth imaging (EDI), axial length, and genome-wide single nucleotide polymorphism (SNP) genotyping data were included. We calculated subfoveal choroidal thickness (SFCT), choroidal vascularity index (CVI), normalized choroidal intensity (NCI), and vertical asymmetry of choroidal thickness. Genome-wide quantitative trait locus (QTL) analyses were performed for each parameter. We screened for four CSC susceptibility SNPs: CFH rs800292, TNFRSF10A rs13278062, GATA5 rs6061548, and VIPR2 rs3793217. Whenever an SNP was not included in the genotyping data after quality control, its proxy SNP was selected.

RESULTS

In total, 4586 participants were evaluated. CFH rs800292 was significantly associated with SFCT (P < 0.001) and CVI (P < 0.001). VIPR2 rs3793217 was significantly associated with SFCT (P < 0.001) but not with CVI. Whereas, TNFRSF10A rs13254617 and GATA5 rs6061548 were not significantly associated with SFCT or CVI. None of these SNPs was associated with NCI_EDI and asymmetry of choroidal thickness.

CONCLUSION

CFH, VIPR2, TNFRSF10A, and GATA5 showed different association patterns with choroidal parameters. Although the mechanism of CSC pathogenesis by choroidal changes is not fully understood, this finding suggests that each gene may be involved in different mechanisms of CSC development. Our genetic study provides a basis for understanding the role of CSC susceptibility genes.

Retinal Transcriptomics Analysis Reveals the Underlying Mechanism of Disturbed Emmetropization Induced by Wavelength Defocus

PURPOSE

Wavelength signals play a vital role in refractive development. This study aimed to explore the retinal transcriptome signature in these processes.

METHODS

Guinea pigs were randomly divided into three groups exposed to white, blue, or green environmental light for eight weeks. Refraction and axial length were evaluated every 4 weeks, and the retinal transcriptome was profiled at 8 weeks.

RESULTS

Compared with the white group, ocular refraction significantly decreased and ocular axial length significantly extended in the green group whereas these parameters showed opposite trends in the blue group. RNA-sequencing showed that, compared with the white group, 184 and 171 differentially expressed genes (DEGs) were found in the blue and green groups, respectively. Among these DEGs, only 31 overlapped. These two sets of DEGs were enriched in distinct biological processes and pathways. There were 268 DEGs between the blue and green groups, which were primarily enriched in the extracellular matrix, and metabolism, receptor activity, and ion binding processes. In addition, nine human genes, including ECEL1, CHRND, SHBG, PRSS56, OVOL1, RDH5, WNT7B, PEBP4, CA12, were identified to be related to myopia development and wavelength response, indicating the potential role of these genes in human wavelength-induced myopia.

CONCLUSIONS

In this study, we identified retinal targets and pathways involved in the response to wavelength signals in emmetropization.

Candidate pathways for retina to scleral signaling in refractive eye growth

The global prevalence of myopia, or nearsightedness, has increased at an alarming rate over the last few decades. An eye is myopic if incoming light focuses prior to reaching the retinal photoreceptors, which indicates a mismatch in its shape and optical power. This mismatch commonly results from excessive axial elongation. Important drivers of the myopia epidemic include environmental factors, genetic factors, and their interactions, e.g., genetic factors influencing the effects of environmental factors. One factor often hypothesized to be a driver of the myopia epidemic is environmental light, which has changed drastically and rapidly on a global scale. In support of this, it is well established that eye size is regulated by a homeostatic process that incorporates visual cues (emmetropization). This process allows the eye to detect and minimize refractive errors quite accurately and locally over time by modulating the rate of elongation of the eye via remodeling its outermost coat, the sclera. Critically, emmetropization is not dependent on post-retinal processing. Thus, visual cues appear to influence axial elongation through a retina-to-sclera, or retinoscleral, signaling cascade, capable of transmitting information from the innermost layer of the eye to the outermost layer. Despite significant global research interest, the specifics of retinoscleral signaling pathways remain elusive. While a few pharmacological treatments have proven to be effective in slowing axial elongation (most notably topical atropine), the mechanisms behind these treatments are still not fully understood. Additionally, several retinal neuromodulators, neurotransmitters, and other small molecules have been found to influence axial length and/or refractive error or be influenced by myopigenic cues, yet little progress has been made explaining how the signal that originates in the retina crosses the highly vascular choroid to affect the sclera. Here, we compile and synthesize the evidence surrounding three of the major candidate pathways receiving significant research attention - dopamine, retinoic acid, and adenosine. All three candidates have both correlational and causal evidence backing their involvement in axial elongation and have been implicated by multiple independent research groups across diverse species. Two hypothesized mechanisms are presented for how a retina-originating signal crosses the choroid - via 1) all-trans retinoic acid or 2) choroidal blood flow influencing scleral oxygenation. Evidence of crosstalk between the pathways is discussed in the context of these two mechanisms.

Myopia Genetics and Heredity

Myopia is the most common eye condition leading to visual impairment and is greatly influenced by genetics. Over the last two decades, more than 400 associated gene loci have been mapped for myopia and refractive errors via family linkage analyses, candidate gene studies, genome-wide association studies (GWAS), and next-generation sequencing (NGS). Lifestyle factors, such as excessive near work and short outdoor time, are the primary external factors affecting myopia onset and progression. Notably, besides becoming a global health issue, myopia is more prevalent and severe among East Asians than among Caucasians, especially individuals of Chinese, Japanese, and Korean ancestry. Myopia, especially high myopia, can be serious in consequences. The etiology of high myopia is complex. Prediction for progression of myopia to high myopia can help with prevention and early interventions. Prediction models are thus warranted for risk stratification. There have been vigorous investigations on molecular genetics and lifestyle factors to establish polygenic risk estimations for myopia. However, genes causing myopia have to be identified in order to shed light on pathogenesis and pathway mechanisms. This report aims to examine current evidence regarding (1) the genetic architecture of myopia; (2) currently associated myopia loci identified from the OMIM database, genetic association studies, and NGS studies; (3) gene-environment interactions; and (4) the prediction of myopia via polygenic risk scores (PRSs). The report also discusses various perspectives on myopia genetics and heredity.

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.

Genome-wide Association Study of Axial Length in Population-based Cohorts in Japan

Purpose

To elucidate the differences in ocular biometric parameters by generation and gender and to identify axial length (AL)-associated genetic variants in Japanese individuals, we analyzed Tohoku Medical Megabank Organization (ToMMo) Eye Study data.

Design

We designed the ToMMo Eye Study, examined AL variations, and conducted genome-wide association studies (GWASs).

Participants

In total, 33 483 participants aged > 18 years who were recruited into the community-based cohort (CommCohort) and the birth and three-generation cohort (BirThree Cohort) of the ToMMo Eye Study were examined.

Methods

Each participant was screened with an interview, ophthalmic examinations, and a microarray analysis. The GWASs were performed in 22 379 participants in the CommCohort (discovery stage) and 11 104 participants in the BirThree Cohort (replication stage). We evaluated the associations of single nucleotide polymorphisms (SNPs) with AL using a genome-wide significance threshold (5 × 10^-8) in each stage of the study and in the subsequent meta-analysis.

Main Outcome Measures

We identified the association of SNPs with AL and distributions of AL in right and left eyes and individuals of different sexes and ages.

Results

In the discovery stage, the mean AL of the right eye (23.99 mm) was significantly greater than that of the left eye (23.95 mm). This difference was reproducible across sexes and ages. The GWASs revealed 703 and 215 AL-associated SNPs with genome-wide significance in the discovery and validation stages, respectively, and many of the SNPs in the discovery stage were replicated in the validation stage. Validated SNPs and their associated loci were meta-analyzed for statistical significance (P < 5 × 10^-8). This study identified 1478 SNPs spread over 31 loci. Of the 31 loci, 5 are known AL loci, 15 are known refractive-error loci, 4 are known corneal-curvature loci, and 7 loci are newly identified loci that are not known to be associated with AL. Of note, some of them shared functional relationships with previously identified loci.

Conclusions

Our large-scale GWASs exploiting ToMMo Eye Study data identified 31 loci linked to variations in AL, 7 of which are newly reported in this article. The results revealed genetic heterogeneity and similarity in SNPs related to ethnic variations in AL.

Identification of the Potential Key Genes and Pathways Involved in Lens Changes of High Myopia

Aim

Methods

Results

Conclusion

Retinal Nerve Fiber Layer Thickness and Rim Area Profiles in Asians: Pooled Analysis from the Asian Eye Epidemiology Consortium

PURPOSE

To evaluate ethnic variations, ocular and systemic determinants of retinal nerve fibre layer (RNFL) thickness and neuroretinal rim area among Asians, using a large consortium of population-based eye studies.

DESIGN

Cross-sectional pooled-analysis.

PARTICIPANTS

22436 participants (22436 eyes) from 10 population-based studies (China, Hong Kong, India, Japan, Russia and Singapore) of the Asian Eye Epidemiology Consortium.

METHODS

Participants aged ≥40 years without glaucoma were included. All participants underwent spectral domain optical coherence tomography (OCT) imaging, systemic and ocular examinations. Data were pooled from each study. Multivariable regression analysis was performed to evaluate inter-ethnic, inter-machine variations, ocular and systemic factors associated with RNFL thickness and rim area, adjusting for age, gender, diabetes, intraocular pressure (IOP), spherical equivalent (SE), ethnicity, OCT model, and study group. When evaluating body mass index, smoking, and hypertension as exposures, these factors were additionally adjusted in the model.

MAIN OUTCOME MEASURE

Average RNFL thickness (μm) and rim area (mm²) RESULTS: Indian and Japanese eyes showed thinner RNFL, compared to other Asian ethnicities (β values ranging 7.31-12.76μm, P<0.001 for all pair-wise comparisons). Compared to measurements by Cirrus HD-OCT, RNFL was on average 7.29μm thicker when measured by Spectralis, 12.85μm thicker by Nidek, and 17.48μm thicker by Optovue (all P<0.001). Additionally, older age (per decade, β=-2.70; 95% confidence interval [CI], -2.85 to -2.55), diabetes (β=-0.72; 95%CI, -1.20 to -0.24), ), higher IOP (per mmHg, β=-0.07; 95% CI, -0.10 to -0.04), more myopic SE (per dioptre, β=-1.13; 95% CI, -1.19 to -1.07), cardiovascular disease (CVD, β=-0.94; 95% CI, -1.49 to -0.40), and hypertension (β=-0.68; 95% CI, -1.04 to -0.32), were associated with thinner RNFL (all P≤0.003). Similarly, older age (β=-0.019; 95% CI, -0.028 to -0.009), higher IOP (β=-0.010; 95% CI, -0.013 to -0.008) and more myopic SE (β=-0.025; 95% CI, -0.029 to -0.021) were associated with smaller rim area (all P<0.001).

CONCLUSIONS

In this large pooled-analysis of multiple Asian population studies, Indian and Japanese eyes were observed to have thinner RNFL profiles. In addition to previously known determinants, hypertension and CVD were associated with thinner RNFL. These findings further suggest the need of ethnic-specific normative database to improve glaucoma detection.

Genetic Variants Associated With Human Eye Size Are Distinct From Those Conferring Susceptibility to Myopia

Purpose

Emmetropization requires coordinated scaling of the major ocular components, corneal curvature and axial length. This coordination is achieved in part through a shared set of genetic variants that regulate eye size. Poorly coordinated scaling of corneal curvature and axial length results in refractive error. We tested the hypothesis that genetic variants regulating eye size in emmetropic eyes are distinct from those conferring susceptibility to refractive error.

Methods

A genome-wide association study (GWAS) for corneal curvature in 22,180 adult emmetropic individuals was performed as a proxy for a GWAS for eye size. A polygenic score created using lead GWAS variants was tested for association with corneal curvature and axial length in an independent sample: 437 classified as emmetropic and 637 as ametropic. The genetic correlation between eye size and refractive error was calculated using linkage disequilibrium score regression for approximately 1 million genetic variants.

Results

The GWAS for corneal curvature in emmetropes identified 32 independent genetic variants (P < 5.0e-08). A polygenic score created using these 32 genetic markers explained 3.5% (P < 0.001) and 2.0% (P = 0.001) of the variance in corneal curvature and axial length, respectively, in the independent sample of emmetropic individuals but was not predictive of these traits in ametropic individuals. The genetic correlation between eye size and refractive error was close to zero (r_g = 0.00; SE = 0.06; P = 0.95).

Conclusions

These results support the hypothesis that genetic variants regulating eye size in emmetropic eyes do not overlap with those conferring susceptibility to myopia. This suggests that distinct biological pathways regulate normal eye growth and myopia development.

Corneal curvature-associated MTOR variant differentiates mild myopia from high myopia in Han Chinese population

BACKGROUND

Myopia is the most prevalent ocular disorder in the world, and corneal parameters have been regarded as key ocular biometric parameters determining the refractive status. Here, we aimed to determine the association of genome-wide association study-identified corneal curvature (CC)-related gene variants with different severity of myopia and ocular biometric parameters in Chinese population.

METHODS

Total 2,101 unrelated Han Chinese subjects were recruited, including 1,649 myopia and 452 control subjects. Five previously reported CC-associated gene variants (PDGFRA, MTOR, WNT7B, CMPK1 and RBP3) were genotyped by TaqMan assay, and their association with different myopia severity and ocular biometric parameters were evaluated.

RESULTS

Joint additive effect analysis showed that MTOR rs74225573 paired with PDGFRA rs2114039 (P = .009, odds ratio (OR) = 4.91) or CMPK1 rs17103186 (P = .002, OR = 13.03) were significantly associated with higher risk in mild myopia. Critically, mild myopia subjects had significantly higher frequency in MTOR rs74225573 C allele than high myopia subjects (P = .003), especially in male subjects (P = .001, OR = 0.49). High myopia subjects carrying MTOR rs74225573 C allele have significant flatter CC (P = .035) and longer corneal radius (P = .044) than those carrying TT genotype.

CONCLUSION

This study revealed that male high myopia subjects are more prone to carry CC-related MTOR rs74225573 T allele, whereas mild myopia subjects are prone to carry the C allele. MTOR rs7422573 variant could be a genetic marker to differentiate mild from high myopia in risk assessment.

ABBREVIATIONS

ACD: anterior chamber depth; AL: axial length; AL/CR: axial length/corneal radius ratio; ANOVA: analysis of variance; CC: corneal curvature; CCT: central corneal thickness; C.I.: confidence interval; CMPK1: cytidine/uridine monophosphate kinase 1; CR: corneal radius; D: diopter; GWAS: genome-wide association studies; HWE: Hardy-Weinberg equilibrium; LT: lens thickness; MIPEP: mitochondrial intermediate peptidase; MTOR: mechanistic target of rapamycin kinase; OR: odds ratio; PDGFRA: platelet-derived growth factor receptor-α; RBP3: retinol-binding protein 3; SD: standard deviation; SE: spherical equivalence; SNTB1: syntrophin beta 1; VCD: vitreous chamber depth; VIPR2: vasoactive intestinal peptide receptor 2; WNT7B: wingless/integrated family member 7B.

Canonical Wnt Signaling Drives Myopia Development and Can Be Pharmacologically Modulated

Purpose

The purpose of this study was to evaluate the role of the canonical Wnt signaling in the development of the myopia.

Methods

Plasma from adult patients with myopia, myopic animal models including the adenomatous polyposis coli (APC) gene mutation mouse model, and the form deprivation (FD) induced mouse model of myopia were used. Niclosamide, a canonical Wnt pathway inhibitor, was orally administrated in animal models. Plasma levels of DKK-1 were determined by using enzyme-linked immunosorbent assay. Refraction, vitreous chamber depth (VCD), axial length (AL), and other parameters, were measured at the end of the FD treatment. Canonical Wnt signaling changes were evaluated by Western blot analysis and immunostaining analysis.

Results

Plasma level of Wnt inhibitor DKK-1 was markedly decreased in patients with myopia. Meanwhile, the canonical Wnt pathway was progressively activated during myopia development in mice. Moreover, inhibition of canonical Wnt signaling by niclosamide in mouse models markedly reduced lens thickness (LT), VCD, and AL elongation, resulting in myopia inhibition.

Conclusions

Dysregulation of canonical Wnt signaling is a characteristic of myopia and targeting Wnt signaling pathways has potential as a therapeutic strategy for myopia.

Distribution of Choroidal Thickness and Choroidal Vessel Dilation in Healthy Japanese Individuals

Purpose

To report fundamental epidemiologic data for choroidal parameters such as choroidal thickness and index of choroidal vascularity in Japanese individuals and to evaluate their correlations with age, sex, systemic parameters, and other ocular parameters.

Design

Population-based cohort study.

Participants

A total of 9850 individuals participated in the first follow-up of the Nagahama Prospective Cohort for Comprehensive Human Bioscience (the Nagahama Study) conducted between 2013 and 2016.

Methods

All participants underwent standardized ophthalmic examinations, including OCT with enhanced depth imaging (EDI; RS-3000 Advance; Nidek). We manually segmented the choroidoscleral interface to measure subfoveal choroidal thickness (SFCT) and calculated the normalized choroidal intensity obtained with EDI (NCI_EDI) and choroidal vascularity index (CVI). These are indices of choroidal brightness in OCT and reportedly represent the dilation of choroidal vessels. After summarizing the age-sex stratified distributions of SFCT, NCI_EDI, and CVI, their associations with age, sex, axial length (AL), and spherical equivalent (SE) were evaluated using linear regression analysis with adjustments for possible confounders.

Main Outcome Measures

Distribution of SFCT, NCI_EDI, and CVI in the healthy Japanese population and their characteristics.

Results

Age-sex standardized SFCT, NCI_EDI, and CVI were 291.2 μm, 0.653, and 66.88%, respectively. In both men and women, SFCT was associated negatively with age (P < 0.001) and NCI_EDI was associated positively with age (P < 0.001). Although both SFCT and NCI_EDI did not differ significantly between men and women overall (P = 0.87 and P = 0.21, respectively), among younger participants (35–50 years of age), men showed significantly greater SFCT than women (P < 0.001). Only in men was CVI associated positively with age (P < 0.001). In the multivariable analysis, SFCT was associated significantly with age, sex, AL, SE, and the interaction term of age and sex (P < 0.001). Independent of SFCT, NCI_EDI and CVI were associated significantly with age (P < 0.001).

Conclusions

We report normative Japanese SFCT, NCI_EDI, and CVI data using a large general Japanese cohort. The association analysis of SFCT with NCI_EDI and CVI suggested that younger individuals have a more lumen-rich choroid for their choroidal thickness than older individuals.

The effect of Wnt/β-catenin pathway on the scleral remolding in the mouse during form deprivation

BACKGROUND

Many reports have shown that Wnt/β-Catenin pathway is associated with a variety of diseases, but its role in the pathogenesis of myopia is still unknown. In order to clarify the role of Wnt/β-catenin pathway in the development of form deprivation myopia (FDM), this study investigated the expression of scleral Wls, β-catenin and TCF4 in mice model of form deprivation (FD) myopia.

METHODS

Three parallel experimental groups, including FD, monocular exposure (SC) and binocular exposure (NC) group, were designed to investigate the effects of Wnt/β-Catenin pathway on scleral remodeling mouse during form deprivation in three-week-old C57BL/6 mice. Diopters and axial lengths (AL) in each sample were measured with an infrared eccentric refractometer or spectral-domain optical coherence tomography. The expression of scleral Wls, β-catenin and TCF4 were detected with Western blot. Morphological changes of posterior sclera were observed with a transmission electron microscope (TEM). The above characterization and analysis were performed on the 0th, 7th, 14th, 21st and 28th days, respectively.

RESULTS

The difference of diopter and AL between the three groups (SC, NC and FD group) gradually increased with the prolongation of FD time (except AL between SC and NC groups). The changes of diopter and AL gradually increased with the prolongation of FD time. Especially, the diopter and AL will increase sharply after FD lasts for a long time, such as the measurement on the 21st for diopter and 28th days for AL. Most notably, the AL of FD eyes significantly increased after 28 days of deprivation. Thinning and disordered arrangement of collagen fibers and a decrease of extracellular matrix were observed with TEM. The expression of scleral Wls, β-catenin and TCF4 increased with age in the NC and SC group. In FD group, they increased significantly on the 7th, 14th and 21st days but decreased on the 28th day.

CONCLUSIONS

The expression of Wls, β-Catenin and TCF4 to FD were more sensitive indicators than that of diopter and AL. Within the first 7 days of FD, the expression of Wls, β-Catenin and TCF4 in sclera increased sharply. With the extension of FD duration, it gradually decreased. It is suggested that the Wnt/β-catenin pathway might be involved in the scleral remodeling induced in FDM mice.

Association of Myopia and Genetic Variants of TGFB2-AS1 and TGFBR1 in the TGF-β Signaling Pathway: A Longitudinal Study in Chinese School-Aged Children

Background

Myopia is a complex multifactorial condition which involves several overlapping signaling pathways mediated by distinct genes. This prospective cohort study evaluated the associations of two genetic variants in the TGF-β signaling pathway with the onset and progression of myopia and ocular biometric parameters in Chinese school-aged children.

Methods

A total of 556 second grade children were examined and followed up for 3.5 years. Non-cycloplegic refraction and ocular biometric parameters were measured annually. Multivariate regression analysis was used to assess the effect of the TGFBR1 rs10760673 and TGFB2-AS1 rs7550232 variants on the occurrence and progression of myopia. A 10,000 permutations test was used to correct for multiple testing. Functional annotation of single nucleotide polymorphisms (SNPs) was performed using RegulomeDB, HaploReg, and rVarBase.

Results

A total of 448 children were included in the analysis. After adjustments for gender, age, near work time and outdoor time with 10,000 permutations, the results indicated that the C allele and the AC or CC genotypes of rs7550232 adjacent to TGFB2-AS1 were associated with a significantly increased risk of the onset of myopia in two genetic models (additive: P’ = 0.022; dominant: P’ = 0.025). Additionally, the A allele and the AA or AG genotypes of rs10760673 of TGFBR1 were associated with a significant myopic shift (additive: P’ = 0.008; dominant: P’ = 0.028; recessive: P’ = 0.027). Furthermore, rs10760673 was associated with an increase in axial length (AL) (P’ = 0.013, β = 0.03) and a change in the ratio of AL to the corneal radius of curvature (AL/CRC) (P’ = 0.031, β = 0.003). Analysis using RegulomeDB, HaploReg, and rVarBase indicated that rs7550232 is likely to affect transcription factor binding, any motif, DNase footprint, and DNase peak.

Conclusion

The present study indicated that rs10760673 and rs7550232 may represent susceptibility loci for the progression and onset of myopia, respectively, in school-aged children. Associations of the variants of the TGFBR1 and TGFB2-AS1 genes with myopia may be mediated by the TGF-β signaling pathway; this hypothesis requires validation in functional studies. This trial was registered as ChiCTR1900020584 at www.Chictr.org.cn.

Aberrant TGF-β1 signaling activation by MAF underlies pathological lens growth in high myopia

High myopia is a leading cause of blindness worldwide. Myopia progression may lead to pathological changes of lens and affect the outcome of lens surgery, but the underlying mechanism remains unclear. Here, we find an increased lens size in highly myopic eyes associated with up-regulation of β/γ-crystallin expressions. Similar findings are replicated in two independent mouse models of high myopia. Mechanistic studies show that the transcription factor MAF plays an essential role in up-regulating β/γ-crystallins in high myopia, by direct activation of the crystallin gene promoters and by activation of TGF-β1-Smad signaling. Our results establish lens morphological and molecular changes as a characteristic feature of high myopia, and point to the dysregulation of the MAF-TGF-β1-crystallin axis as an underlying mechanism, providing an insight for therapeutic interventions.

High myopia is associated with lens changes, but the underlying mechanisms are unclear. Here, the authors show increased equatorial diameter of the lens in subjects affected by high myopia, and find that these changes are associated with an increase in crystallin expression driven by the transcription factor MAF and TGF-β1 signaling.

RNA-seq and GSEA identifies suppression of ligand-gated chloride efflux channels as the major gene pathway contributing to form deprivation myopia

Currently there is no consensus regarding the aetiology of the excessive ocular volume that characterizes high myopia. Thus, we aimed to test whether the gene pathways identified by gene set enrichment analysis of RNA-seq transcriptomics refutes the predictions of the Retinal Ion Driven Efflux (RIDE) hypothesis when applied to the induction of form-deprivation myopia (FDM) and subsequent recovery (post-occluder removal). We found that the induction of profound FDM led to significant suppression in the ligand-gated chloride ion channel transport pathway via suppression of glycine, GABA_A and GABA_C ionotropic receptors. Post-occluder removal for short term recovery from FDM of 6 h and 24 h, induced significant upregulation of the gene families linked to cone receptor phototransduction, mitochondrial energy, and complement pathways. These findings support a model of form deprivation myopia as a Cl⁻ ion driven adaptive fluid response to the modulation of the visual signal cascade by form deprivation that in turn affects the resultant ionic environment of the outer and inner retinal tissues, axial and vitreal elongation as predicted by the RIDE model. Occluder removal and return to normal light conditions led to return to more normal upregulation of phototransduction, slowed growth rate, refractive recovery and apparent return towards physiological homeostasis.

Relationship between Intraocular Pressure and Coffee Consumption in a Japanese Population without Glaucoma: The Nagahama Study

PURPOSE

To evaluate the association between daily coffee consumption and intraocular pressure (IOP) in healthy persons without glaucoma and the association between daily coffee consumption and history of glaucoma.

DESIGN

Cross-sectional study.

PARTICIPANTS

A total of 9850 individuals participated in the first follow-up of the Nagahama Prospective Cohort for Comprehensive Human Bioscience (the Nagahama Study) conducted between 2013 and 2016.

METHODS

All participants underwent a standardized ophthalmic examination. Self-reporting questionnaires were completed by all participants. First, the association between habitual coffee consumption and IOP among nonglaucoma individuals was evaluated by a multivariate linear regression analysis, adjusting for possible confounders. Second, the association between habitual coffee consumption and history of glaucoma also was evaluated using a multivariate logistic regression analysis.

MAIN OUTCOME MEASURES

The association between habitual coffee consumption and IOP among nonglaucoma individuals.

RESULTS

Of 9850 participants, 9418 did not have history of glaucoma. Among these participants, the mean ± standard deviation IOP of both eyes was 14.7 ± 2.9 mmHg. The multivariate regression analysis revealed that habitual coffee consumption was associated significantly with IOP (P < 0.001): the higher the consumption of coffee, the lower the IOP of an individual. The IOP of the group who consumed coffee most frequently (3 times daily or more) was 0.4 mmHg lower (95% confidence interval, 0.2-0.5 mmHg lower) than that of the group that consumed coffee least frequently (less than once daily). However, the logistic regression analysis showed that habitual coffee consumption was not associated significantly with history of glaucoma (P = 0.53).

CONCLUSIONS

Frequent coffee consumption was associated with a slightly lower IOP in people without glaucoma but was not associated with a decreased risk of glaucoma developing. Additional experimental studies are needed to examine the effects of coffee on IOP and glaucoma risk.

Association of WNT7B and RSPO1 with Axial Length in School Children

Purpose

To evaluate the association between single-nucleotide polymorphisms (SNPs) in the ZC3H11B, RSPO1, C3orf26, GJD2, ZNRF3, and WNT7B genes and myopia endophenotypes in children.

Methods

Seven SNPs identified in previous genome-wide association studies of axial length (AL) were genotyped in 2883 Southern Han Chinese children. Multiple linear regression analyses were conducted to evaluate the genotype association with AL, spherical equivalent (SE), corneal curvature (CC), and central corneal thickness (CCT).

Results

Two SNPs-namely, rs12144790 in RSPO1 (allele T, P = 0.0066, β = 0.062) and rs10453441 in WNT7B (allele A, P = 8.03 × 10-6, β = 0.103)-were significantly associated with AL. The association of rs4373767 in ZC3H11B (allele C, P = 0.030, β = -0.053) could not withstand the correction for multiple testing. WNT7B rs10453441 showed a strong association with CC (P = 1.17 × 10-14, β = 0.053) and with CCT (P = 0.0026, β = 2.65). None of the tested SNPs was significantly associated with SE. The C allele of SNP rs12321 in ZNRF3 was associated with CC (P = 0.0060, β = -0.018).

Conclusions

This study revealed that the RSPO1 SNP rs12144790 was associated with AL, whereas WNT7B rs10453441 was associated with AL, CC, and CCT in children. A novel association between ZNRF3 rs12321 and CC was discovered. Our data suggest that the RSPO1 and WNT7B genes might exert their effects on multiple aspects of eye growth during childhood. Potential differences in the genetic profiles of AL between children and adults should be explored in larger cohorts.

Genome-wide association study of individual differences of human lymphocyte profiles using large-scale cytometry data

Human immune systems are very complex, and the basis for individual differences in immune phenotypes is largely unclear. One reason is that the phenotype of the immune system is so complex that it is very difficult to describe its features and quantify differences between samples. To identify the genetic factors that cause individual differences in whole lymphocyte profiles and their changes after vaccination without having to rely on biological assumptions, we performed a genome-wide association study (GWAS), using cytometry data. Here, we applied computational analysis to the cytometry data of 301 people before receiving an influenza vaccine, and 1, 7, and 90 days after the vaccination to extract the feature statistics of the lymphocyte profiles in a nonparametric and data-driven manner. We analyzed two types of cytometry data: measurements of six markers for B cell classification and seven markers for T cell classification. The coordinate values calculated by this method can be treated as feature statistics of the lymphocyte profile. Next, we examined the genetic basis of individual differences in human immune phenotypes with a GWAS for the feature statistics, and we newly identified seven significant and 36 suggestive single-nucleotide polymorphisms associated with the individual differences in lymphocyte profiles and their change after vaccination. This study provides a new workflow for performing combined analyses of cytometry data and other types of genomics data.

Genetics vs chronic corneal mechanical trauma in the etiology of keratoconus

Both genetic and environmental factors have been considered to play a role in the etiology keratoconus. Eye rubbing, and more recently eye compression due to sleeping position, have been identified to be highly related to the condition, and are present in a high percentage of patients. Today, the predominant model is that these factors can provide the "second hit" necessary to generate the condition in a genetically susceptible individual. In addition, the extremely high prevalence in Arab populations, where endogamy could play a role, the high concordance rate in monozygotic twins, and the presence of family history of the condition between 5 and 23% of cases, support a genetic influence. Segregation analysis studies suggest that keratoconus is a complex non-Mendelian disease. Results from linkage analysis, next generation sequencing studies and genome-wide association studies also have suggested that genetic factors are involved in the condition. Recently, it has been proposed that mechanical trauma (i.e. eye rubbing or eye compression at night), is a sine quanon condition for the onset of keratoconus, and quite possibly its only cause. There are various arguments for and against this hypothesis. Indeed, it is possible, as initially suggested around 55 years ago, that the term "keratoconus" include diverse phenotypically similar conditions, which are actually of different etiology.

Cohort Profile: The Ganka-Ekigaku Network (GEN), a Network of Japanese Ophthalmological Epidemiology Studies

PURPOSE

Japan has been known as a super-aged society, and ageing is a well-known risk factor for blinding eye diseases. However, epidemiological studies in ophthalmology are still scarce in Japan, and the sizes of the cohorts are relatively small. "Ganka-Ekigaku Network" (GEN, an acronym for the epidemiological network in ophthalmology in Japanese) is established to develop a capacity to boost each epidemiological study and enrich a potential inter-study collaboration to identify risk factors of visual impairment in aged society.

METHODS

We reviewed cohort studies in Japan with the inclusion criteria as: (1) at least n = 1000 at baseline, (2) multiple modalities of ophthalmic data, and (3) diagnosis reviewed by ophthalmologist(s), and (4) ophthalmologists are involved in the investigators group. As of January 2020, GEN includes 4 individual Japanese epidemiological studies namely, Hisayama study, Yamagata Study (Funagata), Tsuruoka Metabolomics Cohort study, and the Nagahama Prospective Genome Cohort for Comprehensive Human Bioscience.

RESULTS

GEN includes approximately 25,000 Japanese participants. The baseline surveys started from 1998 to 2012, and since then the data has been prospectively collected approximately every 5 years. A variety of ophthalmic measurements and other factors have been collected in each study in GEN: ophthalmic measurements (fundus photography, optical coherence tomography, etc.), systemic conditions (laboratory data, etc.), and others (DNA, etc.).

CONCLUSION

GEN is an open platform for observational ophthalmic epidemiological studies to share standardized methodologies. While each study in GEN pursues specific and original research questions, standardization of the methods will enable us to conduct reliable meta-analysis/pooled data analyses.

Myopia Prevalence and Ocular Biometry Features in a General Japanese Population: The Nagahama Study

PURPOSE

To describe the distribution of ocular biometry and refraction in Japanese adults.

DESIGN

Cross-sectional analysis of a prospective cohort study.

PARTICIPANTS

A total of 9850 individuals participated in the first follow-up of the Nagahama Prospective Cohort for Comprehensive Human Bioscience (the Nagahama Study) conducted between 2013 and 2016. Participants were between 34 and 80 years of age.

METHODS

All participants underwent axial length (AL; in millimeters), anterior chamber depth (ACD; in millimeters), corneal diameter (white to white; in millimeters), and central corneal thickness (CCT; in micrometers) measurement (IOL Master; Carl Zeiss Meditec, Dublin, CA) and refraction (spherical equivalent [SE]; in diopters [D]) and corneal curvature (CC; in millimeters) measurement (ARK-530A; Nidek, Aichi, Japan). Distribution of these ocular biometric parameters and prevalence of myopia, high myopia, and extreme myopia were summarized.

MAIN OUTCOME MEASURES

Distribution of ocular biometry and refraction.

RESULTS

After standardization to the national population of 2015, estimates of mean AL and SE were 24.21 mm and -1.44 D, respectively. Estimates of mean CC, corneal diameter, CCT, and ACD were 7.69 mm, 12.01 mm, 543.96 μm, and 3.21 mm, respectively. After standardization of age and gender, the prevalence of myopia (SE, ≤-0.5 D) and high myopia (SE, ≤-6.0 D) were 49.97% and 7.89%, respectively. Approximately 70% of the younger participants (34-59 years of age) showed myopia, whereas high myopia was observed in approximately 10%. Although the number of individuals with myopia or high myopia was higher in the younger age groups, the prevalence of more extreme phenotypes remained stable across all ages, especially in women. Axial length of more than 30 mm was observed only in older women (n = 5 [0.05%]).

CONCLUSIONS

We showed detailed distributions of various ocular biometry and refraction parameters using a large general Japanese cohort. Prevalences of myopia and high myopia from 2013 through 2016 were higher than those in earlier studies, which reflects recent environmental change. However, constant prevalence of extreme myopia across all ages suggests high genetic predisposition of the extreme phenotype.

Genetic associations of myopia severities and endophenotypes in children

OBJECTIVE

To investigate the associations of multiple single-nucleotide polymorphisms (SNPs) with the severities and endophenotypes of myopia in children.

METHODS

A total of 3300 children aged 5-10 years were recruited: 137 moderate and high myopia (SE≤-3.0D), 670 mild myopia (-3.0D<SE≤-0.5D) and 2493 controls (SE>-0.5D). 13 SNPs in 13 genes/loci were selected for genotyping in all subjects using TaqMan assays. Associations between each SNP with myopia severities and ocular traits (spherical equivalent (SE), axial length (AL) and corneal radius (CR)) were analysed.

RESULTS

When compared with controls, SNPs ZC3H11B rs4373767 (OR=1.15, p=0.038), BICC1 rs7084402 (OR=1.18, p=0.005) and GJD2 rs524952 (OR=1.14, p=0.025) showed nominal associations with overall myopia. ZC3H11B rs4373767 and BICC1 rs7084402 showed stronger associations with moderate and high myopia (rs4373767: OR=1.42, p=0.018; rs7084402: OR=1.33, p=0.025), while GJD2 rs524952 had a stronger association with mild myopia (OR=1.14, p=0.025). GJD2 rs524952 also showed a difference between emmetropia and hyperopia (p=0.018). In quantitative trait locus analysis, ZC3H11B rs4373767, KCNQ5 rs7744813 and GJD2 rs524952 were correlated with both myopic SE (β=-0.09, p=0.03; β=-0.12, p=0.007; β=-0.13, p=0.0006, respectively) and AL (β=0.07, p=0.002; β=0.09, p=0.0008; β=0.07, p=0.0003, respectively). SNTB1 rs7839488 was correlated with both AL (β=0.07, p=0.005) and CR (β=0.02, p=0.006). Moreover, ZC3H11B rs4373767-T (β=0.006; p=0.018), KCNQ5 rs7744813-A (β=0.007; p=0.015) and GJD2 rs524952-T (β=0.009; p=0.0006) were correlated with AL-CR ratio.

CONCLUSIONS AND RELEVANCE

ZC3H11B and BICC1 are genetic risk factors for moderate and high myopia, while ZC3H11B, KCNQ5, SNTB1 and GJD2 confer risk to excessive AL in children.

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.

Meta-analysis of 542,934 subjects of European ancestry identifies new genes and mechanisms predisposing to refractive error and myopia

Refractive errors, in particular myopia, are a leading cause of morbidity and disability worldwide. Genetic investigation can improve understanding of the molecular mechanisms that underlie abnormal eye development and impaired vision. We conducted a meta-analysis of genome-wide association studies (GWAS) that involved 542,934 European participants and identified 336 novel genetic loci associated with refractive error. Collectively, all associated genetic variants explain 18.4% of heritability and improve the accuracy of myopia prediction (area under the curve (AUC) = 0.75). Our results suggest that refractive error is genetically heterogeneous, driven by genes that participate in the development of every anatomical component of the eye. In addition, our analyses suggest that genetic factors controlling circadian rhythm and pigmentation are also involved in the development of myopia and refractive error. These results may enable the prediction of refractive error and the development of personalized myopia prevention strategies in the future.

Genome-Wide Association of Genetic Variants With Refraction, Axial Length, and Corneal Curvature: A Longitudinal Study of Chinese Schoolchildren

Background

Myopia is a common eye disorder that is approaching epidemic proportions worldwide. A genome-wide association study identified AREG (rs12511037), GABRR1 (rs13215566), and PDE10A (rs12206610) as being associated with refractive error in Asian populations. The present study investigated the associations between these three genetic variants and the occurrence and development of myopia, spherical equivalent refraction (SER), axial length (AL), and corneal curvature (CC) in a cohort of southeastern Chinese schoolchildren.

Methods

We examined and followed 550 children in grade 1 enrolled in the Wenzhou Epidemiology of Refractive Error (WERE) project. During the 4-year follow-up, non-cycloplegic refraction was evaluated twice each year, and the AL and CC were measured once every year. Age, sex, and the amounts of time spent on near work and outdoors were documented with a questionnaire. Sanger DNA sequencing was used to genotype single nucleotide polymorphisms (SNPs). SNPtest software was used to identify potential genetic variants associated with myopia, SER, AL, and CC. Ten thousand permutations were used to correct for multiple testing.

Results

In total, 469 children, including 249 (53.1%) boys and 220 (46.9%) girls, were included in analyses. The mean age of all the children was 6.33 ± 0.48 years. After adjusting for age, sex, time spent on near work and time spent outdoors, neither the genotypes nor the allele frequencies of the three SNPs were significantly associated with myopic shift, incident myopia or the change in SER. After adjusting for age, sex, near-work time and outdoor time with 10,000 permutations, the genotype AREG (rs12511037) was associated with an increase in AL (P′-values for the dominant, recessive, additive and general models were 0.0032, 0.0275, 0.0045, and 0.0099, respectively); the genotype PDE10A (rs12206610) was associated with a change in CC in the additive (P′ = 0.0096), dominant (P′ = 0.0096), and heterozygous models (P′ = 0.0096).

Conclusion

These findings preliminarily indicate that AREG SNP rs12511037 and PDE10A SNP rs12206610 are etiologically relevant for ocular traits, providing a basis for further exploration of the development of myopia and its molecular mechanism. However, elucidating the role of AREG and PDE10A in the pathogenesis of myopia requires further animal model and human genetic epidemiology studies. This trial is registered as ChiCTR1900020584 at www.Chictr.org.cn.

Genome-wide association meta-analysis of corneal curvature identifies novel loci and shared genetic influences across axial length and refractive error

Corneal curvature, a highly heritable trait, is a key clinical endophenotype for myopia - a major cause of visual impairment and blindness in the world. Here we present a trans-ethnic meta-analysis of corneal curvature GWAS in 44,042 individuals of Caucasian and Asian with replication in 88,218 UK Biobank data. We identified 47 loci (of which 26 are novel), with population-specific signals as well as shared signals across ethnicities. Some identified variants showed precise scaling in corneal curvature and eye elongation (i.e. axial length) to maintain eyes in emmetropia (i.e. HDAC11/FBLN2 rs2630445, RBP3 rs11204213); others exhibited association with myopia with little pleiotropic effects on eye elongation. Implicated genes are involved in extracellular matrix organization, developmental process for body and eye, connective tissue cartilage and glycosylation protein activities. Our study provides insights into population-specific novel genes for corneal curvature, and their pleiotropic effect in regulating eye size or conferring susceptibility to myopia.

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.

Association of the ZC3H11B, ZFHX1B and SNTB1 genes with myopia of different severities

OBJECTIVE

To investigate the associations of single-nucleotide polymorphisms (SNPs) in the ZC3H11B, ZFHX1B, VIPR2, SNTB1 and MIPEP genes with severities of myopia in Chinese populations.

METHODS

Based on previous myopia genome-wide association studies, five SNPs (ZC3H11B rs4373767, ZFHX1B rs13382811, VIPR2 rs2730260, SNTB1 rs7839488 and MIPEP rs9318086) were selected for genotyping in a Chinese cohort of 2079 subjects: 252 extreme myopia, 277 high myopia, 393 moderate myopia, 366 mild myopia and 791 non-myopic controls. Genotyping was performed by TaqMan assays. Allelic frequencies of the SNPs were compared with myopia severities and ophthalmic biometric measurements.

RESULTS

The risk allele T of ZC3H11B SNP rs4373767 was significantly associated with high myopia (OR=1.39, p=0.007) and extreme myopia (OR=1.34, p=0.013) when compared with controls, whereas ZFHX1B rs13382811 (allele T, OR=1.33, p=0.018) and SNTB1 rs7839488 (allele G, OR=1.71, p=8.44E-05) were significantly associated with extreme myopia only. In contrast, there was no significant association of these SNPs with moderate or mild myopia. When compared with mild myopia, subjects carrying T allele of rs4373767 had a risk of progressing to high myopia (spherical equivalent ≤-6 dioptres) (OR=1.29, p=0.017). Similarly, the T allele of rs13382811 also imposed a significant risk to high myopia (OR=1.36, p=0.007). In quantitative traits analysis, SNPs rs4373767, rs13382811 and rs7839488 were correlated with axial length and refractive errors.

CONCLUSIONS

We confirmed ZC3H11B as a susceptibility gene for high and extreme myopia, and ZFHX1B and SNTB for extreme myopia in Chinese populations. Instead of myopia onset, these three genes were more likely to impose risks of progressing to high and extreme myopia.

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.

Family-Based Genome-Wide Association Study of South Indian Pedigrees Supports WNT7B as a Central Corneal Thickness Locus

Purpose

To identify genetic risk factors contributing to central corneal thickness (CCT) in individuals from South India, a population with a high prevalence of ocular disorders.

Methods

One hundred ninety-five individuals from 15 large South Indian pedigrees were genotyped using the Omni2.5 bead array. Family-based association for CCT was conducted using the score test in MERLIN.

Results

Genome-wide association study (GWAS) identified strongest association for single nucleotide polymorphisms (SNPs) in the first intron of WNT7B and CCT (top SNP rs9330813; β = −0.57, 95% confidence interval [CI]: −0.78 to −0.36; P = 1.7 × 10⁻⁷). We further investigated rs9330813 in a Latino cohort and four independent European cohorts. A meta-analysis of these data sets demonstrated statistically significant association between rs9330813 and CCT (β = −3.94, 95% CI: −5.23 to −2.66; P = 1.7 × 10⁻⁹). WNT7B SNPs located in the same genomic region that includes rs9330813 have previously been associated with CCT in Latinos but with other ocular quantitative traits related to myopia (corneal curvature and axial length) in a Japanese population (rs10453441 and rs200329677). To evaluate the specificity of the observed WNT7B association with CCT in the South Indian families, we completed an ocular phenome-wide association study (PheWAS) for the top WNT7B SNPs using 45 ocular traits measured in these same families including corneal curvature and axial length. The ocular PheWAS results indicate that in the South Indian families WNT7B SNPs are primarily associated with CCT.

Conclusions

The results indicate robust evidence for association between WNT7B SNPs and CCT in South Indian pedigrees, and suggest that WNT7B SNPs can have population-specific effects on ocular quantitative traits.

HOXA9 is a novel myopia risk gene

Purpose

A recent meta-analysis revealed PAX6 as a risk gene for myopia. There is a link between PAX6 and HOXA9. Furthermore, HOXA9 has been reported to activate TGF-β that is a risk factor for myopia. We speculate HOXA9 may participate in myopia development.

Methods

The Singapore GUSTO birth cohort provides data on children’s cycloplegic refraction measured at age of 3 years and their methylation profile based on the umbilical cord DNA. The HOXA9 expression levels were measured in the eyes of mono-ocular form deprivation myopia in mice. The plasmid with the mouse HOXA9 cDNA was constructed and then transfected to mouse primary retinal pigment epithelial (RPE) cells. The expression levels of myopia-related genes and cell proliferation were measured in the HOXA9-overexpressed RPE cells.

Results

A total of 519 children had data on methylation profile and cycloplegic refraction. The mean spherical equivalent refraction (SE) was 0.90D. Among 8 SE outliers (worse than -2D), 7 children had HOXA9 hypomethylation. The HOXA9 levels in the retina of myopic eyes was 2.65-fold (p = 0.029; paired t-test) higher than the uncovered fellow eyes. When HOXA9 was over-expressed in the RPE cells, TGF-β, MMP2, FGF2 and IGF1R expression levels were dose-dependently increased by HOXA9. However, over-expression of HOXA9 had no significant influence on IGF1 or HGF expression. In addition, HOXA9 also increased RPE proliferation.

Conclusion

Based on the human, animal and cellular data, the transcription factor HOXA9 may promote the expression of pro-myopia genes and RPE proliferation, which eventually contribute to myopia development.

Genome-wide DNA hypermethylation and homocysteine increase a risk for myopia

AIM

To test for the association between genome-wide methylation and myopia in human and mice.

METHODS

Long interspersed nucleotide element 1 (LINE-1) methylation levels were used to surrogate genome-wide methylation level. We first tested for the association between high myopia (<-6 D) and LINE-1 methylation in leukocytes in 220 cases and 220 control subjects. Secondly, we validated the results of LINE-1 methylation in eyes from the form deprivation myopia (FDM) mice. Furthermore, we calculated the correlation of LINE-1 methylation levels between leukocyte DNA and ocular DNA in the mice. We also tested whether dopamine can alter LINE-1 methylation levels.

RESULTS

The LINE-1 methylation level was significantly higher in the myopic human subjects than controls. The upper and middle tertiles of the methylation levels increased an approximately 2-fold (P≤0.002) risk for myopia than the lower tertile. Similarly, FDM mice had high LINE-1 methylation levels in the leukocyte, retina and sclera, and furthermore the methylation levels detected from these three tissues were significantly correlated. Immunohistochemical staining revealed higher levels of homocysteine and methionine in the rodent myopic eyes than normal eyes. Dopamine treatment to the cells reduced both LINE-1 methylation and DNA methyltransferase levels.

CONCLUSION

LINE-1 hypermethylation may be associated with high myopia in human and mice. Homocysteine and methionine are accumulated in myopic eyes, which may provide excess methyl group for genome-wide methylation.

Molecular genetic aspects of complicated myopia pathogenesis

Complicated myopia (CM) is not only a refractive error but a complex, multifactorial disorder characterized by a mismatch between the optical power of the eye and the axial length that causes the image to be focused off the retina. Genetic factors in progressive myopia play a key role in determining the impact of ecologic factors on refraction development. The majority of genetic variants underlying CM are characterized by modest effect and/or low frequency, which makes them difficult to identify using classic genetic approaches. The genes identified to date account for less than 10% of all myopia cases, suggesting the existence of a large number of yet unidentified low-frequency and/or small-effect variants, which underlie the majority of myopia cases. Genome analysis revealed dozens of loci associated with non-syndromic myopia, and showed that refractive errors are associated with mutations in genes that are involved in the growth and development of the eye by regulating ion transport, neurotransmission, remodeling of extracellular matrix of the retina and other ocular structures. Genetic study of refractive error provides a unique opportunity to detect key molecules that may play important roles in the development of refractive error. Identifying the molecular basis of refractive error helps to understand mechanisms, and subsequently to design rational therapeutic intervention for this condition.

CFH and VIPR2 as susceptibility loci in choroidal thickness and pachychoroid disease central serous chorioretinopathy

Although central serous chorioretinopathy (CSC) presumptively shares pathophysiological basis with age-related macular degeneration (AMD), the CFH risk alleles for AMD are reportedly protective against CSC development. Our finding, that the CFH risk allele for AMD is protective against choroidal thickening in a Japanese cohort, indicates that CFH affects CSC development through its choroid-thickening effects rather than its association with AMD, highlighting the need for a new AMD classification, with CSC/pachychoroid-associated choroidal neovascularization as a distinct disease. Furthermore, our genome-wide association study (GWAS) addressing choroidal thickness successfully discovered a susceptibility gene for CSC: VIPR2. Future GWASs on choroidal thickness will likely discover additional CSC susceptibility genes and provide key molecules to elucidate the pathophysiological difference between CSC and AMD.

Central serous chorioretinopathy (CSC) is a common disease affecting younger people and may lead to vision loss. CSC shares phenotypic overlap with age-related macular degeneration (AMD). As recent studies have revealed a characteristic increase of choroidal thickness in CSC, we conducted a genome-wide association study on choroidal thickness in 3,418 individuals followed by TaqMan assays in 2,692 subjects, and we identified two susceptibility loci: CFH rs800292, an established AMD susceptibility polymorphism, and VIPR2 rs3793217 (P = 2.05 × 10⁻¹⁰ and 6.75 × 10⁻⁸, respectively). Case–control studies using patients with CSC confirmed associations between both polymorphisms and CSC (P = 5.27 × 10⁻⁵ and 5.14 × 10⁻⁵, respectively). The CFH rs800292 G allele is reportedly a risk allele for AMD, whereas the A allele conferred risk for thicker choroid and CSC development. This study not only shows that susceptibility genes for CSC could be discovered using choroidal thickness as a defining variable but also, deepens the understanding of differences between CSC and AMD pathophysiology.

CCDC102B confers risk of low vision and blindness in high myopia

The incidence of high myopia is increasing worldwide with myopic maculopathy, a complication of myopia, often progressing to blindness. Our two-stage genome-wide association study of myopic maculopathy identifies a susceptibility locus at rs11873439 in an intron of CCDC102B (P = 1.77 × 10⁻¹² and P_corr = 1.61 × 10⁻¹⁰). In contrast, this SNP is not significantly associated with myopia itself. The association between rs11873439 and myopic maculopathy is further confirmed in 2317 highly myopic patients (P = 2.40 × 10⁻⁶ and P_corr = 1.72 × 10⁻⁴). CCDC102B is strongly expressed in the retinal pigment epithelium and choroids, where atrophic changes initially occur in myopic maculopathy. The development of myopic maculopathy thus likely exhibits a unique background apart from the development of myopia itself; elucidation of the roles of CCDC102B in myopic maculopathy development may thus provide insights into preventive methods for blindness in patients with high myopia.

Myopic maculopathy is a complication of myopia that often progresses to blindness. Here, in a genome-wide association study, Hosoda et al. find that rs11873439 intronic to CCDC102B is associated with myopic maculopathy, but not with myopia, thus representing a risk factor independent of myopia.

Prevalence of posterior staphyloma and factors associated with its shape in the Japanese population

Myopia is increasing rapidly worldwide. We performed a cross-sectional study to investigate the prevalence of posterior staphyloma, a complication of myopia, and its shape characteristics in relation to age, sex, and axial length (AL) in a Japanese community-based cohort. The right eyes of 3748 participants who underwent fundus photography and optical coherence tomography (OCT) examination were evaluated. Posterior staphyloma prevalence was evaluated using fundus photographs and OCT images. Furthermore, fundus shapes were analyzed by measuring local fundus curvatures on 6 mm cross-line OCT images at intervals of 1 µm. The mean and variance of the curvatures were calculated to represent the fundus shape of each eye for investigation of the relationship between fundus curvature and age, sex, and AL. Seventy-seven eyes (2.05%) had posterior staphyloma. The mean and variance of the fundus curvatures were significantly greater in women than in men and became greater with age, suggesting that the shape of the staphyloma was steeper and less smooth in women and elderly subjects. AL and mean curvature showed a significant correlation (P = 2 × 10⁻¹⁶, R = 0.480), which was significantly affected by age (P < 2 × 10⁻¹⁶). Quantitative analysis of fundus shapes was useful for statistical analysis of posterior staphyloma in relation to age, sex, and AL.

Genome-wide association study identifies WNT7B as a novel locus for central corneal thickness in Latinos

The cornea is the outermost layer of the eye and is a vital component of focusing incoming light on the retina. Central corneal thickness (CCT) is now recognized to have a significant role in ocular health and is a risk factor for various ocular diseases, such as keratoconus and primary open angle glaucoma. Most previous genetic studies utilized European and Asian subjects to identify genetic loci associated with CCT. Minority populations, such as Latinos, may aid in identifying additional loci and improve our understanding of the genetic architecture of CCT. In this study, we conducted a genome-wide association study (GWAS) in Latinos, a traditionally understudied population in genetic research, to further identify loci contributing to CCT. Study participants were genotyped using either the Illumina OmniExpress BeadChip (∼730K markers) or the Illumina Hispanic/SOL BeadChip (∼2.5 million markers). All study participants were 40 years of age and older. We assessed the association between individual single nucleotide polymorphisms (SNPs) and CCT using linear regression, adjusting for age, gender and principal components of genetic ancestry. To expand genomic coverage and to interrogate additional SNPs, we imputed SNPs from the 1000 Genomes Project reference panels. We identified a novel SNP, rs10453441 (P = 6.01E-09), in an intron of WNT7B that is associated with CCT. Furthermore, WNT7B is expressed in the human cornea. We also replicated 11 previously reported loci, including IBTK, RXRA-COL5A1, COL5A1, FOXO1, LRRK1 and ZNF469 (P < 1.25E-3). These findings provide further insight into the genetic architecture of CCT and illustrate that the use of minority groups in GWAS will help identify additional loci.

Novel evidence for complement system activation in chick myopia and hyperopia models: a meta-analysis of transcriptome datasets

Myopia (short-sightedness) and hyperopia (long-sightedness) occur when the eye grows too long or short, respectively, for its refractive power. There are currently approximately 1.45 billion myopes worldwide and prevalence is rising dramatically. Although high myopia significantly increases the risk of developing a range of sight-threatening disorders, the molecular mechanisms underlying ocular growth regulation and its relationship to these secondary complications remain poorly understood. Thus, this study meta-analyzed transcriptome datasets collected in the commonly used chick model of optically-induced refractive error. Fifteen datasets (collected across five previous studies) were obtained from GEO, preprocessed in Bioconductor, and divided into 4 conditions representing early (≤1 day) and late (>1 day) myopia and hyperopia induction. Differentially expressed genes in each condition were then identified using Rank Product meta-analysis. The results provide novel evidence for transcriptional activation of the complement system during both myopia and hyperopia induction, and confirm existing literature implicating cell signaling, mitochondrial, and structural processes in refractive error. Further comparisons demonstrated that the meta-analysis results also significantly improve concordance with broader omics data types (i.e., human genetic association and animal proteomics studies) relative to previous transcriptome studies, and show extensive similarities with the genes linked to age-related macular degeneration, choroidal neovascularization, and cataract.

Altered Refractive Development in Mice With Reduced Levels of Retinal Dopamine

Purpose

The neuromodulator dopamine (DA) has been implicated in the prevention of excessive ocular elongation and myopia in various animal models. This study used retina-specific DA knockout mice to investigate the role of retinal DA in refractive development and susceptibility to experimental myopia.

Methods

Measurements of refractive error, corneal curvature, and ocular biometrics were obtained as a function of age for both untreated and form-deprived (FD) groups of retina-specific tyrosine hydroxylase knockout (rTHKO) and control (Ctrl) mice. Retinas from each group were analyzed by HPLC for levels of DA and its primary metabolite (DOPAC).

Results

Under normal visual conditions, rTHKO mice showed significantly myopic refractions (F(1,188) = 7.602, P < 0.001) and steeper corneas (main effect of genotype F(1,180) = 5.1, P < 0.01) at 4 and 6 weeks of age compared with Ctrl mice. Retina-specific THKO mice also had thinner corneas (main effect of genotype F(1,181) = 37.17, P < 0.001), thinner retinas (F(6,181) = 6.07, P < 0.001), and shorter axial lengths (F(6,181) = 3.78, P < 0.01) than Ctrl mice. Retina-specific THKO retinas contained less than 15% of DA and DOPAC compared with Ctrl retinas, and the remaining DA had a significantly higher turnover, as indicated by DOPAC/DA ratios (Student's t-test, P < 0.05). Retina-specific THKO mice showed similar, yet more variable, responses to 6 weeks of FD compared with Ctrl mice.

Conclusions

Diminished retinal DA induced spontaneous myopia in mice raised under laboratory conditions without form deprivation. The relative myopic shift in rTHKO mice may be explained by steeper corneas, an unexpected finding. The chronic loss of DA did not significantly alter the FD myopia response in rTHKO mice.

Association of IGF1 polymorphism rs6214 with high myopia: A systematic review and meta-analysis

PURPOSE

To conduct a comprehensive evaluation of the association of Insulin-like growth factor 1 (IGF1) polymorphism rs6214 with high myopia through a systematic review and meta-analysis of candidate genetic association study.

METHODS

All case-control association studies on IGF1 and high myopia reported up to 15 June 2016 in PubMed, Embase, Web of Science, and the Chinese Biomedical Database were retrieved. Odds ratios (ORs) and 95% confidence intervals (95% CIs) were calculated for single-nucleotide polymorphism (SNP) using fixed and random effects models according to between study heterogeneity. Publication bias analyses were conducted using Begg's test.

RESULTS

A total of eight studies from published articles were included in our analysis. The meta-analyses for IGF1 rs6214, composed of 4242 high myopia patients and 4430 controls, showed low heterogeneity for the included populations in all the genetic models, except that of the allelic genetic model in the pooled populations. The analyses of all the genetic models in Chinese, Japanese, and overall pooled populations did not identify any significant association between high myopia and IGF1 rs6214.

CONCLUSIONS

This meta-analysis showed there was no association detected between IGF1 rs6214 and high myopia. Given the limited sample size, further investigations including more ethnic groups are required to validate the association.