Loss of Sip1 leads to migration defects and retention of ectodermal markers during lens development

Transcriptome Meta-Analysis Uncovers Cell-Specific Regulatory Relationships in Embryonic, Juvenile, Adult, and Aged Mouse Lens Epithelium and Fibers

Purpose

The lens transcriptome has been examined using microarrays and RNA-sequencing (RNA-seq). These omics data are the basis of the bioinformatics web-resource iSyTE that has identified new genes involved in lens development and cataract. The lens predominantly contains epithelial and fiber cells, and yet, presently, iSyTE is based on whole lens data. To gain cell-specific regulatory insights, we meta-analyzed isolated epithelium and fiber transcriptomes from embryonic/postnatal, adult and aged lenses.

Methods

Mouse lens epithelium and fiber transcriptome public datasets at embryonic (E) and postnatal (P) stages E12.5, E14.5, E16.5, E18.5, P0.5, P0, P5, P13, and age one month, three months, six months, and two years were analyzed. Microarray or RNA-seq data were analyzed by appropriate methods and compared to other resources (e.g., Cat-Map, CompBio).

Results

Across all RNA-seq datasets examined, 2466 genes are differentially expressed between epithelium and fibers, of which 106 are cataract-linked. Gene ontology enrichment validates epithelial and fiber expression, corroborating the meta-analysis. Whole embryonic-body-in silico subtraction and other analyses identify several new high-priority epithelial- and/or fiber-enriched genes (e.g., Casz1, Ell2). Furthermore, new insights into cell-specific regulatory processes at distinct stages are identified (e.g., ribonucleoprotein regulation in E12.5 epithelium). Finally, this data is made accessible at iSyTE (https://research.bioinformatics.udel.edu/iSyTE/).

Conclusions

This spatiotemporal transcriptome meta-analysis comprehensively informs on epithelium- and fiber-specific regulatory processes in developing, adult and aged lenses. Notably, it includes the first description of an embryonic stage (i.e., E12.5) representing early primary fiber differentiation, thus informing on the initial transcriptome changes as lens cell-types are readily distinguishable.

Unilateral progressive anterior iris adhesions in Mowat-Wilson syndrome: a new ocular finding

Ocular associations in Mowat-Wilson syndrome (MWS) are rare. Those involving the anterior segment are scarce in the literature. We describe a child with genetic confirmation of MWS that presented with acquired onset of unilateral anterior iris adhesions with no known trauma.

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.

ATF4 May Be Essential for Adaption of the Ocular Lens to Its Avascular Environment

The late embryonic mouse lens requires the transcription factor ATF4 for its survival although the underlying mechanisms were unknown. Here, RNAseq analysis revealed that E16.5 Atf4 null mouse lenses downregulate the mRNA levels of lens epithelial markers as well as known markers of late lens fiber cell differentiation. However, a comparison of this list of differentially expressed genes (DEGs) with other known transcriptional regulators of lens development indicated that ATF4 expression is not directly controlled by the previously described lens gene regulatory network. Pathway analysis revealed that the Atf4 DEG list was enriched in numerous genes involved in nutrient transport, amino acid biosynthesis, and tRNA charging. These changes in gene expression likely result in the observed reductions in lens free amino acid and glutathione levels, which would result in the observed low levels of extractable lens protein, finally leading to perinatal lens disintegration. These data demonstrate that ATF4, via its function in the integrated stress response, is likely to play a crucial role in mediating the adaption of the lens to the avascularity needed to maintain lens transparency.

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.

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.

First Case Report of Developmental Bilateral Cataract with a Novel Mutation in the ZEB2 Gene Observed in Mowat-Wilson Syndrome

Background: Mowat-Wilson syndrome (MWS) is extremely rare multisystemic autosomal dominant disorder caused by mutations in the Zinc Finger E-Box Binding Homeobox 2 (ZEB2) gene. Ocular pathologies are one of the symptoms that appear in the clinical picture of MWS individuals, but not many have been described so far. Pathologies such as optic nerve or retinal epithelium atrophy, iris or optic disc coloboma as well as congenital cataracts have been most frequently described until now. Therefore, we would like to report the first case of bilateral developmental cataract in a 9-year-old girl with MWS who underwent successful cataract surgery with intraocular lens implantation. Case Presentation: A 9-year-old girl, diagnosed with p.Gln694Ter mutation in ZEB2 gene and suspicion of MWS was referred to the Children's Outpatient Ophthalmology Clinic for ophthalmological evaluation. Her previous assessments revealed abnormalities of the optic nerve discs. The patient was diagnosed with atrophy of the optic nerves, convergent strabismus, and with-the-rule astigmatism. One year later, during the follow-up visit, the patient was presented with decreased visual acuity (VA), developmental total cataract in the right eye and a developmental partial cataract in the left eye. This resulted in decreased VA confirmed by deteriorated responses in visual evoked potential (VEP) test. The girl underwent a two-stage procedure of cataract removal, first of one eye and then of the other eye with artificial lens implants. In the 2 years following the operation, no complications were observed and, most remarkably, VA improved significantly. Conclusions: The ZEB2 gene is primarily responsible for encoding the Smad interaction protein 1 (SIP1), which is involved in the proper development of various eye components. When mutated, it results in multilevel abnormalities, also in the proper lens formation, that prevent the child from normal vision development. This typically results in the formation of congenital cataracts in children with MWS syndrome, however, our case shows that it also leads to the formation of developmental cataracts. This is presumably due to the effect of the lack of SIP1 on other genes, altering their downstream expression and is a novel insight into the importance of the SIP1 in the occurrence of ocular pathologies. To the best of our knowledge, this is the first case of bilateral developmental cataract in the context of MWS. Moreover, a novel mutation (p.Gln694Ter) in the ZEB2 gene was found corresponding to this syndrome entity. This report allows us to gain a more comprehensive insight into the genetic spectrum and the corresponding phenotypic features in MWS syndrome patients.

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.

ZEB2, the Mowat-Wilson Syndrome Transcription Factor: Confirmations, Novel Functions, and Continuing Surprises

After its publication in 1999 as a DNA-binding and SMAD-binding transcription factor (TF) that co-determines cell fate in amphibian embryos, ZEB2 was from 2003 studied by embryologists mainly by documenting the consequences of conditional, cell-type specific Zeb2 knockout (cKO) in mice. In between, it was further identified as causal gene causing Mowat-Wilson Syndrome (MOWS) and novel regulator of epithelial-mesenchymal transition (EMT). ZEB2's functions and action mechanisms in mouse embryos were first addressed in its main sites of expression, with focus on those that helped to explain neurodevelopmental and neural crest defects seen in MOWS patients. By doing so, ZEB2 was identified in the forebrain as the first TF that determined timing of neuro-/gliogenesis, and thereby also the extent of different layers of the cortex, in a cell non-autonomous fashion, i.e., by its cell-intrinsic control within neurons of neuron-to-progenitor paracrine signaling. Transcriptomics-based phenotyping of Zeb2 mutant mouse cells have identified large sets of intact-ZEB2 dependent genes, and the cKO approaches also moved to post-natal brain development and diverse other systems in adult mice, including hematopoiesis and various cell types of the immune system. These new studies start to highlight the important adult roles of ZEB2 in cell-cell communication, including after challenge, e.g., in the infarcted heart and fibrotic liver. Such studies may further evolve towards those documenting the roles of ZEB2 in cell-based repair of injured tissue and organs, downstream of actions of diverse growth factors, which recapitulate developmental signaling principles in the injured sites. Evident questions are about ZEB2's direct target genes, its various partners, and ZEB2 as a candidate modifier gene, e.g., in other (neuro)developmental disorders, but also the accurate transcriptional and epigenetic regulation of its mRNA expression sites and levels. Other questions start to address ZEB2's function as a niche-controlling regulatory TF of also other cell types, in part by its modulation of growth factor responses (e.g., TGFβ/BMP, Wnt, Notch). Furthermore, growing numbers of mapped missense as well as protein non-coding mutations in MOWS patients are becoming available and inspire the design of new animal model and pluripotent stem cell-based systems. This review attempts to summarize in detail, albeit without discussing ZEB2's role in cancer, hematopoiesis, and its emerging roles in the immune system, how intense ZEB2 research has arrived at this exciting intersection.

Fibronectin has multifunctional roles in posterior capsular opacification (PCO)

Fibrotic posterior capsular opacification (PCO), one of the major complications of cataract surgery, occurs when lens epithelial cells (LCs) left behind post cataract surgery (PCS) undergo epithelial to mesenchymal transition, migrate into the optical axis and produce opaque scar tissue. LCs left behind PCS robustly produce fibronectin, although its roles in fibrotic PCO are not known. In order to determine the function of fibronectin in PCO pathogenesis, we created mice lacking the fibronectin gene (FN conditional knock out -FNcKO) from the lens. While animals from this line have normal lenses, upon lens fiber cell removal which models cataract surgery, FNcKO LCs exhibit a greatly attenuated fibrotic response from 3 days PCS onward as assessed by a reduction in surgery-induced cell proliferation, and fibrotic extracellular matrix (ECM) production and deposition. This is correlated with less upregulation of Transforming Growth Factor β (TGFβ) and integrin signaling in FNcKO LCs PCS concomitant with sustained Bone Morphogenetic Protein (BMP) signaling and elevation of the epithelial cell marker E cadherin. Although the initial fibrotic response of FNcKO LCs was qualitatively normal at 48 h PCS as measured by the upregulation of fibrotic marker protein αSMA, RNA sequencing revealed that the fibrotic response was already quantitatively attenuated at this time, as measured by the upregulation of mRNAs encoding molecules that control, and are controlled by, TGFβ signaling, including many known markers of fibrosis. Most notably, gremlin-1, a known regulator of TGFβ superfamily signaling, was upregulated sharply in WT LCs PCS, while this response was attenuated in FNcKO LCs. As exogenous administration of either active TGFβ1 or gremlin-1 to FNcKO lens capsular bags rescued the attenuated fibrotic response of fibronectin null LCs PCS including the loss of SMAD2/3 phosphorylation, this suggests that fibronectin plays multifunctional roles in fibrotic PCO development.

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.

Cerebral ischemia induces TRPC6 via HIF1α/ZEB2 axis in the glomerular podocytes and contributes to proteinuria

Podocytes are specialized cells of the glomerulus and key component of the glomerular filtration apparatus (GFA). GFA regulates the permselectivity and ultrafiltration of blood. The mechanism by which the integrity of the GFA is compromised and manifest in proteinuria during ischemic stroke remains enigmatic. We investigated the mechanism of ischemic hypoxia-induced proteinuria in a middle cerebral artery occlusion (MCAO) model. Ischemic hypoxia resulted in the accumulation of HIF1α in the podocytes that resulted in the increased expression of ZEB2 (Zinc finger E-box-binding homeobox 2). ZEB2, in turn, induced TRPC6 (transient receptor potential cation channel, subfamily C, member 6), which has increased selectivity for calcium. Elevated expression of TRPC6 elicited increased calcium influx and aberrant activation of focal adhesion kinase (FAK) in podocytes. FAK activation resulted in the stress fibers reorganization and podocyte foot process effacement. Our study suggests overactive HIF1α/ZEB2 axis during ischemic-hypoxia raises intracellular calcium levels via TRPC6 and consequently altered podocyte structure and function thus contributes to proteinuria.

The crucial role of ZEB2: From development to epithelial-to-mesenchymal transition and cancer complexity

Zinc finger E-box binding homeobox 2 (ZEB2) is a DNA-binding transcription factor, which is mainly involved in epithelial-to-mesenchymal transition (EMT). EMT is a conserved process during which mature and adherent epithelial-like state is converted into a mobile mesenchymal state. Emerging data indicate that ZEB2 plays a pivotal role in EMT-induced processes such as development, differentiation, and malignant mechanisms, for example, drug resistance, cancer stem cell-like traits, apoptosis, survival, cell cycle arrest, tumor recurrence, and metastasis. In this regard, the understanding of mentioned subjects in the development of normal and cancerous cells could be helpful in cancer complexity of diagnosis and therapy. In this study, we review recent findings about the biological properties of ZEB2 in healthy and cancerous states to find new approaches for cancer treatment.

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.

A screen for deeply conserved non-coding GWAS SNPs uncovers a MIR-9-2 functional mutation associated to retinal vasculature defects in human

Thousands of human disease-associated single nucleotide polymorphisms (SNPs) lie in the non-coding genome, but only a handful have been demonstrated to affect gene expression and human biology. We computationally identified risk-associated SNPs in deeply conserved non-exonic elements (CNEs) potentially contributing to 45 human diseases. We further demonstrated that human CNE1/rs17421627 associated with retinal vasculature defects showed transcriptional activity in the zebrafish retina, while introducing the risk-associated allele completely abolished CNE1 enhancer activity. Furthermore, deletion of CNE1 led to retinal vasculature defects and to a specific downregulation of microRNA-9, rather than MEF2C as predicted by the original genome-wide association studies. Consistent with these results, miR-9 depletion affects retinal vasculature formation, demonstrating MIR-9-2 as a critical gene underpinning the associated trait. Importantly, we validated that other CNEs act as transcriptional enhancers that can be disrupted by conserved non-coding SNPs. This study uncovers disease-associated non-coding mutations that are deeply conserved, providing a path for in vivo testing to reveal their cis-regulated genes and biological roles.

Lens Epithelial Cells Initiate an Inflammatory Response Following Cataract Surgery

Purpose

Lens epithelial cell (LEC) conversion to myofibroblast is responsible for fibrotic cataract surgery complications including posterior capsular opacification. While transforming growth factor beta (TGFβ) signaling is important, the mechanisms by which the TGFβ pathway is activated post cataract surgery (PCS) are not well understood.

Methods

RNA-seq was performed on LECs obtained from a mouse cataract surgery model at the time of surgery and 24 hours later. Bioinformatic analysis was performed with iPathwayGuide. Expression dynamics were determined by immunofluorescence.

Results

The LEC transcriptome is massively altered by 24 hours PCS. The differentially expressed genes included those important for lens biology, and fibrotic markers. However, the most dramatic changes were in the expression of genes regulating the innate immune response, with the top three altered genes exhibiting greater than 1000-fold upregulation. Immunolocalization revealed that CXCL1, S100a9, CSF3, COX-2, CCL2, LCN2, and HMOX1 protein levels upregulate in LECs between 1 hour and 6 hours PCS and peak at 24 hours PCS, while their levels sharply attenuate by 3 days PCS. This massive upregulation of known inflammatory mediators precedes the infiltration of neutrophils into the eye at 18 hours PCS, the upregulation of canonical TGFβ signaling at 48 hours PCS, and the infiltration of macrophages at 3 days PCS.

Conclusions

These data demonstrate that LECs produce proinflammatory cytokines immediately following lens injury that could drive postsurgical flare, and suggest that inflammation may be a major player in the onset of lens-associated fibrotic disease PCS.

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.

Spatiotemporal dynamics of canonical Wnt signaling during embryonic eye development and posterior capsular opacification (PCO)

The appropriate spatial and temporal regulation of canonical Wnt signaling is vital for eye development. However, the literature often conflicts on the distribution of canonical Wnt signaling in the eye. Here, using a sensitive mouse transgenic reporter line, we report a detailed re-evaluation of the spatiotemporal dynamics of canonical Wnt signaling in the developing eye. Canonical Wnt activity was dynamic in the optic vesicle and later in the retina, while it was absent from the ectodermal precursors of the lens and corneal epithelium. However, later in corneal development, canonical Wnt reporter activity was detected in corneal stroma and endothelium precursors as they form from the neural crest, although this was lost around birth. Interestingly, while no canonical Wnt signaling was detected in the corneal limbus or basal cells at any developmental stage, it was robust in adult corneal wing and squamous epithelial cells. While canonical Wnt reporter activity was also absent from the postnatal lens, upon lens injury intended to model cataract surgery, it upregulated within 12 h in remnant lens epithelial cells, and co-localized with alpha smooth muscle actin in fibrotic lens epithelial cells from 48 h post-surgery onward. This pattern correlated with downregulation of the inhibitor of canonical Wnt signaling, Dkk3. These data demonstrate that canonical Wnt signaling is dynamic within the developing eye and upregulates in lens epithelial cells in response to lens injury. As canonical Wnt signaling can collaborate with TGFβ to drive fibrosis in other systems, these data offer the first evidence in a lens-injury model that canonical Wnt may synergize with TGFβ signaling to drive fibrotic posterior capsular opacification (PCO).

Transcriptome analysis of developing lens reveals abundance of novel transcripts and extensive splicing alterations

Lens development involves a complex and highly orchestrated regulatory program. Here, we investigate the transcriptomic alterations and splicing events during mouse lens formation using RNA-seq data from multiple developmental stages, and construct a molecular portrait of known and novel transcripts. We show that the extent of novelty of expressed transcripts decreases significantly in post-natal lens compared to embryonic stages. Characterization of novel transcripts into partially novel transcripts (PNTs) and completely novel transcripts (CNTs) (novelty score ≥ 70%) revealed that the PNTs are both highly conserved across vertebrates and highly expressed across multiple stages. Functional analysis of PNTs revealed their widespread role in lens developmental processes while hundreds of CNTs were found to be widely expressed and predicted to encode for proteins. We verified the expression of four CNTs across stages. Examination of splice isoforms revealed skipped exon and retained intron to be the most abundant alternative splicing events during lens development. We validated by RT-PCR and Sanger sequencing, the predicted splice isoforms of several genes Banf1, Cdk4, Cryaa, Eif4g2, Pax6, and Rbm5. Finally, we present a splicing browser Eye Splicer (http://www.iupui.edu/~sysbio/eye-splicer/), to facilitate exploration of developmentally altered splicing events and to improve understanding of post-transcriptional regulatory networks during mouse lens development.

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.

Association and interaction of myopia with SNP markers rs13382811 and rs6469937 at ZFHX1B and SNTB1 in Han Chinese and European populations

PURPOSE

Previously, a genome-wide association study (GWAS) identified rs13382811 (near ZFHX1B) and rs6469937 (near SNTB1) to be associated with high myopia. The present study evaluates the association of these two single nucleotide polymorphisms (SNPs) with moderate to high myopia in two Chinese cohorts and two cohorts of European populations.

METHODS

Two Chinese university student cohorts, including one with 300 unrelated subjects with high myopia and 308 emmetropic controls from Guangzhou and a second with 96 unrelated individuals with moderate to high myopia and 96 emmetropic controls of Chaoshanese origin in Guangzhou, were enrolled in this study. Two SNPs, rs6469937 and rs13382811, were selected for genotyping based on their reported associations with severe myopia. The SNPs were genotyped via DNA sequencing. In addition, association analysis of both SNPs was performed using genotype data from the database of Genotypes and Phenotypes (dbGaP) involving a total of 2,423 samples in two independent cohorts of European-derived populations, as follows: Kooperative Gesundheitsforschung in der Region Augsburg (KORA) and TwinsUK. The allelic and genotypic distribution among cases and controls were analyzed using the Chi-square test. Logistic regression was used to evaluate the SNP-SNP interaction. Fisher's exact test was used for two-SNP comparisons.

RESULTS

In the Guangzhou cohort, SNP rs13382811 near ZFHX1B showed significant association with high myopia (p_allelic = 0.0001, p_genotypic = 4.07 × 10^-5), with the minor T allele showing an increased risk of high myopia (odds ratio [OR] = 1.68, 95% confidence interval [CI] = 1.28-2.20). SNP rs6469937 near SNTB1 showed nominal evidence of association (p_allelic = 0.0085, p_genotypic = 0.0166), which did not withstand correction for multiple testing. No significant association was detected in the smaller Chaoshan cohort alone. The association of SNPs rs13382811 and rs6469937 remained significant when both Han Chinese cohorts were combined (p_allelic = 0.0033 and 0.0016, respectively), and it was also significant under the genotypic test (p_genotypic = 0.0036 and 0.0053, respectively). When both SNPs were considered together under a recessive model, their significance increased (p = 8.37 × 10^-4), as did their effect (OR = 4.09, 95%CI = 1.7-9.8). The association between either of these two SNPs alone and myopia did not replicate significantly in the combined cohorts of European descent, providing only suggestive results (p_allelic = 0.0088 for rs13382811 and p_allelic = 0.0319 for rs6469937). However, the effects of the combined SNPs showed significant association (p = 8.2 × 10^-4; OR = 1.56, 95%CI = 1.2-2.0). While the risk for myopia increased with risk alleles from both SNPs, the increase was additive rather representing a multiplicative interaction in both populations.

CONCLUSIONS

Our study confirms that the two susceptibility loci ZFHX1B and SNTB1 are associated with moderate to high myopia in a Han Chinese population, as well as in a European population, when both SNPs are combined. These results confirm previous reports of their associations, extend these observations to a European population, and suggest that additional interactive and possibly population-specific genetic or environmental factors may affect their contribution to myopia.

β1-Integrin Deletion From the Lens Activates Cellular Stress Responses Leading to Apoptosis and Fibrosis

Purpose

Previous research showed that the absence of β1-integrin from the mouse lens after embryonic day (E) 13.5 (β1MLR10) leads to the perinatal apoptosis of lens epithelial cells (LECs) resulting in severe microphthalmia. This study focuses on elucidating the molecular connections between β1-integrin deletion and this phenotype.

Methods

RNA sequencing was performed to identify differentially regulated genes (DRGs) in β1MLR10 lenses at E15.5. By using bioinformatics analysis and literature searching, Egr1 (early growth response 1) was selected for further study. The activation status of certain signaling pathways (focal adhesion kinase [FAK]/Erk, TGF-β, and Akt signaling) was studied via Western blot and immunohistochemistry. Mice lacking both β1-integrin and Egr1 genes from the lenses were created (β1MLR10/Egr1-/-) to study their relationship.

Results

RNA sequencing identified 120 DRGs that include candidates involved in the cellular stress response, fibrosis, and/or apoptosis. Egr1 was investigated in detail, as it mediates cellular stress responses in various cell types, and is recognized as an upstream regulator of numerous other β1MLR10 lens DRGs. In β1MLR10 mice, Egr1 levels are elevated shortly after β1-integrin loss from the lens. Further, pErk1/2 and pAkt are elevated in β1MLR10 LECs, thus providing the potential signaling mechanism that causes Egr1 upregulation in the mutant. Indeed, deletion of Egr1 from β1MLR10 lenses partially rescues the microphthalmia phenotype.

Conclusions

β1-integrin regulates the appropriate levels of Erk1/2 and Akt phosphorylation in LECs, whereas its deficiency results in the overexpression of Egr1, culminating in reduced cell survival. These findings provide insight into the molecular mechanism underlying the microphthalmia observed in β1MLR10 mice.

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.

β1-integrin controls cell fate specification in early lens development

Integrins are heterodimeric cell surface molecules that mediate cell-extracellular matrix (ECM) adhesion, ECM assembly, and regulation of both ECM and growth factor induced signaling. However, the developmental context of these diverse functions is not clear. Loss of β1-integrin from the lens vesicle (mouse E10.5) results in abnormal exit of anterior lens epithelial cells (LECs) from the cell cycle and their aberrant elongation toward the presumptive cornea by E12.5. These cells lose expression of LEC markers and initiate expression of the Maf (also known as c-Maf) and Prox1 transcription factors as well as other lens fiber cell markers. β1-integrin null LECs also upregulate the ERK, AKT and Smad1/5/8 phosphorylation indicative of BMP and FGF signaling. By E14.5, β1-integrin null lenses have undergone a complete conversion of all lens epithelial cells into fiber cells. These data suggest that shortly after lens vesicle closure, β1-integrin blocks inappropriate differentiation of the lens epithelium into fibers, potentially by inhibiting BMP and/or FGF receptor activation. Thus, β1-integrin has an important role in fine-tuning the response of the early lens to the gradient of growth factors that regulate lens fiber cell differentiation.

Sip1 regulates the generation of the inner nuclear layer retinal cell lineages in mammals

The transcription factor Sip1 (Zeb2) plays multiple roles during CNS development from early acquisition of neural fate to cortical neurogenesis and gliogenesis. In humans, SIP1 (ZEB2) haploinsufficiency leads to Mowat-Wilson syndrome, a complex congenital anomaly including intellectual disability, epilepsy and Hirschsprung disease. Here we uncover the role of Sip1 in retinogenesis. Somatic deletion of Sip1 from mouse retinal progenitors primarily affects the generation of inner nuclear layer cell types, resulting in complete loss of horizontal cells and reduced numbers of amacrine and bipolar cells, while the number of Muller glia is increased. Molecular analysis places Sip1 downstream of the eye field transcription factor Pax6 and upstream of Ptf1a in the gene network required for generating the horizontal and amacrine lineages. Intriguingly, characterization of differentiation dynamics reveals that Sip1 has a role in promoting the timely differentiation of retinal interneurons, assuring generation of the proper number of the diverse neuronal and glial cell subtypes that constitute the functional retina in mammals.

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.

The molecular mechanisms underlying lens fiber elongation

Lens fiber cells are highly elongated cells with complex membrane morphologies that are critical for the transparency of the ocular lens. Investigations into the molecular mechanisms underlying lens fiber cell elongation were first reported in the 1960s, however, our understanding of the process is still poor nearly 50 years later. This review summarizes what is currently hypothesized about the regulation of lens fiber cell elongation along with the available experimental evidence, and how this information relates to what is known about the regulation of cell shape/elongation in other cell types, particularly neurons.

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.

Prox1 and fibroblast growth factor receptors form a novel regulatory loop controlling lens fiber differentiation and gene expression

Lens epithelial cells differentiate into lens fibers (LFs) in response to a fibroblast growth factor (FGF) gradient. This cell fate decision requires the transcription factor Prox1, which has been hypothesized to promote cell cycle exit in differentiating LF cells. However, we find that conditional deletion of Prox1 from mouse lenses results in a failure in LF differentiation despite maintenance of normal cell cycle exit. Instead, RNA-seq demonstrated that Prox1 functions as a global regulator of LF cell gene expression. Intriguingly, Prox1 also controls the expression of fibroblast growth factor receptors (FGFRs) and can bind to their promoters, correlating with decreased downstream signaling through MAPK and AKT in Prox1 mutant lenses. Further, culturing rat lens explants in FGF increased their expression of Prox1, and this was attenuated by the addition of inhibitors of MAPK. Together, these results describe a novel feedback loop required for lens differentiation and morphogenesis, whereby Prox1 and FGFR signaling interact to mediate LF differentiation in response to FGF.

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.

Rap1 GTPase is required for mouse lens epithelial maintenance and morphogenesis

Rap1, a Ras-like small GTPase, plays a crucial role in cell-matrix adhesive interactions, cell-cell junction formation, cell polarity and migration. The role of Rap1 in vertebrate organ development and tissue architecture, however, remains elusive. We addressed this question in a mouse lens model system using a conditional gene targeting approach. While individual germline deficiency of either Rap1a or Rap1b did not cause overt defects in mouse lens, conditional double deficiency (Rap1 cKO) prior to lens placode formation led to an ocular phenotype including microphthalmia and lens opacification in embryonic mice. The embryonic Rap1 cKO mouse lens exhibited striking defects including loss of E-cadherin- and ZO-1-based cell-cell junctions, disruption of paxillin and β1-integrin-based cell adhesive interactions along with abnormalities in cell shape and apical-basal polarity of epithelium. These epithelial changes were accompanied by increased levels of α-smooth muscle actin, vimentin and N-cadherin, and expression of transcriptional suppressors of E-cadherin (Snai1, Slug and Zeb2), and a mesenchymal metabolic protein (Dihydropyrimidine dehydrogenase). Additionally, while lens differentiation was not overtly affected, increased apoptosis and dysregulated cell cycle progression were noted in epithelium and fibers in Rap1 cKO mice. Collectively these observations uncover a requirement for Rap1 in maintenance of lens epithelial phenotype and morphogenesis.

Molecular characterization of mouse lens epithelial cell lines and their suitability to study RNA granules and cataract associated genes

The discovery of cytosolic RNA granule (RG) component proteins associated with human cataract has initiated investigations on post-transcriptional mechanisms of gene expression control in the lens. Application of established mouse lens epithelial cell lines (LECs) can provide rapid insights on RG function in lens cells, especially because mouse mutants in several RG components are not available. However, although these LECs represent potential reagents for such analyses, they are uncharacterized for lens gene expression or RG formation. Therefore, a detailed molecular and cellular characterization of three permanent mouse LECs 17EM15, 21EM15 and αTN4 is performed in this study. Comparative analysis between microarray gene expression datasets on LEC 21EM15 and iSyTE lens tissue demonstrates that 30% of top 200 iSyTE identified lens-enriched genes are expressed in these cells. Majority of these candidates are independently validated to either have lens expression, function or linkage to cataract. Moreover, analysis of microarray data with genes described in Cat-Map, an online database of cataract associated genes and loci, demonstrates that 131 genes linked to cataract loci are expressed in 21EM15 cells. Furthermore, gene expression in LECs is compared to isolated lens epithelium or fiber cells by qRT-PCR and by comparative analyses with publically available epithelium or fiber-specific microarray and RNA-seq (sequencing) datasets. Expression of select candidate genes was validated by regular and real-time quantitative RT-PCR. Expression of lens epithelium-enriched genes Foxe3, Pax6, Anxa4 and Mcm4 is up-regulated in LEC lines, compared to isolated lens fiber cells. Moreover, similar to isolated lens epithelium, all three LECs exhibit down-regulation of fiber cell-expressed genes Crybb1, Mip and Prox1 when compared to fiber cells. These data indicate that the LEC lines exhibit greater similarity to lens epithelium than to fiber cells. Compared to non-lens cell line NIH3T3, LECs exhibit significantly enriched expression of transcription factors with important function in the lens, namely Pax6, Foxe3 and Prox1. In addition to these genes, all three LECs also express key lens- and cataract-associated genes, namely Dkk3, Epha2, Hsf4, Jag1, Mab21l1, Meis1, Pknox1, Pou2f1, Sfrp1, Sparc, Tdrd7 and Trpm3. Additionally, 21EM15 microarrays indicate expression of Chmp4b, Cryab and Tcfap2a among others important genes. Immunostaining with makers for Processing bodies (P-bodies) and Stress granules (SGs) demonstrates that these classes of RGs are robustly expressed in all three LECs. Moreover, under conditions of stress, 17EM15 and αTN4 exhibit significantly higher numbers of P-bodies and SGs compared to NIH3T3 cells. In sum, these data indicate that mouse LECs 21EM15, 17EM15 and αTN4 express key lens or cataract genes, are similar to lens epithelium than fiber cells, and exhibit high levels of P-bodies and SGs, indicating their suitability for investigating gene expression control and RG function in lens-derived cells.

Development of novel filtering criteria to analyze RNA-sequencing data obtained from the murine ocular lens during embryogenesis

Next-generation sequencing of the transcriptome (RNA-Seq) is a powerful method that allows for the quantitative determination of absolute gene expression, and can be used to investigate how these levels change in response to an experimental manipulation or disease condition. The sensitivity of this method allows one to analyze transcript levels of all expressed genes, including low abundance transcripts that encode important regulatory molecules, providing valuable insights into the global effects of experimental manipulations. However, this increased sensitivity can also make it challenging to ascertain which expression changes are biologically significant. Here, we describe a novel set of filtering criteria - based on biological insights and computational approaches - that were applied to prioritize genes for further study from an extensive number of differentially expressed transcripts in lenses lacking Smad interacting protein 1 (Sip1) obtained via RNA-Seq by Manthey and colleagues in Mechanisms of Development (Manthey et al., 2014). Notably, this workflow allowed an original list of over 7,100 statistically significant differentially expressed genes (DEGs) to be winnowed down to 190 DEGs that likely play a biologically significant role in Sip1 function during lens development. Focusing on genes whose expression was upregulated or downregulated in a manner opposite to what normally occurs during lens development, we identified 78 genes that appear to be strongly dependent on Sip1 function. From these data (GEO accession number GSE49949), it appears that Sip1 regulates multiple genes in the lens that are generally distinct from those regulated by Sip1 in other cellular contexts, including genes whose expression is prominent in the early head ectoderm, from which the lens differentiates. Further, the analysis criteria outlined here represent a filtering scheme that can be used to prioritize genes in future RNA-Seq investigations performed at this stage of ocular lens development.

The cellular and molecular mechanisms of vertebrate lens development

The ocular lens is a model system for understanding important aspects of embryonic development, such as cell specification and the spatiotemporally controlled formation of a three-dimensional structure. The lens, which is characterized by transparency, refraction and elasticity, is composed of a bulk mass of fiber cells attached to a sheet of lens epithelium. Although lens induction has been studied for over 100 years, recent findings have revealed a myriad of extracellular signaling pathways and gene regulatory networks, integrated and executed by the transcription factor Pax6, that are required for lens formation in vertebrates. This Review summarizes recent progress in the field, emphasizing the interplay between the diverse regulatory mechanisms employed to form lens progenitor and precursor cells and highlighting novel opportunities to fill gaps in our understanding of lens tissue morphogenesis.

Comparative transcriptome analysis of epithelial and fiber cells in newborn mouse lenses with RNA sequencing

PURPOSE

The ocular lens contains only two cell types: epithelial cells and fiber cells. The epithelial cells lining the anterior hemisphere have the capacity to continuously proliferate and differentiate into lens fiber cells that make up the large proportion of the lens mass. To understand the transcriptional changes that take place during the differentiation process, high-throughput RNA-Seq of newborn mouse lens epithelial cells and lens fiber cells was conducted to comprehensively compare the transcriptomes of these two cell types.

METHODS

RNA from three biologic replicate samples of epithelial and fiber cells from newborn FVB/N mouse lenses was isolated and sequenced to yield more than 24 million reads per sample. Sequence reads that passed quality filtering were mapped to the reference genome using Genomic Short-read Nucleotide Alignment Program (GSNAP). Transcript abundance and differential gene expression were estimated using the Cufflinks and DESeq packages, respectively. Gene Ontology enrichment was analyzed using GOseq. RNA-Seq results were compared with previously published microarray data. The differential expression of several biologically important genes was confirmed using reverse transcription (RT)-quantitative PCR (qPCR).

RESULTS

Here, we present the first application of RNA-Seq to understand the transcriptional changes underlying the differentiation of epithelial cells into fiber cells in the newborn mouse lens. In total, 6,022 protein-coding genes exhibited differential expression between lens epithelial cells and lens fiber cells. To our knowledge, this is the first study identifying the expression of 254 long intergenic non-coding RNAs (lincRNAs) in the lens, of which 86 lincRNAs displayed differential expression between the two cell types. We found that RNA-Seq identified more differentially expressed genes and correlated with RT-qPCR quantification better than previously published microarray data. Gene Ontology analysis showed that genes upregulated in the epithelial cells were enriched for extracellular matrix production, cell division, migration, protein kinase activity, growth factor binding, and calcium ion binding. Genes upregulated in the fiber cells were enriched for proteosome complexes, unfolded protein responses, phosphatase activity, and ubiquitin binding. Differentially expressed genes involved in several important signaling pathways, lens structural components, organelle loss, and denucleation were also highlighted to provide insights into lens development and lens fiber differentiation.

CONCLUSIONS

RNA-Seq analysis provided a comprehensive view of the relative abundance and differential expression of protein-coding and non-coding transcripts from lens epithelial cells and lens fiber cells. This information provides a valuable resource for studying lens development, nuclear degradation, and organelle loss during fiber differentiation, and associated diseases.

The Zeb proteins δEF1 and Sip1 may have distinct functions in lens cells following cataract surgery

PURPOSE

Posterior capsular opacification (PCO), the most prevalent side effect of cataract surgery, occurs when residual lens epithelial cells (LECs) undergo fiber cell differentiation or epithelial-to-mesenchymal transition (EMT). Here, we used a murine cataract surgery model to investigate the role of the Zeb proteins, Smad interacting protein 1 (Sip1) and δ-crystallin enhancer-binding factor 1 (δEF1), during PCO.

METHODS

Extracapsular extraction of lens fiber cells was performed on wild-type and Sip1 knockout mice. Protein expression patterns were assessed at multiple time points after surgery using confocal immunofluorescence. βB1-Crystallin mRNA levels were measured using quantitative RT-PCR. We used Transfac searches to identify δEF1 binding sites in the βB1-crystallin promoter and transfection analysis to test the ability of δEF1 to regulate βB1-crystallin expression.

RESULTS

δEF1, which, in other systems, can activate fibrotic genes (e.g., α-smooth muscle actin) and repress epithelial genes, upregulates by 48 hours after fiber cell removal. In culture, δEF1 repressed βB1-crystallin promoter activity, suggesting that it may also turn off lens gene expression following surgery, contributing to "fibrotic PCO" development. Sip1 also upregulates in LECs by 48 hours, but analysis of Sip1 knockout lenses demonstrated that Sip1 does not play a major role in EMT or fiber cell differentiation after surgery. However, Sip1 knockout LECs do express the ectodermal marker keratin 8, suggesting that Sip1 may limit the reprogramming of residual LECs to an embryonic state.

CONCLUSIONS

Zeb transcription factors likely play important, but distinct roles in PCO development after cataract surgery.

Beta-1 integrin is important for the structural maintenance and homeostasis of differentiating fiber cells

β1-Integrin is a heterodimeric transmembrane protein that has roles in both cell-extra-cellular matrix and cell-cell interactions. Conditional deletion of β1-integrin from all lens cells during embryonic development results in profound lens defects, however, it is less clear whether this reflects functions in the lens epithelium alone or whether this protein plays a role in lens fibers. Thus, a conditional approach was used to delete β1-integrin solely from the lens fiber cells. This deletion resulted in two distinct phenotypes with some lenses exhibiting cataracts while others were clear, albeit with refractive defects. Analysis of "clear" conditional knockout lenses revealed that they had profound defects in fiber cell morphology associated with the loss of the F-actin network. Physiological measurements found that the lens fiber cells had a twofold increase in gap junctional coupling, perhaps due to differential localization of connexins 46 and 50, as well as increased water permeability. This would presumably facilitate transport of ions and nutrients through the lens, and may partially explain how lenses with profound structural abnormalities can maintain transparency. In summary, β1-integrin plays a role in maintaining the cellular morphology and homeostasis of the lens fiber cells.

Lens extrusion from Laminin alpha 1 mutant zebrafish

We report analysis of the ocular lens phenotype of the recessive, larval lethal zebrafish mutant, lama1^a69/a69. Previous work revealed that this mutant has a shortened body axis and eye defects including a defective hyaloid vasculature, focal corneal dysplasia, and loss of the crystalline lens. While these studies highlight the importance of laminin α1 in lens development, a detailed analysis of the lens defects seen in these mutants was not reported. In the present study, we analyze the lenticular anomalies seen in the lama1^a69/a69 mutants and show that the lens defects result from the anterior extrusion of lens material from the eye secondary to structural defects in the lens capsule and developing corneal epithelium associated with basement membrane loss. Our analysis provides further insights into the role of the lens capsule and corneal basement membrane in the structural integrity of the developing eye.