Disruption of PIKFYVE causes congenital cataract in human and zebrafish
Abstract
Congenital cataract, an ocular disease predominantly occurring within the first decade of life, is one of the leading causes of blindness in children. However, the molecular mechanisms underlying the pathogenesis of congenital cataract remain incompletely defined. Through whole-exome sequencing of a Chinese family with congenital cataract, we identified a potential pathological variant (p.G1943E) in PIKFYVE, which is located in the PIP kinase domain of the PIKFYVE protein. We demonstrated that heterozygous/homozygous disruption of PIKFYVE kinase domain, instead of overexpression of PIKFYVEG1943E in zebrafish mimicked the cataract defect in human patients, suggesting that haploinsufficiency, rather than dominant-negative inhibition of PIKFYVE activity caused the disease. Phenotypical analysis of pikfyve zebrafish mutants revealed that loss of Pikfyve caused aberrant vacuolation (accumulation of Rab7+Lc3+ amphisomes) in lens cells, which was significantly alleviated by treatment with the V-ATPase inhibitor bafilomycin A1 (Baf-A1). Collectively, we identified PIKFYVE as a novel causative gene for congenital cataract and pinpointed the potential application of Baf-A1 for the treatment of congenital cataract caused by PIKFYVE deficiency.
Editor's evaluation
This manuscript will be of interest to readers in the field of eye development and pathology and ocular geneticists. Using a zebrafish model, the authors have identified a new autosomal dominant mutation in the kinase domain of the human phosphoinositide kinase gene PIKFYVE and assessed its impact on congenital cataract (clouding of the lens). They demonstrated that heterozygous or homozygous disruption of the PIKFYVE kinase domain, but not the overexpression of PIKFYVE G1943E, led to a cataract defect, suggesting the mechanism of haploinsufficiency rather than dominant-negative inhibition of PIKFYVE activity. The authors have presented strong evidence that mutations in the PIKFYVE kinase domain could lead to congenital cataract. Overall, the same pipeline of in vitro and in vivo genetic approaches could be applied to other human genetic disorders.
https://doi.org/10.7554/eLife.71256.sa0Introduction
Congenital cataract is partial or complete opacification of the lens that occurs at birth or during the first decade of life. It is a common ocular abnormality that causes visual impairment and blindness during infancy (Khokhar et al., 2017), accounting for 5–20% blindness in children worldwide (Sheeladevi et al., 2016). According to the location and shape of the lens opacities, congenital cataract can be divided into seven clinical types: nuclear cataract, polar cataract, lamellar cataract, nuclear with cortical cataract, cortical cataract, sutural cataract, and total cataract (Zhai et al., 2017). Different types of cataracts cause different levels of visual impairment in patients. A variety of factors, including gene variants that affect lens metabolism. are closely associated with cataract (Li et al., 2020), and autosomal dominant congenital cataract is the most common mode of inheritance (Berry et al., 2020a). To date, the occurrence of congenital cataract has been linked with genetic variants in at least 50 genes involved in lens structure and development (Berry et al., 2020b; Shiels and Hejtmancik, 2017), including crystallin genes (CRYAA, CRYAB, CRYBB1, CRYBB2, CRYBB3, CRYBA1/A3, CRYBA2, CRYBA4, CRYGC, CRYGD, and CRYGS) (Bhat, 2003; Zhuang et al., 2019), membrane protein genes (GJA3, GJA8, MIP, and LIM2) (Berry et al., 2000; Beyer et al., 2013; Pei et al., 2020), growth and transcription factor genes (PITX3, MAF, and HSF4) (Anand et al., 2018), beaded filament structural protein genes (BFSP1 and BFSP2) (Song et al., 2009), and other genes (CHMP4B and EPHA2) (Dave et al., 2016; Shiels et al., 2007). Although these findings have made tremendous contributions to our understanding of the genetic etiology of congenital cataract, new causative genes as well as the underlying molecular mechanisms remain to be discovered (Berry et al., 2020a).
In this study, we identified a missense genetic variant (p.G1943E) in the phosphatidylinositol phosphate kinase (PIPK) domain of PIKFYVE (phosphoinositide kinase, FYVE-type zinc finger containing) from a Chinese cataract family. The human PIKFYVE gene encodes a large protein consisting of 2098 amino acids, which contains six evolutionarily conserved domains (Kawasaki et al., 2012), including zinc finger phosphoinositide kinase (FYVE) domain, plextrin homology domain, β-sheet winged helix DNA/RNA-binding motif, cytosolic chaperone CCTγ apical domain-like motif, spectrin repeats (SPEC), and C-terminal fragment PIPK domain (Li et al., 2005; Shisheva, 2008). The N-terminal FYVE domain can bind to PtdIns3P on the membrane, while the C-terminal is the kinase domain which has the catalytic function for the production of PtdIns(3,5)P2 (Shisheva et al., 1999). The in-vivo function of PIKFYVE in the development of lens and its association with congenital cataract have not been documented. Different from previously reported PIKFYVE variants affecting the chaperonin-like domain that cause fleck corneal dystrophy (CFD), a disease characterized by tiny white flecks scattered in corneal stroma (Gee et al., 2015; Kawasaki et al., 2012; Li et al., 2005), the p.G1943E variant in the kinase domain of PIKFYVE was revealed for the first time to cause congenital cataract in human patients. Using zebrafish models, we demonstrated that disruption of the PIPK domain in Pikfyve caused early-onset cataract defect (aberrant accumulation of Rab7+Lc3+ amphisomes in the lens) and treating the pikfyve-deficient embryos with V-ATPase inhibitor bafilomycin A1 (Baf-A1) significantly alleviated the vacuolation defect in the lens of zebrafish mutants.
Results
The congenital cataract family and clinical characteristics
We recruited a large Chinese Korean family affected with congenital cataract, which consists of 31 family members including 12 affected individuals across four generations (Figure 1A). All the family members were Chinese Korean except one married-in individual (III-8) who was Han Chinese. There was no consanguineous marriage in this family. Bilateral cataract was observed in every generation of the pedigree and affected parents transmitted the disease to both males and females, suggesting an autosomal dominant inheritance. No other ocular and systemic abnormalities were found in these patients. Clinical information was obtained from all patients except two deceased individuals (Table 1). Seven of the 10 living patients were diagnosed as cataract in their first decade of life, with the smallest age of diagnosis at 2 years old. The best corrected visual acuity (BCVA) of these patients ranged from 20/125 to 20/20. The proband (III-9) showed nuclear pulverulent cataract in both eyes (Figure 1B) with vision loss (OD: 20/40, OS: 20/66) since childhood. The mother (II-4) of the proband also had nuclear pulverulent cataract in both eyes, while his father (II-3) had no cataract in either eye. The son (IV-5) of the proband had Y-sutural cataract. Among the other affected members of this family, individuals II-8, III-5, III-6, and IV-3 presented nuclear pulverulent opacity in both eyes, whereas individual III-1 and his son (IV-1) showed peripheral cortical punctate opacity in both eyes. All patients of this family did not show any abnormities in fundus photography and optical coherence tomography (OCT) examinations.
Identification of a pathogenic genetic variant in PIKFYVE in the cataract family
To identify the pathogenic variant responsible for congenital cataract in this family, we performed whole-exome sequencing (WES) analysis of four selected family members including three affected members (II-1, III-9, and IV-5) and one unaffected member (II-5; Figure 1A). Overall, 101,852–102,401 sequence variants were identified in each individual. Nonsynonymous or splicing-site variants, of which the minor allele frequencies were less than 1% in the genome aggregate database (gnomAD; https://gnomad.broadinstitute.org/) and were shared among all three affected members but were absent from the unaffected member were filtered and retained. Then we focused on variants in known cataract genes but did not find any disease causative variants in this family. Further filtering based on in-silico prediction of pathogenicity refined the candidates to three genes (PIKFYVE, NPHS1, and FPR1; Table 2). We then verified and evaluated these candidate variants in the whole family using Sanger sequencing. Among them, only the PIKFYVE variant completely segregated with the disease in the family (Figure 1A). As shown in Figure 2A, in contrast to the healthy control, the PIKFYVE gene in the cataract patients carried a G to A heterozygous substitution in the 39th exon, which leads to the replacement of a highly conserved amino acid glycine (G) by glutamate (E) at position 1943 in the PIPK domain (Figure 2B and C). To investigate the possible roles of PIKFYVE in the development of cataracts, we first confirmed its expression in human anterior lens capsules by quantitative RT-PCR (Figure 2—figure supplement 1A and B). Then to interpret how p.G1943E variant altered the function of PIKFYVE, we determined the protein stability of PIKFYVEWT and PIKFYVEG1943E by Western blot analysis of pCS2(+)-CMV-PIKFYVEWT and pCS2(+)-CMV-PIKFYVEG1943E transfected HEK293T cells at 34 and 44 hr post-transfection. As shown in Figure 2D, the expression level of PIKFYVEG1943E was comparable with that of PIKFYVEWT, suggesting that the p.G1943E variant has little or no effect on the protein stability of PIKFYVE. Recently, a middle-to-low resolution structure of PIKFYVE was resolved by cryo-EM (PDB ID:7K2V) (Lees et al., 2020). Using this structure as template, we predicted the structure of the p.G1943E variant form of PIKFYVE PIPK domain. As shown in Figure 2E and F, the p.G1943E sits in a linker connecting the N-lobe (gold) and the C-lobe (cyan) of PIPK domain. This loop (red) acts as a hinge between N-lobe and C-lobe (Figure 2F) and forms a highly restricted conformation. However, the p.G1943E variant introduces a negatively charged residue to the loop. This negatively charged residue not only changes the surface electrostatic potential of PIPK domain (Figure 2G), but also generates charge repulsions between the loop and N-lobe (D1872 and D1871) (Figure 2F), thereby potentially affecting the kinase activity of PIKFYVE.
To further strengthen the correlation of PIKFYVE variants with cataracts, we screened another congenital cataract family, 10 sporadic cases with congenital cataract and 200 patients with age-related or complicated cataract. In this experiment, six affected loci from seven age-related or complicated cataract patients were identified. Among them, two variants in exon 36 (c.5399A>C, p.Q1800P) and exon 37 (c.5594C>T, p.A1865V) are located in the PIP kinase domain (Figure 2—figure supplement 2, Figure 2—figure supplement 2—source data 1). These findings further suggested that genetic variants in PIKFYVE might be causal/risk factors for cataract formation.
Disruption of PIPK domain of PIKFYVE caused early onset cataract in zebrafish
Bioinformatic analysis suggests that the p.G1943E variant introduces a negatively charged residue in the loop, which potentially affects the kinase activity of PIKFYVE. To directly demonstrate that heterozygosity of the p.G1943E variant caused haploinsufficiency rather than dominant-negative inhibition of PIKFYVE activity in human patients, we generated pikfyve-deficient zebrafish mutants using the CRISPR/Cas9-directed gene editing method (Shankaran et al., 2018). To maximally mimic the p.G1943E variant in human patients, we designed single-guide RNA (sgRNA) in exon 40 to target the C terminal PIPK domain and screened out the pikfyveΔ8 allele, which harbored an 8 bp deletion and a 112 bp insertion in the cDNA of pikfyve. As shown in Figure 3A, this insertion/deletion in pikfyveΔ8 allele introduced a premature stop codon immediately downstream of the insertion site, thereby presumably generating a truncated Pikfyve protein that lacks the kinase activity. Phenotypical analysis revealed that pikfyveΔ8 homozygous mutants developed normally before 5 days post-fertilization (dpf), but later showed severe developmental defects and all died by 9 dpf (Figure 3—figure supplement 1A, B). Under stereomicroscope, we noticed that compared to siblings (sib), the lens of pikfyveΔ8 mutants was less transparent at 5 dpf (Figure 3B). To characterize this phenotype in details, we followed up the development of lens in all genotypes at different stages under differential interference contrast (DIC) microscope. Our results showed that bubble-like vacuoles first appeared in the developing lens of pikfyveΔ8 homozygous mutants at 3 dpf (Figure 3C and D). Strikingly, the number and size of vacuoles both increased drastically by 5 dpf and almost dominated the lens of mutants. Notably, compared to the wild-type (WT) control, pikfyveΔ8 heterozygous mutants also manifested significantly higher number of vacuoles in their lens (Figure 3C and D), suggesting that normal development of lens highly depends on the activity of Pikfyve and heterozygous disruption of Pikfyve is sufficient to cause cataract phenotype. Furthermore, pharmacological inhibition of PIKFYVE activity by inhibitor YM201636 also induced aberrant vacuolation in human lens epithelial cells HLEB3 (Figure 3E and F), suggesting that the function of PIKFYVE in the regulation of cataract formation is highly conserved across species and our zebrafish cataract model would potentially be an excellent system for further preclinical study.
In addition to loss of function study of Pikfyve, we also utilized a ubiquitously expressed ubiquitin (ubi) promoter to generate transgenic lines, that is, Tg(ubi:PIKFYVEWT) and Tg(ubi:PIKFYVEG1943E) to overexpress PIKFYVEWT and PIKFYVEG1943E in zebrafish, respectively (Figure 3—figure supplement 2A). As expected, both of these two zebrafish lines showed normal development of the lens during early stages (Figure 3—figure supplement 2B), suggesting that indeed haploinsufficiency, rather than dominant-negative inhibition of PIKFYVE activity by the p.G1943E variant caused cataract in human patients. Moreover, the vacuolation phenotype could be partially rescued in both pikfyveΔ8;Tg(ubi:PIKFYVEWT) and pikfyveΔ8;Tg(ubi:PIKFYVEG1943E) zebrafish (Figure 3—figure supplement 3A). However, while the vacuole numbers decreased to less than 15 in all pikfyveΔ8;Tg(ubi:PIKFYVEWT) embryos (n=6/6), about one-third of pikfyveΔ8;Tg(ubi:PIKFYVEG1943E) mutants still harbored more than 15 vacuoles (n=5/17) (Figure 3—figure supplement 3A). As whole-mount in situ hybridization (WISH) results revealed that the expression of PIKFYVE in Tg(ubi:PIKFYVEG1943E) was at least as strong as or even stronger than that of PIKFYVEWT (Figure 3—figure supplement 3B). These results indicated that the rescue effect of PIKFYVEG1943E is less efficient than that of PIKFYVEWT.
Detailed characterization of cataract phenotypes in PikfyveΔ8 mutants
To delineate the details of cataract phenotype in pikfyveΔ8 mutants, we crossed pikfyveΔ8 mutant with the reporter transgenic line Tg(cryaa:DsRed) (Nguyen-Chi et al., 2014), in which expression of DsRed was controlled by promoter of the crystalline gene cryaa, the major constitutive components of lens fibers (Runkle et al., 2002). Intriguingly, in comparison with evenly distributed DsRed signals in the lens of siblings, vacuoles in pikfyveΔ8 mutants almost occupied the surface of the lens and all of them were DsRed negative (Figure 4A), indicating that these vacuoles did not contain lens fibers. Meanwhile, hematoxylin-eosin (HE) staining, together with ZL-1 antibody and DAPI co-staining revealed that while disruption of Pikfyve seemingly had no effects on the enucleation process of lens, lens fibers in pikfyveΔ8 mutants were less organized than those in siblings (Figure 4B and C). To further characterize the vacuoles and lens structure in high resolution, we utilized transmission electron microscope (TEM) to visualize the ultrastructure of lens in both siblings and pikfyveΔ8 mutants. We could detect large vacuoles in the mutant lens at 3 dpf and 5 dpf, while only several tiny vacuoles were observed in the sibling lens (Figure 4D). Lens fiber cells of siblings were mildly edema, and their nuclei were oval with uniform chromatin. The lens fibers (white arrows) were arranged neatly and tightly. Furthermore, mitochondria (black triangles) of siblings were slightly swollen without vacuoles. In contrast, lens fiber cells of the pikfyveΔ8 mutants were obviously edematous and the nuclei were irregularly shaped. Besides, the arrangement of lens fibers was loose and deformed and lipid droplets were formed in pikfyveΔ8 mutants. Compared to the siblings, mitochondria (black triangles) of pikfyveΔ8 mutants were obviously swollen and enlarged. Mitochondrial crests were reduced, and most of them were aberrantly vacuolated. In addition, autophagic lysosomes also appeared swollen and vacuolated in pikfyveΔ8 mutants.
Vacuoles in PikfyveΔ8 mutants were amphisomes
To define the nature of vacuoles in the lens of pikfyveΔ8 mutants, we conducted time-lapse imaging to monitor their behaviors during zebrafish development. Our results showed that vacuoles in pikfyveΔ8 mutants were highly dynamic and small vacuoles were frequently found to be fused with each other to form larger ones (Figure 5A, white arrows). This feature, together with previous findings showing that PIKFYVE was an essential regulator of endomembrane homeostasis (Hasegawa et al., 2017), prompted us to further investigate whether vacuoles in pikfyveΔ8 mutants were virtually endocytic vesicles. To probe this issue, we generated fusion mRNAs encoding GFP and the markers of endocytic vesicles (i.e., the small GTPase Rab5c, Rab7, and Rab11a) (Langemeyer et al., 2018) and injected them into zebrafish embryos to specifically label early, late, and recycling endosomes in pikfyveΔ8 mutants. Our results showed that in pikfyveΔ8 mutants, almost all vacuoles were positive for the late endosome maker Rab7-GFP (Figure 5C, white arrows). By contrast, the early endosome marker Rab5c-GFP and recycling marker Rab11a-GFP showed no co-localization with vacuoles (Figure 5B-D). As autophagy was also closely related to the organelle membrane system and had been shown to be regulated by PIKFYVE (Vicinanza et al., 2015), we therefore checked the status of autophagosomes in pikfyveΔ8 mutants by injecting lc3b-mcherry fusion mRNA. As shown in Figure 5 and a proportion of vacuoles in pikfyveΔ8 mutants were also positive for the autophagosome marker Lc3b. Taken together, these data implied that vacuoles in pikfyveΔ8 mutants were amphisomes formed by the fusion of autophagosomes and late endosomes.
Baf-A1 partially rescued the vacuole defect in the lens of PikfyveΔ8 mutant zebrafish
Baf-A1 is a specific inhibitor of V-ATPase. Previous work has shown that the vacuole phenotype induced by PIKFYVE deficiency in COS-7 cells could be rescued by Baf-A1 (Compton et al., 2016). We thus investigated whether the vacuolation defect in the lens of pikfyveΔ8 mutant zebrafish could also be attenuated by Baf-A1. Indeed, we found that vacuole number in the lens of Baf-A1-treated pikfyveΔ8 mutants was significantly lower than that in the control group treated with dimethyl sulfoxide (DMSO) (Figure 6A and B). To further validate that Baf-A1 could directly alleviate cataract defect, rather than just delay the phenotype, we imaged the lens of the same mutant zebrafish before and after Baf-A1 treatment. As shown in Figure 6C and D, while vacuole numbers of all the mutant embryos showed a slight increase after DMSO treatment or without treatment, we observed significant decrease of vacuole numbers in all pikfyveΔ8 mutants. On the other hand, a previous study also showed overexpression of transient receptor potential mucolipin 1 (TRPML1) could partially rescue the vacuole phenotype in PIKfyve-deficient macrophages (Krishna et al., 2016). However, in our system, ectopic expression of Trpml1 by mRNA injection failed to rescue vacuolation phenotype in pikfyveΔ8 mutants (Figure 6—figure supplement 1). Thus, V-ATPase rather than Trpml1 might be a key component for the formation or maintenance of the large vacuoles in the lens of pikfyve-deficient mutants, and inhibition of V-ATPase by Baf-A1 could directly alleviate the cataract phenotype.
Discussion
In this study, we identified a missense variant (p.G1943E) in PIKFYVE responsible for congenital cataract in a Chinese Korean family (Figures 1 and 2). This variant is very rare in the general population. In the gnomAD database, only four individuals (including three East Asians and one Latino/Admixed American) carry the p.G1943E variant, representing an extremely low allele frequency of 0.00002 (Table 2). Consistent with previous findings (Berry et al., 2020b), the degree and morphology of lens opacity in the cataract family were phenotypically heterogeneous (Table 1). While most of the patients in this family developed cataract and vision loss in their childhood, some cases were diagnosed quite late. Meanwhile, the majority of patients had nuclear pulverulent cataract, while the others had nuclear Y-sutural cataract or peripheral cortical punctate cataract (Table 1). This heterogeneity in clinical manifestations might be ascribed to interactions between genetic and environmental factors during lens development (Berry et al., 2020b).
Variants in PIKFYVE have been reported to be associated with CFD (Gee et al., 2015; Kawasaki et al., 2012; Kotoulas et al., 2011; Li et al., 2005). However, except that two of the CFD patients showed cataract formation at 66 and 58 years old, respectively, none of these PIKFYVE variants in these studies caused congenital cataract (Kotoulas et al., 2011). Interestingly, we also did not observe any corneal defect in patients with congenital cataract in this study. It is noted that all CFD-related mutations in PIKFYVE, either homozygous or heterozygous, are distributed in two regions (i.e., amino acids 667–843 and 1490–1538), corresponding to cytosolic chaperone CCTγ apical domain-like motif and SPEC domains, respectively. By contrast, the p.G1943E variant identified in this study is located in the C-terminal PIPK domain. Based on these findings, we hypothesized that different domains of PIKFYVE might exert different functions by binding with different partners. This hypothesis is further supported by other evidences. First, the predicted protein structure suggested that, instead of gross disruption of PIKFYVE protein structure, the p.G1943E variant only mildly altered the 3-D conformation of PIPK domain by changing its surface electrostatic potential and generating charge repulsions between the loop and N-lobe (Figure 2G and F). Therefore, it is reasonable to argue that this mutation in PIPK domain may specifically alter its kinase activity or change the binding surface with other proteins, while keep the N-terminal cytosolic chaperone CCTγ apical domain-like motif and SPEC domains functional intact. Second, in our pikfyveΔ8 zebrafish mutants, while the C-terminal PIPK domain was truncated in PikfyveΔ8 protein, the N-terminal cytosolic chaperone CCTγ apical domain-like motif and SPEC domains were not affected. Accordingly, we did not observe any cornea defect in the zebrafish mutants either.
Inhibition of PIKFYVE in COS-7 cells had been shown to induce large vacuoles through promoting the enlargement of both early and late endosomes (Ikonomov et al., 2006; Rutherford et al., 2006). Further studies revealed that the Ca2+ releasing channel endolysosome-localized mucolipin TRPML1 acts downstream of PIKFYVE to trigger membrane fusion/fission process (Dong et al., 2010) or promote lysosome/phagosome maturation (Dayam et al., 2015; Kim et al., 2014). Interestingly, while overexpression of TRPML1 partially alleviated the vacuole phenotype in PIKFYVE-deficient macrophages (Krishna et al., 2016), the rescue effect was not observed in the lens of pikfyveΔ8 zebrafish mutants (Figure 6—figure supplement 1), suggesting that PIKFYVE might function in a context-dependent manner. Consistent with this idea, we also noted that vacuoles in microglia (the brain-resident macrophages) and lens cells of pikfyveΔ8 mutants were differently stained by lysosome marker (Figure 5—figure supplement 1). On the other hand, vacuole formation in PIKFYVE-deficient COS-7 cells and macrophages could also be inhibited by a drug called Baf-A1 (Compton et al., 2016; Isobe et al., 2019). In this study, we found that Baf-A1 could partially rescue the vacuole defect in the lens of pikfyveΔ8 mutant zebrafish (Figure 6A–D). Baf-A1 is a macrolide antibiotic that inhibits V-ATPase, the ATP-dependent proton pump located on the membrane of organelles. The cellular acidification process mediated by V-ATPase may affect many basic biological processes, including membrane trafficking (in particular endosome maturation and fusion between autophagosomes and lysosomes; Hammond et al., 1998; Yamamoto et al., 1998), protein degradation, and autophagy (Bowman et al., 1988; Yamamoto et al., 1998; Yoshimori et al., 1991). Interestingly, several works in yeast also identified mutations in V-ATPase that did not affect proton pump function, but indeed caused defects in vacuole fusion (Strasser et al., 2011). Moreover, another study in Drosophila demonstrated that inhibition of vesicle fusion by Baf-A1 did not depend on V-ATPases, but relied on Ca2+ sarco/endoplasmic reticulum Ca2+-ATPase (SERCA) pump, the secondary target of Baf-A1 (Mauvezin et al., 2015). Therefore, the target of Baf-A1 and underlying mechanisms still need to be further investigated. These mechanisms will finally contribute to the potential clinical application of the drug in the treatment of congenital cataract caused by PIKFYVE deficiency.
Materials and methods
Pedigree and patients
Request a detailed protocolA four-generation pedigree consisting of 31 family members (Figure 1A) was recruited from the eye clinic of Shenzhen Eye Hospital, Shenzhen, Guangdong, China. All the living family members underwent a complete ocular examination. Twelve individuals in this family were diagnosed with congenital cataract (Figure 1B, Table 1), while the other 19 individuals were unaffected. The diagnosis of congenital cataract was based on (1) lens opacity detected at birth or during the first decade of life; and (2) no known other causes (e.g., trauma, iatrogenic, or inflammatory disease) (Berry et al., 2020a). Meanwhile, we also recruited another congenital cataract family, 10 sporadic cases with congenital cataract and 200 patients with age-related cataract or complicated cataract. The ethical protocol in this study was approved by the Independent Ethics Committee of Shenzhen Eye Hospital, in accordance with the tenets of the Declaration of Helsinki (1975). Written informed consent and consent to publish were received from all study participants prior to obtaining the human tissue samples. Samples were de-identified by removing the identifiable individual information such as name and identity number before researching.
Identification of genetic variants in PIKFYVE
Request a detailed protocolWES and Sanger sequencing were used in this study to identify the p.G1943E genetic variant in the congenital cataract family. In brief, DNA was extracted from blood or hair samples from cataract patients. Exome capture was performed using a SureSelect Human All Exon V6 Capture Kit (Agilent Technologies, Santa Clara, CA). WES was performed on four selected family members (II-1, II-5, III-9, and IV-5) of Chinese Korean family. Samples were sequenced on a HiSeq 4000 next-generation sequencing system (Illumina, San Diego, CA). Sequence reads were aligned to the human reference genome (hg19) using the Burrows-Wheeler Aligner (BWA) (Li and Durbin, 2009). Sequence variants were called using the Genome Analysis Toolkit (GATK) with the Best Practices for SNP and Indel discovery in germline DNA (McKenna et al., 2010). The variants were annotated using Annotate Variation (ANNOVAR) (Wang et al., 2010) and then filtered based on the following criteria: (1) nonsynonymous SNPs or indels in the exon region or the splice site region; (2) novel or minor allele frequency <1% in gnomAD; (3) PolyPhen score >0.2; (4) genomic evolutionary rate profiling (GERP) score >2.5; and (5) combined annotation dependent depletion score (CADD) >10. Polymerase chain reaction (PCR) and Sanger sequencing were used to verify the identified mutations from WES and screen the other family members for the mutations. The primer sequences used in PCR and Sanger sequencing were: (1) forward primer: 5′-TTTTGACCTTCTCTTGATTAGAGG-3′ and (2) reverse primer: 5′-AAATATGGCCTAGTAACCAAAGTTAAA-3′.
For another congenital cataract family and 10 sporadic cases with congenital cataract, gene panel testing of ophthalmology (testing of known pathological genes in cataract and PIKFYVE) were used, while for 200 patients with age-related cataract or complicated cataract, fast-target sequencing of PIKFYVE exons were utilized to screen for other potential pathological variants in PIKFYVE. The variants were filtered based on the following criteria: (1) sorting intolerant from tolerant (SIFT) score ≤ 0.05; (2) deleterious annotation of genetic variants using neural networks (DANN) score >0.93; (3) variant effect scoring tool (VEST) score from 0 to 1, the higher score, the more likely mutation might cause a functional change; (4) ClinPred score >0.5; and (5) MutationTaster score >0.5.
Cell culture
Request a detailed protocolThe HLEB3 cells (ATCC CRL-11421) were ordered from ATCC. HEK293T cells were obtained from the Laboratory cell bank of Shenzhen PKU-HKUST Medical Center, which were initially ordered from ATCC by colleagues in HKUST. Cell culture was performed according to the instructions (Andley et al., 1994). The mycoplasma contamination was prevented by examining the integrity of cell nucleus via DAPI staining.
YM201636 treatment
Request a detailed protocolThe HLEB3 cells were plated into six-well plates and grown for 4 days. The cells were then treated with YM201636 (Selleck, China) at a working concentration of 800 nM for 4 hr before imaging. Cells only treated with DMSO were used as control.
Quantitative RT-PCR
Request a detailed protocolAnterior lens capsules from three individuals with cataract were obtained during the surgery. Total RNA was extracted with FastPure Cell/Tissue Total RNA Isolation Kit (RC112-01, Vazyme, China) and subjected to reverse transcription with SuperScript III First-Strand Synthesis System (18080051, Thermo Fisher Scientific, USA). Quantitative RT-PCR was conducted with iTaq Universal SYBR Green Supermix (1725124, Bio-Rad, USA) in CFX96 Real-Time PCR Detection System.
Protein detection by western blot
Request a detailed protocolTransfection of the HEK293T cells was performed with Lipofectamine RNAiMAX Reagent (13778030, Thermo Fisher Scientific, USA) according to the manufacturer’s instructions. Cells were transfected with 2 µg pCS2(+)-CMV-PIKFYVEWT or pCS2(+)-CMV-PIKFYVEG1943E and harvested after 34 and 44 hr for Western blot analysis. Expression of the housekeeping gene GAPDH was used as the loading control. The ImageJ software (National Institutes of Health, Bethesda, MD, https://imagej.nih.gov/ij/) was used to measure the intensity of the bands and quantify the expression of PIKFYVE protein. The anti-PIKFYVE (ab137907, Abcam, Cambridge, UK), anti-GAPDH (ab9485, Abcam) primary antibodies, and goat anti-rabbit (ab205718, Abcam) secondary antibodies were used in this study.
Zebrafish husbandry
Request a detailed protocolZebrafish were raised and maintained according to the standard protocols (Westerfield, 1993). The WT AB strain, pikfyveΔ8mutants, and transgenic lines including Tg(cryaa:DsRed;il-1b:GFP-F) (Nguyen-Chi et al., 2014), Tg(ubi:PIKFYVEWT), and Tg(ubi:PIKFYVEG1943E) were used in this study.
Generation of zebrafish Pikfyve mutant
Request a detailed protocolThe sgRNA was designed with the online tool Crisprscan (Vejnar et al., 2016). The sgRNA was synthesized by in vitro transcription with the MEGAshortscript T7 Transcription Kit (AM1354, Ambion, Austin, TX) and purified by the MEGAclear Kit (AM1908, Ambion). The sgRNA, together with the Cas9 protein (EnGen Cas9 NLS, S. pyogenes, M0646M, NEB, Ipswich, MA), was injected into zebrafish embryos at one-cell stage. The final concentration was 100 ng/μl for each sgRNA and 800 ng/μl for Cas9 protein. Stable mutant line was screened by the T7 endonuclease 1 (T7E1) digestion and Sanger sequencing.
Generation of transgenic zebrafish lines
Request a detailed protocolThe coding sequence of WT and G1943E variant form of human PIKFYVE was placed downstream of the ubi promoter (Mosimann et al., 2011), which were cloned to a modified PBSK vector containing two arms of Tol2 elements (Kawakami et al., 2000). The plasmid and transposon mRNA were then injected into embryos at one-cell stage. Injected embryos were raised into adulthood and screened for germline transmission by PCR and WISH (Zhang and Liu, 2013) of F1 embryos.
Histology
Request a detailed protocolZebrafish embryos were fixed in 4% paraformaldehyde at 5 dpf and then dehydrated in 30% sucrose. The whole embryos were mounted in optimal cutting temperature compound and frozen in –80°C. Cryo-section was performed at 10 µM intervals in transverse planes from the head. HE staining and antibody staining were conducted according to the standard protocols. ZL-1 antibody (ab185979, Abcam) was used in this study. TEM analysis of the lens was completed by Servicebio Technology (Wuhan, China).
Imaging of zebrafish lens
Request a detailed protocolImages were taken under ZEISS LSM980 Confocal Laser Scanning Microscope (Carl Zeiss Microscopy, NY). The 63× objective (Plan-Apochromat 63×/1.40 Oil DIC M27) was used in this study. For time-lapse imaging experiment, 4-dpf zebrafish was mounted in 1.0% low-melting agarose with 0.02% tricaine. The lens was directly imaged under ZEISS Celldiscoverer 7 microscope (Carl Zeiss Microscopy). A 50× water lens (Plan-Apochromat 50×/1.2) was used and images were taken every 3 min. Images were analyzed by ImageJ (National Institutes of Health).
Detection of endosomes and autophagosomes in zebrafish lens
Request a detailed protocolCoding sequences of rab5c, rab7, or rab11a were fused with GFP sequence and cloned into the pCS2+ vector. mRNAs were synthesized in vitro using the mMESSAGE mMACHINE SP6 Kit (AM1908, Ambion). About 1–2 nl mRNA (100 ng/μl) was injected into zebrafish embryos at one-cell stage. At 3.5 dpf, embryos were mounted in 1% low-melting agarose and imaged under the confocal microscope.
Baf-A1 treatment
Request a detailed protocolBaf-A1 powder (HY-100558, MedChemExpress, NJ) was dissolved to the concentration of 1 mM in DMSO for stock, and diluted to 1 μM in egg water before use. Zebrafish embryos were treated with either 1‰ DMSO or 1 μM Baf-A1 solution at 28.5°C for 4.5 hr and washed out by egg water before confocal imaging.
Statistical analysis
Request a detailed protocolTwo-tailed Student’s t-tests were used for comparisons between two groups. Data were represented as mean ± SD. Statistical significance was shown as n.s., p>0.05, *p<0.05, ***p<0.001, and ****p<0.0001.
Data availability
All data generated or analysed during this study are included in the manuscript and supporting files. Source data files have been provided for Table 2.
-
Genome Aggregation DatabaseID r2_1. The genetic landscape of crystallins in congenital cataract.
References
-
Propagation and immortalization of human lens epithelial cells in cultureInvestigative Ophthalmology & Visual Science 35:3094–3102.
-
Inherited cataracts: molecular genetics, clinical features, disease mechanisms and novel therapeutic approachesThe British Journal of Ophthalmology 104:1331–1337.https://doi.org/10.1136/bjophthalmol-2019-315282
-
The genetic landscape of crystallins in congenital cataractOrphanet Journal of Rare Diseases 15:333.https://doi.org/10.1186/s13023-020-01613-3
-
Connexin mutants and cataractsFrontiers in Pharmacology 4:43.https://doi.org/10.3389/fphar.2013.00043
-
Crystallins, genes and cataractProgress in Drug Research. Fortschritte Der Arzneimittelforschung 60:205–262.https://doi.org/10.1007/978-3-0348-8012-1_7
-
Active vacuolar H+ ATPase and functional cycle of Rab5 are required for the vacuolation defect triggered by PtdIns(3,5)P2 loss under PIKfyve or Vps34 deficiencyAmerican Journal of Physiology. Cell Physiology 311:C366–C377.https://doi.org/10.1152/ajpcell.00104.2016
-
EPHA2 MUTATIONS CONTRIBUTE TO CONGENITAL CATARACT THROUGH DIVERSE MECHANISMSMolecular Vision 22:18–30.
-
Identification of novel PIKFYVE gene mutations associated with Fleck corneal dystrophyMolecular Vision 21:1093–1100.
-
Membrane potential mediates H(+)-ATPase dependence of “degradative pathway” endosomal fusionThe Journal of Membrane Biology 162:157–167.https://doi.org/10.1007/s002329900353
-
PI5P and PI(3,5)P2: Minor, but Essential PhosphoinositidesCell Structure and Function 42:49–60.https://doi.org/10.1247/csf.17003
-
Localized PtdIns 3,5-P2 synthesis to regulate early endosome dynamics and fusionAmerican Journal of Physiology-Cell Physiology 291:C393–C404.https://doi.org/10.1152/ajpcell.00019.2006
-
A novel mutation (p.Glu1389AspfsX16) of the phosphoinositide kinase, FYVE finger containing gene found in a Japanese patient with fleck corneal dystrophyMolecular Vision 18:2954–2960.
-
Pediatric cataractIndian Journal of Ophthalmology 65:1340–1349.https://doi.org/10.4103/ijo.IJO_1023_17
-
A novel PIKFYVE mutation in fleck corneal dystrophyMolecular Vision 17:2776–2781.
-
Rab GTPase Function in Endosome and Lysosome BiogenesisTrends in Cell Biology 28:957–970.https://doi.org/10.1016/j.tcb.2018.06.007
-
Mutations in PIP5K3 are associated with François-Neetens mouchetée fleck corneal dystrophyAmerican Journal of Human Genetics 77:54–63.https://doi.org/10.1086/431346
-
Molecular genetics of congenital cataractsExperimental Eye Research 191:107872.https://doi.org/10.1016/j.exer.2019.107872
-
Transient infection of the zebrafish notochord with E. coli induces chronic inflammationDisease Models & Mechanisms 7:871–882.https://doi.org/10.1242/dmm.014498
-
A novel mutation of LIM2 causes autosomal dominant membranous cataract in a Chinese familyInternational Journal of Ophthalmology 13:1512–1520.https://doi.org/10.18240/ijo.2020.10.02
-
Sequence and spatial expression of zebrafish (Danio rerio) alphaA-crystallinMolecular Vision 8:45–50.
-
The mammalian phosphatidylinositol 3-phosphate 5-kinase (PIKfyve) regulates endosome-to-TGN retrograde transportJournal of Cell Science 119:3944–3957.https://doi.org/10.1242/jcs.03153
-
CRISPR/Cas9‐Directed Gene Editing for the Generation of Loss‐of‐Function Mutants in High‐Throughput Zebrafish F0 ScreensCurrent Protocols in Molecular Biology 119:31.https://doi.org/10.1002/cpmb.42
-
CHMP4B, a novel gene for autosomal dominant cataracts linked to chromosome 20qAmerican Journal of Human Genetics 81:596–606.https://doi.org/10.1086/519980
-
Mutations and mechanisms in congenital and age-related cataractsExperimental Eye Research 156:95–102.https://doi.org/10.1016/j.exer.2016.06.011
-
PIKfyve: Partners, significance, debates and paradoxesCell Biology International 32:591–604.https://doi.org/10.1016/j.cellbi.2008.01.006
-
Functions of the intermediate filament cytoskeleton in the eye lensThe Journal of Clinical Investigation 119:1837–1848.https://doi.org/10.1172/JCI38277
-
Optimized CRISPR–Cas9 System for Genome Editing in ZebrafishCold Spring Harbor Protocols 2016:prot086850.https://doi.org/10.1101/pdb.prot086850
-
PI(5)P regulates autophagosome biogenesisMolecular Cell 57:219–234.https://doi.org/10.1016/j.molcel.2014.12.007
-
ANNOVAR: functional annotation of genetic variants from high-throughput sequencing dataNucleic Acids Research 38:e164.https://doi.org/10.1093/nar/gkq603
-
BookThe Zebrafish Book: A Guide for the Laboratory Use of Zebrafish (Brachydanio RerioEugene: M. Westerfield.
-
Bafilomycin A1, a specific inhibitor of vacuolar-type H(+)-ATPase, inhibits acidification and protein degradation in lysosomes of cultured cellsThe Journal of Biological Chemistry 266:17707–17712.https://doi.org/10.1016/S0021-9258(19)47429-2
-
Mutation screening of crystallin genes in Chinese families with congenital cataractsMolecular Vision 25:427–437.
Article and author information
Author details
Funding
Science, Technology and Innovation Commission of Shenzhen Municipality (GJHZ20180420180937076)
- Jun Zhao
Science, Technology and Innovation Commission of Shenzhen Municipality (JCYJ20180228164400218)
- Jun Zhao
Sanming Project of Medicine in Shenzhen (SZSM201812090)
- Jun Zhao
The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.
Acknowledgements
The authors would like to thank all family members who participated in this study. The authors thank Dr. Georges Lutfalla from Institut Pasteur for providing Tg(cryaa:DsRed;il-1b:GFP-F) zebrafish. The authors also thank Dr. Keyu Chen from Shenzhen Bay Laboratory for the structure analysis on PIKFYVE. This study was supported by Science, Technology and Innovation Commission of Shenzhen Municipality Grants (GJHZ 20180420180937076 and JCYJ20180228164400218) and Sanming Project of Medicine in Shenzhen Grant (SZSM201812090).
Ethics
Ethical approval was provided by the Medical Sciences Ethics Committee at Shenzhen Eye Hospital Affiliated to Jinan University, China (No.: 20161103005). The ethical protocol in this study was approved by the Independent Ethics Committee of Shenzhen Eye Hospital, in accordance with the tenets of the Declaration of Helsinki (1975). Written informed consent and consent to publish were received from all study participants prior to obtaining the human tissue samples. Samples were de-identified by removing the identifiable individual information such as name and identity number before researching.
Copyright
© 2022, Mei et al.
This article is distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use and redistribution provided that the original author and source are credited.
Metrics
-
- 1,019
- views
-
- 173
- downloads
-
- 15
- citations
Views, downloads and citations are aggregated across all versions of this paper published by eLife.
Download links
Downloads (link to download the article as PDF)
Open citations (links to open the citations from this article in various online reference manager services)
Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)
Further reading
-
- Genetics and Genomics
Recent studies have shown that, in human cancer cells, the tetrameric Shieldin complex (comprising REV7, SHLD1, SHLD2, and SHLD3) facilitates non-homologous end-joining (NHEJ) while blocking homologous recombination (HR). Surprisingly, several eukaryotic species lack SHLD1, SHLD2, and SHLD3 orthologs, suggesting that Rev7 may leverage an alternative mechanism to regulate the double-strand break (DSB) repair pathway choice. Exploring this hypothesis, we discovered that Saccharomyces cerevisiae Rev7 physically interacts with the Mre11–Rad50–Xrs2 (MRX) subunits, impedes G-quadruplex DNA synergized HU-induced toxicity, and facilitates NHEJ, while antagonizing HR. Notably, we reveal that a 42-amino acid C-terminal fragment of Rev7 binds to the subunits of MRX complex, protects rev7∆ cells from G-quadruplex DNA-HU-induced toxicity, and promotes NHEJ by blocking HR. By comparison, the N-terminal HORMA domain, a conserved protein–protein interaction module, was dispensable. We further show that the full-length Rev7 impedes Mre11 nuclease and Rad50’s ATPase activities without affecting the latter’s ATP-binding ability. Combined, these results provide unanticipated insights into the functional interaction between the MRX subunits and Rev7 and highlight a previously unrecognized mechanism by which Rev7 facilitates DSB repair via NHEJ, and attenuation of HR, by blocking Mre11 nuclease and Rad50’s ATPase activities in S. cerevisiae.
-
- Genetics and Genomics
Over the last three decades, human genetics has gone from dissecting high-penetrance Mendelian diseases to discovering the vast and complex genetic etiology of common human diseases. In tackling this complexity, scientists have discovered the importance of numerous genetic processes – most notably functional regulatory elements – in the development and progression of these diseases. Simultaneously, scientists have increasingly used multiplex assays of variant effect to systematically phenotype the cellular consequences of millions of genetic variants. In this article, we argue that the context of genetic variants – at all scales, from other genetic variants and gene regulation to cell biology to organismal environment – are critical components of how we can employ genomics to interpret these variants, and ultimately treat these diseases. We describe approaches to extend existing experimental assays and computational approaches to examine and quantify the importance of this context, including through causal analytic approaches. Having a unified understanding of the molecular, physiological, and environmental processes governing the interpretation of genetic variants is sorely needed for the field, and this perspective argues for feasible approaches by which the combined interpretation of cellular, animal, and epidemiological data can yield that knowledge.