MS/MS in silico subtraction-based proteomic profiling as an approach to facilitate disease gene discovery: application to lens development and cataract

Variant in EZR leads to defects in lens development

BACKGROUND

Congenital cataract is a common cause of blindness. Genetic factors always play important role.

MATERIAL AND METHODS

This study identified a novel missense variant (c.1412C>T (p.P471L)) in the EZR gene in a four-generation Chinese family with nuclear cataract by linkage analysis and whole-exome sequencing. A knockout study in zebrafish using transcription activator-like effector nucleases was carried out to gain insight into candidate gene function.

RESULTS

Conservative and functional prediction suggests that the P-to-L substitution may impair the function of the human ezrin protein. Histology showed developmental delays in the ezrin-mutated zebrafish, manifesting as multilayered lens epithelial cells. Immunohistochemistry revealed abnormal proliferation patterns in mutant fish.

CONCLUSIONS

The study suggests that ezrin may be involved in the enucleation and differentiation of lens epithelial cells.

Proteomic profiling of retina and retinal pigment epithelium combined embryonic tissue to facilitate ocular disease gene discovery

To expedite gene discovery in eye development and its associated defects, we previously developed a bioinformatics resource-tool iSyTE (integrated Systems Tool for Eye gene discovery). However, iSyTE is presently limited to lens tissue and is predominantly based on transcriptomics datasets. Therefore, to extend iSyTE to other eye tissues on the proteome level, we performed high-throughput tandem mass spectrometry (MS/MS) on mouse embryonic day (E)14.5 retina and retinal pigment epithelium combined tissue and identified an average of 3300 proteins per sample (n = 5). High-throughput expression profiling-based gene discovery approaches-involving either transcriptomics or proteomics-pose a key challenge of prioritizing candidates from thousands of RNA/proteins expressed. To address this, we used MS/MS proteome data from mouse whole embryonic body (WB) as a reference dataset and performed comparative analysis-termed "in silico WB-subtraction"-with the retina proteome dataset. In silico WB-subtraction identified 90 high-priority proteins with retina-enriched expression at stringency criteria of ≥ 2.5 average spectral counts, ≥ 2.0 fold-enrichment, false discovery rate < 0.01. These top candidates represent a pool of retina-enriched proteins, several of which are associated with retinal biology and/or defects (e.g., Aldh1a1, Ank2, Ank3, Dcn, Dync2h1, Egfr, Ephb2, Fbln5, Fbn2, Hras, Igf2bp1, Msi1, Rbp1, Rlbp1, Tenm3, Yap1, etc.), indicating the effectiveness of this approach. Importantly, in silico WB-subtraction also identified several new high-priority candidates with potential regulatory function in retina development. Finally, proteins exhibiting expression or enriched-expression in the retina are made accessible in a user-friendly manner at iSyTE ( https://research.bioinformatics.udel.edu/iSyTE/ ), to allow effective visualization of this information and facilitate eye gene discovery.

High-Throughput Transcriptomics of Celf1 Conditional Knockout Lens Identifies Downstream Networks Linked to Cataract Pathology

Defects in the development of the ocular lens can cause congenital cataracts. To understand the various etiologies of congenital cataracts, it is important to characterize the genes linked to this developmental defect and to define their downstream pathways that are relevant to lens biology and pathology. Deficiency or alteration of several RNA-binding proteins, including the conserved RBP Celf1 (CUGBP Elav-like family member 1), has been described to cause lens defects and early onset cataracts in animal models and/or humans. Celf1 is involved in various aspects of post-transcriptional gene expression control, including regulation of mRNA stability/decay, alternative splicing and translation. Celf1 germline knockout mice and lens conditional knockout (Celf1cKO) mice develop fully penetrant cataracts in early postnatal stages. To define the genome-level changes in RNA transcripts that result from Celf1 deficiency, we performed high-throughput RNA-sequencing of Celf1cKO mouse lenses at postnatal day (P) 0. Celf1cKO lenses exhibit 987 differentially expressed genes (DEGs) at cut-offs of >1.0 log2 counts per million (CPM), ≥±0.58 log2 fold-change and <0.05 false discovery rate (FDR). Of these, 327 RNAs were reduced while 660 were elevated in Celf1cKO lenses. The DEGs were subjected to various downstream analyses including iSyTE lens enriched-expression, presence in Cat-map, and gene ontology (GO) and representation of regulatory pathways. Further, a comparative analysis was done with previously generated microarray datasets on Celf1cKO lenses P0 and P6. Together, these analyses validated and prioritized several key genes mis-expressed in Celf1cKO lenses that are relevant to lens biology, including known cataract-linked genes (e.g., Cryab, Cryba2, Cryba4, Crybb1, Crybb2, Cryga, Crygb, Crygc, Crygd, Cryge, Crygf, Dnase2b, Bfsp1, Gja3, Pxdn, Sparc, Tdrd7, etc.) as well as novel candidates (e.g., Ell2 and Prdm16). Together, these data have defined the alterations in lens transcriptome caused by Celf1 deficiency, in turn uncovering downstream genes and pathways (e.g., structural constituents of eye lenses, lens fiber cell differentiation, etc.) associated with lens development and early-onset cataracts.

Proteomic profiling of retina and retinal pigment epithelium combined embryonic tissue to facilitate ocular disease gene discovery

To expedite gene discovery in eye development and its associated defects, we previously developed a bioinformatics resource-tool iSyTE (integrated Systems Tool for Eye gene discovery). However, iSyTE is presently limited to lens tissue and is predominantly based on transcriptomics datasets. Therefore, to extend iSyTE to other eye tissues on the proteome level, we performed high-throughput tandem mass spectrometry (MS/MS) on mouse embryonic day (E)14.5 retina and retinal pigment epithelium combined tissue and identified an average of 3,300 proteins per sample (n=5). High-throughput expression profiling-based gene discovery approaches-involving either transcriptomics or proteomics-pose a key challenge of prioritizing candidates from thousands of RNA/proteins expressed. To address this, we used MS/MS proteome data from mouse whole embryonic body (WB) as a reference dataset and performed comparative analysis-termed "in silico WB-subtraction"-with the retina proteome dataset. In silico WB-subtraction identified 90 high-priority proteins with retina-enriched expression at stringency criteria of ³2.5 average spectral counts, ³2.0 fold-enrichment, False Discovery Rate <0.01. These top candidates represent a pool of retina-enriched proteins, several of which are associated with retinal biology and/or defects (e.g., Aldh1a1, Ank2, Ank3, Dcn, Dync2h1, Egfr, Ephb2, Fbln5, Fbn2, Hras, Igf2bp1, Msi1, Rbp1, Rlbp1, Tenm3, Yap1, etc.), indicating the effectiveness of this approach. Importantly, in silico WB-subtraction also identified several new high-priority candidates with potential regulatory function in retina development. Finally, proteins exhibiting expression or enriched-expression in the retina are made accessible in a user-friendly manner at iSyTE (https://research.bioinformatics.udel.edu/iSyTE/), to allow effective visualization of this information and facilitate eye gene discovery.

Biological properties of the BCL-2 family protein BCL-RAMBO, which regulates apoptosis, mitochondrial fragmentation, and mitophagy

Mitochondria play an essential role in the regulation of cellular stress responses, including cell death. Damaged mitochondria are removed by fission and fusion cycles and mitophagy, which counteract cell death. BCL-2 family proteins possess one to four BCL-2 homology domains and regulate apoptosis signaling at mitochondria. BCL-RAMBO, also known as BCL2-like 13 (BCL2L13), was initially identified as one of the BCL-2 family proteins inducing apoptosis. Mitophagy receptors recruit the ATG8 family proteins MAP1LC3/GABARAP via the MAP1LC3-interacting region (LIR) motif to initiate mitophagy. In addition to apoptosis, BCL-RAMBO has recently been identified as a mitophagy receptor that possesses the LIR motif and regulates mitochondrial fragmentation and mitophagy. In the 20 years since its discovery, many important findings on BCL-RAMBO have been increasingly reported. The biological properties of BCL-RAMBO are reviewed herein.

A possible connection between reactive oxygen species and the unfolded protein response in lens development: From insight to foresight

The lens is a relatively special and simple organ. It has become an ideal model to study the common developmental characteristics among different organic systems. Lens development is a complex process influenced by numerous factors, including signals from the intracellular and extracellular environment. Reactive oxygen species (ROS) are a group of highly reactive and oxygen-containing molecules that can cause endoplasmic reticulum stress in lens cells. As an adaptive response to ER stress, lens cells initiate the unfolded protein response (UPR) to maintain normal protein synthesis by selectively increasing/decreasing protein synthesis and increasing the degradation of misfolded proteins. Generally, the UPR signaling pathways have been well characterized in the context of many pathological conditions. However, recent studies have also confirmed that all three UPR signaling pathways participate in a variety of developmental processes, including those of the lens. In this review, we first briefly summarize the three stages of lens development and present the basic profiles of ROS and the UPR. We then discuss the interconnections between lens development and these two mechanisms. Additionally, the potential adoption of human pluripotent stem-cell-based lentoids in lens development research is proposed to provide a novel perspective on future developmental studies.

Deficiency of the bZIP transcription factors Mafg and Mafk causes misexpression of genes in distinct pathways and results in lens embryonic developmental defects

Deficiency of the small Maf proteins Mafg and Mafk cause multiple defects, namely, progressive neuronal degeneration, cataract, thrombocytopenia and mid-gestational/perinatal lethality. Previous data shows Mafg ^-/-:Mafk ^+/- compound knockout (KO) mice exhibit cataracts age 4-months onward. Strikingly, Mafg ^-/-:Mafk ^-/- double KO mice develop lens defects significantly early in life, during embryogenesis, but the pathobiology of these defects is unknown, and is addressed here. At embryonic day (E)16.5, the epithelium of lens in Mafg ^-/-:Mafk ^-/- animals appears abnormally multilayered as demonstrated by E-cadherin and nuclear staining. Additionally, Mafg ^-/-:Mafk ^-/- lenses exhibit abnormal distribution of F-actin near the "fulcrum" region where epithelial cells undergo apical constriction prior to elongation and reorientation as early differentiating fiber cells. To identify the underlying molecular changes, we performed high-throughput RNA-sequencing of E16.5 Mafg ^-/-:Mafk ^-/- lenses and identified a cohort of differentially expressed genes that were further prioritized using stringent filtering criteria and validated by RT-qPCR. Several key factors associated with the cytoskeleton, cell cycle or extracellular matrix (e.g., Cdk1, Cdkn1c, Camsap1, Col3a1, Map3k12, Sipa1l1) were mis-expressed in Mafg ^-/-:Mafk ^-/- lenses. Further, the congenital cataract-linked extracellular matrix peroxidase Pxdn was significantly overexpressed in Mafg ^-/-:Mafk ^-/- lenses, which may cause abnormal cell morphology. These data also identified the ephrin signaling receptor Epha5 to be reduced in Mafg ^-/-:Mafk ^-/- lenses. This likely contributes to the Mafg ^-/-:Mafk ^-/- multilayered lens epithelium pathology, as loss of an ephrin ligand, Efna5 (ephrin-A5), causes similar lens defects. Together, these findings uncover a novel early function of Mafg and Mafk in lens development and identify their new downstream regulatory relationships with key cellular factors.

RNA-binding proteins and post-transcriptional regulation in lens biology and cataract: Mediating spatiotemporal expression of key factors that control the cell cycle, transcription, cytoskeleton and transparency

Development of the ocular lens - a transparent tissue capable of sustaining frequent shape changes for optimal focusing power - pushes the boundaries of what cells can achieve using the molecular toolkit encoded by their genomes. The mammalian lens contains broadly two types of cells, the anteriorly located monolayer of epithelial cells which, at the equatorial region of the lens, initiate differentiation into fiber cells that contribute to the bulk of the tissue. This differentiation program involves massive upregulation of select fiber cell-expressed RNAs and their subsequent translation into high amounts of proteins, such as crystallins. But intriguingly, fiber cells achieve this while also simultaneously undergoing significant morphological changes such as elongation - involving about 1000-fold length-wise increase - and migration, which requires modulation of cytoskeletal and cell adhesion factors. Adding further to the challenges, these molecular and cellular events have to be coordinated as fiber cells progress toward loss of their nuclei and organelles, which irreversibly compromises their potential for harnessing genetically hardwired information. A long-standing question is how processes downstream of signaling and transcription, which may also participate in feedback regulation, contribute toward orchestrating these cellular differentiation events in the lens. It is now becoming clear from findings over the past decade that post-transcriptional gene expression regulatory mechanisms are critical in controlling cellular proteomes and coordinating key processes in lens development and fiber cell differentiation. Indeed, RNA-binding proteins (RBPs) such as Caprin2, Celf1, Rbm24 and Tdrd7 have now been described in mediating post-transcriptional control over key factors (e.g. Actn2, Cdkn1a (p21Cip1), Cdkn1b (p27Kip1), various crystallins, Dnase2b, Hspb1, Pax6, Prox1, Sox2) that are variously involved in cell cycle, transcription, cytoskeleton maintenance and differentiation in the lens. Furthermore, deficiencies of these RBPs have been shown to result in various eye and lens defects and/or cataract. Because fiber cell differentiation in the lens occurs throughout life, the underlying regulatory mechanisms operational in development are expected to also be recruited for the maintenance of transparency in aged lenses. Indeed, in support of this, TDRD7 and CAPRIN2 loci have been linked to age-related cataract in humans. Here, I will review the role of key RBPs in the lens and their importance in understanding the pathology of lens defects. I will discuss advances in RBP-based gene expression control, in general, and the important challenges that need to be addressed in the lens to define the mechanisms that determine the epithelial and fiber cell proteome. Finally, I will also discuss in detail several key future directions including the application of bioinformatics approaches such as iSyTE to study RBP-based post-transcriptional gene expression control in the aging lens and in the context of age-related cataract.

A Novel Mutation in Cse1l Disrupts Brain and Eye Development with Specific Effects on Pax6 Expression

Forward genetics in the mouse continues to be a useful and unbiased approach to identifying new genes and alleles with previously unappreciated roles in mammalian development and disease. Here, we report a new mouse allele of Cse1l that was recovered from an ENU mutagenesis screen. Embryos homozygous for the anteater allele of Cse1l display a number of variable phenotypes, with craniofacial and ocular malformations being the most obvious. We provide evidence that Cse1l is the causal gene through complementation with a novel null allele of Cse1l generated by CRISPR-Cas9 editing. While the variability in the anteater phenotype was high enough to preclude a detailed molecular analysis, we demonstrate a very penetrant reduction in Pax6 levels in the developing eye along with significant ocular developmental phenotypes. The eye gene discovery tool iSyTE shows Cse1l to be significantly expressed in the lens from early eye development stages in embryos through adulthood. Cse1l has not previously been shown to be required for organogenesis as homozygosity for a null allele results in very early lethality. Future detailed studies of Cse1l function in craniofacial and neural development will be best served with a conditional allele to circumvent the variable phenotypes we report here. We suggest that human next-generation (whole genome or exome) sequencing studies yielding variants of unknown significance in CSE1L could consider these findings as part of variant analysis.

Crystallin gene expression: Insights from studies of transcriptional bursting

Cellular differentiation is marked by temporally and spatially regulated gene expression. The ocular lens is one of the most powerful mammalian model system since it is composed from only two cell subtypes, called lens epithelial and fiber cells. Lens epithelial cells differentiate into fiber cells through a series of spatially and temporally orchestrated processes, including massive production of crystallins, cellular elongation and the coordinated degradation of nuclei and other organelles. Studies of transcriptional and posttranscriptional gene regulatory mechanisms in lens provide a wide range of opportunities to understand global molecular mechanisms of gene expression as steady-state levels of crystallin mRNAs reach very high levels comparable to globin genes in erythrocytes. Importantly, dysregulation of crystallin gene expression results in lens structural abnormalities and cataracts. The mRNA life cycle is comprised of multiple stages, including transcription, splicing, nuclear export into cytoplasm, stabilization, localization, translation and ultimate decay. In recent years, development of modern mRNA detection methods with single molecule and single cell resolution enabled transformative studies to visualize the mRNA life cycle to generate novel insights into the sequential regulatory mechanisms of gene expression during embryogenesis. This review is focused on recent major advancements in studies of transcriptional bursting in differentiating lens fiber cells, analysis of nascent mRNA expression from bi-directional promoters, transient nuclear accumulation of specific mRNAs, condensation of chromatin prior lens fiber cell denucleation, and outlines future studies to probe the interactions of individual mRNAs with specific RNA-binding proteins (RBPs) in the cytoplasm and regulation of translation and mRNA decay.

The cataract-linked RNA-binding protein Celf1 post-transcriptionally controls the spatiotemporal expression of the key homeodomain transcription factors Pax6 and Prox1 in lens development

The homeodomain transcription factors (TFs) Pax6 (OMIM: 607108) and Prox1 (OMIM: 601546) critically regulate gene expression in lens development. While PAX6 mutations in humans can cause cataract, aniridia, microphthalmia, and anophthalmia, among other defects, Prox1 deletion in mice causes severe lens abnormalities, in addition to other organ defects. Furthermore, the optimal dosage/spatiotemporal expression of these key TFs is essential for development. In lens development, Pax6 expression is elevated in cells of the anterior epithelium compared to fiber cells, while Prox1 exhibits the opposite pattern. Whether post-transcriptional regulatory mechanisms control these precise TF expression patterns is unknown. Here, we report the unprecedented finding that the cataract-linked RNA-binding protein (RBP), Celf1 (OMIM: 601074), post-transcriptionally regulates Pax6 and Prox1 protein expression in lens development. Immunostaining shows that Celf1 lens-specific conditional knockout (Celf1^cKO) mice exhibit abnormal elevation of Pax6 protein in fiber cells and abnormal Prox1 protein levels in epithelial cells-directly opposite to their normal expression patterns in development. Furthermore, RT-qPCR shows no change in Pax6 and Prox1 transcript levels in Celf1^cKO lenses, suggesting that Celf1 regulates these TFs on the translational level. Indeed, RNA-immunoprecipitation assays using Celf1 antibody indicate that Celf1 protein binds to Pax6 and Prox1 transcripts. Furthermore, reporter assays in Celf1 knockdown and Celf1-overexpression cells demonstrate that Celf1 negatively controls Pax6 and Prox1 translation via their 3' UTRs. These data define a new mechanism of RBP-based post-transcriptional regulation that enables precise control over spatiotemporal expression of Pax6 and Prox1 in lens development, thereby uncovering a new etiological mechanism for Celf1 deficiency-based cataract.

Generation and proteome profiling of PBMC-originated, iPSC-derived lentoid bodies

Here, we report proteome profiling of peripheral blood mononuclear cell (PBMC)-originated, induced pluripotent stem cell (iPSC)-derived, lens-like organoids termed lentoid bodies at two differentiation time points. A small aliquot of the blood sample was ascertained to collect PBMCs that were reprogrammed to iPSCs. The PBMC-originated, iPSCs were differentiated to lentoid bodies employing the "fried egg" method. Quantitative real-time PCR (qRT-PCR) analysis revealed increased expression levels of lens-associated markers in lentoid bodies while transmission electron microscopy identified closely packed lens epithelial- and differentiating fiber-like cells in lentoid bodies. Total cellular protein was extracted from lentoid bodies at differentiation day 25 and mass spectrometry identified a total of 9,473 proteins. The low counts of crystallin proteins at differentiation day 25 prompted us to re-examine the proteome at differentiation day 35 as we reasoned that 10 additional days of differentiation will increase the crystallin count. However, we did not detect any substantial increase in crystallin protein counts at differentiation day 35. In conclusion, we report generation and proteome profiles of PBMC-originated, iPSC-derived lentoid bodies at multiple differentiation time points.

Modeling ocular lens disease in Xenopus

BACKGROUND

Ocular lens clouding is termed as cataract, which depending on the onset, is classified as congenital or age-related. Developing new cataract treatments requires new models. Thus far, Xenopus embryos have not been evaluated as a system for studying cataract.

RESULTS

We characterized the developmental process of lens formation in Xenopus laevis tailbuds and tadpoles, and we disrupted the orthologues of three mammalian cataract-linked genes in F0 by CRISPR/Cas9. We assessed the consequences of gene inactivation by combining external examination with histochemical analyses and functional vision assays. Inactivating the key metazoan eye development transcription factor gene pax6 produces a strong eye phenotype including an absence of eye tissue. Inactivating the genes for gap-junction protein and a nuclease, gja8 and dnase2b, produces lens defects that share several features of human cataracts, including impaired vision acuity, nuclei retention in lens fiber cells, and actin fibers disorganization. We tested the potential improvement of the visual acuity of gja8 crispant tadpoles upon treatment with the molecular chaperone 4-phenylbutyrate.

CONCLUSION

Xenopus is a valuable model organism to understand the molecular pathology of congenital eye defects, including cataracts, and to screen molecules with a potential to prevent or reverse cataracts.