Drosophila PSI controls circadian period and the phase of circadian behavior under temperature cycle via tim splicing
Abstract
The Drosophila circadian pacemaker consists of transcriptional feedback loops subjected to post-transcriptional and post-translational regulation. While post-translational regulatory mechanisms have been studied in detail, much less is known about circadian post-transcriptional control. Thus, we targeted 364 RNA binding and RNA associated proteins with RNA interference. Among the 43 hits we identified was the alternative splicing regulator P-element somatic inhibitor (PSI). PSI regulates the thermosensitive alternative splicing of timeless (tim), promoting splicing events favored at warm temperature over those increased at cold temperature. Psi downregulation shortens the period of circadian rhythms and advances the phase of circadian behavior under temperature cycle. Interestingly, both phenotypes were suppressed in flies that could produce TIM proteins only from a transgene that cannot form the thermosensitive splicing isoforms. Therefore, we conclude that PSI regulates the period of Drosophila circadian rhythms and circadian behavior phase during temperature cycling through its modulation of the tim splicing pattern.
Introduction
Circadian rhythms are the organism’s physiological and behavioral strategies for coping with daily oscillations in environment conditions. Inputs such as light and temperature feed into a molecular clock via anatomical and molecular input pathways and reset it every day. Light is the dominant cue for entraining the molecular clock, but temperature is also a pervasive resetting signal in natural environments. Paradoxically, clocks must be semi-resistant to temperature: they should not hasten in warm summer months or lag in the winter cold (this is called temperature compensation), but they can synchronize to the daily rise and fall of temperature (temperature entrainment) (Pittendrigh, 1960). Not only can temperature entrain the clock, it also has a role in seasonal adaptation by affecting the phase of behavior (see for example Majercak et al., 1999).
Molecular circadian clocks in eukaryotes are made up of negative transcriptional feedback loops (Dunlap, 1999). In Drosophila, the transcription factors CLOCK (CLK) and CYCLE (CYC) bind to E-boxes in the promoters of the clock genes period (per) and timeless (tim) and activate their transcription. PER and TIM proteins accumulate in the cytoplasm where they heterodimerize and enter the nucleus to feedback and repress the activity of CLK and CYC and thus downregulate their own transcription (Hardin, 2011). This main loop is strengthened by a scaffolding of interlocked feedback loops involving the transcription factors vrille (vri), PAR domain protein 1 (Pdp1) and clockwork orange (cwo). Post-translational modifications are well-established mechanisms for adjusting the speed and timing of the clock (Tataroglu and Emery, 2015).
Increasing evidence indicates that post-transcriptional mechanisms controlling gene expression are also critical for the proper function of circadian clocks in many organisms. In Drosophila, the post-transcriptional regulation of per mRNA has been best studied. per mRNA stability changes as a function of time (So and Rosbash, 1997). In addition, per contains an intron in its 3’UTR (dmpi8) that is alternatively spliced depending on temperature and lighting conditions (Majercak et al., 1999; Majercak et al., 2004). On cold days, the spliced variant is favored, causing an advance in the accumulation of per transcript levels as well as an advance of the evening activity peak. This behavioral shift means that the fly is more active during the day when the temperature would be most tolerable in their natural environment. The temperature sensitivity of dmpi8 is due to the presence of weak non-canonical splice sites. However, the efficiency of the underlying baseline splicing is affected by four single nucleotide polymorphisms (SNPs) in the per 3’UTR that vary in natural populations and form two distinct haplotypes (Low et al., 2012; Cao and Edery, 2017). Also, while this splicing is temperature-sensitive in two Drosophila species that followed human migration, two species that remained in Africa lack temperature sensitivity of dmpi8 splicing, (Low et al., 2008). Furthermore, Zhang et al. (2018) recently demonstrated that the the trans-acting splicing factor B52 enhances dmpi8 splicing efficiency, and this effect is stronger with one of the two haplotypes. per is also regulated post-transcriptionally by the TWENTYFOUR-ATAXIN2 translational activation complex (Zhang et al., 2013; Lim et al., 2011; Lim and Allada, 2013a; Lee et al., 2017). This complex works by binding to per mRNA as well as the cap-binding complex and poly-A binding protein. This may enable more efficient translation by promoting circularization of the transcript. Interestingly, this mechanism appears to be required only in the circadian pacemaker neurons. Non-canonical translation initiation has also been implicated in the control of PER translation (Bradley et al., 2012). Regulation of PER protein translation has also been studied in mammals, with RBM4 being a critical regulator of mPER1 expression (Kojima et al., 2007). In flies however, the homolog of RBM4, LARK, regulates the translation of DBT, a PER kinase (Huang et al., 2014). miRNAs have emerged as important critical regulators of circadian rhythms in Drosophila and mammals, affecting the circadian pacemaker itself, as well as input and output pathways controlling rhythmic behavioral and physiological processes (Tataroglu and Emery, 2015; Lim and Allada, 2013b).
RNA-associated proteins (RAPs) include proteins that either bind directly or indirectly to RNAs. They mediate post-transcriptional regulation at every level. Many of these regulated events – including alternative splicing, splicing efficiency, mRNA stability, and translation – have been shown to function in molecular clocks. Thus, to obtain a broad view of the Drosophila circadian RAP landscape and its mechanism of action, we performed an RNAi screen targeting 364 of these proteins. This led us to discover a role for the splicing factor P-element somatic inhibitor (PSI) in regulating the pace of the molecular clock through alternative splicing of tim.
Results
An RNAi screen for RNA-associated proteins controlling circadian behavioral rhythms
Under constant darkness conditions (DD) flies have an intrinsic period length of about 24 hr. To identify novel genes that act at the post-transcriptional level to regulate circadian locomotor behavior, we screened 364 genes, which were annotated in either Flybase (FB2014_03, Thurmond et al., 2019) or the RNA Binding Protein Database (Cook et al., 2011) as RNA binding or involved in RNA associated processes, using period length as a readout of clock function (Supplementary file 1: RAP Screen Dataset). We avoided many, but not all, genes with broad effects on gene expression, such as those encoding essential splicing or translation factors. When possible, we used at least two non-overlapping RNAi lines from the TRiP and VDRC collections. RNAi lines were crossed to two different GAL4 drivers: tim-GAL4 (Kaneko et al., 2000) and Pdf-GAL4 (Renn et al., 1999) each combined with a UAS-dicer-2 transgene to enhance the strength of the knockdown (Dietzl et al., 2007). These combinations will be abbreviated as TD2 and PD2, respectively. tim-GAL4 drives expression in all cells with circadian rhythms in the brain and body (Kaneko et al., 2000), while Pdf-GAL4 drives expression in a small subset of clock neurons in the brain: the PDF-positive small (s) and large (l) LNvs (Renn et al., 1999). Among them, the sLNvs are critical pacemaker neurons that drive circadian behavior in DD (Renn et al., 1999; Stoleru et al., 2005). In the initial round of screening, we tested the behavior of 4–8 males for each RNAi line crossed to both TD2 and PD2 (occasionally, fewer males were tested if a cross produced little progeny). We also crossed some RNAi lines to w1118 (+) flies (most were lines selected for retest, see below). We noticed that RNAi/+ control flies for the TRiP collection were 0.3 hr shorter than those of the VDRC collection (Figure 1A). Furthermore, the mean period from all RNAi lines crossed to either PD2 or TD2 was significantly shorter for the TRiP collection than for the VDRC collection (Figure 1A) (0.2 hr, TD2 crosses; 0.5 hr, PD2 crosses). We also found that many of the VDRC KK lines that resulted in long period phenotypes when crossed to both drivers contained insertions in the 40D locus (VDRC annotation), although this effect was stronger with PD2 than TD2. It has been shown that this landing site is in the 5’UTR of tiptop (tio) and can lead to non-specific effects in combination with some GAL4 drivers, likely due to misexpression of tio (Vissers et al., 2016; Green et al., 2014). Indeed, when we crossed a control line that contains a UAS insertion at 40D (40D-UAS) to PD2, the progeny also had a ca. 0.6 hr longer period relative to the PD2 control (Figure 1B). Thus, in order to determine a cutoff for candidates to further investigate, we analyzed the data obtained in our screen from the TRiP, VDRC, and the 40D KK VDRC lines independently (Figure 1C). These data are represented in two overlaid histograms that show period distributions: one for the TD2 crosses (blue) and one for the PD2 crosses (magenta). We chose a cutoff of two standard deviations (SD) from the mean period length for each RNAi line set. RNAi lines were selected for repeat if knockdown resulted in period lengths above or below the 2-SD cutoff. We also chose to repeat a subset of lines that did not pass the cutoff but were of interest and showed period lengthening or shortening, as well as lines that were highly arrhythmic in constant darkness (DD) or had an abnormal pattern of behavior in a light-dark cycle (LD). After a total of three independent experiments, we ended up with 43 candidates (Table 1) that passed the period length cutoffs determined by the initial screen; 31 showed a long period phenotype, while 12 had a short period. One line showed a short period phenotype with PD2 but was long with TD2 (although just below the 2-SD cutoff). Although loss of rhythmicity was also observed in many lines (Supplementary file 1), we decided to focus the present screen on period alterations to increase the probability of identifying proteins that regulate the circadian molecular pacemaker. Indeed, a change in the period length of circadian behavior is most likely caused by a defect in the molecular pacemaker of circadian neurons, while an increase in arrhythmicity can also originate from disruption of output pathways, abnormal development of the neuronal circuits underlying circadian behavioral rhythms, or cell death in the circadian neural network, for example.
Among the 43 candidate genes (Tables 1 and 2), we noticed a high proportion of genes involved or presumed to be involved in splicing (17), including five suspected or known to impact alternative splicing. Perhaps not surprisingly, several genes involved in snRNP assembly were identified in our screen. Their downregulation caused long period phenotypes. We also noticed the presence of four members of the CCR4-NOT complex, which can potentially regulate different steps of mRNA metabolism, including deadenylation, and thus mediate translational repression. Their downregulation mostly caused short period phenotypes and tended to result in high levels of arrhythmicity. Rga downregulation, however, resulted in a long period phenotype, suggesting multiple functions for the CCR4-NOT complex in the regulation of circadian rhythms. Interestingly, two genes implicated in mRNA decapping triggered by deadenylation, were also identified, with long periods observed when these genes were downregulated. Moreover, POP2, a CCR4-NOT component, was recently shown to regulate tim mRNA and protein levels (Grima et al., 2019). Another gene isolated in our screen, SMG5, was also recently found to impact circadian behavior (Ri et al., 2019). This validates our screen.
Knockdown of Psi shortens the period of behavioral rhythms
A promising candidate to emerge from our screen was the alternative splicing regulator PSI (Labourier et al., 2001; Siebel et al., 1992). Knockdown of Psi with both TD2 and PD2 crossed to two non-overlapping RNAi lines from the VDRC collection (GD14067 and KK101882) caused a significant period shortening, compared to the TD2/+ and PD2/+ controls (Figure 2A–E, Table 3), which the experimental flies need to be compared to since the GAL4 drivers in the TD2 and PD2 combination cause a previously reported dominant ca. 0.8 hr period lengthening (Figure 2C, left panel (TD2/+ compared to w1118); Kaneko et al., 2000; Renn et al., 1999; Zhang and Emery, 2013; Zhang et al., 2013). Importantly, the RNAi lines did not cause period shortening on their own (Figure 2C left panel, Table 3). While most experiments were performed at 25°C, we noticed that at 30°C, TD2/+ control had a period of ca. 24 hr (Figure 2C). We could thus meaningfully compare TD2/RNAi flies to both RNAi/+ and TD2/+ control at that temperature. The period of the experimental flies was significantly shorter than both controls (Figure 2C). Two additional lines from the TRiP collection (JF01476 and HMS00140) also caused period shortening when crossed to TD2 (Table 1). Interestingly, HMS00140 targets only the Psi-RA isoform, indicating that the RA isoform is important for the control of circadian period (Figure 2A). Since four RNA lines caused a similar phenotype and only two of them partially overlapped (Figure 2A), we are confident that the period shortening was not caused by off-target effects. Moreover, both the KK101882 and GD14067 lines have been shown to efficiently downregulate Psi (Guo et al., 2016), and we confirmed by quantitative Real-Time PCR (qPCR) that the RNAi line KK101882, which gave the shortest period phenotype with TD2, significantly reduced Psi mRNA levels in heads (Figure 2—figure supplement 1). This line was selected for most of the experiments described below as it gave the strongest period phenotype.
The phenotype caused by Psi downregulation was more pronounced with TD2 than with PD2 (Figure 2C–D, Table 3). This was unexpected since the sLNvs - targeted quite specifically by PD2 - determine circadian behavior period in DD (Stoleru et al., 2005; Renn et al., 1999). This could happen if PD2 is less efficient at downregulating Psi in sLNvs than TD2, or if the short period phenotype is not solely caused by downregulation of Psi in the sLNvs. To distinguish between these two possibilities, we used Pdf-GAL80 combined with TD2 to inhibit GAL4 activity specifically in the LNvs (Stoleru et al., 2004), while allowing RNAi expression in all other circadian tissues. With this combination, we also observed a significant period shortening compared to TD2/+; Pdf-GAL80/+ controls, but the period shortening was not as pronounced as with TD2 (Figure 2E, Table 3). We therefore conclude that both the sLNvs and non-PDF cells contribute to the short period phenotype caused by Psi downregulation (see discussion).
Psi overexpression disrupts circadian behavior
Since we observed that downregulating Psi leads to a short period, we wondered whether overexpression would have an inverse effect and lengthen the period of circadian behavior. Indeed, when we overexpressed Psi by driving a UAS-Psi transgene (Labourier et al., 2001) with the tim-GAL4 (TG4) driver, the period length of circadian behavior increased significantly by about 0.7 hr compared to the TG4/+ control (Figure 2F, Table 3). Interestingly, we also observed a severe decrease in the number of rhythmic flies. When we overexpressed Psi with Pdf-GAL4 (PG4), period was not statistically different from control (PG4/+), and rhythmicity was not reduced compared to the UAS-Psi/+ control (Figure 2G). Overexpression of Psi with the tim-GAL4; Pdf-GAL80 combination caused a severe decrease in rhythmicity but caused only a subtle period lengthening compared to TG4/+; Pdf-GAL80/+ controls (Figure 2H, Table 3). The effect of Psi overexpression on period is in line with the knockdown results, indicating that PSI regulates circadian behavioral period through both PDF+ LNvs and non-PDF circadian neurons. However, the increase in arrhythmicity observed with Psi overexpression is primarily caused by non-PDF cells.
Psi downregulation also shortens the period of body clocks
We wanted to further examine the effect of Psi knockdown on the molecular rhythms of two core clock genes: period (per) and timeless (tim). To do this, we took advantage of two luciferase reporter transgenes. We downregulated Psi with the TD2 driver in flies expressing either a TIM-LUCIFERASE (ptim-TIM-LUC) or a PER-LUCIFERASE (BG-LUC) fusion protein under the control of the tim or per promoter, respectively. We estimated period of luciferase activity rhythms over the first two days in DD, because oscillations rapidly dampened. Fully consistent with our behavioral results, the period of LUC activity was significantly shortened by about 1–1.5 hr compared to controls when Psi was downregulated in ptim-TIM-LUC flies (Figure 2—figure supplement 2A and B). Knockdown of Psi in BG-LUC flies resulted in a similar trend, although differences did not reach statistical significance (Figure 2—figure supplement 2C and D). Period was however shorter in experimental flies compared to both control genotypes in all four independent experiments performed with BG-LUC (and all six with ptim-TIM-LUC). Since the luciferase signal in these flies is dominated by light from the abdomen (Lamba et al., 2018; Stanewsky et al., 1997), this indicates that Psi knockdown, shortens the period of circadian clocks in peripheral tissues as well as in the brain neural network that controls circadian behavior.
Alternative splicing of two clock genes, cwo and tim, is altered in Psi knockdown flies
PSI has been best characterized for its role in alternative splicing of the P element transposase gene in somatic cells (Labourier et al., 2001; Siebel et al., 1992). However, it was recently reported that PSI has a wider role in alternative splicing (Wang et al., 2016). Wang et al. reported an RNA-seq dataset of alternative splicing changes that occur when a lethal Psi-null allele is rescued with a copy of Psi in which the AB domain has been deleted (PSIΔAB). This domain is required for the interaction of PSI with the U1 snRNP, which is necessary for PSI to mediate alternative splicing of P element transposase (Labourier et al., 2002). Interestingly, Wang et al. (2016) found that PSIΔAB affects alternative splicing of genes involved in complex behaviors such as learning, memory and courtship. Intriguingly, we found four core clock genes listed in this dataset: tim, cwo, sgg and Pdp1. We decided to focus on cwo and tim, since only one specific splicing isoform of Pdp1 is involved in the regulation circadian rhythm, (Pdp1e) (Zheng et al., 2009), and since the sgg gene produces a very complex set of alternative transcripts. After three days of LD entrainment, we collected RNA samples at four time points on the first day of DD and determined the relative expression of multiple isoforms of cwo and tim in Psi knockdown heads compared to driver and RNAi controls.
CWO is a basic helix-loop-helix (bHLH) transcriptional factor and is part of an interlocked feedback loop that reinforces the main loop by competing with CLK/CYC for E-box binding (Matsumoto et al., 2007; Lim et al., 2007; Kadener et al., 2007; Richier et al., 2008). There are three mRNA isoforms of cwo predicted in Flybase (Figure 3—figure supplement 1A) (Thurmond et al., 2019). Of the three, only cwo-RA encodes a full-length CWO protein. Exon two is skipped in cwo-RB, and in cwo-RC there is an alternative 3’ splice site in the first intron that lengthens exon 2. Translation begins from a downstream start codon in cwo-RB and -RC, because exon two skipping or lengthening, respectively, causes a frameshift after the start codon used in cwo-RA. The predicted start codon in both cwo-RB and cwo-RC would produce an N-terminal truncation of the protein, which would thus be missing the basic region of the bHLH domain and should not be able to bind DNA. The cwo-RB and cwo-RC isoforms may therefore encode endogenous dominant negatives.
We found that the level of the cwo-RB isoform was significantly reduced compared to both controls at CT 9 (Figure 3—figure supplement 1C). The cwo-RA isoform was also reduced compared to both controls at CT9 (Figure 3—figure supplement 1B). This reduction was significant compared to the TD2/+ control (p=0.0002) but was just above the significance threshold compared to the PsiRNAiKK/+ control (p=0.0715). Conversely, cwo-RC isoform expression was significantly increased at CT 15 (Figure 3—figure supplement 1D). The overall expression of all cwo mRNAs in Psi knockdown fly heads was significantly reduced at both CT 9 and CT 15, indicating that the RC isoform’s contribution to total cwo mRNA levels is quite modest (Figure 3—figure supplement 1E).
We then analyzed alternative splicing of tim in Psi knockdown heads compared to controls. Specifically, we looked at the expression of three temperature-sensitive intron inclusion events in tim that all theoretically lead to C-terminal truncations of the protein (Figure 3A). The tim-cold isoform, which is not annotated in Flybase (Thurmond et al., 2019), is dominant at low temperature (18°C) and arises when the last intron is retained (Boothroyd et al., 2007). We found that tim-cold is elevated in Psi knockdown heads at peak levels under 25°C conditions (CT15, Figure 3D). Similarly, we found that another intron inclusion event, tim-sc (tim-short and cold) which has also been shown to be elevated at 18°C and is present in the tim-RN and -RO isoforms (Martin Anduaga et al., 2019), is significantly increased at 25°C in Psi knockdown heads at CT15 (Figure 3B). Thus, interestingly, two intron inclusion events that are upregulated by cold temperature are also both upregulated in Psi knockdown heads at 25°C. In contrast, we found that an intron included in the tim-RM and -RS isoforms (tim-M, for tim-Medium) and shown to be increased at high temperature (29°C, Martin Anduaga et al., 2019; Shakhmantsir et al., 2018) is significantly decreased at CT 9, 15 and 21 in Psi knockdown heads at 25°C (Figure 3F). In the case of tim-sc, it should be noted that the intron is only partially retained, because a cleavage and poly-adenylation signal is located within this intron, thus resulting in a much shorter mature transcript (Martin Anduaga et al., 2019). Based on PSI function, the most parsimonious explanation is that PSI reduces production of tim-sc by promoting splicing of the relevant intron. However, we cannot entirely exclude that PSI regulates the probability of premature cleavage causing the RNA polymerase to undergo transcription termination soon after passing the poly-adenylation signal.
Collectively, these results indicate that, in wild-type flies, PSI shifts the balance of tim alternative splicing events toward a warm temperature tim RNA isoform profile at an intermediate temperature (25°C). This could be achieved either by altering the temperature sensitivity of tim introns, or by promoting a ‘warm temperature splicing pattern’ independently of temperature. We therefore also measured tim splicing isoforms at 18°C and 29°C (Figure 3C,E,G). We entrained flies for 3 days in LD at 25°C to maintain similar levels of GAL4 expression and thus of Psi knockdown (the GAL4/UAS system’s activity increases with temperature, Duffy, 2002). We then shifted them to either 18°C or 29°C at CT 0 on the first day of DD and collected samples at CT 3, 9, 15 and 21. We found that both the tim-cold intron and the tim-sc introns were elevated at 18°C in both Psi knockdown heads and controls (Figure 3C and E). Thus, Psi knockdown does not block the temperature sensitivity of these introns. tim-M levels were unexpectedly variable in DD, particularly in the Psi knockdown flies, perhaps because of the temperature change. Nevertheless, we observed a trend for the tim-M intron retention to be elevated at 29°C (Figure 3G), further supporting our conclusion that Psi knockdown does not affect the temperature sensitivity of tim splicing, but rather determines the ratio of tim mRNA isoforms, and it does this at all temperatures.
As expected from these results, Psi downregulation did not affect the ability of flies to adjust the phase of their evening and morning peak to changes in temperature (Figure 3—figure supplement 2). We also tested whether Psi knockdown flies responded normally to short light pulses, since TIM is the target of the circadian photoreceptor CRY (Emery et al., 1998; Stanewsky et al., 1998; Lin et al., 2001; Busza et al., 2004; Koh et al., 2006). These flies could both delay or advance the phase of their circadian behavior in response to early or late-night light pulses, respectively (Figure 3—figure supplement 3). We noticed however a possible slight shift of the whole Phase Response Curve toward earlier times. This would be expected since the pace of the circadian clock is accelerated.
PSI controls the phase of circadian behavior under temperature cycle
Since PSI regulates thermosensitive tim splicing events, we wondered whether it might have an impact on circadian behavioral responses to temperature. As mentioned above, Psi downregulation does not affect Drosophila’s ability to adjust the phase of their behavior to different constant ambient temperatures, under a LD cycle (Figure 3—figure supplement 2). Psi knockdown did not appear to affect temperature compensation, as these flies essentially responded to temperature in a similar way as their TD2/+ control, with shorter period at 29°C (Figure 4—figure supplement 1). However, we found a striking phenotype in flies with Psi downregulation under temperature cycle (29/20°C). Once flies had reached a stable phase relationship with the entraining temperature cycle (Busza et al., 2007), the phase of the evening peak of activity was advanced by about 2.5 hr in TD2/PsiRNAi, compared to controls, and this with two non-overlapping dsRNAs (Figure 4). Controls included TD2/+ or TD2/VIE-260B (KK host strain), RNAi/+, as well as TD2 crossed to a KK or GD RNAi line that did not produce a circadian phenotype. Importantly, no such phase advance was observed under LD (Figure 3—figure supplement 2), indicating that the short period phenotype does not account for the evening-peak advanced phase under temperature cycle. Rather, the phase advance is specific to temperature entrainment. The morning peak was difficult to quantify as it tended to be of low amplitude.
tim splicing is required for PSI’s regulation of circadian period and circadian behavior phase under temperature cycle
Because tim is a key element of the circadian transcriptional feedback loop and its splicing pattern is determined by the ambient temperature, we wondered whether PSI might be regulating the speed of the clock and the phase of the evening peak through its effects on tim splicing. We therefore rescued the amorphic tim allele (tim0) with a tim transgene that lacks the known temperature sensitive alternatively spliced introns as well as most other introns (timHA) (Figure 5A) (Rutila et al., 1998). Importantly, the tim0 mutation is a frame-shifting deletion located upstream of the temperature-sensitive alternative splicing events (Myers et al., 1995), and would thus truncate any TIM protein produced from the splice variants we studied. Strikingly, we found that knockdown of Psi in timHA rescued tim0 flies had no impact on the period of circadian behavior (Figure 5B–C, Table 3). Likewise, the evening peak phase under temperature cycles was essentially insensitive to Psi knockdown in timHA rescued tim0 flies (Figure 5D–E). This indicates that PSI controls circadian period in DD and the phase of the evening peak under temperature cycle through tim splicing.
Discussion
Our results identify a novel post-transcriptional regulator of the circadian clock: PSI. PSI is required for the proper pace of both brain and body clock, and for proper phase-relationship with ambient temperature cycles. When Psi is downregulated, the circadian pacemaker speeds up and behavior phase under temperature cycles is advanced by 3 hr, and these phenotypes appear to be predominantly caused by an abnormal tim splicing pattern. Indeed, the circadian period and behavior phase of flies that can only produce functional TIM protein from a transgene missing most introns is insensitive to Psi downregulation. We note however that cwo’s splicing pattern is also affected by Psi downregulation, and we did not study sgg splicing pattern, although it might also be controlled by PSI (Wang et al., 2016). We therefore cannot exclude a small contribution of non-tim splicing events to PSI downregulation phenotypes, or that in specific tissues these other splicing events play a greater role than in the brain.
Interestingly, Psi downregulation results in an increase in intron inclusion events that are favored under cold conditions (tim-sc and tim-cold), while an intron inclusion event favored under warm conditions is decreased (tim-M). However, the ability of tim splicing to respond to temperature changes is not abolished when Psi is downregulated (Figure 3C,E,G). This could imply that an as yet unknown factor specifically promotes or represses tim splicing events in a temperature-dependent manner. Another possibility is that the strength of splice sites or tim’s pre-mRNA structure impacts splicing efficiency in a temperature–dependent manner. For example, suboptimal per splicing signals explain the lower efficiency of per’s most 3’ splicing event at warm temperature (Low et al., 2008).
How would the patterns of tim splicing affect the pace of the circadian clock, or advance the phase of circadian behavior under temperature cycles? In all splicing events that we studied, intron retention results in a truncated TIM protein. It is therefore possible that the balance of full length and truncated TIM proteins, which may function as endogenous dominant-negatives, determines circadian period. For example, truncated TIM might be less efficient at protecting PER from degradation, thus accelerating the pacemaker, or affecting its phase. Consistent with this idea, overexpression of the shorter cold-favored tim isoform (tim-sc) shortens period (Martin Anduaga et al., 2019). Strikingly, Psi downregulation increases this isoform’s levels and also results in a short phenotype. Shakhmantsir et al. (2018) also proposed that production of tim-M transcripts (called tim-tiny in their study) delays the rate of TIM accumulation. Such a mechanism could also contribute to the short period we observed when Psi is downregulated, since this reduces tim-M levels, which may accelerate TIM accumulation. Another interesting question is how PSI differentially affects specific splice isoforms of tim. One possibility is that the execution of a specific tim splicing event negatively influences the probability of the occurrence of other splicing events. For example, PSI could downregulate tim-sc and tim-cold by enhancing splicing and removal of the introns whose retention is necessary for production of these isoforms. This could indirectly reduce splicing of the intron that is retained in the warm tim-M isoform and result in tim-M upregulation. Conversely, PSI could directly promote tim-M intron retention and indirectly downregulate production of tim-sc and tim-cold.
Other splicing factors have been shown to be involved in the control of circadian rhythms in Drosophila. SRm160 contributes to the amplitude of circadian rhythms by promoting per expression (Beckwith et al., 2017), while B52/SMp55 and PRMT5 regulate per’s most 3’ splicing, which is temperature sensitive (Zhang et al., 2018; Sanchez et al., 2010). Loss of PRMT5 results in essentially arrhythmic behavior (Sanchez et al., 2010), but this is unlikely to be explained by its effect on per’s thermosensitive splicing. B52/SMp55 knockdown flies show a reduced siesta, which is controlled by the same per splicing (Zhang et al., 2018). With the identification of Psi, we uncover a key regulator of tim alternative splicing pattern and show that this pattern determines circadian period length, while per alternative splicing regulates the timing and amplitude of the daytime siesta. Interestingly, a recent study identified PRP4 kinase and other members of tri-snRNP complexes as regulators of circadian rhythms (Shakhmantsir et al., 2018). Downregulation of prp4 caused excessive retention of the tim-M intron. PSI and PRP4 might thus have complementary functions in tim mRNA splicing regulation, working together to maintain the proper balance of tim isoform expression.
An unexpected finding is the role played by both PDF neurons and other circadian neurons in the short period phenotype observed with circadian locomotor rhythms when we knocked-down Psi. Indeed, it is quite clear from multiple studies that under constant darkness, the PDF-positive sLNvs dictate the pace of circadian behavior (Stoleru et al., 2005; Yao and Shafer, 2014). Why, in the case of Psi downregulation, do PDF negative neurons also play a role in period determination? The explanation might be that PSI alters the hierarchy between circadian neurons, promoting the role of PDF negative neurons. This could be achieved by weakening PDF/PDFR signaling, for example.
While we focused our work on PSI, several other interesting candidates were identified in our screen (Tables 1 and 2). We note the presence of a large number of splicing factors. This adds to the emerging notion that alternative splicing plays a critical role in the control of circadian rhythms. We have already mentioned above several per splicing regulators that can impact circadian behavior. In addition, a recent study demonstrated that specific classes of circadian neurons express specific alternative splicing variants, and that rhythmic alternative splicing is widespread in these neurons (Wang et al., 2018). Interestingly, in this study, the splicing regulator barc, which was identified in our screen and which has been shown to causes intron retention in specific mRNAs (Abramczuk et al., 2017), was found to be rhythmically expressed in LNds. Moreover, in mammals, alternative splicing appears to be very sensitive to temperature, and could explain how body temperature rhythms synchronize peripheral clocks (Preußner et al., 2017). Another intriguing candidate is cg42458, which was found to be enriched in circadian neurons (LNvs and Dorsal Neurons 1) (Wang et al., 2018). In addition to emphasizing the role of splicing, our screen suggests that regulation of polyA tail length is important for circadian rhythmicity, since we identified several members of the CCR4-NOT complex and deadenylation-dependent decapping enzymes. Future work will be required to determine whether these factors directly target mRNAs encoding for core clock components, or whether their effect on circadian period is indirect. Interestingly, the POP2 deadenylase, which is part of the CCR4-NOT complex, was recently shown to regulate tim mRNA levels post-transcriptionally (Grima et al., 2019). It should be noted that while our screen targeted 364 proteins binding or associated with RNA, it did not include all of them. For example, LSM12, which was recently shown to be a part of the ATXN2/TYF complex (Lee et al., 2017), was not included in our screen because it had not been annotated as a potential RAP when we initiated our screen.
In summary, our work provides an important resource for identifying RNA associated proteins regulating circadian rhythms in Drosophila. It identifies PSI is an important regulator of circadian period and circadian phase in response to thermal cycles, and points at additional candidates and processes that determine the periodicity of circadian rhythms.
Materials and methods
Fly stocks
Request a detailed protocolFlies were raised on a standard cornmeal/agar medium at 25°C under a 12 hr:12 hr light:dark (LD) cycle. The following Drosophila strains were used: w1118 -- w; tim-GAL4, UAS-dicer2/CyO (TD2) (Dubruille et al., 2009) -- y w; Pdf-GAL4, UAS-dicer2/CyO (PD2) (Dubruille et al., 2009) -- y w; Tim-GAL4/CyO (TG4) (Kaneko et al., 2000) -- y w; Pdf-GAL4 (PG4) (Renn et al., 1999) -- w;; UAS-dcr2 (Dietzl et al., 2007) -- y w;; timHA (Rutila et al., 1998) -- yw; TD2; Pdf-Gal80, Pdf-GAL80 (Zhang and Emery, 2013). The following combinations were generated for this study: y w; TG4; Pdf-GAL80, Pdf-GAL80 -- w; tim-GAL4/CyO; UAS-dicer2/TM6B -- tim0,TG4/CyO; UAS-Dcr2/TM6B -- tim0, PsiRNAiKK/CyO; timHA/TM6B. TD2, ptim-TIM-LUC and TD2, BG-LUC transgenic flies expressing a tim-luciferase and per-luciferase fusion gene respectively, combined with the TD2 driver, were used for luciferase experiments. The TIM-LUC fusion is under the control of the tim promoter (ca. 5 kb) and 1st intron (Lamba et al., 2018), BG-LUC contains per genomic DNA encoding the N-terminal two-thirds of PER and is under the control of the per promoter (Stanewsky et al., 1997). RNAi lines (names beginning with JF, GL, GLV, HM or HMS) were generated by the Transgenic RNAi Project at Harvard Medical School (Boston, MA) and obtained from the Bloomington Drosophila Stock Center (Indiana University, USA). RNAi lines (names beginning with GD or KK) and control lines (host strain for the KK library containing landing sites for the RNAi transgenes, VIE-260B, and tio misexpression control strain, 40D-UAS) were obtained from the Vienna Drosophila Stock Center. UAS-Psi flies were kindly provided by D. Rio (Labourier et al., 2001).
Behavioral monitoring and analysis
Request a detailed protocolThe locomotor activity of individual male flies (2–5 days old at start of experiment) was monitored in Trikinetics Activity Monitors (Waltham, MA). Flies were entrained to a 12:12 LD cycle for 3–4 days at 25°C (unless indicated) using I-36LL Percival incubators (Percival Scientific, Perry IA). After entrainment, flies were released into DD for five days. Rhythmicity and period length were analyzed using the FaasX software (courtesy of F. Rouyer, Centre National de la Recherche Scientifique, Gif-sur-Yvette, France) (Grima et al., 2002). Rhythmicity was defined by the criteria – power >20, width >1.5 using the χ2 periodogram analysis. Actograms were generated using a signal-processing toolbox implemented in MATLAB (MathWorks), (Levine et al., 2002). For phase-shifting experiments, groups of 16 flies per genotype were entrained to a 12:12 LD cycle for 5–6 days at 25°C exposed to a 5 min pulse of white fluorescent light (1500 lux) at different time points on the last night of the LD cycle. A separate control group of flies was not light-pulsed. Following the light pulse, flies were released in DD for six days. To determine the amplitude of photic phase shifts, data analysis was done in MS Excel using activity data from all flies, including those that were arrhythmic according to periodogram analysis. Activity was averaged within each group, plotted in Excel, and then fitted with a 4 hr moving average. A genotype-blind observer quantified the phase shifts. The peak of activity was found to be the most reliable phase marker for all genotypes. Phase shifts were calculated by subtracting the average peak phase of the light-pulsed group from the average peak phase of non-light pulsed group of flies. Temperature entrainment was performed essentially as described in Busza et al. (2007). Flies were entrained for 4–5 days in LD followed by 11 days in an 8 hr phase advanced temperature cycle. Behavior was analyzed between day 7 and day 10 of the temperature cycle. Actograms were used to ensure that all genotypes had reached – as expected from Busza et al. (2007) – a stable phase relationship with the temperature cycle. The phase of the evening peak of activity was determined as described for the phase response curve above. Because, under a LD cycle, the evening peak tend to be truncated by the light off transition, we used the approach described in Harrisingh et al. (2007), which compares the percent of activity between ZT17.5–23.5 that occurs between ZT20.5–23.5 (Morning anticipation phase score), or the percent of activity between ZT5.5–11.5 that occurs between ZT8.5–11.5 (Evening anticipation phase score). If phase is advanced, and activity increases earlier than normal, this percent will decrease.
Statistical analysis
Request a detailed protocolFor the statistical analysis of behavioral and luciferase period length, Student’s t-test was used to compare means between two groups, and one-way analysis of variance (ANOVA), coupled to post hoc tests, was used for multiple comparisons. Tukey’s post hoc test was used when comparing three or more genotypes and Dunnett’s post hoc test was used when comparing two experimental genotypes to one control. For the statistical analysis of qPCR and the behavioral phase-shifting experiments, two-way ANOVA, coupled to Tukey’s post hoc test, was used for multiple comparisons. Statistical analyses were performed using GraphPad Prism version 7.0 c for Mac OS X, GraphPad Software, La Jolla California USA, www.graphpad.com. P values and 95% Confidence Intervals are reported in data source files ‘Figure statistics’.
Luciferase experiments
Request a detailed protocolThe luciferase activity of whole male flies on Luciferin (Gold-biotech) containing agar/sucrose medium (170 μl volume, 1% agar, 2% sucrose, 25 mM luciferin), was monitored in Berthold LB960 plate reader (Berthold technologies, USE) in l-36LL Percival incubators with 90% humidity (Percival Scientific, Perry IA). Three flies per well were covered with needle-poked Pattern Adhesive PTFE Sealing Film (Analytical sales and services 961801). The distance between the agar and film was such that the flies were not able to move vertically. Period length was determined from light measurements taken during the first two days of DD. The analysis was limited to this window because TIM-LUC and BG-LUC oscillations severely dampened after the second day of DD. Period was estimated by an exponential dampened cosinor fit using the least squares method in MS Excel (Solver function).
Real-time quantitative PCR
Request a detailed protocolTotal RNA from about 30 or 60 fly heads collected at CT 3, CT9, CT15 and CT21 on the first day of DD were prepared using Trizol (Invitrogen) and Zymo Research Direct-zol RNA MiniPrep kit (R2050) following manufacturer’s instructions. 1 μg of total RNA was reverse transcribed using Bio-RAD iSCRIPT cDNA synthesis kit (1708891) following manufacturer’s instructions. Real-time PCR analysis was performed in triplicate (three technical replicates per sample) using Bio-RAD iTaq Universal SYBR Green Supermix (1725121) in a Bio-RAD C1000 Touch Thermal Cycler instrument. A standard curve was generated for each primer pair, using RNA extracted from wild-type fly heads, to verify amplification efficiency. Data were normalized to RpL32 (Dubruille et al., 2009) using the 2-ΔΔCt method. Primers used: RpL32-forward ATGCTAAGCTGTCGCACAAA; RpL32-reverse GTTCGATCCGTAACCGATGT; psi-forward GGTGCCTTGAATGGGTGAT; psi-reverse CGATTTATCCGGGTCCTCG; tim-M-forward TGGGAATCTCGCCCGAAAC; tim-M-reverse AGAAGGAGGAGAAGGAGAGAGG; tim-sc-forward ACTGTGCGATGACTGGTCTG; tim-sc-reverse TGCTTCAAGGAAATCTTCTG; tim-cold-forward CCTCCATGAAGTCCTCGTTCG; tim-cold-reverse ATTGAGCTGGGACACCAGG; cwo-foward TTCCGCTGTCCACCAACTC; cwo-reverse CGATTGCTTTGCTTTACCAGCTC; cwoRA-forward TCAAGTATGAGAGCGAAGCAGC; cwoRA-reverse TGTCTTATTACGTCTTCCGGTGG; cwoRB-forward GTATGAGAGCAAGATCCACTTTCC; cwoRB-reverse GATGATCTCCGTCTTCTCGATAC; cwoRC-forward GTATGAGAGCCAAGCGACCAC; cwoRC-reverse CCAAATCCATCTGTCTGCCTC.
Data availability
All source data are included in this submission.
-
NCBI Gene Expression OmnibusID GSE79916. The PSI-U1 snRNP interaction regulates male mating behavior in Drosophila.
References
-
Rhythmic behavior is controlled by the SRm160 splicing factor in Drosophila melanogasterGenetics 207:genetics.300139.2017.https://doi.org/10.1534/genetics.117.300139
-
Interactions between circadian neurons control temperature synchronization of Drosophila behaviorJournal of Neuroscience 27:10722–10733.https://doi.org/10.1523/JNEUROSCI.2479-07.2007
-
RBPDB: a database of RNA-binding specificitiesNucleic Acids Research 39:D301–D308.https://doi.org/10.1093/nar/gkq1069
-
Molecular genetic analysis of circadian timekeeping in DrosophilaAdvances in Genetics 74:141–173.https://doi.org/10.1016/B978-0-12-387690-4.00005-2
-
Intracellular Ca2+regulates free-running circadian clock oscillation in vivoJournal of Neuroscience 27:12489–12499.https://doi.org/10.1523/JNEUROSCI.3680-07.2007
-
Neural network interactions modulate CRY-Dependent photoresponses in DrosophilaThe Journal of Neuroscience 38:6161–6171.https://doi.org/10.1523/JNEUROSCI.2259-17.2018
-
Emerging roles for post-transcriptional regulation in circadian clocksNature Neuroscience 16:1544–1550.https://doi.org/10.1038/nn.3543
-
Photic signaling by cryptochrome in the Drosophila circadian systemMolecular and Cellular Biology 21:7287–7294.https://doi.org/10.1128/MCB.21.21.7287-7294.2001
-
Splicing of the period gene 3'-terminal intron is regulated by light, circadian clock factors, and phospholipase CMolecular and Cellular Biology 24:3359–3372.https://doi.org/10.1128/MCB.24.8.3359-3372.2004
-
Circadian rhythms and the circadian organization of living systemsCold Spring Harbor Symposia on Quantitative Biology 25:159–184.https://doi.org/10.1101/SQB.1960.025.01.015
-
The clockwork orange Drosophila protein functions as both an activator and a repressor of clock gene expressionJournal of Biological Rhythms 23:103–116.https://doi.org/10.1177/0748730407313817
-
The timSL mutant affects a restricted portion of the Drosophila Melanogaster circadian cycleJournal of Biological Rhythms 13:380–392.https://doi.org/10.1177/074873098129000200
-
The molecular ticks of the Drosophila circadian clockCurrent Opinion in Insect Science 7:51–57.https://doi.org/10.1016/j.cois.2015.01.002
-
FlyBase 2.0: the next generationNucleic Acids Research 47:D759–D765.https://doi.org/10.1093/nar/gky1003
-
A Drosophila RNAi library modulates hippo pathway-dependent tissue growthNature Communications 7:10368.https://doi.org/10.1038/ncomms10368
-
An isoform-specific mutant reveals a role of PDP1 epsilon in the circadian oscillatorJournal of Neuroscience 29:10920–10927.https://doi.org/10.1523/JNEUROSCI.2133-09.2009
Article and author information
Author details
Funding
National Institute of General Medical Sciences (1R35GM118087)
- Patrick Emery
National Institute of General Medical Sciences (1R01GM125859)
- Sebastian Kadener
The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.
Acknowledgements
We are particularly grateful to Vincent van der Vinne for his help with the analysis of the luciferase recording. We also express our gratitude to Monika Chitre for help with qPCR, Pallavi Lamba for help with PRC, Elaine Chang for help with luciferase recordings, and Diana Bilodeau-Wentworth, Dianne Szydlik, Chunyan Yuan and Vinh Phan for technical assistance. We also thank Dr. Donald Rio as well as the Bloomington and Vienna Drosophila Resource Centers for fly stocks. This work was supported by MIRA award 1R35GM118087 from the National Institute of General Medicine Sciences (NIGMS) to PE, and NIGMS grant 1R01GM125859 to SK.
Copyright
© 2019, Foley 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,612
- views
-
- 253
- downloads
-
- 25
- 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
-
- Neuroscience
Sound location coding has been extensively studied at the central nucleus of the mammalian inferior colliculus (CNIC), supporting a population code. However, this population code has not been extensively characterized on the single-trial level with simultaneous recordings or at other anatomical regions like the dorsal cortex of inferior colliculus (DCIC), which is relevant for learning-induced experience dependent plasticity. To address these knowledge gaps, here we made in two complementary ways large-scale recordings of DCIC populations from awake mice in response to sounds delivered from 13 different frontal horizontal locations (azimuths): volumetric two-photon calcium imaging with ~700 cells simultaneously recorded at a relatively low temporal resolution, and high-density single-unit extracellular recordings with ~20 cells simultaneously recorded at a high temporal resolution. Independent of the method, the recorded DCIC population responses revealed substantial trial-to-trial variation (neuronal noise) which was significantly correlated across pairs of neurons (noise correlations) in the passively listening condition. Nevertheless, decoding analysis supported that these noisy response patterns encode sound location on the single-trial basis, reaching errors that match the discrimination ability of mice. The detected noise correlations contributed to minimize the error of the DCIC population code of sound azimuth. Altogether these findings point out that DCIC can encode sound location in a similar format to what has been proposed for CNIC, opening exciting questions about how noise correlations could shape this code in the context of cortico-collicular input and experience-dependent plasticity.
-
- Neuroscience
Complexin determines magnitude and kinetics of synchronized secretion, but the underlying molecular mechanisms remained unclear. Here, we show that the hydrophobic face of the amphipathic helix at the C-terminus of Complexin II (CpxII, amino acids 115–134) binds to fusion-promoting SNARE proteins, prevents premature secretion, and allows vesicles to accumulate in a release-ready state in mouse chromaffin cells. Specifically, we demonstrate that an unrelated amphipathic helix functionally substitutes for the C-terminal domain (CTD) of CpxII and that amino acid substitutions on the hydrophobic side compromise the arrest of the pre-fusion intermediate. To facilitate synchronous vesicle fusion, the N-terminal domain (NTD) of CpxII (amino acids 1–27) specifically cooperates with synaptotagmin I (SytI), but not with synaptotagmin VII. Expression of CpxII rescues the slow release kinetics of the Ca2+-binding mutant Syt I R233Q, whereas the N-terminally truncated variant of CpxII further delays it. These results indicate that the CpxII NTD regulates mechanisms which are governed by the forward rate of Ca2+ binding to Syt I. Overall, our results shed new light on key molecular properties of CpxII that hinder premature exocytosis and accelerate synchronous exocytosis.