The opsin repertoire of Jenynsia onca: a new perspective on gene duplication and divergence in livebearers
BMC Research Notes volume 2, Article number: 159 (2009)
Jenynsia onca, commonly known as the one sided livebearer, is a member of the family Anablepidae. The opsin gene repertoires of J. onca's close relatives, the four-eyed fish (Anableps anableps) and the guppy (Poecilia reticulata), have been characterized and each found to include one unique LWS opsin. Currently, the relationship among LWS paralogs and orthologs in these species are unclear, making it difficult to test the hypotheses that link vision to morphology or life history traits. The phylogenetic signal appears to have been disrupted by gene conversion. Here we have sequenced the opsin genes of J. onca in order to resolve these relationships.
We identified nine visual opsins; LWS S180r, LWS S180, LWS P180, SWS1, SWS2A, SWS2B, RH1, RH2-1, and RH2-2. Key site analysis revealed only one unique haplotype, RH2-2, although this is unlikely to shift λmax significantly. LWS P180 was found to be a product of a gene conversion event with LWS S180, followed by convergence to a proline residue at the 180 site.
Jenynsia onca has at least 9 visual opsins: three LWS, one RH1, two RH2, one SWS1 and two SWS2. The presence of LWS P180 moves the location of the LWS P180-S180 tandem duplication event back to the base of the Poeciliidae-Anablepidae clade, expanding the number of species possessing this unusual blue shifted LWS opsin. The presence of the LWS P180 gene also confirms that gene conversion events have homogenized opsin paralogs in fish, just as they have in humans.
Jenynsia is the sister group to the genus Anableps, together they form the subfamily Anablepinae; the genus Oxyzygonectes is the sister group of the subfamily Anablepinae and jointly they form the family Anablepidae. The genus Jenynsia is comprised of thirteen species . Jenynsia onca is distributed in freshwater lakes and rivers throughout southern Brazil . J. onca is a pelagic freshwater fish, whose eyes are morphologically normal in appearance, unlike A. anableps a member of it's sister genus . Distinguishing features of a male J. onca include a tubular unscaled gonopodium, which is either dextral or sinistral. Females of this species correspondingly have either a dextral or sinistral genital opening and mate only with complementary sided males. Additionally, J. onca has dark circular spots on the ventral portion of the flank [1, 2]. Here we have used PCR to characterize the opsin genes of this species.
Vision is an interesting and dynamic sensory modality, particularly in teleost fish, a group which possess some of the greatest morphological and habitat diversity of any animal group. The first step of vision is light absorption, which occurs via opsins expressed in the rods and cones. There are five sub families of vertebrate opsins each with their maximal absorption (λmax) focused on a different area of the visual spectrum : RH1 a rod specific class (490-500 nm), LWS a long-wave sensitive class (490-570 nm), RH2 a middle-wave sensitive class (480-533 nm), SWS1 and SWS2 short-wave sensitive classes (355-440 nm and 410-490 nm respectively) . Within these five classes of conserved opsin proteins, there are variable amino acid residues that give each opsin its unique spectral sensitivity [4–10]. These residues are termed key sites, they are often found at positions in which contact is made with the chromophore, and each has a disproportionate effect on the λmax [5, 11].
Opsin repertoires in fish are particularly interesting due to the extensive pattern of opsin gene duplication and divergence found throughout teleostei. Furthermore, the expanded repertoire is often correlated to differential expression, both spatially across the retina and developmentally over time. For example, in zebrafish (Danio rerio) it has been demonstrated that LWS duplicates localize differentially across the topology of the retina, and that this pattern of expression also changes during development . This has been hypothesized to be a response to the heterogeneous environment the fish live in, with spectral properties differing between the light that hits the dorsal region of the retina and that which hits the ventral region. Another example of differential expression is in cichlids, in this system a subset of the available opsin repertoire is used by each particular species to tune sensitivity in a habitat dependent way . The first step to characterizing these fascinating expressional patterns is through the elucidation of the opsin repertoire itself.
Particularly interesting are the opsin repertoires of the livebearers; ongoing studies of opsin gene duplication and divergence in guppies (Poecilia reticulata) and the four-eyed fish (Anableps anableps) have shown that both species have an expanded LWS repertoire. Both species have recent species-specific duplicates and a repertoire of ten visual opsins [14–16]. Interestingly both also have remarkable morphology, in the case of the guppy this is a variable male pigmentation pattern, while A. anableps have unusual four-eyed morphology. Based on its phylogenetic position, J. onca functions as a useful out-group for comparison of opsin repertoire, particularly between A. anableps and the Poeciliids . Out-groups are used to identify synapomorphies (shared derived characters) and can indicate whether one apomorphic trait evolved before another. Here our goal is to determine whether or not particular opsin genes and gene sequences are associated with coloration in guppies and the four-eye morphology in A. anableps.
Results and Discussion
PCR screening using gene specific primers (Table 1) identified nine visual opsins: LWS S180r, LWS S180, LWS P180, SWS1, SWS2A, SWS2B, RH1, RH2-1, & RH2-2. All opsins except LWS S180r and S180 are expressed, having been amplified from eye cDNA derived from one adult male and one adult female J. onca. LWS S180r and S180 opsins were amplified only from genomic DNA. However the reason that LWS S180 and S180r could not be retrieved from cDNA could be attributed to life stage, as only adults were used in this study. A particularly interesting finding in this repertoire was the LWS P180 opsin, as it's presence moves the point of LWS duplication farther back within the Cyprinodontiformes order. We did not find any species-specific gene duplication events in J. onca.
The J. onca opsins were aligned to orthologous sequences from other fish species. Sequences in the alignment were 573 to 930 base pairs (bp) long. We used PAUP* 4.0B10 to calculate genetic distances based on the modeltest best-fit model of sequence evolution and to reconstruct Neighbour joining (NJ) (Figure 1) and Maximum likelihood (ML) (see Additional file 1) trees [18–21]. Sequences from each opsin subfamily formed well-supported monophyletic groups, with bootstrap support ≥ 97% (1000 replicates) .
The phylogenetic analysis indicated that gene identity was consistent with regard to both the subfamily that the genes fell into and species taxonomy. The LWS subfamily can be further grouped into haplotype clades, however, J. onca LWS S180 and LWS P180 did not fall out on the tree where they would be predicted to based on haplotype identity (Figure 1). LWS P180 was found to be highly similar to A. anableps LWS S180γ and Poeciliidae LWS P180 only over its 3' end (Figure 2). When only the 3' region of these two genes and the Poeciliidae LWS are used in phylogenetic analysis two distinct clades with 60% bootstrap are observed separating LWS S180 from LWS P180 (including A. anableps LWS S180γ) (Figure 3).
We hypothesize that two recent gene conversion events explain the observed pattern. One event would have occurred between the two J. onca genes, from LWS S180 to LWS P180 leaving only the original 3' end of LWS P180, masking the true orthologous relationship. This conversion event in J. onca also explains why the sequence surrounding the proline at amino acid postion 180, which occurs before the 3' end, is very similar to the J. onca LWS S180. Indeed, based on the pattern of similarity it seems likely that the P180 substitution occurred independently in the Anablepid lineage. The second conversion event we suggest occurred recently in the A. anableps LWS S180γ, with the majority of the 5' region of this gene having been over-written by LWS S180α. This event explains the grouping of LWS S180γ with Poeciliid LWS P180 in the 3' phylogenetic tree (Figure 3) and why LWS S180α and LWS S180γ are nearly identical until the last portion of the 3' end. This conversion event also clarifies why the A. anableps LWS S180γ does not have a proline at the 180 amino acid site. Our hypothesis that gene conversion has influenced sequence variation among opsin paralogs is supported by an analysis using GENECONV 1.81A, which detects sequence pairs that have abnormally long continuous regions of high sequence similarity found within regions of lower similarity overall [23, 24]. GENECONV detected a gene conversion event between LWS S180 and LWS P180/S180γ in A. anableps from the beginning of the sequence to 629 bp (p = 0.012) and in J. onca from the start of the sequence to 567 bp (p = 0.036), which correspond to our predicted conversion events.
The only alternative explanation for this pattern of evolution would be convergent evolution of this 3' region (Figure 2), however given the previous observations of gene conversion in opsins it is not the most likely explanation. Homogenized key site haplotypes have been observed not only in fish, such as the guppy where the LWS A180 is a product of duplication followed by partial gene conversion, but also in humans where conversion events are often detrimental [14, 25, 26]. The conversion events we propose in livebearers may have been facilitated between LWS P180 and LWS S180 due to their position as tandem inverted duplicates, something that has been confirmed in both Poecilia and Xiphophorous .
As mentioned above, there are key sites within each opsin subclass, which may be used to estimate the opsin λmax. We investigated the key site haplotypes of these nine opsins and based on the LWS haplotypes we have estimated the LWS λmax values  [see Additional file 2]. The key site substitutions seen in the LWS P180 although found in other species, confer a significant change in λmax; the serine to proline substitution alone results in a -19 nm shift . Only one key site substitution is not seen in J. onca's relatives; RH2-2 deviates from the isoleucine consensus residue to a valine residue in amino acid site 65 (based on human LWS amino acid numbering ), although this would likely not significantly shift the λmax . The most interesting point to glean from the haplotype comparison is that despite significant nucleotide divergence between J. onca and its relatives there is a large degree of amino acid conservation at key sites between orthologs.
Jenynsia onca has nine visual opsins: three LWS, one RH1, two RH2, one SWS1 and two SWS2. Despite nucleotide level divergence between related orthologs it appears there is significant phenotypic (λmax) conservation with only one instance of amino acid key site residue substitution. J. onca in the future could help us identify differences in opsin expression that are associated with the unusual eyes of A. anableps and with the remarkable coloration of guppies by acting as an out-group.
PCR primers were designed to amplify nine opsins from the five visual opsins subfamilies (Table 1). These primers were complementary to regions in each opsin gene or subfamily that were conserved in guppies (Poecilia reticulata), and A. anableps. Two primer pairs were engaged for each gene.
Each primer pair was used to survey cDNA or genomic DNA in PCR reactions using Bio-Rad iProof high-fidelity DNA polymerase in an Eppendorf™ Mastercycler® EP Grad S thermocycler using the following conditions: Initial denaturation at 98°C for 30 s, 35 cycles with denaturation at 98°C for 5 s, annealing at 55 - 65°C (in 5°C intervals) for 12 s, extension at 72°C for 25 s and a final extension at 72°C for 5 min, additional primers (1 μl at 10 mM) were added, at the beginning of the last PCR cycle to prevent heteroduplex formation.
Amplicons of the predicted size were excised using QIAquick® Gel Extraction Kit or purified using QIAquick® PCR Purification Kit. Purified products were A-tailed using Invitrogen™ Taq polymerase and cloned using the Promega® pGEM™ - T Easy Vector System II kit. Clones containing inserts of the correct size were sequenced using labelled M13 forward and reverse primers and a Licor sequencer at the University of Victoria Centre for Biomedical Research.
Live J. onca were obtained from a commercial supplier (The Afishionados, Winnipeg, Manitoba, Canada). Two adult J. onca (one male and one female) were euthanized in buffered MS222. Total RNA was isolated from the eyes using Aurum™ Total RNA Fatty and Fibrous Tissue Pack, immediately after euthanasia and enucleation cDNA was synthesized using BioRad® iScript Select cDNA Synthesis Kit from total RNA. DNA was isolated from the fish carcass using QIAquick® DNeasy Blood & Tissue Kit.
Two phylogenetic trees were reconstructed for the complete set of opsin sequences. The partial coding sequence tree included sequences from Jenynsia onca, Anableps anableps, Poecilia reticulata, Xiphophorus pygmaeus, Lucania goodei, Oryzias latipes and Danio rerio, Poecilia picta, Poecilia parae, Poecilia bifurca and Tomeurus gracilis sequence files used were 412 to 819 bp long. The second tree was based on 243 bp of LWS 3' coding sequence from J. onca, A. anableps, P. reticulata, X. pygmaeus, P. picta, P. parae, P. bifurca and T. gracilis [Accession numbers see Additional file 3]. The aligned sequences for the partial coding sequence phylogenetic tree were first used to obtain the best-fit model of evolution using Modeltest . The phylogenetic reconstruction was done using ML and NJ (1000 bootstrap reanalyses) in PAUP* 4.8B10 and utilized the optimal model parameters [18, 20–22]. The root of the partial coding sequence tree was positioned along the branch separating the LWS opsins from all others . The 3' coding tree was constructed using MEGA4 using the Jukes-Cantor algorithm, NJ, and support for nodes were estimated using 500 bootstrap reanalyses [20, 22, 31, 32]. Pair-wise deletion was used in the case of missing nucleotides for the analysis.
Gene conversion detection was undertaken using GENECONV version 1.81A . We used the program's default values with the exception of gscale, which we set to 2. The input for analysis was a coding sequence nucleotide alignment containing A. anableps LWS S180α and LWS S180γ as well as J. onca, X. pygmaeus, P. picta, P. parae, P. reticulata, and P. bifurca LWS S180 and LWS P180.
Lucinda PHF, Reis RE, Quevedo R: Jenynsia onca, a new species of anablepid fish (Teleostei: Cyprinodontiformes) from southern Brazil and its phylogenetic position. Ichtyological Exploratoin of Freshwaters. 2002, 13: 33-40.
Nelson JS: Fishes of the World. 1984, New York: John Wiley & Sons
Applebury ML: Relationships of G-protein-coupled receptors. A survey with the photoreceptor opsin subfamily. Soc Gen Physiol Ser. 1994, 49: 235-248.
Yokoyama S: Molecular evolution of vertebrate visual pigments. Prog Retin Eye Res. 2000, 19: 385-419. 10.1016/S1350-9462(00)00002-1.
Yokoyama S, Radlwimmer FB: The molecular genetics of red and green color vision in mammals. Genetics. 1999, 153: 919-932.
Yokoyama S, Takenaka N, Blow N: A novel spectral tuning in the short wavelength-sensitive (SWS1 and SWS2) pigments of bluefin killifish (Lucania goodei). Gene. 2007, 396: 196-202. 10.1016/j.gene.2007.03.019.
Chinen A, Matsumoto Y, Kawamura S: Reconstitution of ancestral green visual pigments of zebrafish and molecular mechanism of their spectral differentiation. Mol Biol Evol. 2005, 22: 1001-1010. 10.1093/molbev/msi086.
Takenaka N, Yokoyama S: Mechanisms of spectral tuning in the RH2 pigments of Tokay gecko and American chameleon. Gene. 2007, 399: 26-32. 10.1016/j.gene.2007.04.036.
Yokoyama S: Molecular evolution of color vision in vertebrates. Gene. 2002, 300: 69-78. 10.1016/S0378-1119(02)00845-4.
Yokoyama S, Tada T, Zhang H, Britt L: Elucidation of phenotypic adaptations: Molecular analyses of dim-light vision proteins in vertebrates. Proc Natl Acad Sci USA. 2008, 105: 13480-13485. 10.1073/pnas.0802426105.
Yokoyama S: Amino acid replacements and wavelength absorption of visual pigments in vertebrates. Mol Biol Evol. 1995, 12: 53-61.
Takechi M, Kawamura S: Temporal and spatial changes in the expression pattern of multiple red and green subtype opsin genes during zebrafish development. J Exp Biol. 2005, 208: 1337-1345. 10.1242/jeb.01532.
Parry JW, Carleton KL, Spady T, Carboo A, Hunt DM, Bowmaker JK: Mix and match color vision: tuning spectral sensitivity by differential opsin gene expression in Lake Malawi cichlids. Curr Biol. 2005, 15: 1734-1739. 10.1016/j.cub.2005.08.010.
Ward MN, Churcher AM, Dick KJ, Laver CR, Owens GL, Polack MD, Ward PR, Breden F, Taylor JS: The molecular basis of color vision in colorful fish: four long wave-sensitive (LWS) opsins in guppies (Poecilia reticulata) are defined by amino acid substitutions at key functional sites. BMC Evol Biol. 2008, 8: 210-10.1186/1471-2148-8-210.
Owens GL, Windsor DJ, Mui J, Taylor JS: A fish eye out of water: ten visual opsins in the four-eyed fish, Anableps anableps. PLoS ONE. 2009, 4: e5970-10.1371/journal.pone.0005970.
Hoffmann M, Tripathi N, Henz SR, Lindholm AK, Weigel D, Breden F, Dreyer C: Opsin gene duplication and diversification in the guppy, a model for sexual selection. Proc Biol Sci. 2007, 274: 33-42. 10.1098/rspb.2006.3707.
Hrbek T, Seckinger J, Meyer A: A phylogenetic and biogeographic perspective on the evolution of poeciliid fishes. Mol Phylogenet Evol. 2007, 43: 986-998. 10.1016/j.ympev.2006.06.009.
Swofford D: PAUP* 4b10. Phylogenetic analysis using parsimony (* and other methods). 1996, Sinauer, Sunderland, Massachusetts
Posada D, Crandall KA: MODELTEST: testing the model of DNA substitution. Bioinformatics. 1998, 14: 817-818. 10.1093/bioinformatics/14.9.817.
Saitou N, Nei M: The neighbor-joining method: a new method for reconstructing phylogenetic trees. Mol Biol Evol. 1987, 4: 406-425.
Felsenstein J: Evolutionary trees from DNA sequences: a maximum likelihood approach. J Mol Evol. 1981, 17: 368-376. 10.1007/BF01734359.
Felsenstein J: Confidence limits on phylogenies: an approach using the bootstrap. Evolution. 1985, 39: 783-791. 10.2307/2408678.
Sawyer SA: GENECONV: a computer package for the statistical detection of gene conversion. 1999, Department of Mathematics, Washington University, St. Louis, Mo
Li WH: Unbiased estimation of the rates of synonymous and nonsynonymous substitution. J Mol Evol. 1993, 36: 96-99. 10.1007/BF02407308.
Winderickx J, Battisti L, Hibiya Y, Motulsky AG, Deeb SS: Haplotype diversity in the human red and green opsin genes: evidence for frequent sequence exchange in exon 3. Hum Mol Genet. 1993, 2: 1413-1421. 10.1093/hmg/2.9.1413.
Reyniers E, Van Thienen MN, Meire F, De Boulle K, Devries K, Kestelijn P, Willems PJ: Gene conversion between red and defective green opsin gene in blue cone monochromacy. Genomics. 1995, 29: 323-328. 10.1006/geno.1995.9998.
Yokoyama S: Evolution of dim-light and color vision pigments. Annu Rev Genomics Hum Genet. 2008, 9: 259-282. 10.1146/annurev.genom.9.081307.164228.
Davies WL, Collin SP, Hunt DM: Adaptive gene loss reflects differences in the visual ecology of basal vertebrates. Mol Biol Evol. 2009, 8: 1803-9. 10.1093/molbev/msp089.
Nathans J, Thomas D, Hogness DS: Molecular genetics of human color vision: the genes encoding blue, green, and red pigments. Science. 1986, 232: 193-202. 10.1126/science.2937147.
Okano T, Kojima D, Fukada Y, Shichida Y, Yoshizawa T: Primary structures of chicken cone visual pigments: vertebrate rhodopsins have evolved out of cone visual pigments. Proc Natl Acad Sci USA. 1992, 89: 5932-5936. 10.1073/pnas.89.13.5932.
Tamura K, Dudley J, Nei M, Kumar S: MEGA4: Molecular Evolutionary Genetics Analysis (MEGA) software version 4.0. Mol Biol Evol. 2007, 24: 1596-1599. 10.1093/molbev/msm092.
Jukes TH, Cantor CR: Evolution of protein molecules. Mammalian protein metabolism. 1969, 3: 21-132.
This research was carried out in Dr. John S. Taylor's lab at the University of Victoria. We thank Dr. Taylor for discussion and help editing the manuscript. We also thank Spencer Jack from Afishionados for importing Jenynsia Onca, Roderick Haesevoets at the University of Victoria Centre for Biomedical Research for sequencing assistance and the two anonymous reviewers for their helpful comments. This work was supported by the Natural Sciences and Engineering Research Council of Canada.
The authors declare that they have no competing interests.
DJW & GLO carried out the PCR survey of cDNA and genomic DNA, completed the sequence alignment, and created the phylogenies. DJW drafted the manuscript. Both authors read and approved the final manuscript.
Electronic supplementary material
Additional file 1: Phylogenetic analysis of J. onca opsins using a Maximum-likelihood approach. A maximum likelihood tree, which uses opsin-coding sequence from J. onca and relatives. PAUP* 4.0B10 was used to estimate genetic distances, based on modeltest's best-fit model of evolution, and complete phylogenetic analysis [18, 19] [accession numbers see Additional file 3]. All codon positions were used and pair-wise deletion was used in the case of missing nucleotides for the analysis. (JPEG 86 KB)
Additional file 2: Amino acid key site haplotypes and LWS lambda max prediction for Jenynsia Onca and relatives. Key site haplotypes for J. onca, A. anableps, P. reticulata and L. goodei [Accession numbers see Additional file 3]. Amino acid site numbering is based on human LWS . Lambda max values estimated from mutational analysis of other LWS opsins [27, 28]. An 'X' indicates that the sequence obtained in our PCR survey did not encompass the respective site. (PDF 65 KB)
Additional file 3: Accession numbers for all sequences used. Accession numbers of sequences generated from a degenerate PCR survey of J. onca and those used for phylogenetic analysis. (DOC 86 KB)
About this article
Cite this article
Windsor, D.J., Owens, G.L. The opsin repertoire of Jenynsia onca: a new perspective on gene duplication and divergence in livebearers. BMC Res Notes 2, 159 (2009). https://doi.org/10.1186/1756-0500-2-159
- Gene Conversion
- Neighbour Join
- Conversion Event
- Opsin Gene
- Gene Conversion Event