Molecular evolution of the betagamma lens crystallin superfamily: evidence for a retained ancestral function in gamma N crystallins?

The sensory shark: high-quality morphological, genomic and transcriptomic data for the small-spotted catshark Scyliorhinus canicula reveal the molecular bases of sensory organ evolution in jawed vertebrates

Cartilaginous fishes (chimaeras and elasmobranchs -sharks, skates and rays) hold a key phylogenetic position to explore the origin and diversifications of jawed vertebrates. Here, we report and integrate reference genomic, transcriptomic and morphological data in the small-spotted catshark Scyliorhinus canicula to shed light on the evolution of sensory organs. We first characterise general aspects of the catshark genome, confirming the high conservation of genome organisation across cartilaginous fishes, and investigate population genomic signatures. Taking advantage of a dense sampling of transcriptomic data, we also identify gene signatures for all major organs, including chondrichthyan specializations, and evaluate expression diversifications between paralogs within major gene families involved in sensory functions. Finally, we combine these data with 3D synchrotron imaging and in situ gene expression analyses to explore chondrichthyan-specific traits and more general evolutionary trends of sensory systems. This approach brings to light, among others, novel markers of the ampullae of Lorenzini electro-sensory cells, a duplication hotspot for crystallin genes conserved in jawed vertebrates, and a new metazoan clade of the Transient-receptor potential (TRP) family. These resources and results, obtained in an experimentally tractable chondrichthyan model, open new avenues to integrate multiomics analyses for the study of elasmobranchs and jawed vertebrates.

Proteome Analysis of Whole-Body Responses in Medaka Experimentally Exposed to Fish-Killing Dinoflagellate Karenia mikimotoi

Karenia mikimotoi is a well-known harmful algal bloom species. Blooms of this dinoflagellate have become a serious threat to marine life, including fish, shellfish, and zooplanktons and are usually associated with massive fish death. Despite the discovery of several toxins such as gymnocins and gymnodimines in K. mikimotoi, the mechanisms underlying the ichthyotoxicity of this species remain unclear, and molecular studies on this topic have never been reported. The present study investigates the fish-killing mechanisms of K. mikimotoi through comparative proteomic analysis. Marine medaka, a model fish organism, was exposed to K. mikimotoi for a three-part time period (LT25, LT50 and LT90). Proteins extracted from the whole fish were separated by using two-dimensional gel electrophoresis, and differentially expressed proteins were identified with reference to an untreated control. The change in fish proteomes over the time-course of exposure were analyzed. A total of 35 differential protein spots covering 19 different proteins were identified, of which most began to show significant change in expression levels at the earliest stage of intoxication. Among the 19 identified proteins, some are closely related to the oxidative stress responses, energy metabolism, and muscle contraction. We propose that oxidative stress-mediated muscle damage might explain the symptoms developed during the ichthyotoxicity test, such as gasping for breath, loss of balance, and body twitching. Our findings lay the foundations for more in-depth studies of the mechanisms of K. mikimotoi's ichthyotoxicity.

Protein refractive index increment is determined by conformation as well as composition

The refractive index gradient of the eye lens is controlled by the concentration and distribution of its component crystallin proteins, which are highly enriched in polarizable amino acids. The current understanding of the refractive index increment ([Formula: see text]) of proteins is described using an additive model wherein the refractivity and specific volume of each amino acid type contributes according to abundance in the primary sequence. Here we present experimental measurements of [Formula: see text] for crystallins from the human lens and those of aquatic animals under uniform solvent conditions. In all cases, the measured values are much higher than those predicted from primary sequence alone, suggesting that structural factors also contribute to protein refractive index.

Characterization of Embryo Transcriptome of Gynogenetic Olive Flounder Paralichthys olivaceus

Olive flounder Paralichthys olivaceus is an important commercially cultured marine flatfish in China, Korea, and Japan. Gynogenesis, via meiogynogenesis and mitogynogenesis, shows advantages in breeding and sex control, but the low survival rate, especially for mitogynogenesis, limits its application. In this study, we sequenced the embryo transcriptomes of gynogenetic haploid, meiogynogenetic diploid, mitogynogenetic diploid, and common diploid flounder and investigated their respective genetic characteristics by analyzing differentiated expressed genes. Compared with common diploid, the gynogenetic haploid showed significant downregulation in notch signaling and wingless-related integration site (Wnt) signaling pathways, which may be the source of haploid syndrome. In both meiogynogenesis and mitogynogenesis, several upregulated genes including complement C3, formin-2, and intelectin may be related to increased survival compared to the haploid. The downregulation of immune system and energy metabolism-related genes caused retarded development of gynogenetic diploids compared with the common diploid. These data provided new and important information for application of artificially induced gynogenesis to aquaculture.

Molecular evolutionary analysis of vertebrate transducins: a role for amino acid variation in photoreceptor deactivation

Transducin is a heterotrimeric G protein that plays a critical role in phototransduction in the rod and cone photoreceptor cells of the vertebrate retina. Rods, highly sensitive cells that recover from photoactivation slowly, underlie dim-light vision, whereas cones are less sensitive, recover more quickly, and underlie bright-light vision. Transducin deactivation is a critical step in photoreceptor recovery and may underlie the functional distinction between rods and cones. Rods and cones possess distinct transducin α subunits, yet they share a common deactivation mechanism, the GTPase activating protein (GAP) complex. Here, we used codon models to examine patterns of sequence evolution in rod (GNAT1) and cone (GNAT2) α subunits. Our results indicate that purifying selection is the dominant force shaping GNAT1 and GNAT2 evolution, but that GNAT2 has additionally been subject to positive selection operating at multiple phylogenetic scales; phylogeny-wide analysis identified several sites in the GNAT2 helical domain as having substantially elevated dN/dS estimates, and branch-site analysis identified several nearby sites as targets of strong positive selection during early vertebrate history. Examination of aligned GNAT and GAP complex crystal structures revealed steric clashes between several positively selected sites and the deactivating GAP complex. This suggests that GNAT2 sequence variation could play an important role in adaptive evolution of the vertebrate visual system via effects on photoreceptor deactivation kinetics and provides an alternative perspective to previous work that focused instead on the effect of GAP complex concentration. Our findings thus further the understanding of the molecular biology, physiology, and evolution of vertebrate visual systems.

The molecular refractive function of lens γ-Crystallins

γ-Crystallins constitute the major protein component in the nucleus of the vertebrate eye lens. Present at very high concentrations, they exhibit extreme solubility and thermodynamic stability to prevent scattering of light and formation of cataracts. However, functions beyond this structural role have remained mostly unclear. Here, we calculate molecular refractive index increments of crystallins. We show that all lens γ-crystallins have evolved a significantly elevated molecular refractive index increment, which is far above those of most proteins, including nonlens members of the βγ-crystallin family from different species. The same trait has evolved in parallel in crystallins of different phyla, including S-crystallins of cephalopods. A high refractive index increment can lower the crystallin concentration required to achieve a suitable refractive power of the lens and thereby reduce their propensity to aggregate and form cataracts. To produce a significant increase in the refractive index increment, a substantial global shift in amino acid composition is required, which can naturally explain the highly unusual amino acid composition of γ-crystallins and their functional homologues. This function provides a new perspective for interpreting their molecular structure.

Dynamic programming procedure for searching optimal models to estimate substitution rates based on the maximum-likelihood method

The substitution rate in a gene can provide valuable information for understanding its functionality and evolution. A widely used method to estimate substitution rates is the maximum-likelihood method implemented in the CODEML program in the PAML package. A limited number of branch models, chosen based on a priori information or an interest in a particular lineage(s), are tested, whereas a large number of potential models are neglected. A complementary approach is also needed to test all or a large number of possible models to search for the globally optional model(s) of maximum likelihood. However, the computational time for this search even in a small number of sequences becomes impractically long. Thus, it is desirable to explore the most probable spaces to search for the optimal models. Using dynamic programming techniques, we developed a simple computational method for searching the most probable optimal branch-specific models in a practically feasible computational time. We propose three search methods to find the optimal models, which explored O(n) (method 1) to O(n(2)) (method 2 and method 3) models when the given phylogeny has n branches. In addition, we derived a formula to calculate the number of all possible models, revealing the complexity of finding the optimal branch-specific model. We show that in a reanalysis of over 50 previously published studies, the vast majority obtained better models with significantly higher likelihoods than the conventional hypothesis model methods.

Explosive expansion of betagamma-crystallin genes in the ancestral vertebrate

In jawed vertebrates, βγ-crystallins are restricted to the eye lens and thus excellent markers of lens evolution. These βγ-crystallins are four Greek key motifs/two domain proteins, whereas the urochordate βγ-crystallin has a single domain. To trace the origin of the vertebrate βγ-crystallin genes, we searched for homologues in the genomes of a jawless vertebrate (lamprey) and of a cephalochordate (lancelet). The lamprey genome contains orthologs of the gnathostome βB1-, βA2- and γN-crystallin genes and a single domain γN-crystallin-like gene. It contains at least two γ-crystallin genes, but lacks the gnathostome γS-crystallin gene. The genome also encodes a non-lenticular protein containing βγ-crystallin motifs, AIM1, also found in gnathostomes but not detectable in the uro- or cephalochordate genome. The four cephalochordate βγ-crystallin genes found encode two-domain proteins. Unlike the vertebrate βγ-crystallins but like the urochordate βγ-crystallin, three of the predicted proteins contain calcium-binding sites. In the cephalochordate βγ-crystallin genes, the introns are located within motif-encoding region, while in the urochordate and in the vertebrate βγ-crystallin genes the introns are between motif- and/or domain encoding regions. Coincident with the evolution of the vertebrate lens an ancestral urochordate type βγ-crystallin gene rapidly expanded and diverged in the ancestral vertebrate before the cyclostomes/gnathostomes split. The β- and γN-crystallin genes were maintained in subsequent evolution, and, given the selection pressure imposed by accurate vision, must be essential for lens function. The γ-crystallin genes show lineage specific expansion and contraction, presumably in adaptation to the demands on vision resulting from (changes in) lifestyle.