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

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

Integrative transcriptomic profiling of ncRNAs and mRNAs in developing mouse lens

In recent years, burgeoning research has underscored the pivotal role of non-coding RNA in orchestrating the growth, development, and pathogenesis of various diseases across organisms. However, despite these advances, our understanding of the specific contributions of long non-coding RNAs (lncRNAs) and circular RNAs (circRNAs) to lens development remains notably limited. Clarifying the intricate gene regulatory networks is imperative for unraveling the molecular underpinnings of lens-related disorders. In this study, we aimed to address this gap by conducting a comprehensive analysis of the expression profiles of messenger RNAs (mRNAs), lncRNAs, and circRNAs at critical developmental time points of the mouse lens, encompassing both embryonic (E10.5, E12.5, and E16.5) and postnatal stages (P0.5, P10.5, and P60). Leveraging RNA-sequencing technology, we identified key transcripts pivotal to lens development. Our analysis revealed differentially expressed (DE) mRNAs, lncRNAs, and circRNAs across various developmental stages. Particularly noteworthy, there were 1831 co-differentially expressed (CO-DE) mRNAs, 150 CO-DE lncRNAs, and 13 CO-DE circRNAs identified during embryonic stages. Gene Ontology (GO) enrichment analysis unveiled associations primarily related to lens development, DNA conformational changes, and angiogenesis among DE mRNAs and lncRNAs. Furthermore, employing protein-protein interaction networks, mRNA-lncRNA co-expression networks, and circRNA-microRNA-mRNA networks, we predicted candidate key molecules implicated in lens development. Our findings underscore the pivotal roles of lncRNAs and circRNAs in this process, offering fresh insights into the pathogenesis of lens-related disorders and paving the way for future exploration in this field.

Multi-tissue transcriptome-wide association study identifies novel candidate susceptibility genes for cataract

Introduction

Cataract is the leading cause of blindness among the elderly worldwide. Twin and family studies support an important role for genetic factors in cataract susceptibility with heritability estimates up to 58%. To date, 55 loci for cataract have been identified by genome-wide association studies (GWAS), however, much work remains to identify the causal genes. Here, we conducted a transcriptome-wide association study (TWAS) of cataract to prioritize causal genes and identify novel ones, and examine the impact of their expression.

Methods

We performed tissue-specific and multi-tissue TWAS analyses to assess associations between imputed gene expression from 54 tissues (including 49 from the Genotype Tissue Expression (GTEx) Project v8) with cataract using FUSION software. Meta-analyzed GWAS summary statistics from 59,944 cataract cases and 478,571 controls, all of European ancestry and from two cohorts (GERA and UK Biobank) were used. We then examined the expression of the novel genes in the lens tissue using the iSyTE database.

Results

Across tissue-specific and multi-tissue analyses, we identified 99 genes for which genetically predicted gene expression was associated with cataract after correcting for multiple testing. Of these 99 genes, 20 (AC007773.1, ANKH, ASIP, ATP13A2, CAPZB, CEP95, COQ6, CREB1, CROCC, DDX5, EFEMP1, EIF2S2, ESRRB, GOSR2, HERC4, INSRR, NIPSNAP2, PICALM, SENP3, and SH3YL1) did not overlap with previously reported cataract-associated loci. Tissue-specific analysis identified 202 significant gene-tissue associations for cataract, of which 166 (82.2%), representing 9 unique genes, were attributed to the previously reported 11q13.3 locus. Tissue-enrichment analysis revealed that gastrointestinal tissues represented one of the highest proportions of the Bonferroni-significant gene-tissue associations (21.3%). Moreover, this gastrointestinal tissue type was the only anatomical category significantly enriched in our results, after correcting for the number of tissue donors and imputable genes for each reference panel. Finally, most of the novel cataract genes (e.g., Capzb) were robustly expressed in iSyTE lens data.

Discussion

Our results provide evidence of the utility of imputation-based TWAS approaches to characterize known GWAS risk loci and identify novel candidate genes that may increase our understanding of cataract etiology. Our findings also highlight the fact that expression of genes associated with cataract susceptibility is not necessarily restricted to lens tissue.

A Cataract-Causing Mutation in the TRPM3 Cation Channel Disrupts Calcium Dynamics in the Lens

TRPM3 belongs to the melastatin sub-family of transient receptor potential (TRPM) cation channels and has been shown to function as a steroid-activated, heat-sensitive calcium ion (Ca2+) channel. A missense substitution (p.I65M) in the TRPM3 gene of humans (TRPM3) and mice (Trpm3) has been shown to underlie an inherited form of early-onset, progressive cataract. Here, we model the pathogenetic effects of this cataract-causing mutation using 'knock-in' mutant mice and human cell lines. Trpm3 and its intron-hosted micro-RNA gene (Mir204) were strongly co-expressed in the lens epithelium and other non-pigmented and pigmented ocular epithelia. Homozygous Trpm3-mutant lenses displayed elevated cytosolic Ca2+ levels and an imbalance of sodium (Na+) and potassium (K+) ions coupled with increased water content. Homozygous TRPM3-mutant human lens epithelial (HLE-B3) cell lines and Trpm3-mutant lenses exhibited increased levels of phosphorylated mitogen-activated protein kinase 1/extracellular signal-regulated kinase 2 (MAPK1/ERK2/p42) and MAPK3/ERK1/p44. Mutant TRPM3-M65 channels displayed an increased sensitivity to external Ca2+ concentration and an altered dose response to pregnenolone sulfate (PS) activation. Trpm3-mutant lenses shared the downregulation of genes involved in insulin/peptide secretion and the upregulation of genes involved in Ca2+ dynamics. By contrast, Trpm3-deficient lenses did not replicate the pathophysiological changes observed in Trpm3-mutant lenses. Collectively, our data suggest that a cataract-causing substitution in the TRPM3 cation channel elicits a deleterious gain-of-function rather than a loss-of-function mechanism in the lens.

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.

Eye Lens Organoids Made Simple: Characterization of a New Three-Dimensional Organoid Model for Lens Development and Pathology

Cataract, the opacification of the lens, is the leading cause of blindness worldwide. Although effective, cataract surgery is costly and can lead to complications. Toward identifying alternate treatments, it is imperative to develop organoid models relevant for lens studies and drug screening. Here, we demonstrate that by culturing mouse lens epithelial cells under defined three-dimensional (3D) culture conditions, it is possible to generate organoids that display optical properties and recapitulate many aspects of lens organization and biology. These organoids can be rapidly produced in large amounts. High-throughput RNA sequencing (RNA-seq) on specific organoid regions isolated via laser capture microdissection (LCM) and immunofluorescence assays demonstrate that these lens organoids display a spatiotemporal expression of key lens genes, e.g., Jag1, Pax6, Prox1, Hsf4 and Cryab. Further, these lens organoids are amenable to the induction of opacities. Finally, the knockdown of a cataract-linked RNA-binding protein encoding gene, Celf1, induces opacities in these organoids, indicating their use in rapidly screening for genes that are functionally relevant to lens biology and cataract. In sum, this lens organoid model represents a compelling new tool to advance the understanding of lens biology and pathology and can find future use in the rapid screening of compounds aimed at preventing and/or treating cataracts.

Eye lens organoids going simple: characterization of a new 3-dimensional organoid model for lens development and pathology

The ocular lens, along with the cornea, focuses light on the retina to generate sharp images. Opacification of the lens, or cataract, is the leading cause of blindness worldwide. Presently, the best approach for cataract treatment is to surgically remove the diseased lens and replace it with an artificial implant. Although effective, this is costly and can have post-surgical complications. Toward identifying alternate treatments, it is imperative to develop organoid models relevant for lens studies and anti-cataract drug screening. Here, we demonstrate that by culturing mouse lens epithelial cells under defined 3-dimensional (3D) culture conditions, it is possible to generate organoids that display optical properties and recapitulate many aspects of lens organization at the tissue, cellular and transcriptomic levels. These 3D cultured lens organoids can be rapidly produced in large amounts. High-throughput RNA-sequencing (RNA-seq) on specific organoid regions isolated by laser capture microdissection (LCM) and immunofluorescence assays demonstrate that these lens organoids display spatiotemporal expression of key lens genes, e.g. , Jag1 , Pax6 , Prox1 , Hsf4 and Cryab . Further, these lens organoids are amenable to induction of opacities. Finally, knockdown of a cataract-linked RNA-binding protein encoding gene, Celf1 , induces opacities in these organoids, indicating their use in rapidly screening for genes functionally relevant to lens biology and cataract. In sum, this lens organoid model represents a compelling new tool to advance the understanding of lens biology and pathology, and can find future use in the rapid screening of compounds aimed at preventing and/or treating cataract.

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

Ocular lens development entails epithelial to fiber cell differentiation, defects in which cause congenital 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.

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

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

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

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

Lens Epithelial Explants Treated with Vitreous Humor Undergo Alterations in Chromatin Landscape with Concurrent Activation of Genes Associated with Fiber Cell Differentiation and Innate Immune Response

Lens epithelial explants are comprised of lens epithelial cells cultured in vitro on their native basement membrane, the lens capsule. Biologists have used lens epithelial explants to study many different cellular processes including lens fiber cell differentiation. In these studies, fiber differentiation is typically measured by cellular elongation and the expression of a few proteins characteristically expressed by lens fiber cells in situ. Chromatin and RNA was collected from lens epithelial explants cultured in either un-supplemented media or media containing 50% bovine vitreous humor for one or five days. Chromatin for ATAC-sequencing and RNA for RNA-sequencing was prepared from explants to assess regions of accessible chromatin and to quantitatively measure gene expression, respectively. Vitreous humor increased chromatin accessibility in promoter regions of genes associated with fiber differentiation and, surprisingly, an immune response, and this was associated with increased transcript levels for these genes. In contrast, vitreous had little effect on the accessibility of the genes highly expressed in the lens epithelium despite dramatic reductions in their mRNA transcripts. An unbiased analysis of differentially accessible regions revealed an enrichment of cis-regulatory motifs for RUNX, SOX and TEAD transcription factors that may drive differential gene expression in response to vitreous.

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

BACKGROUND

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

RESULTS

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

CONCLUSIONS

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

The lens epithelium as a major determinant in the development, maintenance, and regeneration of the crystalline lens

The crystalline lens is a transparent and refractive biconvex structure formed by lens epithelial cells (LECs) and lens fibers. Lens opacity, also known as cataracts, is the leading cause of blindness in the world. LECs are the principal cells of lens throughout human life, exhibiting different physiological properties and functions. During the embryonic stage, LECs proliferate and differentiate into lens fibers, which form the crystalline lens. Genetics and environment are vital factors that influence normal lens development. During maturation, LECs help maintain lens homeostasis through material transport, synthesis and metabolism as well as mitosis and proliferation. If disturbed, this will result in loss of lens transparency. After cataract surgery, the repair potential of LECs is activated and the structure and transparency of the regenerative tissue depends on postoperative microenvironment. This review summarizes recent research advances on the role of LECs in lens development, homeostasis, and regeneration, with a particular focus on the role of cholesterol synthesis (eg., lanosterol synthase) in lens development and homeostasis maintenance, and how the regenerative potential of LECs can be harnessed to develop surgical strategies and improve the outcomes of cataract surgery (Fig. 1). These new insights suggest that LECs are a major determinant of the physiological and pathological state of the lens. Further studies on their molecular biology will offer possibility to explore new approaches for cataract prevention and treatment.

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

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

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

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

Cell atlas of the human ocular anterior segment: Tissue-specific and shared cell types

The anterior segment of the eye consists of the cornea, iris, ciliary body, crystalline lens, and aqueous humor outflow pathways. Together, these tissues are essential for the proper functioning of the eye. Disorders of vision have been ascribed to defects in all of them; some disorders, including glaucoma and cataract, are among the most prevalent causes of blindness in the world. To characterize the cell types that compose these tissues, we generated an anterior segment cell atlas of the human eye using high-throughput single-nucleus RNA sequencing (snRNAseq). We profiled 195,248 nuclei from nondiseased anterior segment tissues of six human donors, identifying >60 cell types. Many of these cell types were discrete, whereas others, especially in the lens and cornea, formed continua corresponding to known developmental transitions that persist in adulthood. Having profiled each tissue separately, we performed an integrated analysis of the entire anterior segment, revealing that some cell types are unique to a single structure, whereas others are shared across tissues. The integrated cell atlas was then used to investigate cell type-specific expression patterns of more than 900 human ocular disease genes identified through either Mendelian inheritance patterns or genome-wide association studies.

Aged Lens Epithelial Cells Suppress Proliferation and Epithelial–Mesenchymal Transition-Relevance for Posterior Capsule Opacification

Posterior capsule opacification (PCO) is a frequent complication after cataract surgery, and advanced PCO requires YAG laser (Nd: YAG) capsulotomy, which often gives rise to more complications. Lens epithelial cell (LEC) proliferation and transformation (i.e., epithelial–mesenchymal transition (EMT)) are two critical elements in PCO initiation and progression pathogenesis. While PCO marginally impacts aged cataract surgery patients, PCO incidences are exceptionally high in infants and children undergoing cataract surgery. The gene expression of lens epithelial cell aging and its role in the discrepancy of PCO prevalence between young and older people have not been fully studied. Here, we conducted a comprehensive differentially expressed gene (DEG) analysis of a cell aging model by comparing the early and late passage FHL124 lens epithelial cells (LECs). In vitro, TGFβ2, cell treatment, and in vivo mouse cataract surgical models were used to validate our findings. We found that aged LECs decelerated rates of cell proliferation accompanied by dysregulation of cellular immune response and cell stress response. Surprisingly, we found that LECs systematically downregulated epithelial–mesenchymal transition (EMT)-promoting genes. The protein expression of several EMT hallmark genes, e.g., fibronectin, αSMA, and cadherin 11, were gradually decreased during LECs aging. We then confirmed these findings in vitro and found that aged LECs markedly alleviated TGFβ2-mediated EMT. Importantly, we explicitly confirmed the in vitro findings from the in vivo mouse cataract surgery studies. We propose that both the high proliferation rate and EMT-enriched young LECs phenotypic characteristics contribute to unusually high PCO incidence in infants and children.

Patterns of Crystallin Gene Expression in Differentiation State Specific Regions of the Embryonic Chicken Lens

Purpose

Transition from lens epithelial cells to lens fiber cell is accompanied by numerous changes in gene expression critical for lens transparency. We identify expression patterns of highly prevalent genes including ubiquitous and enzyme crystallins in the embryonic day 13 chicken lens.

Methods

Embryonic day 13 chicken lenses were dissected into central epithelial cell (EC), equatorial epithelial cell (EQ), cortical fiber cell (FP), and nuclear fiber cell (FC) compartments. Total RNA was prepared, subjected to high-throughput unidirectional mRNA sequencing, analyzed, mapped to the chicken genome, and functionally grouped.

Results

A total of 77,097 gene-specific transcripts covering 17,450 genes were expressed, of which 10,345 differed between two or more lens subregions. Ubiquitous crystallin gene expression increased from EC to EQ and was similar in FP and FC. Highly expressed crystallin genes fell into three coordinately expressed groups with R² ≥ 0.93: CRYAA, CRYBB2, CRYAB, and CRYBA2; CRYBB1, CRYBA4, CRYGN, ASL1, and ASL; and CRYBB3 and CRYBA1. The highly expressed transcription factors YBX1, YBX3, PNRC1, and BASP1 were coordinately expressed with the second group of crystallins (r² > 0.88).

Conclusions

Although it is well known that lens crystallin gene expression changes during the epithelial to fiber cell transition, these data identify for the first time three distinct patterns of expression for specific subsets of crystallin genes, each highly correlated with expression of specific transcription factors. The results provide a quantitative basis for designing functional experiments pinpointing the mechanisms governing the landscape of crystallin expression during fiber cell differentiation to attain lens transparency.

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.

Oxidative stress induces inflammation of lens cells and triggers immune surveillance of ocular tissues

Recent reports have challenged the notion that the lens is immune-privileged. However, these studies have not fully identified the molecular mechanism(s) that promote immune surveillance of the lens. Using a mouse model of targeted glutathione (GSH) deficiency in ocular surface tissues, we have investigated the role of oxidative stress in upregulating cytokine expression and promoting immune surveillance of the eye. RNA-sequencing of lenses from postnatal day (P) 1-aged Gclcf/f;Le-CreTg/- (KO) and Gclcf/f;Le-Cre-/- control (CON) mice revealed upregulation of many cytokines (e.g., CCL4, GDF15, CSF1) and immune response genes in the lenses of KO mice. The eyes of KO mice had a greater number of cells in the aqueous and vitreous humors at P1, P20 and P50 than age-matched CON and Gclcw/w;Le-CreTg/- (CRE) mice. Histological analyses revealed the presence of innate immune cells (i.e., macrophages, leukocytes) in ocular structures of the KO mice. At P20, the expression of cytokines and ROS content was higher in the lenses of KO mice than in those from age-matched CRE and CON mice, suggesting that oxidative stress may induce cytokine expression. In vitro administration of the oxidant, hydrogen peroxide, and the depletion of GSH (using buthionine sulfoximine (BSO)) in 21EM15 lens epithelial cells induced cytokine expression, an effect that was prevented by co-treatment of the cells with N-acetyl-l-cysteine (NAC), a antioxidant. The in vivo and ex vivo induction of cytokine expression by oxidative stress was associated with the expression of markers of epithelial-to-mesenchymal transition (EMT), α-SMA, in lens cells. Given that EMT of lens epithelial cells causes posterior capsule opacification (PCO), we propose that oxidative stress induces cytokine expression, EMT and the development of PCO in a positive feedback loop. Collectively these data indicate that oxidative stress induces inflammation of lens cells which promotes immune surveillance of ocular structures.

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.

An Atlas of Heparan Sulfate Proteoglycans in the Postnatal Rat Lens

Purpose

The arrangement of lens cells is regulated by ocular growth factors. Although the effects of these inductive molecules on lens cell behavior (proliferation, survival, and fiber differentiation) are well-characterized, the precise mechanisms underlying the regulation of growth factor-mediated signaling in lens remains elusive. Increasing evidence highlights the importance of heparan sulfate proteoglycans (HSPGs) for the signaling regulation of growth factors; however, the identity of the different lens HSPGs and the specific roles they play in lens biology are still unknown.

Methods

Semiquantitative real-time (RT)‐PCR and immunolabeling were used to characterize the spatial distribution of all known HSPG core proteins and their associated glycosaminoglycans (heparan and chondroitin sulfate) in the postnatal rat lens. Fibroblast growth factor (FGF)-2-treated lens epithelial explants, cultured in the presence of Surfen (an inhibitor of heparan sulfate [HS]-growth factor binding interactions) were used to investigate the requirement for HS in FGF-2-induced proliferation, fiber differentiation, and ERK1/2-signaling.

Results

The lens expresses all HSPGs. These HSPGs are differentially localized to distinct functional regions of the lens. In vitro, inhibition of HS-sulfation with Surfen blocked FGF-2-mediated ERK1/2-signaling associated with lens epithelial cell proliferation and fiber differentiation, highlighting that these cellular processes are dependent on HS.

Conclusions

These findings support a requirement for HSPGs in FGF-2 driven lens cell proliferation and fiber differentiation. The identification of specific HSPG core proteins in key functional lens regions, and the divergent expression patterns of closely related HSPGs, suggests that different HSPGs may differentially regulate growth factor signaling networks leading to specific biological events involved in lens growth and maintenance.

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.

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

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

Crystallin gene expression: Insights from studies of transcriptional bursting

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

Lens regeneration in humans: using regenerative potential for tissue repairing

The crystalline lens is an important optic element in human eyes. It is transparent and biconvex, refracting light and accommodating to form a clear retinal image. The lens originates from the embryonic ectoderm. The epithelial cells at the lens equator proliferate, elongate and differentiate into highly aligned lens fiber cells, which are the structural basis for maintaining the transparency of the lens. Cataract refers to the opacity of the lens. Currently, the treatment of cataract is to remove the opaque lens and implant an intraocular lens (IOL). This strategy is inappropriate for children younger than 2 years, because a developing eyeball is prone to have severe complications such as inflammatory proliferation and secondary glaucoma. On the other hand, the absence of the crystalline lens greatly affects visual function rehabilitation. The researchers found that mammalian lenses possess regenerative potential. We identified lens stem cells through linear tracking experiments and designed a minimally invasive lens-content removal surgery (MILS) to remove the opaque lens material while preserving the lens capsule, stem cells and microenvironment. In infants with congenital cataract, functional lens regeneration in situ can be observed after MILS, and the prognosis of visual function is better than that of traditional surgery. Because of insufficient regenerative ability in humans, the morphology and volume of the regenerated lens cannot reach the level of a normal lens. The activation, proliferation and differentiation of lens stem cells and the alignment of lens fibers are regulated by epigenetic factors, growth factors, transcription factors, immune system and other signals and their interactions. The construction of appropriate microenvironment can accelerate lens regeneration and improve its morphology. The therapeutic concept of MILS combined with microenvironment manipulation to activate endogenous stem cells for functional regeneration of organs in situ can be extended to other tissues and organs with strong self-renewal and repair ability.

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

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

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

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

Chicken lens development: complete signature of expression of galectins during embryogenesis and evidence for their complex formation with α-, β-, δ-, and τ-crystallins, N-CAM, and N-cadherin obtained by affinity chromatography

The emerging multifunctionality of galectins by specific protein-glycan/protein interactions explains the interest to determine their expression during embryogenesis. Complete network analysis of all seven chicken galectins (CGs) is presented in the course of differentiation of eye lens that originates from a single type of progenitor cell. It answers the questions on levels of expression and individual patterns of distribution. A qualitative difference occurs in the CG-1A/B paralogue pair, underscoring conspicuous divergence. Considering different cell phenotypes, lens fiber and also epithelial cells can both express the same CG, with developmental upregulation for CG-3 and CG-8. Except for expression of the lens-specific CG (C-GRIFIN), no other CG appeared to be controlled by the transcription factors L-Maf and Pax6. Studying presence and nature of binding partners for CGs, we tested labeled galectins in histochemistry and in ligand blotting. Mass spectrometric (glyco)protein identification after affinity chromatography prominently yielded four types of crystallins, N-CAM, and, in the cases of CG-3 and CG-8, N-cadherin. Should such pairing be functional in situ, it may be involved in tightly packing intracellular lens proteins and forming membrane contact as well as in gaining plasticity and stability of adhesion processes. The expression of CGs throughout embryogenesis is postulated to give meaning to spatiotemporal alterations in the local glycome.

Lens differentiation is characterized by stage-specific changes in chromatin accessibility correlating with differentiation state-specific gene expression

Changes in chromatin accessibility regulate the expression of multiple genes by controlling transcription factor access to key gene regulatory sequences. Here, we sought to establish a potential function for altered chromatin accessibility in control of key gene expression events during lens cell differentiation by establishing genome-wide chromatin accessibility maps specific for four distinct stages of lens cell differentiation and correlating specific changes in chromatin accessibility with genome-wide changes in gene expression. ATAC sequencing was employed to generate chromatin accessibility profiles that were correlated with the expression profiles of over 10,000 lens genes obtained by high-throughput RNA sequencing at the same stages of lens cell differentiation. Approximately 90,000 regions of the lens genome exhibited distinct changes in chromatin accessibility at one or more stages of lens differentiation. Over 1000 genes exhibited high Pearson correlation coefficients (r ​> ​0.7) between altered expression levels at specific stages of lens cell differentiation and changes in chromatin accessibility in potential promoter (-7.5kbp/+2.5kbp of the transcriptional start site) and/or other potential cis-regulatory regions ( ±10 ​kb of the gene body). Analysis of these regions identified consensus binding sequences for multiple transcription factors including members of the TEAD, FOX, and NFAT families of transcription factors as well as HIF1a, RBPJ and IRF1. Functional mapping of genes with high correlations between altered chromatin accessibility and differentiation state-specific gene expression changes identified multiple families of proteins whose expression could be regulated through changes in chromatin accessibility including those governing lens structure (BFSP1,BFSP2), gene expression (Pax-6, Sox 2), translation (TDRD7), cell-cell communication (GJA1), autophagy (FYCO1), signal transduction (SMAD3, EPHA2), and lens transparency (CRYBB1, CRYBA4). These data provide a novel relationship between altered chromatin accessibility and lens differentiation and they identify a wide-variety of lens genes and functions that could be regulated through altered chromatin accessibility. The data also point to a large number of potential DNA regulatory sequences and transcription factors whose functional analysis is likely to provide insight into novel regulatory mechanisms governing the lens differentiation program.

Bidirectional Analysis of Cryba4-Crybb1 Nascent Transcription and Nuclear Accumulation of Crybb3 mRNAs in Lens Fibers

Purpose

Crystallin gene expression during lens fiber cell differentiation is tightly spatially and temporally regulated. A significant fraction of mammalian genes is transcribed from adjacent promoters in opposite directions ("bidirectional" promoters). It is not known whether two proximal genes located on the same allele are simultaneously transcribed.

Methods

Mouse lens transcriptome was analyzed for paired genes whose transcriptional start sites are separated by less than 5 kbp to identify coexpressed bidirectional promoter gene pairs. To probe these transcriptional mechanisms, nascent transcription of Cryba4, Crybb1, and Crybb3 genes from gene-rich part of chromosome 5 was visualized by RNA fluorescent in situ hybridizations (RNA FISH) in individual lens fiber cell nuclei.

Results

Genome-wide lens transcriptome analysis by RNA-seq revealed that the Cryba4-Crybb1 pair has the highest Pearson correlation coefficient between their steady-state mRNA levels. Analysis of Cryba4 and Crybb1 nascent transcription revealed frequent simultaneous expression of both genes from the same allele. Nascent Crybb3 transcript visualization in "early" but not "late" differentiating lens fibers show nuclear accumulation of the spliced Crybb3 transcripts that was not affected in abnormal lens fiber cell nuclei depleted of chromatin remodeling enzyme Snf2h (Smarca5).

Conclusions

The current study shows for the first time that two highly expressed lens crystallin genes, Cryba4 and Crybb1, can be simultaneously transcribed from adjacent bidirectional promoters and do not show nuclear accumulation. In contrast, spliced Crybb3 mRNAs transiently accumulate in early lens fiber cell nuclei. The gene pairs coexpressed during lens development showed significant enrichment in human "cataract" phenotype.

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

BACKGROUND

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

RESULTS

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

CONCLUSIONS

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

Proteome-transcriptome analysis and proteome remodeling in mouse lens epithelium and fibers

Epithelial cells and differentiated fiber cells represent distinct compartments in the ocular lens. While previous studies have revealed proteins that are preferentially expressed in epithelial vs. fiber cells, a comprehensive proteomics library comparing the molecular compositions of epithelial vs. fiber cells is essential for understanding lens formation, function, disease and regenerative potential, and for efficient differentiation of pluripotent stem cells for modeling of lens development and pathology in vitro. To compare protein compositions between the lens epithelium and fibers, we employed tandem mass spectrometry (2D-LC/MS) analysis of microdissected mouse P0.5 lenses. Functional classifications of the top 525 identified proteins into gene ontology categories by molecular processes and subcellular localizations, were adapted for the lens. Expression levels of both epithelial and fiber proteomes were compared with whole lens proteome and mRNA levels using E14.5, E16.5, E18.5, and P0.5 RNA-Seq data sets. During this developmental time window, multiple complex biosynthetic and catabolic processes generate the molecular and structural foundation for lens transparency. As expected, crystallins showed a high correlation between their mRNA and protein levels. Comprehensive data analysis confirmed and/or predicted roles for transcription factors (TFs), RNA-binding proteins (e.g. Carhsp1), translational apparatus including ribosomal heterogeneity and initiation factors, microtubules, cytoskeletal [e.g. non-muscle myosin IIA heavy chain (Myh9) and βB2-spectrin (Sptbn2)] and membrane proteins in lens formation and maturation. Our data highlighted many proteins with unknown functions in the lens that were preferentially enriched in epithelium or fibers, setting the stage for future studies to further dissect the roles of these proteins in fiber cell differentiation vs. epithelial cell maintenance. In conclusion, the present proteomic datasets represent the first mouse lens epithelium and fiber cell proteomes, establish comparative analyses of protein and RNA-Seq data, and characterize the major proteome remodeling required to form the mature lens fiber cells.

Promoter-enhancer looping and shadow enhancers of the mouse αA-crystallin locus

Gene regulation by enhancers is important for precise temporal and spatial gene expression. Enhancers can drive gene expression regardless of their location, orientation or distance from the promoter. Changes in chromatin conformation and chromatin looping occur to bring the promoter and enhancers into close proximity. αA-crystallin ranks among one of the most abundantly expressed genes and proteins in the mammalian lens. The αA-crystallin locus is characterized by a 16 kb chromatin domain marked by two distal enhancers, 5′ DCR1 and 3′ DCR3. Here we used chromatin conformation capture (3C) analysis and transgenic approaches to analyze temporal control of the mouse αA-crystallin gene. We find that DCR1 is necessary, but not sufficient alone to drive expression at E10.5 in the mouse lens pit. Chromatin looping revealed interaction between the promoter and the region 3′ to DCR1, identifying a novel enhancer region in the αA-crystallin locus. We determined that this novel enhancer region, DCR1S, recapitulates the temporal control by DCR1. Acting as shadow enhancers, DCR1 and DCR1S are able to control expression in the lens vesicle at E11.5. It remains to be elucidated however, which region of the αA-crystallin locus is responsible for expression in the lens pit at E10.5.

Summary: The αA-crystallin ranks amongst the most highly expressed tissue-specific genes. It is an advantageous model system to probe both promoter-enhancer looping and to identify distal enhancers and their temporal/spatial activities.

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.

Identification of Novel Gata3 Distal Enhancers Active in Mouse Embryonic Lens

BACKGROUND

The tissue-specific transcriptional programs during normal development require tight control by distal cis-regulatory elements, such as enhancers, with specific DNA sequences recognized by transcription factors, coactivators, and chromatin remodeling enzymes. Gata3 is a sequence-specific DNA-binding transcription factor that regulates formation of multiple tissues and organs, including inner ear, lens, mammary gland, T-cells, urogenital system, and thyroid gland. In the eye, Gata3 has a highly restricted expression domain in the posterior part of the lens vesicle; however, the underlying regulatory mechanisms are unknown.

RESULTS

Here we describe the identification of a novel bipartite Gata3 lens-specific enhancer located ∼18 kb upstream from its transcriptional start site. We also found that a 5-kb Gata3 promoter possesses low activity in the lens. The bipartite enhancer contains arrays of AP-1, Ets-, and Smad1/5-binding sites as well as binding sites for lens-associated DNA-binding factors. Transient transfection studies of the promoter with the bipartite enhancer showed enhanced activation by BMP4 and FGF2.

CONCLUSIONS

These studies identify a novel distal enhancer of Gata3 with high activity in lens and indicate that BMP and FGF signaling can up-regulate expression of Gata3 in differentiating lens fiber cells through the identified Gata3 enhancer and promoter elements. Developmental Dynamics 247:1186-1198, 2018. © 2018 The Authors. Developmental Dynamics published by Wiley Periodicals, Inc. on behalf of American Association of Anatomists.