- Research article
- Open Access
Evolution and expression of the phosphodiesterase 6 genes unveils vertebrate novelty to control photosensitivity
BMC Evolutionary Biology volume 16, Article number: 124 (2016)
Phosphodiesterase 6 (PDE6) is a protein complex that hydrolyses cGMP and acts as the effector of the vertebrate phototransduction cascade. The PDE6 holoenzyme consists of catalytic and inhibitory subunits belonging to two unrelated gene families. Rods and cones express distinct genes from both families: PDE6A and PDE6B code for the catalytic and PDE6G the inhibitory subunits in rods while PDE6C codes for the catalytic and PDE6H the inhibitory subunits in cones. We performed phylogenetic and comparative synteny analyses for both gene families in genomes from a broad range of animals. Furthermore, gene expression was investigated in zebrafish.
We found that both gene families expanded from one to three members in the two rounds of genome doubling (2R) that occurred at the base of vertebrate evolution. The PDE6 inhibitory subunit gene family appears to be unique to vertebrates and expanded further after the teleost-specific genome doubling (3R). We also describe a new family member that originated in 2R and has been lost in amniotes, which we have named pde6i. Zebrafish has retained two additional copies of the PDE6 inhibitory subunit genes after 3R that are highly conserved, have high amino acid sequence identity, are coexpressed in the same photoreceptor type as their amniote orthologs and, interestingly, show strikingly different daily oscillation in gene expression levels.
Together, these data suggest specialisation related to the adaptation to different light intensities during the day-night cycle, most likely maintaining the regulatory function of the PDE inhibitory subunits in the phototransduction cascade.
Phosphodiesterases (PDEs) are a large family of proteins with 21 identified human members divided into 11 subfamilies: PDE1-PDE11 . The present work focuses on the photoreceptor cell-specific PDE6 subfamily, which consist of α, β and α’ catalytic subunits encoded by the PDE6A, PDE6B and PDE6C genes. PDE6 catalytic activity is regulated by the γ and γ’ inhibitory subunits encoded by PDE6G and PDE6H, respectively .
Phylogenetic analyses complemented with chromosomal analyses in human previously published by our lab [3, 4] and phylogenetic data from other investigators  suggested that the PDE6 catalytic gene family with the PDE6A, PDE6B and PDE6C genes expanded in the two rounds of whole genome duplications (2R) that occurred early in vertebrate evolution [6, 7]. Other analyses of phylogeny and exon-intron organisation have shown that the PDE6 subfamily is most closely related to the PDE5 and PDE11 subfamilies [5, 8]. 2R was also proposed to explain the duplicates found in the PDE6 inhibitory subunit gene family, whose genes were found to be located in the same paralogon (set of related chromosomes) as the developmentally important homeobox gene clusters  that are known to have been duplicated in 2R .
PDE6 proteins are expressed in vertebrate rod and cone photoreceptor cells. As for the other components of the phototransduction cascade (see  for references), rods and cones use distinct but related PDE6 subunits: rods express the PDE6A and PDE6B genes, which give rise to a catalytic heterodimer, and the PDE6G inhibitory subunit gene, whereas cones express PDE6C, resulting in a catalytic homodimer, and the PDE6H inhibitory subunit gene [2, 4].
The PDE6 catalytic subunit proteins have two GAF domains (GAF domains are named after the proteins that contain them; cGMP-activated PDEs, adenylyl cyclases and FhlA) followed by one catalytic domain, a structure that is shared with the PDE2, PDE5, PDE10 and PDE11 subfamilies . The catalytic domain is present in all of the different PDE subfamilies, with differences in substrate specificity . The PDE6 enzymes are distinguished from the other PDE subfamilies by a much higher catalytic activity and two accessory inhibitory subunits that interact with a GAF domain and the catalytic domain of the catalytic subunits and thus block activity during dark conditions [1, 10]. The emergence of the inhibitory subunits has been proposed as one of the events that made it possible for a higher catalytic rate to evolve in PDE6 compared to the other PDEs, resulting in a fast photoreceptor response .
The PDE6 holoenzyme function can be summarised as follows. The cascade is initiated by an opsin, which is activated by a photon. There are different opsins with distinct spectral selectivity mediating colour vision in various types of cones and a rhodopsin mediating dim-light vision in rods. Activated opsin acts as a G nucleotide exchange factor for the heterotrimeric G-protein transducin. A GTP molecule replaces GDP at the active site of the alpha subunit of transducin leading to dissociation of the transducin heterotrimer into the activated alpha subunit and a heterodimer of the beta and gamma subunits. The alpha transducin then activates the PDE6, a cGMP phosphodiesterase. Activation takes place when transducin alpha subunits remove the two PDE6 inhibitory subunits . The activated PDE6 hydrolyses cGMP into GMP, which reduces the cGMP levels in the cell and leads to a closure of cyclic nucleotide-gated channels and hyperpolarisation of the photoreceptor cell. The regulation of activity by the PDE6 inhibitory subunits involves specific regions within these small proteins responsible for the binding to the PDE6 catalytic subunits and the alpha subunit of transducin but also for the stabilisation of the binding affinity, the stimulation of non-catalytic cGMP binding, the increase in cGMP exchange at the GAF domains and the stimulation of GTPase activity of transducin, in complex with RGS9-1 .
In this study, an extensive repertoire of vertebrate species was used to resolve the evolutionary history of the PDE6 subunit gene families and their chromosomal regions more precisely. This gives us a robust dataset, which demonstrates the expansion of both the catalytic and inhibitory subunit gene families in 2R as well as further duplication of the inhibitory subunit gene family in the teleost-specific whole genome duplication (3R) as well as through local duplications. Previous studies have described the same set of catalytic subunit gene repertoire in zebrafish as in human: pde6a, pde6b and pde6c [13, 14]. Here we demonstrate that zebrafish has retained, as a result of 3R, two paralogous genes for each of the two PDE6 inhibitory genes found in human: pde6ga, pde6gb, pde6ha and pde6hb. Additionally, we also show that zebrafish has retained an extra PDE6 inhibitory subunit gene from 2R that we named pde6i which is not present in amniotes. Due to the importance of zebrafish as a model for both evolution and visual function and disorders, we have analysed the expression of all PDE6 genes in zebrafish and found striking differences related to the adaptation to different light intensities during the day-night cycle.
Three PDE6 catalytic subunit genes expanded in 2R and were subsequently retained in most vertebrate lineages, with minimal to no expansion in 3R
An alignment of the PDE6 catalytic subunit amino acid sequences used in the phylogenetic maximum likelihood (PhyML) analysis is provided in Additional file 1 and the sequence identifiers and the genomic locations of the included sequences are provided in Additional file 2. The PhyML tree, rooted with the human PDE5A and PDE11A genes, of this family forms three well-supported clades with each of the human PDE6A, PDE6B and PDE6C genes and their orthologues (Fig. 1). All vertebrate genomes investigated have three PDE6 catalytic subunit genes, with a few exceptions. The first exception is the green spotted pufferfish, which has two PDE6C paralogs; one located on chromosome 17, that clusters basally to the other teleost sequences, and the other on chromosome 2, that clusters together with the Nile tilapia, three-spined stickleback and medaka sequences. Interestingly, two PDE6C genes were also found in the fugu genome (data not shown), indicating either a pufferfish-specific duplication (with rapid evolution of the duplicate) or that the pufferfishes have retained 3R duplicates that have been lost in other teleost fish, a scenario supported by the phylogenetic analysis (Fig. 1). The second exception is the absence of the PDE6A gene in birds and non-avian reptiles (Additional file 3: Figure S1).
We identified a single gene in the arctic lamprey (Lethenteron camtschaticum) genome assembly (data not shown). This gene has previously been identified in the sea lamprey , indicating that northern hemisphere lampreys have retained only one catalytic subunit gene. The sea lamprey sequence clusters basally to all other PDE6 catalytic subunit clusters in the phylogenetic analysis (Fig. 1), thus we cannot confidently assign orthology to any of the gnathostome sequences. The invertebrate chordate Ciona intestinalis has one gene that has been assigned the same Ensembl protein family ID as the vertebrate PDE6 catalytic subunit genes. This sequence was added to the analysis for relative dating of the duplications of the vertebrate genes and it clusters basally to the three vertebrate clusters, suggesting that it is ortholog to the ancestor of the vertebrate PDE6A, PDE6B and PDE6C genes. This is in agreement with the previously suggested expansion of this gene family in 2R, based on much fewer sequences [3, 4].
The chromosomal region housing the PDE6 catalytic subunit genes belongs to a paralogon that arose in 2R
The human PDE6A, PDE6B and PDE6C genes are located on chromosomes 5, 4 and 10, respectively, like the well-studied neuropeptide Y receptor genes and many neighbouring gene families. These have been previously demonstrated to have expanded in the early vertebrate tetraploidisations . As the PDE6 catalytic subunit genes are located some distance away from those reported previously, we analysed additional families in these regions to see whether the duplications could have taken place in the timeframe of 2R.
A total of eleven neighbouring gene families were identified as having members in the same chromosomal regions as the three PDE6 catalytic subunit genes (Additional file 3: Table S1 and Additional file 2). Three of the neighbouring gene families (CP, SLC26A and ZNF) were excluded at an initial stage of the analysis due to their complex phylogenetic tree topology or their multitude of members. For the remaining eight families four have annotated invertebrate and vertebrate members, DPYS, PPP2R2, STK32 and TBC1D; while the other four, ABLIM, AFAP, JAKMIP and SH3TC, only have vertebrate members. The sequence identifiers and genomic locations of the included neighbouring gene families are provided in Additional file 2.
The phylogenetic trees and the species representation of the STK32 gene family support an expansion in 2R giving rise to three vertebrate genes (Additional file 3: Figure S2). However, teleost fish appear to have lost the STK32B. Early origin of this gene is shown by its presence outside tetrapods, in the genomes of the coelacanth and the spotted gar. The phylogenetic tree of the TBC1D family shows two clusters consistent with duplication in the same timeframe as 2R (Additional file 3: Figure S3). The PPP2R2 phylogenetic tree supports an expansion in 2R and a possible expansion of the PPP2R2A genes in the teleost specific 3R event (Additional file 3: Figure S4). The phylogenetic tree of the DPYS family shows six gene clusters, four of which, namely CRMP1, DPYSL2, DPYSL3 and DPYSL4, seem to form a vertebrate subfamily clustering with a basal putative Ciona intestinalis ortholog (Additional file 3: Figure S5). This is consistent with an expansion in 2R. The topology also supports a duplication of the teleost DPYSL2 genes in 3R. The two other gene clusters that include the human DPYS and DPYSL5 genes with their vertebrate homologs seem to have branched off from the previously mentioned vertebrate subfamily before the divergence of the lineages leading to tunicates and vertebrates.
The phylogenies of the ABLIM, AFAP, JAKMIP and SH3TC families show that these families expanded before the separation of lobe-finned and ray-finned vertebrates (Additional file 3: Figures S6-S9). However, if this expansion occurred due to 2R is not clear due to the lack of invertebrate family members for relative dating. The ABLIM family shows evidence for a duplication of both the ABLIM1 gene and the ABLIM3 gene in 3R (Additional file 3: Figure S6).
A summary of the localisation of the identified neighbouring gene family genes as well as the PDE6 catalytic subunit genes on human, chicken, spotted gar and zebrafish chromosomes are presented in Fig. 2.
An ancestral PDE6 inhibitory subunit gene emerged in the vertebrate ancestor and duplicated in 2R and 3R
An amino acid sequence alignment of the identified PDE6 inhibitory subunit genes is provided in Additional file 4 and the sequence identifiers and genomic locations for the included sequences are listed in Additional file 2. Due to the relatively high level of sequence conservation, 72–81 % between lampreys and humans (data not shown), ≥64 % between human and zebrafish (see Table 1) and short sequence lengths, no reliable phylogenetic signal can be obtained for the PDE6 inhibitory subunit genes. Therefore, their chromosomal locations were thoroughly investigated and used for assigning orthology. As a result, we discovered that this gene family has three paralogs (PDE6G, PDE6H and PDE6I) that are located in the same chromosomal regions as the somatostatin receptor 2, 3 and 5 genes (SSTR2, SSTR3 and SSTR5) and the urotensin receptor genes, previously shown to reside in a paralogon [16, 17]. Thus, now when more genomic resources are available we can show that the PDE6 inhibitory subunit genes are located in a different paralogon than the previously suggested HOX paralogon .
Previously, only PDE6G and PDE6H have been described in different vertebrate species, thus our identification of PDE6I reveals that this family has a third member (Fig. 3). The three PDE6 inhibitory subunit genes were found on linkage groups (LG) 10, 12 and 13 of the spotted gar genome, which also carries SSTR2, SSTR3 and SSTR5 genes, respectively. A recent analysis of the eye transcriptome of the Florida gar (Lepisosteus platyrhincus) showed that the PDE6I is expressed (T.D. Lamb et al., personal communication).
The coelacanth has three different PDE6 inhibitory subunit genes in its genome assembly, one fragmented (of unknown orthology, located on scaffold JH132110.1), one full length (PDE6G located on scaffold JH126581.1) and one that appears to be a pseudogene due to a frame-shift mutation in the first exon disrupting the translation (PDE6I located on scaffold JH127264.1) (Fig. 3 and Additional file 3: Figure S10). The Western clawed frog has four full-length and one fragmented PDE6 inhibitory subunit genes all located on different genomic scaffolds. We could confidently assign the identity as PDE6G and PDE6I for two of the frog genes (located on scaffolds GL172940.1 and GL172810.1 respectively), whereas the others did not share synteny with any of the spotted gar chromosomal regions carrying a PDE6 inhibitory subunit gene (Fig. 3 and Additional file 3: Figure S10).
We identified orthologs for PDE6G and PDE6H in the teleost species investigated but only zebrafish and Mexican cave tetra has retained PDE6I (Additional file 3: Figures S11-S12 and Additional file 5). Most of the teleost species investigated have retained putative 3R duplicates for PDE6H, while only zebrafish and medaka have retained 3R paralogs for PDE6G (pde6ga and pde6gb) (Fig. 3). The teleost PDE6H genes have experienced intron loss at different time points during teleost evolution. The zebrafish and Mexican cave tetra have lost the first intron of the coding region in the pde6ha gene (located on chromosome 6 in zebrafish), thus it probably took place in their common ancestor. The other investigated teleosts have lost the last intron of the coding region in their PDE6H genes (see Fig. 3). These intron losses along with the conservation of synteny allowed assignment of orthology of the teleost PDE6 inhibitory subunit genes. For the complete repertoire of the PDE6 inhibitory subunit genes identified in teleost genomes see Fig. 3. For synteny comparisons of teleost scaffolds or chromosomes with PDE6 inhibitory subunit genes to spotted gar linkage groups, see Additional file 3: Figures S11-S17.
The PDE6I gene could not be identified in the elephant shark (tentatively only PDE6G and PDE6H) or in amniotes. In chicken, the PDE6G and PDE6H genes are located close to the SSTR2 gene on chromosome 18 and SSTR3 on chromosome 1, respectively. The human PDE6G gene is located on chromosome 17 as is the SSTR2 gene. In contrast, the human PDE6H gene on chromosome 12 is not located on the same chromosome as the SSTR3 gene, which is located on chromosome 22 (Additional file 5). However, the opossum PDE6H gene is located on the same chromosome as the SSTR3 gene, chromosome 8, which suggests that a translocation of PDE6H took place in the primate lineage.
Lampreys are a special case. Two PDE6 inhibitory subunit sequences have previously been cloned in the sea lamprey: GenBank accession number ABO16480.2 and NCBI accession number ABO64650.1 . Due to the lack of synteny data we could not assign their orthology to the other vertebrate genes (Fig. 3). In the arctic lamprey genome we could only identify sequence fragments: one first exon, one middle exon and three last exons, all located on different scaffolds or contigs that prevented us to assign any orthology.
In order to analyse the evolutionary origin of these genes, we performed TBLASTN searches in all invertebrate genomes present in the NCBI whole genome shotgun reference database, using the human PDE6G and PDE6H amino acid sequences as templates. We also performed nhmmer searches in the invertebrate genomes present in Table 2. Despite these efforts, not a single homologous sequence for the PDE6 inhibitory subunit genes was identified in invertebrates.
Comparison of zebrafish PDE6 inhibitory subunit paralog sequences reveal highly conserved regions
Using zebrafish as a model, we investigated the possible specialisations of the PDE6 subunit genes; pde6a, pde6b and pde6c for the catalytic subunits and the PDE6 inhibitory subunit genes retained after 3R; pde6ga, pde6gb, pde6ha and pde6hb. A comparison of the amino acid sequence identity between human and zebrafish for the PDE6 inhibitory subunits confirmed high level of conservation: human PDE6G shares 75 % and 93 % identity with zebrafish pde6ga and pde6gb, respectively; while human PDE6H displays approximately 70 % identity with both zebrafish pde6ha and pde6hb (Table 1).
Additionally, an amino acid sequence alignment of the human and zebrafish PDE6 inhibitory subunit sequences showed some differences in each gene, but the regions that are known to be involved in specific functions, previously described in , are highly conserved (Fig. 4). Some exceptions were observed, like an asparagine instead of threonine at position 65 in pde6ga, a phenylalanine instead of a tyrosine at position 84 in human PDE6H and one alteration that could have a functional effect in the pde6i carrying a glutamic acid instead of the lysine at position 29 (Fig. 4). These results, together with the gene expression data provided in the sections that follow, suggest that the specialisations did not involve any neofunctionalisation at the protein structure level.
The percentage of identity for the 3′UTR sequences of the zebrafish PDE6 inhibitory subunit genes were analysed in order to find regions suitable for use in the design of probes able to distinguish between the 3R duplicates. We found that, pde6ga and pde6gb share 48 % nucleotide sequence identity, pde6ha and pde6hb share 59 %, and pde6i shares between 48–52 % identity with the others (Additional file 3: Table S2).
The zebrafish PDE6 inhibitory subunit gene paralogs are expressed in the same cell types as their amniote orthologs
The expression pattern of the three PDE6 catalytic subunit genes present in zebrafish, has been previously reported . In this investigation, we have analysed the expression levels for all zebrafish PDE6 genes by qRT-PCR, as well as the expression pattern of the PDE6 inhibitory subunit genes; pde6ga, pde6gb, pde6ha, pde6hb and pde6i, by in situ hybridisation (ISH). Pde6i was excluded from the ISH experiments due to the extremely low expression levels found in adult eyes (see next section).
There is a strikingly different morphology of the outer retina between day and night due to the retinomotor movements, previously described in several species, including zebrafish [18, 19]. Interestingly, both the intensity and the location of the staining observed for the four studied genes were substantially different in retinae from animals sacrificed in the afternoon (17:00, Fig. 5a-d´´) versus the middle of the night (4:00, Fig. 5e-h´´). Due to this, the photoreceptor naming will be done only in the day retinae (Fig. 5a-d´´). ISH experiments revealed staining only in photoreceptor cells, throughout the entire retina. Expression of pde6ga (Fig. 5a-a´´) and pde6gb (Fig. 5b-b´´) were observed in the myoid and cytoplasm surrounding the nuclei of rods. Staining in the SSC seems to be present in A´-A´´ and B´-B´´ and we cannot discard that pde6ga and pde6gb are expressed in rods and also SSC. However, if true, this will not be in the entire retina (see Fig. 6). Expression of pde6ha (Fig. 5c-c´´) and pde6hb (Fig. 5d-d´´) was observed in all cones. The yellowish colour observed in 5g´´ could suggest double staining; however its absence in the day retina, combined with the around 100 fold lower expression levels (see next section) and the loose morphology of the night retinae led us to conclude that this is just an artefact. Additionally, double ISH experiments revealed coexpression of pde6ga and pde6gb in rods (Fig. 6a-c) and pde6ha and pde6hb in cones (Fig. 6d-f), throughout the retina. Expression of the four genes was also observed in the pineal complex (Fig. 5i-l) without noticeable differences between day and night (data not shown). However, to reveal any differential expression in the pineal complex, qRT-PCR must be performed.
In addition to the analyses on adults, ontogenetic analyses using 1–6 dpf larvae were performed by whole-mount ISH (WISH). The results reveal a synchronised onset of expression in the pineal complex at 1 dpf for all genes, except pde6ha that starts at 2 dpf, and delayed onset of expression for all of them until around 2–3 dpf in the retina (Additional file 3: Figure S18). The expression in the retina starts in the ventral part, in accordance with the differentiation process that starts in this region in zebrafish, different to most vertebrates, in which it is initiated in central locations .
The zebrafish PDE6 inhibitory subunit paralogs show differential oscillatory patterns in expression levels for both rods and cones in adult retinae
qRT-PCR experiments were performed to analyse the expression levels of all zebrafish PDE6 genes in adult eyes from animals collected at six different time points during 24 h. The results were normalised to the mean of the housekeeping genes (actb1 and tuba1b), previously reported to be appropriate for gene expression analysis in zebrafish . For detailed p-values see Additional file 3: Table S3. Amplification of pde6i gene started between Cq 28 and 34, when using ten times higher amount of adult eye cDNA template (100 ng) for the 5 p.m. time point, indicating low expression levels. Pde6i expression levels were also investigated in 3dpf, 6dpf and 10dpf embryos, showing amplification between Cq 29 and 34, either using 10 ng or 50 ng of cDNA per reaction. Additionally, when searching zebrafish EST sequences using tblastn with the mRNA sequence for pde6i (GenBank accession no. XM_003198119.3), no identical sequences were found. Therefore, pde6i was excluded from expression pattern analyses.
In rods (Fig. 7a), significant (p < 0.05) but minor oscillations in gene expression (~2 fold) were observed over 24 h for the 2R duplicates pde6a and pde6b, which code for the rod-specific catalytic subunits. Interestingly, the 3R paralogs pde6ga and pde6gb, which code for the inhibitory subunits, behave quite differently: pde6ga has the highest expression in the early morning (08:00, p < 0.05), after a constant ~20 fold increase since late evening (20:00). Its expression decreases rapidly in the morning (08:00–12:00), to levels comparable to pde6a (Fig. 7a). The expression levels of pde6gb also oscillate, but differently in time and levels: its highest expression is early at night (24:00), after a quick ~10 fold increase since late evening (20:00, p < 0.05). After 24:00, its expression diminishes until noon (12:00). During the day, there is a fairly constant expression for both paralogs; pde6gb being ~10 fold higher than pde6ga.
In cones (Fig. 7b) the pde6c gene, which codes for the catalytic subunit, shows a minor oscillation in its expression during the day (~5 fold), with a significant peak of expression in the early morning (08:00) (Fig. 7b, p < 0.05). However, this oscillation is minimal compared to the oscillation in the expression pde6ha, coding for an inhibitory subunit, which gradually increases its expression (~120 fold) from the middle of the night (04:00) to reaching its peak in the evening (17:00, p < 0.05), to decrease again after 17:00. An oscillation of ~20 fold in the expression of pde6hb was also observed, increasing expression during the night (24:00–08:00) and decreasing during the day (08:00–20:00) (Fig. 7b).
This report describes the evolutionary history of all PDE6 genes as deduced from analyses of a wide range of vertebrate genomes. A résumé of their evolution, including the number and type of genes present in each vertebrate group, is shown in Fig. 8. In addition, zebrafish was used as a model to investigate the possible sub- or neofunctionalisations of PDE6 inhibitory subunit gene 3R paralogs.
The PDE6 catalytic subunit gene repertoire is well conserved within the subphylum Vertebrata
The phylogenetic analysis of the PDE6 catalytic subunit gene family and our data on conserved synteny (including phylogenies for the neighbouring gene families) support an expansion in 2R for this gene family (Figs. 1 and 2). The three-gene repertoire (PDE6A, PDE6B and PDE6C) is well conserved across vertebrates, with the exception of northern hemisphere lampreys, that have only one gene and reptiles and birds that lack PDE6A. Additionally, the teleost-specific genome doubling did not result in any extra duplicates in this gene family, with the possible exception of the pufferfishes, which have an additional PDE6C.
Within the phototransduction cascade, our group has previously reported similar cases of gene repertoire conservation for the visual opsin gene family  and the alpha transducin gene family . Altogether, these results reflect the high degree of conservation of the catalytic function, with no or little selection pressure to retain 3R duplicates; there might even have been selection against their retention.
Functionally, the three PDE6 catalytic subunits have similar catalytic activity and interaction with the PDE6 inhibitory subunits . The difference is the efficiency of activation by transducin, in that rod PDE6 is activated much less efficiently than cone PDE6 . The lack of a PDE6A gene in birds and non-avian reptiles could be explained by a loss in their common ancestor before the divergence of these lineages. Such a gene loss would be in line with a previous study of the chicken retina  which suggests that chicken uses a PDE6ββ homodimer instead of the PDE6αβ heterodimer commonly used by other vertebrates [2, 25].
The appearance of an inhibitory subunit that regulates the PDE6 activity facilitated the success of the vertebrate visual system
The absence of PDE6 inhibitory subunit genes in invertebrates and their location in the SSTR2, 3, 5 and UTSR paralogon [16, 17], has led us to propose the following scenario for the evolution of this gene family. An ancestral PDE6 inhibitory subunit gene emerged de novo in the vertebrate lineage closely before 1R facilitating the shutting off of the PDE6 enzyme during dark conditions . The appearance of de novo genes from non-coding genomic regions has been suggested to be important mechanism contributing to the origin of genes, where they can evolve important and even essential functions rapidly . Subsequently, this gene was quadrupled in 2R followed by a loss of one of the four resulting genes.
Within gnathostomes, we here show for the first time that actinopterygians and non-amniote sarcopterygians have retained a third gene, that we named PDE6I to mark its paralogous relationship with PDE6G and PDE6H. We also found extra duplicates of PDE6G and PDE6H in teleosts, both resulting from 3R and through local duplications. In zebrafish we found that both the pde6g and pde6h genes have retained 3R duplicates. We have named the genes pde6ga, pde6gb, pde6ha and pde6hb. Additionally, zebrafish has the pde6i gene for which we detected very low expression in adult eyes and whole embryos up to 10 dpf, therefore we suggest that it might not have a function or, in adults, it is expressed in another organ.
The different intron losses that occurred in the PDE6H genes within teleosts can be the result of incomplete retrotranscription of the mRNA or retrotranscription of a partially spliced pre-mRNA, followed by replacement of the gene resulting in loss of one or more introns . Similar processes has been observed in many genes [28–31]. It has been observed that the hatching enzyme genes with the highest expression levels are more prone to intron loss . This is in line with our analyses that the expression of the pde6ha gene is much higher than the other genes (Fig. 7).
Among non-amniote sarcopterygians, we found PDE6I in the Western clawed frog and the coelacanth. In coelacanth PDE6I has probably been pseudogenised due to a frame-shift mutation that disrupts normal translation. However, due to the high overall sequence conservation of the coelacanth gene, the mutation is either recent or due to a sequencing error.
The sea lamprey genome assembly (see Table 2) has low coverage (5X) and is based on somatic cell DNA. This poses a problem because it is known that the lamprey genome goes through programmed loss of large amounts (~20 %) of somatic DNA in different cell lineages during embryonic development . Thus, genes might be present in the genome but missing in the genome assembly. When searching the germline genome assembly of the arctic lamprey we found fragmented PDE6 inhibitory subunit gene sequences. The lifestyle of the northern hemisphere lampreys, mostly in pelagic and benthic waters , probably has led to losses in genes related to vision. Thus, to fully understand the visual gene repertoire, a genome assembly from a southern hemisphere lamprey, such as Geotria australis, living predominantly in surface waters  and known to have retained all ancestral subtypes of visual opsins [22, 34, 35], is needed.
What did the visual system of zebrafish gain by retaining PDE6 inhibitory subunit duplicates after 3R?
The retention of duplicated genes may have two possible outcomes; subfunctionalisation, whereby the functions and/or expression of the ancestral gene are partitioned between the duplicates, or neofunctionalisation, whereby one or both of the copies gain novel functions [36, 37].
Here we demonstrate that the paralogs retained after 3R in zebrafish for the PDE6 inhibitory subunits, pde6ga, pde6gb, pde6ha and pde6hb, did not undergo subfunctionalisation involving topographical or temporal specialisations during development, i.e. the four genes are expressed in the retina and the pineal complex from early development. Instead, both paralogous pairs are coexpressed in the same retinal photoreceptor cell types as their amniote orthologs; pde6ga and pde6gb in rods and pde6ha and pde6hb in cones. An exception could be the SSC (UV-opsin containing cones) which possess some rod-like physiology features (see ) that could explain the possible expression of the four PDE6 inhibitory subunits, in addition to the previously suggested gnb1a and gnb1b  and rcv1a . Further analyses are planned to be conducted in this sense in order to solve the transcriptome of each zebrafish cone type.
We also observed a striking difference in expression levels between day and night. Altogether, the high sequence identity, coexpression and different expression levels of the paralogs led us to consider neofunctionalisation as unlikely to have happened. Instead we suggest that the overall regulatory function by the PDE6 inhibitory subunits of the phototransduction cascade is most likely maintained, while specialisation of one or several of mechanisms of this regulation may be affected , with a direct effect on gene expression levels.
Similar coexpression in rods of paralogs involved in the phototransduction cascade has also been found for the gnb1a and gnb1b transducin subunit genes  and has been suggested for the arrSa and arrSb arrestin genes  and the grk1a and grk1b opsin GPCR kinase genes . Additionally, subfunctionalisation of the zebrafish cone-specific arrestins has been reported previously: arr3a expressed in double cones and arr3b in single cones . The current study shows a case of retention of two 3R paralogs in all cones similarly to the recently reported recoverin genes  which also might have been duplicated in 3R.
Altogether, zebrafish has retained several 3R duplicates of phototransduction cascade components that are coexpressed in the same photoreceptor cell types, indicating a dosage effect in regulation of the visual function. In this article, we show that the coexpression of PDE6 inhibitory subunit genes, inhibiting the main effector of the phototransduction cascade, is subfunctionalised by showing a daily oscillation in gene expression (see next section).
The pineal complex of non-mammalian vertebrates has rod- and cone-like photoreceptor cells [43, 44]. The expression of the four PDE6 inhibitory subunit genes in the pineal complex agrees with previous observations of the presence of several other components of the phototransduction cascade [39, 40, 45–47]. Altogether, their expression further supports the idea that the pinealocytes share a common origin with the retinal photoreceptors [38, 39, 48]. However, the functional role for these duplicated genes remains unclear.
The onset of expression for the PDE6 inhibitory subunit genes is earlier in the pineal complex than in the retina, similar to other components of the teleost phototransduction cascade like opsins [49–51], transducins  or cyclic nucleotide gated channels (CNGs) (unpublished results) and recoverins . This is in accordance with the involvement of the pineal complex in regulation of hatching in Atlantic halibut . The delay of expression in the retina is consistent with previous reports that have shown the zebrafish retina to be responsive to light stimuli around 3 dpf .
Opposite rhythmic oscillations in expression for the PDE6 inhibitory subunit genes are related to differential light sensitivity
Circadian changes in the zebrafish and cichlid retinae have been previously reported regarding morphology [18, 19], gene expression [53, 54] and physiology [54, 55]. How those circadian changes influence visual function is not well understood. Here we describe gene expression data of regulatory genes that are likely to directly influence the sensitivity to light.
Zebrafish have been reported to show circadian oscillations in visual sensitivity , which can be partially explained by the striking morphological differences between day and night in the outer retina ; rod outer segments (rOS) are located in outer locations during the day and move inwards at night, while cone outer segments (cOS) are located in the innermost outer retina during the day and move outwards at night (see Fig. 9). These retinomotor movements, present in most vertebrate groups except mammals , facilitate an optimal light exposure for each photoreceptor type at the different light conditions, which increases their sensitivity, albeit losing resolution at night .
Putting together data about phototransduction cascade proteins, cone opsins  and pde6ha (present results) have their highest expression in the late afternoon, when the cones mediating photopic vision have high sensitivity . The cone transducin levels, on the contrary, do not oscillate significantly (unpublished results). High levels of opsins would increase sensitivity while pde6ha would be a key factor for the non-saturating effect, exclusive for cones , tentatively by the binding of the excess of free PDE6 inhibitory subunits to activated alpha transducins in a desensitisation-like effect. For scotopic vision, rhodopsin (unpublished results) and rod transducin  expression levels do not oscillate significantly, while expression of both PDE6G paralogs is higher in rods at night, where dim-light conditions activate rods.
Oscillations in the expression of opsins have been reported in teleosts [53, 54], frog  and chicken , mice , the Sudanian grass rat  with no obvious common denominator among all these results. These circadian changes in gene expression have been attributed to the circadian cycle of outer segment disc shedding of the photoreceptors  or correlated with behavioural sensitivity to light . This is the first report showing gene expression oscillations in a phototransduction cascade protein except opsins. Similarly, oscillations in gene expression has been observed for the interphotoreceptor retinoid binding protein, a gene specific for photoreceptors but not strictly phototransduction cascade related protein, in both zebrafish  and chicken .
In this study we investigated the mRNA expression pattern and expression levels due to the difficulty to generate specific antibodies able to target the 3R paralogs. Therefore, we have to assume that there is no or little delay in the translation of proteins or in the trafficking to and throughout the outer segment. Contrarily, it has been previously reported that transducin is stored in the photoreceptor inner segments in rats before being transported to the outer segments during the night . However, this has not been demonstrated in other species or for other proteins of the phototransduction cascade . In order to investigate the light-dependence of the PDE6 inhibitory subunits oscillations in gene expression levels and its effect on visual sensitivity, our group has initiated functional and semi-functional studies.
We report here that both catalytic and inhibitory subunit gene families of PDE6 were duplicated in 2R, resulting in cone- and rod-specific copies. Together with duplications of several other gene families expressed in the eye [3, 4, 22, 23, 41] this emphasises further the importance of 2R for vertebrate evolution in general and particularly for vision. The 2R event facilitated the differentiation of rods and cones by subfunctionalisation and/or neofunctionalisation of gene duplicates. Additionally, we report that the important regulatory PDE6 inhibitory subunit genes seem to have arisen in conjunction with the origin of the vertebrates, shortly before the genome quadruplication in 2R, presumably facilitating the rapid regulation of visual responses of the vertebrate eye mediated by the PDE6 enzyme. After 2R, the inhibitory subunit gene duplicates became specialised on rods and cones, respectively, and a third ancient copy presumably acquired distinct functions and still exists in some vertebrate lineages, although not amniotes. Finally, we show that the 3R-generated PDE6 inhibitory subunit gene duplicates evolved dramatically different daily expression profiles in zebrafish. Thus, the PDE6 inhibitory gene duplicates display novel specialisations that warrant further functional experiments.
Collection of amino acid sequences for phylogenetic analyses
Amino acid sequence predictions of the PDE6 subunit genes were downloaded from the Ensembl genome browser version 69 (October 2012) and the Pre! Ensembl genome browser from the species listed in Table 2. In addition human PDE6 subunit amino acid sequences were used as TBLASTN queries for searches against the genome assemblies of the species where genes were missing. Regions in the genome assemblies with BLAST hits, but without any Ensembl gene prediction, were downloaded and predicted using GENSCAN  (available at the web server: http://genes.mit.edu/GENSCAN.html) or predicted manually following splice donor-acceptor sites and sequence homology. The sequences for the human PDE5A and PDE11A genes were included in the analysis to serve as outgroup to root the trees.
In order to identify invertebrate orthologs of the PDE6 inhibitory subunit genes, human PDE6G and PDE6H amino acid sequences were used as queries in TBLASTN searches against the whole genome shotgun (WGS) sequences database from the National Centre for Biotechnology Information (NCBI) where all vertebrate sequences was excluded to try to identify invertebrate homologs. Hidden Markov Model (HMM) searches were performed against protein and nucleotide databases both on http://hmmer.janelia.org and locally . A HMM based on a nucleotide alignment of human, chicken and sea lamprey PDE6 inhibitory subunit gene sequences was generated using a local installation of hmmbuild and used as input to a local installation of nhmmer against several invertebrate genomes listed in Table 2.
Alignments and phylogenetic analyses
Amino acid sequences of each protein family were aligned using ClustalO with standard settings within the Seaview 4.5.3 program . Stretches of poorly aligned sequences, that seemed to be the result bad sequence predictions, were checked towards the genomic sequence. To identify the correct prediction of these parts genscan was used for the region of the gene. If genscan did not find the sequence, manual annotation was tried following sequence homology as well as consensus donor and acceptor splice sites in the genomic sequence in combination with TBLASTN searches.
Alignments used for phylogenetic analysis were tested for the best amino acid substitution model using ProtTest 3.2  with the following settings; Starting topology: Fixed BioNJ JTT, Matrices: JTT, LG, DCMut, Dayhoff, WAG, RtREV, CpREV, Blosum62 and VT.
PhyML trees were created using the PhyML 3.0 web server  (available at: http://www.atgc-montpellier.fr/phyml/) or locally using Seaview 4.5.3 using the following settings; JTT was used as substitution model for all of the neighbouring gene families while WAG was used for the PDE6 catalytic subunit gene family. Equilibrium frequencies were set to empirical, proportion of invariable sites estimated, number of substitution rate categories eight, gamma shape parameter estimated, starting tree BIONJ, type of tree improvement SPR & NNI, number of random starting trees no, optimise topology yes, optimise branch lengths yes, compute aLRT no and finally bootstrap was used with 100 replicates.
Analysis of conserved synteny
To analyse the conservation of synteny for the chromosomal regions harbouring the PDE6 catalytic subunit genes, information on associated gene names, Ensembl gene IDs, Ensembl transcript IDs, genomic location and Ensembl protein family IDs were downloaded from Ensembl version 61 in a region spanning approximately 10 Mb upstream and downstream of the PDE6A, PDE6B and PDE6C genes in the human genome. The chosen region for human chromosome 4 was 1 bp to 10.7 Mb (PDE6B), 139.2 Mb to 159.3 Mb for chromosome 5 (PDE6A) and 85.4 Mb to 105.4 Mb for chromosome 10 (PDE6C). Ensembl protein families with members in at least two of the selected chromosomal regions were selected for subsequent phylogenetic analysis and sequences were retrieved from Ensembl version 65.
The PDE6 inhibitory subunit genes were upon initial inspections in tetrapod genomes found to be located in the SSTR2, 3 and 5 paralogon previously described by . Thus, we were able to use their supplementary chromosomal tables to assign the orthology of most of the PDE6 inhibitory subunit genes. To investigate the conservation of the paralogon in the spotted gar genome, a species not included in their analysis, we investigated a selection of neighbouring gene families in regions spanning 5 mega base pairs (mb) upstream and downstream of the spotted gar PDE6 inhibitory subunit genes. We identified and analysed several of the neighbouring gene families they described, indicating that the PDE6 inhibitory subunit genes indeed are located in the same paralogon in the spotted gar (data not shown). Additionally, to be able to assign the orthology of the PDE6 inhibitory subunit genes located in difficult regions (too low resolution in the  analysis) and on short scaffolds, lists of spotted gar genes, also containing Ensembl’s predicted orthologs and their locations in other species, located in regions 1 Mb upstream and downstream of the spotted gar PDE6 inhibitory subunit genes were downloaded from Ensembl version 75. The included species in this orthology analysis were the Western clawed frog, coelacanth, zebrafish, three-spined stickleback, Nile tilapia, Atlantic cod and medaka.
Animals used for gene expression analyses
The animals used in this study were from three zebrafish lines: AB strain zebrafish (n = 45) purchased from the Science for Life Laboratory Zebrafish Technology Platform (Uppsala University, Sweden), Tg(gnat2:EGFP) zebrafish (n = 7)  Tg(rho:EGFP) (n = 8) . All adult zebrafish used in the experiments were from 6 to 18 months old and were maintained under standard conditions with lights on at 8 a.m. and off at 10 p.m.in a 14–10 light–dark cycle and the change of lighting conditions was abrupt.
Quantitative reverse-transcriptase PCR
For quantitative reverse-transcriptase PCR (qRT-PCR) experiments, the fish were anesthetised using Tricaine (0.04 % in system water) and decapitated. The heads were immersed in RNAlater® (Qiagen: cat. no. 76104) prior to dissection of the eyes. Both eyes of each individual were pooled together for RNA extraction. Total RNA was extracted using the RNeasy mini kit (Qiagen: cat. no. 74104). Three to four individuals were collected at six time points during the day: 08:00, 12:00, 17:00, 20:00, 00:00, 04:00. Following RNA extraction, the RNA was precipitated in ethanol for increased purity, treated with DNase I (Thermo Scientific: cat. no. EN0521) and used as template for reverse transcription using iScript™ cDNA synthesis kit (Bio-Rad: cat. no. 170–8890) in both RT and no-RT reactions.
Each qRT-PCR reaction was set up as follows: 10 μl iQ™ SYBR® Green Supermix (Bio-Rad: cat. no. 170–8880), 1,25 μl forward- and 1,25 μl reverse primers (10 pmol/μl), 6,5 μl MilliQ H2O and finally 1 μl 10 ng/μl cDNA or MilliQ H2O. Each reaction was run in triplicates, together with no-RT and no-template controls in duplicates. Primers were designed using primer-blast (http://blast.ncbi.nlm.nih.gov) spanning either exon-exon boundaries or placed in separate exons (Table 3). Actin beta 1 (actb1) and alpha tubulin 1b (tuba1b) were used as housekeeping genes according to  with our own primer pairs. qRT-PCR efficiency and amplification data was analysed using LinRegPCR version 2014.6  and presented as relative to the mean of actb1 and tuba1b (the 2-ΔCq method). Statistics were analysed for each gene using one-way ANOVA with a Tukey’s Multiple Comparison post-hoc test in GraphPad Prism 5.00. The mean relative expression for each gene and time-point were plotted in double plots (Fig. 7).
For the ontogenetic analysis of pde6i expression levels, 3 pools of 20 embryos were isolated for 3dpf, 6dpf and 10dpf embryos. Subsequently, their mRNA was extracted, cDNA was synthesised, the RT-qPCR was performed and the results were analysed as is described above.
Probe design and synthesis
In order to study the expression pattern of the four zebrafish PDE6 inhibitory subunit genes (pde6ga, pde6gb, pde6ha and pde6hb), antisense riboprobes were designed targeting their 3′ untranslated regions (3′UTR). PCR primers were designed using the Primer-BLAST tool , available at the NCBI webpage (http://blast.ncbi.nlm.nih.gov/Blast.cgi) (see Table 3) and PCR reactions were performed using genomic DNA from zebrafish. The resulting amplicons were cloned into pCR®II-TOPO® vectors (Invitrogen: cat. no. K4650-01), sequenced and used for probe synthesis with either T7 or SP6 RNA polymerase using DIG RNA labelling kit (Roche: cat. no. 11175025910) according to manufacturer’s instructions. Sense and antisense probes were synthesised and the former was used to control for specificity.
In situ hybridisation (ISH)
The fish used in the ISH experiments were anesthetized using Tricaine (0.04 % in system water), either during the day (17:00) or during the night (04:00). The heads were dissected and fixed by immersion in 4 % paraformaldehyde (PFA) diluted in phosphate buffer 0.1 M pH 7.4 (PB) for 7 h and washed in phosphate buffered saline 0.1 M pH 7.4 (PBS) overnight, both at 4 °C. Subsequently the heads were cryoprotected in 30 % sucrose and sectioned in a cryostat (Microm Cryo-Star HM 560) obtaining transversal sections 12–20 μm thick stuck on positively charged slides.
The ISH was performed according to  with minor adaptations. The final staining reaction was carried out using different substrates for the AP enzyme bounded to the Fab fragments: NBT/BCIP or Fast Red tablets (Roche: cat. no. 11681451001 and 11496549001, respectively). Sense probes were incubated in parallel with the antisense, as a specificity control, with no staining as a result. In order to preserve the fluorescence of the EGFP when using transgenic animals, the hybridisation reaction time was reduced from overnight to 6 h. Double ISH experiments according to  were also performed in specific cases. Lastly, all slides were mounted using VectaShield® mounting medium with DAPI incorporated.
Cell assignment was based on overall cell morphology, topological location of the nuclei and the mitochondria dense ellipsoids  of the different photoreceptor cell types, as well as immunohistochemistry on rod- or cone-specific EGFP fluorescence transgenic lines. Rod-specific rabbit anti-GNB1 (Nordic BioSite: cat. no. LS-C90703; 1:500) and double cones-specific mouse anti-zpr1 (ZIRC; 1:400) primary antibodies were used with donkey anti-rabbit coupled Alexa 555 (Life technologies cat. no. A31572; 1:1000) and donkey anti-mouse coupled Alexa 488 (Life Technologies cat. no. A21202; 1:1000) as secondary antibodies. We have used the nomenclature proposed by  for the different cone types: double cones (DC, middle and long wavelength), long single cones (LSC, short wavelength) and short single cones (SSC, ultraviolet).
Bright-field, fluorescence and Nomarski contrast photomicrographs, as well as their combinations were taken using a Zeiss Axioplan 2 microscope equipped with a Zeiss AxioCam camera or an inverted LSM510 Zeiss confocal microscope. The figures merged using CorelDRAW Graphics Suite X6.
Whole mount in situ hybridisation (WISH)
All embryos and larvae (1–6 dpf) were collected as described previously . The WISH experiments were performed in an InsituPro VSi (Intavis AG, Köln, Germany) ISH robot at the Science for Life Laboratory Zebrafish Technology Platform (Uppsala Universitet, Sweden). The staining was done outside of the robot using NBP/BCIP at 37 °C. The images were acquired using a stereomicroscope Nikon SMZ1500 with a Nikon DS-Vi1 camera and the figures merged using CorelDRAW Graphics Suite X6.
Conti M, Beavo J. Biochemistry and physiology of cyclic nucleotide phosphodiesterases: essential components in cyclic nucleotide signaling. Annu Rev Biochem. 2007;76(c):481–511.
Cote RH. Characteristics of photoreceptor PDE (PDE6): similarities and differences to PDE5. Int J Impot Res. 2004;16:S28–33.
Nordström K, Larsson TA, Larhammar D. Extensive duplications of phototransduction genes in early vertebrate evolution correlate with block (chromosome) duplications. Genomics. 2004;83:852–72.
Larhammar D, Nordström K, Larsson TA. Evolution of vertebrate rod and cone phototransduction genes. Philos Trans R Soc Lond B Biol Sci. 2009;364:2867–80.
Muradov H, Boyd KK, Kerov V, Artemyev NO. PDE6 in lamprey Petromyzon marinus: implications for the evolution of the visual effector in vertebrates. Biochemistry. 2007;46:9992–10000.
Nakatani Y, Takeda H, Kohara Y, Morishita S. Reconstruction of the vertebrate ancestral genome reveals dynamic genome reorganization in early vertebrates. Genome Res. 2007;17:1254–65.
Putnam NH, Butts T, Ferrier DEK, Furlong RF, Hellsten U, Kawashima T, Robinson-Rechavi M, Shoguchi E, Terry A, Yu J-K, Benito-Gutiérrez EL, Dubchak I, Garcia-Fernàndez J, Gibson-Brown JJ, Grigoriev I V, Horton AC, de Jong PJ, Jurka J, Kapitonov V V, Kohara Y, Kuroki Y, Lindquist E, Lucas S, Osoegawa K, Pennacchio LA, Salamov AA, Satou Y, Sauka-Spengler T, Schmutz J, Shin-I T, et al. The amphioxus genome and the evolution of the chordate karyotype. Nature. 2008;453:1064–71.
Yuasa K, Kanoh Y, Okumura K, Omori K. Genomic organization of the human phosphodiesterase PDE11A gene. Evolutionary relatedness with other PDEs containing GAF domains. Eur J Biochem. 2001;268:168–78.
Sundström G, Larsson TA, Larhammar D. Phylogenetic and chromosomal analyses of multiple gene families syntenic with vertebrate Hox clusters. BMC Evol Biol. 2008;8:254.
Guo L-W, Muradov H, Hajipour AR, Sievert MK, Artemyev NO, Ruoho AE. The inhibitory gamma subunit of the rod cGMP phosphodiesterase binds the catalytic subunits in an extended linear structure. J Biol Chem. 2006;281:15412–22.
Arshavsky VY, Lamb TD, Pugh EN. G proteins and phototransduction. Annu Rev Physiol. 2002;64:153–87.
Zhang XJ, Gao XZ, Yao W, Cote RH. Functional mapping of interacting regions of the photoreceptor phosphodiesterase (PDE6) γ-subunit with PDE6 catalytic dimer, transducin, and regulator of G-protein signaling9-1 (RGS9-1). J Biol Chem. 2012;287:26312–20.
Nishiwaki Y, Komori A, Sagara H, Suzuki E, Manabe T, Hosoya T, Nojima Y, Wada H, Tanaka H, Okamoto H, Masai I. Mutation of cGMP phosphodiesterase 6alpha’-subunit gene causes progressive degeneration of cone photoreceptors in zebrafish. Mech Dev. 2008;125:932–46.
Vihtelic TS, Fadool JM, Gao J, Thornton KA, Hyde DR, Wistow G. Expressed sequence tag analysis of zebrafish eye tissues for NEIBank. Mol Vis. 2005;11:1083–100.
Larsson TA, Olsson F, Sundstrom G, Lundin L-G, Brenner S, Venkatesh B, Larhammar D. Early vertebrate chromosome duplications and the evolution of the neuropeptide Y receptor gene regions. BMC Evol Biol. 2008;8:184.
Ocampo Daza D, Sundström G, Bergqvist CA, Larhammar D. The evolution of vertebrate somatostatin receptors and their gene regions involves extensive chromosomal rearrangements. BMC Evol Biol. 2012;12:231.
Tostivint H, Ocampo Daza D, Bergqvist CA, Quan FB, Bougerol M, Lihrmann I, Larhammar D. Molecular evolution of GPCRs: Somatostatin/urotensin II receptors. J Mol Endocrinol. 2014;52:T61–86.
Menger GJ, Koke JR, Cahill GM. Diurnal and circadian retinomotor movements in zebrafish. Vis Neurosci. 2005;22:203–9.
Hodel C, Neuhauss SCF, Biehlmaier O. Time course and development of light adaptation processes in the outer zebrafish retina. Anat Rec A: Discov Mol Cell Evol Biol. 2006;288:653–62.
Schmitt EA, Dowling JE. Early retinal development in the zebrafish, Danio rerio: light and electron microscopic analyses. J Comp Neurol. 1999;404:515–36.
McCurley AT, Callard GV. Characterization of housekeeping genes in zebrafish: male–female differences and effects of tissue type, developmental stage and chemical treatment. BMC Mol Biol. 2008;9:102.
Lagman D, Daza DO, Widmark J, Abalo XM, Sundström G, Larhammar D. The vertebrate ancestral repertoire of visual opsins, transducin alpha subunits and oxytocin/vasopressin receptors was established by duplication of their shared genomic region in the two rounds of early vertebrate genome duplications. BMC Evol Biol. 2013;13:238.
Lagman D, Sundström G, Ocampo Daza D, Abalo XM, Larhammar D. Expansion of Transducin Subunit Gene Families in Early Vertebrate Tetraploidizations. Genomics. 2012;100:203–11.
Muradov H, Boyd KK, Artemyev NO. Rod phosphodiesterase-6 PDE6A and PDE6B subunits are enzymatically equivalent. J Biol Chem. 2010;285:39828–34.
Huang D, Hinds TR, Martinez SE, Doneanu C, Beavo JA. Molecular determinants of cGMP binding to chicken cone photoreceptor phosphodiesterase. J Biol Chem. 2004;279:48143–51.
Wu DD, Zhang YP. Evolution and function of de novo originated genes. Mol Phylogenet Evol. 2013;67:541–5.
Jeffares DC, Mourier T, Penny D. The biology of intron gain and loss. Trends Genet. 2006;22:16–22.
da Lage JL, Maczkowiak F, Cariou ML. Phylogenetic distribution of intron positions in alpha-amylase genes of bilateria suggests numerous gains and losses. PLoS One. 2011;6:e19673.
Kawaguchi M, Hiroi J, Miya M, Nishida M, Iuchi I, Yasumasu S. Intron-loss evolution of hatching enzyme genes in Teleostei. BMC Evol Biol. 2010;10:260.
Zhou M, Yan J, Ma Z, Zhou Y, Abbood NN, Liu J, Su L, Jia H, Guo AY. Comparative and evolutionary analysis of the HES/HEY gene family reveal exon/intron loss and teleost specific duplication events. PLoS One. 2012;7:e40649.
Bellingham J, Tarttelin E, Foster R, Wells DJ. Structure and Evolution of the Teleost Extraretinal Rod-Like Opsin ( errlo ) and Ocular Rod Opsin ( rho ) Genes: Is Teleost rho a Retrogene? J Exp Zool Part B Mol Dev Evol. 2003;297B:1–10.
Smith JJ, Antonacci F, Eichler EE, Amemiya CT. Programmed loss of millions of base pairs from a vertebrate genome. Proc Natl Acad Sci U S A. 2009;106:11212–7.
Renaud CB. Lampreys of the World: An Annotated and Illustrated Catalogue of Lamprey Species Known to Date. Volume 1. 5th ed. 2011.
Collin SP, Knight MA, Davies WL, Potter IC, Hunt DM, Trezise AEO. Ancient colour vision: multiple opsin genes in the ancestral vertebrates. Curr Biol. 2003;13:R864–5.
Collin SP, Davies WL, Hart NS, Hunt DM. The evolution of early vertebrate photoreceptors. Philos Trans R Soc Lond B Biol Sci. 2009;364:2925–40.
Lynch M, Conery JS. The evolutionary fate and consequences of duplicate genes. Science. 2000;290:1151–5.
Lynch M, Force A. The probability of duplicate gene preservation by subfunctionalization. Genetics. 2000;154:459–73.
Lamb TD. Evolution of phototransduction, vertebrate photoreceptors and retina. Prog Retin Eye Res. 2013;36:52–119.
Lagman D, Callado-Pérez A, Franzén IE, Larhammar D, Abalo XM. Transducin Duplicates in the Zebrafish Retina and Pineal Complex: Differential Specialisation after the Teleost Tetraploidisation. PLoS One. 2015;10:e0121330.
Zang J, Keim J, Kastenhuber E, Gesemann M, Neuhauss SCF. Recoverin depletion accelerates cone photoresponse recovery. Open Biol. 2015;5:150086.
Renninger SL, Gesemann M, Neuhauss SCF. Cone arrestin confers cone vision of high temporal resolution in zebrafish larvae. Eur J Neurosci. 2011;33:658–67.
Wada Y, Sugiyama J, Okano T, Fukada Y. GRK1 and GRK7: Unique cellular distribution and widely different activities of opsin phosphorylation in the zebrafish rods and cones. J Neurochem. 2006;98:824–37.
Tamotsu S, Oishi T, Nakao K, Fukada Y, Shichida Y, Yoshizawa T, Morita Y. Localization of iodopsin and rod-opsin immunoreactivity in the retina and pineal complex of the river lamprey, Lampetra japonica. Cell Tissue Res. 1994;278:1–10.
Allwardt BA, Dowling JE. The pineal gland in wild-type and two zebrafish mutants with retinal defects. J Neurocytol. 2001;30:493–501.
Mano H, Kojima D, Fukada Y. Exo-rhodopsin: a novel rhodopsin expressed in the zebrafish pineal gland. Brain Res Mol Brain Res. 1999;73:110–8.
Rinner O, Makhankov YV, Biehlmaier O, Neuhauss SCF. Knockdown of cone-specific kinase GRK7 in larval zebrafish leads to impaired cone response recovery and delayed dark adaptation. Neuron. 2005;47:231–42.
Li X, Montgomery J, Cheng W, Noh JH, Hyde DR, Li L. Pineal photoreceptor cells are required for maintaining the circadian rhythms of behavioral visual sensitivity in zebrafish. PLoS One. 2012;7:e40508.
Falcón J, Coon SL, Besseau L, Cazaméa-Catalan D, Fuentès M, Magnanou E, Paulin C-H, Boeuf G, Sauzet S, Jørgensen EH, Mazan S, Wolf YI, Koonin EV, Steinbach PJ, Hyodo S, Klein DC. Drastic neofunctionalization associated with evolution of the timezyme AANAT 500 Mya. Proc Natl Acad Sci U S A. 2014;111:314–9.
Van Veen T, Ekström P, Nyberg L. Serotonin and opsin immunoreactivities in the developing pineal organ of the three-spined stickleback, Gasterosteus aculeatus L. Cell Tissue Res. 1984;237:559–64.
Negishi K, Wagner H. Differentiation of photoreceptors, glia, and neurons in the retina of the cichlid fish Aequidens pulcher; an immunocytochemical study. Dev Brain Res. 1995;89:87–102.
Forsell J, Holmqvist B, Helvik JV, Ekstrom P. Role of the pineal organ in the photoregulated hatching of the Atlantic halibut. Int J Dev Biol. 1997;41:591–5.
Easter Jr S, Nicola G. The Development of Vision in the Zebrafish ( Danio rerio). Dev Biol. 1996;663:646–63.
Halstenberg S, Lindgren KM, Samagh SPS, Nadal-Vicens M, Balt S, Fernald RD. Diurnal rhythm of cone opsin expression in the teleost fish Haplochromis burtoni. Vis Neurosci. 2005;22:135–41.
Li P, Temple S, Gao Y, Haimberger TJ, Hawryshyn CW, Li L. Circadian rhythms of behavioral cone sensitivity and long wavelength opsin mRNA expression: a correlation study in zebrafish. J Exp Biol. 2005;208(Pt 3):497–504.
Li L, Dowling JE. Zebrafish visual sensitivity is regulated by a circadian clock. Vis Neurosci. 1998;15:851–7.
Burnside B. Light and circadian regulation of retinomotor movement. Prog Brain Res. 2001;131:477–85. Review.
Haeri M, Calvert PD, Solessio E, Pugh EN, Knox BE. Regulation of rhodopsin-eGFP distribution in transgenic xenopus rod outer segments by light. PLoS One. 2013;8:e80059.
Pierce ME, Sheshberadaran H, Zhang Z, Fox LE, Applebury M, Takahashi JS. Circadian Regulation of lodopsin Gene Expression in Embryonic Photoreceptors in ‘ Retinal Cell Cult & e. Neuron. 1993;10:579–84.
Von Schantz M, Lucas R, Foster R. Circadian oscillation of photopigment transcript levels in the mouse retina. Mol Brain Res. 1999;72:108–14.
Bobu C, Sandu C, Laurent V, Felder-Schmittbuhl MP, Hicks D. Prolonged light exposure induces widespread phase shifting in the circadian clock and visual pigment gene expression of the Arvicanthis ansorgei retina. Mol Vis. 2013;19:1060–73.
Rajendran RR, Van Niel EE, Stenkamp DL, Cunningham LL, Raymond PA, Gonzalez-Fernandez F. Zebrafish interphotoreceptor retinoid-binding protein: differential circadian expression among cone subtypes. J Exp Biol. 1996;199(Pt 12):2775–87.
Stenkamp DL, Calderwood JL, Van Niel EE, Daniels LM, Gonzalez-Fernandez F. The interphotoreceptor retinoid-binding protein (IRBP) of the chicken (Gallus gallus domesticus). Mol Vis. 2005;11:833–45.
Brann MR, Cohen LV. Diumal Expression of Transducin mRNA and Translocation of Transducin in Rods of Rat Retina. Sciencen (80- ). 1987;235:585–7.
Burge C, Karlin S. Prediction of Complete Gene Structures in Human Genomic DNA. J Mol Biol. 1997;268:78–94.
Finn RD, Clements J, Eddy SR. HMMER web server: interactive sequence similarity searching. Nucleic Acids Res. 2011;39(Web Server issue):W29–37.
Gouy M, Guindon S, Gascuel O. SeaView version 4: A multiplatform graphical user interface for sequence alignment and phylogenetic tree building. Mol Biol Evol. 2010;27:221–4.
Abascal F, Zardoya R, Posada D. ProtTest: selection of best-fit models of protein evolution. Bioinformatics. 2005;21:2104–5.
Guindon S, Dufayard JF, Lefort V, Anisimova M, Hordijk W, Gascuel O. New Algorithms and Methods to Estimate Maximum-Likelihood Phylogenies: Assessing the Performance of PhyML 3.0. Syst Biol. 2010;59:307–21.
Smyth VA, Di Lorenzo D, Kennedy BN. A novel, evolutionarily conserved enhancer of cone photoreceptor-specific expression. J Biol Chem. 2008;283:10881–91.
Williams A, Sarkar S, Cuddon P, Ttofi EK, Saiki S, Siddiqi FH, Jahreiss L, Fleming A, Pask D, Goldsmith P, O’Kane CJ, Floto RA, Rubinsztein DC. Novel targets for Huntington’s disease in an mTOR-independent autophagy pathway. Nat Chem Biol. 2008;4:295–305.
Ruijter JM, Ramakers C, Hoogaars WMH, Karlen Y, Bakker O, van den Hoff MJB, Moorman AFM. Amplification efficiency: linking baseline and bias in the analysis of quantitative PCR data. Nucleic Acids Res. 2009;37:e45.
Ye J, Coulouris G, Zaretskaya I, Cutcutache I, Rozen S, Madden TL. Primer-BLAST: a tool to design target-specific primers for polymerase chain reaction. BMC Bioinformatics. 2012;13:134.
Hauptmann G, Gerster T. Multicolor whole-mount in situ hybridization. Methods Mol Biol. 2000;137:139–48.
Hauptmann G. One-, two-, and three-color whole-mount in situ hybridization to Drosophila embryos. Methods. 2001;23:359–72.
Tarboush R, Novales Flamarique I, Chapman GB, Connaughton VP. Variability in mitochondria of zebrafish photoreceptor ellipsoids. Vis Neurosci. 2014;31:11–23.
Raymond PA, Barthel LK. A moving wave patterns the cone photoreceptor mosaic array in the zebrafish retina. Int J Dev Biol. 2004;48:935–45.
Lagman D, Franzén IE, Eggert J, Larhammar D, Abalo XM. Data from: Evolution and expression of the phosphodiesterase 6 genes unveils vertebrate novelty to control photosensitivity. Dryad Digital Repository. 2016. doi:10.5061/dryad.k2f04.
We thank Katarzyna Rogóz for valuable help with probe synthesis, Kateryna Shebanits and Aniruddha Todkar for valuable advice on qRT-PCR assays and analysis, Christina Bergqvist for help with figures. This work was funded by the Swedish Research Council and the Carl Trygger’s and Olle Engkvist Byggmästare Foundations.
Availability of data and materials
The datasets supporting the conclusions of this article are available in the Dryad Digital repository, at doi:10.5061/dryad.k2f04 .
DLag, DLar and XMA conceived and designed the experiments and drafted the manuscript. DLag carried out all the phylogenetic and synteny analyses. DLag, IEF, JE and XMA performed the gene expression analyses. All authors read and approved the final manuscript.
The authors declare that do not have any financial political, personal, religious, ideological, academic, intellectual, commercial or any other competing interests in relation to this manuscript.
Ethics approval and consent to participate
All experiments using animals were performed in accordance to the recommendations of the Federation of Laboratory Animal Science Associations and approved by the Uppsala ethical committee on animal experiments (Uppsala djurförsöksetiska nämnd), permit numbers C294/12 and C315/12.
Multiple sequence alignment in FASTA format for the amino acid sequences of the PDE6 catalytic subunit genes. This alignment was done using ClustalO in the Seaview 4.5.3 program and was one used for the phylogenetic analysis. Sequence names are provided in Additional file 2. (DOC 53 kb)
Tables containing sequence identifiers, genomic locations and sequence names. The tables include information for genes used in the phylogenetic analyses, both from the PDE6 subunit gene families and the identified neighbouring gene families. Each family has its own sheet in the excel file. (XLSX 58 kb)
This file lists all the figures (S1-S18) and tables (S1-S3). Figure S1: Phylogenetic maximum likelihood tree with additional PDE6 catalytic subunit gene amino acid sequences from bird and non-avian reptile genomes. No PDE6A genes could be identified in any of the investigated species, further supporting the loss of PDE6A in this lineage. Figure S2-S9: Phylogenetic maximum likelihood trees of the investigated neighbouring gene families of the PDE6 catalytic subunit genes. Figure S10-S17: Synteny comparisons between spotted gar linkage groups carrying PDE6 inhibitory subunit genes and chromosomes and scaffolds in other species where comparisons with tables from Ocampo Daza et al., 2012 were uninformative. Figure S18: Whole mount in situ hybridisation with PDE6 inhibitory subunit gene expression during zebrafish development (1dpf – 3dpf). Table S1: List of the identified neighbouring gene families of the PDE6 catalytic subunit gene paralogon. Table S2: Table listing the pair-wise percentage nucleotide sequence identity between the PDE6 inhibitory subunit gene 3ʹUTR sequences. Table S3: Table displaying the results from one-way ANOVA of the relative expression for each PDE6 subunit gene between the six time-points. (PDF 1990 kb)
Multiple sequence alignment in FASTA format for the amino acid sequences of the identified PDE6 inhibitory subunit genes. The alignment was done using ClustalO in the Seaview 4.5.3 program. Sequence names are provided in Additional file 2. (RTF 9 kb)
Synteny tables from Ocampo Daza et al. 2012 with our identified PDE6 inhibitory subunit genes added. Tables showing genomic locations for the PDE6 neighbouring gene families organized in chromosomal blocks from Ocampo Daza et al. 2012, to which we have added the genomic locations of the PDE6 inhibitory subunit genes we identified in the current study. (XLSX 87 kb)
About this article
Cite this article
Lagman, D., Franzén, I.E., Eggert, J. et al. Evolution and expression of the phosphodiesterase 6 genes unveils vertebrate novelty to control photosensitivity. BMC Evol Biol 16, 124 (2016). https://doi.org/10.1186/s12862-016-0695-z
- Genome duplication
- Retinomotor movements