Circadian regulation of vertebrate cone photoreceptor function
Abstract
Eukaryotes generally display a circadian rhythm as an adaption to the reoccurring day/night cycle. This is particularly true for visual physiology that is directly affected by changing light conditions. Here we investigate the influence of the circadian rhythm on the expression and function of visual transduction cascade regulators in diurnal zebrafish and nocturnal mice. We focused on regulators of shut-off kinetics such as Recoverins, Arrestins, Opsin kinases, and Regulator of G-protein signaling that have direct effects on temporal vision. Transcript as well as protein levels of most analyzed genes show a robust circadian rhythm-dependent regulation, which correlates with changes in photoresponse kinetics. Electroretinography demonstrates that photoresponse recovery in zebrafish is delayed in the evening and accelerated in the morning. Functional rhythmicity persists in continuous darkness, and it is reversed by an inverted light cycle and disrupted by constant light. This is in line with our finding that orthologous gene transcripts from diurnal zebrafish and nocturnal mice are often expressed in an anti-phasic daily rhythm.
Introduction
Circadian rhythms serve as endogenous clocks that molecularly support the daily occurring oscillations of physiology and ensuing behavior (Brown et al., 2019; Cahill, 2002; Frøland Steindal and Whitmore, 2019; Golombek et al., 2014; Idda et al., 2012; Ukai and Ueda, 2010; Vatine et al., 2011). It has long been recognized that the central pacemaker of circadian rhythms resides in dedicated brain regions, either the suprachiasmatic nucleus in mammals or the pineal gland in non-mammalian vertebrates. The rhythm is entrained by external stimuli (eg, light) that directly act on the core circadian transcriptional feedback loop. Multiple studies have shown that autonomous circadian clocks also exist in other brain regions and in peripheral tissues (Frøland Steindal and Whitmore, 2019; Idda et al., 2012; Vatine et al., 2011). This is particularly true for the retina, which generates its own circadian rhythm (Gladys, 2020). In zebrafish, this rhythmicity is reflected in a number of circadian adaptations, such as a higher response threshold in the morning (Li and Dowling, 1998), photoreceptor retinomotor movement in constant darkness (Menger et al., 2005), and cone photoreceptor synaptic ribbon disassembly at night (Emran et al., 2010). Such adaptations are also found in other animals such as mice, where stronger electrical retinal coupling during the night (Jin et al., 2015; Li et al., 2009; Ribelayga et al., 2008), as well as slower dark adaptation of rods during the day, was observed (Xue et al., 2015). The molecular mechanisms underlying these circadian-dependent retinal regulations are still largely unknown.
In the vertebrate retina, there are two different types of photoreceptors, namely rods and cones (Burns and Baylor, 2001; Fu and Yau, 2007). Rods function mainly during dim light conditions, whereas cones are characterized by lower sensitivity but faster response kinetics, being important for daylight and color vision. About 92% of larval and 60% of adult photoreceptors in the zebrafish retina are cones (Allison et al., 2010; Fadool, 2003; Zimmermann et al., 2018). Although rods and cones generally use the same visual transduction cascade components, the individual reactions are typically mediated by photoreceptor type-specific proteins.
Visual transduction commences by an opsin chromophore-mediated absorption of photons, which triggers the activation of a second messenger cascade including the trimeric G-protein transducin. Activated transducin stimulates the effector enzyme phosphodiesterase (PDE), which leads to a reduction in intracellular cyclic guanosine monophosphate (cGMP) levels, subsequently leading to the closure of cyclic nucleotide -gated (CNG) cation channels, resulting in a membrane potential change (Fain et al., 2001; Lamb and Pugh, 2006).
High-temporal resolution requires a tightly regulated termination of visual transduction (Chen et al., 2012; Matthews and Sampath, 2010; Zang and Matthews, 2012). This depends on the highly effective quenching of both the activated visual pigment (R*) and the PDE-transducin complex (PDE*). R* is phosphorylated by a G-protein receptor kinase (GRK) before being completely deactivated by binding to arrestin. While GRK activity itself is controlled by recoverin (RCV) in a Ca2+-dependent manner (Zang and Neuhauss, 2018), the quenching of PDE* depends on the GTPase activity of its γ-subunit that is regulated by activator protein RGS9 (Regulator of G-protein Signaling 9) (Krispel et al., 2006).
We now show that the expression levels of these important regulators of cone visual transduction decay are modulated by the circadian clock. Moreover, these periodic fluctuations are reflected in oscillating protein levels that correlate with the rhythmicity in visual physiology and behavior observed in zebrafish. Interestingly, we have found that the expression of a selection of mouse orthologs of the investigated regulatory genes is also modulated by the circadian clock. However, the periodicity was opposite to that of zebrafish, fitting the nocturnal lifestyle of mice.
Results
Expression levels of key genes involved in shaping visual transduction decay are regulated by the circadian clock
To determine the influence of the circadian clock on visual behavior, we analyzed gene expression levels of key visual transduction regulators over a 24 hr period using quantitative real-time polymerase chain reacion (qRT-PCR). Eyes from larval (5 days post fertilization [dpf]) and adult zebrafish that were kept under a normal light cycle (LD 14:10, light on at 8 o’clock in the morning), as well as eyes from 5 dpf larvae kept in continuous darkness (DD), were collected every 3 hr over a period of 24 hr and subsequently analyzed. Apart from rcv2a, which seems to have no or weak fluctuating transcript levels in larvae (Figure 1G), expression levels of the other recoverins (rcv1a, rcv1b, which is absent from larval retina, and rcv2b), G-protein receptor kinases (grk7a and grk7b), arrestins (arr3a and arr3b), and regulator of G-protein signaling 9 (rgs9a) were clearly oscillating (statistical information in Supplementary file 1). In many cases, transcripts were most abundant at ZT1 or ZT4 (grk7a, grk7b, rcv2b, arr3a, and arr3b), subsequently declined throughout the day, and recovered during the night. For instance, in adult zebrafish eyes, grk7a expression levels decreased by around 98% from the peak to the lowest expression level (Figure 1A). In the case of adult rgs9a, transcripts reached the highest level at ZT22, with the value very close to ZT1. In situ hybridization (ISH) analysis using digoxigenin-labeled RNA probes validated our qRT-PCR results (Figure 1—figure supplement 2 and Figure 1—figure supplement 3).
Interestingly, two genes, namely rcv1a and rcv2a, displayed different expression profiles in larval and adult eyes (Figure 1E&G). While larval rcv1a mRNA transcript levels peaked around ZT19, larval rcv2a transcript expression was weak/non-cyclic. However, this is in contrast to adult retinas where rcv1a and rcv2a transcripts were highest at ZT7 (Figure 1G). An anti-phasic expression profile between larval and adult stages can also be observed for rod arrestins (arras) (Figure 1—figure supplement 4).
In order to establish that the daily expression changes of these transcripts are indeed regulated by the intrinsic circadian clock, we repeated our experiments in larvae kept in complete darkness (DD), eliminating light as an external factor. Under normal LD, as well as DD conditions, we obtained largely comparable results (Figure 1). Exceptions were arr3a and arr3b, showing a 3-hr phase shift, and rcv1a, showing an almost anti-phase relationship (see ‘Discussion’ section).
Corresponding retinal genes in nocturnal mice display an anti-phasic expression pattern
As zebrafish are diurnal animals having a cone-dominant retina, we wondered if the observed circadian regulation of visual transduction gene transcripts is also seen in the rod-dominant retina of nocturnal mice. We selected mouse Grk1, the only visual grk gene in mice (Chen et al., 1999; Wada et al., 2006), the sole recoverin (Chen et al., 2012) and Rgs9 (Krispel et al., 2006) genes, and the two arrestins Arrb1 and Arrb3, as the counterparts for the above-mentioned zebrafish genes for our analysis.
Expression of all five regulators fluctuated in a 24 hr period (Figure 2), being highest at the beginning of the dark period (ZT13) for the two arrestins (Figure 2A&B), or around midnight (ZT17) for Grk1, Rgs1, and Recvrn (Figure 2C–E). All of them displayed minimal transcript levels early during the day. This oscillation pattern shows a clear anti-phasic relationship with the cyclic fluctuation of the corresponding zebrafish transcripts. Curiously, the amplitude of gene fluctuation in adult zebrafish retina was generally larger than that in the mouse retina (Figures 1 and 2).
Levels of key visual transduction regulator proteins fluctuate in the zebrafish retina
While mRNA half-life is typically in the range of minutes, protein turnover rates can range from minutes to days, explaining why fluctuation of mRNA levels is not always reflected in time-shifted oscillations at the protein level (Cunningham and Gonzalez-Fernandez, 2000; Stenkamp et al., 2005). However, as regulatory proteins often have turnover rates of only a few hours, we were examining whether RNA oscillations are mirrored by corresponding protein level fluctuations. In order to assess protein levels, we generated paralog-specific antibodies against GRK7a and ARR3a. Quantitative western blot analysis indicated periodic changes in protein levels for both proteins. Peak expression was shifted 6 - 12 hr between RNA and protein level (Figure 3A&B). ARR3a reached its highest and lowest levels at ZT7 and ZT22, respectively, whereas GRK7a maintained relatively high levels throughout the day, having the lowest concentrations around midnight. Hence, mRNA circadian oscillations in the zebrafish retina are largely conserved at the protein level with a time shift.
Larval cone response recovery is delayed in the evening
We next asked whether the observed protein and RNA level fluctuations have an impact on functional aspects of visual transduction. Photoresponses at larval zebrafish stages are dominated by cone photoreceptors (Bilotta et al., 2001). In the electroretinogram (ERG), the a-wave directly represents photoreceptor responses. Since in the zebrafish ERG, it is largely masked by the larger b-wave, reflecting the depolarization of ON-bipolar cells, we used the b-wave amplitude as an indirect measure of the cone photoresponse (Figure 4A1). The protein products of the genes analyzed in our study are known to affect photoresponse recovery in zebrafish (Renninger et al., 2011; Rinner et al., 2005; Zang et al., 2015). Therefore, we assessed their function by using the ERG double-flash paradigm. In this experimental setup, the retina receives a conditioning flash, followed by a probing flash of the same light intensity (Figure 4A1). The b-wave amplitude ratio of probing to conditioning response in relation to the interstimulus interval is a normalized read-out for the visual transduction recovery time (Figure 4A2; full example in Figure 4—figure supplement 1). Photoreceptor recovery is complete when the two flashes evoke responses of equal amplitudes. ERG responses are predicted to be contributed by all cone subtypes, given the light source spectrum.
Response recovery was significantly delayed in the evening in comparison to the morning (Figure 4A2). However, as the ERG b-wave is only an indirect measure of the photoreceptor response, we also measured the photoreceptor-induced a-wave by blocking the masking ERG b-wave (Figure 4B1). This was achieved by administering a pharmacological cocktail containing the excitatory amino acid transporter inhibitor DL-threo-beta-benzyloxyaspartate (DL-TBOA) and metabotropic glutamate receptor inhibitor L-2-amino-4-phosphonobutyric acid (L-AP4) (Wong et al., 2004). Consistently, the double-flash paradigm demonstrated that the a-wave response recovery in the evening was delayed (Figure 4B2). According to the light spectrum (Figure 4—figure supplement 2), the a-wave was contributed by all cone subtypes.
In order to prove that increased response recovery times measured in the evening are a bonafide circadian event, we repeated the above experiments on larvae that were kept in constant darkness. At corresponding time points, the decrease in response recovery was comparable (Figure 4C1&C2), verifying that the observed changes are regulated by an intrinsic circadian clock.
As photoresponse recovery is affected by the circadian rhythm, we hypothesized that this should also be apparent in temporal aspects of vision. Therefore, we recorded ERG responses generated by the flickering stimuli with different stimulus frequencies (Figure 5, 5 Hz, 8 Hz, 10 Hz, 12 Hz, and 15 Hz). Fast Fourier transform (FFT) algorithm in MATLAB was used to extract the power at stimulus frequency. This power was then normalized against the power at 50 Hz (line noise), which is far from the stimulus frequencies. In line with our hypothesis, we found that the normalized power at each stimulus frequency was significantly weaker in the evening compared with the power in the morning. This clearly indicates that the cone visual temporal resolution is under circadian control. Note here, the flicker ERG was mainly contributed by double-cone responses because of the spectral content of the stimulus light (Figure 4—figure supplement 2).
Manipulation of gene expression by light is mirrored by functional changes
Next we measured larvae reared in a reversed light cycle (DL) where the night turns into a day. Under this condition, gene expression levels stayed in the fish’s time. ISH for the genes of interest (Figure 6A) reflected this, with a stronger staining intensity in LD fish at 9 o’clock in the morning compared to DL fish at the same time. Consequently, when both groups were recorded at 120 hr post fertilization, a prolonged response recovery time was obtained in the fish maintained in reversed light cycle, reflecting the situation in fish kept in the normal light and recorded in the evening (Figure 6D).
While the intrinsic circadian clock is maintained in the absence of light, continuous light exposure has been shown to disrupt this intrinsic rhythm (Laranjeiro and Whitmore, 2014). We therefore evaluated if the circadian regulation of mRNA expression persists in larvae kept under constant light (LL). Strikingly, the gene expression differences between morning and evening detected under LD conditions were completely lost in LL larvae (Figure 6B&C). This was also reflected on a functional level with no delay of photoresponse recovery in the evening, as measured by ERG.
Taken together, these results demonstrate that changes in the light cycle are reflected in changes of transcript levels of phototransduction regulators that subsequently lead to altered visual performance at different times during the day.
Circadian clock-dependent expressions of key regulator genes tune the single-cone photoresponse kinetics
We applied a computational model of visual transduction to predict how the relative gene expression changes between morning and evening influence the single-cone photoresponse (Invergo et al., 2013; Invergo et al., 2014). The default model was set as morning value (ZT1). We then put the measured gene expression ratio data (arr3a, grk7a, rcv2b and rgs9) between ZT1 and ZT13 into the model for evening simulation. These four genes have been selected due to their pan-cone expression (grk7a, rcv2b and rgs9) and double-cone expression (arr3a), respectively. Running the model with the relative value of arr3b (blue and ultraviolet [UV] cones) produced comparable results to arr3a (data not shown). Detailed parameters are listed in Supplementary file 4. The computed morning and evening values were then compared.
As predicted by our experimental results, the decay of photoresponse to different light intensities in the model was largely prolonged in the evening (Figure 7A–E). The unsaturating response amplitude was slightly elevated in the evening, which may indicate the prolonged lifetime of the visual pigment (Figure 7F).
Discussion
Circadian rhythms have been shown to regulate many biological aspects of vision. An early study demonstrated that zebrafish visual sensitivity is lower before light on and higher prior to light off (Li and Dowling, 1998). Later, another study linked the rhythmic expression of long-wavelength cone opsin to the core clock component CLOCK (Li et al., 2008). A particularly striking finding showed that synaptic ribbons of larval zebrafish photoreceptors disassemble at night. This peculiar phenomenon may save energy in fast-growing larvae (Emran et al., 2010). Our study now demonstrates that regulators of photoresponse decay are not only influenced by the circadian clock but in addition have a clear effect on the varying visual performances throughout a 24 hr cycle. Moreover, kinetics of cone visual transduction quenching is under the control of the circadian clock, which allows the fish to see with better temporal resolution in the morning than in the evening.
It is commonly assumed that circadian gene regulation helps the organism to optimally adapt to its preferential lifestyle and/or environment. Therefore, one would expect that the circadian systems of diurnal and nocturnal animals adapt differently. Our study indeed demonstrates that orthologous zebrafish and mouse genes involved in regulating cone visual transduction decay display an anti-phasic circadian expression pattern, supporting the functional relevance of the oscillating gene expression. While the visual temporal resolution of diurnal species is reduced in the evening, the visual system of nocturnal species is tuned to be most effective during these hours. Zebrafish, therefore, is an interesting model to study the physiology of circadian rhythms of diurnal animals, such as humans.
We would like to point out several additional interesting observations. Although many ohnologs (paralogs generated in a whole-genome duplication event), such as grk7a and grk7b, share a similar circadian phase or oscillatory amplitude, others, such as rcv1a and rcv1b, show an almost anti-phasic relationship. This is remarkable, since these ohnologs have been generated by a teleost-specific whole-genome duplication event (Glasauer and Neuhauss, 2014), implying that initially all ohnologs should have been in synchronicity. Interestingly, these ohnologs also adapted different expression profiles, with rcv1a being expressed in rods and UV cones, while rcv1b being expressed in all cone types in the adult retina (UV, blue, red, and green) (Zang et al., 2015).
While the circadian rhythmicity of most genes persists throughout all developmental stages, some genes do show markedly different expression profiles between larval and adult stages. This may be related to the fact that the larval retina is functionally cone dominant, while the adult retina is a duplex retina with rod and cone contribution. In the case of rcv2 ohnologs, rcv2b displays an in-phase cyclic expression pattern throughout all stages. Conversely, rcv2a did not show an overt cyclic expression pattern at larval stages, but being clearly under circadian control at adult stages (Figure 1). In contrast to the rcv1 ohnologs, both rcv2 genes are expressed in all cone subtypes, and depletion of either one acts to speed up the photoresponse termination (Zang et al., 2015). Other examples of ohnolog-specific cycling have been found for arrs and rgs genes (Figure 1, Figure 1—figure supplement 4). These observations strongly indicate that the transcription of clock-controlled genes (CCGs) is not uniformly regulated.
Interestingly, it has been previously demonstrated that the circadian clock seems to be desynchronized in larvae raised in darkness (Dekens and Whitmore, 2008; Kaneko and Cahill, 2005; Kazimi and Cahill, 1999; Lahiri et al., 2014). The circadian expression of some core clock genes and melatonin rhythms are lost when whole larvae were used as the experimental material in the absence of environmental entrainments. We did not observe this phenomenon in our study of visual transduction genes when only eye tissue was used, consistent with an inheritable maternal clock in the eye. We took care to avoid inadvertent environmental entrainment as described in detail in the ‘Materials and methods’ section. The different experimental results may come from the fact that different experiment materials were used. For example, all the analyzed genes in our study are also expressed in the photoreceptors of the pineal gland, but the transcript fluctuations may not necessarily be synchronized between the eye and pineal gland (eg, rcv1a in Figure 1—figure supplement 1). The use of whole larvae in our qRT-PCR study may have masked the cycling of retinal genes.
Furthermore, in many cases, the DD cycling is in phase with the fluctuations of the transcripts under the LD cycle. Endogenous circadian periods are around, but not exactly, 24 hr, and within 5 days in constant darkness, the peaks may shift relative to the LD cycle. In our experiments, the tissue was collected every 3 hr, and the shift within 5 days may be too small to be visible in the current experimental setting. The observation that the cycle of other genes (arr3a, arr3b, and rcv1a) in DD condition did diverge from LD condition indicates that these genes may be driven by different transcription factors. Furthermore, the DD condition led to the overall upregulation of some (eg, rcv1a, rgs9a) while caused downregulation of other (eg, arr3a, arr3b) genes as calculated from data in Figure 1. This strongly argues against a systematic error.
Among the studied genes in zebrafish, grk7a expression level increased by around 50 times in 1 day (Figure 1A), whereas Grk7a protein level increased by about two times in a 24 hr period (Figure 3A). arr3a transcript increased about 10 times (Figure 1C), while its protein level only grew less than 50% throughout the day (Figure 3A). Therefore, these mRNA expression levels reflect proportionally to protein levels, indicative of a rather fast turnover rate for these proteins.
In the end, we asked whether the observed ERG adaptations between morning and evening directly influence visual behavior. Therefore, we measured the optokinetic response (OKR)(Figure 7—figure supplement 1) and the visual motor response (VMR)(Figure 7—figure supplement 2) . Both behavioral assays showed some changes between the different recording time points, but the direct contribution by visual transduction is hard to assign. Confounding factors, not related to vision, may for instance be circadian regulation of overall activity.
In conclusion, we have shown that key regulators of cone visual transduction at both the mRNA and protein level are under circadian control. Moreover, expression levels of these regulators in diurnal and nocturnal species are anti-phasic, suggesting that circadian changes influencing physiological and behavioral properties of vision are reflected in adaptation to different visual ecologies.
Materials and methods
Zebrafish care
Request a detailed protocolZebrafish (Danio rerio) were maintained at a standard 14 hr light:10 hr dark cycle (LD) with light on at 8 am and light off at 10 pm. Water temperatures were kept between 26 and 28°C (Amores et al., 1998). Fish from the WIK wildtype strain were used in our study. Embryos were raised in E3 medium (5 mM NaCl, 0.17 mM KCl, 0.33 mM CaCl2, and 0.33 mM MgSO4) containing either 0.01% methylene blue to suppress fungal growth and/or 0.2 mM 1-phenyl-2-thiourea (PTU; Sigma-Aldrich) to prevent pigment development. Embryos were collected directly after laying. LD condition embryos were then transferred to the incubation room with normal light cycle (14:10). DD-conditioned embryos were placed in a black box before being transferred to the incubation room. Hence, all larvae (LD and DD) grew in the same environment with a stable temperature at 28°C. LL-conditioned fish were raised under constant light. DL condition was light on at 8 pm and light off at 10 am.
Adult zebrafish were sacrificed using ice water following decapitation. All animal experiments were carried out in line with the ARVO Statement for the Use of Animals in Ophthalmic and Vision Research and were approved by the Veterinary Authorities of Kanton Zurich, Switzerland (TV4206).
Zebrafish quantitative real-time PCR
Request a detailed protocolAround thirty 5 dpf larvae or five eyeballs from adult zebrafish were collected per time point (ZT1, 4, 7, 10, 13, 16, 19 and 22) and the tissue stored in RNAlater (Sigma) at 4°C. Dark adapted tissue was collected under dim red light. Only eyeballs were used for RNA extraction using the NucleoSpin RNA kit (Macherey-Nagel). Complementary DNA (cDNA) was produced using 110 ng total RNA as template for reverse transcription with SuperScript III (Invitrogen, Life Technologies; Zug, Switzerland). The samples collected from different time points were masked during RNA extraction and cDNA generation. qRT-PCR (Applied Biosystems Prism SDS 7900HT; Life Technologies) was performed using the MESA Green qPCR Mastermix Plus for SYBR Assay (Eurogentec, Seraing, Belgium) on a liquid handling robot platform (Tecan Genesis). Three technical replicates were conducted. Primers (Sigma-Aldrich) for qRT-PCR were intron-spanning to avoid amplification of non-digested genomic DNA fragments and were designed by online Universal ProbeLibrary Assay Design Center (Roche). Standard housekeeping genes (elongation factor 1, ef1; β-actin 2, actb2 and ribosomal protein L 13, rpl13) were used as reference (Tang et al., 2007). Primer pairs used are listed in Table 1.Expression levels were normalized to 1. Statistical analysis was performed in R 4.1.0 with ‘rain’ package (Thaben and Westermark, 2014).
Mouse care and gene expression analysis
Request a detailed protocolMice were maintained at the Laboratory Animal Services Center (LASC) of the University of Zurich in a 12 hr light:12 hr dark cycle with lights on at 7 am. All animal experiments were performed according to the ARVO Statement for the Use of Animals in Ophthalmic and Vision Research and the regulations of Veterinary Authorities of Kanton Zurich, Switzerland.
Ten 12-week-old wildtype mice (129S6; Taconic, Ejby, Denmark) were used in our experiments. Dark-phase mice were killed under red light and retinas were processed further under normal light conditions. Three mice at each time point (ZT1, 5, 9, 13, 17 and 21) were sacrificed and RNA was extracted (Macherey-Nagel, Oensingen, Switzerland) according to the manufacturer’s instructions. cDNA synthesized using oligo-dT was done as previously described (Storti et al., 2019). The samples collected from different time points were masked during RNA extraction and cDNA generation. qRT-PCR was performed by ABI QuantStudio3 machine (Thermo Fisher Scientific) with the PowerUp Sybr Green master mix (Thermo Fisher Scientific). Two technical replicates were conducted. Primer pairs used are listed in Table 2 for each gene of interest. Beta-actin (Actb) was used as a housekeeping gene to normalize gene expression with the comparative threshold cycle method (DDCt) using the Relative Quantification software (Thermo Fisher Scientific). The highest expression level was normalized to 1. Statistical analysis was performed in R 4.1.0 with ‘rain’ package (Thaben and Westermark, 2014).
In situ hybridization
Request a detailed protocolPrimers used to generate in situ probes are listed in Table 3. Probes were digoxigenin‐labeled using the DIG RNA Labeling Mix purchased from Roche.
For whole-mount ISH, embryos were treated with E3 containing 0.2 mM PTU (Sigma-Aldrich) to avoid pigmentation. 5 dpf larvae were fixed in 4% paraformaldehyde (PFA; Sigma) in phosphate-buffered saline (PBS) overnight at 4°C. Time points with maximal differences were chosen according to qRT-PCR results. Embryos were washed three times in PBS containing 1% Tween (PBST), dehydrated step wise (25, 50, and 70% methyl alcohol (MeOH) in PBST), and stored in 100% MeOH at –20°C. When comparing two groups of samples fixed at different time points, the tails of the group that may produce weaker staining were cut and mixed with the other group during staining.
For slide ISH, eyeballs were removed from adult zebrafish at different time points and fixed overnight at 4°C using 4% PFA. Detailed ISH processes have been previously described (Haug et al., 2015). When comparing two groups of samples fixed at different time points, both samples were placed on the same slide.
Infrared western blotting
View detailed protocolFive to six eyeballs from adult zebrafish were homogenized in ice-cold 150 ml RIPA buffer (150 mM NaCl, 1% Triton-X, 0.5% sodiumdeoxycholate, 50 mM Tris (pH 8), 1 mM ethylenediaminetetraacetic acid [EDTA], 0.1% sodium dodecyl sulfate [SDS]) containing cOmplete Protease Inhibitor Cocktail ([Roche]). After 2 hr of incubation on a 4°C shaker, lysates were centrifuged for 30 min at 4°C. During this procedure, all the samples were masked. Supernatants were stored at –80°C. Nitrocellulose membranes with 0.45 µm pore size were used. Primary antibodies were diluted to the following concentrations: rabbit anti-Arr3a: 1:4000; rabbit anti-Grk7a: 1:3000; mouse anti-β-actin: 1:6000 (Renninger et al., 2011; Rinner et al., 2005). Anti-arr3a and anti-β-actin antibodies or anti-Grk7a and anti-β-actin antibodies were applied simultaneously. Secondary antibodies IRDye 800CW Goat anti-Rabbit IgG and IRDye 680RD Goat anti-Mouse IgG (LI-COR) were diluted in 1:20,000 ratio in blocking buffer (1% bovine serum albumin [BSA] in PBST). Signal was detected by the Odyssey CLx Imaging System (LI-COR) and data were normalized to the internal loading control β-actin by IMAGEJ (Schindelin et al., 2012).
Electroretinography
Request a detailed protocolERG was recorded as previously described (Zang et al., 2015). Light intensity (light source: Zeiss XBO 75 W) was measured using a spectrometer (Ocean Optics, USB2000b; software Spectra Suite, Ocean Optics) with a spectral range described previously (Supplemental Material 2A in Zang et al., 2015). Pairs of two light flashes with equal intensity and duration (500 ms) were applied (Rinner et al., 2005). Intervals between two flashes were either 100, 200, 300, 500, 1000, 2000, 3000, or 5000 ms. The interval between two pairs was 20 s. b-wave recovery is defined as the ratio of the second b-wave amplitude to the first one in the same pair.
To measure ERG a-wave, 5 dpf larval eyeballs were treated with 400 µM L-AP4 and 200 µM TBOA in Ringer’s solution (111 mM NaCl, 2.5 mM KCl, 1 mM CaCl2, 1.6 mM MgCl2, 10 μm EDTA as a chelator for heavy metal ions, 10 mM glucose, and 3 mM 4- (2-hydroxyethyl) -1-piperazineethanesulfonic acid [HEPES] buffer, adjusted to pH 7.7–7.8 with NaOH). A HPX-2000 Xenon light source (Ocean Optics) was used and its light spectrum was measured by a spectrometer (Ocean Optics, USB2000b; software Spectra Suite, Ocean Optics; Figure 4—figure supplement 2). Electronic signals were amplified 1000 times by a pre-amplifier (P55 AC Preamplifier; Astro-Med. Inc, Grass Technology), digitized by DAQ Board (SCC-68; National Instruments), and recorded by a self-written Labview program (National Instruments). Intervals between two flashes were 300 ms, 500 ms, 1000 ms, and 1500 ms, respectively. a-wave recovery is defined as the ratio of the second a-wave amplitude to the first one in the same pair.
Flicker-fusion ERGs were measured with a white light emitting diode (LED) light source (Ocean Optics; LSM serie) controlled by LDC-1 controller (Ocean Optics). The spectrum of this light source was was measured by a spectrometer (Ocean Optics, USB2000b; software Spectra Suite, Ocean Optics; Figure 4—figure supplement 2). Except for the light source, flicker ERG was performed in the same setup as a-wave ERG. The flicker frequencies of 5 Hz, 8 Hz, 10 Hz, 12 HZ, and 15 Hz at 50% duty cycle were used. Flicker-fusion ERG data were analyzed by MATLAB (R2020b).
Phototransduction modeling
Request a detailed protocolThe computational model of vertebrate phototransduction was introduced and verified previously (Invergo et al., 2014; Invergo et al., 2013). We simulated the photoresponse to different light intensities of 1.7, 4.8, 15.2, 39.4, 125, 444, 1406, and 4630 photons µm–2 with a flash duration of 500 ms. Default parameters in the model were kept for morning (ZT1) simulation. For evening (ZT13) simulation, the relative gene expression change between ZT1 and ZT13 of larvae LD conditions was applied. Parameters for each gene are listed in Supplementary file 4. The simulation was performed in COPASI (Hoops et al., 2006).
Visual motor response
Request a detailed protocolThe VMR was measured using a Zebrabox (ViewPoint Life Science, Lyon, France). 5 dpf larvae were placed in a 96-well plate, subjected to dark adaptation for 10 min inside the Zebrabox, and the larval movement recorded with light off, on, and off for 5 min each. The distance that a single larva moved was measured every 2 s. Baseline activity was calculated as the average movement 1 min before light on or off.
Optokinetic response
Request a detailed protocolThe OKR was recorded as previously described (Rinner et al., 2005). Briefly, 5 dpf larvae were tested with sinusoidal gratings at different time points (ZT1, 4, 7, 10 and 13). To determine the contrast sensitivity, a spatial frequency of 20 cycles/360° and an angular velocity of 7.5°/s were used with different contrast settings (5, 10, 20, 40, 70, and 100%). To explore the spatial sensitivity, an angular velocity of 7.5°/s and 70% of maximum contrast were applied with a varying spatial frequency (7, 14, 21, 28, 42, and 56 cycles/360°). Figures were prepared by SPSS (version 23.0; Armonk, NY: IBM Corp).
Appendix 1
Data availability
All data generated and analysed during this study are included in the manuscript and supporting files. The dataset has been uploaded to dryad at https://doi.org/10.5061/dryad.0cfxpnw26.
-
Dryad Digital RepositoryCircadian Regulation of Vertebrate Cone Photoreceptor Function.https://doi.org/10.5061/dryad.0cfxpnw26
References
-
Ontogeny of cone photoreceptor mosaics in zebrafishThe Journal of Comparative Neurology 518:4182–4195.https://doi.org/10.1002/cne.22447
-
Rod contributions to the electroretinogram of the dark-adapted developing zebrafishDevelopmental Dynamics 222:564–570.https://doi.org/10.1002/dvdy.1188
-
Peripheral Circadian OscillatorsThe Yale Journal of Biology and Medicine 92:327–335.
-
Activation, deactivation, and adaptation in vertebrate photoreceptor cellsAnnual Review of Neuroscience 24:779–805.https://doi.org/10.1146/annurev.neuro.24.1.779
-
Clock mechanisms in zebrafishCell and Tissue Research 309:27–34.https://doi.org/10.1007/s00441-002-0570-7
-
Modulation of mouse rod response decay by rhodopsin kinase and recoverinThe Journal of Neuroscience 32:15998–16006.https://doi.org/10.1523/JNEUROSCI.1639-12.2012
-
Coordination between production and turnover of interphotoreceptor retinoid-binding protein in zebrafishInvestigative Ophthalmology & Visual Science 41:3590–3599.
-
Autonomous onset of the circadian clock in the zebrafish embryoThe EMBO Journal 27:2757–2765.https://doi.org/10.1038/emboj.2008.183
-
Development of a rod photoreceptor mosaic revealed in transgenic zebrafishDevelopmental Biology 258:277–290.https://doi.org/10.1016/s0012-1606(03)00125-8
-
Adaptation in vertebrate photoreceptorsPhysiological Reviews 81:117–151.https://doi.org/10.1152/physrev.2001.81.1.117
-
Phototransduction in mouse rods and conesPflugers Archiv 454:805–819.https://doi.org/10.1007/s00424-006-0194-y
-
Circadian regulation in the retina: From molecules to networkThe European Journal of Neuroscience 51:194–216.https://doi.org/10.1111/ejn.14185
-
Whole-genome duplication in teleost fishes and its evolutionary consequencesMolecular Genetics and Genomics 289:1045–1060.https://doi.org/10.1007/s00438-014-0889-2
-
Minutes, days and years: molecular interactions among different scales of biological timingPhilos Trans R Soc Lond B Biol Sci 369:20120465.https://doi.org/10.1098/rstb.2012.0465
-
Eumetazoan cryptochrome phylogeny and evolutionGenome Biology and Evolution 7:601–619.https://doi.org/10.1093/gbe/evv010
-
COPASI--a COmplex PAthway SImulatorBioinformatics 22:3067–3074.https://doi.org/10.1093/bioinformatics/btl485
-
Circadian clocks: lessons from fishProgress in Brain Research 199:41–57.https://doi.org/10.1016/B978-0-444-59427-3.00003-4
-
Exploring the rate-limiting steps in visual phototransduction recovery by bottom-up kinetic modelingCell Communication and Signaling 11:36.https://doi.org/10.1186/1478-811X-11-36
-
A comprehensive model of the phototransduction cascade in mouse rod cellsMolecular BioSystems 10:1481–1489.https://doi.org/10.1039/c3mb70584f
-
Rod electrical coupling is controlled by a circadian clock and dopamine in mouse retinaThe Journal of Physiology 593:1597–1631.https://doi.org/10.1113/jphysiol.2014.284919
-
Development of a circadian melatonin rhythm in embryonic zebrafishBrain Research. Developmental Brain Research 117:47–52.https://doi.org/10.1016/s0165-3806(99)00096-6
-
Developmental stage-specific regulation of the circadian clock by temperature in zebrafishBioMed Research International 2014:930308.https://doi.org/10.1155/2014/930308
-
Phototransduction, dark adaptation, and rhodopsin regeneration the proctor lectureInvestigative Ophthalmology & Visual Science 47:5137–5152.https://doi.org/10.1167/iovs.06-0849
-
Zebrafish visual sensitivity is regulated by a circadian clockVisual Neuroscience 15:851–857.https://doi.org/10.1017/s0952523898155050
-
CLOCK is required for maintaining the circadian rhythms of Opsin mRNA expression in photoreceptor cellsThe Journal of Biological Chemistry 283:31673–31678.https://doi.org/10.1074/jbc.M803875200
-
Photoreceptor coupling is controlled by connexin 35 phosphorylation in zebrafish retinaThe Journal of Neuroscience 29:15178–15186.https://doi.org/10.1523/JNEUROSCI.3517-09.2009
-
Photopigment quenching is Ca2+ dependent and controls response duration in salamander L-cone photoreceptorsThe Journal of General Physiology 135:355–366.https://doi.org/10.1085/jgp.200910394
-
Diurnal and circadian retinomotor movements in zebrafishVisual Neuroscience 22:203–209.https://doi.org/10.1017/S0952523805222083
-
Cone arrestin confers cone vision of high temporal resolution in zebrafish larvaeThe European Journal of Neuroscience 33:658–667.https://doi.org/10.1111/j.1460-9568.2010.07574.x
-
Fiji: an open-source platform for biological-image analysisNature Methods 9:676–682.https://doi.org/10.1038/nmeth.2019
-
The interphotoreceptor retinoid-binding protein (IRBP) of the chicken (Gallus gallus domesticusMolecular Vision 11:833–845.
-
Detecting rhythms in time series with RAINJournal of Biological Rhythms 29:391–400.https://doi.org/10.1177/0748730414553029
-
Systems biology of mammalian circadian clocksAnnual Review of Physiology 72:579–603.https://doi.org/10.1146/annurev-physiol-073109-130051
-
It’s time to swim! Zebrafish and the circadian clockFEBS Letters 585:1485–1494.https://doi.org/10.1016/j.febslet.2011.04.007
-
Circadian and light-driven regulation of rod dark adaptationScientific Reports 5:17616.https://doi.org/10.1038/srep17616
-
Origin and control of the dominant time constant of salamander cone photoreceptorsThe Journal of General Physiology 140:219–233.https://doi.org/10.1085/jgp.201110762
-
The Binding Properties and Physiological Functions of RecoverinFrontiers in Molecular Neuroscience 11:473.https://doi.org/10.3389/fnmol.2018.00473
Article and author information
Author details
Funding
Schweizerischer Nationalfonds zur Förderung der Wissenschaftlichen Forschung (310030_200376)
- Marijana Samardzija
The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.
Acknowledgements
We would like to thank Kara Kristiansen and Martin Walther for expert animal maintenance.
Ethics
All animal experiments were carried out in the line with the ARVO Statement for the Use of Animals in Ophthalmic and Vision Research and were approved by the Veterinary Authorities of Kanton Zurich, Switzerland (TV4206).
Copyright
© 2021, Zang et al.
This article is distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use and redistribution provided that the original author and source are credited.
Metrics
-
- 1,394
- views
-
- 286
- downloads
-
- 7
- citations
Views, downloads and citations are aggregated across all versions of this paper published by eLife.
Download links
Downloads (link to download the article as PDF)
Open citations (links to open the citations from this article in various online reference manager services)
Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)
Further reading
-
- Computational and Systems Biology
- Neuroscience
Diffusional kurtosis imaging (DKI) is a methodology for measuring the extent of non-Gaussian diffusion in biological tissue, which has shown great promise in clinical diagnosis, treatment planning, and monitoring of many neurological diseases and disorders. However, robust, fast, and accurate estimation of kurtosis from clinically feasible data acquisitions remains a challenge. In this study, we first outline a new accurate approach of estimating mean kurtosis via the sub-diffusion mathematical framework. Crucially, this extension of the conventional DKI overcomes the limitation on the maximum b-value of the latter. Kurtosis and diffusivity can now be simply computed as functions of the sub-diffusion model parameters. Second, we propose a new fast and robust fitting procedure to estimate the sub-diffusion model parameters using two diffusion times without increasing acquisition time as for the conventional DKI. Third, our sub-diffusion-based kurtosis mapping method is evaluated using both simulations and the Connectome 1.0 human brain data. Exquisite tissue contrast is achieved even when the diffusion encoded data is collected in only minutes. In summary, our findings suggest robust, fast, and accurate estimation of mean kurtosis can be realised within a clinically feasible diffusion-weighted magnetic resonance imaging data acquisition time.
-
- Neuroscience
The retina transforms patterns of light into visual feature representations supporting behaviour. These representations are distributed across various types of retinal ganglion cells (RGCs), whose spatial and temporal tuning properties have been studied extensively in many model organisms, including the mouse. However, it has been difficult to link the potentially nonlinear retinal transformations of natural visual inputs to specific ethological purposes. Here, we discover a nonlinear selectivity to chromatic contrast in an RGC type that allows the detection of changes in visual context. We trained a convolutional neural network (CNN) model on large-scale functional recordings of RGC responses to natural mouse movies, and then used this model to search in silico for stimuli that maximally excite distinct types of RGCs. This procedure predicted centre colour opponency in transient suppressed-by-contrast (tSbC) RGCs, a cell type whose function is being debated. We confirmed experimentally that these cells indeed responded very selectively to Green-OFF, UV-ON contrasts. This type of chromatic contrast was characteristic of transitions from ground to sky in the visual scene, as might be elicited by head or eye movements across the horizon. Because tSbC cells performed best among all RGC types at reliably detecting these transitions, we suggest a role for this RGC type in providing contextual information (i.e. sky or ground) necessary for the selection of appropriate behavioural responses to other stimuli, such as looming objects. Our work showcases how a combination of experiments with natural stimuli and computational modelling allows discovering novel types of stimulus selectivity and identifying their potential ethological relevance.