Molecular
Vision 2012; 18:937-944
<http://www.molvis.org/molvis/v18/a97>
Received 10 January 2012 | Accepted 9 April 2012 | Published 12
April 2012
Khanh-Nhat Tran-Viet,1 Elizabeth St.Germain,1 Vincent Soler,2 Caldwell Powell,1 Sing-Hui Lim,1,4 Thomas Klemm,4 Seang Mei Saw,3,4,5 Terri L. Young1,4,6
1The Center for Human Genetics, Duke University Medical Center, Durham, NC; 2UMRS 563, Centre de Physiopathologie de Toulouse Purpan, Université Paul Sabatier, Toulouse, France; 3Saw Swee Hock School of Public Health, National University of Singapore, Singapore; 4Duke-National University of Singapore Graduate Medical School, Singapore; 5Singapore Eye Research Institute, Singapore; 6The Department of Ophthalmology, Duke University Eye Center, Durham, NC
Correspondence to: Terri L. Young, M.D., Professor of Ophthalmology, Pediatrics, and Medicine, Duke University Medical Center, The Center for Human Genetics, 905 S. LaSalle Street, Durham, NC, 27710; Phone: (919) 684-1808; FAX: (919) 684-0906; email: terri.young@duke.edu
Purpose: Myopia, or nearsightedness, is highly prevalent in Asian countries and is considered a serious public health issue globally. High-grade myopia can predispose individuals to myopic maculopathy, premature cataracts, retinal detachment, and glaucoma. A recent study implicated zinc finger protein 644 isoform 1 (ZNF644) variants with non-syndromic high-grade myopia in a Chinese-Asian population. Herein we focused on investigating the role for ZNF644 variants in high-grade myopia in a United States (US) cohort.
Methods: DNA from a case cohort of 131 subject participants diagnosed with high-grade myopia was screened for ZNF644 variants. Spherical refractive error of -≤-6.00 diopters (D) in at least one eye was defined as affected. All coding, intron/exon boundaries were screened using Sanger sequencing. Single nucleotide allele frequencies were determined by screening 672 ethnically matched controls.
Results: Sequencing analysis did not detect previously reported mutations. However, our analysis identified 2 novel single nucleotide variants (c.725C>T, c.821A>T) in 2 high-grade myopia individuals- one Caucasian and one African American, respectively. These variants were not found in normal controls. A rare variant - dbsSNP132 (rs12117237→c.2119A>G) - with a minor allele frequency of 0.2% was present in 6 additional cases, but was also present in 5 controls.
Conclusions: Our study has identified two novel variants in ZNF644 associated with high-grade myopia in a US cohort. Our results suggest that ZNF644 may play a role in myopia development.
Myopia is a common ocular disorder resulting from excessive axial elongation of the globe [1,2]. High myopia in particular can predispose one to several ocular complications such as myopic maculopathy, premature cataracts, retinal detachment, and glaucoma, and is considered a public health concern in numerous countries around the world [3]. Myopia prevalence rates vary world-wide. A recent study by Vitale et al. [4] reported that 33.1% of the USA population greater than the age of 20 has some degree of myopia. Studies have shown that prevalence rates in Asia countries are higher. For example, 84% of school children in Taiwan develop myopia by the age of 18 [5,6]. Additionally, more than 85% of Hong Kong Chinese school children aged between 13 and 15 years are myopic [7]. The economic impact of refractive error management can be substantial- vision impairment correction costs account for 3.8 to 7.2 billion dollars annually in the USA alone [8].
Accepted as a common complex disorder, myopia is thought to be influenced by both genetic and environmental factors [9,10]. To date, over 20 genetic loci have been mapped using linkage analysis for both common myopia and high myopia [11-17]. A recent plethora of genome wide association studies have shown positive associations with refractive error [14,18-25]. Of late, large parallel sequencing techniques have been used to identify causal genes for ocular disorders including myopia. Rare ocular diseases such as retinitis pigmentosa and familial exudative vitreoretinopathy were among the first where causative genes were successfully identified using exome sequencing [26-28]. More recently, Mordechai et al. [29] identified a leprecan-like 1 (LEPREL1) mutation as causal for autosomal recessive high-grade axial myopia in a large, consanguinous Bedouin Israeli kindred. The mutation identified in their study demonstrated an autosomal recessive mode of inheritance with variable expressivity. As the lone finding for high myopia and autosomal dominant inheritance, Shi et al. [30] recently used exome sequencing to identify mutations in zinc finger protein 644 isoform 1 (ZNF644) – first in a large pedigree with autosomal dominant high myopia, and then replicated in a Chinese cohort. ZNF644 is a transcription factor gene, and is expressed in the retina and retinal pigment epithelium (RPE) [30]. To our knowledge, the role of ZNF644 has not been studied in a myopic USA cohort composed primarily of Caucasians. We therefore screened for ZNF644 mutations in a high-grade myopia USA population data set. The identification of mutations in ZNF644 in other ethnicities supplements existing knowledge of the etiology and heterogeneity of this debilitating eye disorder.
Informed consent was obtained from all participants before entering the study, with approval by the Institutional Review Board according to the principles of the Declaration of Helsinki. High-grade myopia cases were defined as individuals with a spherical refractive error of greater than or equal to −6.00 diopters (D) in at least one eye. From venous blood samples, genomic DNA was extracted using AutoPure LS® DNA Extractor and PUREGENE™ reagents (Gentra Systems Inc., Minneapolis, MN). DNA for 672 ethnically matched Caucasian healthy control participants was purchased commercially (The Centre for Applied Genomics, The Hospital for Sick Children, Toronto, Canada). Although the controls had no documented ocular malformations, refractive error information was not available for these individuals. Additionally, 50 ethnically matched African American participant controls with refractive error data were ascertained, and blood for DNA extraction was collected for this study.
ZNF644 (NM_201269.1) encodes for a zinc finger transcription factor which maps to chromosome 1p22.2 (Chromosome1:91,380,860–91,487,671 – GRCh37.p5). The gene comprises 6 exons, five of which are coding. Polymerase chain reaction (PCR) and sequencing primers were designed to cover all coding and untranslated gene regions (UTR) including intron-exon boundaries using the ExonPrimer program (Helmholtz Center, Munich, Germany). Primers were selected to produce amplification products not to exceed 850 base pairs (bp) in size for optimal sequence output and analysis. A total of thirteen primer sets were designed to ensure full coverage of the exons and the flanking intronic regions (Table 1).
Samples were amplified using standard PCR protocol and amplicons were visualized after agarose gel (2%) electrophoresis. Sequencing of the amplicons was then completed on an Applied Biosystems ABI3730 xlrobotics using BigDye™ Terminator 3.1 technology (Applied Biosystems, Inc. [ABI], Foster City, CA). Sequences were analyzed using the Sequencher 5.0™ program (Gene Codes, Ann Arbor, MI), and were compared against the known reference sequence (GRCh37.p5) and analyzed for sequence variation. Single nucleotide variants (SNVs) that were novel and/or coding non-synonymous with a minor allele frequency (MAF) of less than 1% were checked for co-segregation in remaining family member samples.
Allelic discrimination assays were employed to measure the allelic frequencies in 672 Caucasian matched control DNA samples using the TaqMan® SNP Genotyping system (Applied Biosystems). Assays were designed according to Applied Biosystems specifications using a combination of unlabeled primers and minor groove binding (MGB) probes with fluorescently labeled dyes (FAM and VIC) to interrogate the base pair of interest. Reactions were completed and ABI 7900 robotics was used to read the allelic calls for each control sample (Applied Biosystems). SDS v2.4 Software provided by ABI was used to analyze each sample and accurately analyze each genotype call.
Full ophthalmologic exams were performed on all subject participants with one or more individual(s) with high grade myopia. The average cycloplegic spherical refractive error for 131 high myopia cases was −11.22 D for the right eye (OD) and −11.48 for the left eye (OS), and the range was from −6.00 D to −50.00 D (OD) and −5.25 D to −50.00 D (OS) across all cases. Of the 131 case participants screened for ZNF644, 74% (97/131) were Caucasians, 12.2% (16/131) were African Americans, 10% (13/131) were Asians, while 5 individuals (3.8%) were of Hispanic descent or declined to self-identify.
The entire coding and untranslated DNA sequence of ZNF644 was sequenced in 131 high myopia patients. In all, we identified a total of 31 heterozygous SNVs in ZNF644, of which 10 were missense, 7 silent, 9 untranslated, and 5 intronic (Table 2 and Table 3).
From seven synonymous sequence variants detected in our cohort, we identified 4 novel SNVs in four separate individuals, which were either within exon 3 or 6. We also found 9 SNVs in the UTR, where 6 of 7novel variants were unique to an individual. In addition to the exonic variants, 4 intronic variants were previously reported in the dbSNP132 database, while 1 was novel and unique to an individual. All intronic variants observed had high MAFs, and/or were not located near splice-site junctions (Table 2).
Of the 10 missense variants identified, 4 were novel and 6 were reported in the dbSNP132 database (Table 3). One novel missense variant was identified in Caucasian individual IND0603809. This adenine to thymine substitution at position c.821A>T alters glutamic acid to valine (Glu274Val), and was not found in other high myopic case samples. A second novel missense variant was identified on exon 3 in African American individual IND0603564. This cytosine to thymine substitution (c.725C>T) results in a threonine to methionine (Thr242Met) amino acid change was only present in this individual. Due to the lack of additional family members for these individuals, segregation analysis was not possible (Figure 1). Novel variants c.2116C>T (His706Tyr) and c.3299G>A (Arg1100His) were identified in Caucasian individuals IND0519772 and IND0519764, respectively. Neither variant segregated with disease after screening available family members.
To confirm the rarity of c.821A>T, a TaqMan genotyping assay was employed to screen a DNA sample data set of Caucasian controls. The variant (c.821A>T) was not replicated in 672 Caucasian controls by allelic discrimination genotyping method. The novel variant found in African American individual IND0603564 (c.725C>T) was not replicated in 50 African American control samples via Sanger sequencing (Figure 1).
Four missense variants in dbSNP132 were identified in African American individuals IND0603564 and IND0603416 but were ruled out due to minor allele frequency not meeting criteria below 1% (Table 3). Missense variant rs12117237, the only variant in dbSNP132 with a minor allele frequency of less than 1% in public databases was identified in 6 Caucasian high myopia cases (0.2%). The minor allele was also seen in 0.75% (5/672) of controls. Fisher’s exact two-sided test for rs12117237 demonstrated a p-value of 0.0015 when comparing the prevalence rates of the allele between cases and controls (SAS Institute Inc., Cary, NC). Other family members of the cases were sequenced to determine co-segregation. Two cases had uninformative family members or were simplex cases, and thus were not useful for determining segregation of the variant. Of the 6 myopia cases identified with rs12117237, 4 families (MYP19, MYP83, MYP89, and MYP113) had at least 3 additional family members. These members were sequenced, and MYP19 and MYP89 demonstrated phenotype co-segregation with the SNV (Figure 2).
Mutation screening of ZNF644 has successfully identified 2 novel missense variants (c.725C>T, c.821A>T) in two ethnicities within our US population, thus supporting the original report of a causal candidate gene for high-grade myopia (Table 2). The two novel variants from our study (c.725C>T, c.821A>T) localize to a conserved region on the third exon, where Shi et al. [30] previously identified several mutations. The variants were not seen in 1344 and 100 ethnically matched chromosomes, respectively, confirming their rarity. PolyPhen-2 software predicted c.725C>T to be possibly damaging, while c.821A>T was tolerated. With the additional variant information from our study, exon 3 may be a hotspot for susceptibility for high-grade myopia in multiple ethnicities. The clustering of variants in exon 3 may depict the importance in protein domain structures and gene regulatory functions of ZNF644.
The single nucleotide variant rs12117237 was present in 6 Caucasian case samples, but also present in 5 ethnically matched control DNA samples. Public database dbSNP132 reports this variant’s MAF to be 0.2% in a European population, and our p-value of 0.0015 by Fisher’s exact test suggests presence of the alleles likely did not occur by chance. As a result of sequencing additional family members in three high myopia cases for the variant, we discovered 100% co-segregation of the mutation with the disease phenotype in two of four families. PolyPhen-2 software predicts the amino acid change to be tolerated. The variant could be a common rare variant, and proper functional validation would be important to determine whether it is a risk or a causal allele [31]. However, the MAF is in line with prevalence rates for high-grade myopia in a general population and a founder mutation effect may be plausible. Clinical information was missing refractive error from the 5 ethnically matched controls who also presented with the variant in the heterozygous state. The prevalence rate of high myopia is estimated to be 4.5% in the US population [32,33]. Moreover, a recent study suggests that the rarer a variant, the higher the likelihood that the variant functionally alters the protein – thus rs12117237 may be associated with high-grade myopia in a Caucasian population [34]. A larger case-control population data set should be tested to understand the true significance of the association of this variant to high-grade myopia in Caucasians.
ZNF644 is a transcription factor that may play a role in protein domain structures or regulatory functions [30]. Expressed in all tissue types, it follows the trend of ubiquitously expressed genes pathogenic to ocular diseases [35]. To date, it is widely accepted that transcription factor genes in both humans and mouse can play an important role in mammalian eye growth and development. For example, paired-like homeodomain transcription factors 2 and 3 (PITX2, PITX3) have been implicated in Axenfeld-Rieger syndrome and cataracts, respectively [36-38]. Microphthalmia transcription factor (Mitf) is associated with ocular albinism, and paired box 6 (PAX6) has been associated with retinal degeneration, extreme myopia, and corneal innervation [20,39-46]. However, the true function of ZNF644 remains unclear, and molecular characterization of ZNF644 is necessary.
To the best of our knowledge, this is the first successful study confirming a gene implicated with non-syndromic high-grade myopia determined by exome sequencing. We identified two novel variants in ZNF644 in our cohort, in addition to a known variant that demonstrated association. The discovery of previously unidentified variants is not expected due to genetic heterogeneity present in rare complex disease across multiple populations [47]. Ethnic group specific alleles due to founder effect may explain why the previously reported variants were only present in the Asian population whereas the variants discovered in this report appear to be specific to the Caucasian and African American ethnicities. Determining pathogenic rare missense variants remains a challenge for complex diseases. Identification of novel variants in separate ethnicities emphasizes the importance and demonstrates the power of current research approaches.
The authors thank all subjects who participated in this study. This research effort was supported by the National Institutes of Health Grant EY014685, The Lew Wasserman Award from Research To Prevent Blindness Inc., and the Duke- National University of Singapore core grant. (T.L.Y.)