Distinct Evolutionary Patterns Between Two Duplicated Color Vision Genes Within Cyprinid Fishes

Evolution, inactivation and loss of short wavelength-sensitive opsin genes during the diversification of Neotropical cichlids

Natural variation in the number, expression and function of sensory genes in an organism's genome is often tightly linked to different ecological and evolutionary forces. Opsin genes, which code for the first step in visual transduction, are ideal models for testing how ecological factors such as light environment may influence visual system adaptation. Neotropical cichlid fishes are a highly ecologically diverse group that evolved in a variety of aquatic habitats, including black (stained), white (opaque) and clear waters. We used cross-species exon capture to sequence Neotropical cichlid short wavelength-sensitive (SWS) opsins, which mediate ultraviolet (UV) to blue visual sensitivity. Neotropical cichlid SWS1 opsin (UV-sensitive) underwent a relaxation of selective constraint during the early phases of cichlid diversification in South America, leading to pseudogenization and loss. Conversely, SWS2a (blue-sensitive) experienced a burst of episodic positive selection at the base of the South American cichlid radiation. This burst coincides with SWS1 relaxation and loss, and is consistent with findings in ecomorphological studies characterizing a period of extensive ecological divergence in Neotropical cichlids. We use ancestral sequence reconstruction and protein modelling to investigate mutations along this ancestral branch that probably modified SWS2a function. Together, our results suggest that variable light environments played a prominent early role in shaping SWS opsin diversity during the Neotropical cichlid radiation. Our results also illustrate that long-term evolution under light-limited conditions in South America may have reduced visual system plasticity; specifically, early losses of UV sensitivity may have constrained the evolutionary trajectory of Neotropical cichlid vision.

Phylogenetic analysis and expression differences of eye-related genes in cavefish genus Sinocyclocheilus

The adaptive evolution of visual systems has been observed in many cavefish. However, little is known about the molecular mechanisms underlying these adaptations, which include regressive changes such as eye degeneration. Here, we analyzed phylogenetic and expression patterns of 6 eye-related genes (crx, foxg1b, opn1sw2, otx2, rho and sox2) in 12 Sinocyclocheilus species from China, including 8 stygobionts and 4 stygophiles, and examined photoreceptor cell morphology of these species. Those eye-degenerated species of Sinocyclocheilus were polyphyletic and showed different degrees of photoreceptor defects in responses to cave environments. The eye loss and degeneration are the result of convergent evolution. Although S. anophthalmus grouped with the eye-normal species, it displayed not only a high degree of eye degeneration but also significant expression differences in eye-related genes compared with the eye-normal species. The gene foxg1b, which was determined to be under positive selection, might play an important role in the process of eye degeneration in S. anophthalmus based on differential expression. Eye-related gene expression and selection may have contributed to the polyphyly of the cave species. We examined gene expression and duplication in 6 eye-related genes and revealed that these genes displayed considerable diversity in relative expression in Sinocyclocheilus fishes. Otx2 and sox2 were significantly up-regulated in individual cave species, while the other 4 genes (crx, foxg1b, opn1sw2 and rho) were significantly down-regulated. These findings provide a valuable resource for elucidating molecular mechanisms associated with visual system evolution in cavefish.

Differentially-expressed opsin genes identified in Sinocyclocheilus cavefish endemic to China

Eye degeneration is a common troglomorphic character of cave-dwelling organisms. Comparing the morphology and molecular biology of cave species and their close surface relatives is a powerful tool for studying regressive eye evolution and other adaptive phenotypes. We compared two co-occurring and closely-related species of the fish genus Sinocyclocheilus, which is endemic to China and includes both surface- and cave-dwelling species. Sinocyclocheilus tileihornes, a cave species, had smaller eyes than Sinocyclocheilus angustiporus, a surface species. Histological and immunohistochemical analyses revealed that the cavefish had shorter cones and more disorderly rods than did the surface-dwelling species. Using quantitative PCR and in situ hybridization, we found that rhodopsin and a long-wavelength sensitive opsin had significantly lower expression levels in the cavefish. Furthermore, one of two short-wavelength-sensitive opsins was expressed at significantly higher levels in the cavefish. Changes in the expression of opsin genes may have played a role in the degeneration of cavefish eyes.

Evidence for repeated loss of selective constraint in rhodopsin of amblyopsid cavefishes (Teleostei: Amblyopsidae)

The genetic mechanisms underlying regressive evolution-the degeneration or loss of a derived trait--are largely unknown, particularly for complex structures such as eyes in cave organisms. In several eyeless animals, the visual photoreceptor rhodopsin appears to have retained functional amino acid sequences. Hypotheses to explain apparent maintenance of function include weak selection for retention of light-sensing abilities and its pleiotropic roles in circadian rhythms and thermotaxis. In contrast, we show that there has been repeated loss of functional constraint of rhodopsin in amblyopsid cavefishes, as at least three cave lineages have independently accumulated unique loss-of-function mutations over the last 10.3 Mya. Although several cave lineages still possess functional rhodopsin, they exhibit increased rates of nonsynonymous mutations that have greater effect on the structure and function of rhodopsin compared to those in surface lineages. These results indicate that functionality of rhodopsin has been repeatedly lost in amblyopsid cavefishes. The presence of a functional copy of rhodopsin in some cave lineages is likely explained by stochastic accumulation of mutations following recent subterranean colonization.

Complex patterns of divergence among green-sensitive (RH2a) African cichlid opsins revealed by Clade model analyses

BACKGROUND

Gene duplications play an important role in the evolution of functional protein diversity. Some models of duplicate gene evolution predict complex forms of paralog divergence; orthologous proteins may diverge as well, further complicating patterns of divergence among and within gene families. Consequently, studying the link between protein sequence evolution and duplication requires the use of flexible substitution models that can accommodate multiple shifts in selection across a phylogeny. Here, we employed a variety of codon substitution models, primarily Clade models, to explore how selective constraint evolved following the duplication of a green-sensitive (RH2a) visual pigment protein (opsin) in African cichlids. Past studies have linked opsin divergence to ecological and sexual divergence within the African cichlid adaptive radiation. Furthermore, biochemical and regulatory differences between the RH2aα and RH2aβ paralogs have been documented. It thus seems likely that selection varies in complex ways throughout this gene family.

RESULTS

Clade model analysis of African cichlid RH2a opsins revealed a large increase in the nonsynonymous-to-synonymous substitution rate ratio (ω) following the duplication, as well as an even larger increase, one consistent with positive selection, for Lake Tanganyikan cichlid RH2aβ opsins. Analysis using the popular Branch-site models, by contrast, revealed no such alteration of constraint. Several amino acid sites known to influence spectral and non-spectral aspects of opsin biochemistry were found to be evolving divergently, suggesting that orthologous RH2a opsins may vary in terms of spectral sensitivity and response kinetics. Divergence appears to be occurring despite intronic gene conversion among the tandemly-arranged duplicates.

CONCLUSIONS

Our findings indicate that variation in selective constraint is associated with both gene duplication and divergence among orthologs in African cichlid RH2a opsins. At least some of this variation may reflect an adaptive response to differences in light environment. Interestingly, these patterns only became apparent through the use of Clade models, not through the use of the more widely employed Branch-site models; we suggest that this difference stems from the increased flexibility associated with Clade models. Our results thus bear both on studies of cichlid visual system evolution and on studies of gene family evolution in general.

Opsin gene duplication and divergence in ray-finned fish

Opsin gene sequences were first reported in the 1980s. The goal of that research was to test the hypothesis that human opsins were members of a single gene family and that variation in human color vision was mediated by mutations in these genes. While the new data supported both hypotheses, the greatest contribution of this work was, arguably, that it provided the data necessary for PCR-based surveys in a diversity of other species. Such studies, and recent whole genome sequencing projects, have uncovered exceptionally large opsin gene repertoires in ray-finned fishes (taxon, Actinopterygii). Guppies and zebrafish, for example, have 10 visual opsin genes each. Here we review the duplication and divergence events that have generated these gene collections. Phylogenetic analyses revealed that large opsin gene repertories in fish have been generated by gene duplication and divergence events that span the age of the ray-finned fishes. Data from whole genome sequencing projects and from large-insert clones show that tandem duplication is the primary mode of opsin gene family expansion in fishes. In some instances gene conversion between tandem duplicates has obscured evolutionary relationships among genes and generated unique key-site haplotypes. We mapped amino acid substitutions at so-called key-sites onto phylogenies and this exposed many examples of convergence. We found that dN/dS values were higher on the branches of our trees that followed gene duplication than on branches that followed speciation events, suggesting that duplication relaxes constraints on opsin sequence evolution. Though the focus of the review is opsin sequence evolution, we also note that there are few clear connections between opsin gene repertoires and variation in spectral environment, morphological traits, or life history traits.

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.

Of cancer and cave fish

We propose that the drivers of carcinogenesis lie more in the adaptive changes that are enabled by local or systemic alterations of tissue architecture than in the genetic changes observed in cancer cells. A full understanding of cancer biology and therapy through a cataloguing of the cancer genome is unlikely unless it is integrated into an evolutionary and ecological context.