Abstract

A growing number of human diseases have been linked to defects in protein glycosylation that affects a wide range of organs. Among them, O-mannosylation is an unusual type of protein glycosylation that is largely restricted to the muscular and nerve system. Consistently, mutations in genes involved in the O-mannosylation pathway result in infantile-onset, severe developmental defects involving skeleton muscle, brain and eye, such as the muscle−eye−brain disease (MIM no. 253280). However, the functional importance of O-mannosylation in these tissues at later stages remains largely unknown. In our study, we have identified recessive mutations in POMGNT1, which encodes an essential component in O-mannosylation pathway, in three unrelated families with autosomal recessive retinitis pigmentosa (RP), but without extraocular involvement. Enzymatic assay of these mutant alleles demonstrate that they greatly reduce the POMGNT1 enzymatic activity and are likely to be hypomorphic. Immunohistochemistry shows that POMGNT1 is specifically expressed in photoreceptor basal body. Taken together, our work identifies a novel disease-causing gene for RP and indicates that proper protein O-mannosylation is not only essential for early organ development, but also important for maintaining survival and function of the highly specialized retinal cells at later stages.

Introduction

Retinitis pigmentosa (RP) is a Mendelian disorder featured by night blindness, tunnel vision, progressive retinal degeneration and eventually complete blindness. Its prevalence worldwide is ∼1 in 3000–7000 (1). The genetic etiology of RP is highly heterogeneous, with at least 79 disease-causing genes identified (RetNet, https://sph.uth.edu/retnet/). However, there are ∼40% RP cases with unknown causative mutations, strongly suggesting novel genes involved in RP pathogenesis (2). RP can be linked to syndromic disorders since genes important for retina also play roles in other tissues (3). For example, defects of various lysosomal proteins in patients with neuronal ceroid lipofuscinosis (MIM no. 256730) have a wide range of neurological and behavioral abnormalities as well as retinal degeneration (4). Also, the loss of pantothenate kinase leads to HARP syndrome (MIM no. 607236) featured by central nervous system problems, hematological disorders and RP (5). Similar cases can be seen in other genetic disorders such as Wolfram syndrome (MIM no. 222300), hyper-IgD syndrome (MIM no. 260920) (6) and a series of ciliopathies (7).

With the advent of whole-exome sequencing (WES) era, large-scale WES in multiple cohorts have greatly improved our understanding about the molecular etiology of RP. Specifically, a number of novel genotype–phenotype correlations have been identified. In these studies, the phenotype spectrum of the mutations in syndrome-causing genes was expanded to non-syndromic RP. For instance, MVK mutations were known to cause hyper-IgD syndrome (8), which is a systemic inflammatory disorder. However, recent reports show that MVK defects also lead to non-syndromic RP (9). Another example is CLRN1, whose mutations were majorly associated with Usher syndrome (MIM no. 276902) (10), is also linked to non-syndromic RP (11). This phenomenon was observed in a long list of recent reports involving genes such as WDR19 (12), IFT172 (13), IFT140 (14), NEUROD1 (15) and HGSNAT (16). Although these reports were sporadic findings during the genetic screening of RP, they together account for a considerable proportion of RP cases. Furthermore, the identification of these genes will improve our ability of RP molecular diagnosis and our understanding of the basic biology of retina.

Glycosylation is a highly diverse category of protein post-translational modifications. It plays essential roles in regulating the protein subcellular localization, transport and quality control (17). Mutations in DHDDS, encoding an enzyme that plays a role in N-glycosylation, were reported to cause non-syndromic RP (18,19). Another specific type of glycosylation, protein O-mannosylation, was originally detected in the brains and muscles of mammals (20–22) and later found to be functionally important in these tissues (23,24). Correspondingly, disruption of enzymes in the O-mannosylation pathway will lead to a spectrum of genetic disorders majorly affecting human neuronal and muscular systems (25). One of the key member in the pathway is POMGNT1, a β-1,2-N-acetylglucosaminyltransferase enzyme that participates in O-mannosyl glycan synthesis (26). Recessive loss-of-function POMGNT1 mutations can cause muscle–eye–brain disease (MEB, MIM no. 253280) (27), featured by extensive muscular, neurological and ocular abnormalities (28,29). Specifically, a series of ocular symptoms are seen in patients including retinal degeneration, optic nerve atrophy, impaired electroretinogram (ERG) signals, lens opacification, myopia and strabismus (30). The retinal abnormalities documented in MEB patients suggest the possibility that like other ‘syndrome-causing genes with retinal involvement’, POMGNT1 mutations might also be associated with non-syndromic retinal disease.

In this study, we identified recessive POMGNT1 mutations in three unrelated non-syndromic RP families by WES. Further biochemical studies showed that the mutations significantly impaired POMGNT1 enzymatic activity. Immunohistochemical analysis documented POMGNT1 expression in the photoreceptor layer. Our results greatly expand the phenotype spectrum of POMGNT1 mutations and provide the first link between the protein O-mannosylation pathway and non-syndromic RP.

Results

Three non-syndromic RP families with POMGNT1 mutations were identified

The index case we have investigated belongs to a medium-size autosomal recessive RP family of Italian origin diagnosed at the McGill University Health Centre (MUHC), Canada. The proband (MOGL2063) is a 69-year-old male. He was aware of visual problems at the age of 10 years old. He currently has a visual acuity of 20/100 and very restricted visual fields in both eyes (Goldmann visual field test shows 5°). His fundus images show bone spicules, narrow vessels, peripapillary atrophy and a tigroid appearance. The proband's sister is 78 years old (MOGL2064) and also affected with RP. She was aware of nyctalopia and vision problems when she was 12 years old. She currently has visual acuities of 20/80 OD and light perception OS, constricted visual fields (5°) (Fig. 1A and B) and peripapillary atrophy. Funduscopy featured bone spicules and attenuated retinal vessels. Optical coherence tomography (OCT) revealed chorioretinal atrophy, retinal thinning with absent inner segment/outer segment (IS/OS) junctions and cystoid macular edema in both eyes (Fig. 1C and D).

Clinical findings of non-syndromic RP patients in this study. Visual field test results of MOGL2064, right eye (A) and left eye (B), showing restricted visual field in both eyes. OCT images of MOGL2064, right eye (C) and left eye (D), featuring chorioretinal atrophy, retinal thinning with absent IS/OS junctions. Fundus images of SRF1105, right eye (E) and left eye (F), showing the tigroid appearance, thinning of mid-peripheral RPE and loss of choroid-capillaries. OCT images of SRF1105, right eye (G) and left eye (H), showing thinning of outer retinal layer with preserved foveal photoreceptors. (I) ERG results of SRF1105, showing no detectable ERG signals. Fluorescein angiography results of SRF1105, right eye (J) and left eye (K), showing oval shaped hypofluorescence with a hyperfluorescent ring in the macular and mottled hypofluorescence area in the peripheral retina. OCT images of SRF1630, right eye (L) and left eye (M), showing thinning of whole retina layers, a flat fovea and disappearance of IS/OS (right eye more severe than left eye). Visual field test results of patient SRF1630, right eye (N) and left eye (O), showing tunnel vision in both eyes.
Figure 1.

Clinical findings of non-syndromic RP patients in this study. Visual field test results of MOGL2064, right eye (A) and left eye (B), showing restricted visual field in both eyes. OCT images of MOGL2064, right eye (C) and left eye (D), featuring chorioretinal atrophy, retinal thinning with absent IS/OS junctions. Fundus images of SRF1105, right eye (E) and left eye (F), showing the tigroid appearance, thinning of mid-peripheral RPE and loss of choroid-capillaries. OCT images of SRF1105, right eye (G) and left eye (H), showing thinning of outer retinal layer with preserved foveal photoreceptors. (I) ERG results of SRF1105, showing no detectable ERG signals. Fluorescein angiography results of SRF1105, right eye (J) and left eye (K), showing oval shaped hypofluorescence with a hyperfluorescent ring in the macular and mottled hypofluorescence area in the peripheral retina. OCT images of SRF1630, right eye (L) and left eye (M), showing thinning of whole retina layers, a flat fovea and disappearance of IS/OS (right eye more severe than left eye). Visual field test results of patient SRF1630, right eye (N) and left eye (O), showing tunnel vision in both eyes.

We first performed target capture sequencing on patient MOGL2064 to see if she has mutations in known RP-causing genes. The capture sequencing data did not identify any putative causative variants (data not shown). Then we performed WES on both affected individuals and particularly searched for biallelic variants shared by these two individuals. After data analysis, variants in POMGNT1 were prioritized as top candidate due to its known link with a retinal phenotype. Both patients have a stop-gain variant (c.187C>T; p.R63*) and a missense variant (c.860T>G; p.I287S) in POMGNT1. The R63* variant is absent in the ExAC database (http://exac.broadinstitute.org) and the I287S variant has a frequency of 1 in 40 000, indicating they are rare in the population (Table 1). The stop-gain variant completely abolishes the catalytic domain of POMGNT1 and should confer no enzymatic activity according to previous structure–function studies (31). The I287S variant affects an amino acid (AA) site highly conserved in vertebrates (Fig. 2B) and it is considered damaging by all prediction algorithms (Table 2). The missense variant was not reported in MEB patients before. We performed Sanger sequencing and confirmed that the two variants are in trans and two unaffected siblings are wild-type (WT) or single heterozygous at these two loci (Fig. 2A, Supplementary Material, Fig. S1). We revisited the patients to see if they have muscular or neurological abnormalities as shown in MEB patients. Interestingly, both patients (now 69 and 78 years old) have no extraocular disease. These results strongly suggest that POMGNT1 mutations are associated with non-syndromic RP.

Table 1.

Clinical features and POMGNT1 variants of RP patients in this study

ID/ethnicityDxSexAgeMSBROtherPOMGNT1 variantsO FreqEM Freq
MOGL2063/CauRPM69NoNoNoc.187C>T (p.R63*)
c.860T>G (p.I287S)
Absent
3 in 121 400
Absent
1 in 66 738
MOGL2064/CauRPF78NoNoNoc.187C>T (p.R63*)
c.860T>G (p.I287S)
Absent
3 in 121 400
Absent
1 in 66 738
SRF1105/ChiRPF32NoNoNoc.466G>A (p.E156K) (homozygous)AbsentAbsent
SRF1630/ChiRPM52NoNoNoc.1895 + 1G>A (p.V633*)
c.1505G>C (p.G502A)
Absent
Absent
Absent
Absent
ID/ethnicityDxSexAgeMSBROtherPOMGNT1 variantsO FreqEM Freq
MOGL2063/CauRPM69NoNoNoc.187C>T (p.R63*)
c.860T>G (p.I287S)
Absent
3 in 121 400
Absent
1 in 66 738
MOGL2064/CauRPF78NoNoNoc.187C>T (p.R63*)
c.860T>G (p.I287S)
Absent
3 in 121 400
Absent
1 in 66 738
SRF1105/ChiRPF32NoNoNoc.466G>A (p.E156K) (homozygous)AbsentAbsent
SRF1630/ChiRPM52NoNoNoc.1895 + 1G>A (p.V633*)
c.1505G>C (p.G502A)
Absent
Absent
Absent
Absent

Chi, Han Chinese; Cau, Caucasian; Dx, primary diagnosis; MS, if the patient has muscular abnormalities; BR, if the patient has brain abnormalities; Other, if the patient has other disorders. O Freq, overall frequency, occurrence in ExAC database (all population groups); EM Freq, ethnically matched frequency, occurrence in ExAC database within the ethnically matched (EM) population group. Specifically in the ExAC database, East Asian is the EM population group for Han Chinese patients. Non-Finnish European is the EM population group for Caucasian patients; cDNA and protein changes are based on reference cDNA sequence NM_017739.

Table 1.

Clinical features and POMGNT1 variants of RP patients in this study

ID/ethnicityDxSexAgeMSBROtherPOMGNT1 variantsO FreqEM Freq
MOGL2063/CauRPM69NoNoNoc.187C>T (p.R63*)
c.860T>G (p.I287S)
Absent
3 in 121 400
Absent
1 in 66 738
MOGL2064/CauRPF78NoNoNoc.187C>T (p.R63*)
c.860T>G (p.I287S)
Absent
3 in 121 400
Absent
1 in 66 738
SRF1105/ChiRPF32NoNoNoc.466G>A (p.E156K) (homozygous)AbsentAbsent
SRF1630/ChiRPM52NoNoNoc.1895 + 1G>A (p.V633*)
c.1505G>C (p.G502A)
Absent
Absent
Absent
Absent
ID/ethnicityDxSexAgeMSBROtherPOMGNT1 variantsO FreqEM Freq
MOGL2063/CauRPM69NoNoNoc.187C>T (p.R63*)
c.860T>G (p.I287S)
Absent
3 in 121 400
Absent
1 in 66 738
MOGL2064/CauRPF78NoNoNoc.187C>T (p.R63*)
c.860T>G (p.I287S)
Absent
3 in 121 400
Absent
1 in 66 738
SRF1105/ChiRPF32NoNoNoc.466G>A (p.E156K) (homozygous)AbsentAbsent
SRF1630/ChiRPM52NoNoNoc.1895 + 1G>A (p.V633*)
c.1505G>C (p.G502A)
Absent
Absent
Absent
Absent

Chi, Han Chinese; Cau, Caucasian; Dx, primary diagnosis; MS, if the patient has muscular abnormalities; BR, if the patient has brain abnormalities; Other, if the patient has other disorders. O Freq, overall frequency, occurrence in ExAC database (all population groups); EM Freq, ethnically matched frequency, occurrence in ExAC database within the ethnically matched (EM) population group. Specifically in the ExAC database, East Asian is the EM population group for Han Chinese patients. Non-Finnish European is the EM population group for Caucasian patients; cDNA and protein changes are based on reference cDNA sequence NM_017739.

Table 2.

Predicted damaging effect of POMGNT1 missense variants identified in this study.

POMGNT1 missense variantsSIFTPolyphen2_HDIVPolyphen2_HVARLRTMutation TasterMutation AssessorFATHMM
c.860T>G (p.I287S)DPDPDDeDiMD
c.466G>A (p.E156K)DPDPDDeDiMD
c.1505G>C (p.G502A)DPDPDDeDiMD
POMGNT1 missense variantsSIFTPolyphen2_HDIVPolyphen2_HVARLRTMutation TasterMutation AssessorFATHMM
c.860T>G (p.I287S)DPDPDDeDiMD
c.466G>A (p.E156K)DPDPDDeDiMD
c.1505G>C (p.G502A)DPDPDDeDiMD

SIFT, scale-invariant feature transform; Polyphen2, polymorphism phenotyping v2; LRT, likelihood-ratio test; FATHMM, functional analysis through hidden Markov models; D, damaging; PD, probably damaging; De, deleterious; Di, disease-causing; M, medium damaging.

cDNA and protein changes are based on reference cDNA sequence NM_017739.

Table 2.

Predicted damaging effect of POMGNT1 missense variants identified in this study.

POMGNT1 missense variantsSIFTPolyphen2_HDIVPolyphen2_HVARLRTMutation TasterMutation AssessorFATHMM
c.860T>G (p.I287S)DPDPDDeDiMD
c.466G>A (p.E156K)DPDPDDeDiMD
c.1505G>C (p.G502A)DPDPDDeDiMD
POMGNT1 missense variantsSIFTPolyphen2_HDIVPolyphen2_HVARLRTMutation TasterMutation AssessorFATHMM
c.860T>G (p.I287S)DPDPDDeDiMD
c.466G>A (p.E156K)DPDPDDeDiMD
c.1505G>C (p.G502A)DPDPDDeDiMD

SIFT, scale-invariant feature transform; Polyphen2, polymorphism phenotyping v2; LRT, likelihood-ratio test; FATHMM, functional analysis through hidden Markov models; D, damaging; PD, probably damaging; De, deleterious; Di, disease-causing; M, medium damaging.

cDNA and protein changes are based on reference cDNA sequence NM_017739.

Genetic findings of non-syndromic RP patients in this study. (A) Four non-syndromic RP patients from three unrelated families with POMGNT1 biallelic variants were identified. Genotypes (cDNA and protein changes based on reference cDNA sequence NM_017739) are labeled for each family. (B) POMGNT1 protein sequence alignment among selected vertebrates at the loci of POMGNT1 missense variants identified in our study. All the disrupted amino acids are highly conserved in vertebrates.
Figure 2.

Genetic findings of non-syndromic RP patients in this study. (A) Four non-syndromic RP patients from three unrelated families with POMGNT1 biallelic variants were identified. Genotypes (cDNA and protein changes based on reference cDNA sequence NM_017739) are labeled for each family. (B) POMGNT1 protein sequence alignment among selected vertebrates at the loci of POMGNT1 missense variants identified in our study. All the disrupted amino acids are highly conserved in vertebrates.

To test our hypothesis that POMGNT1 is a novel gene for non-syndromic RP, we then screened POMGNT1 mutations in the WES data of 308 unsolved RP patients. Strikingly, we identified another two cases. Both patients were Han Chinese and diagnosed in Peking Union Medical College Hospital (PUMCH). The first patient is a 32-year-old female born in a consanguineous family (Fig. 2A). She suffered night blindness since childhood. Her best corrected visual acuity is 20/25 OD and 20/50 OS. Her diopter is 6.00–2.75 × 169° OD and −9.75 to 3.00 × 6° OS. Her fundus showed a tigroid appearance, which possibly indicates thinning of mid-peripheral retinal pigment epithelium (RPE) and loss of choroid-capillaries (Fig. 1E and F). OCT revealed thinning of outer retinal layer with preserved foveal photoreceptor IS/OS junctions (Fig. 1G and H). Visual field testing demonstrated constricted visual field of both eyes. The ERG examination displayed non-detectable waves OU (Fig. 1I). Fundus auto-fluorescence result showed an oval shaped hypofluorescence with a hyperfluorescent ring in the macula and mottled hypofluorescence area in the peripheral retina (Fig. 1J and K). This patient was re-visited, but no syndromic abnormalities were identified. Her serum creatine kinase level is normal (60 U/l, reference: 24∼170 U/l) and there is no evidence of periventricular white matter abnormality, ventriculomegaly, pontocerebellar hypoplasia or cerebellar cyst in magnetic resonance images. WES data show that she possesses a homozygous missense variant (c.466G>A; p.E156K) in POMGNT1 (Fig. 2A). The missense variant is not found in any control database suggesting it is extremely rare (Table 1). It disrupts a highly conserved AA site (Fig. 2B) and considered damaging by all prediction algorithms (Table 2). The variant is also not found in the list of known POMGNT1 mutations identified in MEB patients.

The other patient is a 52-year-old male. He presented with night blindness since his youth. Currently he has a visual acuity of 20/400 OD, 20/40 OS. OCT features thinning of the whole retina with a flat fovea and absent IS/OS junctions (Fig. 1L and M). Bone-spicule and salt-and-pepper pigment proliferation were scattered in the whole retina with macular involvement in both eyes. Visual field testing demonstrated tunnel vision in both eyes (Fig. 1N and O). He has a splicing (c.1895 + 1G>A; p.V633*) and a missense (c.1505G>C; p.G502A) POMGNT1 variant inherited from one parent each. The splicing variant was previously reported in MEB patients and it causes intron retention and immediately generates a premature stop codon (32). The missense variant is absent in control databases, conserved among vertebrates and also predicated to be damaging by all algorithms (Fig. 2B, Tables 1 and 2).

In addition, we performed homozygosity mapping using SNP arrays and identified five index cases with homozygous regions that include POMGNT1 (Supplementary Material, Table S1). However, mutation analysis of the whole-coding sequence using Sanger sequencing did not reveal any possible pathogenic mutation in these patients.

Functional analysis suggested the hypomorphic nature of the POMGNT1 missense variants

We performed biochemical analyses to further test the pathogenicity of POMGNT1 missense variants we identified. Specifically, we did site-directed mutagenesis on POMGNT1-expressing plasmids and transfected them into HEK293T cells. The culture supernatants were used for western blot analysis and POMGNT1 enzymatic activity measurement. Western blot results showed the successful expression of POMGNT1 proteins and their approximately equal molecular weights (Fig. 3A). During the enzymatic assay, we repeated the experiments five times. Compared with the WT POMGNT1 protein control, the two mutant POMGNT1 proteins we assayed (E156K and I287S) have significantly reduced (P < 10−5) enzymatic activity, with only 30 and 10% of the WT level retained, respectively (Fig. 3B). This indicated that the POMGNT1 variants identified in non-syndromic RP cases were hypomorphic mutations instead of complete loss-of-function ones.

Biochemical studies to test the pathogenicity of two POMGNT1 missense variants. (A) Western blot analysis of POMGNT1 WT and mutant proteins, showing their expression and approximately equal molecular weights. (B) Enzymatic activity of POMGNT1 WT and mutant proteins. Both mutants retained significantly lower activity compared with the WT level. Two-tailed Welch's t-test. ***P < 10−5. Error bars indicate sample standard deviations.
Figure 3.

Biochemical studies to test the pathogenicity of two POMGNT1 missense variants. (A) Western blot analysis of POMGNT1 WT and mutant proteins, showing their expression and approximately equal molecular weights. (B) Enzymatic activity of POMGNT1 WT and mutant proteins. Both mutants retained significantly lower activity compared with the WT level. Two-tailed Welch's t-test. ***P < 10−5. Error bars indicate sample standard deviations.

POMGNT1 is expressed in photoreceptor cells

We further performed immunohistochemistry to explore the localization of POMGNT1 protein in the retina. In adult mouse retina cryo-sections, POMGNT1 shows no expression in nuclear layers or plexiform layers (Supplementary Material, Fig. S2), but clear expression in photoreceptor cells (Fig. 4). Specifically, we found that POMGNT1 is localized at the basal body and daughter centriole of photoreceptor cells, suggesting functional involvement at these sites (Fig. 4).

POMGNT1 is expressed in photoreceptor cells. For each individual photoreceptor, the POMGNT1 signal is next to the cilia marker (acetylated-α-tubulin) signal. Two-dot pattern (labeled by double arrows) indicates the localization of POMGNT1 at both the basal body and the daughter centriole of photoreceptor cells. Ac-α-Tub, acetylated-α-tubulin; OS, outer segments; CC, connecting cilium; IS, inner segments; ONL, outer nuclear layer.
Figure 4.

POMGNT1 is expressed in photoreceptor cells. For each individual photoreceptor, the POMGNT1 signal is next to the cilia marker (acetylated-α-tubulin) signal. Two-dot pattern (labeled by double arrows) indicates the localization of POMGNT1 at both the basal body and the daughter centriole of photoreceptor cells. Ac-α-Tub, acetylated-α-tubulin; OS, outer segments; CC, connecting cilium; IS, inner segments; ONL, outer nuclear layer.

Discussion

POMGNT1 catalyzes the formation of GlcNAcβ1-2Man linkage of O-mannosyl glycans (27). O-mannosyl glycans play roles in the assembly and organization of the basal membranes in muscular and nervous systems, thus crucial to motor activity and neurodevelopment (33,34). As expected, Pomgnt1-null mice show typical MEB-like phenotype including muscular dystrophy, brain developmental defects and eye anomalies (35,36). Specifically, in the retina of Pomgnt1-null mice, the retinal layers are thinner than those in WT mice (36) and there is a reactive gliosis of astrocytes and Müller glial cells (35), which is also seen in other neurodegenerative disorders with retinal involvement (37). These mice also show abnormal retinal vasculature organization and reduced a- and b-wave ERG response levels (35,36), which are typical features of RP. The phenotypes in Pomgnt1-null mice indicated the essential role of POMGNT1 in the maintenance of retinal structures and functions. Previous studies have identified O-mannosylated dystroglycan and its binding partner Pikachurin are required for photoreceptor ribbon synapse formation (38,39), which is crucial for the neuronal connections between photoreceptors and second-order neurons. Here, we first showed that POMGNT1 is indeed expressed in photoreceptor cells, supporting its functional importance in this cell type. Interestingly, we found POMGNT1 is specifically localized near photoreceptor cilium basal body, a subcellular structure critical for protein transport from inner segments to outer segments, suggesting that POMGNT1 is not only critical for dystroglycan O-mannosylation, but also may modify outer segment proteins. Further proteomic studies would discover outer segment protein substrates of POMGNT1 and unravel the pathological mechanism of POMGNT1-associated retinal disease.

Recessive POMGNT1 mutations have been reported to cause MEB (27). The phenotype spectrum of MEB is highly variable, ranging from very severe Walker–Warburg syndrome-like phenotype with early lethality (40), to mild phenotype mimicking limb-girdle muscular dystrophy without mental retardation (41). In general, congenital or early-onset muscular dystrophy are the universal features presented, with most patients also developing neurological and ocular symptoms. Interestingly, the discoveries in our study marked a clearly different phenotype from the canonical phenotypes caused by POMGNT1 mutations.

The non-syndromic feature of RP patients in our study probably resulted from two reasons. First, the missense POMGNT1 mutations in our study are hypomorphic rather than the complete loss-of-function. The majority of the reported POMGNT1 causal alleles in MEB patients were protein-truncating mutations (42,43), which will lead to the loss of protein activity since it was shown that the deletion of 19 AAs in the C-terminal is sufficient to abolish POMGNT1 enzymatic activity (31). Another two studies showed that nearly all missense POMGNT1 MEB-causing variants also completely lose their activities (44,45). These results strongly argue that the total loss of POMGNT1 activity underlies typical MEB phenotype. While in our study, both missense variants we tested still maintains subnormal activity (30 and 10%), which seems plausible to correspond to the relatively milder non-syndromic phenotype. Interestingly, another report identified a homozygous hypomorphic POMGNT1 variant in a patient diagnosed with limb-girdle muscular dystrophy, a phenotype milder than MEB (41), further supporting our contention that the POMGNT1 enzymatic activity level play an important role in determining the phenotype. In fact, studies in other genetic disease have also observed a correlation between disease severity and residual enzyme activity (46–48). The specialized structure, highly dynamic protein turnover and high-energy consumption in retina (49) probably confer higher threshold of POMGNT1 activity than other tissues for maintaining normal structure and functions. This model was also proposed recently in other similar studies identifying RP-related genotype–phenotype correlation (2,16,50). The unique biological profiles in the retina also explains why retinal defects often lie in the mild extremity of the spectra of syndromic phenotypes (51,52).

Secondly, genetic backgrounds or even non-genetic factors probably affect the phenotypes of individuals with POMGNT1 variants. This has been proposed before (43,53,54) since remarkable phenotypic difference was observed in MEB patients, even with the same causative alleles in the same family. In addition, C57BL/6 Pomgnt1-null mice, which presumably possess the same genetic background, also show phenotypic variability in the eyes (35), further suggesting the notion that non-genetic factors have an impact on phenotype. Future studies on Pomgnt1-null mice in different genetic backgrounds or under various conditions might help to reveal the complicated mechanisms underlying phenotypic variability.

It should be noted that a number of MEB patients with POMGNT1 mutations show no ocular symptoms or only non-retinal ocular symptoms including myopia, anterior chamber malformation, microphthalmia, nystagmus, strabismus and cataract (30,43,55). One reason may be that those MEB patients were too young to develop progressive retinal degeneration. MEB patients have a lifespan between 10 and 30 years (40) while the patients in our study are at older ages. In addition, early developmental defects in the eye would hinder the presentation and comprehensive clinical examination of late-onset retinal phenotype. As for the non-retinal ocular symptoms, since O-mannosyl glycan exist exclusively in neuronal and muscular tissues (20–22), they could probably be caused by secondary effects of neurodevelopmental defects or extraocular muscle anomalies as suggested by previous literature (30,56–60).

In summary, the present study identified POMGNT1 as a novel disease-causing gene for non-syndromic RP. The significance of POMGNT1 specifically in the retina has been underappreciated before since muscular dystrophy and neurodevelopmental defects are the predominant symptoms and major cause for the short lifespan in MEB patients. Our results greatly expand the phenotype spectrum of POMGNT1 mutations and highlighted the functional importance of POMGNT1 and O-mannosyl glycan in the retina. WES analysis in our cohort suggested POMGNT1 mutations as a rare cause of non-syndromic RP. Nevertheless, future genetic screening of mutations in other genes involving O-mannosylation pathway might improve the molecular diagnosis of non-syndromic RP cohorts and enhance our understanding of glycobiology in the retina.

Materials and Methods

Clinical diagnosis of patients and DNA sample collection

The patients were diagnosed by ophthalmic examinations including best corrected visual acuity testing, fundus examination, optical coherence topography (OCT, 3D OCT-2000 Spectral Domain; Topcon, Tokyo, Japan), visual field tests (Octopus, Interzeag, Schlieren, Switzerland), autofluorescence (Spectralis HRA + OCT; Heidelberg, Germany) and ERG (RetiPort ERG system, Roland Consult, Wiesbaden, Germany). RP was diagnosed according to medical and family history, typical fundus and OCT features, visual field defects and attenuated or abolished ERG responses. Blood samples were obtained from all patients and their family members if available. DNA was extracted using QIAamp DNA Blood Midi Kit as instructed by the manufacturer (QIAGEN, Hilden, Germany). This study adhered to the Declaration of Helsinki and was approved by the Institutional Review Board at the MUHC and PUMCH. We obtained informed written consent from all participating individuals for this study.

DNA library preparation and next generation sequencing

Approximately 1 μg of the genomic DNA sample was sheared into fragments of 200–500 bp in length. The sheared fragments were blunt-end repaired and a single-adenine base was added to the 3′ends using Klenow exonuclease. Illumina adapters were ligated to the repaired ends and DNA fragments were polymerase chain reaction (PCR) amplified for eight cycles after ligation. In each capture reaction, 50 pre-capture DNA libraries were pooled together. The targeted DNA was captured by customized retinal disease gene panel for retinal capture sequencing to screen for variants in known disease-causing genes. The detailed information of the retinal capture panel was described previously (2). If no plausible causative variants were identified, the DNA was then captured by NimbleGenSeqCap EZ Hybridization and Wash kit (NimblegenSeqCap EZ Human Exome Library v.2.0) following the manufacturer's protocols for WES. Captured libraries were sequenced on Illumina HiSeq 2000 (Illumina, San Diego, CA, USA) as 100 bp paired-end reads according to the manufacturer's protocol.

Bioinformatics analysis

Paired-end sequencing reads were obtained for each sample. Reads were mapped to human reference genome hg19 using Burrows–Wheeler aligner (61). Base-quality recalibration, local realignment and variant calling were performed as previously described (2). Since RP are rare Mendelian disorders, variants with a frequency >1/200 (for a recessive model) or 1/1000 (for a dominant model) in a series of public databases and internal control databases were filtered out. The list of databases can be found in previous literature (2). We also retrieved variant frequencies from the Exome Aggregation Consortium (ExAC, http://exac.broadinstitute.org/) database. After frequency-based filtering, we filtered out synonymous variants, identified known retinal disease-causing variants and predicted the pathogenicity of variants using Scale-invariant feature transform (SIFT) (62), polymorphism phenotyping v2 (Polyphen2) (63), likelihood-ratio test (LRT) (64), MutationTaster (65), MutationAssessor (66) and functional analysis through hidden Markov models (FATHMM) (67) as previously described (2,68).

Sanger sequencing

For each suspected causative mutation, a 500 bp flanking sequence at both sides was obtained from the UCSC genome browser (hg19 assembly). RepeatMasker (http://www.repeatmasker.org/) was used to mask the repetitive sequences in human genome. Primer 3 (69) was used to design a pair of primers for generating a 400–600 bp PCR product to sequence the mutation site and at least 50 bp region surrounding it. After PCR amplification, the amplicons were sequenced on ABI 3730xl. Family members of patients were also Sanger-sequenced when available for confirming allele segregation. In addition, five cases with RP and a significant homozygous region covering POMGNT1 were screened for mutations using primers (Supplementary Material, Table S2) covering all coding exons that were designed using the UCSC ExonPrimer tool (https://genome.ucsc.edu/index.html).

POMGNT1 enzymatic assay

HEK293T cells were maintained in Dulbecco's modified Eagle's medium (Invitrogen, Carlsbad, CA, USA) supplemented with 10% fetal bovine serum (FBS), 2 mml-glutamine, 100 U/ml penicillin and 100 µg/ml streptomycin at 37°C with 5% CO2. The expression plasmids of soluble POMGNT1 were transfected into HEK293T cells, and the media were replaced with 10 ml fresh FBS-free Opti-MEM (Invitrogen) 2 days after transfection, and incubated for one more day. The FBS-free culture supernatants were subjected to western blot analysis and assay for POMGNT1 activity.

The proteins were separated by sodium dodecyl sulfate–polyacrylamide gel electrophoresis (10% gel) and transferred to a PVDF membrane. The membrane was blocked in PBS containing 5% skim milk and 0.5% Tween-20, incubated with anti-POMGNT1 C-terminus antibody (70), and subsequently treated with anti-rabbit IgG conjugated with horseradish peroxidase (GE Healthcare Bio-sciences Corp., Piscataway, NJ, USA). Blots were developed using an ECL kit (GE Healthcare). Purified sPOMGNT1 (70) was used as a mass standard to determine the amount of each protein. Optical density measurement of the bands (Image J software) produced a standard curve, which was used to calculate the concentration of each mutant protein.

POMGNT1 activity was obtained from the amount of [3H]GlcNAc transferred to an O-mannosyl peptide (26). The enzyme reaction were carried out in 20 µl reaction buffer containing 140 mm MES (pH 7.0), 0.5 mm UDP-[3H]GlcNAc (100 000 dpm/nmol), 0.5 mm O-mannosyl peptide [Ac-AAPT(Man)PVAAP-NH2], 10 mm MnCl2, 2% Triton X-100, 5 mm AMP, 200 mm GlcNAc, 10% glycerol and 2 μl culture supernatant at 37°C for 2 h. After boiling for 3 min, the mixture was analyzed by reversed phase high-pressure (or high performance) liquid chromatography with a Wakopak 5C18–200 column (4.6 × 250 mm). Solvent A was 0.085% trifluoroacetic acid in distilled water and solvent B was 0.085% trifluoroacetic acid in acetonitrile. The peptide was eluted at a flow rate of 1 ml/min using a linear gradient of 0–25% solvent B. The peptide separation was monitored continuously at 214 nm, and the radioactivity of each fraction was measured using a liquid scintillation counter.

POMGNT1 immunohistochemistry

Immunohistochemistry was carried out as described in previous literature (71). Briefly, unfixed eyes of adult WT mice were harvested and frozen in melting isopentane,16 μm cryo-sections were cut, post-fixed in 0.5% PFA for 10 min and treated with 0.1% Triton X-100 in PBS for 10 min and subsequently blocked in 10% normal goat serum with 0.1% Triton X-100 in PBS (blocking buffer) for 30 min to 1 h. Cryo-sections were then incubated for overnight at 4°C with primary antibodies (rabbit POMGNT1; 1:500; GeneTex) and anti-mouse acetylated α-tubulin (1:200; Santa-Cruz Biotechnology) diluted in blocking buffer. Slides were then washed in PBS, incubated with secondary antibody (Alexa 488 anti-rabbit; 1:500 dilution; Invitrogen) (Cy3 anti-mouse; 1:500; Jackson Immunochemicals) diluted in hybridization buffer at room temperature for 1 h, washed in PBS, DAPI (1:1000 dilution) was used for nuclear counter staining, slides were mounted with anti-fade medium (Prolong gold, Invitrogen) to reduce bleaching, and cover-slipped. Fluorescent images were captured with a Zeiss-Apotome.2 microscope, images were processed using ZEN and Adobe Photoshop CS4.

Funding

This work was supported by the National Institute of Health (1S10RR026550 to R.C.); National Eye Institute (R01EY022356 and R01EY018571 to R.C.); Foundation Fighting Blindness (BR-GE-0613-0618-BCM to R.C. and CD-CL-0214-0631-PUMCH to R.S.); Canadian Foundation for Fighting Blindness (FFB Canada) (to R.K.K); Foundation of Canadian Institutes for Health Research (CIHR) (to R.K.K.); the Reseau Vision, the Fondation de Recherche en Santee du Quebec (FRSQ) (to R.K.K.); the Montreal Children’s Hospital Foundation (MCHF) (to R.K.K.); National Natural Science Foundation of China (81470669 to R.S); Beijing Natural Science Foundation (7152116 to R.S.); Grants-in-Aid for Scientific Research (B) (25293016 to T.E.); Grants-in-Aid for Intramural Research grant (26-8 to T.E.); Grants-in-Aid for Neurological and Psychiatric Disorders of NCNP (to T.E.); Grants-in-Aid for Scientific Research on Innovative Areas (26110727 to H.M.); Mizutani Foundation for Glycoscience (150171 to H.M.) and Cullen Foundation endowment (to M.X.).

Acknowledgements

We thank all the patients and their family members for participating in this study. We thank the Exome Aggregation Consortium and the groups that provided exome variant data.

Conflict of Interest statement. None declared.

References

1

Ferrari
S.
,
Di Iorio
E.
,
Barbaro
V.
,
Ponzin
D.
,
Sorrentino
F.S.
,
Parmeggiani
F.
(
2011
)
Retinitis pigmentosa: genes and disease mechanisms
.
Curr. Genomics
,
12
,
238
249
.

2

Wang
F.
,
Wang
H.
,
Tuan
H.F.
,
Nguyen
D.H.
,
Sun
V.
,
Keser
V.
,
Bowne
S.J.
,
Sullivan
L.S.
,
Luo
H.
,
Zhao
L.
et al. . (
2014
)
Next generation sequencing-based molecular diagnosis of retinitis pigmentosa: identification of a novel genotype-phenotype correlation and clinical refinements
.
Hum. Genet.
,
133
,
331
345
.

3

Rachel
R.A.
,
Li
T.
,
Swaroop
A.
(
2012
)
Photoreceptor sensory cilia and ciliopathies: focus on CEP290, RPGR and their interacting proteins
.
Cilia
,
1
,
22
.

4

Haltia
M.
(
2006
)
The neuronal ceroid-lipofuscinoses: from past to present
.
Biochim. Biophys. Acta
,
1762
,
850
856
.

5

Hayflick
S.J.
(
2003
)
Unraveling the Hallervorden-Spatz syndrome: pantothenate kinase-associated neurodegeneration is the name
.
Curr. Opin. Pediatr.
,
15
,
572
577
.

6

Haas
D.
,
Hoffmann
G.F.
(
2006
)
Mevalonate kinase deficiencies: from mevalonic aciduria to hyperimmunoglobulinemia D syndrome
.
Orphanet J. Rare Dis.
,
1
,
13
.

7

Hildebrandt
F.
,
Benzing
T.
,
Katsanis
N.
(
2011
)
Ciliopathies
.
N. Eng. J. Med.
,
364
,
1533
1543
.

8

Houten
S.M.
,
Kuis
W.
,
Duran
M.
,
de Koning
T.J.
,
van Royen-Kerkhof
A.
,
Romeijn
G.J.
,
Frenkel
J.
,
Dorland
L.
,
de Barse
M.M.
,
Huijbers
W.A.
et al. . (
1999
)
Mutations in MVK, encoding mevalonate kinase, cause hyperimmunoglobulinaemia D and periodic fever syndrome
.
Nat. Genet.
,
22
,
175
177
.

9

Siemiatkowska
A.M.
,
van den Born
L.I.
,
van Hagen
P.M.
,
Stoffels
M.
,
Neveling
K.
,
Henkes
A.
,
Kipping-Geertsema
M.
,
Hoefsloot
L.H.
,
Hoyng
C.B.
,
Simon
A.
et al. . (
2013
)
Mutations in the mevalonate kinase (MVK) gene cause nonsyndromic retinitis pigmentosa
.
Ophthalmology
,
120
,
2697
2705
.

10

Joensuu
T.
,
Hamalainen
R.
,
Yuan
B.
,
Johnson
C.
,
Tegelberg
S.
,
Gasparini
P.
,
Zelante
L.
,
Pirvola
U.
,
Pakarinen
L.
,
Lehesjoki
A.E.
et al. . (
2001
)
Mutations in a novel gene with transmembrane domains underlie Usher syndrome type 3
.
Am. J. Hum. Genet.
,
69
,
673
684
.

11

Khan
M.I.
,
Kersten
F.F.
,
Azam
M.
,
Collin
R.W.
,
Hussain
A.
,
Shah
S.T.
,
Keunen
J.E.
,
Kremer
H.
,
Cremers
F.P.
,
Qamar
R.
et al. . (
2011
)
CLRN1 mutations cause nonsyndromic retinitis pigmentosa
.
Ophthalmology
,
118
,
1444
1448
.

12

Coussa
R.G.
,
Otto
E.A.
,
Gee
H.Y.
,
Arthurs
P.
,
Ren
H.
,
Lopez
I.
,
Keser
V.
,
Fu
Q.
,
Faingold
R.
,
Khan
A.
et al. . (
2013
)
WDR19: an ancient, retrograde, intraflagellar ciliary protein is mutated in autosomal recessive retinitis pigmentosa and in Senior-Loken syndrome
.
Clin. Genet.
,
84
,
150
159
.

13

Bujakowska
K.M.
,
Zhang
Q.
,
Siemiatkowska
A.M.
,
Liu
Q.
,
Place
E.
,
Falk
M.J.
,
Consugar
M.
,
Lancelot
M.E.
,
Antonio
A.
,
Lonjou
C.
et al. . (
2014
)
Mutations in IFT172 cause isolated retinal degeneration and Bardet-Biedl syndrome
.
Hum. Mol. Genet.
,
24
,
230
242
.

14

Xu
M.
,
Yang
L.
,
Wang
F.
,
Li
H.
,
Wang
X.
,
Wang
W.
,
Ge
Z.
,
Wang
K.
,
Zhao
L.
,
Li
H.
et al. . (
2015
)
Mutations in human IFT140 cause non-syndromic retinal degeneration
.
Hum. Genet.
,
134
,
1069
1078
.

15

Wang
F.
,
Li
H.
,
Xu
M.
,
Li
H.
,
Zhao
L.
,
Yang
L.
,
Zaneveld
J.E.
,
Wang
K.
,
Li
Y.
,
Sui
R.
et al. . (
2015
)
A homozygous missense mutation in NEUROD1 is associated with nonsyndromic autosomal recessive retinitis pigmentosa
.
Invest. Ophthalmol. Vis. Sci.
,
56
,
150
155
.

16

Haer-Wigman
L.
,
Newman
H.
,
Leibu
R.
,
Bax
N.M.
,
Baris
H.N.
,
Rizel
L.
,
Banin
E.
,
Massarweh
A.
,
Roosing
S.
,
Lefeber
D.J.
et al. . (
2015
)
Non-syndromic retinitis pigmentosa due to mutations in the mucopolysaccharidosis type IIIC gene, heparan-alpha-glucosaminide N-acetyltransferase (HGSNAT)
.
Hum. Mol. Genet.
,
24
,
3742
3751
.

17

Spiro
R.G.
(
2002
)
Protein glycosylation: nature, distribution, enzymatic formation, and disease implications of glycopeptide bonds
.
Glycobiology
,
12
,
43r
56r
.

18

Zelinger
L.
,
Banin
E.
,
Obolensky
A.
,
Mizrahi-Meissonnier
L.
,
Beryozkin
A.
,
Bandah-Rozenfeld
D.
,
Frenkel
S.
,
Ben-Yosef
T.
,
Merin
S.
,
Schwartz
S.B.
et al. . (
2011
)
A missense mutation in DHDDS, encoding dehydrodolichyl diphosphate synthase, is associated with autosomal-recessive retinitis pigmentosa in Ashkenazi Jews
.
Am. J. Hum. Genet.
,
88
,
207
215
.

19

Zuchner
S.
,
Dallman
J.
,
Wen
R.
,
Beecham
G.
,
Naj
A.
,
Farooq
A.
,
Kohli
M.A.
,
Whitehead
P.L.
,
Hulme
W.
,
Konidari
I.
et al. . (
2011
)
Whole-exome sequencing links a variant in DHDDS to retinitis pigmentosa
.
Am. J. Hum. Genet.
,
88
,
201
206
.

20

Chiba
A.
,
Matsumura
K.
,
Yamada
H.
,
Inazu
T.
,
Shimizu
T.
,
Kusunoki
S.
,
Kanazawa
I.
,
Kobata
A.
,
Endo
T.
(
1997
)
Structures of sialylated O-linked oligosaccharides of bovine peripheral nerve alpha-dystroglycan. The role of a novel O-mannosyl-type oligosaccharide in the binding of alpha-dystroglycan with laminin
.
J. Biol. Chem.
,
272
,
2156
2162
.

21

Sasaki
T.
,
Yamada
H.
,
Matsumura
K.
,
Shimizu
T.
,
Kobata
A.
,
Endo
T.
(
1998
)
Detection of O-mannosyl glycans in rabbit skeletal muscle alpha-dystroglycan
.
Biochim. Biophys. Acta
,
1425
,
599
606
.

22

Smalheiser
N.R.
,
Haslam
S.M.
,
Sutton-Smith
M.
,
Morris
H.R.
,
Dell
A.
(
1998
)
Structural analysis of sequences O-linked to mannose reveals a novel Lewis X structure in cranin (dystroglycan) purified from sheep brain
.
J. Biol. Chem.
,
273
,
23698
23703
.

23

Williamson
R.A.
,
Henry
M.D.
,
Daniels
K.J.
,
Hrstka
R.F.
,
Lee
J.C.
,
Sunada
Y.
,
Ibraghimov-Beskrovnaya
O.
,
Campbell
K.P.
(
1997
)
Dystroglycan is essential for early embryonic development: disruption of Reichert's membrane in Dag1-null mice
.
Hum. Mol. Genet.
,
6
,
831
841
.

24

Willer
T.
,
Prados
B.
,
Falcon-Perez
J.M.
,
Renner-Muller
I.
,
Przemeck
G.K.
,
Lommel
M.
,
Coloma
A.
,
Valero
M.C.
,
de Angelis
M.H.
,
Tanner
W.
et al. . (
2004
)
Targeted disruption of the Walker-Warburg syndrome gene Pomt1 in mouse results in embryonic lethality
.
Proc. Natl Acad. Sci. USA
,
101
,
14126
14131
.

25

Godfrey
C.
,
Foley
A.R.
,
Clement
E.
,
Muntoni
F.
(
2011
)
Dystroglycanopathies: coming into focus
.
Curr. Opin. Genet. Dev.
,
21
,
278
285
.

26

Takahashi
S.
,
Sasaki
T.
,
Manya
H.
,
Chiba
Y.
,
Yoshida
A.
,
Mizuno
M.
,
Ishida
H.
,
Ito
F.
,
Inazu
T.
,
Kotani
N.
et al. . (
2001
)
A new beta-1,2-N-acetylglucosaminyltransferase that may play a role in the biosynthesis of mammalian O-mannosyl glycans
.
Glycobiology
,
11
,
37
45
.

27

Yoshida
A.
,
Kobayashi
K.
,
Manya
H.
,
Taniguchi
K.
,
Kano
H.
,
Mizuno
M.
,
Inazu
T.
,
Mitsuhashi
H.
,
Takahashi
S.
,
Takeuchi
M.
et al. . (
2001
)
Muscular dystrophy and neuronal migration disorder caused by mutations in a glycosyltransferase, POMGnT1
.
Dev. Cell
,
1
,
717
724
.

28

Godfrey
C.
,
Clement
E.
,
Mein
R.
,
Brockington
M.
,
Smith
J.
,
Talim
B.
,
Straub
V.
,
Robb
S.
,
Quinlivan
R.
,
Feng
L.
et al. . (
2007
)
Refining genotype phenotype correlations in muscular dystrophies with defective glycosylation of dystroglycan
.
Brain
,
130
,
2725
2735
.

29

Yis
U.
,
Uyanik
G.
,
Rosendahl
D.M.
,
Carman
K.B.
,
Bayram
E.
,
Heise
M.
,
Comertpay
G.
,
Kurul
S.H.
(
2014
)
Clinical, radiological, and genetic survey of patients with muscle-eye-brain disease caused by mutations in POMGNT1
.
Pediatr. Neurol.
,
50
,
491
497
.

30

Pihko
H.
,
Lappi
M.
,
Raitta
C.
,
Sainio
K.
,
Valanne
L.
,
Somer
H.
,
Santavuori
P.
(
1995
)
Ocular findings in muscle-eye-brain (MEB) disease: a follow-up study
.
Brain Dev.
,
17
,
57
61
.

31

Akasaka-Manya
K.
,
Manya
H.
,
Kobayashi
K.
,
Toda
T.
,
Endo
T.
(
2004
)
Structure-function analysis of human protein O-linked mannose beta1,2-N-acetylglucosaminyltransferase 1, POMGnT1
.
Biochem. Biophys. Res. Commun.
,
320
,
39
44
.

32

Saredi
S.
,
Ardissone
A.
,
Ruggieri
A.
,
Mottarelli
E.
,
Farina
L.
,
Rinaldi
R.
,
Silvestri
E.
,
Gandioli
C.
,
D'Arrigo
S.
,
Salerno
F.
et al. . (
2012
)
Novel POMGNT1 point mutations and intragenic rearrangements associated with muscle-eye-brain disease
.
J. Neurol. Sci.
,
318
,
45
50
.

33

Barresi
R.
,
Campbell
K.P.
(
2006
)
Dystroglycan: from biosynthesis to pathogenesis of human disease
.
J. Cell Sci.
,
119
,
199
207
.

34

Lommel
M.
,
Strahl
S.
(
2009
)
Protein O-mannosylation: conserved from bacteria to humans
.
Glycobiology
,
19
,
816
828
.

35

Takahashi
H.
,
Kanesaki
H.
,
Igarashi
T.
,
Kameya
S.
,
Yamaki
K.
,
Mizota
A.
,
Kudo
A.
,
Miyagoe-Suzuki
Y.
,
Takeda
S.
,
Takahashi
H.
(
2011
)
Reactive gliosis of astrocytes and Muller glial cells in retina of POMGnT1-deficient mice
.
Mol. Cell Neurosci.
,
47
,
119
130
.

36

Liu
J.
,
Ball
S.L.
,
Yang
Y.
,
Mei
P.
,
Zhang
L.
,
Shi
H.
,
Kaminski
H.J.
,
Lemmon
V.P.
,
Hu
H.
(
2006
)
A genetic model for muscle-eye-brain disease in mice lacking protein O-mannose 1,2-N-acetylglucosaminyltransferase (POMGnT1)
.
Mech. Dev.
,
123
,
228
240
.

37

Mirza
M.
,
Volz
C.
,
Karlstetter
M.
,
Langiu
M.
,
Somogyi
A.
,
Ruonala
M.O.
,
Tamm
E.R.
,
Jagle
H.
,
Langmann
T.
(
2013
)
Progressive retinal degeneration and glial activation in the CLN6 (nclf) mouse model of neuronal ceroid lipofuscinosis: a beneficial effect of DHA and curcumin supplementation
.
PLoS One
,
8
,
e75963
.

38

Kanagawa
M.
,
Omori
Y.
,
Sato
S.
,
Kobayashi
K.
,
Miyagoe-Suzuki
Y.
,
Takeda
S.
,
Endo
T.
,
Furukawa
T.
,
Toda
T.
(
2010
)
Post-translational maturation of dystroglycan is necessary for pikachurin binding and ribbon synaptic localization
.
J. Biol. Chem.
,
285
,
31208
31216
.

39

Sato
S.
,
Omori
Y.
,
Katoh
K.
,
Kondo
M.
,
Kanagawa
M.
,
Miyata
K.
,
Funabiki
K.
,
Koyasu
T.
,
Kajimura
N.
,
Miyoshi
T.
et al. . (
2008
)
Pikachurin, a dystroglycan ligand, is essential for photoreceptor ribbon synapse formation
.
Nat. Neurosci.
,
11
,
923
931
.

40

Teber
S.
,
Sezer
T.
,
Kafali
M.
,
Manzini
M.C.
,
Konuk Yuksel
B.
,
Tekin
M.
,
Fitoz
S.
,
Walsh
C.A.
,
Deda
G.
(
2008
)
Severe muscle-eye-brain disease is associated with a homozygous mutation in the POMGnT1 gene
.
Eur. J. Paediatr. Neurol.
,
12
,
133
136
.

41

Clement
E.M.
,
Godfrey
C.
,
Tan
J.
,
Brockington
M.
,
Torelli
S.
,
Feng
L.
,
Brown
S.C.
,
Jimenez-Mallebrera
C.
,
Sewry
C.A.
,
Longman
C.
et al. . (
2008
)
Mild POMGnT1 mutations underlie a novel limb-girdle muscular dystrophy variant
.
Arch. Neurol.
,
65
,
137
141
.

42

Hehr
U.
,
Uyanik
G.
,
Gross
C.
,
Walter
M.C.
,
Bohring
A.
,
Cohen
M.
,
Oehl-Jaschkowitz
B.
,
Bird
L.M.
,
Shamdeen
G.M.
,
Bogdahn
U.
et al. . (
2007
)
Novel POMGnT1 mutations define broader phenotypic spectrum of muscle-eye-brain disease
.
Neurogenetics
,
8
,
279
288
.

43

Taniguchi
K.
,
Kobayashi
K.
,
Saito
K.
,
Yamanouchi
H.
,
Ohnuma
A.
,
Hayashi
Y.K.
,
Manya
H.
,
Jin
D.K.
,
Lee
M.
,
Parano
E.
et al. . (
2003
)
Worldwide distribution and broader clinical spectrum of muscle-eye-brain disease
.
Hum. Mol. Genet.
,
12
,
527
534
.

44

Manya
H.
,
Sakai
K.
,
Kobayashi
K.
,
Taniguchi
K.
,
Kawakita
M.
,
Toda
T.
,
Endo
T.
(
2003
)
Loss-of-function of an N-acetylglucosaminyltransferase, POMGnT1, in muscle-eye-brain disease
.
Biochem. Biophys. Res. Commun.
,
306
,
93
97
.

45

Voglmeir
J.
,
Kaloo
S.
,
Laurent
N.
,
Meloni
M.M.
,
Bohlmann
L.
,
Wilson
I.B.
,
Flitsch
S.L.
(
2011
)
Biochemical correlation of activity of the alpha-dystroglycan-modifying glycosyltransferase POMGnT1 with mutations in muscle-eye-brain disease
.
Biochem. J.
,
436
,
447
455
.

46

Bunge
S.
,
Clements
P.R.
,
Byers
S.
,
Kleijer
W.J.
,
Brooks
D.A.
,
Hopwood
J.J.
(
1998
)
Genotype-phenotype correlations in mucopolysaccharidosis type I using enzyme kinetics, immunoquantification and in vitro turnover studies
.
Biochim. Biophys. Acta
,
1407
,
249
256
.

47

Fu
R.
,
Jinnah
H.A.
(
2012
)
Genotype-phenotype correlations in Lesch-Nyhan disease: moving beyond the gene
.
J. Biol. Chem.
,
287
,
2997
3008
.

48

Zurutuza
L.
,
Muller
F.
,
Gibrat
J.F.
,
Taillandier
A.
,
Simon-Bouy
B.
,
Serre
J.L.
,
Mornet
E.
(
1999
)
Correlations of genotype and phenotype in hypophosphatasia
.
Hum. Mol. Genet.
,
8
,
1039
1046
.

49

Wright
A.F.
,
Chakarova
C.F.
,
Abd El-Aziz
M.M.
,
Bhattacharya
S.S.
(
2010
)
Photoreceptor degeneration: genetic and mechanistic dissection of a complex trait
.
Nat. Rev. Genet.
,
11
,
273
284
.

50

Roosing
S.
,
van den Born
L.I.
,
Sangermano
R.
,
Banfi
S.
,
Koenekoop
R.K.
,
Zonneveld-Vrieling
M.N.
,
Klaver
C.C.
,
van Lith-Verhoeven
J.J.
,
Cremers
F.P.
,
den Hollander
A.I.
et al. . (
2015
)
Mutations in MFSD8, encoding a lysosomal membrane protein, are associated with nonsyndromic autosomal recessive macular dystrophy
.
Ophthalmology
,
122
,
170
179
.

51

Estrada-Cuzcano
A.
,
Roepman
R.
,
Cremers
F.P.
,
den Hollander
A.I.
,
Mans
D.A.
(
2012
)
Non-syndromic retinal ciliopathies: translating gene discovery into therapy
.
Hum. Mol. Genet.
,
21
,
R111
R124
.

52

Tobin
J.L.
,
Beales
P.L.
(
2009
)
The nonmotile ciliopathies
.
Genet. Med.
,
11
,
386
402
.

53

Diesen
C.
,
Saarinen
A.
,
Pihko
H.
,
Rosenlew
C.
,
Cormand
B.
,
Dobyns
W.B.
,
Dieguez
J.
,
Valanne
L.
,
Joensuu
T.
,
Lehesjoki
A.E.
(
2004
)
POMGnT1 mutation and phenotypic spectrum in muscle-eye-brain disease
.
J. Med. Genet.
,
41
,
e115
.

54

Mercuri
E.
,
Messina
S.
,
Bruno
C.
,
Mora
M.
,
Pegoraro
E.
,
Comi
G.P.
,
D'Amico
A.
,
Aiello
C.
,
Biancheri
R.
,
Berardinelli
A.
et al. . (
2009
)
Congenital muscular dystrophies with defective glycosylation of dystroglycan: a population study
.
Neurology
,
72
,
1802
1809
.

55

Cormand
B.
,
Pihko
H.
,
Bayes
M.
,
Valanne
L.
,
Santavuori
P.
,
Talim
B.
,
Gershoni-Baruch
R.
,
Ahmad
A.
,
van Bokhoven
H.
,
Brunner
H.G.
et al. . (
2001
)
Clinical and genetic distinction between Walker-Warburg syndrome and muscle-eye-brain disease
.
Neurology
,
56
,
1059
1069
.

56

Abadi
R.V.
(
2002
)
Mechanisms underlying nystagmus
.
J. R. Soc. Med.
,
95
,
231
234
.

57

Dobbs
M.D.
,
Mawn
L.A.
,
Donahue
S.P.
(
2011
)
Anomalous extraocular muscles with strabismus
.
AJNR Am. J. Neuroradiol.
,
32
,
E167
E168
.

58

Ghasia
F.F.
,
Shaikh
A.G.
(
2013
)
Pattern strabismus: where does the brain's role end and the muscle's begin?
J. Ophthalmol.
,
2013
,
301256
.

59

Sowden
J.C.
(
2007
)
Molecular and developmental mechanisms of anterior segment dysgenesis
.
Eye (Lond)
,
21
,
1310
1318
.

60

Verma
A.S.
,
Fitzpatrick
D.R.
(
2007
)
Anophthalmia and microphthalmia
.
Orphanet J. Rare Dis.
,
2
,
47
.

61

Li
H.
,
Durbin
R.
(
2009
)
Fast and accurate short read alignment with Burrows-Wheeler transform
.
Bioinformatics
,
25
,
1754
1760
.

62

Ng
P.C.
,
Henikoff
S.
(
2003
)
SIFT: predicting amino acid changes that affect protein function
.
Nucleic Acids Res.
,
31
,
3812
3814
.

63

Adzhubei
I.A.
,
Schmidt
S.
,
Peshkin
L.
,
Ramensky
V.E.
,
Gerasimova
A.
,
Bork
P.
,
Kondrashov
A.S.
,
Sunyaev
S.R.
(
2010
)
A method and server for predicting damaging missense mutations
.
Nat. Methods
,
7
,
248
249
.

64

Chun
S.
,
Fay
J.C.
(
2009
)
Identification of deleterious mutations within three human genomes
.
Genome Res.
,
19
,
1553
1561
.

65

Schwarz
J.M.
,
Rodelsperger
C.
,
Schuelke
M.
,
Seelow
D.
(
2010
)
MutationTaster evaluates disease-causing potential of sequence alterations
.
Nat. Methods
,
7
,
575
576
.

66

Reva
B.
,
Antipin
Y.
,
Sander
C.
(
2011
)
Predicting the functional impact of protein mutations: application to cancer genomics
.
Nucleic Acids Res.
,
39
,
e118
.

67

Shihab
H.A.
,
Gough
J.
,
Cooper
D.N.
,
Stenson
P.D.
,
Barker
G.L.
,
Edwards
K.J.
,
Day
I.N.
,
Gaunt
T.R.
(
2013
)
Predicting the functional, molecular, and phenotypic consequences of amino acid substitutions using hidden Markov models
.
Hum. Mutat.
,
34
,
57
65
.

68

Salvo
J.
,
Lyubasyuk
V.
,
Xu
M.
,
Wang
H.
,
Wang
F.
,
Nguyen
D.
,
Wang
K.
,
Luo
H.
,
Wen
C.
,
Shi
C.
et al. . (
2015
)
Next-generation sequencing and novel variant determination in a cohort of 92 familial exudative vitreoretinopathy patients
.
Invest. Ophthalmol. Vis. Sci.
,
56
,
1937
1946
.

69

Untergasser
A.
,
Cutcutache
I.
,
Koressaar
T.
,
Ye
J.
,
Faircloth
B.C.
,
Remm
M.
,
Rozen
S.G.
(
2012
)
Primer3--new capabilities and interfaces
.
Nucleic Acids Res.
,
40
,
e115
.

70

Xin
X.
,
Akasaka-Manya
K.
,
Manya
H.
,
Furukawa
J.
,
Kuwahara
N.
,
Okada
K.
,
Tsumoto
H.
,
Higashi
N.
,
Kato
R.
,
Shinohara
Y.
et al. . (
2015
)
POMGNT1 is glycosylated by mucin-type O-glycans
.
Biol. Pharm. Bull.
,
38
,
1389
1394
.

71

Eblimit
A.
,
Nguyen
T.M.
,
Chen
Y.
,
Esteve-Rudd
J.
,
Zhong
H.
,
Letteboer
S.
,
Van Reeuwijk
J.
,
Simons
D.L.
,
Ding
Q.
,
Wu
K.M.
et al. . (
2015
)
Spata7 is a retinal ciliopathy gene critical for correct RPGRIP1 localization and protein trafficking in the retina
.
Hum. Mol. Genet.
,
24
,
1584
1601
.

Author notes

The authors wish it to be known that, in their opinion, the last four authors should be regarded as joint corresponding authors.

Supplementary data