Gene duplication and spectral diversification of cone visual pigments of zebrafish

Expansion and Functional Diversification of Long-Wavelength-Sensitive Opsin in Anabantoid Fishes

Gene duplication is one of the most important sources of novel genotypic diversity and the subsequent evolution of phenotypic diversity. Determining the evolutionary history and functional changes of duplicated genes is crucial for a comprehensive understanding of adaptive evolution. The evolutionary history of visual opsin genes is very dynamic, with repeated duplication events followed by sub- or neofunctionalization. While duplication of the green-sensitive opsins rh2 is common in teleost fish, fewer cases of multiple duplication events of the red-sensitive opsin lws are known. In this study, we investigate the visual opsin gene repertoire of the anabantoid fishes, focusing on the five lws opsin genes found in the genus Betta. We determine the evolutionary history of the lws opsin gene by taking advantage of whole-genome sequences of nine anabantoid species, including the newly assembled genome of Betta imbellis. Our results show that at least two independent duplications of lws occurred in the Betta lineage. The analysis of amino acid sequences of the lws paralogs of Betta revealed high levels of diversification in four of the seven transmembrane regions of the lws protein. Amino acid substitutions at two key-tuning sites are predicted to lead to differentiation of absorption maxima (λmax) between the paralogs within Betta. Finally, eye transcriptomics of B. splendens at different developmental stages revealed expression shifts between paralogs for all cone opsin classes. The lws genes are expressed according to their relative position in the lws opsin cluster throughout ontogeny. We conclude that temporal collinearity of lws expression might have facilitated subfunctionalization of lws in Betta and teleost opsins in general.

Light wavelength modulates search behavior performance in zebrafish

Visual systems have evolved to discriminate between different wavelengths of light. The ability to perceive color, or specific light wavelengths, is important as color conveys crucial information about both biotic and abiotic features in the environment. Indeed, different wavelengths of light can drive distinct patterns of activity in the vertebrate brain, yet what remains incompletely understood is whether distinct wavelengths can invoke etiologically relevant behavioral changes. To address how specific wavelengths in the visible spectrum modulate behavioral performance, we use larval zebrafish and a stereotypic light-search behavior. Prior work has shown that the cessation of light triggers a transitional light-search behavior, which we use to interrogate wavelength-dependent behavioral modulation. Using 8 narrow spectrum light sources in the visible range, we demonstrate that all wavelengths induce motor parameters consistent with search behavior, yet the magnitude of search behavior is spectrum sensitive and the underlying motor parameters are modulated in distinct patterns across short, medium, and long wavelengths. However, our data also establishes that not all motor features of search are impacted by wavelength. To define how wavelength modulates search performance, we performed additional assays with alternative wavelengths, dual wavelengths, and variable intensity. Last, we also tested blind larvae to resolve which components of wavelength dependent behavioral changes potentially include signaling from non-retinal photoreception. These findings have important implications as organisms can be exposed to varying wavelengths in laboratory and natural settings and therefore impose unique behavioral outputs.

Ancient origin of the rod bipolar cell pathway in the vertebrate retina

Vertebrates rely on rod photoreceptors for vision in low-light conditions. The specialized downstream circuit for rod signalling, called the primary rod pathway, is well characterized in mammals, but circuitry for rod signalling in non-mammals is largely unknown. Here we demonstrate that the mammalian primary rod pathway is conserved in zebrafish, which diverged from extant mammals ~400 million years ago. Using single-cell RNA sequencing, we identified two bipolar cell types in zebrafish that are related to mammalian rod bipolar cell (RBCs), the only bipolar type that directly carries rod signals from the outer to the inner retina in the primary rod pathway. By combining electrophysiology, histology and ultrastructural reconstruction of the zebrafish RBCs, we found that, similar to mammalian RBCs, both zebrafish RBC types connect with all rods in their dendritic territory and provide output largely onto amacrine cells. The wiring pattern of the amacrine cells postsynaptic to one RBC type is strikingly similar to that of mammalian RBCs and their amacrine partners, suggesting that the cell types and circuit design of the primary rod pathway emerged before the divergence of teleost fish and mammals. The second RBC type, which forms separate pathways, was either lost in mammals or emerged in fish.

Knockout of lws1 in zebrafish (Danio rerio) reveals its role in regulating feeding and vision-guided behavior

Long-wave sensitive (LWS) is a G protein-coupled receptor expressed in the retina, and zebrafish is a better model organism for studying vision, but the role of LWS1 in vision-guided behavior of larvae fish has rarely been reported. In this study, we found that zebrafish lws1 and lws2 are tandemly replicated genes, both with six exons, with lws1 being more evolutionarily conserved. The presence of Y277F in the amino acid sequence of lws2 may have contributed to the shift of λmax to green light. We established a lws1 knockout zebrafish model using CRISPR/Cas9 technology. Lws1-/- larvae showed significantly higher levels of feeding and appetite gene (agrp) expression than WT, and significantly lower levels of anorexia gene (pomc, cart) expression. In addition, green light gene compensation was observed in lws1-/- larvae with significantly increased expression levels of rh2-1. The light-dark movement test showed that lws1-/- larvae were more active under light-dark transitions or vibrational stimuli, and the expression of phototransduction-related genes was significantly up-regulated. This study reveals the important role of lws1 gene in the regulation of vision-guided behavior in larvae.

Color vision evolution in egg-laying mammals: insights from visual photoreceptors and daily activities of Australian echidnas

Egg-laying mammals (monotremes) are considered "primitive" due to traits such as oviparity, cloaca, and incomplete homeothermy, all of which they share with reptiles. Two groups of monotremes, the terrestrial echidna (Tachyglossidae) and semiaquatic platypus (Ornithorhynchidae), have evolved highly divergent characters since their emergence in the Cenozoic era. These evolutionary differences, notably including distinct electrosensory and chemosensory systems, result from adaptations to species-specific habitat conditions. To date, very few studies have examined the visual adaptation of echidna and platypus. In the present study, we show that echidna and platypus have different light absorption spectra in their dichromatic visual sensory systems at the molecular level. We analyzed absorption spectra of monotreme color opsins, long-wavelength sensitive opsin (LWS) and short-wavelength sensitive opsin 2 (SWS2). The wavelength of maximum absorbance (λmax) in LWS was 570.2 in short-beaked echidna (Tachyglossus aculeatus) and 560.6 nm in platypus (Ornithorhynchus anatinus); in SWS2, λmax was 451.7 and 442.6 nm, respectively. Thus, the spectral range in echidna color vision is ~ 10 nm longer overall than in platypus. Natural selection analysis showed that the molecular evolution of monotreme color opsins is generally functionally conserved, suggesting that these taxa rely on species-specific color vision. In order to understand the usage of color vision in monotremes, we made 24-h behavioral observations of captive echidnas at warm temperatures and analyzed the resultant ethograms. Echidnas showed cathemeral activity and various behavioral repertoires such as feeding, traveling, digging, and self-grooming without light/dark environment selectivity. Halting (careful) behavior is more frequent in dark conditions, which suggests that echidnas may be more dependent on vision during the day and olfaction at night. Color vision functions have contributed to dynamic adaptations and dramatic ecological changes during the ~ 60 million years of divergent monotreme evolution. The ethogram of various day and night behaviors in captive echidnas also contributes information relevant to habitat conservation and animal welfare in this iconic species, which is locally endangered.

Retinoic acid signaling regulates spatiotemporal specification of human green and red cones

Trichromacy is unique to primates among placental mammals, enabled by blue (short/S), green (medium/M), and red (long/L) cones. In humans, great apes, and Old World monkeys, cones make a poorly understood choice between M and L cone subtype fates. To determine mechanisms specifying M and L cones, we developed an approach to visualize expression of the highly similar M- and L-opsin mRNAs. M-opsin was observed before L-opsin expression during early human eye development, suggesting that M cones are generated before L cones. In adult human tissue, the early-developing central retina contained a mix of M and L cones compared to the late-developing peripheral region, which contained a high proportion of L cones. Retinoic acid (RA)-synthesizing enzymes are highly expressed early in retinal development. High RA signaling early was sufficient to promote M cone fate and suppress L cone fate in retinal organoids. Across a human population sample, natural variation in the ratios of M and L cone subtypes was associated with a noncoding polymorphism in the NR2F2 gene, a mediator of RA signaling. Our data suggest that RA promotes M cone fate early in development to generate the pattern of M and L cones across the human retina.

Ancient origin of the rod bipolar cell pathway in the vertebrate retina

Vertebrates rely on rod photoreceptors for vision in low-light conditions. Mammals have a specialized downstream circuit for rod signaling called the primary rod pathway, which comprises specific cell types and wiring patterns that are thought to be unique to this lineage. Thus, it has been long assumed that the primary rod pathway evolved in mammals. Here, we challenge this view by demonstrating that the mammalian primary rod pathway is conserved in zebrafish, which diverged from extant mammals ~400 million years ago. Using single-cell RNA-sequencing, we identified two bipolar cell (BC) types in zebrafish that are related to mammalian rod BCs (RBCs) of the primary rod pathway. By combining electrophysiology, histology, and ultrastructural reconstruction of the zebrafish RBCs, we found that, like mammalian RBCs, both zebrafish RBC types connect with all rods in their dendritic territory, and provide output largely onto amacrine cells. The wiring pattern of the amacrine cells post-synaptic to one RBC type is strikingly similar to that of mammalian RBCs, suggesting that the cell types and circuit design of the primary rod pathway have emerged before the divergence of teleost fish and amniotes. The second RBC type, which forms separate pathways, is either lost in mammals or emerged in fish.

Ancient origin of the rod bipolar cell pathway in the vertebrate retina

Vertebrates rely on rod photoreceptors for vision in low-light conditions1. Mammals have a specialized downstream circuit for rod signaling called the primary rod pathway, which comprises specific cell types and wiring patterns that are thought to be unique to this lineage2-6. Thus, it has been long assumed that the primary rod pathway evolved in mammals3,5-7. Here, we challenge this view by demonstrating that the mammalian primary rod pathway is conserved in zebrafish, which diverged from extant mammals ~400 million years ago. Using single-cell RNA-sequencing, we identified two bipolar cell (BC) types in zebrafish that are related to mammalian rod BCs (RBCs) of the primary rod pathway. By combining electrophysiology, histology, and ultrastructural reconstruction of the zebrafish RBCs, we found that, like mammalian RBCs8, both zebrafish RBC types connect with all rods and red-cones in their dendritic territory, and provide output largely onto amacrine cells. The wiring pattern of the amacrine cells post-synaptic to one RBC type is strikingly similar to that of mammalian RBCs. This suggests that the cell types and circuit design of the primary rod pathway may have emerged before the divergence of teleost fish and amniotes (mammals, bird, reptiles). The second RBC type in zebrafish, which forms separate pathways from the first RBC type, is either lost in mammals or emerged in fish to serve yet unknown roles.

Plasticity of cone photoreceptors in adult zebrafish revealed by thyroid hormone exposure

Vertebrate color vision is predominantly mediated by the presence of multiple cone photoreceptor subtypes that are each maximally sensitive to different wavelengths of light. Thyroid hormone (TH) has been shown to be essential in the spatiotemporal patterning of cone subtypes in many species, including cone subtypes that express opsins that are encoded by tandemly replicated genes. TH has been shown to differentially regulate the tandemly replicated lws opsin genes in zebrafish, and exogenous treatments alter the expression levels of these genes in larvae and juveniles. In this study, we sought to determine whether gene expression in cone photoreceptors remains plastic to TH treatment in adults. We used a transgenic lws reporter line, multiplexed fluorescence hybridization chain reaction in situ hybridization, and qPCR to examine the extent to which cone gene expression can be altered by TH in adults. Our studies revealed that opsin gene expression, and the expression of other photoreceptor genes, remains plastic to TH treatment in adult zebrafish. In addition to retinal plasticity, exogenous TH treatment alters skin pigmentation patterns in adult zebrafish after 5 days. Taken together, our results show a remarkable level of TH-sensitive plasticity in the adult zebrafish.

Long wavelength-sensing cones of zebrafish retina exhibit multiple layers of transcriptional heterogeneity

Introduction

Understanding how photoreceptor genes are regulated is important for investigating retinal development and disease. While much is known about gene regulation in cones, the mechanism by which tandemly-replicated opsins, such as human long wavelength-sensitive and middle wavelength-sensitive opsins, are differentially regulated remains elusive. In this study, we aimed to further our understanding of transcriptional heterogeneity in cones that express tandemly-replicated opsins and the regulation of such differential expression using zebrafish, which express the tandemly-replicated opsins lws1 and lws2.

Methods

We performed bulk and single cell RNA-Seq of LWS1 and LWS2 cones, evaluated expression patterns of selected genes of interest using multiplex fluorescence in situ hybridization, and used exogenous thyroid hormone (TH) treatments to test selected genes for potential control by thyroid hormone: a potent, endogenous regulator of lws1 and lws2 expression.

Results

Our studies indicate that additional transcriptional differences beyond opsin expression exist between LWS1 and LWS2 cones. Bulk RNA-Seq results showed 95 transcripts enriched in LWS1 cones and 186 transcripts enriched in LWS2 cones (FC > 2, FDR < 0.05). In situ hybridization results also reveal underlying heterogeneity within the lws1- and lws2-expressing populations. This heterogeneity is evident in cones of mature zebrafish, and further heterogeneity is revealed in transcriptional responses to TH treatments.

Discussion

We found some evidence of coordinate regulation of lws opsins and other genes by exogenous TH in LWS1 vs. LWS2 cones, as well as evidence of gene regulation not mediated by TH. The transcriptional differences between LWS1 and LWS2 cones are likely controlled by multiple signals, including TH.

Transcription factors underlying photoreceptor diversity

During development, retinal progenitors navigate a complex landscape of fate decisions to generate the major cell classes necessary for proper vision. Transcriptional regulation is critical to generate diversity within these major cell classes. Here, we aim to provide the resources and techniques required to identify transcription factors necessary to generate and maintain diversity in photoreceptor subtypes, which are critical for vision. First, we generate a key resource: a high-quality and deep transcriptomic profile of each photoreceptor subtype in adult zebrafish. We make this resource openly accessible, easy to explore, and have integrated it with other currently available photoreceptor transcriptomic datasets. Second, using our transcriptomic profiles, we derive an in-depth map of expression of transcription factors in photoreceptors. Third, we use efficient CRISPR-Cas9 based mutagenesis to screen for null phenotypes in F0 larvae (F0 screening) as a fast, efficient, and versatile technique to assess the involvement of candidate transcription factors in the generation of photoreceptor subtypes. We first show that known phenotypes can be easily replicated using this method: loss of S cones in foxq2 mutants and loss of rods in nr2e3 mutants. We then identify novel functions for the transcription factor Tbx2, demonstrating that it plays distinct roles in controlling the generation of all photoreceptor subtypes within the retina. Our study provides a roadmap to discover additional factors involved in this process. Additionally, we explore four transcription factors of unknown function (Skor1a, Sall1a, Lrrfip1a, and Xbp1), and find no evidence for their involvement in the generation of photoreceptor subtypes. This dataset and screening method will be a valuable way to explore the genes involved in many other essential aspects of photoreceptor biology.

Bisphenol S disrupts opsins gene expression and impairs the light-sensing function via antagonizing TH-TRβ signaling pathway in zebrafish larvae

Bisphenol S (BPS) is extensively used in "bisphenol A-free" products such as baby bottles. Although the visual toxicity of BPS has been reported, the underlying mechanism was largely unknown. In the present study, zebrafish were exposed to 0, 4 and 400 nM BPS from 2 h post-fertilization (hpf) to 120 hpf to further explore the thyroid disruption mechanism underlying the BPS induced impairment of visual function. The results showed that BPS decreased T3 levels in larval eyes, induced retinal expression of thyroid hormone receptor β (TRβ), and thereby down-regulated the expression of TH-mediated opsin genes (opn1lw1, opn1lw2, opn1mw1, opn1mw2, opn1mw3, and opn1sw2) and impaired subsequent phototransduction pathways, leading to decreased visually mediated phototactic response and body color adaptation but stimulated visual motor response (VMR). Combining exposure of exogenous T3 or 1-850 (antagonist for TRβ) with BPS could partly compensate the inhibited expression of opsin genes (opn1mw2, opn1lw1, and opn1lw2) and alleviate the hyperactivity of larval VMR caused by BPS alone, suggesting that BPS disrupted the opsins expression and also light-sensing function via antagonizing TH-TRβ signaling pathway. This study underlined the importance of TH signaling in regulating the proper vision and proposed a novel mechanism for the visual toxicity of BPS.

From perception to behavior: The neural circuits underlying prey hunting in larval zebrafish

A key challenge for neural systems is to extract relevant information from the environment and make appropriate behavioral responses. The larval zebrafish offers an exciting opportunity for studying these sensing processes and sensory-motor transformations. Prey hunting is an instinctual behavior of zebrafish that requires the brain to extract and combine different attributes of the sensory input and form appropriate motor outputs. Due to its small size and transparency the larval zebrafish brain allows optical recording of whole-brain activity to reveal the neural mechanisms involved in prey hunting and capture. In this review we discuss how the larval zebrafish brain processes visual information to identify and locate prey, the neural circuits governing the generation of motor commands in response to prey, how hunting behavior can be modulated by internal states and experience, and some outstanding questions for the field.

The evolutionary history and spectral tuning of vertebrate visual opsins

Many animals depend on the sense of vision for survival. In eumetazoans, vision requires specialized, light-sensitive cells called photoreceptors. Light reaches the photoreceptors and triggers the excitation of light-detecting proteins called opsins. Here, we describe the story of visual opsin evolution from the ancestral bilaterian to the extant vertebrate lineages. We explain the mechanisms determining color vision of extant vertebrates, focusing on opsin gene losses, duplications, and the expression regulation of vertebrate opsins. We describe the sequence variation both within and between species that has tweaked the sensitivities of opsin proteins towards different wavelengths of light. We provide an extensive resource of wavelength sensitivities and mutations that have diverged light sensitivity in many vertebrate species and predict how these mutations were accumulated in each lineage based on parsimony. We suggest possible natural and sexual selection mechanisms underlying these spectral differences. Understanding how molecular changes allow for functional adaptation of animals to different environments is a major goal in the field, and therefore identifying mutations affecting vision and their relationship to photic selection pressures is imperative. The goal of this review is to provide a comprehensive overview of our current understanding of opsin evolution in vertebrates.

Optic flow in the natural habitats of zebrafish supports spatial biases in visual self-motion estimation

Animals benefit from knowing if and how they are moving. Across the animal kingdom, sensory information in the form of optic flow over the visual field is used to estimate self-motion. However, different species exhibit strong spatial biases in how they use optic flow. Here, we show computationally that noisy natural environments favor visual systems that extract spatially biased samples of optic flow when estimating self-motion. The performance associated with these biases, however, depends on interactions between the environment and the animal's brain and behavior. Using the larval zebrafish as a model, we recorded natural optic flow associated with swimming trajectories in the animal's habitat with an omnidirectional camera mounted on a mechanical arm. An analysis of these flow fields suggests that lateral regions of the lower visual field are most informative about swimming speed. This pattern is consistent with the recent findings that zebrafish optomotor responses are preferentially driven by optic flow in the lateral lower visual field, which we extend with behavioral results from a high-resolution spherical arena. Spatial biases in optic-flow sampling are likely pervasive because they are an effective strategy for determining self-motion in noisy natural environments.

Optic flow in the natural habitats of zebrafish supports spatial biases in visual self-motion estimation

Animals benefit from knowing if and how they are moving. Across the animal kingdom, sensory information in the form of optic flow over the visual field is used to estimate self-motion. However, different species exhibit strong spatial biases in how they use optic flow. Here, we show computationally that noisy natural environments favor visual systems that extract spatially biased samples of optic flow when estimating self-motion. The performance associated with these biases, however, depends on interactions between the environment and the animal's brain and behavior. Using the larval zebrafish as a model, we recorded natural optic flow associated with swimming trajectories in the animal's habitat with an omnidirectional camera mounted on a mechanical arm. An analysis of these flow fields suggests that lateral regions of the lower visual field are most informative about swimming speed. This pattern is consistent with the recent findings that zebrafish optomotor responses are preferentially driven by optic flow in the lateral lower visual field, which we extend with behavioral results from a high-resolution spherical arena. Spatial biases in optic-flow sampling are likely pervasive because they are an effective strategy for determining self-motion in noisy natural environments.

Optic flow in the natural habitats of zebrafish supports spatial biases in visual self-motion estimation

Animals benefit from knowing if and how they are moving. Across the animal kingdom, sensory information in the form of optic flow over the visual field is used to estimate self-motion. However, different species exhibit strong spatial biases in how they use optic flow. Here, we show computationally that noisy natural environments favor visual systems that extract spatially biased samples of optic flow when estimating self-motion. The performance associated with these biases, however, depends on interactions between the environment and the animal's brain and behavior. Using the larval zebrafish as a model, we recorded natural optic flow associated with swimming trajectories in the animal's habitat with an omnidirectional camera mounted on a mechanical arm. An analysis of these flow fields suggests that lateral regions of the lower visual field are most informative about swimming speed. This pattern is consistent with the recent findings that zebrafish optomotor responses are preferentially driven by optic flow in the lateral lower visual field, which we extend with behavioral results from a high-resolution spherical arena. Spatial biases in optic-flow sampling are likely pervasive because they are an effective strategy for determining self-motion in noisy natural environments.

The Developmental Progression of Eight Opsin Spectral Signals Recorded from the Zebrafish Retinal Cone Layer Is Altered by the Timing and Cell Type Expression of Thyroxin Receptor β2 (trβ2) Gain-Of-Function Transgenes

Zebrafish retinal cone signals shift in spectral shape through larval, juvenile, and adult development as expression patterns of eight cone-opsin genes change. An algorithm extracting signal amplitudes for the component cone spectral types is developed and tested on two thyroxin receptor β2 (trβ2) gain-of-function lines crx:mYFP-2A-trβ2 and gnat2:mYFP-2A-trβ2, allowing correlation between opsin signaling and opsin immunoreactivity in lines with different developmental timing and cell-type expression of this red-opsin-promoting transgene. Both adult transgenics became complete, or nearly complete, "red-cone dichromats," with disproportionately large long-wavelength-sensitive (LWS)1 opsin amplitudes as compared with controls, where LWS1 and LWS2 amplitudes were about equal, and significant signals from SWS1, SWS2, and Rh2 opsins were detected. But in transgenic larvae and juveniles of both lines it was LWS2 amplitudes that increased, with LWS1 cone signals rarely encountered. In gnat2:mYFP-2A-trβ2 embryos at 5 d postfertilization (dpf), red-opsin immunoreactive cone density doubled, but red-opsin amplitudes (LWS2) increased <10%, and green-opsin, blue-opsin, and UV-opsin signals were unchanged, despite co-expressed red opsins, and the finding that an sws1 UV-opsin reporter gene was shut down by the gnat2:mYFP-2A-trβ2 transgene. By contrast both LWS2 red-cone amplitudes and the density of red-cone immunoreactivity more than doubled in 5-dpf crx:mYFP-2A-trβ2 embryos, while UV-cone amplitudes were reduced 90%. Embryonic cones with trβ2 gain-of-function transgenes were morphologically distinct from control red, blue or UV cones, with wider inner segments and shorter axons than red cones, suggesting cone spectral specification, opsin immunoreactivity and shape are influenced by the abundance and developmental timing of trβ2 expression.

Azole-Induced Color Vision Deficiency Associated with Thyroid Hormone Signaling: An Integrated In Vivo, In Vitro, and In Silico Study

Azoles that are used in pesticides, pharmaceuticals, and personal care products can have toxic effects on fish. However, there is no information regarding azole-induced visual disorder associated with thyroid disruption. We evaluated changes in retinal morphology, optokinetic response, transcript abundance of the genes involved in color perception and hypothalamic-pituitary-thyroid (HPT) axis, and thyroid hormone (TH) levels in zebrafish larvae exposed to common azoles, such as climbazole (CBZ, 0.1 and 10 μg/L) and triadimefon (TDF, 50 and 500 μg/L), at environmentally relevant and predicted worst-case environmental concentrations. Subsequently, the effect of azoles on TH-dependent GH3 cell proliferation and thyroid receptor (TR)-regulated transcriptional activity, as well as the in silico binding affinity between azoles and TR isoforms, was investigated. Azole exposure decreased cell densities of the ganglion cell layer, inner nuclear layer, and photoreceptor layer. Zebrafish larvae exposed to environmentally relevant concentrations of CBZ and TDF showed a decrease in optokinetic response to green-white and red-white stripes but not blue-white stripes, consistent with disturbance in the corresponding opsin gene expression. Azole exposure also reduced triiodothyronine levels and concomitantly increased HPT-related gene expression. Molecular docking analysis combined with in vitro TR-mediated transactivation and dual-luciferase reporter assays demonstrated that CBZ and TDF exhibited TR antagonism. These results are comparable to those obtained from a known TR antagonist, namely, TR antagonist 1, as a positive control. Therefore, damage to specific color perception by azoles appears to result from lowered TH signaling, indicating the potential threat of environmental TH disruptors to the visual function of fish.

An EvoDevo Study of Salmonid Visual Opsin Dynamics and Photopigment Spectral Sensitivity

Salmonids are ideal models as many species follow a distinct developmental program from demersal eggs and a large yolk sac to hatching at an advanced developmental stage. Further, these economically important teleosts inhabit both marine- and freshwaters and experience diverse light environments during their life histories. At a genome level, salmonids have undergone a salmonid-specific fourth whole genome duplication event (Ss4R) compared to other teleosts that are already more genetically diverse compared to many non-teleost vertebrates. Thus, salmonids display phenotypically plastic visual systems that appear to be closely related to their anadromous migration patterns. This is most likely due to a complex interplay between their larger, more gene-rich genomes and broad spectrally enriched habitats; however, the molecular basis and functional consequences for such diversity is not fully understood. This study used advances in genome sequencing to identify the repertoire and genome organization of visual opsin genes (those primarily expressed in retinal photoreceptors) from six different salmonids [Atlantic salmon (Salmo salar), brown trout (Salmo trutta), Chinook salmon (Oncorhynchus tshawytcha), coho salmon (Oncorhynchus kisutch), rainbow trout (Oncorhynchus mykiss), and sockeye salmon (Oncorhynchus nerka)] compared to the northern pike (Esox lucius), a closely related non-salmonid species. Results identified multiple orthologues for all five visual opsin classes, except for presence of a single short-wavelength-sensitive-2 opsin gene. Several visual opsin genes were not retained after the Ss4R duplication event, which is consistent with the concept of salmonid rediploidization. Developmentally, transcriptomic analyzes of Atlantic salmon revealed differential expression within each opsin class, with two of the long-wavelength-sensitive opsins not being expressed before first feeding. Also, early opsin expression in the retina was located centrally, expanding dorsally and ventrally as eye development progressed, with rod opsin being the dominant visual opsin post-hatching. Modeling by spectral tuning analysis and atomistic molecular simulation, predicted the greatest variation in the spectral peak of absorbance to be within the Rh2 class, with a ∼40 nm difference in λ max values between the four medium-wavelength-sensitive photopigments. Overall, it appears that opsin duplication and expression, and their respective spectral tuning profiles, evolved to maximize specialist color vision throughout an anadromous lifecycle, with some visual opsin genes being lost to tailor marine-based vision.

Correlated Expression of the Opsin Retrogene LWS-R and its Host Gene in Two Poeciliid Fishes

The important role of retrogenes in genome evolution and species differentiation is becoming increasingly accepted. One synapomorphy among cyprinodontoid fish is a retrotransposed version of a long-wavelength sensitive (LWS) opsin gene, LWS-R, within an intron of the gephyrin (GPHN) gene. These two genes display opposing orientations. It had been speculated that LWS-R hijacks the cis-regulatory elements of GPHN for transcription, but whether their expression is correlated had remained unclear. Here, in silico predictions identified putative promoters upstream of the translation start site of LWS-R, indicating that its transcription is driven by its own promoter rather than by the GPHN promoter. However, consistent expression ratios of LWS-R:GPHN in the eyeball and brain of fishes indicate that the respective gene transcriptions are correlated. Co-expression is potentially modulated by histone exchange during GPHN transcription. Two isoforms were detected in this study, i.e., intron-free and intron-retaining. Intron-free LWS-R was only expressed in the eyeball of fishes, whereas intron-retaining LWS-R occurred in both eyeball and brain. Expression of vision-associated LWS-R beyond the eyeball supports that it is co-expressed with more ubiquitous GPHN.

Bis (2-ethylhexyl)-2,3,4,5-tetrabromophthalate showed poor penetrability but increased the permeability of blood brain barrier: Evidences from in vitro and in vivo studies

Bis(2ethylhexyl)-2,3,4,5-tetrabromophthalate (TBPH), a replacement for restricted flame retardants, has become ubiquitous in the environment. To reveal the neurotoxicity and underlying mechanism of TBPH, we first evaluated its penetrability through the blood-brain barrier (BBB) using hCMEC/D3 cells as in vitro model, and found TBPH had poor penetrability through BBB with a maximum Papp of 14.8 × 10^-6 cms^-1. Further study using transgenetic zebrafish (Tg flk1: EGFP) as in vivo model confirmed that TBPH could affect the BBB permeability, probably via affecting the transcription of genes encoding tight junction proteins. Finally, wild type zebrafish embryos/larvae were exposed to TBPH to evaluate the neurotoxicity. The neurodevelopment, neurotransmitters and locomotor activity of zebrafish larvae did not changed, which may be because TBPH can hardly cross the BBB to pose direct exposure to the central nervous system. However, the transcription of opsins genes and visual response to light stimulation in zebrafish larvae were inhibited, pointing to additional mechanism that may cause visual impairment indirectly. Above all, these results can help further understand the neurotoxicity and underlying mechanism by TBPH, and also pointed out potential risk of this chemical to aquatic organisms.

Analysis of opsin gene family of Crimson snapper (Lutjanus erythropterus)

Opsin is a fellow of the G protein-coupled receptors (GPCRs) superfamily. It can be divided into visual and non-visual opsin according to whether it is directly involved in visual imaging. Opsin plays an important role in visual image formation and the regulation of non-image forming functions such as circadian entrainment in the growth, development and evolution of fish. Crimson snapper belongs to Perciforme mainly found in the Indo-West Pacific and the South China Sea. It is one of the most influential economic fishes in the South China Sea. In order to study the existence and expression of opsin gene in Crimson snapper, we sequenced the genome and tissue sample transcriptome of Crimson snapper. In this study, 32 opsin genes were identified from the genome of Crimson snapper. The length of these genes ranged from 1061 bp to 86203 bp and were distributed on 15 different chromosomes. The analysis of opsin gene family of Crimson snapper showed that the sws2 had two extra copies as compared with that of Zebrafish. Domain and motif analysis revealed that all the 32 opsin genes have seven-(pass)-transmembrane domain receptors (7TM receptors) each, and the opsin family contained 10 common motifs. The expression level of opsin gene, confirmed by RT-qPCR, was analyzed by using nine tissues transcriptome databases of Crimson snapper. The results showed that almost all opsin genes were highly expressed in the retina and brain, except opn7a and opn7b which were expressed in intestine and red skin, and almost no expression in other tissues. Our results provide a comprehensive basic knowledge for the opsin gene family of Crimson snapper, which has significance for the study of the function of opsin in Lutjanidaes.

Molecular mechanisms of developmental toxicity induced by BBP in zebrafish embryos

Butylbenzyl phthalate (BBP) has been shown to negatively affect the development of zebrafish embryos, however, its underlying mechanisms remain unclear. Therefore, our study aims to reveal the molecular mechanisms of developmental toxicity on zebrafish embryos. Zebrafish embryos were exposed to BBP (0, 0.6, and 1.2 mg/L) from 4 to 72 h post-fertilization (hpf). The adverse effects on zebrafish embryos were evaluated and the transcriptional profiles of zebrafish embryos were analyzed at 72hpf. Exposure to BBP decreased hatching and survival rates and induced obvious morphology abnormalities in zebrafish embryos. The activities of superoxide dismutase (SOD), catalase (CAT), and glutathione peroxidase (GSH-Px) significantly decrease and the malondialdehyde (MDA) content significantly increased with 1.2 mg/L BBP exposure. Global transcriptome profiling analysis demonstrated that 578 and 1257 genes were differentially expressed in zebrafish embryos in the 0.6 and 1.2 mg/L groups, respectively. Gene Ontology (GO) term enrichment analysis demonstrated that DEGs are related to many aspects of cell composition, biological processes, and molecular functions. Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis demonstrated that 13 and 22 pathways were significantly enriched in the 0.6 and 1.2 mg/L groups, respectively. DEGs were primarily concentrated in the metabolism of the 0.6 mg/L group and in the organismal systems and particularly affected vision and digestion in the 1.2 mg/L group. Our results contribute to a better understanding of the underlying mechanisms of developmental toxicity induced by phthalates.

Use of Zebrafish Models to Boost Research in Rare Genetic Diseases

Rare genetic diseases are a group of pathologies with often unmet clinical needs. Even if rare by a single genetic disease (from 1/2000 to 1/more than 1,000,000), the total number of patients concerned account for approximatively 400 million peoples worldwide. Finding treatments remains challenging due to the complexity of these diseases, the small number of patients and the challenge in conducting clinical trials. Therefore, innovative preclinical research strategies are required. The zebrafish has emerged as a powerful animal model for investigating rare diseases. Zebrafish combines conserved vertebrate characteristics with high rate of breeding, limited housing requirements and low costs. More than 84% of human genes responsible for diseases present an orthologue, suggesting that the majority of genetic diseases could be modelized in zebrafish. In this review, we emphasize the unique advantages of zebrafish models over other in vivo models, particularly underlining the high throughput phenotypic capacity for therapeutic screening. We briefly introduce how the generation of zebrafish transgenic lines by gene-modulating technologies can be used to model rare genetic diseases. Then, we describe how zebrafish could be phenotyped using state-of-the-art technologies. Two prototypic examples of rare diseases illustrate how zebrafish models could play a critical role in deciphering the underlying mechanisms of rare genetic diseases and their use to identify innovative therapeutic solutions.

Whole-genome resequencing of three Coilia nasus population reveals genetic variations in genes related to immune, vision, migration, and osmoregulation

Background

Coilia nasus is an important anadromous fish, widely distributed in China, Japan, and Korea. Based on morphological and ecological researches of C. nasus, two ecotypes were identified. One is the anadromous population (AP). The sexually mature fish run thousands of kilometers from marine to river for spawning. Another one is the resident population which cannot migrate. Based on their different habitats, they were classified into landlocked population (LP) and sea population (SP) which were resident in the freshwater lake and marine during the entire lifetime, respectively. However, they have never been systematically studied. Moreover, C. nasus is declining sharply due to overfishing and pollution recently. Therefore, further understandings of C. nasus populations are needed for germplasm protection.

Results

Whole-genome resequencing of AP, LP, and SP were performed to enrich the understanding of different populations of C. nasus. At the genome level, 3,176,204, 3,307,069, and 3,207,906 single nucleotide polymorphisms (SNPs) and 1,892,068, 2,002,912, and 1,922,168 insertion/deletion polymorphisms (InDels) were generated in AP, LP, and SP, respectively. Selective sweeping analysis showed that 1022 genes were selected in AP vs LP; 983 genes were selected in LP vs SP; 116 genes were selected in AP vs SP. Among them, selected genes related to immune, vision, migration, and osmoregulation were identified. Furthermore, their expression profiles were detected by quantitative real-time PCR. Expression levels of selected genes related to immune, and vision in LP were significantly lower than AP and SP. Selected genes related to migration in AP were expressed significantly more highly than LP. Expression levels of selected genes related to osmoregulation were also detected. The expression of NKAα and NKCC1 in LP were significantly lower than SP, while expression of NCC, SLC4A4, NHE3, and V-ATPase in LP was significantly higher than SP.

Conclusions

Combined to life history of C. nasus populations, our results revealed that the molecular mechanisms of their differences of immune, vision, migration, and osmoregulation. Our findings will provide a further understanding of different populations of C. nasus and will be beneficial for wild C. nasus protection.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-08182-0.

Ancestral circuits for vertebrate color vision emerge at the first retinal synapse

For color vision, retinal circuits separate information about intensity and wavelength. In vertebrates that use the full complement of four “ancestral” cone types, the nature and implementation of this computation remain poorly understood. Here, we establish the complete circuit architecture of outer retinal circuits underlying color processing in larval zebrafish. We find that the synaptic outputs of red and green cones efficiently rotate the encoding of natural daylight in a principal components analysis–like manner to yield primary achromatic and spectrally opponent axes, respectively. Blue cones are tuned to capture most remaining variance when opposed to green cones, while UV cone present a UV achromatic axis for prey capture. We note that fruitflies use essentially the same strategy. Therefore, rotating color space into primary achromatic and chromatic axes at the eye’s first synapse may thus be a fundamental principle of color vision when using more than two spectrally well-separated photoreceptor types.

Foxq2 determines blue cone identity in zebrafish

Most vertebrate lineages retain a tetrachromatic visual system, which is supported by a functional combination of spectrally distinct multiple cone photoreceptors, ultraviolet (UV), blue, green, and red cones. The blue cone identity is ensured by selective expression of blue (sws2) opsin, and the mechanism is poorly understood because sws2 gene has been lost in mammalian species such as mouse, whose visual system has been extensively studied. Here, we pursued loss-of-function studies on transcription factors expressed predominantly in zebrafish cone photoreceptors and identified Foxq2 as a blue cone–specific factor driving sws2 gene expression. Foxq2 has dual functions acting as an activator of sws2 transcription and as a suppressor of UV (sws1) opsin transcription in blue cones. A wide range of vertebrate species retain both foxq2 and sws2 genes. We propose that Foxq2-dependent sws2 expression is a prevalent regulatory mechanism that was acquired at the early stage of vertebrate evolution.

Ganglion cells in larval zebrafish retina integrate inputs from multiple cone types

We recently showed the presence of seven physiological cone opsins-R1 (575 nm), R2 (556 nm), G1 (460 nm), G3 (480 nm), B1 (415 nm), B2 (440 nm), and UV (358 nm)-in electroretinogram (ERG) recordings of larval zebrafish (Danio rerio) retina. Larval ganglion cells (GCs) are generally thought to integrate only four cone opsin signals (red, green, blue, and UV). We address the question as to whether they may integrate seven cone spectral signals. Here we examined the 127 possible combinations of seven cone signals to find the optimal representation, as based on impulse discharge data sets from GC axons in the larval optic nerve. We recorded four varieties of light-response waveform, sustained-ON, transient-ON, ON-OFF, and OFF, based on the time course of mean discharge rates to all stimulus wavelengths combined. Modeling of GC responses revealed that each received 1-6 cone opsin signals, with a mean of 3.8 ± 1.3 cone signals/GC. Most onset or offset responses were opponent (ON, 80%; OFF, 100%). The most common cone signals were UV (93%), R2 (50%), G3 (55%), and G1 (60%). Seventy-three percent of cone opsin signals were excitatory, and 27% were inhibitory. UV signals favored excitation, whereas G3 and B2 signals favored inhibition. R1/R2, G1/G3, and B1/B2 opsin signals were selectively associated along a nonsynergistic/opponent axis. Overall, these results suggest that larval zebrafish GC spectral responses are complex and use inputs from the seven expressed opsins.NEW & NOTEWORTHY Ganglion cells in larval zebrafish retina have complex spectral responses driven by seven different cone opsin types. UV cone inputs are significant and excitatory to ganglion cells, whereas green and blue cone inputs favor inhibition. Most dramatic are the pentachromatic cells. These responses were identified at 5-6 days after fertilization, reflecting an impressive level of color processing not seen in older fish or mammals.

Partitioning of gene expression among zebrafish photoreceptor subtypes

Vertebrate photoreceptors are categorized into two broad classes, rods and cones, responsible for dim- and bright-light vision, respectively. While many molecular features that distinguish rods and cones are known, gene expression differences among cone subtypes remain poorly understood. Teleost fishes are renowned for the diversity of their photoreceptor systems. Here, we used single-cell RNA-seq to profile adult photoreceptors in zebrafish, a teleost. We found that in addition to the four canonical zebrafish cone types, there exist subpopulations of green and red cones (previously shown to be located in the ventral retina) that express red-shifted opsin paralogs (opn1mw4 or opn1lw1) as well as a unique combination of cone phototransduction genes. Furthermore, the expression of many paralogous phototransduction genes is partitioned among cone subtypes, analogous to the partitioning of the phototransduction paralogs between rods and cones seen across vertebrates. The partitioned cone-gene pairs arose via the teleost-specific whole-genome duplication or later clade-specific gene duplications. We also discovered that cone subtypes express distinct transcriptional regulators, including many factors not previously implicated in photoreceptor development or differentiation. Overall, our work suggests that partitioning of paralogous gene expression via the action of differentially expressed transcriptional regulators enables diversification of cone subtypes in teleosts.

Molecular evolution of ultraviolet visual opsins and spectral tuning of photoreceptors in anemonefishes (Amphiprioninae)

Many animals including birds, reptiles, insects, and teleost fishes can see ultraviolet (UV) light (shorter than 400 nm), which has functional importance for foraging and communication. For coral reef fishes, shallow reef environments transmit a broad spectrum of light, rich in UV, driving the evolution of diverse spectral sensitivities. However, the identities and sites of the specific visual genes that underly vision in reef fishes remain elusive and are useful in determining how evolution has tuned vision to suit life on the reef. We investigated the visual systems of 11 anemonefish (Amphiprioninae) species, specifically probing for the molecular pathways that facilitate UV-sensitivity. Searching the genomes of anemonefishes, we identified a total of eight functional opsin genes from all five vertebrate visual opsin subfamilies. We found rare instances of teleost UV-sensitive SWS1 opsin gene duplications that produced two functionally coding paralogs (SWS1α and SWS1β) and a pseudogene. We also found separate green sensitive RH2A opsin gene duplicates not yet reported in the family Pomacentridae. Transcriptome analysis revealed false clown anemonefish (Amphiprion ocellaris) expressed one rod opsin (RH1) and six cone opsins (SWS1β, SWS2B, RH2B, RH2A-1, RH2A-2, LWS) in the retina. Fluorescent in situ hybridization highlighted the (co-)expression of SWS1β with SWS2B in single cones, and either RH2B, RH2A, or RH2A together with LWS in different members of double cone photoreceptors (two single cones fused together). Our study provides the first in-depth characterization of visual opsin genes found in anemonefishes and provides a useful basis for the further study of UV-vision in reef fishes.

Behavioural red-light sensitivity in fish according to the optomotor response

Various procedures have been adopted to investigate spectral sensitivity of animals, e.g. absorption spectra of visual pigments, electroretinography, optokinetic response, optomotor response (OMR) and phototaxis. The use of these techniques has led to various conclusions about animal vision. However, visual sensitivity should be evaluated consistently for a reliable comparison. In this study, we retrieved behavioural data of several fish species using a single OMR procedure and compared their sensitivities to near-infrared light. Besides cavefish that lack eyes, some species were not appropriate for the OMR test because they either stayed still or changed swimming direction frequently. Eight of 13 fish species tested were OMR positive. Detailed analyses using medaka, goldfish, zebrafish, guppy, stickleback and cichlid revealed that all the fish were sensitive to light at a wavelength greater than or equal to 750 nm, where the threshold wavelengths varied from 750 to 880 nm. Fish opsin repertoire affected the perception of red light. By contrast, the copy number of long-wavelength-sensitive (LWS) genes did not necessarily improve red-light sensitivity. While the duplication of LWS and other cone opsin genes that has occurred extensively during fish evolution might not aid increasing spectral sensitivity, it may provide some other advantageous ophthalmic function, such as enhanced spectral discrimination.

Orthologous Divergence and Paralogous Anticonvergence in Molecular Evolution of Triplicated Green Opsin Genes in Medaka Fish, Genus Oryzias

Gene duplication of green (RH2) opsin genes and their spectral differentiation are well documented in many teleost fish. However, their evolutionary divergence or conservation patterns among phylogenetically close but ecologically diverse species is not well explored. Medaka fish (genus Oryzias) are broadly distributed in fresh and brackish waters of Asia, with many species being laboratory-housed and feasible for genetic studies. We previously showed that a Japan strain (HNI) of medaka (Oryzias latipes) possessed three RH2 opsin genes (RH2-A, RH2-B, and RH2-C) encoding spectrally divergent photopigments. Here, we examined the three RH2 opsin genes from six Oryzias species representing three species groups: the latipes, the celebensis, and the javanicus. Photopigment reconstitution revealed that the peak absorption spectra (λ_max) of RH2-A were divergent among the species (447–469 nm), whereas those of RH2-B and RH2-C were conservative (516–519 and 486–493 nm, respectively). For the RH2-A opsins, the largest spectral shift was detected in the phylogenetic branch leading to the latipes group. A single amino acid replacement T94C explained most of the spectral shift. For RH2-B and -C opsins, we detected tracts of gene conversion between the two genes homogenizing them. Nevertheless, several amino acid differences were maintained. We showed that the spectral difference between the two opsins was attributed to largely the E/Q amino acid difference at the site 122 and to several sites with individually small spectral effects. These results depict dynamism of spectral divergence of orthologous and paralogous green opsin genes in phylogenetically close but ecologically diverse species exemplified by medaka.

A Chromosome-Level Genome Assembly of the Mandarin Fish (Siniperca chuatsi)

The mandarin fish, Siniperca chuatsi, is an economically important perciform species with widespread aquaculture practices in China. Its special feeding habit, acceptance of only live prey fishes, contributes to its delicious meat. However, little is currently known about related genetic mechanisms. Here, we performed whole-genome sequencing and assembled a 758.78 Mb genome assembly of the mandarin fish, with the scaffold and contig N50 values reaching 2.64 Mb and 46.11 Kb, respectively. Approximately 92.8% of the scaffolds were ordered onto 24 chromosomes (Chrs) with the assistance of a previously established genetic linkage map. The chromosome-level genome contained 19,904 protein-coding genes, of which 19,059 (95.75%) genes were functionally annotated. The special feeding behavior of mandarin fish could be attributable to the interaction of a variety of sense organs (such as vision, smell, and endocrine organs). Through comparative genomics analysis, some interesting results were found. For example, olfactory receptor (OR) genes (especially the beta and delta types) underwent a significant expansion, and endocrinology/vision related npy, spexin, and opsin genes presented various functional mutations. These may contribute to the special feeding habit of the mandarin fish by strengthening the olfactory and visual systems. Meanwhile, previously identified sex-related genes and quantitative trait locis (QTLs) were localized on the Chr14 and Chr17, respectively. 155 toxin proteins were predicted from mandarin fish genome. In summary, the high-quality genome assembly of the mandarin fish provides novel insights into the feeding habit of live prey and offers a valuable genetic resource for the quality improvement of this freshwater fish.

Evolutionary ecology of the visual opsin gene sequence and its expression in turbot (Scophthalmus maximus)

Background

As flatfish, turbot undergo metamorphosis as part of their life cycle. In the larval stage, turbot live at the ocean surface, but after metamorphosis they move to deeper water and turn to benthic life. Thus, the light environment differs greatly between life stages. The visual system plays a great role in organic evolution, but reports of the relationship between the visual system and benthic life are rare. In this study, we reported the molecular and evolutionary analysis of opsin genes in turbot, and the heterochronic shifts in opsin expression during development.

Results

Our gene synteny analysis showed that subtype RH2C was not on the same gene cluster as the other four green-sensitive opsin genes (RH2) in turbot. It was translocated to chromosome 8 from chromosome 6. Based on branch-site test and spectral tuning sites analyses, E122Q and M207L substitutions in RH2C, which were found to be under positive selection, are closely related to the blue shift of optimum light sensitivities. And real-time PCR results indicated the dominant opsin gene shifted from red-sensitive (LWS) to RH2B1 during turbot development, which may lead to spectral sensitivity shifts to shorter wavelengths.

Conclusions

This is the first report that RH2C may be an important subtype of green opsin gene that was retained by turbot and possibly other flatfish species during evolution. Moreover, E122Q and M207L substitutions in RH2C may contribute to the survival of turbot in the bluish colored ocean. And heterochronic shifts in opsin expression may be an important strategy for turbot to adapt to benthic life.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12862-021-01837-2.

Circuit mechanisms for colour vision in zebrafish

The use of spectral information in natural light to inform behaviour is one of the oldest and most fundamental abilities of visual systems. It long-predates animals' venture onto the land, and even the appearance of image-forming eyes. Accordingly, circuits for colour vision evolved under the surface of ancient oceans for hundreds of millions of years. These aquatic beginnings fundamentally underpin, and likely constrain, the organisation of modern visual systems. In contrast to our detailed circuit level understanding from diverse terrestrial vertebrates, however, comparatively little is known about their aquatic counterparts. Here, I summarise some of what is known about neural circuits for colour vision in fish, the most species-diverse group of vertebrates. With a focus on zebrafish, I will explore how their computational strategies are linked to the statistics of natural light in the underwater world, and how their study might help us understand vision in general, including in our own eyes.

Microplastics decrease the toxicity of triphenyl phosphate (TPhP) in the marine medaka (Oryzias melastigma) larvae

Plastics have been recognized as a serious threat to the environment. Besides their own toxicity, microplastics can interact with other environmental pollutants, acting as carriers and potentially modulating their toxicity. In this study, the toxicity of polystyrene (PS) microplastic fragments (plain PS; carboxylated PS, PS-COOH and aminated PS, PS-NH₂) and triphenyl phosphate (TPhP) (an emerging organophosphate flame retardant) at the environmentally relevant concentrations to the marine medaka (Oryzias melastigma) larvae was investigated. Larvae were exposed to 20 μg/L of microplastic fragments or 20 and 100 μg/L of TPhP or a combination of both for 7 days. The results showed that the three microplastics did not affect the larval locomotor activity. For TPhP, the larval moving duration and distance moved were significantly decreased by the TPhP exposure, with a maximum decrease of 43.5% and 59.4% respectively. Exposure to 100 μg/L TPhP respectively down-regulated the expression levels of sine oculis homeobox homologue 3 (six3) and short wavelength-sensitive type 2 (sws2) by 19.1% and 41.7%, suggesting that TPhP might disturb eye development and photoreception and consequently the low locomotor activity in the larvae. Interestingly, during the binary mixture exposure, the presence of PS, PS-COOH or PS-NH₂ reversed the low locomotor activity induced by 100 μg/L TPhP to the normal level. Relative to the larvae from the 100 μg/L TPhP group, the movement duration and distance moved were increased by approximately 60% and 100%, respectively, in the larvae from the TPhP + PS, TPhP + PS-COOH and TPhP + PS-NH₂ groups. However, the gene expression profiles were distinct among the fish from the TPhP + PS, TPhP + PS-COOH and TPhP + PS-NH₂ groups, implying different mechanisms underlying the reversal of the locomotor activity. The findings in this study challenge the general view that microplastics aggravate the toxicity of the adsorbed pollutants, and help better understand the environmental risk of microplastic pollution.

Expression dynamics of visual opsin genes in marbled sole Pseudopleuronectes yokohamae during metamorphosis from the larval to the juvenile stage

We clarified the properties of visual opsin genes in the marbled sole (Pseudopleuronectes yokohamae) by cDNA sequencing, quantification of the opsin gene expression from the larval to the juvenile stage, and measurement of the maximum absorption spectra (λ_max) using photopigment reconstitution. In the marbled sole eye, at least eight visual opsin genes, lws, rh2-a, rh2-b, rh2-c, sws2a, sws2b, sws1, and rh1, were expressed. Quantitative RT-PCR analysis revealed that the expression of opsin genes increased (lws, rh2-c, sws2a, and rh1) or decreased (rh2-a, rh2-b, sws2b, and sws1) from the larval to the juvenile stage. Notably, rh2-a expression was observed only in pre- to mid-metamorphic stage larvae and disappeared after metamorphosis. Thus, pre-metamorphism-specific expression of rh2-a in the marbled sole suggests that its function is restricted to the developmental stage. The reconstituted RH2-A opsin λ_max was 470 nm, which is typical of acanthopterygian species. These results strongly suggest that mid-wavelength-sensitive rh2-a expression was diminished drastically in the marbled sole, probably resulting in a shift of spectral sensitivity during its metamorphosis from the larval to the juvenile stage.

Neurodegeneration, Neuroprotection and Regeneration in the Zebrafish Retina

Neurodegenerative retinal diseases, such as glaucoma and diabetic retinopathy, involve a gradual loss of neurons in the retina as the disease progresses. Central nervous system neurons are not able to regenerate in mammals, therefore, an often sought after course of treatment for neuronal loss follows a neuroprotective or regenerative strategy. Neuroprotection is the process of preserving the structure and function of the neurons that have survived a harmful insult; while regenerative approaches aim to replace or rewire the neurons and synaptic connections that were lost, or induce regrowth of damaged axons or dendrites. In order to test the neuroprotective effectiveness or the regenerative capacity of a particular agent, a robust experimental model of retinal neuronal damage is essential. Zebrafish are being used more often in this type of study because their eye structure and development is well-conserved between zebrafish and mammals. Zebrafish are robust genetic tools and are relatively inexpensive to maintain. The large array of functional and behavioral tests available in zebrafish makes them an attractive model for neuroprotection studies. Some common insults used to model retinal disease and study neuroprotection in zebrafish include intense light, chemical toxicity and mechanical damage. This review covers the existing retinal neuroprotection and regeneration literature in the zebrafish and highlights their potential for future studies.

Thyroid hormone signaling specifies cone photoreceptor subtypes during eye development: Insights from model organisms and human stem cell-derived retinal organoids

Cones are the color-detecting photoreceptors of the vertebrate eye. Cones are specialized into subtypes whose functions are determined by the expression of color-sensitive opsin proteins. Organisms differ greatly in the number and patterning of cone subtypes. Despite these differences, thyroid hormone is an important regulator of opsin expression in most vertebrates. In this chapter, we outline how the timing of thyroid hormone signaling controls cone subtype fates during retinal development. We first examine our current understanding of cone subtype specification in model organisms and then describe advances in human stem cell-derived organoid technology that identified mechanisms controlling development of the human retina.

Biochemistry and physiology of zebrafish photoreceptors

All vertebrates share a canonical retina with light-sensitive photoreceptors in the outer retina. These photoreceptors are of two kinds: rods and cones, adapted to low and bright light conditions, respectively. They both show a peculiar morphology, with long outer segments, comprised of ordered stacks of disc-shaped membranes. These discs host numerous proteins, many of which contribute to the visual transduction cascade. This pathway converts the light stimulus into a biological signal, ultimately modulating synaptic transmission. Recently, the zebrafish (Danio rerio) has gained popularity for studying the function of vertebrate photoreceptors. In this review, we introduce this model system and its contribution to our understanding of photoreception with a focus on the cone visual transduction cascade.

Identification of molecular toxicity pathways across early life-stages of zebrafish exposed to PCB126 using a whole transcriptomics approach

Although withdrawn from the market in the 1980s, polychlorinated biphenyls (PCBs) are still found ubiquitously in the aquatic environment and pose a serious risk to biota due to their teratogenic potential. In fish, early life-stages are often considered most sensitive with regard to their exposure to PCBs and other dioxin-like compounds. However, little is known about the molecular drivers of the frequently observed teratogenic effects. Therefore, the aims of our study were to: (1) characterize the baseline transcriptome profiles at different embryonic life-stages in zebrafish (Danio rerio); and (2) to identify the molecular response to PCB exposure and life-stage specific-effects of the chemical on associated processes. For both objectives, embryos were sampled at 12, 48, and 96 h post-fertilization (hpf) and subjected to Illumina sequence-by-synthesis and RNAseq analysis. Results revealed that with increasing age more genes and related pathways were upregulated both in terms of number and magnitude. Yet, other transcripts followed an opposite pattern with greater transcript abundance at the earlier time points. Additionally, embryos were exposed to PCB126, a potent agonist of the aryl hydrocarbon receptor (AHR). ClueGO network analysis revealed significant enrichment of genes associated with basic cell metabolism, communication, and homeostasis as well as eye development, muscle formation, and skeletal formation. We selected eight genes involved in the affected pathways for an in-depth characterization of their regulation throughout normal embryogenesis and after exposure to PCB126 by quantification of transcript abundances every 12 h until 118 hpf. Among these, fgf7 and c9 stood out because of their strong upregulation by PCB126 exposure at 48 and 96 hpf, respectively. Cyp2aa12 was upregulated from 84 hpf on. Fabp10ab, myhz1.1, col8a1a, sulf1, and opn1sw1 displayed specific regulation depending on the developmental stage. Overall, we demonstrate that (1) the developmental transcriptome of zebrafish is highly dynamic, and (2) dysregulation of gene expression by exposure to PCB126 was significant and in several cases not directly connected to AHR-signaling. Hence, this study improves the understanding of linkages between molecular events and apical outcomes that are of regulatory relevance.

Defect patterns on the curved surface of fish retinae suggest a mechanism of cone mosaic formation

The outer epithelial layer of zebrafish retinae contains a crystalline array of cone photoreceptors, called the cone mosaic. As this mosaic grows by mitotic addition of new photoreceptors at the rim of the hemispheric retina, topological defects, called "Y-Junctions", form to maintain approximately constant cell spacing. The generation of topological defects due to growth on a curved surface is a distinct feature of the cone mosaic not seen in other well-studied biological patterns like the R8 photoreceptor array in the Drosophila compound eye. Since defects can provide insight into cell-cell interactions responsible for pattern formation, here we characterize the arrangement of cones in individual Y-Junction cores as well as the spatial distribution of Y-junctions across entire retinae. We find that for individual Y-junctions, the distribution of cones near the core corresponds closely to structures observed in physical crystals. In addition, Y-Junctions are organized into lines, called grain boundaries, from the retinal center to the periphery. In physical crystals, regardless of the initial distribution of defects, defects can coalesce into grain boundaries via the mobility of individual particles. By imaging in live fish, we demonstrate that grain boundaries in the cone mosaic instead appear during initial mosaic formation, without requiring defect motion. Motivated by this observation, we show that a computational model of repulsive cell-cell interactions generates a mosaic with grain boundaries. In contrast to paradigmatic models of fate specification in mostly motionless cell packings, this finding emphasizes the role of cell motion, guided by cell-cell interactions during differentiation, in forming biological crystals. Such a route to the formation of regular patterns may be especially valuable in situations, like growth on a curved surface, where the resulting long-ranged, elastic, effective interactions between defects can help to group them into grain boundaries.

Daily mitochondrial dynamics in cone photoreceptors

Cone photoreceptors in the retina are exposed to intense daylight and have higher energy demands in darkness. Cones produce energy using a large cluster of mitochondria. Mitochondria are susceptible to oxidative damage, and healthy mitochondrial populations are maintained by regular turnover. Daily cycles of light exposure and energy consumption suggest that mitochondrial turnover is important for cone health. We investigated the three-dimensional (3D) ultrastructure and metabolic function of zebrafish cone mitochondria throughout the day. At night retinas undergo a mitochondrial biogenesis event, corresponding to an increase in the number of smaller, simpler mitochondria and increased metabolic activity in cones. In the daytime, endoplasmic reticula (ER) and autophagosomes associate more with mitochondria, and mitochondrial size distribution across the cluster changes. We also report dense material shared between cone mitochondria that is extruded from the cell at night, sometimes forming extracellular structures. Our findings reveal an elaborate set of daily changes to cone mitochondrial structure and function.

Short-wavelength-sensitive 2 (Sws2) visual photopigment models combined with atomistic molecular simulations to predict spectral peaks of absorbance

For many species, vision is one of the most important sensory modalities for mediating essential tasks that include navigation, predation and foraging, predator avoidance, and numerous social behaviors. The vertebrate visual process begins when photons of the light interact with rod and cone photoreceptors that are present in the neural retina. Vertebrate visual photopigments are housed within these photoreceptor cells and are sensitive to a wide range of wavelengths that peak within the light spectrum, the latter of which is a function of the type of chromophore used and how it interacts with specific amino acid residues found within the opsin protein sequence. Minor differences in the amino acid sequences of the opsins are known to lead to large differences in the spectral peak of absorbance (i.e. the λmax value). In our prior studies, we developed a new approach that combined homology modeling and molecular dynamics simulations to gather structural information associated with chromophore conformation, then used it to generate statistical models for the accurate prediction of λmax values for photopigments derived from Rh1 and Rh2 amino acid sequences. In the present study, we test our novel approach to predict the λmax of phylogenetically distant Sws2 cone opsins. To build a model that can predict the λmax using our approach presented in our prior studies, we selected a spectrally-diverse set of 11 teleost Sws2 photopigments for which both amino acid sequence information and experimentally measured λmax values are known. The final first-order regression model, consisting of three terms associated with chromophore conformation, was sufficient to predict the λmax of Sws2 photopigments with high accuracy. This study further highlights the breadth of our approach in reliably predicting λmax values of Sws2 cone photopigments, evolutionary-more distant from template bovine RH1, and provided mechanistic insights into the role of known spectral tuning sites.

Method Standardization for Conducting Innate Color Preference Studies in Different Zebrafish Strains

The zebrafish has a tetrachromatic vision that is able to distinguish ultraviolet (UV) and visible wavelengths. Recently, zebrafish color preferences have gained much attention because of the easy setup of the instrument and its usefulness to screen behavior-linked stimuli. However, several published papers dealing with zebrafish color preferences have contradicting results that underscore the importance of method standardization in this field. Different laboratories may report different results because of variations in light source, color intensity, and other parameters such as age, gender, container size, and strain of fish. In this study, we aim to standardize the color preference test in zebrafish by measuring light source position, light intensity, gender, age, animal size to space ratio, and animal strain. Our results showed that color preferences for zebrafish are affected by light position, age, strain, and social interaction of the fish, but not affected by fish gender. We validated that ethanol can significantly induce color preference alteration in zebrafish which may be related to anxiety and depression. We also explored the potential use of the optimized method to examine color preference ranking and index differences in various zebrafish strains and species, such as the tiger barb and glass catfish. In conclusion, zebrafish color preference screening is a powerful tool for high-throughput neuropharmacological applications and the standardized protocol established in this study provides a useful reference for the zebrafish research community.

Fovea-like Photoreceptor Specializations Underlie Single UV Cone Driven Prey-Capture Behavior in Zebrafish

In the eye, the function of same-type photoreceptors must be regionally adjusted to process a highly asymmetrical natural visual world. Here, we show that UV cones in the larval zebrafish area temporalis are specifically tuned for UV-bright prey capture in their upper frontal visual field, which may use the signal from a single cone at a time. For this, UV-photon detection probability is regionally boosted more than 10-fold. Next, in vivo two-photon imaging, transcriptomics, and computational modeling reveal that these cones use an elevated baseline of synaptic calcium to facilitate the encoding of bright objects, which in turn results from expressional tuning of phototransduction genes. Moreover, the light-driven synaptic calcium signal is regionally slowed by interactions with horizontal cells and later accentuated at the level of glutamate release driving retinal networks. These regional differences tally with variations between peripheral and foveal cones in primates and hint at a common mechanistic origin.

Exposure route affects the distribution and toxicity of polystyrene nanoplastics in zebrafish

The widespread use of polystyrene (PS) products in a myriad of consumer products has resulted in widespread contamination of PS nanoplastics (PSNPs) in aquatic ecosystems. Fish early life stages are exposed to nanoplastics dermally and via gills. Additional routes of exposure include oral via the ingestion of contaminated prey and maternal transfer. However, there is limited amount of work studying the impact of exposure route in the toxicokinetics and toxicodynamics of PSNPs. The objective of this study was to compare the effects of exposure routes (aqueous and microinjection) on the organ distribution and toxicity of PSNPs. We "mimicked" the maternal exposure of PSNPs to zebrafish by injecting a known concentration of fluorescent particles directly into 2-cell stage embryos. Endpoints were collected starting at 96 h post-fertilization until several weeks post-hatch to evaluate depuration. Although both exposure routes led to the accumulation of PSNPs in the yolk sac followed by brain, eyes, gut and swim bladder, the aqueous exposure caused higher PSNP concentrations in the brain and eyes and the injection exposure caused PSNP accumulation mainly in the trunk area. A waterborne exposure also reduced antioxidant gene expression; increased frequency of developmental abnormalities such as bent tails, jaw deformities and pericardial edema; and resulted in lower growth rates and hypoactivity. Overall, a waterborne exposure to PSNPs resulted in higher transfer to the brain and caused greater toxic effects to zebrafish compared to an injection exposure and highlights the key role of exposure routes in the uptake, localization and subsequent distribution of nanoparticles.

Visual pigment evolution in Characiformes: The dynamic interplay of teleost whole-genome duplication, surviving opsins and spectral tuning

Vision represents an excellent model for studying adaptation, given the genotype-to-phenotype map that has been characterized in a number of taxa. Fish possess a diverse range of visual sensitivities and adaptations to underwater light, making them an excellent group to study visual system evolution. In particular, some speciose but understudied lineages can provide a unique opportunity to better understand aspects of visual system evolution such as opsin gene duplication and neofunctionalization. In this study, we showcase the visual system evolution of neotropical Characiformes and the spectral tuning mechanisms they exhibit to modulate their visual sensitivities. Such mechanisms include gene duplications and losses, gene conversion, opsin amino acid sequence and expression variation, and A₁ /A₂ -chromophore shifts. The Characiforms we studied utilize three cone opsin classes (SWS2, RH2, LWS) and a rod opsin (RH1). However, the characiform's entire opsin gene repertoire is a product of dynamic evolution by opsin gene loss (SWS1, RH2) and duplication (LWS, RH1). The LWS- and RH1-duplicates originated from a teleost specific whole-genome duplication as well as characiform-specific duplication events. Both LWS-opsins exhibit gene conversion and, through substitutions in key tuning sites, one of the LWS-paralogues has acquired spectral sensitivity to green light. These sequence changes suggest reversion and parallel evolution of key tuning sites. Furthermore, characiforms' colour vision is based on the expression of both LWS-paralogues and SWS2. Finally, we found interspecific and intraspecific variation in A₁ /A₂ -chromophores proportions, correlating with the light environment. These multiple mechanisms may be a result of the diverse visual environments where Characiformes have evolved.

Origin and adaptation of green-sensitive (RH2) pigments in vertebrates

One of the critical times for the survival of animals is twilight where the most abundant visible lights are between 400 and 550 nanometres (nm). Green-sensitive RH2 pigments help nonmammalian vertebrate species to better discriminate wavelengths in this blue-green region. Here, evaluation of the wavelengths of maximal absorption (λmax s) of genetically engineered RH2 pigments representing 13 critical stages of vertebrate evolution revealed that the RH2 pigment of the most recent common ancestor of vertebrates had a λmax of 503 nm, while the 12 ancestral pigments exhibited an expanded range in λmax s between 474 and 524 nm, and present-day RH2 pigments have further expanded the range to ~ 450-530 nm. During vertebrate evolution, eight out of the 16 significant λmax shifts (or |Δλmax | ≥ 10 nm) of RH2 pigments identified were fully explained by the repeated mutations E122Q (twice), Q122E (thrice) and M207L (twice), and A292S (once). Our data indicated that the highly variable λmax s of teleost RH2 pigments arose from gene duplications followed by accelerated amino acid substitution.