|Ahead of print publication
A novel G26A variation in 5′ half of TGIF1 gene associates with high myopia in ethnic Kashmiri population from India
Shabhat Rasool1, Rubiya Dar2, Arif Akbar Bhat3, Shiekh Gazalla Ayub1, Muneeb U Rehman4, Sabia Rashid5, Tariq Jan6, Khursheed Iqbal Andrabi2
1 Department of Biotechnology, University of Kashmir, Srinagar, Jammu and Kashmir, India; Department of Biochemistry, Government Medical College, Srinagar, Jammu and Kashmir, India
2 Department of Biotechnology, University of Kashmir, Srinagar, Jammu and Kashmir, India
3 Department of Biochemistry, Government Medical College, Srinagar, Jammu and Kashmir, India
4 Department of Clinical Pharmacy, College of Pharmacy, King Saud University, Riyadh, Saudi Arabia
5 Department of Ophthalmology, Government Medical College, Srinagar, Jammu and Kashmir, India
6 Department of Statistics, University of Kashmir, Srinagar, Jammu and Kashmir, India
|Date of Submission||20-Feb-2019|
|Date of Acceptance||15-Jun-2019|
|Date of Web Publication||17-Oct-2019|
Source of Support: None, Conflict of Interest: None
This study aims to look at novel variations in TGIF1 gene and explores their potential association with high myopia in an ethnic population from Kashmir (India). Genomic DNA was genotyped for polymorphic variations, and allele frequencies were tested for the Hardy–Weinberg disequilibrium in 240 ethnic Kashmiri cases with high myopia with a spherical equivalent of >−6 diopters (D) and compared with emmetropic controls with spherical equivalent within −0.5D in one or both eyes represented by a sample size of 228. In this study, we found a novel sequence variation G26A (GAT to AAT) in 5′ half of TGIF1 gene (p. aspartic acid >asparagine) at a frequency of 62% (148/240, P ≤ 0.0001). Variation appears to associate with high myopia significantly (P ≤ 0.001) as it happens to be present only in high myopia affected individuals. Further, it shows statistical significance for its association with gender and the degree of myopia (P ≤ 0.05). In addition, in silico predictions show that variation likely has an impact on the structure and functional properties of the protein. The assessment of the I-TASSER protein structure showed higher energy for a wild-type protein (−5820.186 kJ/mol) as compared to mutant protein (−6595.593 kJ/mol).
Keywords: Ethnic, gene, myopia, novel, TG1F1, variation
|How to cite this URL:|
Rasool S, Dar R, Bhat AA, Ayub SG, Rehman MU, Rashid S, Jan T, Andrabi KI. A novel G26A variation in 5′ half of TGIF1 gene associates with high myopia in ethnic Kashmiri population from India. Taiwan J Ophthalmol [Epub ahead of print] [cited 2020 Sep 19]. Available from: http://www.e-tjo.org/preprintarticle.asp?id=269486
| Introduction|| |
Myopia, known as nearsightedness, is the most common cause of visual impairment worldwide. It is a condition of the eye where the light that comes in does not directly focus on the retina but in front of it. It is said to be a consequence of mismatch between the power of optical components and the axial length of the eye. This common ocular disorder was earlier thought as a benign refractive error. Nowadays even at low levels, it is associated with high risk for a number of ocular diseases. The high and severe grades of myopia worsen to the stage of “high or pathological myopia.” This condition is often associated with complications, leading to reduced vision or blindness. The ocular refractive components precisely undergo coordinated physical alterations during ocular growth, to attain and maintain normal emmetropic visual acuity, and hence that image focuses directly on the retinal plane. Any discordance between the axial length and other optical refractive components, such as corneal and lenticular curvatures, would result in ametropia and blurred visual acuity.
Myopia is being reported as an epidemic occurring worldwide. Both genetic and environmental components have been associated with the etiology of this potentially blinding ocular condition, as such, making prevention and treatment challenging. Furthermore, studies conducted on humans and animals over the past four decades also provide strong evidence for the involvement of both environmental and genetic factors in the progression of myopia.,,,,
Ethnic diversity plays a great role in the development of myopia reaching as high as 70%–90% in some parts of Asia, 30%–40% in Americans and Europeans, and up to 20% in Africans. It is a highly prevalent and complex disorder involving both genetic and environmental factors or it may be interplay of both genetic and environmental factors. Recent mapping studies have mapped 14 genomic loci to be associated with myopia. MYP2 is a candidate locus of the nonsyndromic autosomal dominant high myopia first identified by Young et al. There are nine known and six hypothetical genes considered to be candidates based on mapped position within the MYP2 interval. Among other genes present in the MYP2 locus is transforming growth factor beta induced factor 1 (TGIF1). TGIF1 is expressed in the sclera, retina, cornea, and optic nerve and competitively inhibits binding of the retinoic acid receptor to a retinoid-responsive promoter. It is possible that mutations in TGIF1 gene may alter its function, and hence, the phase of eye development, thus making it a potential candidate gene to study high myopia.
Two hundred and forty Kashmiri adult participants (from India, above 20 years) with high myopia and ethnically matched 228 healthy adult controls (above 20 years) were enrolled for the study from Government Medical College Hospital (Ophthalmology unit) as well as from our ophthalmologist's clinic. This study was performed with informed consent and following all the guidelines for experimental investigations required by the Institutional Board of Research Studies (BORS), of which all authors are affiliated (No: BORS-BT-20109A-2010). Participants were encouraged to narrate all the details relevant to this study such as age, history of onset of myopia, information regarding close work, and any associated ocular complication. Individuals were excluded if there was a known ocular disease such as retinopathy, cataract or if they had a known genetic disease associated with myopia such as Stickler or Marfan syndrome. The ophthalmic evaluation was done which included measuring visual acquity, keratometry, retinoscopy, slit lamp examination of the anterior segment, fundus examination, and measurement of the axial length.
| Methodology|| |
Genomic DNA was prepared from peripheral blood and subjected to amplification by polymerase chain reaction using following primer sequences: 5'-GGGAATAAGTGAGGGGCTCT-3' (sense) and 5'-CCTGAACCAGTCGCAAAGTT-3' (antisense) generating a 472 bp fragment. Purified PCR products were subjected to denaturation and renaturation procedures for the generation of potential heteroduplexes [Figure 1] and analyzed using conformation-sensitive gel electrophoresis strictly as described by Ganguly et al. Samples that showed unusual mobility during these assays were finally sequenced to confirm the presence of sequence variations along with controls (Scigenom, Cochin). Sequence results obtained in fasta and pdf formats were analyzed using ClustalX version 2 software and by Chromas Pro version 1.49 beta 2 software for the detailed inspection of individual chromatograms allele frequencies were tested for the Hardy–Weinberg disequilibrium. The genotype and allele frequencies were evaluated using the Chi-square test or Fisher's exact test. Odds ratio and confidence interval were also calculated. The amino acid sequence of the protein in fasta format obtained from (NCBI) (www.ncbi.nlm.nih.gov) was submitted to an automated server (I-TASSER) (zhang.bioinformatics.ku.edu/I-TASSER) for three-dimensional (3D) structure prediction. The server furnishes predicted 3D structure in pdb format. Swiss-PdbViewer was used for viewing pdb files and computing the free energy of the predicted 3D structures. Swiss-PdbViewer includes a version of the GROMOS43B1 force field to allow the energy evaluation of a protein structure (macromolecular) as well as repair the distorted geometry through energy minimization. In this implementation, all computations are done in vacuo without reaction field.
|Figure 1: Representative conformation-sensitive gel electrophoresis gel showing heteroduplex bands (as indicated by arrow) in all lanes suggestive of the presence of variation|
Click here to view
| Results|| |
In this study, finally, DNA sequencing was used to confirm the results of heteroduplex assays. In addition to previously reported variations that we have already published the data claiming the association of gene with the disease, a novel missense sequence variation G > A [Figure 2] at codons 26 was identified in 5' half of TGIF1 gene (variant-003, ensemble). An interesting finding of this study was novelty of the variation. The sequence variation G26A (p. aspartic acid >asparagine) which has not been reported till date was present at a frequency of 62% (148/240).
|Figure 2: Representative partial chromatogram of (a) normal sample (b) affected sample showing G > A variation at codon 26 in TGIF1 as indicated by arrow|
Click here to view
A subtle and statistically significant [P ≤ 0.001; [Table 1] difference in the allelic frequency for this sequence variation was indicative of its possible association with high myopia. Furthermore, it could be associated with gender and the degree of myopia [P = 0.01 and P ≤ 0.0001, [Table 2], with the frequency of GA genotype significantly higher in females with degree of myopia >−6 D.
|Table 1: Genotype and allele frequencies of transforming growth factor beta-induced factor 1 gene variation in cases and controls|
Click here to view
|Table 2: Association of transforming growth factor beta-induced factor 1 gene variations with clinical variables in high myopia affected patients|
Click here to view
Insilco prediction results show that calculated energy for the wild-type protein is more (−5820.186 kJ/mol) as compared to the mutant protein (−6595.593 kJ/mol). This change in energy of mutant protein is suggestive of affecting the protein tertiary structure which may, in turn, have some impact on protein function. Therefore, further studies are needed to elucidate the actual role of this mutation on protein structure and function.
| Discussion|| |
Diverse populations have presented inconsistent profile of association data owing largely to heterogeneous nature of the subject populations. Genetic polymorphisms have widely been in use to test the association of a gene with a commonly seen and multifactorial disease instead of single gene disease. Since ethnic differences do exit, it is imperative to substantiate or dispute the relevance of such polymorphism in genetically purer cohorts.
Scavello et al. reported negative association of TGIF1 gene with high myopia, although six polymorphic variations were reported to be associated with high myopia in a Chinese population. However, a second Chinese study excludes the association of these variations with high myopia. A Japanese study group also failed to associate TGIF1 gene with high myopia. Here, it is evident that discrepancy in the association of TGIF1 gene with high myopia might be due to population heterogeneity. In light of these discrepancies, our results can be of immense value as they substantiate the claims of association of TGIF1 gene with high myopia.
Investigating the genetics of common and complex disorders such as myopia remains one of the great challenges in human genetics. Myopia is considered to be a complex and multigenic condition involving several overlapping signaling pathways, each one mediated by a group of distinct genetic profiles. Therefore, studying the genetic polymorphisms of myopia-related genes can further clarify the relationship between genetics and myopia. The association between myopia and various genetic markers has helped increase our knowledge of prevention and treatment of myopia.
Kashmiri population being a pure ethnic group provides an ideal scenario to substantiate the contribution of TGIF1 (if any) in the development of high myopia. Two hundred and forty high myopic and 228 normal controls of Kashmiri ethnicity were recruited for TGIF1 polymorphic studies. Two genotypes GG and GA for codon 26 sequence variations occurred at a frequency of 100%: 0.00% in the control group versus 38%:74% in the high myopia group. A subtle and statistically significant [P ≤ 0.001; [Table 1] difference in the allelic frequency was indicative of its possible association with high myopia. Furthermore, variation was significantly associated with gender and the degree of myopia [P = 0.01 and P ≤ 0.0001, [Table 2] with the frequency of GA genotype significantly higher in females with degree of myopia >−6 D. The calculated energy for the wild-type protein is more (−5820.186 kJ/mol) as compared to the mutant protein (−6595.593 kJ/mol). This change in energy of mutant protein is suggestive of affecting the protein tertiary structure which may, in turn, have some impact on protein function. Therefore, further studies are needed to elucidate the actual role of this mutation on protein structure and function.
In addition, the sequence variation is present in coding sequence of the gene affecting the physicochemical properties of the protein causing change from polar negatively charged aspartic acid to polar and neutral amino acid asparagine.
| Conclusion|| |
Focused investigation is needed to establish the precise role played by TGIF1 in the high myopia development, especially in the context of the above-observed results. Therefore, further studies are needed to rule out the actual role of this novel variation on the protein structure and function.
This study was approved by Board of Research Studies, Departments of Biotechnology, University of Kashmir. Approval No. BORS-BT-20109A-2010.
Financial support and sponsorship
Conflicts of interest
The authors declare that there are no conflicts of interests of this paper.
| References|| |
Wensor M, McCarty CA, Taylor HR. Prevalence and risk factors of Myopia in Victoria, Australia. Arch Ophthalmol 1999;117:658-63.
Flitcroft DI. The complex interactions of retinal, optical and environmental factors in myopia aetiology. Prog Retin Eye Res 2012;31:622-60.
Hosoda Y, Yoshikawa M, Miyake M, Tabara Y, Shimada N, Zhao W, et al.
CCDC102B confers risk of low vision and blindness in high myopia. Nat Commun 2018;9:1782.
Wildsoet CF. Active emmetropization – Evidence for its existence and ramifications for clinical practice. Ophthalmic Physiol Opt 1997;17:279-90.
Zhou J, Rappaport EF, Tobias JW, Young TL. Differential gene expression in mouse sclera during ocular development. Invest Ophthalmol Vis Sci 2006;47:1794-802.
Holden BA. The Charles F. Prentice award lecture 2014: A 50-year research journey: Giants and great collaborators. Optom Vis Sci 2015;92:741-9.
Tkatchenko AV, Tkatchenko TV, Guggenheim JA, Verhoeven VJ, Hysi PG, Wojciechowski R, et al.
APLP2 regulates refractive error and myopia development in mice and humans. PLoS Genet 2015;11:e1005432.
Morgan IG. The biological basis of myopic refractive error. Clin Exp Optom 2003;86:276-88.
Young TL. Molecular genetics of human myopia: An update. Optom Vis Sci 2009;86:E8-22.
Baird PN, Schäche M, Dirani M. The Genes in myopia (GEM) study in understanding the aetiology of refractive errors. Prog Retin Eye Res 2010;29:520-42.
Wojciechowski R. Nature and nurture: The complex genetics of myopia and refractive error. Clin Genet 2011;79:301-20.
Cooper J, Tkatchenko AV. A review of current concepts of the etiology and treatment of myopia. Eye Contact Lens 2018;44:231-47.
Lin HJ, Wan L, Tsai Y, Tsai YY, Fan SS, Tsai CH, et al.
The TGFbeta1 gene codon 10 polymorphism contributes to the genetic predisposition to high myopia. Mol Vis 2006;12:698-703.
Ibay G, Doan B, Reider L, Dana D, Schlifka M, Hu H, et al.
Candidate high myopia loci on chromosomes 18p and 12q do not play a major role in susceptibility to common myopia. BMC Med Genet 2004;5:20.
Young TL, Ronan SM, Drahozal LA, Wildenberg SC, Alvear AB, Oetting WS, et al.
Evidence that a locus for familial high myopia maps to chromosome 18p. Am J Hum Genet 1998;63:109-19.
Scavello GS Jr., Paluru PC, Zhou J, White PS, Rappaport EF, Young TL. Genomic structure and organization of the high grade myopia-2 locus (MYP2) critical region: Mutation screening of 9 positional candidate genes. Mol Vis 2005;11:97-110.
Ganguly A, Rock MJ, Prockop DJ. Conformation-sensitive gel electrophoresis for rapid detection of single-base differences in double-stranded PCR products and DNA fragments: Evidence for solvent-induced bends in DNA heteroduplexes. Proc Natl Acad Sci U S A 1993;90:10325-9.
Thompson JD, Gibson TJ, Plewniak F, Jeanmougin F, Higgins DG. The CLUSTAL_X windows interface: Flexible strategies for multiple sequence alignment aided by quality analysis tools. Nucleic Acids Res 1997;25:4876-82.
Zhang Y. Template-based modeling and free modeling by I-TASSER in CASP7. Proteins 2007;69 Suppl 8:108-17.
Camacho CJ, Gatchell DW, Kimura SR, Vajda S. Scoring docked conformations generated by rigid-body protein-protein docking. Proteins 2000;40:525-37.
Sippl MJ. Calculation of conformational ensembles from potentials of mean force. An approach to the knowledge-based prediction of local structures in globular proteins. J Mol Biol 1990;213:859-83.
Ahmed I, Rasool S, Jan T, Qureshi T, Naykoo NA, Andrabi KI. TGIF1 is a potential candidate gene for high myopia in ethnic Kashmiri population. Curr Eye Res 2014;39:282-90.
Lam DS, Lee WS, Leung YF, Tam PO, Fan DS, Fan BJ, et al.
TGFbeta-induced factor: A candidate gene for high myopia. Invest Ophthalmol Vis Sci 2003;44:1012-5.
Li J, Zhang QJ, Xiao XS, Li JZ, Zhang FS, Li SQ, et al.
The SNPs analysis of encoding sequence of interacting factor gene in Chinese population. Zhonghua Yi Xue Yi Chuan Xue Za Zhi 2003;20:454-6.
Hasumi Y, Inoko H, Mano S, Ota M, Okada E, Kulski JK, et al.
Analysis of single nucleotide polymorphisms at 13 loci within the transforming growth factor-induced factor gene shows no association with high myopia in Japanese subjects. Immunogenetics 2006;58:947-53.
Honda S, Fujii S, Sekiya Y, Yamamoto M. Retinal control on the axial length mediated by transforming growth factor-beta in chick eye. Invest Ophthalmol Vis Sci 1996;37:2519-26.
[Figure 1], [Figure 2]
[Table 1], [Table 2]