Transcriptome Assemblies

De novo Trinity assemblies of the colubrid snake eye transcriptomes

Sequence Alignments

Multiple sequence alignments for each gene

Draft Genome CDS

Sequences extracted from various draft genomes used in the selection analyses

Data from: Shifts in selective pressures on snake phototransduction genes associated with photoreceptor transmutation and dim-light ancestry

The visual systems of snakes are heavily modified relative to other squamates, a condition often thought to reflect their fossorial origins. Further modifications are seen in caenophidian snakes, where evolutionary transitions between rod and cone photoreceptors, termed photoreceptor transmutations, have occurred in many lineages. Little previous work, however, has focused on the molecular evolutionary underpinnings of these morphological changes. To address this, we sequenced seven snake eye transcriptomes and utilized new whole genome and targeted capture sequencing data. We used this data to analyze gene loss and shifts in selection pressures in phototransduction genes that may be associated with snake evolutionary origins and photoreceptor transmutation. We identified the surprising loss of rhodopsin kinase (GRK1), despite a low degree of gene loss overall and a lack of relaxed selection early during snake evolution. These results provide some of the first evolutionary genomic corroboration for a dim-light ancestor that lacks strong fossorial adaptations. Our results also indicate that snakes with photoreceptor transmutation experienced significantly different selection pressures from other reptiles. Significant positive selection was found primarily in cone-specific genes, but not rod-specific genes, contrary to our expectations. These results reveal potential molecular adaptations associated with photoreceptor transmutation, and also highlight unappreciated functional differences between rod- and cone-specific phototransduction proteins. This intriguing example of snake visual system evolution illustrates how the underlying molecular components of a complex system can be reshaped in response to changing selection pressures

Comparative sequence analyses of rhodopsin and RPE65 reveal patterns of selective constraint across hereditary retinal disease mutations

Retinitis pigmentosa (RP) comprises several heritable diseases that involve photoreceptor, and ultimately retinal, degeneration. Currently, mutations in over 50 genes have known links to RP. Despite advances in clinical characterization, molecular characterization of RP remains challenging due to the heterogeneous nature of causal genes, mutations, and clinical phenotypes. In this study, we compiled large datasets of two important visual genes associated with RP: rhodopsin, which initiates the phototransduction cascade, and the retinoid isomerase RPE65, which regenerates the visual cycle. We used a comparative evolutionary approach to investigate the relationship between interspecific sequence variation and pathogenic mutations that lead to degenerative retinal disease. Using codon-based likelihood methods, we estimated evolutionary rates (d N/d S) across both genes in a phylogenetic context to investigate differences between pathogenic and nonpathogenic amino acid sites. In both genes, disease-associated sites showed significantly lower evolutionary rates compared to nondisease sites, and were more likely to occur in functionally critical areas of the proteins. The nature of the dataset (e.g., vertebrate or mammalian sequences), as well as selection of pathogenic sites, affected the differences observed between pathogenic and nonpathogenic sites. Our results illustrate that these methods can serve as an intermediate step in understanding protein structure and function in a clinical context, particularly in predicting the relative pathogenicity (i.e., functional impact) of point mutations and their downstream phenotypic effects. Extensions of this approach may also contribute to current methods for predicting the deleterious effects of candidate mutations and to the identification of protein regions under strong constraint where we expect pathogenic mutations to occur.This work was supported by an NSERC Discovery Grant and a Foundation for Fighting Blindness seed grant to BSWC, Vision Science Research Fellowships to FEH, RKS, GMC, and AVN, and Ontario Graduate Fellowships to FEH and RKS. Amir Sabouhanian provided valuable feedback on the statistical analyses

PAML Results Tables

Detailed results and summary tables of each of the PAML analyses

PAML Files

Alignment and tree files needed to run the PAML analyses

FDR Control Results

Results of analyses to control the false discovery rate (FDR)