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

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.

Multiomic analysis implicates FOXO4 in genetic regulation of chick lens fiber cell differentiation

A classic model for identification of novel differentiation mechanisms and pathways is the eye lens that consists of a monolayer of quiescent epithelial cells that are the progenitors of a core of mature fully differentiated fiber cells. The differentiation of lens epithelial cells into fiber cells follows a coordinated program involving cell cycle exit, expression of key structural proteins and the hallmark elimination of organelles to achieve transparency. Although multiple mechanisms and pathways have been identified to play key roles in lens differentiation, the entirety of mechanisms governing lens differentiation remain to be discovered. A previous study established that specific chromatin accessibility changes were directly associated with the expression of essential lens fiber cell genes, suggesting that the activity of transcription factors needed for expression of these genes could be regulated through binding access to the identified chromatin regions. Sequence analysis of the identified chromatin accessible regions revealed enhanced representation of the binding sequence for the transcription factor FOXO4 suggesting a direct role for FOXO4 in expression of these genes. FOXO4 is known to regulate a variety of cellular processes including cellular response to metabolic and oxidative stress, cell cycle withdrawal, and homeostasis, suggesting a previously unidentified role for FOXO4 in the regulation of lens cell differentiation. To further evaluate the role of FOXO4 we employed a multiomics approach to analyze the relationship between genome-wide FOXO4 binding, the differentiation-specific expression of key genes, and chromatin accessibility. To better identify active promoters and enhancers we also examined histone modification through analysis of H3K27ac. Specific methods included CUT&RUN (FOXO4 binding and H3K27ac modification), RNA-seq (differentiation state specific gene expression), and ATAC-seq (chromatin accessibility). CUT&RUN identified 20,966 FOXO4 binding sites and 33,921 H3K27ac marked regions across the lens fiber cell genome. RNA-seq identified 956 genes with significantly greater expression levels in fiber cells compared to epithelial cells (log2FC > 0.7, q < 0.05) and 2548 genes with significantly lower expression levels (log2FC < -0.7, q < 0.05). Integrated analysis identified 1727 differentiation-state specific genes that were nearest neighbors to at least one FOXO4 binding site, including genes encoding lens gap junctions (GJA1, GJA3), lens structural proteins (BFSP1, CRYBB1, ASL1), and genes required for lens transparency (HSF4, NRCAM). Multiomics analysis comparing the identified FOXO4 binding sites in published ATAC-seq data revealed that chromatin accessibility was associated with FOXO4-dependent gene expression during lens differentiation. The results provide evidence for an important requirement for FOXO4 in the regulated expression of key genes required for lens differentiation and link epigenetic regulation of chromatin accessibility and H3K27ac histone modification with the function of FOXO4 in controlling lens gene expression during lens fiber cell differentiation.

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.

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.

A tamoxifen-inducible Cre knock-in mouse for lens-specific gene manipulation

Mouse models are valuable tools in studying lens biology and biochemistry, and the Cre-loxP system is the most used technology for gene targeting in the lens. However, numerous genes are indispensable in lens development. The conventional knockout method either prevents lens formation or causes simultaneous cataract formation, hindering the studies of their roles in lens structure, growth, metabolism, and cataractogenesis during lens aging. An inducible Cre-loxP mouse line is an excellent way to achieve such a purpose. We established a lens-specific Cre ERT2 knock-in mouse (LCEK), an inducible mouse model for lens-specific gene targeting in a spatiotemporal manner. LCEK mice were created by in-frame infusion of a P2A-CreERT2 at the C-terminus of the last coding exon of the gene alpha A crystallin (Cryaa). LCEK mice express tamoxifen-inducible Cre recombinase uniquely in the lens. Through ROSAmT/mG and two endogenous genes (Gclc and Rbpj) targeting, we found no Cre recombinase leakage in the lens epithelium, but 50-80% leakage was observed in the lens cortex and nucleus. Administration of tamoxifen almost completely abolished target gene expression in both lens epithelium and cortex but only mildly enhanced gene deletion in the lens nucleus. Notably, no overt leakage of Cre activity was detected in developing LCEK lens when bred with mice carrying loxP floxed genes that are essential for lens development. This newly generated LCEK line will be a powerful tool to target genes in the lens for gene functions study in lens aging, posterior capsule opacification (PCO), and other areas requiring precision gene targeting.

Generation of Lens Progenitor Cells and Lentoid Bodies from Pluripotent Stem Cells: Novel Tools for Human Lens Development and Ocular Disease Etiology

In vitro differentiation of human pluripotent stem cells (hPSCs) into specialized tissues and organs represents a powerful approach to gain insight into those cellular and molecular mechanisms regulating human development. Although normal embryonic eye development is a complex process, generation of ocular organoids and specific ocular tissues from pluripotent stem cells has provided invaluable insights into the formation of lineage-committed progenitor cell populations, signal transduction pathways, and self-organization principles. This review provides a comprehensive summary of recent advances in generation of adenohypophyseal, olfactory, and lens placodes, lens progenitor cells and three-dimensional (3D) primitive lenses, "lentoid bodies", and "micro-lenses". These cells are produced alone or "community-grown" with other ocular tissues. Lentoid bodies/micro-lenses generated from human patients carrying mutations in crystallin genes demonstrate proof-of-principle that these cells are suitable for mechanistic studies of cataractogenesis. Taken together, current and emerging advanced in vitro differentiation methods pave the road to understand molecular mechanisms of cataract formation caused by the entire spectrum of mutations in DNA-binding regulatory genes, such as PAX6, SOX2, FOXE3, MAF, PITX3, and HSF4, individual crystallins, and other genes such as BFSP1, BFSP2, EPHA2, GJA3, GJA8, LIM2, MIP, and TDRD7 represented in human cataract patients.

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.

A functional map of genomic HIF1α-DNA complexes in the eye lens revealed through multiomics analysis

BACKGROUND

During eye lens development the embryonic vasculature regresses leaving the lens without a direct oxygen source. Both embryonically and throughout adult life, the lens contains a decreasing oxygen gradient from the surface to the core that parallels the natural differentiation of immature surface epithelial cells into mature core transparent fiber cells. These properties of the lens suggest a potential role for hypoxia and the master regulator of the hypoxic response, hypoxia-inducible transcription factor 1 (HIF1), in the regulation of genes required for lens fiber cell differentiation, structure and transparency. Here, we employed a multiomics approach combining CUT&RUN, RNA-seq and ATACseq analysis to establish the genomic complement of lens HIF1α binding sites, genes activated or repressed by HIF1α and the chromatin states of HIF1α-regulated genes.

RESULTS

CUT&RUN analysis revealed 8375 HIF1α-DNA binding complexes in the chick lens genome. One thousand one hundred ninety HIF1α-DNA binding complexes were significantly clustered within chromatin accessible regions (χ2 test p < 1 × 10- 55) identified by ATACseq. Formation of the identified HIF1α-DNA complexes paralleled the activation or repression of 526 genes, 116 of which contained HIF1α binding sites within 10kB of the transcription start sites. Some of the identified HIF1α genes have previously established lens functions while others have novel functions never before examined in the lens. GO and pathway analysis of these genes implicate HIF1α in the control of a wide-variety of cellular pathways potentially critical for lens fiber cell formation, structure and function including glycolysis, cell cycle regulation, chromatin remodeling, Notch and Wnt signaling, differentiation, development, and transparency.

CONCLUSIONS

These data establish the first functional map of genomic HIF1α-DNA complexes in the eye lens. They identify HIF1α as an important regulator of a wide-variety of genes previously shown to be critical for lens formation and function and they reveal a requirement for HIF1α in the regulation of a wide-variety of genes not yet examined for lens function. They support a requirement for HIF1α in lens fiber cell formation, structure and function and they provide a basis for understanding the potential roles and requirements for HIF1α in the development, structure and function of more complex tissues.

Mechanisms of organelle elimination for lens development and differentiation

A hallmark feature of lens development and differentiation is the complete elimination of organelles from the center of the eye lens. A long unanswered question in lens biology is what are the mechanisms that control the elimination of organelles during the terminal remodeling program to form mature lens fiber cells? Recent advances have expanded our understanding of these mechanisms including newly discovered signaling pathways, proteasomal regulators, autophagy proteins, transcription factors and the hypoxic environment of the lens itself. These recent discoveries suggest that distinct mechanisms coordinate the elimination of the nucleus, mitochondria, endoplasmic reticulum and Golgi apparatus during lens fiber cell differentiation. Since regulation of organelle number and distribution is also a feature of the terminal remodeling programs of more complex cell-types and tissues, these advances are likely to impact a wide-variety of fields.

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.

Ultrasound elastography for evaluating stiffness of the human lens nucleus with aging: a feasibility study

AIM

To investigate the significance of ultrasound elastography for evaluating stiffness of the human lens nucleus in volunteers with different ages.

METHODS

A total of 90 volunteers (lens transparency, uncorrected visual acuity ≥0.5, intraocular pressure: 14-19 mm Hg) were divided into 3 groups according to age: Group A (30 people, median age: 82±3.5y, mean axial lengths 23.7±0.5 mm); Group B (30 people, median age: 46±2.1y, mean axial lengths 23.9±0.4 mm); and Group C (30 people, median age: 22±3.5y, mean axial lengths 24.0±0.4 mm). Lens nuclear stiffness was measured by Free-hand qualitative elastography by independent operators. Strain gray scale and color-coded elastography maps were recorded. In each case, three consecutive detections were performed and strain ratio was used for statistical analysis.

RESULTS

Elastography analysis showed excellent diagnostic performance for lens sclerosis. Lens strain ratio was lowest (0.03±0.01)% in Group A and highest (2.03±0.43)% in Group C. Lens strain ratio was moderate (0.64±0.10)% in Group B. There were significant differences between these three groups (P<0.05). The lens nucleus strain rate changes with age. With aging, the lens nucleus strain rate and resilience decrease, demonstrating harder texture.

CONCLUSION

The relationship between human lens stiffness and age is demonstrated by ultrasound elastography. Older age is associated with lower strain ratio and less resilience of the lens.

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 Tudor-domain protein TDRD7, mutated in congenital cataract, controls the heat shock protein HSPB1 (HSP27) and lens fiber cell morphology

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

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.