- Research article
- Open Access
cryptochrome genes form an oscillatory loop independent of the per/tim loop in the circadian clockwork of the cricket Gryllus bimaculatus
Zoological Letters volume 3, Article number: 5 (2017)
Animals exhibit circadian rhythms with a period of approximately 24 h in various physiological functions, including locomotor activity. This rhythm is controlled by an endogenous oscillatory mechanism, or circadian clock, which consists of cyclically expressed clock genes and their product proteins. cryptochrome (cry) genes are thought to be involved in the clock mechanism, and their functions have been examined extensively in holometabolous insects, but in hemimetabolous insects their role is less well understood.
In the present study, the role of cry genes was investigated using RNAi technology in a hemimetabolous insect, the cricket Gryllus bimaculatus. Using a molecular cloning approach, we obtained cDNAs for two cry genes: Drosophila-type cry1 (Gb’cry1) and mammalian-type cry2 (Gb’cry2). Gb’cry2 has six splicing variants, most of which showed rhythmic mRNA expression. Gb’cry1 RNAi treatment had only a limited effect at the behavioral and molecular levels, while Gb’cry2 RNAi had a significant effect on behavioral rhythms and molecular oscillatory machinery, alone or in combination with Gb’cry1 RNAi. In Gb’cry1/Gb’cry2 double-RNAi crickets, most clock genes showed arrhythmic expression, except for timeless, which retained clear rhythmic expression. Molecular analysis revealed that some combination of Gb’cry1 and Gb’cry2 variants suppressed CLK/CYC transcriptional activity in cultured cells.
Based on these results, we propose a new model of the cricket’s circadian clock, including a molecular oscillatory loop for Gb’cry2, which can operate independent of the Gb’per/Gb’tim loop.
Circadian rhythms are a biological rhythm with a period of approximately 24 h, which can be observed in various physiological functions of animals, including general activity, hormonal secretion, and metabolism. Rhythmicity is generated by an endogenous mechanism, the so-called circadian clock, which consists of cyclic expression of clock genes and their product proteins. In insects, the core clock components are period (per), timeless (tim), Clock (Clk), and cycle (cyc) [1, 2]. Clk and cyc encode the transcription factors CLOCK (CLK) and CYCLE (CYC), which form a heterodimer and activate transcription of per and tim. Transcripts of per and tim are subsequently translated to the product proteins PERIOD (PER) and TIMELESS (TIM), which heterodimerize, enter the nucleus, and inhibit transcriptional activity of CLK/CYC. This negative feedback loop produces the 24 h oscillation. The two cryptochromes, cry1 (or Drosophila-type cry) and cry2 (or mammalian-type cry), are also important components involved in the insect clock. cry proteins are members of the photolyase family, and cry1 functions as a photoreceptor, resetting the clock through TIM degradation upon exposure to blue light in Drosophila and some other insects [3, 4]. cry2 is thought to be a component of the clock, working together with PER to form a feedback machinery . This hypothesis has been supported by a reporter assay using cultured cells [3, 6]. However, to date the roles of cry genes have been investigated mostly in holometabolous insects, meaning their functions in hemimetabolous insects must be investigated if we are to understand their role in these insects and their functional diversification.
The cricket, Gryllus bimaculatus, is well suited to the study of these questions, since the molecular machinery of the clock has been extensively studied and RNAi can be used to dissect the roles of clock genes [7, 8]. The clock is located in the optic lobe, and its oscillatory components include Gb’per, Gb’tim, Gb’Clk, and Gb’cyc, as in Drosophila, but their roles differ substantially from those in Drosophila . Similar to the case in Drosophila, Gb’per and Gb’tim are rhythmically expressed and their transcripts increase at night, peaking from early night to midnight [9, 10]. However, unlike in Drosophila, Gb’cyc, not Gb’Clk, is rhythmically expressed [11, 12]. Double-stranded RNA (dsRNA) treatment of Gb’per and Gb’Clk abolished the locomotor rhythm and the molecular expression rhythms of Gb’per and Gb’tim [9, 11, 13]. However, Gb’tim and Gb’cyc RNAi arrests the rhythmic expression of Gb’per and Gb’tim, but does not prevent behavioral rhythmicity [10, 12], suggesting the existence of other oscillatory components. We suggest that cry is a candidate component as, in some insects, CRY functions as a transcriptional repressor by itself, at least in cultured cells .
In the present study, we obtained two cry genes in G. bimaculatus by molecular cloning and confirmed that they are homologous to cry1 (Gb’cry1) and cry2 (Gb’cry2) by structural analysis. PCR analysis revealed that there are six splicing variants of Gb’cry2. We then examined their expression profiles and functional roles in the clock machinery by RNAi experiments. RNAi of Gb’cry2 substantially changed the free-running period of the locomotor rhythm, and the change was further enhanced by Gb’cry1 RNAi, suggesting that the cry genes may be involved in the determination of the free-running period. Although the molecular oscillation of clock genes was mostly halted by Gb’cry1/Gb’cry2 double RNAi, Gb’tim retained its rhythmic expression. Cellular reporter assays revealed that some combination of Gb’cry1 and Gb’cry2 variants represses CLK/CYC transcriptional activity. We propose a unique model for the molecular clock of the cricket that incorporates these new findings and previously described properties of clock genes.
Third instar nymphs and adult male crickets, Gryllus bimaculatus, were used. These were purchased or obtained from a laboratory colony maintained under standard environmental conditions with a lighting regimen of LD 12:12 (light: 0600-1800; Japan Standard Time) and a constant temperature of 25 ± 1.0 °C. They were fed laboratory chow and water.
Cloning and structural analysis of the cry genes
Total RNA was extracted with ISOGEN (Nippon Gene, Tokyo, Japan) from 20 heads of third instar nymphs collected at ZT 10 (ZT stands for zeitgeber time and ZT0 corresponds to lights-on and ZT12 to lights-off). We used 4.5 μg of total RNA for reverse transcription to obtain cDNA, using SuperScript II (Invitrogen, Carlsbad, CA, USA). Using the single-stranded cDNA as a template, we performed PCR with degenerate primers deduced from the conserved amino acid sequences among insect Drosophila-type cryptochrome (cry1) and mammalian-type cryptochrome (cry2) homologues. The primers used for cry1 and cry2 were 5′-CAGGAACAAAGATGGTGGGATAYAAYMGNATG-3′ for forward and 5′-CCAGTTTCCTGCGCACACNSWCCARTC-3′ for reverse, and 5′-TGTTCGTGATCCGAGGACARCCNGCNGA-3′ for forward and 5′-CACTGGGATGTACTTTCGGATGWARTCNCCRTT-3′ for reverse, respectively. The PCR conditions employed were 30 s for denaturation at 95 °C, 30 s for annealing at 57 °C, and 1 min 30 s for extension at 72 °C for 35 cycles with ExTaq DNA polymerase (Takara, Otsu, Japan). The purified fragment was cloned into TOPO-pCR II vector (Invitrogen) and sequenced with BigDye Terminator v3.1 Cycle Sequencing Kit (Applied Biosystems, Foster City, CA, USA). 5′ and 3′ RACEs for cry1 and cry2 were done with GeneRacer kit (Invitrogen) and SMARTer RACE cDNA Amplification kit (Takara) with gene specific primers, 5′-TGGTTCACAATCCTGCTCAA-3′ and 5′-ATGGAGCGCAATCCTATCTG-3′, and 5′-CCACTTGGCTAAGGCTTCTG-3′ and 5′-GCAGTTCCATCAAGTCAGCA-3′, respectively. RACE fragments were purified, cloned and sequenced as mentioned above. Sequences were analyzed by Genetyx ver. 6 (Genetic Information Processing Software, Tokyo, Japan) and BioEdit ver. 7.0.5 (Biological Sequence Alignment Editor, Ibis Therapeutic, Carlsbad, CA, USA). There were six splicing variants in the cricket cry2 cDNA. For cloning these cDNAs, primer sets of 5′-ATGCAGGTACCATGGGAAGCACACTTGCATT-3′ for forward and 5′-ATGCAGCGGCCGCAATTGTGCTGTTGATTTAAAC-3′ for reverse, or 5′-ATGCAGGTACCAGTGCTCGTGTGTGTGTTTG-3′ for forward and 5′-ATGCAGCGGCCGCAATTGTGCTGTTGATTTAAAC-3′ for reverse were used. Exon/intron structure of cry genes were analyzed with genomic DNA sequence data.
Amino acid sequences of CRY1 and CRY2 were aligned using ClustalW in MEGA 7.0. A phylogenetic tree of CRYs was constructed using maximum likelihood methods based on JTT matrix-based model in MEGA 7.0. Sequences of known insect cry1 and cry2 genes were obtained from GenBank.
Measurement of mRNA levels
Quantitative real-time reverse transcription polymerase chain reaction (qPCR) and reverse transcription polymerase chain reaction (RT-PCR) were used to measure mRNA levels. Total RNA was extracted and purified from six adult male optic lobes with TRIzol Reagent (Invitrogen). To remove contaminating genomic DNA, the total RNA was treated with DNase I. About 500 ng total RNA of each sample was reverse transcribed with random hexamers using PrimeScript RT reagent Kit (Takara). Real-time PCR was performed by Mx3000P Real-Time PCR System (Stratagene, La Jolla, CA, USA) using FastStart Universal SYBR Green Master (Roche, Tokyo, Japan) including SYBR Green with primers designed for Gb’cry1, Gb’cry2, Gb’per (BAG48878), Gb’tim (BAJ16356), and Gb’rpl18a (DC448653) (Table 1). In all cases, a single expected amplicon was confirmed by melting analysis. Quantification was performed based on a standard curve obtained from a known amount of templates. The results were analyzed using the software associated with the instrument.
To estimate expression levels of Gb’cry2 variants, semi-quantitative RT-PCR was performed with Gb’cry2-x2a, Gb’cry2-x4 and Gb’cry2-x10b primer sets that were designed to detect variants with or without Gb’cry2 exon2a, 4 and 10b, respectively (Table 1). Gb’rpl18a was also used as internal reference. The PCR conditions were 10 s for denaturation at 98 °C, 30 s for annealing at 52 °C, and 30 s for extension at 72 °C for 35 cycles with Emeraldamp® MAX PCR Master Mix (Takara, Otsu, Japan) for Gb’cry2 exon2a, 10 s at 98 °C, 30 s at 55 °C, and 30 s at 68 °C for 33 cycles with KOD FX Neo (Toyobo, Osaka, Japan) for exon4 and 10b, 10 s at 98 °C, 30 s at 55 °C, and 30 s at 72 °C for 28 cycles with Emeraldamp® MAX PCR Master Mix for Gb’rpl18a. The PCR samples were electrophoresed on 1.5% agarose gels in TBE buffer (89 mM Tris-base pH 7.6, 89 mM boric acid, 2 mM EDTA). Gels were stained with GelRedTM nucleic acid gel stains (Biotium, CA, USA) and photographed on a 280 nm UV light box (TP-20MP, ATTO, Tokyo, Japan). The gel images were digital-imaged and quantified by imageJ software (available from https://imagej.nih.gov/ij/). The quantity of the PCR samples were estimated relative to the internal reference gene, Gb’rpl18a, amplified with the same PCR protocol.
For both qPCR and RT-PCR, the values were normalized with those of Gb’rpl18a at each time point. Results of 3–5 independent experiments were used to calculate the mean ± SEM.
dsRNA for Gb’cry1, Gb’cry2, Gb’per (GenBank/EMBL/DDBJ Accession No. BAG48878), Gb’tim (BAJ16356), Gb’Clk (AB738083), Gb’cyc (AB762416) and DsRed2 derived from a coral species (Discosoma sp.), were synthesized using MEGAscript High Yield Transcription kit (Ambion, Austin, TX, USA). For Gb’cry1 and Gb’cry2, cDNAs prepared as described above, were used as templates of PCR, which was performed with ExTaq DNA polymerase (Takara). The T7 or T3 containing primers used were listed in Table 1. Amplified Gb’cry1#d1 (575 bp), Gb’cry1#d2 (515 bp), Gb’cry2#d1 (463 bp), Gb’cry2#d2 (423 bp), Gb’per (456 bp), Gb’tim (519 bp), Gb’Clk (407 bp), and Gb’cyc (450 bp) fragments were extracted with phenol/chloroform and precipitated with ethanol then resuspended in Ultra Pure Water (Invitrogen). For DsRed2 dsRNA, linearized DsRed2 fragment was amplified from pDsRed2-N1 (Clontech, Mountain View, CA, USA), with the primers shown in Table 1. With each of these linearized fragments as a template, RNA was synthesized with T7 or T3 RNA polymerase. Synthesized RNAs were extracted with phenol/chloroform, and suspended in 50 μl TE solution after isopropanol precipitation. The yield and quality of RNA were assessed by absorbance using a spectral photometer (Genequant Pro, Amersham Bioscience, Piscataway, NJ, USA) and the same amount of sense and antisense RNA were mixed. The RNAs were denatured for 5 min at 100 °C and annealed by a gradual cool down to room temperature. After ethanol precipitation, the obtained dsRNA was suspended in Ultra Pure Water (Invitrogen) and adjusted to the final concentration of 20 μM. The dsRNA solution was stored at –80 °C until use. 760 nl of dsRNA solution was injected with a nanoliter injector (WPI, Sarasota, FL, USA) into the abdomen of adults anesthetized with CO2.
Assays in cultured cells
The coding sequences of cry genes, Gb’cry1 and Gb’cry2a–f, were subcloned individually into pAc5.1B-V5/His (Invitrogen, Carlsbad, CA). Primers used for cloning are listed in Table 1. Cultured Drosophila S2 cells were plated in 24-well tissue culture plates in Shields and Sang M3 insect medium (Sigma, St Louis, MO) supplemented with 12.5% fetal bovine serum (Biowest, Canada) and antibiotics (12.5 U/ml penicillin, 12.5 mg/ml streptomycin; GIBCO, Grand Island, NY). S2 cells were transfected by Effectene Transfection Reagent (Qiagen, Hilden, Germany) with 100 ng of pAct-Clk as an activator , 100 ng of Drosophila tim-luc , which contains a typical E-box  showing a relatively high expression , and 10 ng of pAc5.1-Rluc  as a positive control for luciferase activity, along with a combination of 500 ng of pAc5.1-Gb’cry1 and pAc5.1-Gb’cry2a ~ f. As needed, the empty vector pAc5.1B-V5/His was used instead of the vector constructs of Gb’cry1 and Gb’cry2 genes to ensure that an equal amount of DNA was used for transfection in each well. Forty-eight hours after transfection, luciferase assays of transfected cells were carried out using the Dual-Luciferase® Reporter Assay System (Promega Corporation, Madison, Wisconsin, USA) and then normalized to Rluc activity as a control for transfection efficiency. Assays were performed at least three times.
Locomotor activities were recorded in the same way as described previously . Briefly, the final instar nymphs or adult crickets were individually housed in a transparent plastic box (18 × 9 × 4.5 cm) with a rocking substratum. The number of substratum rocking was recorded every 6 min by a computerized system. Food and water were provided ad libitum. The actographs were placed in an incubator (MIR-153, Sanyo Biomedica, Osaka, Japan) in which temperature was kept at 25 ± 0.5 °C and lighting conditions were maintained by a cool white fluorescent lamp connected to an electric timer. The light intensity was 600–1000 lux at the animal’s level, varying with the proximity to the lamp. The raw data were displayed as conventional double-plotted actograms to judge activity patterns, and the free-running period was calculated by the chi-square periodogram  with Actogram J (freely available at http://actogramj.neurofly.de/) . If a peak in the periodogram appeared above the 0.05 confidence level (alpha = 0.005), the power value (height of the peak above the confidence level) was greater than or equal to 10, and the width of the peak was greater than or equal to 2, the period for the peak was designated statistically significant .
The one-way analysis of variance (ANOVA) followed by a post hoc Tukey-test was used to compare the differences in means of mRNA levels between the different time points and the differences in means of free-running periods and phase angle differences between insect groups with different treatments. To compare the means of two groups, t-test was used. Dunnett’s-test was also used for multiple comparisons, where applicable. Significance of the daily and circadian rhythmicity in mRNA expression was examined by the single cosinor method , fitting a cosine curve of 24 h period using Time Series Analysis Serial Cosinor 6.3 (Expert Soft Technologie, Richelieu, France). In all statistical tests, the significance level was set at P < 0.05.
Cloning and structural analysis of Gb’cry1 and Gb’cry2
cry1 and cry2 homologues were cloned from the optic lobe of G. bimaculatus using a degenerate PCR strategy with degenerate primers based on conserved amino acid sequences from insects cry1 and cry2 homologues. A Gryllus bimaculatus cry1 (Gb’cry1) fragment of 1116 bps and a cry2 (Gb’cry2) fragment of 876 bps, including the FAD-binding 7 domain and a part of PhrB, was first obtained. By 3′ rapid amplification of cDNA ends (RACE) and 5′ RACE, we obtained a full length of 1937 bp Gb’cry1 cDNA (GenBank/EMBL/DDBJ Accession No. LC202047). The Gb’cry1 cDNA had 189 bp 5′-untranslated regions (UTR) and 131 bp 3′-UTR. Genomic DNA analysis revealed that Gb’cry1 consists of nine exons, and no splicing variants were detected by RT-PCR. The putative product protein Gb’CRY1 consists of 537 amino acid residues (Fig. 1a).
As to Gb’cry2, we obtained six variants of cDNAs, those are Gb’cry2a ~ f (GenBank/EMBL/DDBJ Accession No. LC202048, LC202049, LC202050, LC202051, LC202052, LC202053) (Fig. 1a). The lengths of cDNA varied among the variants. Their structures are illustrated in Fig. 1a with deduced proteins. With genomic sequence analysis, the Gb’cry2 gene was shown to consist of 11 exons and the splicing variants were shown to commonly include exons 1–3 and 5–11 with variable lengths of exons 2, 6 and 10. Exceptionally, a variant Gb’cry2e lacked exon 4. Gb’cry2a ~ c had a short exon 2 lacking the 2a region, Gb’cry2c with short exon 6 lacking 6a region, and Gb’cry2a, d with longer exon 10 including a 10b region. From these mRNAs six different proteins are deduced. Gb’CRY2a ~ c lack the photolyase α/β domain, because the start codon is located in middle of exon 4, Gb’CRY2c has a shorter N-terminal region lacking exon 6a region, and Gb’CRY2a and -d have a shorter C-terminal region because of stop codon in exon 10b. Gb’CRY2e is the smallest lacking most part and consists of only 202 amino acid residues, while Gb’CRY2f is the largest, includes all coding regions, and consists of 587 amino acid residues.
A BLAST database search indicated that the amino acid sequence of Gb’CRY1 has 55–57% identities along the entire length of the protein with those of CRY1 characterized in other insects, including the fruit fly Drosophila melanogaster (GenBank/EMBL/DDBJ accession No. AAF55649.1), the monarch butterfly, Danaus plexippus (EHJ63675.1), the silkmoth, Antheraea pernyi (AAK11644.1), and the mosquito, Anopheles gambiae (ABB29886.1, Table 2). The longest Gb’CRY2f has 62–79% identity along the entire length of the protein with those of known insect CRY2s and mouse CRYs, including the monarch butterfly, D. plexippus (EHJ74426.1), the silkmoth, A. pernyi (ABO38435.1), the honeybee, Apis mellifera (NP_001077099.1), the mosquito, A. gambiae (ABB29887.1), the red flour beetle, Tribolium castaneum (EFA04537.1) and the mouse, Mus musculus (CRY1, NP_031797.1, CRY2, NP_034093.1) (Table 2).
The deduced proteins Gb’CRY1 and Gb’CRY2f commonly had five highly conserved regions that are characteristic of CRY proteins (Fig. 1a): (i) DNA photolyase α/β domain, (ii) FAD-binding 7 domain, (iii) RD-2b domain, which is necessary for nuclear localization and the repression of CLOCK: BMAL-mediated transcription, (iv) inhibition of CLOCK-ARNTL-mediated transcription (ICAT) domain, and (v) PhrB domain. Identity of these domains had 51–69% in Gb’CRY1, while those in Gb’CRY2f had high identity of 52–100% and also share high identity with mouse CRY1 and CRY2 (Table 2). A phylogenetic tree based on the amino acid sequences of CRYs from known insects and some vertebrates revealed that Gb’CRY1 and Gb’CRY2f form separate clades with CRY1 and CRY2 in other insects (Fig. 1b).
Expression profiles of Gb’cry1 and Gb’cry2 mRNA in the optic lobe
To examine whether Gb’cry1 and Gb’cry2 transcript oscillated in the optic lobe, their mRNA levels under LD12:12 were measured using qPCR in the adult male crickets. Primers used for qPCR were shown in Fig. 1a and Table 1. As shown in Fig. 2a, Gb’cry1 appeared to be constitutively expressed, and no clear daily rhythm was observed (ANOVA, F5,23 = 1.03, P > 0.42; cosinor, P > 0.05). For Gb’cry2, rhythmic mRNA expression was found with a peak during the midnight (ZT18, ANOVA, F5,21 = 21.95, P < 0.01; cosinor, P < 0.01) when all the splicing variants were amplified with primers designed at exon 11 (Fig. 2b). The amplitude of the rhythm was approximately 4-fold.
We then examined the expression levels of Gb’cry2a ~ c and Gb’cry2d ~ f with short and long exon 2 using a primer set Gb’cry2-x2a (Table 1), respectively, in the optic lobe at ZT 6 and 18 by a semi-quantitative method, RT-PCR. Both of Gb’cry2a ~ c and Gb’cry2d ~ f mRNAs with short and long exon 2, respectively, showed a higher level at ZT18, but the expression level was much greater in the latter (Fig. 2c; ANOVA, F3,11 = 454.26, P < 0.01). When the levels of Gb’cry2a,d and Gb’cry2b, c, e, f with long and short exon 10, respectively, were examined using a primer set Gb’cry2-x10b (Table 1), both were again higher at ZT18 but the levels of Gb’cry2b, c, e, f were higher (Fig. 2c; ANOVA, F3,11 = 110.3, P < 0.01). These results suggest that the Gb’cry2f was expressed most abundantly. We also compared the levels of cry2e that lacked exon 4 and Gb’cry2a ~ d,f that included exon 4 using a primer set Gb’cry2-x4 (Table 1). The results showed that they again showed higher expression at ZT18 (ANOVA, F3,11 = 167.36, P < 0.01) and that most of expressed Gb’cry2 variants were those including exon 4 and the cry2e expression was only negligible (Fig. 2c).
RNAi suppresses Gb’cry1 and Gb’cry2 expression
We examined the effects of dsRNA of Gb’cry1 (dsGb’cry1) and Gb’cry2 (dsGb’cry2) on respective transcript levels, by measuring their mRNA levels in the optic lobe of the adult male crickets by qPCR. dsGb’cry2#d1 and dsGb’cry2#d2 were synthesized for the region included in all splicing variants (Fig. 1a). Treatment of dsGb’cry1#d1 or dsGb’cry1#d2 significantly reduced Gb’cry1 mRNA levels throughout the day (Fig. 2a). In the crickets injected with dsGb’cry2#d1 or dsGb’cry2#d2, the levels of Gb’cry2 mRNA were knocked down to nearly the basal level of that of intact crickets; although a small but significant fluctuation was observed (ANOVA: F5,19 = 5.44, P < 0.05 for dsGb’cry2#d1; F5,23 = 3.18, P < 0.05 for dsGb’cry2#d2, Fig. 2b), no significant daily rhythm was detected on single cosinor analysis (P > 0.05).
Effects of Gb’cry1 and Gb’cry2 RNAi on the locomotor rhythm
To examine the role of Gb’cry1 and Gb’cry2 in circadian locomotor rhythm regulation, locomotor activity was recorded in adult males injected with dsGb’cry1 (dsGb’cry1#d1, N = 19; dsGb’cry1#d2, N = 22) or dsGb’cry2 (dsGb’cry2#d1, N = 37; dsGb’cry2#d2, N = 22). Since similar results were obtained in crickets treated with two different dsRNAs for Gb’cry1 and Gb’cry2, the results were pooled. We used DsRed2 RNAi crickets (n = 21) as a negative control. As was the case for DsRed2 RNAi crickets, all of the Gb’cry1 RNAi and the Gb’cry2 RNAi crickets exhibited a nocturnal activity rhythm under LD12:12, with a major peak at lights-off and a minor peak at lights-on (Fig. 3). In the ensuing constant darkness (DD), the rhythm free-ran, except in 4 Gb’cry2 RNAi crickets which became arrhythmic, and the free-running periods varied with treatments (Figs. 3 and 4). The Gb’cry1 RNAi (n = 16) and DsRed2 RNAi (n = 15) crickets showed free-running periods shorter than 24 h, averaging 23.58 ± 0.25 h and 23.71 ± 0.24 h, respectively. However, the free-running period varied widely in Gb’cry2 RNAi crickets (n = 33): the average period was 24.12 ± 0.34 h ranging from 23.2 h to 25.0 h. When treated doubly with dsGb’cry1 and dsGb’cry2, the crickets were all rhythmic but showed a wider range of free-running period (22.4–25.8 h) than Gb’cry2 RNAi crickets. Notably, some of the Gb’cry2 RNAi crickets (dsGb’cry2#d1, n = 3; dsGb’cry2#d2, n = 2) showed a rhythm that split into two components in the free-running condition (as exemplified in Fig. 3d) and that individuals with the free-running periods longer than 24 h often became arrhythmic (dsGb’cry2#d1, N = 4) or only weakly rhythmic (dsGb’cry2#d1, n = 8; dsGb’cry2#d2, n = 10) in due course (as exemplified in Fig. 3d, j).
Effects of Gb’cry1 and Gb’cry2 RNAi on the molecular clockwork
We then examined the effects of RNAi of Gb’cry1 and Gb’cry2 on Gb’per and Gb’tim genes, using dsGb’cry1#d1 and dsGb’cry2#d1. In intact crickets, transcripts of Gb’per and Gb’tim showed a rhythmic expression to peak around midnight (ZT18). These expression profiles were similar to those reported previously [9–11]. When Gb’cry1 was knocked-down by RNAi, Gb’cry2 and Gb’tim maintained a rhythmic expression (ANOVA, F5,24 = 15.10, P < 0.01; cosinor, P < 0.05 for Gb’cry2 and ANOVA, F5,23 = 21.38, P < 0.01; cosinor, P < 0.05 for Gb’tim) with a peak similar to that of untreated crickets (Gb’cry2) or advanced by about 4–8 h to peak during late day to early night (Gb’tim); Gb’per also showed a similar rhythmic pattern to Gb’tim, showing a significant fluctuation (ANOVA, F5,21 = 3.59, P < 0.05), but no daily rhythm was detected by the single cosinor method (P > 0.05) (Fig. 5a). The results of Gb’cry2 RNAi are shown in Fig. 5b. No clear effect was observed on Gb’cry1 mRNA levels. Gb’tim showed a rhythm (ANOVA, F5,20 = 3.90, P < 0.05; cosinor, P < 0.05) with a peak at late day (ZT10), while Gb’per fluctuated around a medium range of the control, with no significant daily rhythm (ANOVA, F5,21 = 1.54, P > 0.23; cosinor, P > 0.05).
When Gb’cry1 and Gb’cry2 were simultaneously knocked down, both Gb’cry1 and Gb’cry2 showed no significant daily rhythm (Fig. 5c: ANOVA, F5,22 = 0.54, P > 0.73; cosinor, P > 0.05 for Gb’cry1 and ANOVA, F5,22 = 2.14 P > 0.11; cosinor, P > 0.05 for Gb’cry2). Gb’per showed an arrhythmic expression pattern with mRNA levels around a medium range of control (ANOVA, F5,23 = 2.02, P > 0.12; cosinor, P > 0.05), while Gb’tim maintained a robust oscillation with a peak slightly advanced to peak at ZT14 (Fig. 5c: ANOVA, F5,20 = 9.36, P < 0.01; cosinor, P < 0.05).
Gb’cry1 and Gb'cry2 regulates CLK/CYC transcriptional activity
We examined the role of Gb'cry1 and Gb’cry2 in the molecular clockwork using cell-based assay. Various combinations of Gb’cry1 and Gb’cry2a ~ f were transfected into S2 cells together with Drosophila Clk and Drosophila tim-Luc. When Gb’cry2c was co-transfected with Gb’cry1 or Gb’cry2f, transcriptional activity of CLK/CYC was significantly suppressed (Dunnett’s test, P < 0.05; Fig. 6). CLK/CYC transcriptional activity was enhanced when Gb’cry2e was co-transfected with Gb’cry2a or Gb’cry2d (Dunnett’s test, P < 0.05; Fig. 6). Transfection in other combination had no significant effects on transcriptional activity through the E-box.
Effects of RNAi of Gb’per, Gb’tim, Gb’Clk and Gb’cyc on Gb’cry1 and Gb’cry2 expression
To investigate the relationship of Gb’cry1 and Gb’cry2 with a Gb’per/Gb’tim loop, we examined the effects of RNAi of Gb’per, Gb’tim, Gb’Clk and Gb’cyc on the expression levels of Gb’cry1 and Gb’cry2. These treatments have shown that molecular oscillation of respective gene is eliminated [9–12]. The RNAi of these clock genes had no clear effects on Gb’cry1 mRNA levels except for Gb’cyc RNAi at ZT22, at which a significant increase in Gb’cry1 mRNA was observed (Fig. 7a). The same treatments, however, significantly affected Gb’cry2 mRNA expression (Fig. 7b). Gb’per RNAi upregulated the Gb’cry2 mRNA levels to raise the bottom level to exceed the peak level of control crickets and lead to arrhythmic expression (ANOVA, F5,18 = 1.12, P > 0.39; cosinor, P > 0.05). In Gb’Clk RNAi crickets, Gb’cry2 mRNA showed no expression rhythm, remaining at the median range (ANOVA, F5,17 = 2.01, P > 0.14; cosinor, P > 0.05). Gb’cry2 mRNA retained a daily rhythm with a normal level in Gb’tim RNAi and Gb’cyc RNAi crickets (ANOVA, F5,23 = 4.62, P < 0.01; cosinor, P < 0.05 for Gb’tim RNAi and ANOVA, F5,17 = 7.02, P < 0.01; cosinor P < 0.05 for Gb’cyc RNAi), but the amplitude was slightly attenuated with the increased levels during daytime. These results strongly suggest that rhythmic expression of Gb’cry2 is regulated by a mechanism other than the Gb’per/Gb’tim loop.
Gb’cry1 and Gb’cry2 genes of the cricket
In insects, two cry genes have been identified: one is the Drosophila-type cry1 and the other is the mammalian-type cry2 [3, 22]. The present study showed that the cricket, Gryllus bimaculatus also has two cry genes. Gb’cry2 has six splicing variants with variations in exons 2, 4, 6 and 10. The expression levels were different among those variants and greatest in Gb’cry2f with an exon 4, long exons 2 and 6, and a short exon 10. Their daily expression profiles in the optic lobe are similar to those in other insects; Gb’cry2 showed a rhythmic expression with a peak around mid-night, while Gb’cry1 is rather constitutively expressed throughout a day with a slight reduction at the late day phase. The reduction of Gb’cry1 may be caused by a light-induced suppression mechanism that is known to regulate cry1 in other insects . Gb’cry1 expression is most likely independent of the circadian clock, since RNAi of clock genes had nearly no effect on Gb’cry1 mRNA levels (Fig. 7a).
Role of Gb’cry1 and Gb’cry2 in the cricket’s clock
In flies and butterflies, CRY1 is supposed to be a blue light receptor, and resets the clock by leading to a degradation of TIMELESS protein in a light-dependent manner [4, 6, 23]. Our results, however, clearly show that Gb’cry1 RNAi crickets synchronized to the given light dark cycle as robustly as the control crickets, suggesting that Gb’CRY1 is not the major photoreceptor for entrainment of the clock. This finding is consistent with our previous reports that the compound eye is the only photoreceptor necessary for light entrainment of the cricket clock  and that opsin-long wavelength (opsin-LW) expressed in the compound eye plays a major role in photic entrainment . Similarly Gb’cry2 is apparently not a major photoreceptor for entrainment, because Gb’cry2 RNAi crickets also robustly synchronized to light cycles, similar to the case for Gb’cry1 RNAi crickets. However, we could not exclude the possibility that the two cry genes play subtle roles in photic entrainment. To draw a conclusion on the involvement of these genes in photic entrainment, careful examination of entrainability in Gb’cry1 RNAi and Gb’cry2 RNAi crickets is needed.
On the free-running rhythm in DD, Gb’cry2 RNAi had significant effects; the treated crickets showed variable free-running periods, suggesting that Gb’cry2 plays a significant role in the clock oscillatory mechanism. In contrast, Gb’cry1 RNAi had no clear effect, and the treated crickets showed a free-running rhythm with a period similar to that of control crickets treated with dsDsRed2. However, Gb’cry1 RNAi enhanced the variability of the free-running period of Gb’cry2 RNAi crickets (Figs. 3, 4), suggesting that Gb’cry1 and Gb’cry2 cooperate to determine the free-running period.
Gb’CRYs are transcriptional repressors
The present study showed that Gb’CRYs play in the cricket molecular clockwork as transcriptional repressors (Fig. 6). We have shown for the first time that Gb’cry2 has six splicing variants and that none has a repressor activity alone. The repressor activity was evident only when Gb’cry2c was co-expressed with Gb’cry1 or Gb’cry2f (Fig. 6). The requirement of two types of CRYs contrasts to mammalian and other insect clocks where a single molecular species of CRY can work as transcriptional repressor [3, 26, 27]. Since Gb’cry2c, Gb’cry2f, and Gb’cry1 commonly have ICAT and RD-2b, the former is likely important for repression of CLK/CYC transcriptional activity  and the latter may be required for the nuclear translocation . Gb’cry1’s role as a transcriptional repressor in combination with Gb’cry2 (Fig. 6) is reminiscent of a finding from peripheral tissues of Drosophila, in which cry1 is also known to function as a transcriptional repressor [23, 30]. For example, overexpression of cry1 and per in the Drosophila compound eye represses CLK/CYC-activated transcription . However, the Gb’cry1’s role in transcriptional regulation is contrast to an earlier report on Anopheles gambiae, Antheraea pernyi, and Danaus plexippus, in which cry1 has no repressor activity . The role of cry1 in the core clock mechanism may be lost in these insects, in which cry2 plays a major role as a transcriptional repressor. However, we also suggest that the role of cry1 should be carefully reexamined in those insects, since Yuan et al.  did not examine the role of cry1 in combination with cry2.
The role of Gb’CRYs as a component in the clock machinery partly explains the changes of free-running periods in Gb’cry2 RNAi crickets. Even after RNAi treatment it is expected that a small amount of mRNA will be expressed. The variable amount of expressed Gb’cry2 probably resulted in variable free-running periods among individuals treated by Gb’cry2 RNAi or Gb’cry1/Gb’cry2 double RNAi (Fig. 4) as suggested for CRYs in the mammalian suprachiasmatic nucleus (SCN) clock . In Gb’cry1 RNAi crickets, the lack of Gb’cry1 may be compensated by another complex formed by Gb’CRY2c and Gb’CRY2f, which would account for the lack of any observable significant changes in this period.
Our reporter assay also showed that CLK/CYC transcriptional activity was enhanced when Gb’cry2e was co-expressed with Gb’cry2a or Gb’cry2d. Although this possibility cannot be ruled out, their role as transcriptional activator may be negligible, since Gb’cry2e was expressed at only trace levels and the expression of Gb’cry2a and Gb’cry2d was also quite low (Fig. 2c).
Another possible role of Gb’cry2 may be as a coupling factor between clock neurons. This hypothesis is based on the observation that Gb’cry2 RNAi crickets sometimes show a rhythm dissociation into some components running with different free-running periods in DD (Fig. 3d). The change in coupling strength by Gb’cry2 RNAi may affect in variable degrees the free-running period of the locomotor rhythm. This may represent a parallel case with that described in reports on mammalian clocks in the SCN [32, 33]. Some of the Gb’CRY2 variants might play a role in this coupling and/or period change in individual clock neurons caused by Gb’cry2 RNAi may indirectly affect the coupling. This issue remains to be addressed in future studies.
Gb’cry2 oscillation is regulated by a mechanism different from the Gb’per/Gb’tim loop
We have shown that the cricket circadian clock includes a loop for rhythmic expression of Gb’per and Gb’tim [9, 10], of which transcription is regulated by Gb’CLK/Gb’CYC [11, 12]. It has been suggested that CRY2 forms a complex with PER or the PER/TIM complex to enter the nucleus and repress the CLK/CYC transcriptional activity in other insects [5, 6]. However, the results of the present study suggest that Gb’CRY2 composes a loop to repress Gb’CLK/Gb’CYC activity by forming complex between its variants or with Gb’CRY1, and that the loop can operate independent of the Gb’per/Gb’tim loop to some degree (Fig. 8). This hypothesis is based on the following observations. When Gb’tim was knocked down by RNAi or Gb’tim and Gb’per stopped their oscillation by Gb’cyc RNAi , Gb’cry2 maintains its oscillation within a normal range but with a slightly attenuated amplitude (Fig. 7). Alternatively when Gb’cry2 is knocked down, Gb’tim maintained its oscillation. Our data suggest that the circadian locomotor rhythm is expressed when either Gb’cry2 or Gb’tim is rhythmically expressed.
However, there must be an inter-connection between the assumed Gb’cry2 loop and the Gb’per/Gb’tim loop. When Gb’per or Gb’Clk was knocked down by RNAi, Gb’cry2 lost its rhythmic expression (Fig. 7), and in Gb’per RNAi crickets the mRNA levels of Gb’cry2 were significantly increased, with the trough level nearly equivalent to the peak level of DsRed2 RNAi crickets (Fig. 7). Alternatively when Gb’cry2 was knocked down, Gb’per lost the rhythm but Gb’tim showed a phase advance of its expression rhythm (Fig. 5). Although further detailed study is necessary to understand the underlying mechanism, these data suggest that some coupling mechanism exists between the two loops. For Gb’per, Gb’tim and Gb’cry2, Gb’CLK/Gb’CYC may be a common transcriptional activator. Thus, the two loops may couple with Gb’CLK/Gb’CYC as a hinge (Fig. 8). The transcriptional activity of Gb’CLK/Gb’CYC may be suppressed by both Gb’PER/Gb’TIM complex and Gb’CRY1/Gb’CRY2 or Gb’CRY2 complexes. Knock-down of either Gb’cry1 or Gb’cry2 weakens the repression of CLK/CYC, resulting in the upregulation of Gb’per mRNA. The transcription of Gb’tim and Gb’cry2 may require stronger transactivation than Gb’per, and thus Gb’tim and Gb’cry2 mRNA maintained rhythmic expression even after knockdown of Gb’cry1/Gb’cry2 and Gb’tim, respectively.
The cricket clock consists of at least four loops, each of which regulates rhythmic expression of Gb’per and Gb’tim, Gb’Clk, Gb’cyc, and Gb’cry2 [9, 10, 12] (Fig. 8). The clock system appears to be more complex than those of Drosophila and other insects, which basically include three loops for rhythmic expression of per and tim, Clk or cyc, and clockwork orange [1, 34]. The cricket is a hemimetabolous insect, and therefore more primitive than holometabolous insects, including Diptera (Drosophila), Lepidoptera (Danaus, Antheraea), and Hymenoptera (Apis). It may thus possess a more primitive clock machinery, of which some components may have changed or been lost in the course of phylogenetic development in higher insects.
The cricket’s clock possesses some properties common to mammalian clock. Gb’cyc shows a rhythmic expression similar to that of its mammalian homologue, Bmal1 . Both Gb’cry1 and Gb’cry2 are involved in the clock machinery. The reset mechanism is exclusively relies on retinal photoreceptors [24, 25, 35]. Similarly the firebrat Thermobia domestica, one of the most primitive ametabolous insects, has a mammalian-like clock machinery in which cyc is rhythmically expressed [36, 37]. These findings suggest that the ancestral insect clock resembled the mammalian clock and subsequently diversified into various types .
The evolution and diversification of insect cry genes remains a challenging issue. There is a variation for presence and absence of cry genes. Although many insects have both cry1 and cry2, some dipteran species lack cry2 and some coleopteran and hymenopteran species lack cry1 [1, 38]. The present study revealed that Gb’CRY1 is not a photoreceptor, but rather a component of the clock and that Gb’cry2 has splicing variants, providing another example of diversification of cry1 and cry2. These findings make it difficult to resolve the evolution or diversification of insect cry genes. To deal with this issue, extensive comparative studies are desired across phylogenetically diverse species, including hemimetabola and ametabola.
Analysis of variance
- Clk :
- cry :
- cry1 :
- cry2 :
- cyc :
Double stranded RNA
- Gb :
Inhibition of CLOCK-ARNTL-mediated transcription domain
- opsin-LW :
- per :
Quantitative real-time reverse transcription polymerase chain reaction
Rapid amplification of cDNA ends
- rpl18a :
ribosomal protein L18a
Reverse transcription polymerase chain reaction
- tim :
Tomioka K, Matsumoto A. Circadian molecular clockworks in non-model insects. Curr Opin Insect Sci. 2015;7:58–64.
Tataroglu O, Emery P. The molecular ticks of the Drosophila circadian clock. Curr Opin Insect Sci. 2015;7:51–7.
Yuan Q, Metterville D, Briscoe AD, Reppert SM. Insect cryptochromes: gene duplication and loss define diverse ways to construct insect circadian clocks. Mol Biol Evol. 2007;24:948–55.
Ceriani MF, Darlington TK, Staknis D, Mas P, Petti AA, Weitz CJ, et al. Light-dependent sequentation of TIMELESS by CRYPTOCHROME. Science. 1999;285:553–6.
Rubin EB, Shemesh Y, Cohen M, Elgavish S, Robertson HM, Bloch G. Molecular and phylogenetic analyses reveal mammalian-like clockwork in the honey bee (Apis mellifera) and shed new light on the molecular evolution of the circadian clock. Genome Res. 2006;16:1352–65.
Zhu H, Sauman I, Yuan Q, Casselman A, Emery-Le M, Emery P, et al. Cryptochromes define a novel circadian clock mechanism in monarch butterflies that may underlie sun compass navigation. PLoS Biol. 2008;6:e4.
Tomioka K. Chronobiology of crickets: a review. Zool Sci. 2014;31:624–32.
Tomioka K, Sakamoto T, Moriyama Y. RNA interference is a powerful tool for chronobiological study in the cricket. Sleep Biol Rhythms. 2009;7:144–51.
Moriyama Y, Sakamoto T, Karpova SG, Matsumoto A, Noji S, Tomioka K. RNA interference of the clock gene period disrupts circadian rhythms in the cricket Gryllus bimaculatus. J Biol Rhythms. 2008;23:308–18.
Danbara Y, Sakamoto T, Uryu O, Tomioka K. RNA interference of timeless gene does not disrupt circadian locomotor rhythms in the cricket Gryllus bimaculatus. J Insect Physiol. 2010;56:1738–45.
Moriyama Y, Kamae Y, Uryu O, Tomioka K. Gb'Clock is expressed in the optic lobe and required for the circadian clock in the cricket Gryllus bimaculatus. J Biol Rhythms. 2012;27(6):467–77.
Uryu O, Karpova SG, Tomioka K. The clock gene cycle plays an important role in the circadian clock of the cricket Gryllus bimaculatus. J Insect Physiol. 2013;59:697–704.
Moriyama Y, Sakamoto T, Matsumoto A, Noji S, Tomioka K. Functional analysis of the circadian clock gene period by RNA interference in nymphal crickets Gryllus bimaculatus. J Insect Physiol. 2009;55:396–400.
Darlington TK, Wager-Smith K, Ceriani MF, Staknis D, Gekakis N, Steeves TDL, et al. Closing the circadian loop: CLOCK-induced transcription of its own inhibitors per and tim. Science. 1998;280:1599–603.
McDonald MJ, Rosbash M, Emery P. Wild-type circadian rhythmicity is dependent on closely spaced E boxes in the Drosophila timeless promotor. Mol Cell Biol. 2001;21:1207–17.
Ueda HR, Matsumoto A, Kawamura M, Tanimura T, Hashimoto S. Genome-wide transcriptional orchestration of circadian rhythm in Drosophila. J Biol Chem. 2002;277:14048–52.
Matsumoto A, Ukai-Tadenuma M, Yamada RG, Houl J, Uno KD, Kasukawa T, et al. A functional genomics strategy reveals clockwork orange as a transcriptional regulator in the Drosophila circadian clock. Gene Dev. 2007;21:1687–700.
Sokolove PG, Bushell WN. The chi square periodogram: its utility for analysis of circadian rhythm. J Theor Biol. 1978;72:131–60.
Schmid B, Helfrich-Förster C, Yoshii T. A new ImageJ plug-in “ActogramJ” for chronobiological analyses. J Biol Rhythms. 2011;26:464–7.
Kaneko M, Hall JC. Neuroanatomy of cells expressing clock genes in Drosophila: transgenic manipulation of the period and timeless genes to mark the perikarya of circadian pacemaker neurons and their projections. J Comp Neurol. 2000;422:66–94.
Nelson W, Tong Y, Lee J, Halberg F. Methods for cosinor-rhythmometry. Chronobiologia. 1979;6:305–23.
Zhu H, Yuan Q, Briscoe AD, Froy O, Casselman A, Reppert SM. The two CRYs of the butterfly. Curr Biol. 2005;15:R953–954.
Ivanchenko M, Stanewsky R, Giebultowicz JM. Circadian photoreception in Drosophila: functions of cryptochrome in peripheral and central clocks. J Biol Rhythms. 2001;16:205–15.
Tomioka K, Chiba Y. Effects of nymphal stage optic nerve severance or optic lobe removal on the circadian locomotor rhythm of the cricket, Gryllus bimaculatus. Zool Sci. 1984;1:385–94.
Komada S, Kamae Y, Koyanagi M, Tatewaki K, Hassaneen E, Saifullah A, et al. Green-sensitive opsin is the photoreceptor for photic entrainment of an insect circadian clock. Zool Lett. 2015;1:11.
Ye R, Selby CP, Chiou Y-Y, Ozkan-Dagliyan I, Gaddameedhi S, Sancar A. Dual modes of CLOCK:BMAL1 inhibition mediated by Cryptochrome and Period proteins in the mammalian circadian clock. Gene Dev. 2016;28:1989–98.
Langmesser S, Tallone T, Bordon A, Rusconi S, Albrecht U. Interaction of circadian clock proteins PER2 and CRY with BMAL1 and CLOCK. BMC Mol Biol. 2008;9:41.
Ingram KK, Kutowoi A, Wurm Y, Shoemaker D, Meier R. The molecular clockwork of the fire ant Solenopsis invicta. PLoS One. 2012;7:e45715.
Hirayama J, Nakamura H, Ishikawa T, Kobayashi Y, Todo T. Functional and structural analyses of cryptochrome. Vertebrate CRY regions responsible for interaction with the CLOCK:BMAL1 heterodimer and its nuclear localization. J Biol Chem. 2003;278:35620–8.
Collins B, Massoni EO, Stanewsky R, Blau J. Drosophila CRYPTOCHROME is a circadian transcriptional repressor. Curr Biol. 2006;16:441–9.
Anand SN, Maywood ES, Chesham JE, Joynson G, Banks GT, Hastings MH, et al. Distinct and separable roles for endogenous CRY1 and CRY2 within the circadian molecular clockwork of the suprachiasmatic nucleus, as revealed by the Fbxl3Afh mutation. J Neurosci. 2013;33:7145–53.
Ono D, Honma S, Honma K-I. Cryptochromes are critical for the development of coherent circadian rhythms in the mouse suprachiasmatic nucleus. Nat Commun. 2013;4:1666.
Evans JA, Pan H, Liu AC, Welsh DK. Cry1 −/− circadian rhythmicity depends on SCN intercellular coupling. J Biol Rhythms. 2012;27:443–52.
Hardin PE. Essential and expendable features of the circadian timekeeping mechanism. Curr Opin Neurobiol. 2006;16:686–92.
Tomioka K, Okada Y, Chiba Y. Distribution of circadian photoreceptors in the compound eye of the cricket Gryllus bimaculatus. J Biol Rhythms. 1990;5:131–9.
Kamae Y, Tanaka F, Tomioka K. Molecular cloning and functional analysis of the clock genes, Clock and cycle, in the firebrat Thermobia domestica. J Insect Physiol. 2010;56:1291–9.
Kamae Y, Uryu O, Miki T, Tomioka K. The nuclear receptor genes HR3 and E75 are required for the circadian rhythm in a primitive insect. PLoS One. 2014;9:e114899.
Zhan S, Merlin C, Boore JL, Reppert SM. The monarch butterfly genome yields insights into long-distance migration. Cell. 2011;147:1171–85.
The authors thank Dr. Sumihare Noji of Tokushima University for providing genomic information on the cricket.
This study was supported in part by grant-in-aid No. 15H0440017 from JSPS.
Availability of data and materials
The datasets supporting the conclusions of this article are included within the article.
AT, TQI, SH, OU, YD and MN performed the experiments. AT and KT wrote the paper. KT, AT and YD designed the experiments. TB and TT supported molecular analysis of cry genes. All authors read and approved the final manuscript.
The authors declare they have no competing interests.
Ethics approval and consent to participate
Consent for publication
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Shinryo Hori Deceased
About this article
Cite this article
Tokuoka, A., Itoh, T.Q., Hori, S. et al. cryptochrome genes form an oscillatory loop independent of the per/tim loop in the circadian clockwork of the cricket Gryllus bimaculatus . Zoological Lett 3, 5 (2017) doi:10.1186/s40851-017-0066-7
- Circadian clock
- Clock gene
- Molecular oscillatory mechanism