Birc7: A Late Fiber Gene of the Crystalline Lens

Lens Regeneration: The Application of iSyTE and In Silico Approaches to Evaluate Gene Expression in Lens Organoids

Several protocols have been established for the generation of lens organoids from embryonic stem cells (ESCs), induced pluripotent stem cells (iPSCs), and other cells with regenerative potential in humans or various animal models. It is important to examine how well the regenerated lens organoids reflect lens biology, in terms of its development, homeostasis, and aging. Toward this goal, the iSyTE database (integrated Systems Tool for Eye gene discovery; https://research.bioinformatics.udel.edu/iSyTE/ ), a bioinformatics resource tool that contains meta-analyzed gene expression data in wild-type lens across different embryonic, postnatal, and adult stages, can serve as a resource for comparative analysis. This article outlines the approaches toward effective use of iSyTE to gain insights into normal gene expression in the mouse lens, enriched expression in the lens, and differential gene expression in select mouse gene-perturbation cataract/lens defects models, which in turn can be used to evaluate expression of key lens-relevant genes in lens organoids by transcriptomics (e.g., RNA-sequencing (RNA-seq), microarrays, etc.) or other downstream methods (e.g., RT-qPCR, etc.).

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.

Impaired GSH biosynthesis disrupts eye development, lens morphogenesis and PAX6 function

PURPOSE

The purpose of this study was to elucidate the role and molecular consequences of impaired glutathione (GSH) biosynthesis on eye development.

METHODS

GSH biosynthesis was impaired in surface ectoderm-derived ocular tissues by crossing Gclc^f/f mice with hemizygous Le-Cre transgenic mice to produce Gclc^f/f/Le-Cre^Tg/- (KO) mice. Control mice included Gclc^f/fand Gclc^wt/wt/Le-Cre^Tg/- mice (CRE). Eyes from all mice (at various stages of eye development) were subjected to histological, immunohistochemical, Western blot, RT-qPCR, RNA-seq, and subsequent Gene Ontology, Ingenuity Pathway Analysis and TRANSFAC analyses. PAX6 transactivation activity was studied using a luciferase reporter assay in HEK293T cells depleted of GSH using buthionine sulfoximine (BSO).

RESULTS

Deletion of Gclc diminished GSH levels, increased reactive oxygen species (ROS), and caused an overt microphthalmia phenotype characterized by malformation of the cornea, iris, lens, and retina that is distinct from and much more profound than the one observed in CRE mice. In addition, only the lenses of KO mice displayed reduced crystallin (α, β), PITX3 and Foxe3 expression. RNA-seq analyses at postnatal day 1 revealed 1552 differentially expressed genes (DEGs) in the lenses of KO mice relative to those from Gclc^f/f mice, with Crystallin and lens fiber cell identity genes being downregulated while lens epithelial cell identity and immune response genes were upregulated. Bioinformatic analysis of the DEGs implicated PAX6 as a key upstream regulator. PAX6 transactivation activity was impaired in BSO-treated HEK293T cells.

CONCLUSIONS

These data suggest that impaired ocular GSH biosynthesis may disrupt eye development and PAX6 function.

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.

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

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

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β.

High-throughput transcriptome analysis reveals that the loss of Pten activates a novel NKX6-1/RASGRP1 regulatory module to rescue microphthalmia caused by Fgfr2-deficient lenses

FGFR signaling is critical to development and disease pathogenesis, initiating phosphorylation-driven signaling cascades, notably the RAS-RAF-MEK-ERK and PI3 K-AKT cascades. PTEN antagonizes FGFR signaling by reducing AKT and ERK activation. Mouse lenses lacking FGFR2 exhibit microphakia and reduced ERK and AKT phosphorylation, widespread apoptosis, and defective lens fiber cell differentiation. In contrast, simultaneous deletion of both Fgfr2 and Pten restores ERK and AKT activation levels as well as lens size, cell survival and aspects of fiber cell differentiation; however, the molecular basis of this "rescue" remains undefined. We performed transcriptomic analysis by RNA sequencing of mouse lenses with conditional deletion of Fgfr2, Pten or both Fgfr2 and Pten, which reveal new molecular mechanisms that uncover how FGFR2 and PTEN signaling interact during development. The FGFR2-deficient lens transcriptome demonstrates overall loss of fiber cell identity with deregulated expression of 1448 genes. We find that ~ 60% of deregulated genes return to normal expression levels in lenses lacking both Fgfr2 and Pten. Further, application of customized filtering parameters to these RNA-seq data sets identified 68 high-priority candidate genes. Bioinformatics analyses showed that the cis-binding motif of a high-priority homeodomain transcription factor, NKX6-1, was present in the putative promoters of ~ 78% of these candidates. Finally, biochemical reporter assays demonstrate that NKX6-1 activated the expression of the high-priority candidate Rasgrp1, a RAS-activating protein. Together, these data define a novel regulatory module in which NKX6-1 directly activates Rasgrp1 expression to restore the balance of ERK and AKT activation, thus providing new insights into alternate regulation of FGFR downstream events.

Considerations for the use of Cre recombinase for conditional gene deletion in the mouse lens

Background

Despite a number of different transgenes that can mediate DNA deletion in the developing lens, each has unique features that can make a given transgenic line more or less appropriate for particular studies. The purpose of this work encompasses both a review of transgenes that lead to the expression of Cre recombinase in the lens and a comparative analysis of currently available transgenic lines with a particular emphasis on the Le-Cre and P0-3.9GFPCre lines that can mediate DNA deletion in the lens placode. Although both of these transgenes are driven by elements of the Pax6 P0 promoter, the Le-Cre transgene consistently leads to ocular abnormalities in homozygous state and can lead to ocular defects on some genetic backgrounds when hemizygous.

Result

Although both P0-3.9GFPCre and Le-Cre hemizygous transgenic mice undergo normal eye development on an FVB/N genetic background, Le-Cre homozygotes uniquely exhibit microphthalmia. Examination of the expression patterns of these two transgenes revealed similar expression in the developing eye and pancreas. However, lineage tracing revealed widespread non-ocular CRE reporter gene expression in the P0-3.9GFPCre transgenic mice that results from stochastic CRE expression in the P0-3.9GFPCre embryos prior to lens placode formation. Postnatal hemizygous Le-Cre transgenic lenses express higher levels of CRE transcript and protein than the hemizygous lenses of P0-3.9GFPCre mice. Transcriptome analysis revealed that Le-Cre hemizygous lenses deregulated the expression of 15 murine genes, several of which are associated with apoptosis. In contrast, P0-3.9GFPCre hemizygous lenses only deregulated two murine genes. No known PAX6-responsive genes or genes directly associated with lens differentiation were deregulated in the hemizygous Le-Cre lenses.

Conclusions

Although P0-3.9GFPCre transgenic mice appear free from ocular abnormalities, extensive non-ocular CRE expression represents a potential problem for conditional gene deletion studies using this transgene. The higher level of CRE expression in Le-Cre lenses versus P0-3.9GFPCre lenses may explain abnormal lens development in homozygous Le-Cre mice. Given the lack of deregulation of PAX6-responsive transcripts, we suggest that abnormal eye development in Le-Cre transgenic mice stems from CRE toxicity. Our studies reinforce the requirement for appropriate CRE-only expressing controls when using CRE as a driver of conditional gene targeting strategies.

Electronic supplementary material

The online version of this article (10.1186/s40246-019-0192-8) contains supplementary material, which is available to authorized users.

Spatial Analysis of Single Fiber Cells of the Developing Ocular Lens Reveals Regulated Heterogeneity of Gene Expression

The developing eye lens presents an exceptional paradigm for spatial transcriptomics. It is composed of highly organized long, slender transparent fiber cells, which differentiate from the edges of the anterior epithelium of the lens (equator), attended by high expression of crystallins, which generates transparency. Every fiber cell, therefore, is an optical unit whose refractive properties derive from its gene activity. Here, we probe this tangible relationship between the gene activity and the phenotype by studying the expression of all known 17 crystallins and 77 other non-crystallin genes in single fiber cells isolated from three states/regions of differentiation, allowing us to follow molecular progression at the single-cell level. The data demonstrate highly variable gene activity in cortical fibers, interposed between the nascent and the terminally differentiated fiber cell transcription. These data suggest that the so-called stochastic, highly heterogeneous gene activity is a regulated intermediate in the realization of a functional phenotype.

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.

IAP genes partake weighty roles in the astogeny and whole body regeneration in the colonial urochordate Botryllus schlosseri

Inhibitors of Apoptosis Protein (IAP) genes participate in processes like apoptosis, proliferation, innate immunity, inflammation, cell motility, differentiation and in malignancies. Here we reveal 25 IAP genes in the tunicate Botryllus schlosseri's genome and their functions in two developmental biology phenomena, a new mode of whole body regeneration (WBR) induced by budectomy, and blastogenesis, the four-staged cycles of botryllid ascidian astogeny. IAP genes that were specifically upregulated during these developmental phenomena were identified, and protein expression patterns of one of these genes, IAP28, were followed. Most of the IAP genes upregulation recorded at blastogenetic stages C/D was in concert with the upregulation at 100 μM H₂O₂ apoptotic-induced treatment and in parallel to expressions of AIF1, Bax, Mcl1, caspase 2 and two orthologues of caspase 7. Wnt agonist altered the takeover duration along with reduced IAP expressions, and displacement of IAP28⁺ phagocytes. WBR was initiated solely at blastogenetic stage D, where zooidal absorption was attenuated and regeneration centers were formed either from remains of partially absorbed zooids or from deformed ampullae. Subsequently, bud-bearing zooids developed, in concert with a massive IAP28-dependent phagocytic wave that eliminated the old zooids, then proceeded with the establishment of morphologically normal-looking colonies. IAP4, IAP14 and IAP28 were also involved in WBR, in conjunction with the expression of the pro-survival PI3K-Akt pathway. IAPs function deregulation by Smac mimetics resulted in severe morphological damages, attenuation in bud growth and differentiation, and in destabilization of colonial coordination. Longtime knockdown of IAP functions prior to the budectomy, resulted in colonial death.

Exposure to subthreshold dose of UVR-B induces apoptosis in the lens epithelial cells and does not in the lens cortical fibre cells

PURPOSE

The aim of this study was to investigate in which part of the lens in vivo exposure to subthreshold dose of UVR-B radiation induces apoptosis.

METHODS

Twenty 6-week-old female albino Sprague-Dawley rats were exposed to subthreshold dose (1 kJ/m² ) of UVR-B unilaterally and killed at 120 hr after exposure. Lenses were enucleated and dissected on three regions: the lens epithelium, the cortex and the nucleus. The lens nucleus then was removed. Apoptosis markers p53 and caspase 3 were used to study apoptosis in the lens regions. qRT-PCR and Western blot were utilized to analyse the lenses.

RESULTS

TP53 and CASP3 mRNA expressions are increased in exposed lenses, both in the lens epithelium and in the cortex regions, in relation to non-exposed lenses. Expression of p53 protein is increased in exposed lens epithelium in relation to non-exposed lens epithelium. Caspase 3 protein is expressed in exposed lens epithelial cells, while it is not expressed in non-exposed lens epithelial cells. p53 and caspase 3 proteins are not expressed in either exposed nor non-exposed lens fibre cells.

CONCLUSION

Exposure to UVR-B increases mRNA transcription of apoptosis marker p53 in vivo in both regions of the lens and of apoptosis marker caspase 3 in the lens cortex. Exposure to UVR-B increases p53 and caspase 3 proteins expression just in the lens epithelium. In vivo exposure to subthreshold dose of UVR-B induces apoptosis in the lens epithelial cells and does not in the lens fibre cells.

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.

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.