Novel truncation mutations in MYRF cause autosomal dominant high hyperopia mapped to 11p12–q13.3

High hyperopia is a common and severe form of refractive error. Genetic factors play important roles in the development of high hyperopia but the exact gene responsible for this condition is mostly unknown. We identified a large Chinese family with autosomal dominant high hyperopia. A genome-wide linkage scan mapped the high hyperopia to chromosome 11p12–q13.3, with maximum log of the odds scores of 4.68 at theta = 0 for D11S987. Parallel whole-exome sequencing detected a novel c.3377delG (p.Gly1126Valfs*31) heterozygous mutation in the MYRF gene within the linkage interval. Whole-exome sequencing in other 121 probands with high hyperopia identified additional novel mutations in MYRF within two other families: a de novo c.3274_3275delAG (p.Leu1093Profs*22) heterozygous mutation and a c.3194+2T>C heterozygous mutation. All three mutations are located in the C-terminal region of MYRF and are predicted to result in truncation of that portion. Two patients from two of the three families developed angle-closure glaucoma. These three mutations were present in neither the ExAC database nor our in-house whole-exome sequencing data from 3280 individuals. No other truncation mutations in MYRF were detected in the 3280 individuals. Knockdown of myrf resulted in small eye size in zebrafish. These evidence all support that truncation mutations in the C-terminal region of MYRF are responsible for autosomal dominant high hyperopia in these families. Our results may provide useful clues for further understanding the functional role of the C-terminal region of this critical myelin regulatory factor, as well as the molecular pathogenesis of high hyperopia and its associated angle-closure glaucoma. Electronic supplementary material The online version of this article (10.1007/s00439-019-02039-z) contains supplementary material, which is available to authorized users.


Introduction
High hyperopia, defined as cycloplegic sphere refraction ≥ + 5.00 diopters (D) (Association 2006) or ocular axial length < 21 mm (Carricondo et al. 2018;Fuchs et al. 2005), is a common severe refractive error. High hyperopia is present in approximately 0.8-3.0% of children between 3 and 6 years old, according to population studies (Giordano et al. 2009;Multi-Ethnic Pediatric Eye Disease Study 2010). This condition is frequently associated with blurred vision, asthenopia, accommodative and binocular dysfunction, amblyopia, strabismus, or even primary angle-closure glaucoma (Association 2006;Jonas et al. 2017;Klimek et al. 2004;Shen et al. 2016). Most cases with high hyperopia are characterized by a normal corneal diameter and a short axial length, with no other significant ocular or systemic anomalies. High hyperopia is also observed in other ocular or systemic diseases, such as microphthalmia (Gal et al. 2011;Sundin et al. 2005), Leber congenital amaurosis (Abouzeid et al. 2006), retinal dystrophy (Bifari et al. 2016), Xueshan Xiao and Wenmin Sun contributed equally to this work.

Electronic supplementary material
The online version of this article (https ://doi.org/10.1007/s0043 9-019-02039 -z) contains supplementary material, which is available to authorized users. and Alstrom syndrome (Khan et al. 2015), in which microphthalmia accompanies a small cornea, while others have retinal degeneration.
Two other extreme forms of high hyperopia have been frequently described: posterior microphthalmia and nanophthalmos. Both of these conditions feature axial length < 20 mm and extremely high hyperopia (usually ≥ + 8 D) (Altintas et al. 1997;Nowilaty et al. 2013a;Park et al. 2016;Vingolo et al. 1994). Posterior microphthalmia usually features a normal-appearing anterior segment but is frequently associated with typical papillomacular fold, sclerochoroidal thickening, and/or several other related changes to the posterior fundus (Erdol et al. 2008;Goldblum and Mojon 1999;Khairallah et al. 2002;Nowilaty et al. 2013b;Park et al. 2016;Spitznas et al. 1983). Nanophthalmos is characterized by a structurally normal but generally small eye, with the anterior and posterior segments both being affected. Apart from microcornea and a shallow anterior chamber, nanophthalmos may also feature fundus changes similar to those observed in posterior microphthalmia (Helvacioglu et al. 2014;Yalcindag et al. 2011). Biometric and molecular characterization suggests that posterior microphthalmia and nanophthalmos both belong to a spectrum of high hyperopia but are not distinct phenotypes (Nowilaty et al. 2013a). This characterization is further supported by mutations in MFRP in different families with nanophthalmos (Sundin et al. 2005), high hyperopia (Xu et al. 2016), or posterior microphthalmia (Aldahmesh et al. 2011;Matsushita et al. 2012;Wasmann et al. 2014). Similarly, mutations in PRSS56 (Gal et al. 2011;Nair et al. 2011;Orr et al. 2011;Said et al. 2013) have been identified in patients with either posterior microphthalmia or nanophthalmos.
Genetic factors play important roles in the development of high hyperopia. In the clinic, high hyperopia may demonstrate an autosomal dominant or autosomal recessive pattern of inheritance in families, although most cases are sporadic. Families with high hyperopia have been reported previously (Fledelius et al. 2004;Fuchs et al. 2005). Family aggregation studies, twin studies, and genome-wide association studies all support a genetic contribution to hyperopia (Dirani et al. 2006;Hammond et al. 2001;Lee et al. 2001;Simpson et al. 2014;Wojciechowski et al. 2005). Four loci for nanophthalmos, i.e., NNO1 for autosomal dominant (Othman et al. 1998), NNO2 for autosomal recessive (Sundin et al. 2005), NNO3 for autosomal dominant (Li et al. 2008), and NNO4 for autosomal dominant (Awadalla et al. 2014), have been mapped, in which mutations in MFRP are responsible for NNO2 while those in TMEM98 are responsible for NNO4. Apart from MFRP and TMEM98, high hyperopia is a frequent sign in other ocular or systemic diseases associated with mutations in several following genes, including PRSS56 (Gal et al. 2011), CRB1 (Abouzeid et al. 2006), LCA5 (den Hollander et al. 2007), IFT140 (Bifari et al. 2016), and KERA (Khan et al. 2004). However, most patients with high hyperopia do not have mutations in these genes, according to our in-house unpublished data and our previous studies (Jiang et al. 2013;Wang et al. 2009;Xu et al. 2016).
In the current study, genome-wide linkage analysis mapped high hyperopia in a large family to chromosome 11p12-q13.3. Whole-exome sequencing identified three novel truncation mutations in MYRF located in the linkage interval in three families with high hyperopia, including the family mapped to this locus. Sanger sequencing confirmed the mutations and their cosegregation with high hyperopia in the three families. All patients with high hyperopia were free of other systemic anomalies, including cardiac and urogenital anomalies. Two patients from two families developed angle-closure glaucoma. Additional analyses of mutation prediction, mutations in the database, mutation location, and myrf knockdown in zebrafish all support that truncation mutations in the C-terminal region of this critical transcriptional regulator for central nervous system myelination are responsible for high hyperopia that is associated with angle-closure glaucoma.

Families and clinical data
One hundred and twenty-two probands with high hyperopia and their available family members, including a large family with 14 affected individuals across four generations (Family ZOC710536), were recruited from our outpatient clinic at the Zhongshan Ophthalmic Center as part of our ongoing program for the collection of genomic DNA to investigate hereditary eye diseases. Written informed consent was obtained from the participating individuals or their guardians prior to the study, in compliance with the Declaration of Helsinki and the Guidance of Sample Collection of Human Genetic Diseases (863-Plan) by the Ministry of Public Health of China. This study has been approved by the Institutional Review Board of the Zhongshan Ophthalmic Center. The disease in family ZOC710536 was likely transmitted as an autosomal dominant trait (Fig. 1). The diagnosis of high hyperopia was based on the following criteria as previously described (Carricondo et al. 2018;Fuchs et al. 2005;Xu et al. 2016): (1) bilateral cycloplegic sphere refraction ≥ + 5.00 D based on the Optometric Clinical Practice Guidelines for the Care of the Patient with Hyperopia (Association 2006), or ocular axial length < 21 mm (Carricondo et al. 2018;Fuchs et al. 2005); (2) corneal diameter in the normal range (≥ 11 mm) (Carricondo et al. 2018); (3) no systemic hereditary disease or ocular disease, including retinal degeneration or developmental anomalies of the iris, lens, choroid or optic disc. Available subjects usually received routine ophthalmological examinations, including visual acuity testing, slit-lamp examination and direct ophthalmoscopic visualization. Refractive errors were measured using an autorefractometer after cycloplegia, and ocular axial length was measured using an optical biometer, as described in our previous study (Xu et al. 2016).

Mapping and identification of the novel causative gene
For the largest family in this study with high hyperopia (ZOC710536), parallel genome-wide linkage scans, whole-exome sequencing (WES), and wholegenome sequencing of genomic DNA were carried out as described in our previous study . WES alone was also performed on genomic DNA from 121 probands of additional families with high hyperopia. Sanger sequencing was used to confirm potential pathogenic mutations.
Briefly, a genome-wide linkage scan was initially conducted on genomic DNA from 16 of the 19 individuals in family ZOC710536 (excluding IV:1, IV:3, and V:3; Fig. 1) using 400 5′-fluorescently labeled microsatellite markers according to panels 1-28 of the ABI PRISM linkage Mapping Set Version 2 as previously described (Zhang et al. 2005). Two-point linkage analysis was performed using the MLINK program of the FASTLINK implementation of the LINKAGE program package (Lathrop and Lalouel 1984;Schaffer et al. 1994). The high hyperopia in the family was characterized as an autosomal dominant trait with complete penetrance and a disease allele frequency of 0.0001. Fine mapping of a candidate locus was carried out on all 19 individuals in the family using an M13tailed primer PCR method to genotype additional markers (Barkley et al. 2007). Haplotypes were generated using the Cyrillic 2.1 program (Cyrillic Software, Wallingford, Oxfordshire OX10 8BA) and confirmed by inspection. The criteria for establishing linkage have been previously described (Lander and Kruglyak 1995).
Meanwhile, WES was performed on genomic DNA from four affected individuals (III:5, III:9, IV:2, and V:1 in Fig. 1) in family ZOC710536, as well as from 121 probands in other additional families with high hyperopia, with an Agilent SureSelect Human All Exon Enrichment Kit V4 (Santa Clara, CA, USA) (51189318 base pairs) array on an Illumina HiSeq 2000 101PE (San Diego, CA, USA) as described in our previous study . Variants detected by WES were initially filtered by multistep bioinformatics analysis: the stages were as follows: (1) excluding variants in noncoding regions as well as synonymous variants that did not affect splice sites according to the Berkeley Drosophila Genome Project [BDGP; http://www.fruit fly.org/ (in the public domain)]; (2) excluding variants with minor allele frequency (MAF) ≥ 1%; (3) excluding variants not shared by four affected individuals in the largest family, ZOC710536; and (4) searching for unique variants in genes with similar mutations in 121 additional families with high hyperopia. WES data from the Exome Aggregation Consortium (ExAC) database and our in-house WES data from 3280 unrelated probands with different forms of inherited eye diseases served as reference. Sanger dideoxy sequencing was used to confirm potential pathogenic variants and to evaluate their cosegregation in family members. Wholegenome sequencing on two Individuals (IV:7 and IV:9 in Fig. 1) was carried out using a commercial service to test for mutations in the linkage interval that had been missed by WES.

Knockdown of myrf in zebrafish
Wild-type AB zebrafish (Danio rerio) embryos were obtained from Sun Yat-sen University. Zebrafish embryos were fertilized in vitro at 28.5 °C, maintained on a 14-h light:10-h dark cycle, and staged in hours post-fertilization (hpf) or days post-fertilization (dpf). Two morpholino oligonucleotides (MOs) were purchased from Gene Tools, LLC (Corvallis, OR, USA), including an antisense MO targeting translation initiation (ATG) of myrf (myrf MO, 5′-CAA TCA CGG TTT CCA TTC TAG CGG C-3′) and a standard control MO (std MO, 5′-CCT CTT ACC TCA GTT ACA ATT TAT A-3′). Both of the MOs were dissolved in distilled water and microinjected into the egg yolks of zebrafish embryos in the one-cell stage. To test the efficiency of myrf MO, we amplified the complementary sequence of myrf MO using primers crossing the MO sequence, designed in Primer3 (S1 Table) and inserted in-frame before the green fluorescent protein (GFP) coding sequence of the pEGFPN1 vector.
For the rescue experiment, the full-length coding sequence of myrf in zebrafish was amplified by PCR using cDNA from wild-type zebrafish at 72 hpf as the template and cloned into the pCS2+ vector (primers are listed in S1 Table). The myrf mRNA was synthesized using the mMES-SAGE mMACHINE™ SP6 kit (Invitrogen AM1340) and pCS2-myrf plasmid as a template. The mRNA was diluted to 106 pg and coinjected with myrf MO into the egg yolks during the one-cell stage. Zebrafish larvae were anesthetized with 0.03% tricaine at 72 hpf to observe the phenotypes. The ocular phenotypes of zebrafish larvae were analyzed at 72 hpf by light microscopy. The eye sizes of the larvae, defined as the distance from the nasal side to the temporal side of the eye, were measured using ImageJ software and were compared among different groups by a one-way ANOVA by SPSS23.0 software.
Images were collected with a Nikon C2 confocal microscope and a Zeiss 788 confocal microscope.

Identification of the genetic locus in family ZOC710536
In family ZOC710536, high hyperopia is transmitted as an autosomal dominant trait (Fig. 1), with male-to-male transmission across at least four generations. Clinical data and genomic DNA were obtained from 19 individuals, including 14 affected and 5 unaffected members, comprising 11 males and 8 females. Refraction in the 14 affected individuals ranged from + 5.00 to + 12.00 diopters, while refraction in 5 unaffected members ranged from 0 to − 3.00 diopters. Ocular axial length was available in three affected members: 17.60 mm/OD and 17.47 mm/OS for IV:7, 17.49 mm/ OD and 17.53 mm/OS for V:1, and 18.62 mm/OD and 18.76 mm/OS for V:2. All 19 family members had clear corneas with corneal diameters between 11 and 12 mm. Affected individuals had fundus changes typical of high hyperopia (Fig. 2). Strabismus was not observed in any of the 14 affected members. Angle-closure glaucoma in the right eye, with maximum intraocular pressure of 43 mmHg, had recently developed in one individual (IV:7) at the age of 38 years, 9 years after collection of the family data.
An initial genome-wide linkage scan was performed on 16 of the 19 family members (excluding IV:1, IV:3, and V:3; Fig. 1), to save funds and time (two rows of 16 samples fit well in our ABI3130 sequencer). Two-point linkage analysis revealed only one candidate locus on chromosome 11, with maximum LOD scores of 3.46, 3.17, and 2.72 for markers D11S4191, D11S987, and D11S905, respectively. Fine mapping with all 19 individuals (including IV:1, IV:3, and V:3 in Fig. 1) and additional markers confirmed the locus on chromosome 11, where ten closely spaced microsatellite markers in the region generated positive LOD scores, with D11S987, D11S1889, and D11S4191 yielding LOD scores of 4.68, 4.21, and 4.07, respectively, at theta = 0 (Table 1). Haplotype construction demonstrated that all patients in the family shared a haplotype block of D11S4102-D11S1360-D11S905-D11S1313-D11S4191-D11S4076-D11S4205-D11S1889-D11S987-D11S4113. The proximal and telomeric boundaries were indicated by recombination between D11S935 and D11S4102 in individuals III:7, III:9, IV:5, IV:6, IV:7, V:1, V:2, and V:3 as well as recombination between D11S4113 and D11S4196 in individual IV:5. The data strongly support that the disease gene is situated in a 27.1 cM region between D11S935 and D11S4196 on chromosome 11p12-q13.3 (Fig. 1).

Identification of a novel gene for high hyperopia in family ZOC710536
For parallel WES of genomic DNA from four affected individuals (III:5, III:9, IV:2, and V:1) in family ZOC710536, the average throughput depth, mean read depth, and coverage of the target regions were 195.5X, 49.9X, and 89.1% for III:5; 228.8X, 63.6X, and 90.2% for III:9; 191.7X, 57.3X, and 90.0% for IV:2; and 203.3X, 66.7X, and 90.1% for V:1. Comparison of sequencing data among the four individuals revealed that 14 rare variants affecting coding residues were shared by all four patients, including 12 missense mutations, one non-frameshift insertion, and one frameshift deletion. Of the 14, five were located inside the linkage interval on chromosome 11p12-q13.3. All but one of the 14 were excluded as candidate mutations for high hyperopia in the family after comparison of the 14 variants with data on the corresponding genes from ExAC and our in-house data from 3280 unrelated individuals. The remaining one is a novel variant located at chr11:61553536 in exon 27 (the last exon) of the myelin regulatory factor (MYRF) gene (HGNC ID: 1181; Entrez Gene: 745; Ensembl: ENSG00000124920; OMIM: 608329), also named C11ORF9 or KIAA0954, located at 11q12.2 [genomic location: chr11: 61,752,617-61,788,518 (GRCh38/hg38), or chr11: 61,520,121-61,555,990 (GRCh37/hg19)] that is situated inside the linkage interval. This potential pathogenic variant, a c.3377delG (p.Gly1126Valfs*31) based on NM_001127392, is predicted to result in the loss of the last 26 of the 1151 residues of MYRF and the addition of 30 new residues. Such variants are classified as loss-of-function (LoF) mutations, a category that includes nonsense, frameshift, splicing acceptor, and splicing donor variants; these mutations are extremely rare in MYRF, and the gene is, therefore, considered extremely LoF intolerant (pLI = 1.0) based on the ExAC database (http://exac.broad insti tute.org/ gene/ENSG0 00001 24920 , Sep 25th 2018). This variant is not present in ExAC or in any additional families among the 3280 with WES data. The c.3377delG mutation was further confirmed by Sanger sequencing (Fig. 3) and was completely cosegregated with high hyperopia in family ZOC710536.

Detection of MYRF mutations in two additional families with high hyperopia
Apart from the novel mutation identified in the large family mentioned above, two additional novel truncation variants in MYRF were detected in 2 of the other 121 probands with high hyperopia based on WES. Of the two, a c.3274_3275del (p. Leu1093Profs*22) mutation in exon 25 was detected in a 6-year-old boy with high hyperopia and axial lengths of 20.90 mm and 20.68 mm for the right and left eyes, respectively. The c.3274_3275del mutation was a de novo mutation that was not present in the boy's parents (Fig. 3), who had normal visual acuity and normal axial length. The remaining one, a c.3194+2T>C predicted to result in loss of the splicing donor in intron 24, was detected in a 40-year-old singleton female with high hyperopia and axial lengths of 19.73 mm and 19.83 mm. This individual had recently developed angle-closure glaucoma of both eyes, with a highest intraocular pressure of 42 mmHg.
Similar to the c.3377delG mutation, the c.3274_3275del and c.3194+2T>C mutations were confirmed by Sanger sequencing (Fig. 3) and were present in neither the ExAC database nor the rest of the 3280 probands from our in-house  Fig. 1 at 38 years old, 2 weeks after the onset of angle-closure glaucoma in the right eye (with IOP 43 mmHg before treatment). An enlarged optic disc was observed compared with the fundus photograph in D. f Heidelberg retina tomograph (HRT) results for the affected individual IV:7 from Fig. 1. The HRT result in F demonstrated partial loss of the retinal ganglion cell layer in the temporal half of the retina in the right eye. OD right eye, OS left eye 1 3 WES data. Except for the three novel truncation mutations described in the above three unrelated families, no additional truncation mutations were found in WES data from the 3280 probands with other forms of genetic eye diseases. A number of missense variants were also detected in MYRF by WES, but they were distributed in patients with different genetic eye diseases, as well as controls, without enrichment in any special subgroup.

Knockdown of myrf resulted in small eye size
To further analyze the function of myrf in eye development, we developed a myrf knockdown zebrafish model by injecting of myrf MO into zebrafish eggs to block translation of the gene. The myrf MO showed high efficiency in knocking down expression of the myrf gene in zebrafish, as no GFP signal was detected in embryos at 12 hpf after coinjection of 254 pg of myrf-pEGFPN1 plasmid and 3 ng of myrf MO coinjected embryos, while a strong GFP signal was detected in embryos injected with 254 pg myrf-pEGFPN1 alone or in conjunction with 3 ng of std MO (Fig. 4a, c). The mean (± standard error, SE) was 296.10 (± 11.68) μm of 96 larvae injected with 3 ng of std MO at 72 hpf and 279.07 (± 18.03) μm of 111 larvae injected with 3 ng of myrf MO at 72 hpf, respectively. A significant reduction in eye size was detected in larvae injected with 3 ng of myrf MO at 72 hpf than larvae injected with 3 ng of std MO (p = 7.11 × 10 −14 ) (Fig. 4b,  d). No other obvious bodily abnormalities were found in the zebrafish larvae. To further confirm the specificity of myrf MO against myrf, we synthesized myrf mRNA with a 3-nt silent mutation in the myrf MO targeting sequence and coinjected it with myrf MO for the rescue experiments. Indeed, the mean (± SE) was 310.25 (± 14.81) μm of 66 larvae coinjected with myrf mRNA and myrf MO at 72 hpf, which significantly larger than that of larvae injected with myrf MO alone (p = 5.56 × 10 −25 ) (Fig. 4b, d). These results suggested that the phenotype of small eye size in myrf morphants was specifically caused by myrf gene knockdown in zebrafish. Taken together, myrf MO knockdown and the rescue experiments illustrated that myrf played an important role in zebrafish eye development.
Since a phenotype of small eye size was detected in myrf morphants, we wanted to detect whether myrf affected retinal cell development. Immunofluorescent staining for different retinal cell labels showed that all of the ganglion cells, amacrine cells, bipolar cells, Müller cells, cone bipolar cells, and long double cone outer segments in the retina developed normally in myrf morphants compared with wild-type larvae at 5 pdf (Fig. 5).

Discussion
In this study, autosomal dominant high hyperopia in a large Chinese family was mapped to chromosome 11p12-q13.3 between D11S935 and D11S4196. The exclusion of other regions in the genome, the maximum positive LOD score of 4.68 together with haplotype support, the detection of a novel truncation mutation in MYRF that completely segregated with high hyperopia in the family, the identification of two other novel truncation mutations in the same gene in two additional families, the extreme rarity of truncation mutations in this gene, and the exclusion of potential pathogenic mutations in other genes, and the results from myrf knockdown zebrafish all support the case that truncation mutations in the C-terminal portion in MYRF are responsible for autosomal dominant high hyperopia in these families. Previously, an autosomal dominant form of nanophthalmos (NNO1) with high hyperopia and angle-closure glaucoma was mapped to chromosome 11p between D11S905 and D11S987, with a maximum LOD score of 5.92 (Othman et al. 1998). The linkage interval of NNO1 overlaps with the linkage interval identified in our current study. However, the MYRF gene is located outside the linkage interval of NNO1; thus, the gene responsible for NNO1 has yet to be identified. Therefore, MYRF-associated high hyperopia should be considered to have a separate locus. Nevertheless, MYRF-associated high hyperopia and NNO1 share common clinical features, such as high hyperopia, short axial length, and angle-closure glaucoma. It would be interesting to know the sequence variations of MYRF in the affected individuals of the family with NNO1. Meanwhile, the association of angle-closure glaucoma with MYRF mutations is of special interest for future studies because the molecular basis of angle-closure glaucoma is largely unknown, although this condition is a common cause of blindness, especially in Asian populations (Sun et al. 2017).
The initial MYRF full-length cDNA was cloned from an eye-tissue-specific cDNA library and was expressed in the retina and several other tissues (Stohr et al. 2000). Updated annotation of the MYRF gene shows two variant transcripts, of which transcript variant 2 is longer, with 27 exons coding for 1151 residues, while transcript variant 1 has 26 exons coding for 1111 residues (https ://www.ncbi.nlm.nih. gov/nucco re/NC_00001 1.10?from=61752 617&to=61788 518&repor t=genba nk). Transcript variant 1 has a different 5′ terminal exon, lacks an exon corresponding to exon 19 in transcript variant 2, and has different splicing acceptors for exons 19 and 22 (corresponding to exons 20 and 23 of transcript variant 2). All three novel truncation mutations identified in the current study are located in exons shared by the two variant transcripts (splicing donor of intron 24, exon 25, and exon 27). MYRF is widely expressed in mouse tissue (Nakayama et al. 2018) and is enriched in the central nervous system (Cahoy et al. 2008;Emery et al. 2009). This gene encodes an endoplasmic reticulum membrane-anchored transcription factor that is specifically expressed in oligodendrocytes and is critical for oligodendrocyte differentiation, as well as for the generation and maintenance of myelination in the central nervous system (Bujalka et al. 2013;Duncan et al. 2017;Emery et al. 2009;Kim et al. 2017;Koenning et al. 2012;Li and Richardson 2016). Through autocleavage, the nuclear-targeted N-terminal portion, which contains a DNA-binding domain, is separated from the C-terminal portion and then translocated from the endoplasmic reticulum to the nucleus, where it directly binds the enhancer regions of genes underlying myelination in oligodendrocytes (Bujalka et al. 2013;Li et al. 2013). MYRF is also indispensable for synaptic plasticity in C. elegans (Meng et al. 2017). The role of the C-terminal region is still unknown and may potentially involve MYRF folding, localization, cleavage, or regulation (Bujalka et al. 2013). Interestingly, all three novel truncation mutations identified in the current study are located in the C-terminal region of MYRF (Fig. 3).
Thus far, mutations in MYRF have been reported to have possible associations with several diseases, although evidence is not confirmatory from a genetic point of view. Missense variants in MYRF have been reported to have a significant association (p < 0.05) with late-onset Alzheimer's disease (Vardarajan et al. 2018). Two SNPs near MYRF have been reported to associate with blood systolic pressure based on genotyping of approximately 200,000 SNPs in 1947 black South African individuals (rs11230796, p = 2.16 × 10 −7 ; rs400075, p = 2.88 × 10 −7 ) (Hendry et al. 2018). In addition, three heterozygous de novo truncation variants in MYRF have been detected in three singleton males with cardiac and urogenital anomalies (Chitayat et al. 2018;Pinz et al. 2018). Furthermore, a heterozygous c.1208A>G (p.Gln403Arg) mutation has been detected in two families with encephalopathy with reversible myelin vacuolization, with incomplete penetrance (Kurahashi et al. 2018). Mice with homozygous conditional knockout of MYRF in oligodendrocytes demonstrate severe tremors and ataxia from postnatal day 11, subsequently develop seizures, and die during the third postnatal week (Emery et al. 2009). However, the phenotype of mice with heterozygous conditional knockout has not been described (Emery et al. 2009). The nature of the variants described in our current study is different from most of what is described in the above-mentioned reports, except for the three males with cardiac and urogenital anomalies, and the location of variants is also different. Our patients with high hyperopia do not present systemic signs or symptoms suggesting brain, cardiac or urogenital defects. All the three truncation variants identified in families with high hyperopia in the current study are located at the C-terminal of MYRF and absent in gnomAD database. In contrast, the other truncation variants in MYRF are all located upstream of the three variants, including three variants identified in patients with cardiac and urogenital anomalies (Table S2, Fig. 3b) and additional 19 LoF variants in gnomAD database. Therefore, it might be a gain-of-function model by Fig. 3 Sequence chromatography, distribution, and conservation of the three variants in MYRF. a Sequence chromatography from Sanger sequencing. Variants in MYRF identified by exome sequencing were further confirmed with Sanger sequencing and then validated in available family members. Pedigrees of ZOC710972 and ZOC712464 are shown in the left column. Under each individual, + indicates a wild-type allele, and M indicates a mutant allele. The sequences with or without variants detected in patients with high hyperopia, available parents, as well as normal controls are shown in the middle and right columns. b Distribution of variants in MYRF and exon coverage corresponding to the predicted protein schematic of MYRF. Mutations associated with high hyperopia are indicated in red, while mutations associated with cardiac and urogenital anomalies are indicated in blue. a′, b′ Predicted schematic of MYRF referring to the NM_013279.3 and NM_001127392.2 isoforms, respectively (Pinz et al. 2018). DBD DNA-binding domain, ICA intramolecular chaperone autoprocessing domain, TM transmembrane domain, S74 peptidase S74, chaperone of endosialidase domain, MRF C myelin gene regulatory factor C-terminal domain. c′, d′ Exon coverage of MYRF, referring to the NM_013279.3 and NM_001127392.2 isoforms, respectively. Exons differing between the two isoforms are indicated in red. c Protein sequence alignment of seven MYRF orthologs. The residues differing from the human sequence are indicated with boxes. The three variants and the associated new residue sequences resulted from them are shown above and indicated with red ◂ the three variants in patients with high hyperopia in the current study, especially one (c.3377delG) is located in the last exon which could escape the nonsense-mediated decay.
The phenotypes associated with different MYRF variants may need further clarification, although our current study provides the strongest genetic evidence to date in support Fig. 4 Phenotype of myrf knockdown zebrafish. a The efficiency of myrf MO knockdown of the myrf gene. Strong GFP signals were present in the left and middle images, which represent injection with myrf-pEGFPN1 plasmid alone and coinjected of myrf-pEGFPN1 plasmid with std MO, respectively. No GFP signal was detected in the embryos coinjected with myrf-pEGFPN1 plasmid and myrf MO, as shown in the right image. b Embryos with myrf knockdown showed a phenotype of small eye size (middle image) compared to the std MOinjected embryos (left image). This phenotype could be rescued by coinjection of myrf mRNA (right image). Control is larvae injected with 3 ng of std MO at 72 hpf. Mutant is larvae injected 3 ng of myrf MO at 72 hpf. Rescue is larvae coinjected with 3 ng of myrf mRNA and 106 pg of myrf MO at 72 hpf. c Proportions of embryos with GFP + or GFP − treated with 245 pg pEGFPN1-myrf, co-injection of 3 ng std MO with 245 pg of pEGFPN1-myrf and coinjection of 3 ng myrf MO with 245 pg of pEGFPN1-myrf. The p values by the Chisquared test are shown above. d Box plot of the median and deviations of eye size in three groups. The p values by one-way ANOVA are shown above of the association between MYRF truncation mutation and high hyperopia.
A disease called myelinated retinal nerve fibers (MRNF), consisting of abnormal myelination with ectopic oligodendrocyte-like cells in the retina, is a common congenital anomaly that is observed in approximately 1% of all eyes (Nangia et al. 2014;Tarabishy et al. 2007). This condition is frequently associated with ipsilateral myopia or hyperopia (Lee and Salchow 2008;Nangia et al. 2014;Tarabishy et al. 2007). A recent study demonstrated a significant association between MRNF and hyperopia based on an analysis of 52 adult eyes with MRNF in rural Central India (Nangia et al. 2014). It is still unknown whether variants in MYRF contribute to MRNF, although MYRF plays a critical role in myelination. SOX10 regulates and cooperates with MYRF (Hornig et al. 2013;Lopez-Anido et al. 2015). Mutations in SOX10 cause Waardenburg syndrome, which is characterized by developmental anomalies of the eye (Pingault et al. 1998(Pingault et al. , 2010. These data also imply that mutations in MYRF may contribute to developmental anomalies of the eye, including high hyperopia with short axial length.
The molecular mechanism underlying high hyperopia development with MYRF truncation mutations is not yet known. Retinal progenitor neuron cells and specifically expressed functional genes in these cells may not only play roles in defining specific ocular structure or biological functions but also participate in the developmental process of eye development. Further study on these truncation mutations in the C-terminal portion of MYRF may not only reveal the molecular pathway contributing to high hyperopia but also provide a valuable avenue for investigating the role of MYRF C-terminal portion.
In summary, high hyperopia in a large family is mapped to chromosome 11p12-q13.3, a locus overlapping with NNO1, by genome-wide linkage analysis. WES and wholegenome sequencing identified novel truncation mutations in the C-terminal portion of MYRF, located inside the linkage interval of high hyperopia but outside the linkage interval of NNO1, in three unrelated families including the family whose mutation is mapped to this region. Informatic approaches using mutation prediction and a mutation database, as well as an experimental approach using knockdown of myrf in zebrafish, further support linkage mapping and mutation identification. Our findings may provide useful clues to further elucidate the function of the critical myelin regulatory factor MYRF as well as the molecular pathogenesis of high hyperopia and its associated angle-closure glaucoma.