Evolutionary conservation of VSX2 super-enhancer modules in retinal development.

Super-enhancers (SEs) are expansive regions of genomic DNA that regulate the expression of genes involved in cell identity and cell fate. We recently identified developmental stage- and cell type-specific modules within the murine Vsx2 SE. Here, we show that the human VSX2 SE modules have similar developmental stage- and cell type-specific activity in reporter gene assays. By inserting the human sequence of one VSX2 SE module into a mouse with microphthalmia, eye size was rescued. To understand the function of these SE modules during human retinal development, we deleted individual modules in human embryonic stem cells and generated retinal organoids. Deleting one module results in small organoids, recapitulating the small-eyed phenotype of mice with microphthalmia, while deletion of the other module led to disruptions in bipolar neuron development. This prototypical SE serves as a model for understanding developmental stage- and cell type-specific effects of neurogenic transcription factors with complex expression patterns. Moreover, by elucidating the gene regulatory mechanisms, we can begin to examine how dysregulation of these mechanisms contributes to phenotypic diversity and disease.

Identification and Characterization of ATOH7-Regulated Target Genes and Pathways in Human Neuroretinal Development.

The proneural transcription factor atonal basic helix-loop-helix transcription factor 7 (ATOH7) is expressed in early progenitors in the developing neuroretina. In vertebrates, this is crucial for the development of retinal ganglion cells (RGCs), as mutant animals show an almost complete absence of RGCs, underdeveloped optic nerves, and aberrations in retinal vessel development. Human mutations are rare and result in autosomal recessive optic nerve hypoplasia (ONH) or severe vascular changes, diagnosed as autosomal recessive persistent hyperplasia of the primary vitreous (PHPVAR). To better understand the role of ATOH7 in neuroretinal development, we created ATOH7 knockout and eGFP-expressing ATOH7 reporter human induced pluripotent stem cells (hiPSCs), which were differentiated into early-stage retinal organoids. Target loci regulated by ATOH7 were identified by Cleavage Under Targets and Release Using Nuclease with sequencing (CUT&RUN-seq) and differential expression by RNA sequencing (RNA-seq) of wildtype and mutant organoid-derived reporter cells. Additionally, single-cell RNA sequencing (scRNA-seq) was performed on whole organoids to identify cell type-specific genes. Mutant organoids displayed substantial deficiency in axon sprouting, reduction in RGCs, and an increase in other cell types. We identified 469 differentially expressed target genes, with an overrepresentation of genes belonging to axon development/guidance and Notch signaling. Taken together, we consolidate the function of human ATOH7 in guiding progenitor competence by inducing RGC-specific genes while inhibiting other cell fates. Furthermore, we highlight candidate genes responsible for ATOH7-associated optic nerve and retinovascular anomalies, which sheds light to potential future therapy targets for related disorders.

Automated quantification of photoreceptor outer segments in developing and degenerating retinas on microscopy images across scales.

The functionality of photoreceptors, rods, and cones is highly dependent on their outer segments (POS), a cellular compartment containing highly organized membranous structures that generate biochemical signals from incident light. While POS formation and degeneration are qualitatively assessed on microscopy images, reliable methodology for quantitative analyses is still limited. Here, we developed methods to quantify POS (QuaPOS) maturation and quality on retinal sections using automated image analyses. POS formation was examined during the development and in adulthood of wild-type mice via light microscopy (LM) and transmission electron microscopy (TEM). To quantify the number, size, shape, and fluorescence intensity of POS, retinal cryosections were immunostained for the cone POS marker S-opsin. Fluorescence images were used to train the robust classifier QuaPOS-LM based on supervised machine learning for automated image segmentation. Characteristic features of segmentation results were extracted to quantify the maturation of cone POS. Subsequently, this quantification method was applied to characterize POS degeneration in "cone photoreceptor function loss 1" mice. TEM images were used to establish the ultrastructural quantification method QuaPOS-TEM for the alignment of POS membranes. Images were analyzed using a custom-written MATLAB code to extract the orientation of membranes from the image gradient and their alignment (coherency). This analysis was used to quantify the POS morphology of wild-type and two inherited retinal degeneration ("retinal degeneration 19" and "rhodopsin knock-out") mouse lines. Both automated analysis technologies provided robust characterization and quantification of POS based on LM or TEM images. Automated image segmentation by the classifier QuaPOS-LM and analysis of the orientation of membrane stacks by QuaPOS-TEM using fluorescent or TEM images allowed quantitative evaluation of POS formation and quality. The assessments showed an increase in POS number, volume, and membrane coherency during wild-type postnatal development, while a decrease in all three observables was detected in different retinal degeneration mouse models. All the code used for the presented analysis is open source, including example datasets to reproduce the findings. Hence, the QuaPOS quantification methods are useful for in-depth characterization of POS on retinal sections in developmental studies, for disease modeling, or after therapeutic interventions affecting photoreceptors.

A cis-regulatory module underlies retinal ganglion cell genesis and axonogenesis.

Atoh7 is transiently expressed in retinal progenitor cells (RPCs) and is required for retinal ganglion cell (RGC) differentiation. In humans, a deletion in a distal non-coding regulatory region upstream of ATOH7 is associated with optic nerve atrophy and blindness. Here, we functionally interrogate the significance of the Atoh7 regulatory landscape to retinogenesis in mice. Deletion of the Atoh7 enhancer structure leads to RGC deficiency, optic nerve hypoplasia, and retinal blood vascular abnormalities, phenocopying inactivation of Atoh7. Further, loss of the Atoh7 remote enhancer impacts ipsilaterally projecting RGCs and disrupts proper axonal projections to the visual thalamus. Deletion of the Atoh7 remote enhancer is also associated with the dysregulation of axonogenesis genes, including the derepression of the axon repulsive cue Robo3. Our data provide insights into how Atoh7 enhancer elements function to promote RGC development and optic nerve formation and highlight a key role of Atoh7 in the transcriptional control of axon guidance molecules.

Platelet-activating factor receptor (PAFR) regulates neuronal maturation and synaptic transmission during postnatal retinal development.

Introduction: Platelet-activating factor (PAF), PAF receptor (PAFR), and PAF- synthesis/degradation systems are involved in essential CNS processes such as neuroblast proliferation, differentiation, migration, and synaptic modulation. The retina is an important central nervous system (CNS) tissue for visual information processing. During retinal development, the balance between Retinal Progenitor Cell (RPC) proliferation and differentiation is crucial for proper cell determination and retinogenesis. Despite its importance in retinal development, the effects of PAFR deletion on RPC dynamics are still unknown. Methods: We compared PAFR knockout mice (PAFR-/-) retinal postnatal development proliferation and differentiation aspects with control animals. Electrophysiological responses were analyzed by electroretinography (ERG). Results and discussion: In this study, we demonstrate that PAFR-/- mice increased proliferation during postnatal retinogenesis and altered the expression of specific differentiation markers. The retinas of postnatal PAFR-/- animals decreased neuronal differentiation and synaptic transmission markers, leading to differential responses to light stimuli measured by ERG. Our findings suggest that PAFR signaling plays a critical role in regulating postnatal RPC cell differentiation dynamics during retinal development, cell organization, and neuronal circuitry formation.

Photoreceptor deficits appear at eye opening in Rs1 mutant mouse models of X-linked retinoschisis.

X-linked retinoschisis (XLRS) is an early onset degenerative retinal disease characterized by cystic lesions in the middle layers of the retina. These structural changes are accompanied by a loss of visual acuity and decreased contrast sensitivity. XLRS is caused by mutations in the gene Rs1 which encodes the secreted protein Retinoschisin 1. Young Rs1-mutant mouse models develop key hallmarks of XLRS including intraretinal schisis and abnormal electroretinograms. The electroretinogram (ERG) comprises activity of multiple cellular generators, and it is not known how and when each of these is impacted in Rs1 mutant mice. Here we use an ex vivo ERG system and pharmacological blockade to determine how ERG components generated by photoreceptors, ON-bipolar, and Müller glial cells are impacted in Rs1 mutants and to determine the time course of these changes. We report that ERG abnormalities begin near eye-opening and that all ERG components are involved.

Specific photoreceptor cell fate pathways are differentially altered in NR2E3-associated diseases.

Mutations in NR2E3, a gene encoding an orphan nuclear transcription factor, cause two retinal dystrophies with a distinct phenotype, but the precise role of NR2E3 in rod and cone transcriptional networks remains unclear. To dissect NR2E3 function, we performed scRNA-seq in the retinas of wildtype and two different Nr2e3 mouse models that show phenotypes similar to patients carrying NR2E3 mutations. Our results reveal that rod and cone populations are not homogeneous and can be separated into different sub-classes. We identify a previously unreported cone pathway that generates hybrid cones co-expressing both cone- and rod-related genes. In mutant retinas, this hybrid cone subpopulation is more abundant and includes a subpopulation of rods transitioning towards a cone cell fate. Hybrid photoreceptors with high misexpression of cone- and rod-related genes are prone to regulated necrosis. Overall, our results shed light on the role of NR2E3 in modulating photoreceptor differentiation towards cone and rod fates and explain how different mutations in NR2E3 lead to distinct visual disorders in humans.

Role of PKN1 in Retinal Cell Type Formation.

We recently identified PKN1 as a developmentally active gatekeeper of the transcription factor neuronal differentiation-2 (NeuroD2) in several brain areas. Since NeuroD2 plays an important role in amacrine cell (AC) and retinal ganglion cell (RGC) type formation, we aimed to study the expression of NeuroD2 in the postnatal retina of WT and Pkn1-/- animals, with a particular focus on these two cell types. We show that PKN1 is broadly expressed in the retina and that the gross retinal structure is not different between both genotypes. Postnatal retinal NeuroD2 levels were elevated upon Pkn1 knockout, with Pkn1-/- retinae showing more NeuroD2+ cells in the lower portion of the inner nuclear layer. Accordingly, immunohistochemical analysis revealed an increased amount of AC in postnatal and adult Pkn1-/- retinae. There were no differences in horizontal cell, bipolar cell, glial cell and RGC numbers, nor defective axon guidance to the optic chiasm or tract upon Pkn1 knockout. Interestingly, we did, however, see a specific reduction in SMI-32+ α-RGC in Pkn1-/- retinae. These results suggest that PKN1 is important for retinal cell type formation and validate PKN1 for future studies focusing on AC and α-RGC specification and development.

Maternal exposure to polystyrene nanoplastics causes defective retinal development and function in progeny mice by disturbing metabolic profiles.

Microplastics (MPs) and nanoplastics (NPs) are widely spreading in our living environment, accumulating in the human body and potentially threating human health. The retina, which is a terminally differentiated extension of the central nervous system, is essential for the visual system. However, the effects and molecular mechanisms of MPs/NPs on retina development and function are still unclear. Here, we investigated the effects and modes of action of polystyrene NPs (PS-NPs) on the retina using mice as a mammalian model species. Maternal PS-NP exposure (100 nm) at an environmentally realistic concentration of 10 mg L-1 (or 2.07 *1010 particles mL-1) via drinking water from the first day of pregnancy till the end of lactation (21 days after birth) caused defective neural retinal development in the neonatal mice, by depositing in the retinal tissue and reducing the number of retinal ganglion cells and bipolar cells. Exposure to PS-NPs retarded retinal vascular development, while abnormal electroretinogram (ERG) responses and an increased level of oxidative stress were also observed in the retina of the progeny mice after maternal PS-NP exposure. Metabolomics showed significant dysregulation of amino acids that are pivotal to neuron retinal function, such as glutamate, aspartate, alanine, glycine, serine, threonine, taurine, and serotonin. Transcriptomics identified significantly dysregulated genes, which were enriched in processes of angiogenesis, visual system development and lens development. Regulatory analysis showed that Fos gene mediated pathways could be a potential key target for PS-NP exposure in retinal development and function. Our study revealed that maternal exposure to PS-NPs generated detrimental effects on retinal development and function in progeny mice, offering new insights into the visual toxicity of PS-NPs.

The Interplay between Neurotransmitters and Calcium Dynamics in Retinal Synapses during Development, Health, and Disease.

The intricate functionality of the vertebrate retina relies on the interplay between neurotransmitter activity and calcium (Ca2+) dynamics, offering important insights into developmental processes, physiological functioning, and disease progression. Neurotransmitters orchestrate cellular processes to shape the behavior of the retina under diverse circumstances. Despite research to elucidate the roles of individual neurotransmitters in the visual system, there remains a gap in our understanding of the holistic integration of their interplay with Ca2+ dynamics in the broader context of neuronal development, health, and disease. To address this gap, the present review explores the mechanisms used by the neurotransmitters glutamate, gamma-aminobutyric acid (GABA), glycine, dopamine, and acetylcholine (ACh) and their interplay with Ca2+ dynamics. This conceptual outline is intended to inform and guide future research, underpinning novel therapeutic avenues for retinal-associated disorders.

RNA fusion in human retinal development.

Chimeric RNAs have been found in both cancerous and healthy human cells. They have regulatory effects on human stem/progenitor cell differentiation, stemness maintenance, and central nervous system development. However, whether they are present in human retinal cells and their physiological functions in the retinal development remain unknown. Based on the human embryonic stem cell-derived retinal organoids (ROs) spanning from days 0 to 120, we present the expression atlas of chimeric RNAs throughout the developing ROs. We confirmed the existence of some common chimeric RNAs and also discovered many novel chimeric RNAs during retinal development. We focused on CTNNBIP1-CLSTN1 (CTCL) whose downregulation caused precocious neuronal differentiation and a marked reduction of neural progenitors in human cerebral organoids. CTCL is universally present in human retinas, ROs, and retinal cell lines, and its loss-of-function biases the progenitor cells toward retinal pigment epithelial cell fate at the expense of retinal cells. Together, this work provides a landscape of chimeric RNAs and reveals evidence for their critical role in human retinal development.

Comparison of Retinal Metabolic Activity and Structural Development between rd10 Mice and Normal Mice Using Multiphoton Fluorescence Lifetime Imaging Microscopy.

Fluorescence lifetime imaging microscopy (FLIM) is a technique that analyzes the metabolic state of tissues based on the spatial distribution of fluorescence lifetimes of certain interacting molecules. We used multiphoton FLIM to study the metabolic state of developing C57BL6/J and rd10 retinas based on the fluorescence lifetimes of free versus bound nicotinamide adenine dinucleotide and nicotinamide adenine dinucleotide phosphate (NAD(P)H), with free NAD(P)H percentages suggesting increased glycolysis and bound NAD(P)H percentages indicating oxidative phosphorylation. The mice were sacrificed and enucleated at various time points throughout their first 3 months of life. The isolated eyecups were fixed, sectioned using a polyacrylamide gel embedding technique, and then analyzed with FLIM. The results suggested that in both C57BL6/J mice and rd10 mice, oxidative phosphorylation initially decreased and then increased, plateauing over time. This trend, however, was accelerated in rd10 mice, with its turning point occurring at p10 versus the p30 turning point in C57BL6/J mice. There was also a noticeable difference in oxidative phosphorylation rates between the outer and inner retinas in both strains, with greater oxidative phosphorylation present in the latter. A greater understanding of rd10 and WT metabolic changes during retinal development may provide deeper insights into retinal degeneration and facilitate the development of future treatments.

Expression patterns of CYP26A1, FGF8, CDKN1A, and NPVF in the developing rhesus monkey retina.

The fovea centralis (fovea) is a specialized region of the primate retina that plays crucial roles in high-resolution visual acuity and color perception. The fovea is characterized by a high density of cone photoreceptors and no rods, and unique anatomical properties that contribute to its remarkable visual capabilities. Early histological analyses identified some of the key events that contribute to foveal development, but the mechanisms that direct the specification of this area are not understood. Recently, the expression of the retinoic acid-metabolizing enzyme CYP26A1 has become a hallmark of some of the retinal specializations found in vertebrates, including the primate fovea and the high-acuity area in avian species. In chickens, the retinoic acid pathway regulates the expression of FGF8 to then direct the development of a rod-free area. Similarly, high levels of CYP26A1, CDKN1A, and NPVF expression have been observed in the primate macula using transcriptomic approaches. However, which retinal cells express these genes and their expression dynamics in the developing primate eye remain unknown. Here, we systematically characterize the expression patterns of CYP26A1, FGF8, CDKN1A, and NPVF during the development of the rhesus monkey retina, from early stages of development in the first trimester until the third trimester (near term). Our data suggest that some of the markers previously proposed to be fovea-specific are not enriched in the progenitors of the rhesus monkey fovea. In contrast, CYP26A1 is expressed at high levels in the progenitors of the fovea, while it localizes in a subpopulation of macular Müller glia cells later in development. Together these data provide invaluable insights into the expression dynamics of several molecules in the nonhuman primate retina and highlight the developmental advancement of the foveal region.

Stage-specific dynamic reorganization of genome topology shapes transcriptional neighborhoods in developing human retinal organoids.

We have generated a high-resolution Hi-C map of developing human retinal organoids to elucidate spatiotemporal dynamics of genomic architecture and its relationship with gene expression patterns. We demonstrate progressive stage-specific alterations in DNA topology and correlate these changes with transcription of cell-type-restricted gene markers during retinal differentiation. Temporal Hi-C reveals a shift toward A compartment for protein-coding genes and B compartment for non-coding RNAs, displaying high and low expression, respectively. Notably, retina-enriched genes are clustered near lost boundaries of topologically associated domains (TADs), and higher-order assemblages (i.e., TAD cliques) localize in active chromatin regions with binding sites for eye-field transcription factors. These genes gain chromatin contacts at their transcription start site as organoid differentiation proceeds. Our study provides a global view of chromatin architecture dynamics associated with diversification of cell types during retinal development and serves as a foundational resource for in-depth functional investigations of retinal developmental traits.

Identification of Evolutionarily Conserved VSX2 Enhancers in Retinal Development.

Super-enhancers (SEs) are expansive regions of genomic DNA that regulate the expression of genes involved in cell identity and cell fate. Recently, we found that distinct modules within a murine SE regulate gene expression of master regulatory transcription factor Vsx2 in a developmental stage- and cell-type specific manner. Vsx2 is expressed in retinal progenitor cells as well as differentiated bipolar neurons and Müller glia. Mutations in VSX2 in humans and mice lead to microphthalmia due to a defect in retinal progenitor cell proliferation. Deletion of a single module within the Vsx2 SE leads to microphthalmia. Deletion of a separate module within the SE leads to a complete loss of bipolar neurons, yet the remainder of the retina develops normally. Furthermore, the Vsx2 SE is evolutionarily conserved in vertebrates, suggesting that these modules are important for retinal development across species. In the present study, we examine the ability of these modules to drive retinal development between species. By inserting the human build of one Vsx2 SE module into a mouse with microphthalmia, eye size was rescued. To understand the implications of these SE modules in a model of human development, we generated human retinal organoids. Deleting one module results in small organoids, recapitulating the small-eyed phenotype of mice with microphthalmia, while deletion of the other module leads to a complete loss of ON cone bipolar neurons. This prototypical SE serves as a model for uncoupling developmental stage- and cell-type specific effects of neurogenic transcription factors with complex expression patterns. Moreover, by elucidating the gene regulatory mechanisms, we can begin to examine how dysregulation of these mechanisms contributes to phenotypic diversity and disease.

Effects of prenatal Doppler ultrasound on the retina of the chick embryo in ovo.

Currently, the effect of prenatal ultrasound on foetal development is intensively discussed and the guidelines for prenatal diagnostics have been changed. However, data supporting these concerns are scarce. Therefore, we used an established in ovo model of the chicken embryo to investigate cell proliferation and apoptosis within the retina. A total of 21 chicken eggs were fenestrated on Day 5 and allocated to either the control (n = 8) or exposition group (n = 13). The exposition group was treated with pulsed-wave Doppler ultrasound (PWD) for 10 min while controls remained without treatment. After subsequent incubation (6-48 h), the eggs were sacrificed, and chicken embryos were examined morphologically (HE-staining) and immunohistochemically. Counting of apoptotic and proliferating cells per retina was performed using antibodies specific for phospho-histone-H3 and active caspase-3 in combination with a biotin-labelled secondary antibody and peroxidase conjugated avidin-biotin complex for chromogenic detection. Due to a rather low number of specimens at each time point after ultrasound exposition, we neglected the effects of incubation time and focused on treatment effects. This approach revealed that the median number of proliferating cells is reduced after 10 min of exposure to PWD (569 vs. 766), while the number of apoptotic cells is fairly comparable between groups (5 vs. 6). Our data contribute to a better understanding of prenatal US on foetal development by suggesting that PWD could have an impact on the number of proliferating cells in the developing chicken retina and therefore justify further investigations.

Mitophagy in the retina: Viewing mitochondrial homeostasis through a new lens.

Mitochondrial function is key to support metabolism and homeostasis in the retina, an organ that has one of the highest metabolic rates body-wide and is constantly exposed to photooxidative damage and external stressors. Mitophagy is the selective autophagic degradation of mitochondria within lysosomes, and can be triggered by distinct stimuli such as mitochondrial damage or hypoxia. Here, we review the importance of mitophagy in retinal physiology and pathology. In the developing retina, mitophagy is essential for metabolic reprogramming and differentiation of retina ganglion cells (RGCs). In basal conditions, mitophagy acts as a quality control mechanism, maintaining a healthy mitochondrial pool to meet cellular demands. We summarize the different autophagy- and mitophagy-deficient mouse models described in the literature, and discuss the potential role of mitophagy dysregulation in retinal diseases such as glaucoma, diabetic retinopathy, retinitis pigmentosa, and age-related macular degeneration. Finally, we provide an overview of methods used to monitor mitophagy in vitro, ex vivo, and in vivo. This review highlights the important role of mitophagy in sustaining visual function, and its potential as a putative therapeutic target for retinal and other diseases.

The role of EZH2 in ocular diseases: a narrative review.

EZH2, acting as a catalytic subunit of PRC2 to catalyze lysine 27 in histone H3, induces the suppression of gene expression. EZH2 can regulate cell proliferation and differentiation of retinal progenitors, which are required for physiological retinal development. Meanwhile, an abnormal level of EZH2 has been observed in ocular tumors and other pathological tissues. This review summarizes the current knowledge on EZH2 in retinal development and ocular diseases, including inherited retinal diseases, ocular tumors, corneal injury, cataract, glaucoma, diabetic retinopathy and age-related retinal degeneration. We highlight the potential of targeting EZH2 as a precision therapeutic target in ocular diseases.

EZH2 is a protein that helps to regulate the activity of genes in cells. It works as a part of a complex called PRC2 to control a chemical group called lysine 27 in histone H3 and then inhibit the expression of genes. EZH2 is important for the normal development of the retina. Abnormal levels of EZH2 are associated with various eye diseases. This review summarizes the role of EZH2 in different ocular diseases and the potential mechanisms. Targeting EZH2 may be a novel way to treat or prevent ocular diseases.

Identification of the zebrafish homologues of IMPG2, a retinal proteoglycan.

Photoreceptor outer segments are surrounded by a carbohydrate-rich matrix, the interphotoreceptor matrix, necessary for physiological retinal function. Few roles for molecules characterizing the interphotoreceptor matrix have been clearly defined. Recent studies have found the presence of nonsense mutations in the interphotoreceptor matrix proteoglycan 2 (IMPG2) gene in patients affected by retinal dystrophies. IMPG2 encodes for a proteoglycan synthesized by photoreceptors and secreted in the interphotoreceptor matrix. Little is known about the structure and function of this protein, we thus decided to characterize zebrafish impg2. In zebrafish there are two Impg2 proteins, Impg2a and Impg2b. We generated a phylogenetic tree based on IMPG2 protein sequence similarity among vertebrates, showing a significant similarity between humans and teleosts. The human and zebrafish proteins share conserved domains, as also shown by homology models. Expression analyses of impg2a and impg2b show a continued expression in the photoreceptor layer starting from developmental stages and continuing through adulthood. Between 1 and 6 months post-fertilization, there is a significant shift of Impg2 expression toward the outer segment region, suggesting an increase in secretion. This raises intriguing hypotheses about its possible role(s) during retinal maturation, laying the groundwork for the generation of most needed models for the study of IMPG2-related inherited retinal dystrophies.

Mechanistic insight into the adverse outcome of tire wear and road particle leachate exposure in zebrafish (Danio rerio) larvae.

Tire wear particles (TWP) have become the major microplastic pollution in China. Road runoff containing TWP leachate can decrease the eye size and even induced mortality in the aquatic organisms. However, the toxic mechanism of TWP and road particles (RP) leachate on aquatic organisms is still unclear. In this study, the zebrafish embryos were exposed to TWP or RP leachate for 5 days at both environmental relevant and high concentrations. The adverse outcome pathways (AOPs) were screened from individual to molecular levels. The morphological and behavioral analysis demonstrated that the leachate exposure mainly impaired the eye development of zebrafish larvae and inhibited the larval swim behavior and phototactic response, which are the adverse outcomes. The phototransduction modulated by zebrafish retina was significantly down-regulated through transcriptomics and metabolomics analysis. The eye histopathological analysis showed that the decreased thickness of the retinal outer nuclear layer (ONL) and retinal pigmented epithelium (RPE) after leachate exposure were caused by the decreased photoreceptor cells. Moreover, the expression of NR2E3 and TPO genes showed concentration-dependent down-regulation after leachate exposure. The inhibition of photoreceptor cell proliferation was identified as the main reason for photoreceptor cell decrease in zebrafish larval eye. This study, for the first time, uncovered the underlying toxic mechanism of TWP and RP on zebrafish larval eyes.