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

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.

Interactions between genetic variants and near-work activities in incident myopia in schoolchildren: a 4-year prospective longitudinal study

CLINICAL RELEVANCE

Knowledge of interactions between genetic variants and near-work activities at the onset of myopia can facilitate health education regarding myopia.

BACKGROUND

To investigate the interactions between genetic variants (PDE10A, AREG and GABRR1) and near-work activities in the onset of myopia in southeastern Chinese school children.

METHODS

A total of 458 non-myopic, grade 1 children aged 6-7 years were included in a 4-year follow-up examination; 409 children were assessed further. Manifest (non-cycloplegic) refraction and axial length (AL) were measured every year, and questionnaires were administered annually to assess information regarding the demographic characteristics of children, near-work activities, outdoor exposure and parental myopia. Oral mucosa was collected in the last year of follow-up, and Sanger sequencing was used to genotype single nucleotide polymorphisms (SNPs) in DNA.

RESULTS

The cumulative change in the spherical equivalent refraction (SER) over 4 years was -1.20 ± 1.00 D, and the proportion of children with incident myopia was 42.9%. Multivariate logistic regression analysis showed that an increased amount of time spent doing homework (>2 h/d) was an independent risk factor for incident myopia. The PDE10A rs12206610CT genotype and spending > 5 h/d on near-work activities showed an interaction for incident myopia (OR = 4.29, 95% CI: 1.27-14.53; P_interaction = 0.02); moreover, the rs12206610CT genotype carriers who used electronic devices for > 1 h/d displayed an increased risk of incident myopia (OR = 3.43, 95% CI: 1.07-11.01; P_interaction = 0.043).

CONCLUSIONS

The rs2206610CT genotype carriers with near-work activities of >5 h/d were more likely to show incident myopia, especially those who used electronic devices >1 h/d. However, interactions between the rs12206610 SNP and near-work activities require further verification in animal models and larger sample cohorts.

Gene-environment Interaction in Spherical Equivalent and Myopia: An Evidence-based Review

PURPOSE

Association between gene-environment interaction and myopia/spherical equivalent has not been systematically reported. This paper reviewed nine studies concerning gene-environment interaction in myopia.

METHODS

We obtained relevant studies concerning gene-environment interaction in myopia by systematically searching the MEDLINE(PubMed), Cochrane, Web of Science, CNKI, Wanfang databases before 31 March 2020. Data were analyzed by STATA version 16.0 software, and figures were drawn by ArcGIS V.10.0 software.

RESULTS

Nine studies were included in this review concerning gene-environment interaction. Gene and education interaction in adult cohorts suggested a more significant genetic effect in higher education levels than lower education levels, using both candidate genes and PRS approaches. Several interacted genetic variants, including ZMAT4(rs2137277), GJD2(rs524952), TJP2 (rs11145488) from adult study and ZMAT4(rs7829127) from child study are pinpointed out, but the replication attempts were limited. Besides, the genetic effect was associated with a significant shift at a higher educational level (Pooled β = -0.15,95%CI = -0.19-0.11) towards myopia than that at a lower education level (Pooled β = -0.10,95%CI = -0.11-0.09).

CONCLUSION

This study summarizes the relationship between gene-environment interaction and myopia, and interaction effect of the gene or genetic risk score with the environment could be found in these studies. The effect of gene-environment (higher education) interaction substantially impacts myopia in adult studies. Evidence that environmental factors (Increased near-work time/decreased outdoor activities) increase the genetic risk is still limited, and specific SNPs contributing to gene-environment effect are not determined yet.

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.