Compound mouse mutants of bZIP transcription factors Mafg and Mafk reveal a regulatory network of non-crystallin genes associated with cataract

Integrated single-cell multiomics uncovers foundational regulatory mechanisms of lens development and pathology

Ocular lens development entails epithelial to fiber cell differentiation, defects in which cause congenital cataracts. We report the first single-cell multiomic atlas of lens development, leveraging snRNA-seq, snATAC-seq and CUT&RUN-seq to discover previously unreported mechanisms of cell fate determination and cataract-linked regulatory networks. A comprehensive profile of cis- and trans-regulatory interactions, including for the cataract-linked transcription factor MAF, is established across a temporal trajectory of fiber cell differentiation. Furthermore, we identify an epigenetic paradigm of cellular differentiation, defined by progressive loss of the H3K27 methylation writer Polycomb repressive complex 2 (PRC2). PRC2 localizes to heterochromatin domains across master-regulator transcription factor gene bodies, suggesting it safeguards epithelial cell fate. Moreover, we demonstrate that FGF hyper-stimulation in vivo leads to MAF network activation and the emergence of novel lens cell states. Collectively, these data depict a comprehensive portrait of lens fiber cell differentiation, while defining regulatory effectors of cell identity and cataract formation.

Integrated single-cell multiomics uncovers foundational regulatory mechanisms of lens development and pathology

Ocular lens development entails epithelial to fiber cell differentiation, defects in which cause congenital cataract. We report the first single-cell multiomic atlas of lens development, leveraging snRNA-seq, snATAC-seq, and CUT&RUN-seq to discover novel mechanisms of cell fate determination and cataract-linked regulatory networks. A comprehensive profile of cis- and trans-regulatory interactions, including for the cataract-linked transcription factor MAF, is established across a temporal trajectory of fiber cell differentiation. Further, we divulge a conserved epigenetic paradigm of cellular differentiation, defined by progressive loss of H3K27 methylation writer Polycomb repressive complex 2 (PRC2). PRC2 localizes to heterochromatin domains across master-regulator transcription factor gene bodies, suggesting it safeguards epithelial cell fate. Moreover, we demonstrate that FGF hyper-stimulation in vivo leads to MAF network activation and the emergence of novel lens cell states. Collectively, these data depict a comprehensive portrait of lens fiber cell differentiation, while defining regulatory effectors of cell identity and cataract formation.

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.

Multiomics Analysis Reveals Novel Genetic Determinants for Lens Differentiation, Structure, and Transparency

Recent advances in next-generation sequencing and data analysis have provided new gateways for identification of novel genome-wide genetic determinants governing tissue development and disease. These advances have revolutionized our understanding of cellular differentiation, homeostasis, and specialized function in multiple tissues. Bioinformatic and functional analysis of these genetic determinants and the pathways they regulate have provided a novel basis for the design of functional experiments to answer a wide range of long-sought biological questions. A well-characterized model for the application of these emerging technologies is the development and differentiation of the ocular lens and how individual pathways regulate lens morphogenesis, gene expression, transparency, and refraction. Recent applications of next-generation sequencing analysis on well-characterized chicken and mouse lens differentiation models using a variety of omics techniques including RNA-seq, ATAC-seq, whole-genome bisulfite sequencing (WGBS), chip-seq, and CUT&RUN have revealed a wide range of essential biological pathways and chromatin features governing lens structure and function. Multiomics integration of these data has established new gene functions and cellular processes essential for lens formation, homeostasis, and transparency including the identification of novel transcription control pathways, autophagy remodeling pathways, and signal transduction pathways, among others. This review summarizes recent omics technologies applied to the lens, methods for integrating multiomics data, and how these recent technologies have advanced our understanding ocular biology and function. The approach and analysis are relevant to identifying the features and functional requirements of more complex tissues and disease states.

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.

The role and regulation of Maf proteins in cancer

The Maf proteins (Mafs) belong to basic leucine zipper transcription factors and are members of the activator protein-1 (AP-1) superfamily. There are two subgroups of Mafs: large Mafs and small Mafs, which are involved in a wide range of biological processes, such as the cell cycle, proliferation, oxidative stress, and inflammation. Therefore, dysregulation of Mafs can affect cell fate and is closely associated with diverse diseases. Accumulating evidence has established both large and small Mafs as mediators of tumor development. In this review, we first briefly describe the structure and physiological functions of Mafs. Then we summarize the upstream regulatory mechanisms that control the expression and activity of Mafs. Furthermore, we discuss recent studies on the critical role of Mafs in cancer progression, including cancer proliferation, apoptosis, metastasis, tumor/stroma interaction and angiogenesis. We also review the clinical implications of Mafs, namely their potential possibilities and limitations as biomarkers and therapeutic targets in cancer.

Dynamic changes in whole genome DNA methylation, chromatin and gene expression during mouse lens differentiation

BACKGROUND

Cellular differentiation is marked by temporally and spatially coordinated gene expression regulated at multiple levels. DNA methylation represents a universal mechanism to control chromatin organization and its accessibility. Cytosine methylation of CpG dinucleotides regulates binding of methylation-sensitive DNA-binding transcription factors within regulatory regions of transcription, including promoters and distal enhancers. Ocular lens differentiation represents an advantageous model system to examine these processes as lens comprises only two cell types, the proliferating lens epithelium and postmitotic lens fiber cells all originating from the epithelium.

RESULTS

Using whole genome bisulfite sequencing (WGBS) and microdissected lenses, we investigated dynamics of DNA methylation and chromatin changes during mouse lens fiber and epithelium differentiation between embryos (E14.5) and newborns (P0.5). Histone H3.3 variant chromatin landscapes were also generated for both P0.5 lens epithelium and fibers by chromatin immunoprecipitation followed by next generation sequencing (ChIP-seq). Tissue-specific features of DNA methylation patterns are demonstrated via comparative studies with embryonic stem (ES) cells and neural progenitor cells (NPCs) at Nanog, Pou5f1, Sox2, Pax6 and Six3 loci. Comparisons with ATAC-seq and RNA-seq data demonstrate that reduced methylation is associated with increased expression of fiber cell abundant genes, including crystallins, intermediate filament (Bfsp1 and Bfsp2) and gap junction proteins (Gja3 and Gja8), marked by high levels of histone H3.3 within their transcribed regions. Interestingly, Pax6-binding sites exhibited predominantly DNA hypomethylation in lens chromatin. In vitro binding of Pax6 proteins showed Pax6's ability to interact with sites containing one or two methylated CpG dinucleotides.

CONCLUSIONS

Our study has generated the first data on methylation changes between two different stages of mammalian lens development and linked these data with chromatin accessibility maps, presence of histone H3.3 and gene expression. Reduced DNA methylation correlates with expression of important genes involved in lens morphogenesis and lens fiber cell differentiation.

Mapping the Universe of Eph Receptor and Ephrin Ligand Transcripts in Epithelial and Fiber Cells of the Eye Lens

The eye lens is a transparent, ellipsoid organ in the anterior chamber of the eye that is required for fine focusing of light onto the retina to transmit a clear image. Cataracts, defined as any opacity in the lens, remains the leading cause of blindness in the world. Recent studies in humans and mice indicate that Eph–ephrin bidirectional signaling is important for maintaining lens transparency. Specifically, mutations and polymorphisms in the EphA2 receptor and the ephrin-A5 ligand have been linked to congenital and age-related cataracts. It is unclear what other variants of Ephs and ephrins are expressed in the lens or whether there is preferential expression in epithelial vs. fiber cells. We performed a detailed analysis of Eph receptor and ephrin ligand mRNA transcripts in whole mouse lenses, epithelial cell fractions, and fiber cell fractions using a new RNA isolation method. We compared control samples with EphA2 knockout (KO) and ephrin-A5 KO samples. Our results revealed the presence of transcripts for 12 out of 14 Eph receptors and 8 out of 8 ephrin ligands in various fractions of lens cells. Using specific primer sets, RT-PCR, and sequencing, we verified the variant of each gene that is expressed, and we found two epithelial-cell-specific genes. Surprisingly, we also identified one Eph receptor variant that is expressed in KO lens fibers but is absent from control lens fibers. We also identified one low expression ephrin variant that is only expressed in ephrin-A5 control samples. These results indicate that the lens expresses almost all Ephs and ephrins, and there may be many receptor–ligand pairs that play a role in lens homeostasis.

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.

Changes in DNA methylation hallmark alterations in chromatin accessibility and gene expression for eye lens differentiation

Background

Methylation at cytosines (mCG) is a well-known regulator of gene expression, but its requirements for cellular differentiation have yet to be fully elucidated. A well-studied cellular differentiation model system is the eye lens, consisting of a single anterior layer of epithelial cells that migrate laterally and differentiate into a core of fiber cells. Here, we explore the genome-wide relationships between mCG methylation, chromatin accessibility and gene expression during differentiation of eye lens epithelial cells into fiber cells.

Results

Whole genome bisulfite sequencing identified 7621 genomic loci exhibiting significant differences in mCG levels between lens epithelial and fiber cells. Changes in mCG levels were inversely correlated with the differentiation state-specific expression of 1285 genes preferentially expressed in either lens fiber or lens epithelial cells (Pearson correlation r = − 0.37, p < 1 × 10^–42). mCG levels were inversely correlated with chromatin accessibility determined by assay for transposase-accessible sequencing (ATAC-seq) (Pearson correlation r = − 0.86, p < 1 × 10^–300). Many of the genes exhibiting altered regions of DNA methylation, chromatin accessibility and gene expression levels in fiber cells relative to epithelial cells are associated with lens fiber cell structure, homeostasis and transparency. These include lens crystallins (CRYBA4, CRYBB1, CRYGN, CRYBB2), lens beaded filament proteins (BFSP1, BFSP2), transcription factors (HSF4, SOX2, HIF1A), and Notch signaling pathway members (NOTCH1, NOTCH2, HEY1, HES5). Analysis of regions exhibiting cell-type specific alterations in DNA methylation revealed an overrepresentation of consensus sequences of multiple transcription factors known to play key roles in lens cell differentiation including HIF1A, SOX2, and the MAF family of transcription factors.

Conclusions

Collectively, these results link DNA methylation with control of chromatin accessibility and gene expression changes required for eye lens differentiation. The results also point to a role for DNA methylation in the regulation of transcription factors previously identified to be important for lens cell differentiation.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13072-022-00440-z.

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.

Aged Nrf2-Null Mice Develop All Major Types of Age-Related Cataracts

Purpose

Age-related cataracts affect the majority of older adults and are a leading cause of blindness worldwide. Treatments that delay cataract onset or severity have the potential to delay cataract surgery, but require relevant animal models that recapitulate the major types of cataracts for their development. Unfortunately, few such models are available. Here, we report the lens phenotypes of aged mice lacking the critical antioxidant transcription factor Nfe2l2 (designated as Nrf2 −/−).

Methods

Three independent cohorts of Nrf2 −/− and wild-type C57BL/6J mice were evaluated for cataracts using combinations of slit lamp imaging, photography of freshly dissected lenses, and histology. Mice were fed high glycemic diets, low glycemic diets, regular chow ad libitum, or regular chow with 30% caloric restriction.

Results

Nrf2 −/− mice developed significant opacities between 11 and 15 months and developed advanced cortical, posterior subcapsular, anterior subcapsular, and nuclear cataracts. Cataracts occurred similarly in male mice fed high or low glycemic diets, and were also observed in 21-month male and female Nrf2 −/− mice fed ad libitum or 30% caloric restriction. Histological observation of 18-month cataractous lenses revealed significant disruption to fiber cell architecture and the retention of nuclei throughout the cortical region of the lens. However, fiber cell denucleation and initiation of lens differentiation was normal at birth, with the first abnormalities observed at 3 months.

Conclusions

Nrf2 −/− mice offer a tool to understand how defective antioxidant signaling causes multiple forms of cataract and may be useful for screening drugs to prevent or delay cataractogenesis in susceptible adults.

Transcriptomics Analysis of Lens from Patients with Posterior Subcapsular Congenital Cataract

To gain insight into the aetiology of posterior subcapsular congenital cataract from the perspective of transcriptional changes, we conducted an mRNA sequencing analysis of the lenses in posterior subcapsular congenital cataract patients and in normal children. There were 1533 differentially expressed genes from 19,072 genes in the lens epithelial cells of the posterior subcapsular congenital cataract patients compared to in the normal controls at a cut-off criteria of |log2 fold change| of >1 and a p-value of <0.05, including 847 downregulated genes and 686 upregulated genes. To further narrow down the DEGs, we utilised the stricter criteria of |log2 fold change| of >1 and an FDR value of <0.05, and we identified 551 DEGs, including 97 upregulated genes and 454 downregulated genes. This study also identified 1263 differentially expressed genes of the 18,755 genes in lens cortex and nuclear fibres, including 646 downregulated genes and 617 upregulated genes. The downregulated genes in epithelial cells were significantly enriched in the structural constituent of lenses, lens development and lens fibre cell differentiation. After filtering the DEGs using the databases iSyTE and Cat-Map, several high-priority candidate genes related to posterior subcapsular congenital cataract such as GRIFIN, HTRA1 and DAPL1 were identified. The findings of our study may provide a deeper understanding of the mechanisms of posterior subcapsular congenital cataract and help in the prevention and treatment of this disease.

A large multiethnic GWAS meta-analysis of cataract identifies new risk loci and sex-specific effects

Cataract is the leading cause of blindness among the elderly worldwide and cataract surgery is one of the most common operations performed in the United States. As the genetic etiology of cataract formation remains unclear, we conducted a multiethnic genome-wide association meta-analysis, combining results from the GERA and UK Biobank cohorts, and tested for replication in the 23andMe research cohort. We report 54 genome-wide significant loci, 37 of which were novel. Sex-stratified analyses identified CASP7 as an additional novel locus specific to women. We show that genes within or near 80% of the cataract-associated loci are significantly expressed and/or enriched-expressed in the mouse lens across various spatiotemporal stages as per iSyTE analysis. Furthermore, iSyTE shows 32 candidate genes in the associated loci have altered gene expression in 9 different gene perturbation mouse models of lens defects/cataract, suggesting their relevance to lens biology. Our work provides further insight into the complex genetic architecture of cataract susceptibility.

The aging mouse lens transcriptome

Age is a major risk factor for cataract (ARC). However, the influence of aging on the lens transcriptome is under studied. Lens epithelial (LEC) and fiber cells (LFC) were isolated from young (3 month old) and aged (24 month old) C57BL/6J mice, and the transcriptome elucidated via RNAseq. EdgeR estimated differential gene expression in pairwise contrasts, and Advaita's Ipathway guide and custom R scripts were used to evaluate the potential biological significance of differentially expressed genes (DEGs). This analysis revealed age-dependent decreases in lens differentiation marker expression in both LECs and LFCs, with gamma crystallin transcripts downregulating nearly 50 fold in aged LFCs. The expression of the transcription factors Hsf4 and Maf, which are known to activate lens fiber cell preferred genes, are downregulated, while FoxE3, which represses gamma crystallin expression, is upregulated in aged fibers. Aged LECs upregulate genes controlling the immune response, complement pathways, and cellular stress responses, including glutathione peroxidase 3 (Gpx3). Aged LFCs exhibit broad changes in the expression of genes regulating cell communication, and upregulate genes involved in antigen processing/presentation and cholesterol metabolism, while changes in the expression of mitochondrial respiratory chain genes are consistent with mitochondrial stress, including upregulation of NDufa4l2, which encodes an alternate electron transport chain protein. However, age did not profoundly affect the response of LECs to injury as both young and aged LECs upregulate inflammatory gene signatures at 24 h post injury to similar extents. These RNAseq profiles provide a rich data set that can be mined to understand the genetic regulation of lens aging and how this impinges on the pathophysiology of age related cataract.

The master transcription factor SOX2, mutated in anophthalmia/microphthalmia, is post-transcriptionally regulated by the conserved RNA-binding protein RBM24 in vertebrate eye development

Mutations in the key transcription factor, SOX2, alone account for 20% of anophthalmia (no eye) and microphthalmia (small eye) birth defects in humans-yet its regulation is not well understood, especially on the post-transcription level. We report the unprecedented finding that the conserved RNA-binding motif protein, RBM24, positively controls Sox2 mRNA stability and is necessary for optimal SOX2 mRNA and protein levels in development, perturbation of which causes ocular defects, including microphthalmia and anophthalmia. RNA immunoprecipitation assay indicates that RBM24 protein interacts with Sox2 mRNA in mouse embryonic eye tissue. and electrophoretic mobility shift assay shows that RBM24 directly binds to the Sox2 mRNA 3'UTR, which is dependent on AU-rich elements (ARE) present in the Sox2 mRNA 3'UTR. Further, we demonstrate that Sox2 3'UTR AREs are necessary for RBM24-based elevation of Sox2 mRNA half-life. We find that this novel RBM24-Sox2 regulatory module is essential for early eye development in vertebrates. We show that Rbm24-targeted deletion using a constitutive CMV-driven Cre in mouse, and rbm24a-CRISPR/Cas9-targeted mutation or morpholino knockdown in zebrafish, results in Sox2 downregulation and causes the developmental defects anophthalmia or microphthalmia, similar to human SOX2-deficiency defects. We further show that Rbm24 deficiency leads to apoptotic defects in mouse ocular tissue and downregulation of eye development markers Lhx2, Pax6, Jag1, E-cadherin and gamma-crystallins. These data highlight the exquisite specificity that conserved RNA-binding proteins like RBM24 mediate in the post-transcriptional control of key transcription factors, namely, SOX2, associated with organogenesis and human developmental defects.

Single cell transcriptomics of the developing zebrafish lens and identification of putative controllers of lens development

The vertebrate lens is a valuable model system for investigating the gene expression changes that coordinate tissue differentiation due to its inclusion of two spatially separated cell types, the outer epithelial cells and the deeper denucleated fiber cells that they support. Zebrafish are a useful model system for studying lens development given the organ's rapid development in the first several days of life in an accessible, transparent embryo. While we have strong foundational knowledge of the diverse lens crystallin proteins and the basic gene regulatory networks controlling lens development, no study has detailed gene expression in a vertebrate lens at single cell resolution. Here we report an atlas of lens gene expression in zebrafish embryos and larvae at single cell resolution through five days of development, identifying a number of novel putative regulators of lens development. Our data address open questions about the temperospatial expression of α-crystallins during lens development that will support future studies of their function and provide the first detailed view of β- and γ-crystallin expression in and outside the lens. We describe divergent expression in transcription factor genes that occur as paralog pairs in the zebrafish. Finally, we examine the expression dynamics of cytoskeletal, membrane associated, RNA-binding, and transcription factor genes, identifying a number of novel patterns. Overall these data provide a foundation for identifying and characterizing lens developmental regulatory mechanisms and revealing targets for future functional studies with potential therapeutic impact.

Genome-Wide Analysis of Differentially Expressed miRNAs and Their Associated Regulatory Networks in Lenses Deficient for the Congenital Cataract-Linked Tudor Domain Containing Protein TDRD7

Mutations/deficiency of TDRD7, encoding a tudor domain protein involved in post-transcriptional gene expression control, causes early onset cataract in humans. While Tdrd7 is implicated in the control of key lens mRNAs, the impact of Tdrd7 deficiency on microRNAs (miRNAs) and how this contributes to transcriptome misexpression and to cataracts, is undefined. We address this critical knowledge-gap by investigating Tdrd7-targeted knockout (Tdrd7-/-) mice that exhibit fully penetrant juvenile cataracts. We performed Affymetrix miRNA 3.0 microarray analysis on Tdrd7-/- mouse lenses at postnatal day (P) 4, a stage preceding cataract formation. This analysis identifies 22 miRNAs [14 over-expressed (miR-15a, miR-19a, miR-138, miR-328, miR-339, miR-345, miR-378b, miR-384, miR-467a, miR-1224, miR-1935, miR-1946a, miR-3102, miR-3107), 8 reduced (let-7b, miR-34c, miR-298, miR-382, miR-409, miR-1198, miR-1947, miR-3092)] to be significantly misexpressed (fold-change ≥ ± 1.2, p-value < 0.05) in Tdrd7-/- lenses. To understand how these misexpressed miRNAs impact Tdrd7-/- cataract, we predicted their mRNA targets and examined their misexpression upon Tdrd7-deficiency by performing comparative transcriptomics analysis on P4 and P30 Tdrd7-/- lens. To prioritize these target mRNAs, we used various stringency filters (e.g., fold-change in Tdrd7-/- lens, iSyTE-based lens-enriched expression) and identified 98 reduced and 89 elevated mRNA targets for overexpressed and reduced miRNAs, respectively, which were classified as “top-priority” “high-priority,” and “promising” candidates. For Tdrd7-/- lens overexpressed miRNAs, this approach identified 18 top-priority reduced target mRNAs: Alad, Ankrd46, Ceacam10, Dgat2, Ednrb, H2-Eb1, Klhl22, Lin7a, Loxl1, Lpin1, Npc1, Olfm1, Ppm1e, Ppp1r1a, Rgs8, Shisa4, Snx22 and Wnk2. Majority of these targets were also altered in other gene-specific perturbation mouse models (e.g., Brg1, E2f1/E2f2/E2f3, Foxe3, Hsf4, Klf4, Mafg/Mafk, Notch) of lens defects/cataract, suggesting their importance to lens biology. Gene ontology (GO) provided further insight into their relevance to lens pathology. For example, the Tdrd7-deficient lens capsule defect may be explained by reduced mRNA targets (e.g., Col4a3, Loxl1, Timp2, Timp3) associated with “basement membrane”. GO analysis also identified new genes (e.g., Casz1, Rasgrp1) recently linked to lens biology/pathology. Together, these analyses define a new Tdrd7-downstream miRNA-mRNA network, in turn, uncovering several new mRNA targets and their associated pathways relevant to lens biology and offering molecular insights into the pathology of congenital cataract.

iTRAQ-Based Proteomic Profiling of Potential Biomarkers in Rat Serum for Uranium Tailing Suspension Intratracheal Instillation

Protection against low-dose ionizing radiation is of great significance. Uranium tailings are formed as a byproduct of uranium mining and a potential risk to organisms. In this study, we identified potential biomarkers associated with exposure to low-dose radiation from uranium tailings. We established a Wistar rat model of low dose rate irradiation by intratracheal instillation of a uranium tailing suspension. We observed pathological changes in the liver, lung, and kidney tissues of the rats. Using isobaric tags for relative and absolute quantification, we screened 17 common differentially expressed proteins in three dose groups. We chose alpha-1 antiproteinase (Serpina1), keratin 17 (Krt17), and aldehyde dehydrogenase (Aldh3a1) for further investigation. Our data showed that expression of Serpina1, Krt17, and Aldh3a1 had changed after the intratracheal instillation in rats, which may be potential biomarkers for uranium tailing low-dose irradiation. However, the underlying mechanisms require further investigation.

The Tudor-domain protein TDRD7, mutated in congenital cataract, controls the heat shock protein HSPB1 (HSP27) and lens fiber cell morphology

Mutations of the RNA granule component TDRD7 (OMIM: 611258) cause pediatric cataract. We applied an integrated approach to uncover the molecular pathology of cataract in Tdrd7-/- mice. Early postnatal Tdrd7-/- animals precipitously develop cataract suggesting a global-level breakdown/misregulation of key cellular processes. High-throughput RNA sequencing integrated with iSyTE-bioinformatics analysis identified the molecular chaperone and cytoskeletal modulator, HSPB1, among high-priority downregulated candidates in Tdrd7-/- lens. A protein fluorescence two-dimensional difference in-gel electrophoresis (2D-DIGE)-coupled mass spectrometry screen also identified HSPB1 downregulation, offering independent support for its importance to Tdrd7-/- cataractogenesis. Lens fiber cells normally undergo nuclear degradation for transparency, posing a challenge: how is their cell morphology, also critical for transparency, controlled post-nuclear degradation? HSPB1 functions in cytoskeletal maintenance, and its reduction in Tdrd7-/- lens precedes cataract, suggesting cytoskeletal defects may contribute to Tdrd7-/- cataract. In agreement, scanning electron microscopy (SEM) revealed abnormal fiber cell morphology in Tdrd7-/- lenses. Further, abnormal phalloidin and wheat germ agglutinin (WGA) staining of Tdrd7-/- fiber cells, particularly those exhibiting nuclear degradation, reveals distinct regulatory mechanisms control F-actin cytoskeletal and/or membrane maintenance in post-organelle degradation maturation stage fiber cells. Indeed, RNA immunoprecipitation identified Hspb1 mRNA in wild-type lens lysate TDRD7-pulldowns, and single-molecule RNA imaging showed co-localization of TDRD7 protein with cytoplasmic Hspb1 mRNA in differentiating fiber cells, suggesting that TDRD7-ribonucleoprotein complexes may be involved in optimal buildup of key factors. Finally, Hspb1 knockdown in Xenopus causes eye/lens defects. Together, these data uncover TDRD7's novel upstream role in elevation of stress-responsive chaperones for cytoskeletal maintenance in post-nuclear degradation lens fiber cells, perturbation of which causes early-onset cataracts.

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.

Transcriptomic analysis and novel insights into lens fibre cell differentiation regulated by Gata3

Gata3 is a DNA-binding transcription factor involved in cellular differentiation in a variety of tissues including inner ear, hair follicle, kidney, mammary gland and T-cells. In a previous study in 2009, Maeda et al. (Dev. Dyn.238, 2280–2291; doi:10.1002/dvdy.22035) found that Gata3 mutants could be rescued from midgestational lethality by the expression of a Gata3 transgene in sympathoadrenal neuroendocrine cells. The rescued embryos clearly showed multiple defects in lens fibre cell differentiation. To determine whether these defects were truly due to the loss of Gata3 expression in the lens, we generated a lens-specific Gata3 loss-of-function model. Analogous to the previous findings, our Gata3 null embryos showed abnormal regulation of cell cycle exit during lens fibre cell differentiation, marked by reduction in the expression of the cyclin-dependent kinase inhibitors Cdkn1b/p27 and Cdkn1c/p57, and the retention of nuclei accompanied by downregulation of Dnase IIβ. Comparisons of transcriptomes between control and mutated lenses by RNA-Seq revealed dysregulation of lens-specific crystallin genes and intermediate filament protein Bfsp2. Both Cdkn1b/p27 and Cdkn1c/p57 loci are occupied in vivo by Gata3, as well as Prox1 and c-Jun, in lens chromatin. Collectively, our studies suggest that Gata3 regulates lens differentiation through the direct regulation of the Cdkn1b/p27and Cdkn1c/p57 expression, and the direct/or indirect transcriptional control of Bfsp2 and Dnase IIβ.

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

While the bioinformatics resource-tool iSyTE (integrated Systems Tool for Eye gene discovery) effectively identifies human cataract-associated genes, it is currently based on just transcriptome data, and thus, it is necessary to include protein-level information to gain greater confidence in gene prioritization. Here, we expand iSyTE through development of a novel proteome-based resource on the lens and demonstrate its utility in cataract gene discovery. We applied high-throughput tandem mass spectrometry (MS/MS) to generate a global protein expression profile of mouse lens at embryonic day (E)14.5, which identified 2371 lens-expressed proteins. A major challenge of high-throughput expression profiling is identification of high-priority candidates among the thousands of expressed proteins. To address this problem, we generated new MS/MS proteome data on mouse whole embryonic body (WB). WB proteome was then used as a reference dataset for performing "in silico WB-subtraction" comparative analysis with the lens proteome, which effectively identified 422 proteins with lens-enriched expression at ≥ 2.5 average spectral counts, ≥ 2.0 fold enrichment (FDR < 0.01) cut-off. These top 20% candidates represent a rich pool of high-priority proteins in the lens including known human cataract-linked genes and many new potential regulators of lens development and homeostasis. This rich information is made publicly accessible through iSyTE (https://research.bioinformatics.udel.edu/iSyTE/), which enables user-friendly visualization of promising candidates, thus making iSyTE a comprehensive tool for cataract gene discovery.

Molecular characterization of the human lens epithelium-derived cell line SRA01/04

Cataract-associated gene discovery in human and animal models have informed on key aspects of human lens development, homeostasis and pathology. Additionally, in vitro models such as the culture of permanent human lens epithelium-derived cell lines (LECs) have also been utilized to understand the molecular biology of lens cells. However, these resources remain uncharacterized, specifically regarding their global gene expression and suitability to model lens cell biology. Therefore, we sought to molecularly characterize gene expression in the human LEC, SRA01/04, which is commonly used in lens studies. We first performed short tandem repeat (STR) analysis and validated SRA01/04 LEC for its human origin, as recommended by the eye research community. Next, we used Illumina HumanHT-12 v3.0 Expression BeadChip arrays to gain insights into the global gene expression profile of SRA01/04. Comparative analysis of SRA01/04 microarray data was performed using other resources such as the lens expression database iSyTE (integrated Systems Tool for Eye gene discovery), the cataract gene database Cat-Map and the published lens literature. This analysis showed that SRA01/04 significantly expresses >40% of the top iSyTE lens-enriched genes (313 out of 749) across different developmental stages. Further, SRA01/04 also significantly expresses ~53% (168 out of 318) of cataract-associated genes in Cat-Map. We also performed comparative gene expression analysis between SRA01/04 cells and the previously validated mouse LEC 21EM15. To gain insight into whether SRA01/04 reflects epithelial or fiber cell characteristics, we compared its gene expression profile to previously reported differentially expressed genes in isolated mouse lens epithelial and fiber cells. This analysis suggests that SRA01/04 has reduced expression of several fiber cell-enriched genes. In agreement with these findings, cell culture analysis demonstrates that SRA01/04 has reduced potential to initiate spontaneous lentoid body formation compared to 21EM15 cells. Next, to independently validate SRA01/04 microarray gene expression, we subjected several candidate genes to RT-PCR and RT-qPCR assays. This analysis demonstrates that SRA01/04 supports expression of many key genes associated with lens development and cataract, including CRYAB, CRYBB2, CRYGS, DKK3, EPHA2, ETV5, GJA1, HSPB1, INPPL1, ITGB1, PAX6, PVRL3, SFRP1, SPARC, TDRD7, and VIM, among others, and therefore can be relevant for understanding the mechanistic basis of these factors. At the same time, SRA01/04 cells do not exhibit robust expression of several genes known to be important to lens biology and cataract such as ALDH1A1, COL4A6, CP, CRYBA4, FOXE3, HMX1, HSF4, MAF, MEIS1, PITX3, PRX, SIX3, and TRPM3, among many others. Therefore, the present study offers a rich transcript-level resource for case-by-case evaluation of the potential advantages and limitations of SRA01/04 cells prior to their use in downstream investigations. In sum, these data show that the human LEC, SRA01/04, exhibits lens epithelial cell-like character reflected in the expression of several lens-enriched and cataract-associated genes, and therefore can be considered as a useful in vitro resource when combined with in vivo studies to gain insight into specific aspects of human lens epithelial cells.

iSyTE 2.0: a database for expression-based gene discovery in the eye

Although successful in identifying new cataract-linked genes, the previous version of the database iSyTE (integrated Systems Tool for Eye gene discovery) was based on expression information on just three mouse lens stages and was functionally limited to visualization by only UCSC-Genome Browser tracks. To increase its efficacy, here we provide an enhanced iSyTE version 2.0 (URL: http://research.bioinformatics.udel.edu/iSyTE) based on well-curated, comprehensive genome-level lens expression data as a one-stop portal for the effective visualization and analysis of candidate genes in lens development and disease. iSyTE 2.0 includes all publicly available lens Affymetrix and Illumina microarray datasets representing a broad range of embryonic and postnatal stages from wild-type and specific gene-perturbation mouse mutants with eye defects. Further, we developed a new user-friendly web interface for direct access and cogent visualization of the curated expression data, which supports convenient searches and a range of downstream analyses. The utility of these new iSyTE 2.0 features is illustrated through examples of established genes associated with lens development and pathobiology, which serve as tutorials for its application by the end-user. iSyTE 2.0 will facilitate the prioritization of eye development and disease-linked candidate genes in studies involving transcriptomics or next-generation sequencing data, linkage analysis and GWAS approaches.

Profiling of chromatin accessibility and identification of general cis-regulatory mechanisms that control two ocular lens differentiation pathways

BACKGROUND

Promoters and enhancers are cis-regulatory DNA sequences that control specificity and quantity of transcription. Both are rich on clusters of cis-acting sites that interact with sequence-specific DNA-binding transcription factors (TFs). At the level of chromatin, these regions display increased nuclease sensitivity, reduced nucleosome density, including nucleosome-free regions, and specific combinations of posttranslational modifications of core histone proteins. Together, "open" and "closed" chromatins represent transcriptionally active and repressed states of individual genes, respectively. Cellular differentiation is marked by changes in local chromatin structure. Lens morphogenesis, regulated by TF Pax6, includes differentiation of epithelial precursor cells into lens fibers in parallel with differentiation of epithelial precursors into the mature lens epithelium.

RESULTS

Using ATAC-seq, we investigated dynamics of chromatin changes during mouse lens fibers and epithelium differentiation. Tissue-specific features of these processes are demonstrated via comparative studies of embryonic stem cells, forebrain, and liver chromatins. Unbiased analysis reveals cis-regulatory logic of lens differentiation through known (e.g., AP-1, Ets, Hsf4, Maf, and Pax6 sites) and novel (e.g., CTCF, Tead, and NF1) motifs. Twenty-six DNA-binding TFs, recognizing these cis-motifs, are markedly up-regulated in differentiating lens fibers. As specific examples, our ATAC-seq data uncovered both the regulatory regions and TF binding motifs in Foxe3, Prox1, and Mip loci that are consistent with previous, though incomplete, experimental data. A cross-examination of Pax6 binding with ATAC-seq data demonstrated that Pax6 bound to both open (H3K27ac and P300-enriched) and closed chromatin domains in lens and forebrain.

CONCLUSIONS

Our study has generated the first lens chromatin accessibility maps that support a general model of stage-specific chromatin changes associated with transcriptional activities of batteries of genes required for lens fiber cell formation. Analysis of active (or open) promoters and enhancers reveals important cis-DNA motifs that establish the molecular foundation for temporally and spatially regulated gene expression in lens. Together, our data and models open new avenues for the field to conduct mechanistic studies of transcriptional control regions, reconstruction of gene regulatory networks that govern lens morphogenesis, and identification of cataract-causing mutations in noncoding sequences.

Photoreceptor Survival Is Regulated by GSTO1-1 in the Degenerating Retina

Purpose

Glutathione-S-transferase omega 1-1 (GSTO1-1) is a cytosolic glutathione transferase enzyme, involved in glutathionylation, toll-like receptor signaling, and calcium channel regulation. GSTO1-1 dysregulation has been implicated in oxidative stress and inflammation, and contributes to the pathogenesis of several diseases and neurological disorders; however, its role in retinal degenerations is unknown. The aim of this study was to investigate the role of GSTO1-1 in modulating oxidative stress and consequent inflammation in the normal and degenerating retina.

Methods

The role of GSTO1-1 in retinal degenerations was explored by using Gsto1-/- mice in a model of retinal degeneration. The expression and localization of GSTO1-1 were investigated with immunohistochemistry and Western blot. Changes in the expression of inflammatory (Ccl2, Il-1β, and C3) and oxidative stress (Nox1, Sod2, Gpx3, Hmox1, Nrf2, and Nqo1) genes were investigated via quantitative real-time polymerase chain reaction. Retinal function in Gsto1-/- mice was investigated by using electroretinography.

Results

GSTO1-1 was localized to the inner segment of cone photoreceptors in the retina. Gsto1-/- photo-oxidative damage (PD) mice had decreased photoreceptor cell death as well as decreased expression of inflammatory (Ccl2, Il-1β, and C3) markers and oxidative stress marker Nqo1. Further, retinal function in the Gsto1-/- PD mice was increased as compared to wild-type PD mice.

Conclusions

These results indicate that GSTO1-1 is required for inflammatory-mediated photoreceptor death in retinal degenerations. Targeting GSTO1-1 may be a useful strategy to reduce oxidative stress and inflammation and ameliorate photoreceptor loss, slowing the progression of retinal degenerations.

A novel MAF missense mutation leads to congenital nuclear cataract by impacting the transactivation of crystallin and noncrystallin genes

The transcription factor v-maf avain musculoaponeurotic fibrosarcoma oncogene homolog (MAF) plays an important role in lens development. It contains a unique extended homology region (EHR) in the DNA binding domain. MAF mutations are associated with phenotypically distinct forms of congenital cataract and show different effects on the transactivation of target genes. Mutations in the MAF EHR region were rarely reported and their corresponding phenotype and impact on target genes' transactivation were not evaluated. A three- generation Chinese family with congenital cataract was recruited. The patients in the family present non-syndromic congenital nuclear and lamellar opacities. A novel MAF mutation (c.812 T > A, p.Val271Glu) was identified by targeted next-generation sequencing. The mutation is in highly conserved EHR region of MAF and co-segregates with the cataract in the family. It is predicted to be pathogenic by multiple algorithms and is absent in a control population. Dual luciferase activity assay shows the mutation significantly impair the transcriptional activity of four crystallin genes (CRYAA, CRYBA4, CRYBA1, and CRYGA) and two non-crystallin genes (HMOX1 and KDELR2). Herein, we report a novel missense mutation in the MAF EHR region of the DNA binding domain in a family with congenital cataract. The mutation is associated with non-syndromic bilateral nuclear cataract and impacts the transactivation of cataract associated genes involved in lens structure and stress response.

RNA sequencing-based transcriptomic profiles of embryonic lens development for cataract gene discovery

Isolated or syndromic congenital cataracts are heterogeneous developmental defects, making the identification of the associated genes challenging. In the past, mouse lens expression microarrays have been successfully applied in bioinformatics tools (e.g., iSyTE) to facilitate human cataract-associated gene discovery. To develop a new resource for geneticists, we report high-throughput RNA sequencing (RNA-seq) profiles of mouse lens at key embryonic stages (E)10.5 (lens pit), E12.5 (primary fiber cell differentiation), E14.5 and E16.5 (secondary fiber cell differentiation). These stages capture important events as the lens develops from an invaginating placode into a transparent tissue. Previously, in silico whole-embryo body (WB)-subtraction-based "lens-enriched" expression has been effective in prioritizing cataract-linked genes. To apply an analogous approach, we generated new mouse WB RNA-seq datasets and show that in silico WB subtraction of lens RNA-seq datasets successfully identifies key genes based on lens-enriched expression. At ≥2 counts-per-million expression, ≥1.5 log₂ fold-enrichment (p < 0.05) cutoff, E10.5 lens exhibits 1401 enriched genes (17% lens-expressed genes), E12.5 lens exhibits 1937 enriched genes (22% lens-expressed genes), E14.5 lens exhibits 2514 enriched genes (31% lens-expressed genes), and E16.5 lens exhibits 2745 enriched genes (34% lens-expressed genes). Biological pathway analysis identified genes associated with lens development, transcription regulation and signaling pathways, among other functional groups. Furthermore, these new RNA-seq data confirmed high expression of established cataract-linked genes and identified new potential regulators in the lens. Finally, we developed new lens stage-specific UCSC Genome Brower annotation tracks and made these publicly accessible through iSyTE ( https://research.bioinformatics.udel.edu/iSyTE/ ) for user-friendly visualization of lens gene expression/enrichment to prioritize genes from high-throughput data from cataract cases.

A comprehensive spatial-temporal transcriptomic analysis of differentiating nascent mouse lens epithelial and fiber cells

Elucidation of both the molecular composition and organization of the ocular lens is a prerequisite to understand its development, function, pathology, regenerative capacity, as well as to model lens development and disease using in vitro differentiation of pluripotent stem cells. Lens is comprised of the anterior lens epithelium and posterior lens fibers, which form the bulk of the lens. Lens fibers differentiate from lens epithelial cells through cell cycle exit-coupled differentiation that includes cellular elongation, accumulation of crystallins, cytoskeleton and membrane remodeling, and degradation of organelles within the central region of the lens. Here, we profiled spatiotemporal expression dynamics of both mRNAs and non-coding RNAs from microdissected mouse nascent lens epithelium and lens fibers at four developmental time points (embryonic [E] day 14.5, E16.5, E18.5, and P0.5) by RNA-seq. During this critical time window, multiple complex biosynthetic and catabolic processes generate the molecular and structural foundation for lens transparency. Throughout this developmental window, 3544 and 3518 genes show consistently and significantly greater expression in the nascent lens epithelium and fibers, respectively. Comprehensive data analysis confirmed major roles of FGF-MAPK, Wnt/β-catenin, PI3K/AKT, TGF-β, and BMP signaling pathways and revealed significant novel contributions of mTOR, EIF2, EIF4, and p70S6K signaling in lens formation. Unbiased motif analysis within promoter regions of these genes with consistent expression changes between epithelium and fiber cells revealed an enrichment for both established (e.g. E2Fs, Etv5, Hsf4, c-Maf, MafG, MafK, N-Myc, and Pax6) transcription factors and a number of novel regulators of lens formation, such as Arntl2, Dmrta2, Stat5a, Stat5b, and Tulp3. In conclusion, the present RNA-seq data serves as a comprehensive reference resource for deciphering molecular principles of normal mammalian lens differentiation, mapping a full spectrum of signaling pathways and DNA-binding transcription factors operating in both lens compartments, and predicting novel pathways required to establish lens transparency.

Mutation update of transcription factor genes FOXE3, HSF4, MAF, and PITX3 causing cataracts and other developmental ocular defects

Mutations in the transcription factor genes FOXE3, HSF4, MAF, and PITX3 cause congenital lens defects including cataracts that may be accompanied by defects in other components of the eye or in nonocular tissues. We comprehensively describe here all the variants in FOXE3, HSF4, MAF, and PITX3 genes linked to human developmental defects. A total of 52 variants for FOXE3, 18 variants for HSF4, 20 variants for MAF, and 19 variants for PITX3 identified so far in isolated cases or within families are documented. This effort reveals FOXE3, HSF4, MAF, and PITX3 to have 33, 16, 18, and 7 unique causal mutations, respectively. Loss-of-function mutant animals for these genes have served to model the pathobiology of the associated human defects, and we discuss the currently known molecular function of these genes, particularly with emphasis on their role in ocular development. Finally, we make the detailed FOXE3, HSF4, MAF, and PITX3 variant information available in the Leiden Online Variation Database (LOVD) platform at https://www.LOVD.nl/FOXE3, https://www.LOVD.nl/HSF4, https://www.LOVD.nl/MAF, and https://www.LOVD.nl/PITX3. Thus, this article informs on key variants in transcription factor genes linked to cataract, aphakia, corneal opacity, glaucoma, microcornea, microphthalmia, anterior segment mesenchymal dysgenesis, and Ayme-Gripp syndrome, and facilitates their access through Web-based databases.

The RNA-binding protein Celf1 post-transcriptionally regulates p27Kip1 and Dnase2b to control fiber cell nuclear degradation in lens development

Opacification of the ocular lens, termed cataract, is a common cause of blindness. To become transparent, lens fiber cells undergo degradation of their organelles, including their nuclei, presenting a fundamental question: does signaling/transcription sufficiently explain differentiation of cells progressing toward compromised transcriptional potential? We report that a conserved RNA-binding protein Celf1 post-transcriptionally controls key genes to regulate lens fiber cell differentiation. Celf1-targeted knockout mice and celf1-knockdown zebrafish and Xenopus morphants have severe eye defects/cataract. Celf1 spatiotemporally down-regulates the cyclin-dependent kinase (Cdk) inhibitor p27Kip1 by interacting with its 5' UTR and mediating translation inhibition. Celf1 deficiency causes ectopic up-regulation of p21Cip1. Further, Celf1 directly binds to the mRNA of the nuclease Dnase2b to maintain its high levels. Together these events are necessary for Cdk1-mediated lamin A/C phosphorylation to initiate nuclear envelope breakdown and DNA degradation in fiber cells. Moreover, Celf1 controls alternative splicing of the membrane-organization factor beta-spectrin and regulates F-actin-crosslinking factor Actn2 mRNA levels, thereby controlling fiber cell morphology. Thus, we illustrate new Celf1-regulated molecular mechanisms in lens development, suggesting that post-transcriptional regulatory RNA-binding proteins have evolved conserved functions to control vertebrate oculogenesis.

Identification of OAF and PVRL1 as candidate genes for an ocular anomaly characterized by Peters anomaly type 2 and ectopia lentis

Keratolenticular dysgenesis (KLD) and ectopia lentis are congenital eye defects. The aim of this study is the identification of molecular genetic alterations responsible for those ocular anomalies with neurologic impairment in an individual with a de novo balanced chromosome translocation t(11;18)(q23.3;q11.2)dn. Disruption of OAF, the human orthologue of the Drosophila oaf, by the 11q23.3 breakpoint results in reduced expression of this transcriptional regulator. Furthermore, four most likely nonfunctional chimeric transcripts comprising up to OAF exon 3, derived from the der(11) allele, have also been identified. This locus has been implicated by publicly available genome-wide association data in corneal disease and corneal topography. The expression of the poliovirus receptor-related 1(PVRL1) or nectin cell adhesion molecule 1 (NECTIN1), a paralogue of nectin cell adhesion molecule 3 (PVRL3) associated with congenital ocular defects, situated 500 kb upstream from 11q23.3 breakpoint, is increased. The 18q11.2 breakpoint is localized between cutaneous T-cell lymphoma-associated antigen 1(CTAGE1) and retinoblastoma binding protein 8 (RBBP8) genes. Genomic imbalance that could contribute to the observed phenotype was excluded. Analysis of gene expression datasets throughout normal murine ocular lens embryogenesis suggests that OAF expression is significantly enriched in the lens from early stages of development through adulthood, whereas PVRL1 is lens-enriched until E12.5 and then down-regulated. This contrasts with the observation that the proposita's lymphoblastoid cell lines exhibit low OAF and high PVRL1 expression as compared to control, which offers further support that the alterations described above are most likely responsible for the clinical phenotype. Finally, gene interaction topology data for PVRL1 also agree with our proposal that disruption of OAF by the translocation breakpoint and misregulation of PVRL1 due to a position effect contribute to the observed ocular and neurological phenotype.

miR-128 Is Implicated in Stress Responses by Targeting MAFG in Skeletal Muscle Cells

MAFG (v-Maf avian musculoaponeurotic fibrosarcoma oncogene homolog G) is a bZIP-type transcriptional regulator that belongs to the small MAF (sMAFs) protein family. By interacting with other bZIP transcription factors, sMAFs can form homo- and heterodimers governing either repressive or activating transcriptional functions. As heterodimeric partner of Nrf2, MAFG positively influences the ARE-dependent antioxidant/xenobiotic pathways, at least in condition of a correct MAFG:Nrf2 balance. MicroRNAs (miRs) participate to different regulatory networks being involved as fine-tuning regulators of gene expression. However, the connections between cellular surveillance to stresses mediated by MAFG:Nrf2 and miR regulations are not well understood. Here, we explored the impact of miR-128 in expression of genes related to stress response. Bioinformatic predictions coupled with functional analysis revealed the presence of miR-128 binding site in the 3′UTR of MAFG. Ectopic miR-128 expression correlated with reduced expression of endogenous MAFG-dependent genes and negatively affected ARE-mediated molecular phenotype based on Nrf2 activity. Indeed, miR-128 impairs redox-dependent pathways induced in response to oxidative stress. Moreover, in condition of hypoxia, MAFG induction correlated with reduced levels of miR-128. This lead to increased mRNA levels of HMOX-1 and x-CT for blunting stress. Overall, these findings identify MAFG as novel direct target of miR-128.

Signaling and Gene Regulatory Networks in Mammalian Lens Development

Ocular lens development represents an advantageous system in which to study regulatory mechanisms governing cell fate decisions, extracellular signaling, cell and tissue organization, and the underlying gene regulatory networks. Spatiotemporally regulated domains of BMP, FGF, and other signaling molecules in late gastrula-early neurula stage embryos generate the border region between the neural plate and non-neural ectoderm from which multiple cell types, including lens progenitor cells, emerge and undergo initial tissue formation. Extracellular signaling and DNA-binding transcription factors govern lens and optic cup morphogenesis. Pax6, c-Maf, Hsf4, Prox1, Sox1, and a few additional factors regulate the expression of the lens structural proteins, the crystallins. Extensive crosstalk between a diverse array of signaling pathways controls the complexity and order of lens morphogenetic processes and lens transparency.

Novel phenotypes and loci identified through clinical genomics approaches to pediatric cataract

Pediatric cataract is highly heterogeneous clinically and etiologically. While mostly isolated, cataract can be part of many multisystem disorders, further complicating the diagnostic process. In this study, we applied genomic tools in the form of a multi-gene panel as well as whole-exome sequencing on unselected cohort of pediatric cataract (166 patients from 74 families). Mutations in previously reported cataract genes were identified in 58% for a total of 43 mutations, including 15 that are novel. GEMIN4 was independently mutated in families with a syndrome of cataract, global developmental delay with or without renal involvement. We also highlight a recognizable syndrome that resembles galactosemia (a fulminant infantile liver disease with cataract) caused by biallelic mutations in CYP51A1. A founder mutation in RIC1 (KIAA1432) was identified in patients with cataract, brain atrophy, microcephaly with or without cleft lip and palate. For non-syndromic pediatric cataract, we map a novel locus in a multiplex consanguineous family on 4p15.32 where exome sequencing revealed a homozygous truncating mutation in TAPT1. We report two further candidates that are biallelically inactivated each in a single cataract family: TAF1A (cataract with global developmental delay) and WDR87 (non-syndromic cataract). In addition to positional mapping data, we use iSyTE developmental lens expression and gene-network analysis to corroborate the proposed link between the novel candidate genes and cataract. Our study expands the phenotypic, allelic and locus heterogeneity of pediatric cataract. The high diagnostic yield of clinical genomics supports the adoption of this approach in this patient group.

RNA-binding proteins in eye development and disease: implication of conserved RNA granule components

The molecular biology of metazoan eye development is an area of intense investigation. These efforts have led to the surprising recognition that although insect and vertebrate eyes have dramatically different structures, the orthologs or family members of several conserved transcription and signaling regulators such as Pax6, Six3, Prox1, and Bmp4 are commonly required for their development. In contrast, our understanding of posttranscriptional regulation in eye development and disease, particularly regarding the function of RNA-binding proteins (RBPs), is limited. We examine the present knowledge of RBPs in eye development in the insect model Drosophila as well as several vertebrate models such as fish, frog, chicken, and mouse. Interestingly, of the 42 RBPs that have been investigated for their expression or function in vertebrate eye development, 24 (~60%) are recognized in eukaryotic cells as components of RNA granules such as processing bodies, stress granules, or other specialized ribonucleoprotein (RNP) complexes. We discuss the distinct developmental and cellular events that may necessitate potential RBP/RNA granule-associated RNA regulon models to facilitate posttranscriptional control of gene expression in eye morphogenesis. In support of these hypotheses, three RBPs and RNP/RNA granule components Tdrd7, Caprin2, and Stau2 are linked to ocular developmental defects such as congenital cataract, Peters anomaly, and microphthalmia in human patients or animal models. We conclude by discussing the utility of interdisciplinary approaches such as the bioinformatics tool iSyTE (integrated Systems Tool for Eye gene discovery) to prioritize RBPs for deriving posttranscriptional regulatory networks in eye development and disease. WIREs RNA 2016, 7:527-557. doi: 10.1002/wrna.1355 For further resources related to this article, please visit the WIREs website.

Small Maf proteins (MafF, MafG, MafK): History, structure and function

The small Maf proteins (sMafs) are basic region leucine zipper (bZIP)-type transcription factors. The basic region of the Maf family is unique among the bZIP factors, and it contributes to the distinct DNA-binding mode of this class of proteins. MafF, MafG and MafK are the three vertebrate sMafs, and no functional differences have been observed among them in terms of their bZIP structures. sMafs form homodimers by themselves, and they form heterodimers with cap 'n' collar (CNC) proteins (p45 NF-E2, Nrf1, Nrf2, and Nrf3) and also with Bach proteins (Bach1 and Bach2). Because CNC and Bach proteins cannot bind to DNA as monomers, sMafs are indispensable partners that are required by CNC and Bach proteins to exert their functions. sMafs lack the transcriptional activation domain; hence, their homodimers act as transcriptional repressors. In contrast, sMafs participate in transcriptional activation or repression depending on their heterodimeric partner molecules and context. Mouse genetic analyses have revealed that various biological pathways are under the regulation of CNC-sMaf heterodimers. In this review, we summarize the history and current progress of sMaf studies in relation to their partners.

Structure, function and disease relevance of Omega-class glutathione transferases

The Omega-class cytosolic glutathione transferases (GSTs) have distinct structural and functional attributes that allow them to perform novel roles unrelated to the functions of other GSTs. Mammalian GSTO1-1 has been found to play a previously unappreciated role in the glutathionylation cycle that is emerging as significant mechanism regulating protein function. GSTO1-1-catalyzed glutathionylation or deglutathionylation of a key signaling protein may explain the requirement for catalytically active GSTO1-1 in LPS-stimulated pro-inflammatory signaling through the TLR4 receptor. The observation that ML175 a specific GSTO1-1 inhibitor can block LPS-stimulated inflammatory signaling has opened a new avenue for the development of novel anti-inflammatory drugs that could be useful in the treatment of toxic shock and other inflammatory disorders. The role of GSTO2-2 remains unclear. As a dehydroascorbate reductase, it could contribute to the maintenance of cellular redox balance and it is interesting to note that the GSTO2 N142D polymorphism has been associated with multiple diseases including Alzheimer's disease, Parkinson's disease, familial amyotrophic lateral sclerosis, chronic obstructive pulmonary disease, age-related cataract and breast cancer.

Systems biology of lens development: A paradigm for disease gene discovery in the eye

Over the past several decades, the biology of the developing lens has been investigated using molecular genetics-based approaches in various vertebrate model systems. These efforts, involving target gene knockouts or knockdowns, have led to major advances in our understanding of lens morphogenesis and the pathological basis of cataracts, as well as of other lens related eye defects. In particular, we now have a functional understanding of regulators such as Pax6, Six3, Sox2, Oct1 (Pou2f1), Meis1, Pnox1, Zeb2 (Sip1), Mab21l1, Foxe3, Tfap2a (Ap2-alpha), Pitx3, Sox11, Prox1, Sox1, c-Maf, Mafg, Mafk, Hsf4, Fgfrs, Bmp7, and Tdrd7 in this tissue. However, whether these individual regulators interact or their targets overlap, and the significance of such interactions during lens morphogenesis, is not well defined. The arrival of high-throughput approaches for gene expression profiling (microarrays, RNA-sequencing (RNA-seq), etc.), which can be coupled with chromatin immunoprecipitation (ChIP) or RNA immunoprecipitation (RIP) assays, along with improved computational resources and publically available datasets (e.g. those containing comprehensive protein-protein, protein-DNA information), presents new opportunities to advance our understanding of the lens tissue on a global systems level. Such systems-level knowledge will lead to the derivation of the underlying lens gene regulatory network (GRN), defined as a circuit map of the regulator-target interactions functional in lens development, which can be applied to expedite cataract gene discovery. In this review, we cover the various systems-level approaches such as microarrays, RNA-seq, and ChIP that are already being applied to lens studies and discuss strategies for assembling and interpreting these vast amounts of high-throughput information for effective dispersion to the scientific community. In particular, we discuss strategies for effective interpretation of this new information in the context of the rich knowledge obtained through the application of traditional single-gene focused experiments on the lens. Finally, we discuss our vision for integrating these diverse high-throughput datasets in a single web-based user-friendly tool iSyTE (integrated Systems Tool for Eye gene discovery) - a resource that is already proving effective in the identification and characterization of genes linked to lens development and cataract. We anticipate that application of a similar approach to other ocular tissues such as the retina and the cornea, and even other organ systems, will significantly impact disease gene discovery.

An integrative approach to analyze microarray datasets for prioritization of genes relevant to lens biology and disease

Microarray-based profiling represents an effective method to analyze cellular or tissue-specific gene expression on the genome-level. However, in comparative analyses between control and mutant samples, microarrays often identify a large number of differentially expressed genes, in turn making it challenging to isolate the select “high-priority candidates” that are most relevant to an observed mutant phenotype. Here, we describe an integrative approach for mouse mutant lens microarray gene expression analysis using publically accessible systems-level information such as wild-type mouse lens expression data in iSyTE (integrated Systems Tool for Eye gene discovery), protein–protein interaction data in public databases, gene ontology enrichment data, and transcription factor binding profile data. This strategy, when applied to small Maf Mafg −/−:Mafk +/− mouse lens microarray datasets (deposited in NCBI Gene Expression Omnibus database with accession number GSE65500) in Agrawal et al. 2015 [1], led to the effective prioritization of candidate genes linked to lens defects in these mutants. Indeed, from the original list of genes that are differentially expressed at ± 1.5-fold and p < 0.05 in Mafg −/−:Mafk +/− mutant lenses, this analysis led to the identification of thirty-six high-priority candidates, in turn reducing the number of genes for further study by approximately 1/3 of the total. Moreover, eight of these genes are linked to mammalian cataract in the published literature, validating the efficacy of this approach. Additionally, these high-priority candidates contribute valuable information for the assembly of a gene regulatory network in the lens. In sum, the pipeline outlined in this report represents an effective approach for initial as well as downstream microarray expression data analysis to identify genes important for lens biology and cataracts. We anticipate that this integrative strategy can be extended to prioritize phenotypically relevant candidate genes from microarray data in other cells and tissues.

Deficiency of the RNA binding protein caprin2 causes lens defects and features of Peters anomaly

BACKGROUND

It was recently demonstrated that deficiency of a conserved RNA binding protein (RBP) and RNA granule (RG) component Tdrd7 causes ocular defects including cataracts in human, mouse and chicken, indicating the importance of posttranscriptional regulation in eye development. Here we investigated the function of a second conserved RBP/RG component Caprin2 that is identified by the eye gene discovery tool iSyTE.

RESULTS

In situ hybridization, Western blotting and immunostaining confirmed highly enriched expression of Caprin2 mRNA and protein in mouse embryonic and postnatal lens. To gain insight into its function, lens-specific Caprin2 conditional knockout (cKO) mouse mutants were generated using a lens-Cre deleter line Pax6GFPCre. Phenotypic analysis of Caprin2(cKO/cKO) mutants revealed distinct eye defects at variable penetrance. Wheat germ agglutinin staining and scanning electron microscopy demonstrated that Caprin2(cKO/cKO) mutants have an abnormally compact lens nucleus, which is the core of the lens comprised of centrally located terminally differentiated fiber cells. Additionally, Caprin2(cKO/cKO) mutants also exhibited at 8% penetrance a developmental defect that resembles a human condition called Peters anomaly, wherein the lens and the cornea remain attached by a persistent stalk.

CONCLUSIONS

These data suggest that a conserved RBP Caprin2 functions in distinct morphological events in mammalian eye development.