The Resonance and Adaptation of Neurospora crassa Circadian and Conidiation Rhyth ms to Short Light-Dark Cycles

Circadian clocks control the physiological and behavioral rhythms to adapt to the environment with a period of ~24 h. However, the influences and mechanisms of the extreme light/dark cycles on the circadian clock remain unclear. We showed that, in Neurospora crassa, both the growth and the microconidia production contribute to adaptation in LD12:12 (12 h light/12 h dark, periodically). Mathematical modeling and experiments demonstrate that in short LD cycles, the expression of the core clock protein FREQUENCY was entrained to the LD cycles when LD > 3:3 while it free ran when T ≤ LD3:3. The conidial rhythmicity can resonate with a series of different LD conditions. Moreover, we demonstrate that the existence of unknown blue light photoreceptor(s) and the circadian clock might promote the conidiation rhythms that resonate with the environment. The ubiquitin E3 ligase FWD-1 and the previously described CRY-dependent oscillator system were implicated in regulating conidiation under short LD conditions. These findings shed new light on the resonance of Neurospora circadian clock and conidiation rhythms to short LD cycles, which may benefit the understandings of both the basic regulatory aspects of circadian clock and the adaptation of physiological rhythms to the extreme conditions.


Introduction
The Earth rotates with a period of approximately 24 h, which causes periodic changes in many environmental factors. Circadian clocks are the inner mechanisms that allow organisms to adjust their physiology and behavior to the daily cycling of environmental factors, e.g., light and temperature [1].
In the natural daily cycling conditions, rhythms with the same 24 h period to the environmental factors occur in organisms, which are called diurnal rhythms. Whereas in constant conditions in which any environmental cues are absent, e.g., light, dark, temperature or humidity, organisms display rhythms with periods slightly deviated from 24 h, which are called free-running rhythms [2]. The period length of regular alternation of temporal cues, e.g., light and dark, which is different from 24 h, is defined as T-cycle [3]. In certain ranges, non-24 h light-dark (LD) cycles can entrain the circadian rhythms, and the flexibility of this entrainability dramatically varies between different organisms. Within the range of entrainment, which is also called the limit of stable entrainment, the rhythm changes in accordance with the non-24 h T-cycles and show stable phase angles. Phase angle represents the relationship between the timing of the rhythm and the timing of an points the cultures were harvested by filtration, frozen in liquid nitrogen, and ground in liquid nitrogen.
An incubator (Percival Scientific, USA) was used to simulate different LD cycles, the white light intensity during light was 5000 lux or 1000 lux as indicated. The red light (λ = 660 nm) was generated by an LED lamp, and the intensity was 5000 lux or 1000 lux as indicated. The blue light (λ = 470 nm) was generated by an LED lamp, and the intensity was 5000 lux or 1000 lux as indicated.
The counting of microconidia was conducted on a hemocytometer under an optical microscope. Briefly, since the shape of microconidia is nearly round and they usually contain a single nucleus, we consider those hyphae with length/width ratio <1.5 as a subjective criteria as microconidia. The proportion of microconidia was the ratio of the microconidia quantity divided by the total quantity of both microconidia and macroconidia.

Race Tube Assay
The rhythmicity of Neurospora conidiation is observed in race tubes, which are hollow glass growth tubes, bent up at both ends in order to contain the solid agar medium sufficient for linear growth for~one week. In a race tube assay, the conidia are inoculated on the surface of the solid media at one end of the race tube thus it grows toward the other end. During growth, Neurospora yields conidiation bands periodically, which permit the measurement of the period length by calculating the growth time between two adjacent conidiation bands [21].
In this work, the conidiation banding profiles were assayed on race tubes under standard conditions [14]. The growth front was marked every 24 h under a red safe light in all light conditions in this work. All race tube experiments were carried out at room temperature (25 • C).

Protein Analysis
Protein extraction, quantification and western blot analysis were conducted as described elsewhere [41,42]. For Western blotting, equal amounts of total protein (50 µg) were loaded in each protein lane of SDS-PAGE (7.5%) gels containing a ratio of 37.5:1 acrylamide/bisacrylamide.

Dynamic Modeling
The modeling was conducted based on the following assumptions: (1) frq mRNA is synthesized in the nucleus and transfers to the cytosol where it is degraded; (2) nuclear FRQ protein, and combined the processes of multiple phosphorylation and translocation of the nuclear to a time delay; (3) the rate of synthesis of FRQ is proportional to frq mRNA several h ago; (4) wc-1 mRNA and two states of nuclear WCC complex are involved in this model: phosphorylated and light dependent; (5) to consider the adaptation effect to light, VVD is introduced to the model ( Figure S1; Table S2 and supplemental protocols for details in Supplemental File S1).

RNA-Seq and Analysis
The strains subjected to RNA-seq were grown and the RNA samples were isolated from duplicates. All of the Pearson's correlation values between each duplicate were >0.99. Tophat was used as the aligner to map the reads to the reference genome [N. crassa OR74A (NC12)]. The genes showing ≥1.5-fold increase were selected for further analysis. See detailed RNA-seq methods and parameters in Supplemental File S1.

Statistical Analysis
Data are mean ± SE or mean ± SD where indicated. n ≥ 3. The Student's t test was used for all statistical analyses. * Represents the p-value of the statistical tests is less than the significance level of 0.05 (p ≤ 0.05); ** represents p ≤ 0.01 and *** represents p ≤ 0.001.
The arrhythmicity of Neurospora under extremely short LD cycles (e.g., LD1:1) could be explained by two possibilities: (1) the conidiation rhythms were abolished; (2) the bands were too compact and ambiguous to detect. To validate, we used the non-band strain FGSC4200 which grows very fast in race tube. In constant dark, FGSC4200 shows mild conidiation rhythms [45], while in LD1:1 it showed weak but recognizable conidiation rhythms ( Figure S2). These data suggest that although the free-running circadian period of Neurospora is approximately 22 h [20], the conidiation rhythm can be driven by extremely short LD cycles.
The circadian clock can adjust to fit with the cycling environment, especially the light. We compared the growth rate of two frq mutants frq 2 and frq 7 , which possess shorter and longer free-running periods in constant dark, respectively [46]. The race tube assay results reveal that in LD9:9 the growth ratio of frq 2 vs. frq 7 was higher than that in LD13.5:13.5 ( Figure 1B,C; Table S3), suggesting a resonance of Neurospora circadian period to the environment.
However, when we compared the growth rate of the wild-type strain under various LD cycles, we obtained unexpected results showing that Neurospora grew much faster in constant dark than that in constant light or most LD conditions. Moreover, the growth rate in shorter LD cycles was not necessarily slower than LD12:12. For instance, the growth in constant darkness (DD), LD6:6, LD3:3 and LD1:1 was significantly faster than that in LD12:12. The growth under LD65 min:25 min was also significantly faster than LD12:12 ( Figure 1D). The dramatic differences in growth under LL vs. DD, LD45 min:45 min vs. LD65 min:25 min indicate that the growth is influenced by the daily length of lighting in Neurospora.
Neurospora microconidia are small uninucleate spores that serve as male gametes or as asexual reproductive structures which are important to stress adaptation [47]. When we counted the microconidia produced under LD12:12 and short LD cycles, respectively, and the results show that the proportion of microconidia in LD12:12 was significantly higher than those in other LD conditions. Statistic results revealed a significant correlation between the LD cycles and microconidia proportion while there was no significant correlation between LD cycles and growth rates ( Figures 1E and S3A,B).

Expression of the FRQ Protein in Different LD Cycles
As the circadian clock is involved in controlling the conidiation rhythmicity, we asked whether the molecular rhythmicity of the FRQ protein is subject to changes in response to the length of the LD cycles. As one of the core factors in Neurospora circadian clock, FRQ acts as the negative component in the transcription/translational feedback loop [20]. We conducted mathematical modeling and the results suggested that under LD12:12 and LD6:6, the concentration of FRQ protein matches the LD cycles ( Figure 2A; Table S2; protocol and parameters for modeling in File S1), instead, under LD2:2 the FRQ levels show a free running pattern with some saw-tooth shaped fluctuations. The simulation results show that the mRNA levels have a direct response to light, but the protein curves are smoother, especially for high-frequency light signals. Therefore, the post-translational regulatory and translocation processes may perform a buffering function or high-frequency filtering function. Furthermore, the numerical data also predicted that the direct response of frq mRNA to the light stimuli is lower under shorter LD cycles (Figure 2A), which may be due to the light adaptation module (VVD-WCC loop) and the core negative feedback loop. The arrhythmicity of Neurospora under extremely short LD cycles (e.g., LD1:1) could be explained by two possibilities: (1) the conidiation rhythms were abolished; (2) the bands were too compact and ambiguous to detect. To validate, we used the non-band strain FGSC4200 which grows very fast in race tube. In constant dark, FGSC4200 shows Representative results are shown (n ≥ 3), (C) ratio of the growth rates between frq 2 and frq 7 calculated according to the results of (B). Data are mean ± SD, n = 3, (D) relative growth rate of WT strain calculated according to the results of (B). The daily growth lengths were normalized to the data in LD12:12. Data are mean ± SE, n = 3, (E) proportion of microconidia produced under different LD cycles. Data are mean ± SE, n = 3. The light intensity was 5000 lux. ** represents p ≤ 0.01 and *** represents p ≤ 0.001. We next conducted Western blotting and showed that under LD12:12 and LD6:6, the expression of the FRQ protein and the phosphorylation profile displayed coincident rhythms with the periods of LD cycles. However, under LD3:3 and shorter cycles, the FRQ levels seemed to free run but we also observed some saw-tooth shaped fluctuations under LD3:3 ( Figure 2B-G). These data support the mathematical modeling.

FWO Is Not Required for Adaptation of Conidiation Rhythms to Short LD Cycles in High-Intensity Light
In the Neurospora circadian oscillator, WC-1 and WC-2 are the two positive elements in which WC-1 functions as a blue light sensor [23]. We asked whether WCC proteins are responsible for conidiation rhythmicity under LD cycles, therefore, we grew the wc-1 RIP , wc-2 KO and wcc DKO and frq 10 strains in race tubes under different LD cycles and observed condition rhythms for all of these mutants exhibit conidiation rhythms under LD12:12, LD6:6, LD3:3andLD2:2, but not in DD and LD45 min:45 min (Figure 3). Among these strains, wc-1 RIP is a wc-1-disrupted strain using the RIPing method [23]; wc-2 KO and frq 10 are deletion mutants lacking functional WC-2 and FRQ, respectively; wcc DKO is a double knockout strain of wc-1 and wc-2. The strains were exposed to a light intensity of 5000 lux, during the light periods. These data suggest that under short LD conditions, the conidiation rhythms cycle to the environment even without a functional circadian clock. We have also shown that in WT, the FRQ protein oscillation is masked under LD cycles (6 ≤ T< 24), but free runs when T < 6 ( Figure 2). Taken together, these results suggest that the FWO system is not required for the conidiation banding rhythms driven by the cycling environment.
with a period of approximately 6 h. If the frequency of the light signal is further increased to LD2:2, although the frq mRNA can respond to the induction by the light, the FRQ levels show the free running pattern with some fluctuations, (B-G) Western blot results showing FRQ levels in different LD cycles. The LD routines are LD12:12 for two cycles (B), LD3:3 for 48 h (C), LD2:2 for 48 h (D), LD6:6 for 48 h (E), LD3:3 for two cycles (F) and LD2:2 for two cycles (G), respectively. The black and white bars denote different LD regimes. Membranes (Mem) stained with Amido Black were used as loading control. The densitometric analysis of FRQ levels is shown at the bottom of each panel. Representative data of triplicates are shown. The light intensity was 5000 lux. Data are mean ± SE, n = 3.
We next conducted Western blotting and showed that under LD12:12 and LD6:6, the expression of the FRQ protein and the phosphorylation profile displayed coincident rhythms with the periods of LD cycles. However, under LD3:3 and shorter cycles, the FRQ levels seemed to free run but we also observed some saw-tooth shaped fluctuations under LD3:3 ( Figure 2B-G). These data support the mathematical modeling.

FWO is not Required for Adaptation of Conidiation Rhythms to Short LD Cycles in High-Intensity Light
In the Neurospora circadian oscillator, WC-1 and WC-2 are the two positive elements in which WC-1 functions as a blue light sensor [23]. We asked whether WCC proteins are responsible for conidiation rhythmicity under LD cycles, therefore, we grew the wc-1 RIP , wc-2 KO and wcc DKO and frq 10 strains in race tubes under different LD cycles and observed condition rhythms for all of these mutants exhibit conidiation rhythms under LD12:12, LD6:6, LD3:3andLD2:2, but not in DD and LD45 min:45 min (Figure 3). Among these strains, wc-1 RIP is a wc-1-disrupted strain using the RIPing method [23]; wc-2 KO and frq 10 are deletion mutants lacking functional WC-2 and FRQ, respectively; wcc DKO is a double knockout strain of wc-1 and wc-2. The strains were exposed to a light intensity of 5000 lux, during the light periods. These data suggest that under short LD conditions, the conidiation rhythms cycle to the environment even without a functional circadian clock. We have also shown that in WT, the FRQ protein oscillation is masked under LD cycles (6 ≤ T< 24), but free runs when T < 6 ( Figure 2). Taken together, these results suggest that the FWO system is not required for the conidiation banding rhythms driven by the cycling environment.

Conidiation Rhythms of Photosensor-Related Mutants under Short LD Regimes
First, we investigated the impacts of individual photosensing-associated genes on the conidiation rhythms under short LD cycles in ∆cry; ∆vvd; ∆phy-1; ∆phy-2; ∆nop-1 and in cog-1, a newly identified mutant showing circadian conidiation rhythmicity under constant light [14]. WC-1, VVD and CRY are blue light sensors while PHY-1 and PHY-2 are potential red-light receptors [24,27,35]. Unexpectedly, all of these strains showed conidiation rhythms with periods in accordance with the LD cycling periods, under conditions of LD12:12, LD6:6, LD3:3 and LD2:2 ( Figure 4A-D). Together with the fact that the wc-1 RIP strain can be entrained at LD2:2, these results suggest that either there exists additional photosensor(s), or the lack of an individual photosensor is not sufficient to abolish the adaptation of conidiation to short LD cycling cues.
The temperature increase that accompanied the switch-on of the lighting in the incubator, which might also have affected the conidiation. To rule out this possibility, we measured conidiation rhythms under temperature cycles (25.5 • C 6 h: 24.5 • C 6 h) as the range of incubator temperature is~±0.5 • C in our experiments, as previously described [14]. The results showed no resonance in all tested strains ( Figure S5H), demonstrating this temperature variation is not sufficient for entrainment. Taken together, these results show that all of these tested light sensors and associated factors were not critical for the conidiation rhythms. Instead, an unidentified blue light photoreceptor might be involved in controlling the conidiation rhythms.
The effects of different photosensing-associated genes on the production of carotenoid in all of these strains in constant light were also analyzed. As has been previously shown, the strains containing vvd mutations, including ∆vvd; ∆vvd, cog-1; ∆cry, ∆vvd and ∆cry, ∆vvd, cog-1 displayed bright orange color in their mycelia [14,32]. In contrast, the mycelia color of strains ∆vvd, ∆wc-1, cog-1; ∆cry, ∆vvd, ∆wc-1 and ∆vvd, ∆cry, ∆wc-1, cog-1 was not bright orange ( Figure S6A). These data suggest that some of these light sensors might regulate the functions of others.
In contrast to our results, Nsa et al. showed that the ∆wc-1 strain exhibit no obvious rhythms under a number of T-cycle conditions including LD6:6, LD9:9, LD12:12 and LD14:14 [14]. However, we used a white light intensity (~5000 lux) that was much higher than that used by Nsa et al. (1200 lux), which may account for the inconsistency. In support of this view, when we grew the clock mutants in LD conditions and the light intensity was approximately 1000 lux, and conidiation rhythmicity was abolished at the tested LD conditions ( Figures 5A,B, S5 and S6B). Conversely, the conidiation rhythms in clock gene mutants, wc-1 RIP , wcc DKO , frq 10 and ∆wc-2, were present in 5000 lux blue light but were absent in 1000 lux ( Figure S5). These data suggest that while clock components might act to promote the function of a potential photoreceptor to elicit conidiation rhythms, the response of clock mutants to a high-intensity but not low-intensity light, suggests that higher-intensity light might overcome the repression effect caused by the loss of clock components ( Figure 5C-E).
It is intriguing to note that in BD3:3 (blue light intensity: 1000 lux), the strains cog-1; ∆wc-1, cog-1 and ∆vvd, ∆wc-1, cog-1 exhibited biomodal patterns on the first day ( Figure S5D), reflecting the role of the CDO pathway in mediating the resonance of conidiation rhythmicity although the underlying mechanisms remain unclear.

Implication of FWD-1 in Regulating the Conidiation Rhythms under Short LD Cycles
To probe the possible mechanisms regulating the adaptation to LD cycles, we also checked the conidiation rhythms in the fwd-1 RIP strain that lacks the functional fwd-1 gene, as FWD-1 has been implicated in Neurospora circadian clock through mediating the FRQ ubiquitination and turn over. Neurospora FWD-1 is the homolog of Drosophila Slimb, a ubiquitin E3 ligase [39,49].
The fwd-1 RIP strain exhibited no overt conidiation bands under LD3:3 and shorter LD cycles (Figures 4, 5 and S5). The mutants lacking tested photosensors showed normal responses to short LD cycles, suggesting that the phenotype of fwd-1 RIP is not based on the impacts of FWD-1 on these genes, i.e., additional photosensor(s) controlled by FWD-1 might be implicated in regulating the conidiation rhythms under short LD cycles. Alternatively, a potential factor connecting the specific photosensor to conidiation and, that is controlled by FWD-1, might be implicated ( Figure 5C,D). gene mutants, wc-1 RIP , wcc DKO , frq 10 and wc-2, were present in 5000 lux blue light but were absent in 1000 lux ( Figure S5). These data suggest that while clock components might act to promote the function of a potential photoreceptor to elicit conidiation rhythms, the response of clock mutants to a high-intensity but not low-intensity light, suggests that higher-intensity light might overcome the repression effect caused by the loss of clock components ( Figure 5C-E). It is intriguing to note that in BD3:3 (blue light intensity: 1000 lux), the strains cog-1; wc-1, cog-1 and vvd, wc-1, cog-1 exhibited biomodal patterns on the first day (Figure The white light intensity was 1000 lux. Triangles denote the conidiation bands. (C-E). Schematics of the control of conidiation rhythms under short LD cycles of WT in low-intensity light (C) and clock mutants in low-intensity light (D) and high-intensity light (E). WC-1, VVD and additional photoreceptor(s)are implicated in the regulation of the conidiation rhythms. In 6 h≤ T ≤ 24 h, the FWO system is entrained, and the FRQ rhythms are not endogenous as they cannot be maintained in constant dark. The question masks denote unidentified photoreceptor and putative coupler linking FWO and conidiation, respectively. Transducers denote the factors linking photoreceptor and conidiation.

Transcriptomic Analysis of Genes Driven by High-Intensity Light
The above data suggest that high-intensity light influences conidiation rhythms, thus we conducted RNA-seq analysis to assess the changes of gene expression in response to high-intensity light. We compared the transcriptomic changes of the strain ∆cry, ∆vvd, ∆wc-1, which contains no known blue light photoreceptors, after exposure for 45 min to either 1000 lux or 5000 lux white light ( Figure 6A). The results showed that upon light exposure, a variety of genes involved in metabolism and gene expression were up-regulated (folds ≥ 1.5). For the genes exclusively driven by 5000-lux light, one putative motif (GAXGA) was identified by the XXmotif online tool (http://xxmotif.genzentrum.lmu.de/, accessed on 15 January 2018), which is present in 56 promoter areas of the 73 genes ( Figure 6B,C). Light of 5000-lux driven the exclusive expression of 73 specific genes, many of which were enriched in fatty acid metabolism, ribosome biogenesis, etc. (Figure 6D; supplemental File S2 and S3). The induced expression of three representative genes (NCU00298, NCU00992 and NCU09771) were confirmed in ∆cry, ∆vvd, ∆wc-1; ∆phy-1, ∆phy-2 and ∆nop-1 by qRT-PCR ( Figure 6E-G). Although the pathways of ribosome biogenesis and propanoate metabolism were induced in both 1000-lux and 5000-lux light (File S3), the genes involved were different, suggesting that higher-intensity light plays specific roles in regulating the expression certain genes and metabolism. Taken together, these data confirm that high-intensity light induces the expression of a set of genes in clock mutants that might represent the targets for studying the potential unidentified photoreceptor(s). Taken together, these data confirm that high-intensity light induces the expression of a set of genes in clock mutants that might represent the targets for studying the potential unidentified photoreceptor(s).

Discussion
Within a certain range, the cyclical environmental changes (T-cycles) can drive rhythm with the same periodicity as the environment even when these are non-24-h cycles, which mean T-cycles mask the endogenous circadian periodicity. This range varies in different organisms [50,51]. Usually, entrainment involves the synchronization of the innate circadian clock to the environmental cues. It is important to note that entrainment here means a rhythm displaying a period that equals the length of the environmental cycle and having a stable phase angle, which takes no account of whether the innate circadian system is synchronized or not [9,14]. In Neurospora, we show that under almost all of the short LD cycles, including even very short cycles, such as LD1:1, the conidiation rhythms were driven by the T-cycles. The conidiation rhythms under short T-cycles are not autonomous as they disappeared in constant dark; instead, they are the hourglass-type rhythms.
Similar to the conidiation rhythms, the period of FRQ changed in accordance with the LD cycles when T > LD3:3. In contrast, under LD cycles with T ≤ LD3:3, FRQ level no longer oscillated with the T cycles, and instead exhibit periods close to 24 h, i.e., the masking effects faded while the FRQ rhythm free ran. The dynamic modeling results showed the superimposition of endogenous and the light-driven rhythms, which can also be observed in the western blot results of FRQ protein under very short LD cycles (T≤ LD3:3) (Figure 2). The inconsistency between conidiation rhythms and FRQ rhythms might reflect the decoupling between the oscillator and the clock-controlled output processes when T < LD3:3.
In many species, it has been demonstrated that non-24 h T-cycles decrease their fitness or adaptation [18], suggesting the adaptation of the endogenous rhythm of an organism to external cycling cues is critical for the development and growth. We show that the growth rate of Neurospora was slower under LD12:12 compared to many other conditions. However, the proportion of microconidia produced in LD12:12 prevailed over those in other conditions.
Faster growth of an organism does not necessarily mean optimal adaptation to the environment [18]. Instead, the capability to produce variable and fertile progenies may be more crucial. Neurospora microconidia function predominantly as spermatid [47], suggesting that LD12:12 might be optimal for sexual reproduction. It is likely that in nature, the circadian clock renders Neurospora better prepared to use its sexual reproduction strategies to cope with constantly changing environmental stresses. Compared to the findings in some other organisms, these results reflect that different organisms may employ different strategies for adaptation and propagation [52]. In addition, the faster growth in non-24 h conditions might help Neurospora escape from the abnormal environment.
In Neurospora, it has been claimed that WCC and VVD are responsible for most light-dependent physiological processes, if not all. Surprisingly, we show that upon removal of several photosensors, the conidiation rhythms still exist in short LD cycles. The clock mutants showed resonance to 5000 lux white light but not 1000 lux white light [14] ( Figures 3 and 5). Moreover, the conidiation rhythmicity only occurs in blue light but not red light ( Figure S5). These data suggest the presence of unidentified blue photoreceptor(s) which drives conidiation rhythmicity. The conidiation rhythms in short LD cycles are independent of either FWO or CDO, though CDO and FWD-1 are involved in modulating the conidiation rhythms in short LD cycles.
Under 1000 lux white light, clock mutants displayed no conidiation rhythms, suggesting that the circadian clock functions to maintain or promote conidiation responses to short LD cycles. These data also suggest that the circadian clock might help Neurospora to promptly adjust its conidiation according to changeable light condition in the natural environment; for instance, caused by cloud movement or shadows being cast in the day-time. Under high-intensity light, even the clock mutants showed conidiation rhythms, suggesting the unidentified photoreceptor(s) might only sense very strong light, which can confer the condition rhythmicity in clock mutants.
The light characteristics including wavelength and intensity affects fungal germination, growth, development and metabolism and light initiates adaptation in a comprehensive set of metabolic pathways and [53][54][55]. In this work, we showed that light intensity induces differential gene expression at the transcriptomic level and found that a subset of metabolic pathways was affected upon higher light exposure relative to lower light, suggesting that additional photoreceptor(s) than those that are known (e.g., WC-1, VVD and CRY) play a role in regulating light-induced metabolic adaptation ( Figure 6).
Taken together, these results provide insights into the resonance and adaptation of Neurospora circadian clock and conidiation rhythms to short LD cycles, which may benefit the understandings of the basic regulatory aspects of circadian clock. In addition, these findings might help us understand the changes and underlying mechanisms in circadian rhythms under some special conditions, e.g., the organisms living in the space laboratory exposed to the very short LD cycles (90-120-min), or the welfare of the people who are required to live and work under non-24 h work/rest schedules, e.g., maritime operations, oil mining and some other industries [56,57].
Supplementary Materials: The following are available online at https://www.mdpi.com/article/ 10.3390/jof8010027/s1. File S1 contains Table S1-S3, Figure S1-S6, and supplemental protocol and parameters for modeling and RNA-seq methods. Figure S1: Scheme for the mathematical model of Neurospora circadian rhythm, Figure S2: Conidiation rhythms of the FGSC4200 strain in LD1:1 and DD conditions, Figure S3: Coefficient analysis between the LD cycles and growth rates and between the LD cycles and microconidia proportion, Figure S4: Impacts of light sensors on the adaptation to short LD cycles, Figure S5: Conidiation rhythms of Neurospora strains in short blue light/dark (BD) and red light dark (RD) cycles, Figure S6: Growth phenotype and conidiation rhythms of indicated strains. Table S1: Neurospora crassa strains used in this study, Table S2: Parameters of the mathematical model, Table S3: the growth rates and periods of conidiation rhythms of indicated Neurospora strains. File S2 shows the list of genes induced under different light intensity. File S3 shows KEGG pathway annotation analysis of induced genes under different light conditions.