The loss of previously adaptive traits is typically linked to relaxation in selection, yet the molecular steps leading to such repeated losses are rarely known. Molecular studies of loss have tended to focus on gene sequences alone, but overlooking other aspects of protein expression might underestimate phenotypic diversity. Insights based almost solely on opsin gene evolution, for instance, have made mammalian color vision a textbook example of phenotypic loss. We address this gap by investigating retention and loss of opsin genes, transcripts and proteins across ecologically diverse noctilionoid bats. We find multiple, independent losses of short-wave sensitive opsins. Mismatches between putatively functional DNA sequences, mRNA transcripts, and proteins implicate transcriptional and post-transcriptional processes in the ongoing loss of S-opsins in some noctilionoid bats. Our results provide a snapshot of evolution in progress during phenotypic trait loss, and suggest vertebrate visual phenotypes cannot always be predicted from genotypes alone.
R code and data files for calculating ecological correlates of opsin presence and density
Sadier_et_al_Opsins.zip
RHO_cds_al
RHO sequence alignment used for PAML analyses. Sequences were obtained from novel RNA-Seq data, publicly available genomes (GenBank), and unpublished genomes (Rossiter Lab and the Broad Institute). The final alignment included an additional sequence from Phyllostomus discolor available from https://vgp.github.io/genomeark/Phyllostomus_discolor/GenomeArk. The terminal stop codon has been removed.
OPN1SW_cds_al
OPN1SW sequence alignment used for PAML analyses. Sequences were obtained from novel RNA-Seq data, publicly available genomes (GenBank), and unpublished genomes (Rossiter Lab and the Broad Institute). The final alignment included an additional sequence from Phyllostomus discolor available from https://vgp.github.io/genomeark/Phyllostomus_discolor/GenomeArk. The terminal stop codon has been removed. Sequences were aligned by eye, premature stop codons were masked with ‘Ns’ and columns containing insertions that shifted the translation frame were deleted to keep codons in frame.
OPN1LW_cds_al
OPN1LW sequence alignment used for PAML analyses. Sequences were obtained from novel RNA-Seq data, publicly available genomes (GenBank), and unpublished genomes (Rossiter Lab and the Broad Institute). The final alignment included an additional sequence from Phyllostomus discolor available from https://vgp.github.io/genomeark/Phyllostomus_discolor/GenomeArk. The terminal stop codon has been removed.