Enhanced catabolism of glycine betaine and derivatives provides improved osmotic stress protection in Methylorubrum extorquens PA1.

Integration of metabolites into the overall metabolic network of a cell requires careful coordination dependent upon the ultimate usage of the metabolite. Different stoichiometric needs, and thus pathway fluxes, must exist for compounds destined for diverse uses, such as carbon sources, nitrogen sources, or stress-protective agents. Herein, we expand upon our previous work that highlighted the nature of glycine betaine (GB) metabolism in Methylobacteria to examine the utilization of GB-derivative compounds dimethylglycine (DMG) and sarcosine into Methylorubrum extorquens in different metabolic capacities, including as sole nitrogen and/or carbon sources. We isolated gain-of-function mutations that allowed M. extorquens PA1 to utilize dimethylglycine as a carbon source and dimethylglycine and sarcosine as nitrogen source. Characterization of mutants demonstrated selection for variants of the AraC-like regulator Mext_3735 that confer constitutive expression of the GB metabolic gene cluster, allowing direct utilization of the downstream GB derivatives. Finally, among the distinct isolates examined, we found that catabolism of the osmoprotectant used for selection (GB or dimethylglycine) enhanced osmotic stress resistance provided in the presence of that particular osmolyte. Thus, access to the carbon and nitrogen and osmoprotective effects of GB and DMG are made readily accessible through adaptive mutations. In M. extorquens PA1, the limitations to exploiting this group of compounds appear to exist predominantly at the levels of gene regulation and functional activity, rather than being constrained by transport or toxicity.IMPORTANCEOsmotic stress is a common challenge for bacteria colonizing the phyllosphere, where glycine betaine (GB) can be found as a prevalent osmoprotectant. Though Methylorubrum extorquens PA1 cannot use GB or its demethylation products, dimethylglycine (DMG) and sarcosine, as a sole carbon source, utilization is highly selectable via single nucleotide changes for both GB and DMG growth. The innate inability to use these compounds is due to limited flux through steps in the pathway and regulatory constraints. Herein, the characterization of the transcriptional regulator, Mext_3735 (GbdR), expands our understanding of the various roles in which GB derivatives can be used in M. extorquens PA1. Interestingly, increased catabolism of GB and derivatives does not interfere with, but rather improves, the ability of cells to thrive under increased salt stress conditions, suggesting that metabolic flux improves stress tolerance rather than providing a distinct tension between uses.

EfgA is a conserved formaldehyde sensor that leads to bacterial growth arrest in response to elevated formaldehyde.

Normal cellular processes give rise to toxic metabolites that cells must mitigate. Formaldehyde is a universal stressor and potent metabolic toxin that is generated in organisms from bacteria to humans. Methylotrophic bacteria such as Methylorubrum extorquens face an acute challenge due to their production of formaldehyde as an obligate central intermediate of single-carbon metabolism. Mechanisms to sense and respond to formaldehyde were speculated to exist in methylotrophs for decades but had never been discovered. Here, we identify a member of the DUF336 domain family, named efgA for enhanced formaldehyde growth, that plays an important role in endogenous formaldehyde stress response in M. extorquens PA1 and is found almost exclusively in methylotrophic taxa. Our experimental analyses reveal that EfgA is a formaldehyde sensor that rapidly arrests growth in response to elevated levels of formaldehyde. Heterologous expression of EfgA in Escherichia coli increases formaldehyde resistance, indicating that its interaction partners are widespread and conserved. EfgA represents the first example of a formaldehyde stress response system that does not involve enzymatic detoxification. Thus, EfgA comprises a unique stress response mechanism in bacteria, whereby a single protein directly senses elevated levels of a toxic intracellular metabolite and safeguards cells from potential damage.

Statistical modeling to quantify the uncertainty of FoldX-predicted protein folding and binding stability.

BACKGROUND: Computational methods of predicting protein stability changes upon missense mutations are invaluable tools in high-throughput studies involving a large number of protein variants. However, they are limited by a wide variation in accuracy and difficulty of assessing prediction uncertainty. Using a popular computational tool, FoldX, we develop a statistical framework that quantifies the uncertainty of predicted changes in protein stability. RESULTS: We show that multiple linear regression models can be used to quantify the uncertainty associated with FoldX prediction for individual mutations. Comparing the performance among models with varying degrees of complexity, we find that the model precision improves significantly when we utilize molecular dynamics simulation as part of the FoldX workflow. Based on the model that incorporates information from molecular dynamics, biochemical properties, as well as FoldX energy terms, we can generally expect upper bounds on the uncertainty of folding stability predictions of ± 2.9 kcal/mol and ± 3.5 kcal/mol for binding stability predictions. The uncertainty for individual mutations varies; our model estimates it using FoldX energy terms, biochemical properties of the mutated residue, as well as the variability among snapshots from molecular dynamics simulation. CONCLUSIONS: Using a linear regression framework, we construct models to predict the uncertainty associated with FoldX prediction of stability changes upon mutation. This technique is straightforward and can be extended to other computational methods as well.

Molecular Modeling Predicts Novel Antibody Escape Mutations in the Respiratory Syncytial Virus Fusion Glycoprotein.

Monoclonal antibodies are increasingly used for the prevention and/or treatment of viral infections. One caveat of their use is the ability of viruses to evolve resistance to antibody binding and neutralization. Computational strategies to identify viral mutations that may disrupt antibody binding would leverage the wealth of viral genomic sequence data to monitor for potential antibody-resistant mutations. The respiratory syncytial virus is an important pathogen for which monoclonal antibodies against the fusion (F) protein are used to prevent severe disease in high-risk infants. In this study, we used an approach that combines molecular dynamics simulations with FoldX to estimate changes in free energy in F protein folding and binding to the motavizumab antibody upon each possible amino acid change. We systematically selected 8 predicted escape mutations and tested them in an infectious clone. Consistent with our F protein stability predictions, replication-effective viruses were observed for each selected mutation. Six of the eight variants showed increased resistance to neutralization by motavizumab. Flow cytometry was used to validate the estimated (model-predicted) effects on antibody binding to F. Using surface plasmon resonance, we determined that changes in the on-rate of motavizumab binding were associated with the reduced affinity for two novel escape mutations. Our study empirically validated the accuracy of our molecular modeling approach and emphasized the role of biophysical protein modeling in predicting viral resistance to antibody-based therapeutics that can be used to monitor the emergence of resistant viruses and to design improved therapeutic antibodies. IMPORTANCE Respiratory syncytial virus (RSV) causes severe disease in young infants, particularly those with heart or lung diseases or born prematurely. Because no vaccine is currently available, monoclonal antibodies are used to prevent severe RSV disease in high-risk infants. While it is known that RSV evolves to avoid recognition by antibodies, screening tools that can predict which changes to the virus may lead to antibody resistance are greatly needed.

Paired Simulations and Experimental Investigations into the Calcium-Dependent Conformation of Albumin.

Serum albumin is the most abundant protein in blood plasma, and it is involved in multiple biological processes. Serum albumin has recently been adapted for improving biomaterial integration with bone tissue, and studies have shown the importance of this protein in bone repair and regeneration. However, the mechanism of action is not yet clear. In stark contrast, other studies have demonstrated that albumin blocks cell adhesion to surfaces, which is seen as a limitation to its bone healing role. These apparent contradictions suggest that the conformation of albumin facilitates its bioactivity, leading to enhanced bone repair. Serum albumin is known to play a major role in maintaining the calcium ion concentration in blood plasma. Due to the prevalence of calcium at bone repair and regeneration sites, it has been hypothesized that calcium binding to serum albumin triggers a conformational change, leading to bioactivity. In the current study, molecular modeling approaches including molecular docking, atomic molecular dynamics (MD) simulation, and coarse-grained MD simulation were used to test this hypothesis by investigating the conformational changes induced in bovine serum albumin by interaction with calcium ions. The computational results were qualitatively validated with experimental Fourier-transform infrared spectroscopy analysis. We find that free calcium ions in solution transiently bind with the three major loops in albumin, triggering a conformational change where N-terminal and C-terminal domains separate from each other in a partial unfolding process. The separation distance between these domains was found to correlate with the calcium ion concentration. The experimental data support the simulation results showing that albumin has enhanced conformational heterogeneity upon exposure to intermediate levels of calcium, without any significant secondary structure changes.

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.

Use of an electronic seizure diary in a randomized, controlled trial of natalizumab in adult participants with drug-resistant focal epilepsy.

OBJECTIVE: To analyze electronic diary (e-diary) use in a phase 2, randomized, controlled clinical trial (OPUS; NCT03283371) of natalizumab in adult participants with drug-resistant focal epilepsy. METHODS: We developed an e-diary, which incorporated an episodic seizure diary and a daily diary reminder, for use as the primary source to record participants' daily seizure activity in the OPUS phase 2 clinical trial. Participants and/or their designated caregivers made e-diary entries by selecting seizure descriptions generated in the participants' and/or caregivers' own words at the time of screening. Seizures and seizure-free days were reported for the current day and for up to 5 and 4 retrospective days, respectively. A record of seizure symptoms entered within the prior 5-day period was displayed on accessing the diary. Changes were not permitted in the e-diary once a seizure record was saved unless a data change request was made. A paper backup diary was available. RESULTS: E-diary entries (N = 15,176) from the 6-week baseline period and subsequent 24-week placebo-controlled period were analyzed for 66 adults who were randomized and dosed in the OPUS trial. The overall e-diary compliance, defined as the total number of days with any entry out of the total number of days in the baseline and placebo-controlled periods for all participants combined, was 83.6%. Caregivers made 190 (1.3%) e-diary entries. Day-of-event e-diary entries totaled 11,248 (74.1%). At least one paper backup diary was used by 36 (54.5%) participants. SIGNIFICANCE: Our data highlight that good e-diary compliance can be achieved across participants in randomized clinical trials in adult focal epilepsy. In addition to identifying and addressing any barriers that may prevent a minority of participants from achieving good e-diary compliance, consideration of e-diary elements, such as recall period and reporting of seizure-free days, will facilitate the most accurate data capture in epilepsy clinical trials.

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

Vision is the dominant sensory modality in many organisms for foraging, predator avoidance, and social behaviors including mate selection. Vertebrate visual perception is initiated when light strikes rod and cone photoreceptors within the neural retina of the eye. Sensitivity to individual colors, i.e., peak spectral sensitivities (λmax) of visual pigments, are a function of the type of chromophore and the amino acid sequence of the associated opsin protein in the photoreceptors. Large differences in peak spectral sensitivities can result from minor differences in amino acid sequence of cone opsins. To determine how minor sequence differences could result in large spectral shifts we selected a spectrally-diverse group of 14 teleost Rh2 cone opsins for which sequences and λmax are experimentally known. Classical molecular dynamics simulations were carried out after embedding chromophore-associated homology structures within explicit bilayers and water. These simulations revealed structural features of visual pigments, particularly within the chromophore, that contributed to diverged spectral sensitivities. Statistical tests performed on all the observed structural parameters associated with the chromophore revealed that a two-term, first-order regression model was sufficient to accurately predict λmax over a range of 452-528 nm. The approach was accurate, efficient and simple in that site-by-site molecular modifications or complex quantum mechanics models were not required to predict λmax. These studies identify structural features associated with the chromophore that may explain diverged spectral sensitivities, and provide a platform for future, functionally predictive opsin modeling.

Adaptive Mutations in RNA Polymerase and the Transcriptional Terminator Rho Have Similar Effects on Escherichia coli Gene Expression.

Modifications to transcriptional regulators play a major role in adaptation. Here, we compared the effects of multiple beneficial mutations within and between Escherichia coli rpoB, the gene encoding the RNA polymerase ß subunit, and rho, which encodes a transcriptional terminator. These two genes have harbored adaptive mutations in numerous E. coli evolution experiments but particularly in our previous large-scale thermal stress experiment, where the two genes characterized alternative adaptive pathways. To compare the effects of beneficial mutations, we engineered four advantageous mutations into each of the two genes and measured their effects on fitness, growth, gene expression and transcriptional termination at 42.2 °C. Among the eight mutations, two rho mutations had no detectable effect on relative fitness, suggesting they were beneficial only in the context of epistatic interactions. The remaining six mutations had an average relative fitness benefit of ∼20%. The rpoB mutations affected the expression of ∼1,700 genes; rho mutations affected the expression of fewer genes but most (83%) were a subset of those altered by rpoB mutants. Across the eight mutants, relative fitness correlated with the degree to which a mutation restored gene expression back to the unstressed, 37.0 °C state. The beneficial mutations in the two genes did not have identical effects on fitness, growth or gene expression, but they caused parallel phenotypic effects on gene expression and genome-wide transcriptional termination.

Impact of embryo co-culture with cumulus cells on pregnancy & implantation rate in patients undergoing in vitro fertilization using donor oocyte.

BACKGROUND & OBJECTIVES: Cumulus cell co-culture of embryo had been found to be beneficial for achieving better pregnancy and implantation rate (IR). The present study was aimed to evaluate efficiency of cumulus co-culture technique over simple culture of embryo in terms of pregnancy rate (PR) and IR in patients undergoing treatment for infertility using donor oocytes fertilized by intracytoplasmic sperm injection. METHODS: This was a quasi-experimental study between control and study groups. The primary endpoint was achievement of pregnancy. Control group included 508 women who underwent embryo development without cumulus cell co-culture and study group included 394 women who underwent embryo development with cumulus cell co-culture using donor's cumulus cells. RESULTS: The present study demonstrated a significant increase in the IR (37.2 vs 24.2%, P<0.001) and in PR (45.7 vs 37.8%, P<0.05) in study group than in control group. The PR and IR were found to be higher in study group, among all groups of women, grouped on the basis of different indications for use of donor oocytes. INTERPRETATION & CONCLUSIONS: Cumulus cell co-culture technique was found to be more effective than simple culture technique for embryo development in women undergoing treatment for infertility using donor oocytes fertilized by intracytoplasmic sperm injection.

Meningoencephalitis associated with disseminated sarcocystosis in a free-ranging moose (Alces alces) calf.

A wild moose (Alces alces) calf was presented for necropsy due to severe neurologic signs. Histopathologic examination revealed multisystemic inflammation with intralesional mature and immature schizonts. Schizonts in the brain reacted positively to Sarcocystis spp. polyclonal antibodies. Gene sequencing of PCR-amplified DNA identified the species as Sarcocystis alceslatrans.

Méningo-encéphalite associée avec une sarcocystose disséminée chez un jeune orignal en liberté(Alces alces). Un jeune orignal en liberté (Alces alces) a été présenté pour une nécropsie en raison de signes neurologiques graves. L'examen histopathologique a révélé une inflammation multisystémique avec des schizontes matures et immatures intralésionnels. Les schizontes du cerveau ont réagi positivement aux anticorps polyclonaux de Sarcocystis spp. Le séquençage des gènes de l'ADN par amplification en chaîne par polymérase a identifié l'espèce comme étant Sarcocystis alceslatrans.(Traduit par Isabelle Vallières).

Local generation and propagation of ripples along the septotemporal axis of the hippocampus.

A topographical relationship exists between the septotemporal segments of the hippocampus and their entorhinal-neocortical targets, but the physiological organization of activity along the septotemporal axis is poorly understood. We recorded sharp-wave ripple patterns in rats during sleep from the entire septotemporal axis of the CA1 pyramidal layer. Qualitatively similar ripples emerged at all levels. From the local seed, ripples traveled septally or temporally at a speed of ∼0.35 m/s, and the spatial spread depended on ripple magnitude. Ripples propagated smoothly across the septal and intermediate segments of the hippocampus, but ripples in the temporal segment often remained isolated. These findings show that ripples can combine information from the septal and intermediate hippocampus and transfer integrated signals downstream. In contrast, ripples that emerged in the temporal pole broadcast largely independent information to their cortical and subcortical targets.

Steered molecular dynamics simulations for studying protein-ligand interaction in cyclin-dependent kinase 5.

In this study, we applied steered molecular dynamics (SMD) simulations to investigate the unbinding mechanism of nine inhibitors of the enzyme cyclin-dependent kinase 5 (CDK5). The study had two major objectives: (i) to create a correlation between the unbinding force profiles and the inhibition activities of these compounds expressed as IC50 values; (ii) to investigate the unbinding mechanism and to reveal atomistic insights, which could help identify accessory binding sites and transient interactions. Overall, we carried out 1.35 µs of cumulative SMD simulations. We showed that SMD could qualitatively discriminate binders from nonbinders, while it failed to properly rank series of inhibitors, particularly when IC50 values were too similar. From a mechanistic standpoint, SMD provided useful insights related to transient and dynamical interactions, which could complement static description obtained by X-ray crystallography experiments. In conclusion, the present study represents a further step toward a systematic exploitation of SMD and other dynamical approaches in structure-based drug design and computational medicinal chemistry.

Exploring the ability of the MD+FoldX method to predict SARS-CoV-2 antibody escape mutations using large-scale data.

Antibody escape mutations pose a significant challenge to the effectiveness of vaccines and antibody-based therapies. The ability to predict these escape mutations with computer simulations would allow us to detect threats early and develop effective countermeasures, but a lack of large-scale experimental data has hampered the validation of these calculations. In this study, we evaluate the ability of the MD+FoldX molecular modeling method to predict escape mutations by leveraging a large deep mutational scanning dataset, focusing on the SARS-CoV-2 receptor binding domain. Our results show a positive correlation between predicted and experimental data, indicating that mutations with reduced predicted binding affinity correlate moderately with higher experimental escape fractions. We also demonstrate that better performance can be achieved using affinity cutoffs tailored to distinct antibody-antigen interactions rather than a one-size-fits-all approach. We find that 70% of the systems surpass the 50% precision mark, and demonstrate success in identifying mutations present in significant variants of concern and variants of interest. Despite promising results for some systems, our study highlights the challenges in comparing predicted and experimental values. It also emphasizes the need for new binding affinity methods with improved accuracy that are fast enough to estimate hundreds to thousands of antibody-antigen binding affinities.

Human Milk Protein-Derived Bioactive Peptides from In Vitro-Digested Colostrum Exert Antimicrobial Activities against Common Neonatal Pathogens.

Human milk reduces risk for necrotizing enterocolitis in preterm infants. Necrotizing enterocolitis occurs in the ileocecal region where thousands of milk protein-derived peptides have been released from digestion. Digestion-released peptides may exert bioactivity, such as antimicrobial and immunomodulatory activities, in the gut. In this study, we applied mass spectrometry-based peptidomics to characterize peptides present in colostrum before and after in vitro digestion. Sequence-based computational modeling was applied to predict peptides with antimicrobial activity. We identified more peptides in undigested samples, yet the abundances were much higher in the digested samples. Heatmapping demonstrated highly different peptide profiles between undigested and digested samples. Four peptides (αS1-casein [157-163], αS1-casein [157-165], ß-casein [153-159] and plasminogen [591-597]) were selected, synthesized and tested against common pathogenic bacteria associated with necrotizing enterocolitis. All four exhibited bacteriostatic, though not bactericidal, activities against Klebsiella aerogenes, Citrobacter freundii and Serratia marcescens, but not Escherichia coli.

The effect of mutations on binding interactions between the SARS-CoV-2 receptor binding domain and neutralizing antibodies B38 and CB6.

SARS-CoV-2 is the pathogen responsible for COVID-19 that has claimed over six million lives as of July 2022. The severity of COVID-19 motivates a need to understand how it could evolve to escape potential treatments and to find ways to strengthen existing treatments. Here, we used the molecular modeling methods MD + FoldX and PyRosetta to study the SARS-CoV-2 spike receptor binding domain (S-RBD) bound to two neutralizing antibodies, B38 and CB6 and generated lists of antibody escape and antibody strengthening mutations. Our resulting watchlist contains potential antibody escape mutations against B38/CB6 and consists of 211/186 mutations across 35/22 S-RBD sites. Some of these mutations have been identified in previous studies as being significant in human populations (e.g., N501Y). The list of potential antibody strengthening mutations that are predicted to improve binding of B38/CB6 to S-RBD consists of 116/45 mutations across 29/13 sites. These mutations could be used to improve the therapeutic value of these antibodies.

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.

Defining the HIV Capsid Binding Site of Nucleoporin 153.

The interaction between the HIV-1 capsid and human nucleoporin 153 (NUP153) is vital for delivering the HIV-1 preintegration complex into the nucleus via the nuclear pore complex. The interaction with the capsid requires a phenylalanine/glycine-containing motif in the C-terminus of NUP153 (NUP153C). This study used molecular modeling and biochemical assays to comprehensively determine the amino acids in NUP153 that are important for capsid interaction. Molecular dynamics, FoldX, and PyRosetta simulations delineated the minimal capsid binding motif of NUP153 based on the known structure of NUP153 bound to the HIV-1 capsid hexamer. Computational predictions were experimentally validated by testing the interaction of NUP153 with capsid using an in vitro binding assay and a cell-based TRIM-NUP153C restriction assay. This work identified eight amino acids from P1411 to G1418 that stably engage with capsid, with significant correlations between the interactions predicted by molecular models and empirical experiments. This validated the usefulness of this multidisciplinary approach to rapidly characterize the interaction between human proteins and the HIV-1 capsid. IMPORTANCE The human immunodeficiency virus (HIV) can infect nondividing cells by interacting with the host nuclear pore complex. The host nuclear pore protein NUP153 directly interacts with the HIV capsid to promote viral nuclear entry. This study used a multidisciplinary approach combining computational and experimental techniques to comprehensively map the effect of mutating the amino acids of NUP153 on HIV capsid interaction. This work showed a significant correlation between computational and empirical data sets, revealing that the HIV capsid interacted specifically with only six amino acids of NUP153. The simplicity of the interaction motif suggested other FG-containing motifs could also interact with the HIV-1 capsid. Furthermore, it was predicted that naturally occurring polymorphisms in human and nonhuman primates would disrupt NUP153 interaction with capsid, potentially protecting certain populations from HIV-1 infection.

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.