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.