Convergent Phenotypic Evolution Research Articles

Adaptive Evolution of Pseudomonas aeruginosa in Human Airways Shows Phenotypic Convergence Despite Diverse Patterns of Genomic Changes.

Selective forces in the environment drive bacterial adaptation to novel niches, choosing the fitter variants in the population. However, in dynamic and changing environments, the evolutionary processes controlling bacterial adaptation are difficult to monitor. Here, we follow 9 people with cystic fibrosis chronically infected with Pseudomonas aeruginosa, as a proxy for bacterial adaptation. We identify and describe the bacterial changes and evolution occurring between 15 and 35 yr of within-host evolution. We combine whole-genome sequencing, RNA sequencing, and metabolomics and compare the evolutionary trajectories directed by the adaptation of 4 different P. aeruginosa lineages to the lung. Our data suggest divergent evolution at the genomic level for most of the genes, with signs of convergent evolution with respect to the acquisition of mutations in regulatory genes, which drive the transcriptional and metabolomic program at late time of evolution. Metabolomics further confirmed convergent adaptive phenotypic evolution as documented by the reduction of the quorum-sensing molecules acyl-homoserine lactone, phenazines, and rhamnolipids (except for quinolones). The modulation of the quorum-sensing repertoire suggests that similar selective forces characterize at late times of evolution independent of the patient. Collectively, our data suggest that similar environments and similar P. aeruginosa populations in the patients at prolonged time of infection are associated with an overall reduction of virulence-associated features and phenotypic convergence.

Speciation and repeated origins of hypertrophied lips in parallel adaptive radiations of cyprinid fish from East Africa.

The evolution of convergent phenotypes is one of the most interesting phenomena of repeated adaptive radiations. Here, we examined the repeated patterns of thick-lipped or "rubberlip" phenotype of cyprinid fish of the genus Labeobarbus discovered in riverine environments of the Ethiopian Highlands, East Africa. To test the adaptive value of thickened lips, identify the ecological niche of the thick-lipped ecomorphs, and test whether these ecomorphs are the products of adaptive divergence, we studied six sympatric pairs of ecomorphs with hypertrophied lips and the normal lip structure from different riverine basins. Trophic morphology, diet, stable isotope (δ15N and δ13C) signatures, as well as mtDNA markers and genome-wide SNP variation, were analyzed. Our results show that thick-lipped ecomorphs partition trophic resources with generalized ecomorphs in only one-half of the examined sympatric pairs despite the pronounced divergence in lip structure. In these thick-lipped ecomorphs that were trophically diverged, the data on their diet along with the elevated 15N values suggest an invertivorous specialization different from the basal omnivorous-detritivouros feeding mode of the generalized ecomorphs. Genetic data confirmed an independent and parallel origin of all six lipped ecomorphs. Yet, only one of those six thick-lipped ecomorphs had a notable genetic divergence with sympatric non-lipped ecomorphs based on nuclear SNPs data (F ST = 0.21). Sympatric pairs can be sorted by combinations of phenotypic, ecological, and genetic divergence from an ecologically non-functional mouth polymorphism via ecologically functional polymorphism to a matured speciation stage via divergent evolution.

Evolution of immune genes in island birds: reduction in population sizes can explain island syndrome

Shared ecological conditions encountered by species that colonize islands often lead to the evolution of convergent phenotypes, commonly referred to as “island syndrome”. Reduced immune functions have been previously proposed to be part of this syndrome, as a consequence of the reduced diversity of pathogens on island ecosystems. According to this hypothesis, immune genes are expected to exhibit genomic signatures of relaxed selection pressure in island species. In this study, we used comparative genomic methods to study immune genes in island species (N = 20) and their mainland relatives (N = 14). We gathered public data as well as generated new data on innate (TLR: Toll-Like Receptors, BD: Beta Defensins) and acquired immune genes (MHC: Major Histocompatibility Complex classes I and II), but also on hundreds of genes with various immune functions. As a control, we used a set of 97 genes, not known to be involved in immune functions based on the literature, to account for the increased drift effects of the lower effective population sizes in island species. We used synonymous and non-synonymous variants to estimate the selection pressure acting on immune genes. We found that BDs and TLRs have higher ratios of non-synonymous over synonymous polymorphisms (Pn/Ps) than randomly selected control genes, suggesting that they evolve under a different selection regime. However, simulations show that this is unlikely to be explained by ongoing positive selection or balancing selection. For the MHC genes, which evolve under balancing selection, we used simulations to estimate the impact of population size variation. We found a significant effect of drift on immune genes of island species leading to a reduction in genetic diversity and efficacy of selection. However, the intensity of relaxed selection was not significantly different from control genes, except for MHC class II genes. These genes exhibit a significantly higher level of non-synonymous loss of polymorphism than expected assuming only drift and evolution under frequency dependent selection, possibly due to a reduction of extracellular parasite communities on islands. Overall, our results showed that demographic effects lead to a decrease in the immune functions of island species, but the relaxed selection that is expected to be caused by a reduced parasite pressure may only occur in some categories of immune genes.

Convergent Phenotypic Evolution of Rhodopsin for Dim-Light Sensing across Deep-Diving Vertebrates

Rhodopsin comprises an opsin attached to a retinal chromophore and is the only visual pigment conferring dim-light vision in vertebrates. On activation by photons, the retinal group becomes detached from the opsin, which is then inactive until it is recharged. Of all vertebrate species, those that dive face unique visual challenges, experiencing rapid decreases in light level and hunting in near darkness. Here, we combine sequence analyses with functional assays to show that the rhodopsin pigments of four divergent lineages of deep-diving vertebrates have undergone convergent increases in their retinal release rate. We compare gene sequences and detect parallel amino acids between penguins and diving mammals and perform mutagenesis to show that a single critical residue fully explains the observed increases in retinal release rate in both the emperor penguin and beaked whale. At the same time, we find that other shared sites have no significant effect on retinal release, implying that convergence does not always signify adaptive significance. We propose that accelerated retinal release confers rapid rhodopsin recharging, enabling the visual systems of diving species to adjust quickly to changing light levels as they descend through the water column. This contrasts with nocturnal species, where adaptation to darkness has been attributed to slower retinal release rates.

Developmental Transcriptomics Reveals a Gene Network Driving Mimetic Color Variation in a Bumble Bee.

A major goal of evolutionary genetics and evo-devo is to understand how changes in genotype manifest as changes in phenotype. Bumble bees display remarkable color pattern diversity while converging onto numerous regional Müllerian mimicry patterns, thus enabling exploration of the genetic mechanisms underlying convergent phenotypic evolution. In western North America, multiple bumble bee species converge onto local mimicry patterns through parallel shifts of midabdominal segments from red to black. It was previously demonstrated that a Hox gene, Abd-B, is the key regulator of the phenotypic switch in one of these species, Bombus melanopygus, however, the mechanism by which Abd-B regulates color differentiation remains unclear. Using tissue/stage-specific transcriptomic analysis followed by qRT–PCR validation, this study reveals a suite of genes potentially involved downstream of Abd-B during color pattern differentiation. The data support differential genes expression of not only the first switch gene Abd-B, but also an intermediate developmental gene nubbin, and a whole suite of downstream melanin and redox genes that together reinforce the observed eumelanin (black)-pheomelanin (red) ratios. These include potential genes involved in the production of insect pheomelanins, a pigment until recently not thought to occur in insects and thus lacking known regulatory enzymes. The results enhance understanding of pigmentation gene networks involved in bumble bee color pattern development and diversification, while providing insights into how upstream regulators such as Hox genes interact with downstream morphogenic players to facilitate this adaptive phenotypic radiation.

First approach to pod dehiscence in faba bean: genetic and histological analyses

Pod dehiscence causes important yield losses in cultivated crops and therefore has been a key trait strongly selected against in crop domestication. In spite of the growing knowledge on the genetic basis of dehiscence in different crops, no information is available so far for faba bean. Here we conduct the first comprehensive study for faba bean pod dehiscence by combining, linkage mapping, comparative genomics, QTL analysis and histological examination of mature pods. Mapping of dehiscence-related genes revealed conservation of syntenic blocks among different legumes. Three QTLs were identified in faba bean chromosomes II, IV and VI, although none of them was stable across years. Histological analysis supports the convergent phenotypic evolution previously reported in cereals and related legume species but revealed a more complex pattern in faba bean. Contrary to common bean and soybean, the faba bean dehiscence zone appears to show functional equivalence to that described in crucifers. The lignified wall fiber layer, which is absent in the paucijuga primitive line Vf27, or less lignified and vacuolated in other dehiscent lines, appears to act as the major force triggering pod dehiscence in this species. While our findings, provide new insight into the mechanisms underlying faba bean dehiscence, full understanding of the molecular bases will require further studies combining precise phenotyping with genomic analysis.

Convergent molecular evolution among ash species resistant to the emerald ash borer.

Recent studies show that molecular convergence plays an unexpectedly common role in the evolution of convergent phenotypes. We exploited this phenomenon to find candidate loci underlying resistance to the emerald ash borer (EAB; Agrilus planipennis), the USA’s most costly invasive forest insect to date, within the pan-genome of ash trees (the genus Fraxinus). We show that EAB-resistant taxa occur within three independent phylogenetic lineages. In genomes from these resistant lineages, we detect 53 genes with evidence of convergent amino acid evolution. Gene tree reconstruction indicates that for 48 of these candidates, the convergent amino acids are more likely to have arisen via independent evolution than by another process, such as hybridisation or incomplete lineage sorting. Seven of the candidate genes have putative roles connected to the phenylpropanoid biosynthesis pathway and 17 relate to herbivore recognition, defence signalling or programmed cell death. Evidence for loss-of-function mutations among these candidates is more frequent in susceptible species, than in resistant ones. Our results on evolutionary relationships, variability in resistance, and candidate genes for defence response within the ash genus could inform breeding for EAB resistance, facilitating ecological restoration in areas this beetle has invaded.

Convergent evolution of seasonal camouflage in response to reduced snow cover across the snowshoe hare range.

Determining how different populations adapt to similar environments is fundamental to understanding the limits of adaptation under changing environments. Snowshoe hares (Lepus americanus) typically molt into white winter coats to remain camouflaged against snow. In some warmer climates, hares have evolved brown winter camouflage-an adaptation that may spread in response to climate change. We used extensive range-wide genomic data to (1) resolve broad patterns of population structure and gene flow and (2) investigate the factors shaping the origins and distribution of winter-brown camouflage variation. In coastal Pacific Northwest (PNW) populations, winter-brown camouflage is known to be determined by a recessive haplotype at the Agouti pigmentation gene. Our phylogeographic analyses revealed deep structure and limited gene flow between PNW and more northern Boreal populations, where winter-brown camouflage is rare along the range edge. Genome sequencing of a winter-brown snowshoe hare from Alaska shows that it lacks the winter-brown PNW haplotype, reflecting a history of convergent phenotypic evolution. However, the PNW haplotype does occur at low frequency in a winter-white population from Montana, consistent with the spread of a locally deleterious recessive variant that is masked from selection when rare. Simulations of this population further show that this masking effect would greatly slow the selective increase of the winter-brown Agouti allele should it suddenly become beneficial (e.g., owing to dramatic declines in snow cover). Our findings underscore how allelic dominance can shape the geographic extent and rate of convergent adaptation in response to rapidly changing environments.

Ontological analyses reveal clinically-significant clear cell renal cell carcinoma subtypes with convergent evolutionary trajectories into an aggressive type.

BackgroundClear cell renal cell carcinoma (ccRCC) is a particularly challenging tumor type because of its extensive phenotypic variability as well as intra-tumoral heterogeneity (ITH). Clinically, this complexity has been reduced to a handful of pathological variables such as stage, grade and necrosis, but these variables fail to capture the breadth of the disease. How different phenotypes affect patient prognosis and influence therapeutic response is poorly understood. Extensive ITH illustrates remarkable plasticity, providing a framework to study tumor evolution. While multiregional genomic analyses have shown evolution from an ancient clone that acquires metastatic competency over time, these studies have been conducted agnostic to morphological cues and phenotypic plasticity.MethodsWe established a systematic ontology of ccRCC phenotypic variability by developing a multi-scale framework along three fundamental axes: tumor architecture, cytology and the microenvironment. We defined 33 parameters, which we comprehensively evaluated in 549 consecutive ccRCCs retrospectively. We systematically evaluated the impact of each parameter on patient outcomes, and assessed their contribution through multivariate analyses. We measured therapeutic impact in the context of anti-angiogenic therapies. We applied dimensionality reduction by t-distributed stochastic neighbor embedding (t-SNE) algorithms to tumor architectures for the study of tumor evolution superimposing tumor size and grade vectors. Evolutionary models were refined through empirical analyses of directed evolution of tumor intravascular extensions, and metastatic competency (as determined by tumor reconstitution in a heterologous host).FindingsWe discovered several novel ccRCC phenotypes, developed an integrated taxonomy, and identified features that improve current prognostic models. We identified a subset of ccRCCs refractory to anti-angiogenic therapies. We developed a model of tumor evolution, which revealed converging evolutionary trajectories into an aggressive type.InterpretationThis work serves as a paradigm for deconvoluting tumor complexity and illustrates how morphological analyses can improve our understanding of ccRCC pleiotropy. We identified several subtypes associated with aggressive biology, and differential response to targeted therapies. By analyzing patterns of spatial and temporal co-occurrence, intravascular tumor extensions and metastatic competency, we were able to identify distinct trajectories of convergent phenotypic evolution.

A functional enrichment test for molecular convergent evolution finds a clear protein-coding signal in echolocating bats and whales

Distantly related species entering similar biological niches often adapt by evolving similar morphological and physiological characters. How much genomic molecular convergence (particularly of highly constrained coding sequence) contributes to convergent phenotypic evolution, such as echolocation in bats and whales, is a long-standing fundamental question. Like others, we find that convergent amino acid substitutions are not more abundant in echolocating mammals compared to their outgroups. However, we also ask a more informative question about the genomic distribution of convergent substitutions by devising a test to determine which, if any, of more than 4,000 tissue-affecting gene sets is most statistically enriched with convergent substitutions. We find that the gene set most overrepresented (q-value = 2.2e-3) with convergent substitutions in echolocators, affecting 18 genes, regulates development of the cochlear ganglion, a structure with empirically supported relevance to echolocation. Conversely, when comparing to nonecholocating outgroups, no significant gene set enrichment exists. For aquatic and high-altitude mammals, our analysis highlights 15 and 16 genes from the gene sets most affected by molecular convergence which regulate skin and lung physiology, respectively. Importantly, our test requires that the most convergence-enriched set cannot also be enriched for divergent substitutions, such as in the pattern produced by inactivated vision genes in subterranean mammals. Showing a clear role for adaptive protein-coding molecular convergence, we discover nearly 2,600 convergent positions, highlight 77 of them in 3 organs, and provide code to investigate other clades across the tree of life.

The population genetics of structural variants in grapevine domestication.

Structural variants (SVs) are a largely unexplored feature of plant genomes. Little is known about the type and size of SVs, their distribution among individuals and, especially, their population dynamics. Understanding these dynamics is critical for understanding both the contributions of SVs to phenotypes and the likelihood of identifying them as causal genetic variants in genome-wide associations. Here, we identify SVs and study their evolutionary genomics in clonally propagated grapevine cultivars and their outcrossing wild progenitors. To catalogue SVs, we assembled the highly heterozygous Chardonnay genome, for which one in seven genes is hemizygous based on SVs. Using an integrative comparison between Chardonnay and Cabernet Sauvignon genomes by whole-genome, long-read and short-read alignment, we extended SV detection to population samples. We found that strong purifying selection acts against SVs but particularly against inversion and translocation events. SVs nonetheless accrue as recessive heterozygotes in clonally propagated lineages. They also define outlier regions of genomic divergence between wild and cultivated grapevines, suggesting roles in domestication. Outlier regions include the sex-determination region and the berry colour locus, where independent large, complex inversions have driven convergent phenotypic evolution.

The role of pollination effectiveness on the attributes of interaction networks: from floral visitation to plant fitness.

Network analysis is a powerful tool to understand community-level plant-pollinator interactions. We evaluated the role of floral visitors on plant fitness through a series of pollination exclusion experiments to test the effectiveness of pollinators of an Ipomoea community in the Pacific coast of Mexico, including: (1) all flower visitors, (2) visitors that contact the reproductive organs, (3) visitors that deposit pollen on stigmas, and (4) visitors that mediate fruit and seed production. Our results show that networks built from effective pollination interactions are smaller, less connected, more specialized and modular than floral visitor networks. Modules are associated with pollinator functional groups and they provide strong support for pollination syndromes only when non-effective interactions are excluded. In contrast to other studies, the analyzed networks are not nested. Our results also show that only 59% of floral visitors were legitimate pollinators that contribute to seed production. Furthermore, only 27% of the links in visitation network resulted in seed production. Our study shows that plant-pollination networks that consider effectiveness measures of pollination in addition to floral visitation provide insightful information about the different role floral visitors play in a community, encompassing a large number of commensalistic/antagonistic interactions and the more restricted set of mutualistic relationships that underlie the evolution of convergent floral phenotypes in plants.

Developmental control of convergent floral pigmentation across evolutionary timescales.

Convergent phenotypic evolution has been widely documented across timescales, from populations, to species, to major lineages. The extent to which convergent phenotypes arise from convergent genetic and developmental mechanisms remains an open question, although studies to-date reveal examples of both similar and different underlying mechanisms. This variation likely relates to a range of factors, including the genetic architecture of the trait and selective filtering of mutations over time. Here we focus on floral pigmentation, and examine the degree of developmental convergence between white-flowered lineages and white morphs within pigmented species. Using the model clade Iochrominae, we find that white morphs and white-flowered species are biochemically convergent, sharing an absence of colorful anthocyanin pigments. Regression analyses suggest that the expression levels of upstream genes are the strongest drivers of total pigmentation across species, although white species also show sharp down-regulation of the downstream genes. The white morphs do not share this pattern and present overall expression profiles more similar to the pigmented species. These results suggest that the mechanisms underlying variation within populations differ from those which give rise to fixed differences between species. Future work will aim to uncover the genetic changes responsible for this developmental non-convergence.

Molecular convergence and positive selection associated with the evolution of symbiont transmission mode in stony corals.

Heritable symbioses have been critical for the evolution of life. The genetic consequences of evolving a heritable symbiosis from the perspective of the symbiont are well established, but concomitant changes in the host remain unresolved. In stony corals, heritable, vertical transmission has evolved repeatedly, providing a unique opportunity to investigate the genomic basis of this complex trait. We conducted a comparative analysis of 25 coral transcriptomes to identify orthologous genes exhibiting signatures of positive selection and convergent amino acid substitutions in vertically transmitting lineages. The frequency of convergence events tends to be higher among vertically transmitting lineages, consistent with the proposed role of selection in driving the evolution of convergent transmission mode phenotypes. Of 10 774 orthologous genes, 403 exhibited at least one molecular convergence event and evidence of positive selection in at least one vertically transmitting lineage. Functional enrichments among these top candidate genes include processes previously implicated in symbiosis including endocytosis, immune response, cytoskeletal protein binding and cytoplasmic membrane-bounded vesicles. Finally, several novel candidates were identified among 100 genes showing evidence of positive selection at the particular convergence event, highlighting the value of our approach for generating new insight into host mechanisms associated with the evolution of heritable symbioses.

Convergent regulatory evolution and loss of flight in paleognathous birds.

A core question in evolutionary biology is whether convergent phenotypic evolution is driven by convergent molecular changes in proteins or regulatory regions. We combined phylogenomic, developmental, and epigenomic analysis of 11 new genomes of paleognathous birds, including an extinct moa, to show that convergent evolution of regulatory regions, more so than protein-coding genes, is prevalent among developmental pathways associated with independent losses of flight. A Bayesian analysis of 284,001 conserved noncoding elements, 60,665 of which are corroborated as enhancers by open chromatin states during development, identified 2355 independent accelerations along lineages of flightless paleognaths, with functional consequences for driving gene expression in the developing forelimb. Our results suggest that the genomic landscape associated with morphological convergence in ratites has a substantial shared regulatory component.

Genetic and lutionary mechanisms of convergent phenotypes in animals

It is not uncommon that similar phenotypes independently evolved in species of evolutionarily distant lineages, known as convergent evolution. The occurrence of convergent phenotypes is widely believed to promote adaptations to the same or similar environments of evolutionarily distant species by increasing their fitness. More and more studies show that the similar or same genetic variations play key roles in origination and evolution of convergent phenotypes. As genomic data accumulating rapidly, some researchers systematically explored the relationships between convergent phenotypes and convergent genetic variations on a genome-wide scale, comprehensively revealed molecular bases in genetics and evolution of convergent phenotypes, and identified many genes undergoing convergent evolution that are closely related to convergent phenotypes. Nevertheless, it is difficult to verify strictly whether convergent variations at the molecular level are under Darwinian natural selection. Except for convergent variations satisfying functional convergence of the related genes, it is more important to demonstrate that molecular convergence indeed leads to occurrence of convergent phenotypes and significantly increases the fitness of the species undergoing convergent evolution.

Convergent Evolution of the Seed Shattering Trait.

Loss of seed shattering is a key trait in crop domestication, particularly for grain crops. For wild plants, seed shattering is a crucial mechanism to achieve greater fitness, although in the agricultural context, this mechanism reduces harvesting efficiency, especially under dry conditions. Loss of seed shattering was acquired independently in different monocotyledon and dicotyledon crop species by ‘convergent phenotypic evolution’, leading to similar low dehiscent and indehiscent phenotypes. Here, the main aim is to review the current knowledge about seed shattering in crops, in order to highlight the tissue modifications that underlie the convergent phenotypic evolution of reduced shattering in different types of fruit, from the silique of Brassicaceae species, to the pods of legumes and spikes of cereals. Emphasis is given to legumes, with consideration of recent data obtained for the common bean. The current review also discusses to what extent convergent phenotypes arose from parallel changes at the histological and/or molecular levels. For this reason, an overview is included of the main findings relating to the genetic control of seed shattering in the model species Arabidopsis thaliana and in other important crops.

Genomic dissection of pod shattering in common bean: mutations at non-orthologous loci at the basis of convergent phenotypic evolution under domestication of leguminous species.

The complete or partial loss of shattering ability occurred independently during the domestication of several crops. Therefore, the study of this trait can provide an understanding of the link between phenotypic and molecular convergent evolution. The genetic dissection of 'pod shattering' in Phaseolus vulgaris is achieved here using a population of introgression lines and next-generation sequencing techniques. The 'occurrence' of the indehiscent phenotype (indehiscent versus dehiscent) depends on a major locus on chromosome 5. Furthermore, at least two additional genes are associated with the 'level' of shattering (number of shattering pods per plant: low versus high) and the 'mode' of shattering (non-twisting versus twisting pods), with all of these loci contributing to the phenotype by epistatic interactions. Comparative mapping indicates that the major gene identified on common bean chromosome 5 corresponds to one of the four quantitative trait loci for pod shattering in Vigna unguiculata. None of the loci identified comprised genes that are homologs of the known shattering genes in Glycine max. Therefore, although convergent domestication can be determined by mutations at orthologous loci, this was only partially true for P.vulgaris and V.unguiculata, which are two phylogenetically closely related crop species, and this was not the case for the more distant P.vulgaris and G.max. Conversely, comparative mapping suggests that the convergent evolution of the indehiscent phenotype arose through mutations in different genes from the same underlying gene networks that are involved in secondary cell-wall biosynthesis and lignin deposition patterning at the pod level.

Parallel selection on a dormancy gene during domestication of crops from multiple families.

Domesticated species often exhibit convergent phenotypic evolution, termed the domestication syndrome, of which loss of seed dormancy is a component. To date, dormancy genes that contribute to parallel domestication across different families have not been reported. Here, we cloned the classical stay-green G gene from soybean and found that it controls seed dormancy and showed evidence of selection during soybean domestication. Moreover, orthologs in rice and tomato also showed evidence of selection during domestication. Analysis of transgenic plants confirmed that orthologs of G had conserved functions in controlling seed dormancy in soybean, rice, and Arabidopsis. Functional investigation demonstrated that G affected seed dormancy through interactions with NCED3 and PSY and in turn modulated abscisic acid synthesis. Therefore, we identified a gene responsible for seed dormancy that has been subject to parallel selection in multiple crop families. This may help facilitate the domestication of new crops.

Convergent phenotypic evolution of the visual system via different molecular routes: How Neotropical cichlid fishes adapt to novel light environments.

How predictable is evolution? This remains a fundamental but contested issue in evolutionary biology. When independent lineages colonize the same environment, we are presented with a natural experiment that allows us to ask if genetic and ecological differences promote species‐specific evolutionary outcomes or whether species phenotypically evolve in a convergent manner in response to shared selection pressures. If so, are the molecular mechanisms underlying phenotypic convergence the same? In Nicaragua, seven species of cichlid fishes concurrently colonized two novel photic environments. Hence, their visual system represents a compelling model to address these questions, particularly since the adaptive value of phenotypic changes is well‐understood. By analyzing retinal transcriptomes, we found that differential expression of genes responsible for color vision (cone opsins and cyp27c1) produced rapid and mostly convergent changes of predicted visual sensitivities. Notably, these changes occurred in the same direction in all species although there were differences in underlying gene expression patterns illustrating nonconvergence at the molecular level. Adaptive phenotypes evolved deterministically, even when species differ substantially in ecology and genetic variation. This provides strong evidence that phenotypic evolution of the visual system occurred in response to similar selective forces of the photic environment.