Transcriptome Profiling of Developing Murine Lens Through RNA Sequencing

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.

miR-26 Deficiency Causes Alterations in Lens Transcriptome and Results in Adult-Onset Cataract

Purpose

Despite strong evidence demonstrating that normal lens development requires regulation governed by microRNAs (miRNAs), the functional role of specific miRNAs in mammalian lens development remains largely unexplored.

Methods

A comprehensive analysis of miRNA transcripts in the newborn mouse lens, exploring both differential expression between lens epithelial cells and lens fiber cells and overall miRNA abundance, was conducted by miRNA sequencing. Mouse lenses lacking each of three abundantly expressed lens miRNAs (miR-184, miR-26, and miR-1) were analyzed to explore the role of these miRNAs in lens development.

Results

Mice lacking all three copies of miR-26 (miR-26TKO) developed postnatal cataracts as early as 4 to 6 weeks of age. RNA sequencing analysis of neonatal lenses from miR-26TKO mice exhibited abnormal reduced expression of a cohort of genes found to be lens enriched and linked to cataract (e.g., Foxe3, Hsf4, Mip, Tdrd7, and numerous crystallin genes) and abnormal elevated expression of genes related to neural development (Lhx3, Neurod4, Shisa7, Elavl3), inflammation (Ccr1, Tnfrsf12a, Csf2ra), the complement pathway, and epithelial to mesenchymal transition (Tnfrsf1a, Ccl7, Stat3, Cntfr).

Conclusions

miR-1, miR-184, and miR-26 are each dispensable for normal embryonic lens development. However, loss of miR-26 causes lens transcriptome changes and drives cataract formation.

Conditional Ablation of Spred1 and Spred2 in the Eye Lens Negatively Impacts Its Development and Growth

The development and growth of the eye depends on normal lens morphogenesis and its growth. This growth, in turn, is dependent on coordinated proliferation of the lens epithelial cells and their subsequent differentiation into fiber cells. These cellular processes are tightly regulated to maintain the precise cellular structure and size of the lens, critical for its transparency and refractive properties. Growth factor-mediated MAPK signaling driven by ERK1/2 has been reported as essential for regulating cellular processes of the lens, with ERK1/2 signaling tightly regulated by endogenous antagonists, including members of the Sprouty and related Spred families. Our previous studies have demonstrated the importance of both these inhibitory molecules in lens and eye development. In this study, we build on these findings to highlight the importance of Spreds in regulating early lens morphogenesis by modulating ERK1/2-mediated lens epithelial cell proliferation and fiber differentiation. Conditional loss of both Spred1 and Spred2 in early lens morphogenesis results in elevated ERK1/2 phosphorylation, hyperproliferation of lens epithelia, and an associated increase in the rate of fiber differentiation. This results in transient microphakia and microphthalmia, which disappears, owing potentially to compensatory Sprouty expression. Our data support an important temporal role for Spreds in the early stages of lens morphogenesis and highlight how negative regulation of ERK1/2 signaling is critical for maintaining lens proliferation and fiber differentiation in situ throughout life.

miR-26 deficiency causes alterations in lens transcriptome and results in adult-onset cataract

Purpose

Despite strong evidence demonstrating that normal lens development requires regulation governed by miRNAs, the functional role of specific miRNAs in mammalian lens development remains largely unexplored.

Methods

A comprehensive analysis of miRNA transcripts in the newborn mouse lens, exploring both differential expression between lens epithelial cells and lens fiber cells and overall miRNA abundance was conducted by miRNA-seq. Mouse lenses lacking each of three abundantly expressed lens miRNAs: miR-184, miR-26 and miR-1 were analyzed to explore the role of these miRNAs in lens development.

Results

Mice lacking all three copies of miR-26 (miR-26TKO) developed postnatal cataracts as early as 4-6 weeks of age. RNA-seq analysis of neonatal lenses from miR-26TKO mice exhibited abnormal reduced expression of a cohort of genes found to be lens-enriched and linked to cataract (e.g. Foxe3, Hsf4, Mip, Tdrd7, and numerous crystallin genes), and abnormal elevated expression of genes related to neural development (Lhx3, Neurod4, Shisa7, Elavl3 ), inflammation (Ccr1, Tnfrsf12a, Csf2ra), the complement pathway, and epithelial to mesenchymal transition (Tnfrsf1a, Ccl7, Stat3, Cntfr).

Conclusion

miR-1, miR-184 and miR-26 are each dispensable for normal embryonic lens development. However, loss of miR-26 causes lens transcriptome changes and drives cataract formation.

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

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

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

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

Spatiotemporal Localisation of Heparan Sulphate Proteoglycans throughout Mouse Lens Morphogenesis

Heparan sulphate proteoglycans (HSPGs) consist of a core protein decorated with sulphated HS-glycosaminoglycan (GAG) chains. These negatively charged HS-GAG chains rely on the activity of PAPSS synthesising enzymes for their sulfation, which allows them to bind to and regulate the activity of many positively charged HS-binding proteins. HSPGs are found on the surfaces of cells and in the pericellular matrix, where they interact with various components of the cell microenvironment, including growth factors. By binding to and regulating ocular morphogens and growth factors, HSPGs are positioned to orchestrate growth factor-mediated signalling events that are essential for lens epithelial cell proliferation, migration, and lens fibre differentiation. Previous studies have shown that HS sulfation is essential for lens development. Moreover, each of the full-time HSPGs, differentiated by thirteen different core proteins, are differentially localised in a cell-type specific manner with regional differences in the postnatal rat lens. Here, the same thirteen HSPG-associated GAGs and core proteins as well as PAPSS2, are shown to be differentially regulated throughout murine lens development in a spatiotemporal manner. These findings suggest that HS-GAG sulfation is essential for growth factor-induced cellular processes during embryogenesis, and the unique and divergent localisation of different lens HSPG core proteins indicates that different HSPGs likely play specialized roles during lens induction and morphogenesis.

Heparan sulfate proteoglycans (HSPGs) of the ocular lens

Heparan sulfate proteoglycans (HSPGs) reside in most cells; on their surface, in the pericellular milieu and/or extracellular matrix. In the eye, HSPGs can orchestrate the activity of key signalling molecules found in the ocular environment that promote its development and homeostasis. To date, our understanding of the specific roles played by individual HSPG family members, and the heterogeneity of their associated sulfated HS chains, is in its infancy. The crystalline lens is a relatively simple and well characterised ocular tissue that provides an ideal stage to showcase and model the expression and unique roles of individual HSPGs. Individual HSPG core proteins are differentially localised to eye tissues in a temporal and spatial developmental- and cell-type specific manner, and their loss or functional disruption results in unique phenotypic outcomes for the lens, and other ocular tissues. More recent work has found that different HS sulfation enzymes are also presented in a cell- and tissue-specific manner, and that disruption of these different sulfation patterns affects specific HS-protein interactions. Not surprisingly, these sulfated HS chains have also been reported to be required for lens and eye development, with dysregulation of HS chain structure and function leading to pathogenesis and eye-related phenotypes. In the lens, HSPGs undergo significant and specific changes in expression and function that can drive pathology, or in some cases, promote tissue repair. As master signalling regulators, HSPGs may one day serve as valuable biomarkers, and even as putative targets for the development of novel therapeutics, not only for the eye but for many other systemic pathologies.

In Silico Localization of Perilymph Proteins Enriched in Meńier̀e Disease Using Mammalian Cochlear Single-cell Transcriptomics

Hypothesis

Proteins enriched in the perilymph proteome of Meńier̀e disease (MD) patients may identify affected cell types. Utilizing single-cell transcriptome datasets from the mammalian cochlea, we hypothesize that these enriched perilymph proteins can be localized to specific cochlear cell types.

Background

The limited understanding of human inner ear pathologies and their associated biomolecular variations hinder efforts to develop disease-specific diagnostics and therapeutics. Perilymph sampling and analysis is now enabling further characterization of the cochlear microenvironment. Recently, enriched inner ear protein expression has been demonstrated in patients with MD compared to patients with other inner ear diseases. Localizing expression of these proteins to cochlear cell types can further our knowledge of potential disease pathways and subsequent development of targeted therapeutics.

Methods

We compiled previously published data regarding differential perilymph proteome profiles amongst patients with MD, otosclerosis, enlarged vestibular aqueduct, sudden hearing loss, and hearing loss of undefined etiology (controls). Enriched proteins in MD were cross-referenced against published single-cell/single-nucleus RNA-sequencing datasets to localize gene expression to specific cochlear cell types.

Results

In silico analysis of single-cell transcriptomic datasets demonstrates enrichment of a unique group of perilymph proteins associated with MD in a variety of intracochlear cells, and some exogeneous hematologic and immune effector cells. This suggests that these cell types may play an important role in the pathology associated with late MD, suggesting potential future areas of investigation for MD pathophysiology and treatment.

Conclusions

Perilymph proteins enriched in MD are expressed by specific cochlear cell types based on in silico localization, potentially facilitating development of disease-specific diagnostic markers and therapeutics.

The role of FYCO1-dependent autophagy in lens fiber cell differentiation

FYCO1 (FYVE and coiled-coil domain containing 1) is an adaptor protein, expressed ubiquitously and required for microtubule-dependent, plus-end-directed transport of macroautophagic/autophagic vesicles. We have previously shown that loss-of-function mutations in FYCO1 cause cataracts with no other ocular and/or extra-ocular phenotype. Here, we show fyco1 homozygous knockout (fyco1-/-) mice recapitulate the cataract phenotype consistent with a critical role of FYCO1 and autophagy in lens morphogenesis. Transcriptome coupled with proteome and metabolome profiling identified many autophagy-associated genes, proteins, and lipids respectively perturbed in fyco1-/- mice lenses. Flow cytometry of FYCO1 (c.2206C>T) knock-in (KI) human lens epithelial cells revealed a decrease in autophagic flux and autophagic vesicles resulting from the loss of FYCO1. Transmission electron microscopy showed cellular organelles accumulated in FYCO1 (c.2206C>T) KI lens-like organoid structures and in fyco1-/- mice lenses. In summary, our data confirm the loss of FYCO1 function results in a diminished autophagic flux, impaired organelle removal, and cataractogenesis.Abbreviations: CC: congenital cataracts; DE: differentially expressed; ER: endoplasmic reticulum; FYCO1: FYVE and coiled-coil domain containing 1; hESC: human embryonic stem cell; KI: knock-in; OFZ: organelle-free zone; qRT-PCR: quantitative real-time PCR; PE: phosphatidylethanolamine; RNA-Seq: RNA sequencing; SD: standard deviation; sgRNA: single guide RNA; shRNA: shorthairpin RNA; TEM: transmission electron microscopy; WT: wild type.

Transcriptional profiling of the chick retina identifies down-regulation of VIP and UTS2B genes during early lens-induced myopia

Gene expression of the chick retina was examined during the early development of lens-induced myopia (LIM) using whole transcriptome sequencing. Monocular treatment of the right eyes with −10 diopter (D)...

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.

Defect of LSS Disrupts Lens Development in Cataractogenesis

Congenital cataract is one of the leading causes of blindness in children worldwide. About one-third of congenital cataracts are caused by genetic defects. LSS, which encodes lanosterol synthase, is a causal gene for congenital cataracts. LSS is critical in preventing abnormal protein aggregation of various cataract-causing mutant crystallins; however, its roles in lens development remain largely unknown. In our study, we generated a mouse model harboring Lss G589S mutation, which is homologous to cataract-causing G588S mutation in human LSS. Lss^G589S/G589S mice exhibited neonatal lethality at postal day 0 (P0), whereas these mice showed severe opacity in eye lens. Also, we found that cataract was formed at E17.5 after we examined the opacity of embryonic lens from E13.5 to E18.5. Moreover, disrupted lens differentiation occurred at E14.5 prior to formation of the opacity of eye lens, shown as delayed differentiation of lens secondary fiber and disordered lens fiber organization. In addition, RNA-seq analysis indicated that cholesterol synthesis signaling pathways were significantly downregulated. Overall, our findings provide clear evidence that a mouse model harboring a homozygous Lss G589S mutation can recapitulate human congenital cataract. Our study points out that LSS functions as a critical determinant of lens development, which will contribute to better understanding LSS defects in cataractogenesis and developing therapies for cataracts.

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.

Examining the effects of cigarette smoke on mouse lens through a multi OMIC approach

Here, we report a multi OMIC (transcriptome, proteome, and metabolome) approach to investigate molecular changes in lens fiber cells (FC) of mice exposed to cigarette smoke (CS). Pregnant mice were placed in a whole-body smoke chamber and a few days later pups were born, which were exposed to CS for 5 hours/day, 5 days/week for a total of 3½ months. We examined the mice exposed to CS for CS-related cataractogenesis after completion of the CS exposure but no cataracts were observed. Lenses of CS-exposed and age-matched, untreated control mice were extracted and lens FC were subjected to multi OMIC profiling. We identified 348 genes, 130 proteins, and 14 metabolites exhibiting significant (p < 0.05) differential levels in lens FC of mice exposed to CS, corresponding to 3.6%, 4.3%, and 5.0% of the total genes, protein, and metabolites, respectively identified in this study. Our multi OMIC approach confirmed that only a small fraction of the transcriptome, the proteome, and the metabolome was perturbed in the lens FC of mice exposed to CS, which suggests that exposure of CS had a minimal effect on the mouse lens. It is worth noting that while our results confirm that CS exposure does not have a substantial impact on the molecular landscape of the mouse lens FC, we cannot rule out that CS exposure for longer durations and/or in combination with other morbidities or environmental factors would have a more robust effect and/or result in cataractogenesis.

The DNA methyltransferase DNMT3A contributes to autophagy long-term memory

Macroautophagy/autophagy is a conserved catabolic pathway that targets cytoplasmic components for their degradation and recycling in an autophagosome-dependent lysosomal manner. Under physiological conditions, this process maintains cellular homeostasis. However, autophagy can be stimulated upon different forms of cellular stress, ranging from nutrient starvation to exposure to drugs. Thus, this pathway can be seen as a central component of the integrated and adaptive stress response. Here, we report that even brief induction of autophagy is coupled in vitro to a persistent downregulation of the expression of MAP1LC3 isoforms, which are key components of the autophagy core machinery. In fact, DNA-methylation mediated by de novo DNA methyltransferase DNMT3A of MAP1LC3 loci upon autophagy stimulation leads to the observed long-term decrease of MAP1LC3 isoforms at transcriptional level. Finally, we report that the downregulation of MAP1LC3 expression can be observed in vivo in zebrafish larvae and mice exposed to a transient autophagy stimulus. This epigenetic memory of autophagy provides some understanding of the long-term effect of autophagy induction and offers a possible mechanism for its decline upon aging, pathological conditions, or in response to treatment interventions.Abbreviations: ACTB: actin beta; ATG: autophagy-related; 5-Aza: 5-aza-2'-deoxycytidine; BafA1: bafilomycin A1; CBZ: carbamazepine; CDKN2A: cyclin dependent kinase inhibitor 2A; ChIP: chromatin immunoprecipitation; Clon.: clonidine; CpG: cytosine-guanine dinucleotide: DMSO: dimethyl sulfoxide; DNA: deoxyribonucleic acid; DNMT: DNA methyltransferase; DNMT1: DNA methyltransferase 1; DNMT3A: DNA methyltransferase alpha; DNMT3B: DNA methyltransferase beta; dpf: days post-fertilization; EBSS: Earle's balanced salt solution; EM: Zebrafish embryo medium; GABARAP: GABA type A receptor associated protein; GABARAPL1: GABA type A receptor associated protein like 1; GABARAPL2: GABA type A receptor associated protein like 2; GAPDH: glyceraldehyde-3-phosphate dehydrogenase; GRO-Seq: Global Run-On sequencing; MAP1LC3/LC3: microtubule-associated protein 1 light chain 3; MAP1LC3A: microtubule-associated protein 1 light chain 3 alpha; MAP1LC3B: microtubule-associated protein 1 light chain 3 beta; MAP1LC3B2: microtubule-associated protein 1 light chain 3 beta 2; MEM: minimum essential medium; MEF: mouse embryonic fibroblasts; mRNA: messenger RNA; MTOR: mechanistic target of rapamycin kinase; PBS: phosphate-buffered saline; PIK3C3: phosphatidylinositol 3-kinase catalytic subunit type 3; RB1CC1/FIP200: RB1 inducible coiled-coil 1; RT-qPCR: quantitative reverse transcription polymerase chain reaction; SQSTM1/p62: sequestosome 1; Starv.: starvation; Treh.: trehalose; ULK1: unc-51 like autophagy activating kinase 1.

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

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

A missense allele of PEX5 is responsible for the defective import of PTS2 cargo proteins into peroxisomes

Peroxisomes, single-membrane intracellular organelles, play an important role in various metabolic pathways. The translocation of proteins from the cytosol to peroxisomes depends on peroxisome import receptor proteins and defects in peroxisome transport result in a wide spectrum of peroxisomal disorders. Here, we report a large consanguineous family with autosomal recessive congenital cataracts and developmental defects. Genome-wide linkage analysis localized the critical interval to chromosome 12p with a maximum two-point LOD score of 4.2 (θ = 0). Next-generation exome sequencing identified a novel homozygous missense variant (c.653 T > C; p.F218S) in peroxisomal biogenesis factor 5 (PEX5), a peroxisome import receptor protein. This missense mutation was confirmed by bidirectional Sanger sequencing. It segregated with the disease phenotype in the family and was absent in ethnically matched control chromosomes. The lens-specific knockout mice of Pex5 recapitulated the cataractous phenotype. In vitro import assays revealed a normal capacity of the mutant PEX5 to enter the peroxisomal Docking/Translocation Module (DTM) in the presence of peroxisome targeting signal 1 (PTS1) cargo protein, be monoubiquitinated and exported back into the cytosol. Importantly, the mutant PEX5 protein was unable to form a stable trimeric complex with peroxisomal biogenesis factor 7 (PEX7) and a peroxisome targeting signal 2 (PTS2) cargo protein and, therefore, failed to promote the import of PTS2 cargo proteins into peroxisomes. In conclusion, we report a novel missense mutation in PEX5 responsible for the defective import of PTS2 cargo proteins into peroxisomes resulting in congenital cataracts and developmental defects.

Metabolome profiling of the developing murine lens

Metabolomics is a study of the entire repertoire of metabolites in a cell at a particular time point. Here, we investigate the mouse lens at multiple embryonic and postnatal time points to establish the metabolome profile during early lens development. The lenses were isolated at six time points including embryonic day 15 (E15) and E18 and postnatal day 0 (P0), P3, P6, and P9. A total of four biological replicates of each time point, each consisting of 25 mg of lens tissue were preserved. Sample preparation was performed by protein precipitation followed by centrifugation to remove proteins and recover metabolites. The resulting extract was subjected to reverse phase/ultra-performance liquid chromatography-tandem mass spectrometry. Metabolome profiling identified a total of 353 metabolites in mouse lens, marked with an abundance of collagen, antioxidant, glycosaminoglycans, lipid, amino acid, and energy-related metabolites. A comparative metabolome analysis identified >200 metabolites exhibiting increased levels (p < 0.05) at latter time points relative to E15. Principal component analysis revealed distinct metabolomic signatures running from E15 to P9 while random forest analysis categorized lipid-, amino acid-, and nucleotide-related metabolites contributing significantly to the separation of the time points. To the best of our knowledge, this is the first report investigating the mouse lens metabolome at multiple embryonic and postnatal time points.

Differential miRNAs expression pattern of irradiated breast cancer cell lines is correlated with radiation sensitivity

Radiotherapy is a fundamental step in the treatment of breast cancer patients. The treatment efficiency is however reduced by the possible onset of radiation resistance. In order to develop the effective treatment approach, it is important to understand molecular basis of radiosensitivity in breast cancer. The purpose of the present study was to investigate different radiation response of breast cancer cell lines, and find out if this response may be related to change in the microRNAs expression profile. MDA-MB-231 and T47D cells were subjected to different doses of radiation, then MTT and clonogenic assays were performed to assess radiation sensitivity. Cytofluorometric and western blot analysis were performed to gain insight into cell cycle distribution and protein expression. MicroRNA sequencing and bioinformatics prediction methods were used to identify the difference in microRNAs expression between two breast cancer cells and the related genes and pathways. T47D cells were more sensitive to radiation respect to MDA-MB-231 cells as demonstrated by a remarkable G2 cell cycle arrest followed by a greater reduction in cell viability and colony forming ability. Accordingly, T47D cells showed higher increase in the phosphorylation of ATM, TP53 and CDK1 (markers of radiation response) and faster and more pronounced increase in RAD51 and γH2AX expression (markers of DNA damage), when compared to MDA-MB-231 cells. The two cell lines had different microRNAs expression profiles with a confirmed significant differential expression of miR-16-5p, which targets cell cycle related genes and predicts longer overall survival of breast cancer patients, as determined by bioinformatics analysis. These results suggest a possible role for miR-16-5p as radiation sensitizing microRNA and as prognostic/predictive biomarker in breast cancer.

Environmental and industrial developments in radiation cataractogenesis

Purpose: This review discusses recent developments in our understanding of biological and physiological mechanisms underlying radiation cataractogenesis. The areas discussed include effects of low-dose exposures to the lens including potential relevance of non-targeted effects, the development of new personal-protective equipment (PPE) and standards in clinical and nuclear settings motivated by the updated ICRP recommendations to mitigate exposures to the lens of the eye. The review also looks at evidence from the field linking cataracts in birds and mammals to low dose exposures.Conclusions: The review suggests that there is evidence that cataractogenesis is not a tissue reaction (deterministic effect) but rather is a low dose effect which shows a saturable dose response relationship similar to that seen for non-targeted effects in general. The review concludes that new research is needed to determine the dose response relationship in environmental studies where field data are contradictory and lab studies confined to rodent models for human exposure studies.

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

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

TRPM3_miR-204: a complex locus for eye development and disease

First discovered in a light-sensitive retinal mutant of Drosophila, the transient receptor potential (TRP) superfamily of non-selective cation channels serve as polymodal cellular sensors that participate in diverse physiological processes across the animal kingdom including the perception of light, temperature, pressure, and pain. TRPM3 belongs to the melastatin sub-family of TRP channels and has been shown to function as a spontaneous calcium channel, with permeability to other cations influenced by alternative splicing and/or non-canonical channel activity. Activators of TRPM3 channels include the neurosteroid pregnenolone sulfate, calmodulin, phosphoinositides, and heat, whereas inhibitors include certain drugs, plant-derived metabolites, and G-protein subunits. Activation of TRPM3 channels at the cell membrane elicits a signal transduction cascade of mitogen-activated kinases and stimulus response transcription factors. The mammalian TRPM3 gene hosts a non-coding microRNA gene specifying miR-204 that serves as both a tumor suppressor and a negative regulator of post-transcriptional gene expression during eye development in vertebrates. Ocular co-expression of TRPM3 and miR-204 is upregulated by the paired box 6 transcription factor (PAX6) and mutations in all three corresponding genes underlie inherited forms of eye disease in humans including early-onset cataract, retinal dystrophy, and coloboma. This review outlines the genomic and functional complexity of the TRPM3_miR-204 locus in mammalian eye development and disease.

Autosomal recessive congenital cataracts linked to HSF4 in a consanguineous Pakistani family

PURPOSE

To investigate the genetic basis of autosomal recessive congenital cataracts (arCC) in a large consanguineous Pakistani family.

METHODS

All participating members of family, PKCC074 underwent an ophthalmic examination. Slit-lamp photographs were ascertained for affected individuals that have not been operated for the removal of the cataractous lens. A small aliquot of the blood sample was collected from all participating individuals and genomic DNAs were extracted. A genome-wide scan was performed with polymorphic short tandem repeat (STR) markers and the logarithm of odds (LOD) scores were calculated. All coding exons and exon-intron boundaries of HSF4 were sequenced and expression of Hsf4 in mouse ocular lens was investigated. The C-terminal FLAG-tagged wild-type and mutant HSF4b constructs were prepared to examine the nuclear localization pattern of the mutant protein.

RESULTS

The ophthalmological examinations suggested that nuclear cataracts are present in affected individuals. Genome-wide linkage analyses localized the critical interval to a 10.95 cM (14.17 Mb) interval on chromosome 16q with a maximum two-point LOD score of 4.51 at θ = 0. Sanger sequencing identified a novel missense mutation: c.433G>C (p.Ala145Pro) that segregated with the disease phenotype in the family and was not present in ethnically matched controls. Real-time PCR analysis identified the expression of HSF4 in mouse lens as early as embryonic day 15 with a steady level of expression thereafter. The immunofluorescence tracking confirmed that both wild-type and mutant HSF4 (p.Ala145Pro) proteins localized to the nucleus.

CONCLUSION

Here, we report a novel missense mutation in HSF4 associated with arCC in a familial case of Pakistani descent.

Comparative transcriptome analysis of hESC- and iPSC-derived lentoid bodies

The ocular lens serves as an excellent system to investigate the intricate details of development and differentiation. Generation of lentoid bodies or lens-like structures using pluripotent stem cells is important for understanding the processes critical for lens morphogenesis and the mechanism of cataractogenesis. We previously reported the generation of peripheral blood mononuclear cell (PBMC)-originated, induced pluripotent stem cells (iPSCs). Here, we report generation of lentoid bodies from human embryonic stem cells (hESCs) and (PBMC)-originated, iPSCs employing the “fried egg” method with brief modifications. The ultrastructure analysis of hESC- and iPSC-derived lentoid bodies identified closely packed lens epithelial- and differentiating fiber-like cells. In addition, we performed RNA sequencing (RNA-Seq) based transcriptome profiling of hESC- and iPSC-derived lentoid bodies at differentiation day 25. Next-generation RNA sequencing (RNA-Seq) of hESC- and iPSC-derived lentoid bodies detected expression (≥0.659 RPKM) of 13,975 and 14,003 genes, respectively. Comparative transcriptome analysis of hESC- and iPSC-derived lentoid bodies revealed 13,563 (>96%) genes common in both datasets. Among the genes common in both transcriptome datasets, 12,856 (~95%) exhibited a quantitatively similar expression profile. Next, we compared the mouse lens epithelial and fiber cell transcriptomes with hESC- and iPSC-derived lentoid bodies transcriptomes and identified > 96% overlap with lentoid body transcriptomes. In conclusion, we report first-time comparative transcriptome analysis of hESC- and iPSC-derived lentoid bodies at differentiation day 25.

Cardiomyocyte dedifferentiation and remodeling in 3D scaffolds to generate the cellular diversity of engineering cardiac tissues

The use of engineered cardiac tissues (ECTs) is a new strategy for the repair and replacement of cardiac tissues in patients with myocardial infarction, particularly at late stages. However, the mechanisms underlying the development of ECTs, including cell-scaffold interactions, are not fully understood, although they are closely related to their therapeutic effect. In the present study, we aimed to determine the cellular fate of cardiomyocytes in a 3D scaffold microenvironment, as well as their role in generating the cellular diversity of ECTs by single-cell sequencing analysis. Consistent with the observed plasticity of cardiomyocytes during cardiac regeneration, cardiomyocytes in 3D scaffolds appeared to dedifferentiate, showing an initial loss of normal cytoskeleton organization in the adaptive response to the new scaffold microenvironment. Cardiomyocytes undergoing this process regained their proliferation potential and gradually developed into myocardial cells at different developmental stages, generating heterogeneous regenerative ECTs. To better characterize the remodeled ECTs, high-throughput single-cell sequencing was performed. The ECTs contained a wide diversity of cells related to endogenous classes in the heart, including myocardial cells at different developmental stages and different kinds of interstitial cells. Non-cardiac cells seemed to play important roles in cardiac reconstruction, especially Cajal-like interstitial cells and macrophages. Altogether, our results showed for the first time that cells underwent adaptive dedifferentiation for survival in a 3D scaffold microenvironment to generate heterogeneous tissues. These findings provide an important basis for an improved understanding of the development and assembly of engineered tissues.

Sall1 plays pivotal roles for lens fiber cell differentiation in mouse

Mammals possess four Sall transcription factors that play various roles in organogenesis. Previously, we found that Sall1 is expressed in microglia in the central nervous system, and it plays pivotal roles in microglia maturation. In the eye, Sall1 was also expressed in the developing lens, and we examined its role in lens development. A knock-in mouse harboring the EGFP gene in the Sall1 locus (Sall1-gfp) was used to analyze the Sall1 expression pattern. In Sall1-gfp/wild, EGFP was expressed throughout the presumptive lens at E11.5, and subsequently the expression in the lens epithelium became weaker. After birth, signals were observed in the equator region. The effects of Sall1 knockout on lens development were examined in Sall1-gfp/gfp. Lens sections revealed small vacuole-like holes and gaps in the center of the lens fibers at E14.5. Subsequently, the vacuoles appeared in most regions of the fiber cells. Electron microscopic analysis indicated that the vacuoles were between the fiber cells, leading to huge gaps. In addition, contact between the lens epithelium and apical end of the fiber cell was disrupted, and there were gaps between the adjoining lens epithelial cells. However, gap junction structure was observed by electron microscopic analysis, and immunostaining of Zo1 showed rather appropriate expression pattern. Immunohistochemistry indicated that the major lens transcription factors Prox1 and Pax6 were expressed in relatively normal patterns. However, although the expression of Prox1 and Pax6 decreased in nuclei in the control lens, it remained in Sall1-gfp/gfp. In addition, lower expression level of c-Maf protein was observed. Therefore, Sall1 is strongly expressed in the lens from the early developmental stage and plays an essential role in the maintenance of fiber cell and lens epithelium adhesion.

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.

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.

Dotting the eyes: mouse strain dependency of the lens epithelium to low dose radiation-induced DNA damage

PURPOSE

Epidemiological evidence regarding the radiosensitivity of the lens of the eye and radiation cataract development has led to changes in the EU Basic Safety Standards for protection of the lens against ionizing radiation. However, mechanistic details of lens radiation response pathways and their significance for cataractogenesis remain unclear. Radiation-induced DNA damage and the potential impairment of repair pathways within the lens epithelium, a cell monolayer that covers the anterior hemisphere of the lens, are likely to be involved.

MATERIALS AND METHODS

In this work, the lens epithelium has been analyzed for its DNA double-strand break (DSB) repair response to ionizing radiation. The responses of epithelial cells located at the anterior pole (central region) have been compared to at the very periphery of the monolayer (germinative and transitional zones). Described here are the different responses in the two regions and across four strains (C57BL/6, 129S2, BALB/c and CBA/Ca) over a low dose (0-25 mGy) in-vivo whole body X-irradiation range up to 24 hours post exposure.

RESULTS

DNA damage and repair as visualized through 53BP1 staining was present across the lens epithelium, although repair kinetics appeared non-uniform. Epithelial cells in the central region have significantly more 53BP1 foci. The sensitivities of different mouse strains have also been compared.

CONCLUSIONS

129S2 and BALB/c showed higher levels of DNA damage, with BALB/c showing significantly less inter-individual variability and appearing to be a more robust model for future DNA damage and repair studies. As a result of this study, BALB/c was identified as a suitable radiosensitive lens strain to detect and quantify early low dose ionizing radiation DNA damage effects in the mouse eye lens specifically, as an indicator of cataract formation.

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

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

Proteome Profiling of Developing Murine Lens Through Mass Spectrometry

Purpose

We previously completed a comprehensive profile of the mouse lens transcriptome. Here, we investigate the proteome of the mouse lens through mass spectrometry–based protein sequencing at the same embryonic and postnatal time points.

Methods

We extracted mouse lenses at embryonic day 15 (E15) and 18 (E18) and postnatal day 0 (P0), 3 (P3), 6 (P6), and 9 (P9). The lenses from each time point were preserved in three distinct pools to serve as biological replicates for each developmental stage. The total cellular protein was extracted from the lens, digested with trypsin, and labeled with isobaric tandem mass tags (TMT) for three independent TMT experiments.

Results

A total of 5404 proteins were identified in the mouse ocular lens in at least one TMT set, 4244 in two, and 3155 were present in all three TMT sets. The majority of the proteins exhibited steady expression at all six developmental time points; nevertheless, we identified 39 proteins that exhibited an 8-fold differential (higher or lower) expression during the developmental time course compared to their respective levels at E15. The lens proteome is composed of diverse proteins that have distinct biological properties and functional characteristics, including proteins associated with cataractogenesis and autophagy.

Conclusions

We have established a comprehensive profile of the developing murine lens proteome. This repository will be helpful in identifying critical components of lens development and processes essential for the maintenance of its transparency.

Identification of novel transcripts and peptides in developing murine lens

We previously investigated the transcriptome and proteome profiles of the murine ocular lens at six developmental time points including two embryonic (E15 and E18) and four postnatal time points (P0, P3, P6, and P9). Here, we extend our analyses to identify novel transcripts and peptides in developing  mouse lens. We identified a total of 9,707 novel transcripts and 325 novel fusion genes in developing mouse lens. Additionally, we identified 13,281 novel alternative splicing (AS) events in mouse lens including 6,990 exon skipping (ES), 2,447 alternative 3' splice site (A3SS), 1,900 alternative 5' splice site (A5SS), 1,771 mutually exclusive exons (MXE), and 173 intron retention (IR). Finally, we integrated our OMIC (Transcriptome and Proteome) datasets identifying 20 novel peptides in mouse lens. All 20 peptides were validated through matching MS/MS spectra of synthetic peptides. To the best of our knowledge, this is the first report integrating OMIC datasets to identify novel peptides in developing murine lens.

Probing the changes in gene expression due to α-crystallin mutations in mouse models of hereditary human cataract

The mammalian eye lens expresses a high concentration of crystallins (α, β and γ-crystallins) to maintain the refractive index essential for lens transparency. Crystallins are long-lived proteins that do not turnover throughout life. The structural destabilization of crystallins by UV exposure, glycation, oxidative stress and mutations in crystallin genes leads to protein aggregation and development of cataracts. Several destabilizing mutations in crystallin genes are linked with human autosomal dominant hereditary cataracts. To investigate the mechanism by which the α-crystallin mutations Cryaa-R49C and Cryab-R120G lead to cataract formation, we determined whether these mutations cause an altered expression of specific transcripts in the lens at an early postnatal age by RNA-seq analysis. Using knock-in mouse models previously generated in our laboratory, in the present work, we identified genes that exhibited altered abundance in the mutant lenses, including decreased transcripts for Clic5, an intracellular water channel in Cryaa-R49C heterozygous mutant lenses, and increased transcripts for Eno1b in Cryab-R120G heterozygous mutant lenses. In addition, RNA-seq analysis revealed increased histones H2B, H2A, and H4 gene expression in Cryaa-R49C mutant lenses, suggesting that the αA-crystallin mutation regulates histone expression via a transcriptional mechanism. Additionally, these studies confirmed the increased expression of histones H2B, H2A, and H4 by proteomic analysis of Cryaa-R49C knock-in and Cryaa;Cryab gene knockout lenses reported previously. Taken together, these findings offer additional insight into the early transcriptional changes caused by Cryaa and Cryab mutations associated with autosomal dominant human cataracts, and indicate that the transcript levels of certain genes are affected by the expression of mutant α-crystallin in vivo.

Express: A database of transcriptome profiles encompassing known and novel transcripts across multiple development stages in eye tissues

Advances in sequencing have facilitated nucleotide-resolution genome-wide transcriptomic profiles across multiple mouse eye tissues. However, these RNA sequencing (RNA-seq) based eye developmental transcriptomes are not organized for easy public access, making any further analysis challenging. Here, we present a new database "Express" (http://www.iupui.edu/∼sysbio/express/) that unifies various mouse lens and retina RNA-seq data and provides user-friendly visualization of the transcriptome to facilitate gene discovery in the eye. We obtained RNA-seq data encompassing 7 developmental stages of lens in addition to that on isolated lens epithelial and fibers, as well as on 11 developmental stages of retina/isolated retinal rod photoreceptor cells from publicly available wild-type mouse datasets. These datasets were pre-processed, aligned, quantified and normalized for expression levels of known and novel transcripts using a unified expression quantification framework. Express provides heatmap and browser view allowing easy navigation of the genomic organization of transcripts or gene loci. Further, it allows users to search candidate genes and export both the visualizations and the embedded data to facilitate downstream analysis. We identified total of >81,000 transcripts in the lens and >178,000 transcripts in the retina across all the included developmental stages. This analysis revealed that a significant number of the retina-expressed transcripts are novel. Expression of several transcripts in the lens and retina across multiple developmental stages was independently validated by RT-qPCR for established genes such as Pax6 and Lhx2 as well as for new candidates such as Elavl4, Rbm5, Pabpc1, Tia1 and Tubb2b. Thus, Express serves as an effective portal for analyzing pruned RNA-seq expression datasets presently collected for the lens and retina. It will allow a wild-type context for the detailed analysis of targeted gene-knockout mouse ocular defect models and facilitate the prioritization of candidate genes from Exome-seq data of eye disease patients.

Light-focusing human micro-lenses generated from pluripotent stem cells model lens development and drug-induced cataract in vitro

Cataracts cause vision loss and blindness by impairing the ability of the ocular lens to focus light onto the retina. Various cataract risk factors have been identified, including drug treatments, age, smoking and diabetes. However, the molecular events responsible for these different forms of cataract are ill-defined, and the advent of modern cataract surgery in the 1960s virtually eliminated access to human lenses for research. Here, we demonstrate large-scale production of light-focusing human micro-lenses from spheroidal masses of human lens epithelial cells purified from differentiating pluripotent stem cells. The purified lens cells and micro-lenses display similar morphology, cellular arrangement, mRNA expression and protein expression to human lens cells and lenses. Exposing the micro-lenses to the emergent cystic fibrosis drug Vx-770 reduces micro-lens transparency and focusing ability. These human micro-lenses provide a powerful and large-scale platform for defining molecular disease mechanisms caused by cataract risk factors, for anti-cataract drug screening and for clinically relevant toxicity assays.

A New Long Noncoding RNA ALB Regulates Autophagy by Enhancing the Transformation of LC3BI to LC3BII during Human Lens Development

Autophagy is essential in lens organelle degradation. This study aimed to seek potential autophagy-associated long noncoding RNAs (lncRNAs) and their relative mechanisms in human lens development using the “fried egg” lentoid body (LB) generation system. The expression pattern of LC3B in differentiating LBs was similar to that in developing a mouse lens in vivo. Among the massive lncRNAs expressed with a significant difference between induced pluripotent stem cells (iPSCs) and LBs, lncRNA affecting LC3B (ALB), which was predicted to have a co-relationship with the autophagy marker LC3B, was highly expressed in differentiating lens fibers in LBs. This result was consistent with its high expression in human embryonic lenses compared to those in embryonic stem cells (ESCs). Furthermore, lncRNA ALB knockdown resulted in the downregulation of LC3BII at the protein level, therefore inhibiting the autophagy process in human lens epithelial cells (HLECs). Our results identify lncRNA ALB, a potential autophagy regulator in organelle degradation during human lens development, highlighting the importance of lncRNAs in lens development.

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

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

Deletion at the GCNT2 Locus Causes Autosomal Recessive Congenital Cataracts

PURPOSE

The aim of this study is to identify the molecular basis of autosomal recessive congenital cataracts (arCC) in a large consanguineous pedigree.

METHODS

All participating individuals underwent a detailed ophthalmic examination. Each patient's medical history, particularly of cataracts and other ocular abnormalities, was compiled from available medical records and interviews with family elders. Blood samples were donated by all participating family members and used to extract genomic DNA. Genetic analysis was performed to rule out linkage to known arCC loci and genes. Whole-exome sequencing libraries were prepared and paired-end sequenced. A large deletion was found that segregated with arCC in the family, and chromosome walking was conducted to estimate the proximal and distal boundaries of the deletion mutation.

RESULTS

Exclusion and linkage analysis suggested linkage to a region of chromosome 6p24 harboring GCNT2 (glucosaminyl (N-acetyl) transferase 2) with a two-point logarithm of odds score of 5.78. PCR amplifications of the coding exons of GCNT2 failed in individuals with arCC, and whole-exome data analysis revealed a large deletion on chromosome 6p in the region harboring GCNT2. Chromosomal walking using multiple primer pairs delineated the extent of the deletion to approximately 190 kb. Interestingly, a failure to amplify a junctional fragment of the deletion break strongly suggests an insertion in addition to the large deletion.

CONCLUSION

Here, we report a novel insertion/deletion mutation at the GCNT2 locus that is responsible for congenital cataracts in a large consanguineous family.

Mutation in LIM2 Is Responsible for Autosomal Recessive Congenital Cataracts

Purpose

To identify the molecular basis of non-syndromic autosomal recessive congenital cataracts (arCC) in a consanguineous family.

Methods

All family members participating in the study received a comprehensive ophthalmic examination to determine their ocular phenotype and contributed a blood sample, from which genomic DNA was extracted. Available medical records and interviews with the family were used to compile the medical history of the family. The symptomatic history of the individuals exhibiting cataracts was confirmed by slit-lamp biomicroscopy. A genome-wide linkage analysis was performed to localize the disease interval. The candidate gene, LIM2 (lens intrinsic membrane protein 2), was sequenced bi-directionally to identify the disease-causing mutation. The physical changes caused by the mutation were analyzed in silico through homology modeling, mutation and bioinformatic algorithms, and evolutionary conservation databases. The physiological importance of LIM2 to ocular development was assessed in vivo by real-time expression analysis of Lim2 in a mouse model.

Results

Ophthalmic examination confirmed the diagnosis of nuclear cataracts in the affected members of the family; the inheritance pattern and cataract development in early infancy indicated arCC. Genome-wide linkage analysis localized the critical interval to chromosome 19q with a two-point logarithm of odds (LOD) score of 3.25. Bidirectional sequencing identified a novel missense mutation, c.233G>A (p.G78D) in LIM2. This mutation segregated with the disease phenotype and was absent in 192 ethnically matched control chromosomes. In silico analysis predicted lower hydropathicity and hydrophobicity but higher polarity of the mutant LIM2-encoded protein (MP19) compared to the wild-type. Moreover, these analyses predicted that the mutation would disrupt the secondary structure of a transmembrane domain of MP19. The expression of Lim2, which was detected in the mouse lens as early as embryonic day 15 (E15) increased after birth to a level that was sustained through the postnatal time points.

Conclusion

A novel missense mutation in LIM2 is responsible for autosomal recessive congenital cataracts.

A Common Ancestral Mutation in CRYBB3 Identified in Multiple Consanguineous Families with Congenital Cataracts

PURPOSE

This study was performed to investigate the genetic determinants of autosomal recessive congenital cataracts in large consanguineous families.

METHODS

Affected individuals underwent a detailed ophthalmological examination and slit-lamp photographs of the cataractous lenses were obtained. An aliquot of blood was collected from all participating family members and genomic DNA was extracted from white blood cells. Initially, a genome-wide scan was performed with genomic DNAs of family PKCC025 followed by exclusion analysis of our familial cohort of congenital cataracts. Protein-coding exons of CRYBB1, CRYBB2, CRYBB3, and CRYBA4 were sequenced bidirectionally. A haplotype was constructed with SNPs flanking the causal mutation for affected individuals in all four families, while the probability that the four familial cases have a common founder was estimated using EM and CHM-based algorithms. The expression of Crybb3 in the developing murine lens was investigated using TaqMan assays.

RESULTS

The clinical and ophthalmological examinations suggested that all affected individuals had nuclear cataracts. Genome-wide linkage analysis localized the causal phenotype in family PKCC025 to chromosome 22q with statistically significant two-point logarithm of odds (LOD) scores. Subsequently, we localized three additional families, PKCC063, PKCC131, and PKCC168 to chromosome 22q. Bidirectional Sanger sequencing identified a missense variation: c.493G>C (p.Gly165Arg) in CRYBB3 that segregated with the disease phenotype in all four familial cases. This variation was not found in ethnically matched control chromosomes, the NHLBI exome variant server, or the 1000 Genomes or dbSNP databases. Interestingly, all four families harbor a unique disease haplotype that strongly suggests a common founder of the causal mutation (p<1.64E-10). We observed expression of Crybb3 in the mouse lens as early as embryonic day 15 (E15), and expression remained relatively steady throughout development.

CONCLUSION

Here, we report a common ancestral mutation in CRYBB3 associated with autosomal recessive congenital cataracts identified in four familial cases of Pakistani origin.

FOXE3 contributes to Peters anomaly through transcriptional regulation of an autophagy-associated protein termed DNAJB1

FOXE3 is a lens-specific transcription factor that has been associated with anterior segment ocular dysgenesis. To determine the transcriptional target(s) of FOXE3 that are indispensable for the anterior segment development, we examined the transcriptome and the proteome of cells expressing truncated FOXE3 responsible for Peters anomaly identified through linkage-coupled next-generation whole-exome sequencing. We found that DNAJB1, an autophagy-associated protein, was the only candidate exhibiting differential expression in both screens. We confirmed the candidacy of DNAJB1 through chromatin immunoprecipitation and luciferase assays while knockdown of DNAJB1 in human lens epithelial cells resulted in a mitotic arrest. Subsequently, we targeted dnajb1a in zebrafish through injection of a splice-blocking morpholino. The dnajb1a morphants exhibited underdeveloped cataractous lenses with persistent apoptotic nuclei. In conclusion, here we report DNAJB1 is a transcriptional target of FOXE3 in a novel pathway that is crucial for the development of the anterior segment of the eye.

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

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

Non-coding RNA profiling of the developing murine lens

Non-coding RNAs (ncRNAs) are emerging as an important player in the regulation of genome integrity and gene expression, and they have been implicated in the pathogenesis of many diseases. The aim of the present study is to identify the repertoire of ncRNAs expressed in the developing mouse lens. We previously reported the mouse lens transcriptome, including mRNA and microRNA (miRNA) profiling at two embryonic (E15 and E18) and four postnatal (P0, P3, P6, and P9) time points. We analyzed the data from small RNA-Seq and mRNA-Seq libraries to investigate the ncRNA profile. Our analysis revealed expression of 12 different classes of ncRNA in the murine lens at six developmental time points. Annotation of small RNA data showed expression of 1,756 antisense ncRNA (asncRNA) in the mouse lens transcriptome. Likewise, we identified 82 P-element-induced wimpy testis (PIWI)-interacting RNA (piRNA), 345 transfer RNA (tRNA), 12 small nuclear RNA (snRNA), 167 small nucleolar RNA (snoRNA), 19 small Cajal body-specific RNA (scaRNA), six ribosomal RNA (rRNA), 18 tRNA-like structures, one MALAT1-associated small cytoplasmic RNA (mascRNA), one Vault RNA (vtRNA), and one Y RNA expressed in the developing mouse lens. In parallel, bioinformatic investigation of mRNA-Seq data identified expression of 1,952 long intergenic ncRNA (lincRNA) in the developing mouse lens. In conclusion, we report a comprehensive ncRNA profile in the murine lens at six developmental time points. To the best of our knowledge, this is first report investigating different classes of ncRNAs in the developing mouse lens and will be monumental in elucidating processes essential for the development of the ocular lens and the maintenance of its transparency.