Predicting peak spectral sensitivities of vertebrate cone visual pigments using atomistic molecular simulations

Deep-sea fish reveal alternative pathway for vertebrate visual development

Vertebrate vision is accomplished by two phenotypically distinct types of photoreceptors in the retina: the saturation-resistant cones for the detection of bright light and the highly sensitive rods for dim light conditions [1]. The current dogma is that, during development, all vertebrates initially feature a cone-dominated retina, and rods are added later [2, 3]. By studying the ontogeny of vision in three species of deep-sea fishes, we show that their larvae express cone-specific genes in photoreceptors with rod-like morphologies. Through development, these fishes either retain this rod-like cone retina (Maurolicus mucronatus) or switch to a retina with true rod photoreceptors with expression of rod-specific genes and transcription factors (Vinciguerria mabahiss and Benthosema pterotum). In contrast to the larvae of most marine fishes, which inhabit the bright upper layer of the open ocean, the larvae of deep-sea fishes occur deeper, exposing them to a dimmer light environment [4-7]. Spectral maxima predictions from molecular dynamics simulations and environmental light estimations suggest that using transmuted photoreceptors that combine the characteristics of both cones and rods maximises visual performance in these dimmer light conditions. Our findings provide molecular, morphological, and functional evidence for the evolution of an alternative developmental pathway for vertebrate vision.

Molecular Mechanisms Underlying the Spectral Shift in Zebrafish Cone Opsins

Visual pigments are essential for converting light into electrical signals during vision. Composed of an opsin protein and a retinal-based chromophore, pigments in vertebrate rods (Rh1) and cones (Rh2) have different spectral sensitivities, with distinct peak absorption wavelengths determined by the shape and composition of the chromophore binding pocket. Despite advances in understanding Rh1 pigments such as bovine rhodopsin, the molecular basis of spectral shifts in Rh2 cone opsins has been less studied, particularly the E122Q mutation, which accounts for about half of the observed spectral shift in these pigments. In this study, we employed molecular modeling and quantum mechanical techniques to investigate the molecular mechanisms behind the spectral difference in blue-shifted Rh2-1 (absorption peak = 467 nm, 122Q) and green-shifted Rh2-4 (absorption peak = 505 nm, 122E) zebrafish cone opsins. We modeled the pigments 3D structures based on their sequences and conducted all-atom molecular dynamics simulations totaling 2 microseconds. Distance analysis of the trajectories identified three key sites: E113, E181, and E122. The E122Q mutation, previously known, validates our findings, while E181 and E113 are newly identified contributors. Structural analysis revealed key features with differing values that explain the divergent spectral sensitivities of Rh2-1 and Rh2-4: 1) chromophore atom fluctuations and C5-C6 torsion angle, 2) binding pocket volume, 3) hydration patterns, and 4) E113-chromophore interaction stability. Quantum mechanics further confirms the critical role of residue E181 in Rh2-1 and E122 in Rh2-4 for their spectral behavior. Our study provides new insights into the molecular determinants of spectral shifts in cone opsins, and we anticipate that it will serve as a starting point for a broader understanding of the functional diversity of visual pigments.

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.

Discovering genotype-phenotype relationships with machine learning and the Visual Physiology Opsin Database (VPOD)

BACKGROUND

Predicting phenotypes from genetic variation is foundational for fields as diverse as bioengineering and global change biology, highlighting the importance of efficient methods to predict gene functions. Linking genetic changes to phenotypic changes has been a goal of decades of experimental work, especially for some model gene families, including light-sensitive opsin proteins. Opsins can be expressed in vitro to measure light absorption parameters, including λmax-the wavelength of maximum absorbance-which strongly affects organismal phenotypes like color vision. Despite extensive research on opsins, the data remain dispersed, uncompiled, and often challenging to access, thereby precluding systematic and comprehensive analyses of the intricate relationships between genotype and phenotype.

RESULTS

Here, we report a newly compiled database of all heterologously expressed opsin genes with λmax phenotypes that we call the Visual Physiology Opsin Database (VPOD). VPOD_1.0 contains 864 unique opsin genotypes and corresponding λmax phenotypes collected across all animals from 73 separate publications. We use VPOD data and deepBreaks to show regression-based machine learning (ML) models often reliably predict λmax, account for nonadditive effects of mutations on function, and identify functionally critical amino acid sites.

CONCLUSION

The ability to reliably predict functions from gene sequences alone using ML will allow robust exploration of molecular-evolutionary patterns governing phenotype, will inform functional and evolutionary connections to an organism's ecological niche, and may be used more broadly for de novo protein design. Together, our database, phenotype predictions, and model comparisons lay the groundwork for future research applicable to families of genes with quantifiable and comparable phenotypes.

Molecular advances to study the function, evolution and spectral tuning of arthropod visual opsins

Visual opsins of vertebrates and invertebrates diversified independently and converged to detect ultraviolet to long wavelengths (LW) of green or red light. In both groups, colour vision largely derives from opsin number, expression patterns and changes in amino acids interacting with the chromophore. Functional insights regarding invertebrate opsin evolution have lagged behind those for vertebrates because of the disparity in genomic resources and the lack of robust in vitro systems to characterize spectral sensitivities. Here, we review bioinformatic approaches to identify and model functional variation in opsins as well as recently developed assays to measure spectral phenotypes. In particular, we discuss how transgenic lines, cAMP-spectroscopy and sensitive heterologous expression platforms are starting to decouple genotype-phenotype relationships of LW opsins to complement the classical physiological-behavioural-phylogenetic toolbox of invertebrate visual sensory studies. We illustrate the use of one heterologous method by characterizing novel LW Gq opsins from 10 species, including diurnal and nocturnal Lepidoptera, a terrestrial dragonfly and an aquatic crustacean, expressing them in HEK293T cells, and showing that their maximum absorbance spectra (λmax) range from 518 to 611 nm. We discuss the advantages of molecular approaches for arthropods with complications such as restricted availability, lateral filters, specialized photochemistry and/or electrophysiological constraints. This article is part of the theme issue 'Understanding colour vision: molecular, physiological, neuronal and behavioural studies in arthropods'.

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.

Rhodopsins: An Excitingly Versatile Protein Species for Research, Development and Creative Engineering

The first member and eponym of the rhodopsin family was identified in the 1930s as the visual pigment of the rod photoreceptor cell in the animal retina. It was found to be a membrane protein, owing its photosensitivity to the presence of a covalently bound chromophoric group. This group, derived from vitamin A, was appropriately dubbed retinal. In the 1970s a microbial counterpart of this species was discovered in an archaeon, being a membrane protein also harbouring retinal as a chromophore, and named bacteriorhodopsin. Since their discovery a photogenic panorama unfolded, where up to date new members and subspecies with a variety of light-driven functionality have been added to this family. The animal branch, meanwhile categorized as type-2 rhodopsins, turned out to form a large subclass in the superfamily of G protein-coupled receptors and are essential to multiple elements of light-dependent animal sensory physiology. The microbial branch, the type-1 rhodopsins, largely function as light-driven ion pumps or channels, but also contain sensory-active and enzyme-sustaining subspecies. In this review we will follow the development of this exciting membrane protein panorama in a representative number of highlights and will present a prospect of their extraordinary future potential.

The Characterization of a Novel Virus Discovered in the Yeast Pichia membranifaciens

Mycoviruses are widely distributed across fungi, including the yeasts of the Saccharomycotina subphylum. This manuscript reports the first double-stranded RNA (dsRNA) virus isolated from Pichia membranifaciens. This novel virus has been named Pichia membranifaciens virus L-A (PmV-L-A) and is a member of the Totiviridae. PmV-L-A is 4579 bp in length, with RNA secondary structures similar to the packaging, replication, and frameshift signals of totiviruses that infect Saccharomycotina yeasts. PmV-L-A was found to be part of a monophyletic group within the I-A totiviruses, implying a shared ancestry between mycoviruses isolated from the Pichiaceae and Saccharomycetaceae yeasts. Energy-minimized AlphaFold2 molecular models of the PmV-L-A Gag protein revealed structural conservation with the Gag protein of Saccharomyces cerevisiae virus L-A (ScV-L-A). The predicted tertiary structure of the PmV-L-A Pol and other homologs provided a possible mechanism for totivirus RNA replication due to structural similarities with the RNA-dependent RNA polymerases of mammalian dsRNA viruses. Insights into the structure, function, and evolution of totiviruses gained from yeasts are essential because of their emerging role in animal disease and their parallels with mammalian viruses.

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.

The pesticide chlorpyrifos promotes obesity by inhibiting diet-induced thermogenesis in brown adipose tissue

Obesity results from a caloric imbalance between energy intake, absorption and expenditure. In both rodents and humans, diet-induced thermogenesis contributes to energy expenditure and involves the activation of brown adipose tissue (BAT). We hypothesize that environmental toxicants commonly used as food additives or pesticides might reduce BAT thermogenesis through suppression of uncoupling protein 1 (UCP1) and this may contribute to the development of obesity. Using a step-wise screening approach, we discover that the organophosphate insecticide chlorpyrifos suppresses UCP1 and mitochondrial respiration in BAT at concentrations as low as 1 pM. In mice housed at thermoneutrality and fed a high-fat diet, chlorpyrifos impairs BAT mitochondrial function and diet-induced thermogenesis, promoting greater obesity, non-alcoholic fatty liver disease (NAFLD) and insulin resistance. This is associated with reductions in cAMP; activation of p38MAPK and AMPK; protein kinases critical for maintaining UCP1 and mitophagy, respectively in BAT. These data indicate that the commonly used pesticide chlorpyrifos, suppresses diet-induced thermogenesis and the activation of BAT, suggesting its use may contribute to the obesity epidemic.

The Visual Opsin Gene Repertoires of Teleost Fishes: Evolution, Ecology, and Function

Visual opsin genes expressed in the rod and cone photoreceptor cells of the retina are core components of the visual sensory system of vertebrates. Here, we provide an overview of the dynamic evolution of visual opsin genes in the most species-rich group of vertebrates, teleost fishes. The examination of the rich genomic resources now available for this group reveals that fish genomes contain more copies of visual opsin genes than are present in the genomes of amphibians, reptiles, birds, and mammals. The expansion of opsin genes in fishes is due primarily to a combination of ancestral and lineage-specific gene duplications. Following their duplication, the visual opsin genes of fishes repeatedly diversified at the same key spectral-tuning sites, generating arrays of visual pigments sensitive from the ultraviolet to the red spectrum of the light. Species-specific opsin gene repertoires correlate strongly with underwater light habitats, ecology, and color-based sexual selection. Expected final online publication date for the Annual Review of Cell and Developmental Biology, Volume 37 is October 2021. Please see http://www.annualreviews.org/page/journal/pubdates for revised estimates.

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.

The Chemosensory Repertoire of the Eastern Diamondback Rattlesnake (Crotalus adamanteus) Reveals Complementary Genetics of Olfactory and Vomeronasal-Type Receptors

Pitviper sensory perception incorporates diverse stimuli through the integration of trichromatic color vision, bifocal heat-sensing, and dual-system chemoperception. Chemoperception, or olfaction, is mediated by chemoreceptors in the olfactory bulb and the vomeronasal organ, but the true genomic complexity of the gene families and their relative contributions is unknown. A full genomic accounting of pitviper chemoperception directly complements our current understanding of their venoms by generating a more complete polyphenic representation of their predatory arsenal. To characterize the genetic repertoire of pitviper chemoperception, we analyzed a full-genome assembly for Crotalus adamanteus, the eastern diamondback rattlesnake. We identified hundreds of genes encoding both olfactory receptors (ORs; 362 full-length genes) and type-2 vomeronasal receptors (V2Rs; 430 full-length genes). Many chemoreceptor genes are organized into large tandem repeat arrays. Comparative analysis of V2R orthologs across squamates demonstrates how gene array expansion and contraction underlies the evolution of the chemoreceptor repertoire, which likely reflects shifts in life history traits. Chromosomal assignments of chemosensory genes identified sex chromosome specific chemoreceptor genes, providing gene candidates underlying observed sex-specific chemosensory-based behaviors. We detected widespread episodic evolution in the extracellular, ligand-binding domains of both ORs and V2Rs, suggesting the diversification of chemoreceptors is driven by transient periods of positive selection. We provide a robust genetic framework for studying pitviper chemosensory ecology and evolution.

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.

The Retinal Basis of Vertebrate Color Vision

The jawless fish that were ancestral to all living vertebrates had four spectral cone types that were probably served by chromatic-opponent retinal circuits. Subsequent evolution of photoreceptor spectral sensitivities is documented for many vertebrate lineages, giving insight into the ecological adaptation of color vision. Beyond the photoreceptors, retinal color processing is best understood in mammals, especially the blueON system, which opposes short- against long-wavelength receptor responses. For other vertebrates that often have three or four types of cone pigment, new findings from zebrafish are extending older work on teleost fish and reptiles to reveal rich color circuitry. Here, horizontal cells establish diverse and complex spectral responses even in photoreceptor outputs. Cone-selective connections to bipolar cells then set up color-opponent synaptic layers in the inner retina, which lead to a large variety of color-opponent channels for transmission to the brain via retinal ganglion cells.

Vision using multiple distinct rod opsins in deep-sea fishes

Vertebrate vision is accomplished through light-sensitive photopigments consisting of an opsin protein bound to a chromophore. In dim light, vertebrates generally rely on a single rod opsin [rhodopsin 1 (RH1)] for obtaining visual information. By inspecting 101 fish genomes, we found that three deep-sea teleost lineages have independently expanded their RH1 gene repertoires. Among these, the silver spinyfin (Diretmus argenteus) stands out as having the highest number of visual opsins in vertebrates (two cone opsins and 38 rod opsins). Spinyfins express up to 14 RH1s (including the most blueshifted rod photopigments known), which cover the range of the residual daylight as well as the bioluminescence spectrum present in the deep sea. Our findings present molecular and functional evidence for the recurrent evolution of multiple rod opsin-based vision in vertebrates.