<title>Introduction</title> <p>Plants have adapted their growth and development to make use of the diurnal light/dark cycle. This is manifested at both the physiological level, with leaf and stomatal movements, modulation of growth, and photoperiodic flowering, and, at the molecular level, with diurnal cycles in enzyme and gene activity. The day/night cycling of gene expression is controlled, primarily, by light and temperature and, secondarily, by a free-running internal molecular timekeeper known as the circadian clock. The intimate connection between light signalling pathways and the circadian oscillator allows the anticipation of the environmental transitions and the measurement of day-length as an indicator of changing seasons.</p> <p>Our current understanding of the plant circadian clock derives mostly from genetic studies in <italic>Arabidopsis thaliana</italic> and rice ##REF##15208414##[1]##. Commonly, the circadian clock system is divided into three parts ##REF##9988221##[2]##: an input pathway that entrains the clock, by transmitting light or temperature signals to the core oscillator, the central oscillator (the clock), responsible for driving 24-h rhythms, and the output signals that generate the fluctuation of a wide range of molecular, biochemical and developmental responses.</p> <p>Plants have evolved several classes of photoreceptors to monitor their environmental light signals. These include red and far-red-light–absorbing phytochromes (PHYs) and UV-A/blue light–absorbing cryptochromes (CRYs) and phototropins ##REF##10441063##[3]##, ##REF##10673498##[4]##. Recent evidence shows that UV-B and green light also affect plant development, but the sensing mechanisms underlying these responses have not been elucidated. Green light has been shown to affect plant processes via both cryptochrome-dependent and cryptochrome-independent pathways ##UREF##0##[5]##.</p> <p>Genetic experiments with <italic>Arabidopsis</italic> mutants have established roles for PHYA, PHYB, PHYD, PHYE, CRY1 and CRY2 in the clock input pathway ##REF##11148293##[6]##–##REF##9425143##[8]##. Light-labile PHYA is the predominant photoreceptor at low intensities of red and blue light, whereas PHYB and CRY1 predominate at high intensities of red and blue light, respectively ##REF##9425143##[8]##. Multiple photoreceptor mutants, such as <italic>cry1 cry2</italic> or <italic>phyA phyB cry1 cry2</italic>
##REF##11148293##[6]##, retain rhythmicity and are still able to entrain the clock to a light-dark cycle, suggesting that photoreceptors can provide light input to the clock in a redundant fashion ##REF##10985392##[9]##. Also, novel putative photoreceptor families have been implicated in circadian responses, such as the ZTL/FKF/ LKP2 family ##REF##11260718##[10]##–##REF##10998191##[13]##.</p> <p>Considerable evidence exists for interaction among photoreceptors. For instance, PHYA and CRY1 directly interact at the molecular level, with CRY1 serving as a phosphorylation substrate for PHYA <italic>in</italic> vitro ##REF##9651577##[14]##. <italic>In vivo</italic>, CRY1 is phosphorylated in response to red light in a far-red reversible manner ##REF##9651577##[14]##. A <italic>cry1</italic> null mutant shows lengthened period in low intensity red or white light, and there is no additivity seen in the double <italic>phyA cry1</italic> mutant ##REF##11148293##[6]##. This suggests that CRY1 acts as a signal transduction component downstream from PHYA in the low intensity light input pathway to the clock ##REF##11181972##[15]##.</p> <p>Genetic studies have implicated two other genes, <italic>EARLY FLOWERING 3</italic> (<italic>ELF3</italic>) and <italic>GIGANTEA</italic> (<italic>GI</italic>), in light signalling to the clock. <italic>elf3</italic> loss-of-function alleles result in early flowering, hypocotyl elongation, and conditional arrhythmicity in continuous light ##REF##11402162##[16]##, ##REF##11130072##[17]##. Genetic experiments suggest substantial redundancy in ELF3 and PHYB function ##REF##10759510##[18]##. ELF3 interacts with PHYB and seems to act as a negative modulator of PHYB signalling to the clock, as <italic>ELF3</italic> overexpression both lengthens the circadian period and attenuates the resetting effects of red light pulses, whereas loss of <italic>ELF3</italic> function renders the plant hypersensitive to light signals ##REF##11402162##[16]##, ##REF##11130072##[17]##, ##REF##11402161##[19]##.</p> <p>In <italic>Arabidopsis</italic>, <italic>GI</italic> positively regulates expression of the flowering time genes <italic>CONSTANS</italic> (<italic>CO</italic>) and <italic>FLOWERING LOCUS T</italic> (<italic>FT</italic>). <italic>GI</italic> encodes a nucleoplasmically localized protein that mediates a number of responses, such as photoperiodic flowering, circadian rhythms and phytochrome/cryptochrome signalling ##REF##17098855##[20]##. The key roles played by <italic>GI</italic> are evident when analyzing the effect of <italic>gi</italic> mutants over leaf movement and gene expression rhythms of multiple clock controlled and flowering genes, including <italic>GI</italic> itself ##REF##10469647##[21]##, ##REF##10477524##[22]##.</p> <p>In-depth studies on <italic>Arabidopsis</italic> have begun to shed light on the molecular mechanisms underlying the functioning of the circadian clock. The current best candidates for <italic>Arabidopsis</italic> clock components are CIRCADIAN CLOCK ASSOCIATED 1 (CCA1) and its redundant homolog LATE ELONGATED HYPOCOTYL (LHY), which are transcription factors containing a single MYB domain ##REF##9657154##[23]##–##REF##12015970##[25]##. Furthermore, pseudo-response regulators (<italic>PRR</italic>), <italic>GI</italic>, <italic>ZEITLUPE/ADAGIO</italic> (<italic>ZTL/ADO1</italic>), <italic>LOV KELCH PROTEIN</italic> (<italic>LKP2</italic>), <italic>EARLY FLOWERING 3</italic> and <italic>4</italic> (<italic>ELF3</italic> and <italic>ELF4</italic>), <italic>LUX ARRHYTHMO/PHYTOCLOCK</italic> (<italic>LUX/PCL1</italic>), <italic>TIME FOR COFFEE</italic> (<italic>TIC</italic>), <italic>SENSITIVITY TO RED LIGHT REDUCED</italic> (<italic>SRR1</italic>) and <italic>TEJ</italic>
##REF##16761955##[26]##, ##REF##16595397##[27]##, have also been involved in the circadian machinery.</p> <p>The clock mechanism in <italic>Arabidopsis</italic> was first proposed to comprise a feedback loop, in which two partially redundant genes, <italic>LHY</italic> and <italic>CCA1</italic>, repress the expression of their activator, TOC1 ##REF##11486091##[28]##. This circuit cannot fit all experimental data ##REF##15784272##[29]##, as a short-period rhythm persists for several cycles both in <italic>lhy cca1</italic>
##REF##12007421##[30]##, ##UREF##1##[31]## and in <italic>toc1</italic> mutant plants ##REF##12509533##[32]##. Subsequently, many other clock-associated genes have been identified and incorporated into the simple model, resulting in a somewhat complicated interlocking multiloop model, comprising the feedback loop between LHY, CCA1, and TOC1, and a predicted, interlocking feedback loop involving TOC1 and a hypothetical component Y ##UREF##1##[31]##. The model was recently extended, suggesting GI as a candidate for Y and including a feedback loop between PRR7, PRR9 and LHY\CCA1, giving rise to a three loop circuit ##UREF##2##[33]##. Analysis of the three-loop network suggests that the plant clock consists of morning and evening oscillators, coupled intracellularly, which may be analogous to coupled, morning and evening clock cells in <italic>Drosophila</italic> and mouse ##UREF##2##[33]##.</p> <p>Light signals typically trigger rapid changes in the mRNA levels of transcription factors, but the position that they occupy in a putative transcriptional cascade, and the steps interposed between the photoreceptors and the first row of transcription factors, have not been fully established ##REF##16096960##[34]##. Recent work using an expressed sequence tag (EST)-based DNA microarray has suggested that nearly one-third of the genome is regulated in white light. In addition, the genome expression patterns largely overlap in 6-day-old seedlings grown under white, far-red, red, and blue light. More than 26 cellular pathways, ranging from DNA replication to transcription, metabolism, protein degradation, plant defence, and developmental regulation, have been found to be redundantly regulated by all light signals ##REF##11752374##[35]##.</p> <p>Furthermore, Schaffer and collaborators observed that 11% of genes showed differential expression at one or more of the phases tested during a light/dark cycle ##REF##11158533##[36]##. A large fraction of <italic>Arabidopsis</italic> genes that showed diurnal regulation was also circadian-regulated, as revealed by differential transcript abundance under constant light conditions ##REF##11118138##[37]##, ##REF##11158533##[36]##. Oligo-based microarray experiments on <italic>Arabidopsis</italic>
##REF##11118138##[37]## allowed the detection of circadian oscillations in mRNA abundance of 5–6% of the 8200 genes examined. In tomato, several photosynthesis-associated genes, including <italic>RBCS</italic>, <italic>LHCI</italic> and <italic>LHCII</italic>, <italic>PSAD</italic>, and <italic>OEE1</italic>, were shown to be regulated in a circadian fashion through Northern blot and nuclear run-on experiments ##REF##3208743##[38]##.</p> <p>Here, we report the characterization of temporal transcript oscillations within the tomato genome using the novel, long oligo-based TOM2 microarray. Focused Real Time RT-PCR analyses over photoreceptor gene transcripts in both wild type tomato and genotypes with altered cryptochrome gene expression provided useful information about possible functional interactions between cryptochromes and the circadian clock machinery, as well as on regulatory interactions between different photoreceptors.</p>
<title>Materials and Methods</title> <p>Standard molecular biology protocols were followed as described in Sambrook and colleagues ##UREF##6##[66]##.</p> <p>
<italic>Solanum lycopersicum</italic> (cv <italic>Moneymaker</italic>), <italic>cry1a-</italic> and transgenic <italic>CRY2-OX</italic>
##REF##11260499##[42]##, ##REF##15618424##[43]## were grown in a growth chamber for 28 days in LD conditions (16 h light-25°C/8 h dark-23°C). A light intensity of about 100 µmol m<sup>−2</sup> s<sup>−1</sup> was provided by Osram (Munich) 11–860 daylight lamps. For LL experiments, plants grown as described above for 28 days, were shifted to continuous light at the dawn of the 29th day. The aerial parts of three plants for each genotype (wt, <italic>cry1a-</italic> and <italic>CRY2-OX</italic>) were harvested at the times shown.</p> <title>Microarray analyses</title> <p>Samples were assayed on the tomato TOM2 oligo-arrays printed at the University of Arizona which comprise contain 12,160 70-mer oligonucleotide elements (<ext-link ext-link-type="uri" xlink:href="http://www.operon.com/arrays/oligosets_Tomato.php">http://www.operon.com/arrays/oligosets_Tomato.php</ext-link>).</p> <p>For each experiment, 2 µg of DNA-free total RNA was reverse-transcribed and amplified using the Aminoallyl Message Amp II kit (Ambion) following the manufacturer's instructions. 2 µg of amplified aminoallyl-modified RNA were labeled in the presence of Cy3 and Cy5 for 2 hours at room temperature. Unincorporated dyes were eliminated using RNeasy MinElute column (Qiagen) according to the manufacturer's specifications.</p> <p>200 pmoles of purified Cy3- and Cy5-labelled aRNAs were combined in a buffer containing 2× SSC, 0.08% SDS and Liquid Blocking Reagent (GE Healtcare), and were dispensed over the microarray surface, and incubated at 55°C overnight with agitation. Slides were washed in decreasing SSC concentrations and 0.1% SDS at 55°C and room temperature, respectively. The last wash was carried out in 0.1× SSC at room temperature. The hybridization and post-hybridization washes were performed using an automatic hybridization station (HybArray 12, Perkin-Elmer). Hybridized microarrays were then scanned using ScanArray Lite (Perkin-Elmer) and the resulting Cy3 and Cy5 images were analyzed with the software ScanArray Express (Perkin-Elmer) in order to determine the Cy3/Cy5 spot intensities.</p> <p>Raw hybridization signals were filtered by imposing a minimal signal/noise ratio of 2.0 and flagging the non-passed spots. In order to obtain a homogeneous dataset for all hybridized slides, we filtered microarray data imposing good quality spots to be present in at least three out of four hybridized slides (two dye-swap and two biological replicas, respectively) for each experimental point. Raw values were then normalized with the locally weighted linear regression (LOWESS) method using the 20% of data for smoothing ##UREF##7##[67]## and gene expression analysis of the array data were performed using GeneSpring version 7.3 (Agilent).</p> <p>After quality analysis and normalization (described above), we had a three-point LD time course with four microarrays per time point (2 independent biological replicates and 2 dye-swap experiments). For each of the transcripts which passed quality controls on the microarrays, a single factor ANOVA was performed across all three time points. Each time point was treated as a group, and arrays at each time point were treated as the individuals within that group. A nonadjusted ANOVA p-value of 0.05 or less was required for any particular transcript to pass the screen. After ANOVA-based statistical prescreening, genes showing equal to or more than 3-fold change in at least one of the three time points were considered diurnally regulated.</p> <p>Moreover, in order to identify genes showing major transcript regulation at dawn (ZT0), for each of the transcripts which passed quality controls on the arrays, a Student's t test was performed across ZT8, ZT16, ZT20 together and their common reference ZT0. After above mentioned analysis, of the transcripts with a Student's t test p-value of 0.05 or less, that showing at least 3-fold change simultaneously at all time points were considered differentially regulated at dawn (##TAB##0##Table 1##).</p> <p>Cluster Analyses were performed using Cluster and Treeview algorithms ##REF##9843981##[68]##. Microarray experiments have been deposited to the EBI public repository ArrayExpress (Accession number E-MEXP-1456).</p> <title>Quantitative RT-PCR</title> <p>Total RNA (1 µg) was reverse-transcribed with oligo-dT and Superscript III (Invitrogen), according to the manufacturer's instructions. First strand cDNA (5 ng) was used as template for QRT-PCR. QRT-PCR assays were carried out with gene-specific primers, using an ABI PRISM 7900HT (Applied Biosystems) and the Platinum SYBR Green master mix (Invitrogen), according to manufacturer's instructions. PCR conditions were: 5 min at 95°C, followed by 45 cycles at 95°C for 15 sec, and at 58°C for 60 sec. At the end of the PCR, the thermocycler has been programmed to generate a thermal denaturation curve of the amplified DNA and to measure the melting temperature of the PCR product(s). The shape of the melting curve indicates whether the amplified products are homogeneous and the melting temperature provides confirmation that the correct product has been specifically amplified. Relative template abundance was quantified using the relative standard curve method described in the ABI PRISM 7900HT manual and the data were normalized for the quantity of the β-actin transcript ##REF##11846609##[69]##. A serial dilution of 10-, 100-, 1000-, 10000-, and 100000-fold of each studied gene fragment was used to determine the amplification efficiency of each target and housekeeping gene. At least three PCR runs were carried out for each cDNA to serve as technical replicates and two independent experiments were carried out by using two independent RNAs for each sample. Means from two independent experiments were subjected to SEM calculation, student's t test using PAST (<ext-link ext-link-type="uri" xlink:href="http://folk.uio.no/ohammer/past/">http://folk.uio.no/ohammer/past/</ext-link>).</p>
<title>Results</title> <p>In order to identify transcripts showing temporal rhythmicity and to establish possible functional relations between photoreceptors and the circadian clock machinery in tomato, we performed extensive transcription analyses using the TOM2 microarray and using Quantitative Real-Time PCR (QRT-PCR) of the <italic>PHYA</italic>, <italic>PHYB1</italic>, <italic>PHYB2</italic>, <italic>PHYE</italic>, <italic>PHYF</italic>
##REF##8616214##[39]## and <italic>CRY1a</italic>, <italic>CRY1b</italic>, <italic>CRY2</italic> genes ##REF##10809448##[40]##, ##UREF##3##[41]##. Additional genes, already known to be regulated by the circadian clock in other plant species, including <italic>GI</italic> and <italic>LHC</italic>, were also investigated by QRT-PCR ##REF##10469647##[21]##, ##REF##3208743##[38]##.</p> <p>Wild-type (wt) tomato plants were grown under a light cycle of 16h light/8h darkness (LD), as described in <xref ref-type="sec" rid="s4">Materials and Methods</xref>, and sampled every 4h for 24 hours. Because diurnal changes of gene expression frequently reflect an underlying circadian rhythm, tomato plants entrained in LD were transferred to light constant conditions (LL), and then gene expression was monitored for additional 40 h at 4h intervals. Two genotypes with altered cryptochrome gene expression, <italic>cry1a</italic>- and <italic>CRY2-OX</italic>
##REF##11260499##[42]##, ##REF##15618424##[43]##, were also included in the experimental scheme.</p> <p>To classify the time points at which the sampling was carried out, we used Zeitgeber time (ZT), that is defined as the time in hours from the start of a normal 16 h light–8 h dark cycle ##REF##2117644##[44]##.</p> <title>Transcriptional profiling using the TOM2 microarray</title> <p>We hybridised the TOM2 microarray with target RNAs extracted from ZT0 (presumptive dawn), ZT8 (eight hours after dawn), ZT16 (presumptive dusk) and ZT20 (four hours after dusk), in LD conditions (see <xref ref-type="sec" rid="s4">Materials and Methods</xref>). The experimental design compared three time points to ZT0 used as a common reference: ZT8 vs. ZT0, ZT16 vs. ZT0 and ZT20 vs. ZT0.</p> <p>Transcripts corresponding to microarray spots which passed ANOVA test at ZT8/ZT0, ZT16/ZT0, ZT20/ZT0 and showing an expression difference greater than three-fold in at least one of the time points (see <xref ref-type="sec" rid="s4">Materials and Methods</xref>) were classified as diurnally regulated. According to this criterion, 1016 transcripts showed a diurnally regulated expression pattern, corresponding to 15% of all spots (6953) which passed quality controls (see <xref ref-type="sec" rid="s4">Materials and Methods</xref>). Compared to their expression at dawn (ZT0), the majority of the genes coding for mitochondrial and cytosolic proteins were up-regulated in the middle of the light phase (ZT8) while the genes coding for ribosomal, nuclear and thylakoid proteins were preferentially more expressed at dusk (ZT16) (##FIG##0##Figure 1A##). Many genes coding for cytoplasmic membrane, cell wall and plastid proteins showed an up-regulation in the dark phase (ZT20) (##FIG##0##Figure 1A##). Regarding their molecular function, several genes coding for proteins involved in “transporter and transferase activity” were preferentially more expressed at dusk (ZT16), while two thirds of the genes up-regulated at ZT20 are associated with transcriptional control (“transcription factor activity and DNA or RNA binding” categories) (##FIG##0##Figure 1B##).</p> <p>In order to identify representative expression patterns of day/night regulated genes, we performed cluster analysis according to similarity of gene expression patterns at ZT8/ZT0, ZT16/ZT0 and ZT20/ZT0. Among transcripts with higher relative accumulation in the middle of the light phase (ZT8) (##FIG##1##Figure 2##, cluster 1) or both at ZT8 and at the presumptive dusk (ZT16) (##FIG##1##Figure 2##, cluster 4), we found many stress-responsive genes, such as genes coding for peroxidases, caspases, salt tolerance proteins, oxygenases and some members of the WRKY family ##REF##10785665##[45]##. Other transcripts in these clusters are involved in circadian rhythms, light signal transduction and flowering - <italic>PSEUDO-RESPONSE REGULATOR 7</italic> (<italic>PRR7</italic>), <italic>FLAVIN-BINDING</italic>, <italic>KELCH REPEAT</italic>, <italic>F-BOX 1</italic> (<italic>FKF1</italic>), <italic>CONSTANS-LIKE 1</italic> (<italic>COL1</italic>), and the flowering time gene <italic>GI</italic>
##REF##14563930##[46]##, ##REF##10847687##[11]##, ##REF##10469647##[21]##. Finally, these clusters included a number of genes implicated in the light-harvesting reactions of photosynthesis that, as expected, are more expressed in the light phase.</p> <p>The converse clusters (##FIG##1##Figure 2##, cluster 3 and cluster 2) represent transcripts relatively more abundant in the dark phase (ZT20) or up-regulated at dusk (ZT16) and persisting at ZT20. These clusters contain elements involved in cellular communication and signal transduction and several transcripts involved in glycolysis/glycogenolysis and in the tricarboxylic acid cycle. We have also found some genes involved in the carotenoid pathway and several transcription factors of the WRKY, MYB, bHLH, leucine zipper, and zinc finger families.</p> <p>Transcripts with a peak or a trough at the presumptive dusk (ZT16) are grouped, respectively, in cluster 6 and cluster 5 (##FIG##1##Figure 2##). Like for cluster 1, these clusters contain several genes coding for stress responsive-proteins. In particular protease inhibitors, endotransglycosylases, the cold stress responsive protein DREB1A ##REF##9707537##[47]## and a cell death-associated protein decreased at presumptive dusk (ZT16), while transcripts peaking at dusk comprise the circadian clock related gene, <italic>ELF4</italic>
##REF##12214234##[48]##.</p> <p>Several transcripts are differentially regulated at least three-fold simultaneously at all time points (ZT8, ZT16 and ZT20) (Student's t test P≤0.05) with respect to presumptive dawn (ZT0) (##TAB##0##Table 1##). Among the 27 transcripts with higher expression at ZT0 we found elements related to light signalling and flowering (##TAB##0##Table 1B##). These include a MYB-related transcription factor, <italic>LATE ELONGATED HYPOCOTYL</italic> (<italic>LHY</italic>), a member of the YABBY family (<italic>ABNORMAL FLORAL ORGANS</italic> (<italic>AFO</italic>)), <italic>CONSTANS-LIKE</italic> 1 (<italic>COL1</italic>), <italic>SUPPRESSOR OF PHYTOCHROME A</italic> (<italic>SPA1</italic>), and <italic>EARLY LIGHT INDUCIBLE</italic> (<italic>ELIP</italic>) genes ##REF##11359606##[49]##–##REF##15475376##[52]##. Several of the 37 transcripts with lower expression at ZT0 are related to stress responses and hormone pathways (##TAB##0##Table 1A##). As expected, the <italic>PSEUDO-RESPONSE REGULATOR</italic> 7 (<italic>PRR7</italic>) transcript is also less expressed at dawn ##REF##14563930##[46]##.</p> <p>Transcription of cryptochrome and phytochrome photoreceptor genes was detected on the microarray. These genes were further assayed in depth by QRT-PCR with similar results (##FIG##2##Figure 3##). Given the higher resolution and sensitivity of the latter assay, we report the QRT-PCR data below.</p> <title>Diurnal mRNA oscillations of tomato photoreceptor genes</title> <p>We measured, by QRT-PCR, changes of mRNA accumulation of tomato phytochrome (<italic>PHYA</italic>, <italic>PHYB1</italic>, <italic>PHYB2</italic>, <italic>PHYE</italic>, <italic>PHYF</italic>) and cryptochrome (<italic>CRY1a</italic>, <italic>CRY1b</italic>, <italic>CRY2</italic>) transcripts under LD at 4-h intervals for 24 h. Although with differences in amplitude, most of the tomato photoreceptor transcripts showed diurnal fluctuations, reaching maximum expression between ZT8 and ZT12, and, with the exception of <italic>PHYB1</italic>, declined significantly during the dark period (##FIG##3##Figure 4##). The absolute expression and amplitude of the fluctuations of <italic>PHYF</italic> were significantly reduced relative to the other tomato photoreceptors (##FIG##3##Figure 4D##). The expression pattern of <italic>CRY1b</italic> was somewhat divergent from the other photoreceptors (##FIG##3##Figure 4A##). Only in early-morning (ZT4), <italic>CRY1b</italic> transcripts showed a trough (##FIG##3##Figure 4A##). <italic>CRY1b</italic> was also the most abundant photoreceptor transcript among those tested, remaining high throughout the 24h period.</p> <p>In general, comparative analyses of diurnal expression pattern of phytochrome and cryptochrome genes showed qualitatively comparable oscillation phases, though significant differences in mRNA abundance were detected throughout the 24h period (##FIG##3##Figure 4##). Concerning the overall amplitude of oscillations, photoreceptor transcripts showed modest fold-changes, ranging from about 2× (<italic>PHYF</italic>) to 7× (<italic>PHYB2</italic>) (##FIG##3##Figure 4##), compared to other diurnally regulated and circadian genes like <italic>LHC4</italic> and <italic>GI</italic> (see below).</p> <title>Circadian rhythmicity of tomato photoreceptor transcripts in LL conditions</title> <p>One of the diagnostic features of circadian rhythms is that they persist under constant light or darkness conditions. To determine whether the rhythmic fluctuations of the photoreceptor transcripts observed in LD conditions were maintained also in LL conditions, we measured their transcription in plants transferred to LL, after entraining the clock in LD. Samples were harvested at 4h intervals during a period of 40 h.</p> <p>Under LL, transcript levels of <italic>CRY1b</italic>, <italic>CRY2</italic>, and <italic>PHYB2</italic> continued to cycle, suggesting that the circadian clock controls the expression of these genes (##FIG##4##Figure 5##). <italic>PHYB1</italic> and <italic>PHYE</italic> transcripts lost any detectable oscillation, while <italic>PHYF</italic> showed increased oscillation, compared to LD conditions.</p> <title>Effects of cryptochrome gene perturbation on expression of tomato photoreceptor genes in LD and LL</title> <p>To study possible effects of cryptochrome-mediated light signals on the expression profiles of tomato photoreceptor genes, we compared mRNA levels in LD conditions in wt, in a <italic>cry1a</italic>- mutant ##REF##11260499##[42]## and in a <italic>CRY2</italic> transgenic over-expressor (<italic>CRY2-OX</italic>) ##REF##15618424##[43]##. The results indicated that loss of <italic>CRY1a</italic> as well as the over-expression of <italic>CRY2</italic> influenced the diurnal transcription profiles of several genes (##FIG##5##Figure 6##). In <italic>cry1a-</italic> or <italic>CRY2-OX</italic> plants, most tomato photoreceptor transcripts continued to cycle in LD conditions with the same phase as the wt, although with reduced or increased amplitude.</p> <p>Major alterations involved <italic>CRY1a</italic>, <italic>CRY2</italic>, <italic>PHYA</italic>, <italic>PHYB2</italic> and <italic>PHYF</italic> transcripts. Interestingly, in the <italic>cry1a-</italic> background, the peak of the (non functional) <italic>CRY1a</italic> and of the <italic>CRY2</italic> transcripts was increased by about 2-fold with respect to the wt (##FIG##5##Figure 6A##). In <italic>CRY2-OX</italic> plants, as expected for the presence of the transgene, the <italic>CRY2</italic> mRNA was about 10–15 fold more abundant relative to the wt (##FIG##5##Figure 6A##). Remarkably, the <italic>CRY2</italic> transcript, over-expressed under the control of the <italic>35S</italic> promoter still showed vigorous LD cycles (##FIG##5##Figure 6A##). <italic>PHYA</italic> transcripts were altered in a similar way in both the <italic>cry1a-</italic> and <italic>CRY2-OX</italic> backgrounds, showing a decrease of mRNA abundance especially at ZT12 (##FIG##5##Figure 6B##). A different effect was observed for <italic>PHYB2</italic>, whose oscillation was slightly increased in <italic>cry1a-</italic> and reduced in <italic>CRY2-OX</italic> plants (##FIG##5##Figure 6B##). Finally, <italic>PHYF</italic> mRNA shows a dramatic increase in the <italic>cry1a-</italic> genotype (##FIG##5##Figure 6B##).</p> <p>In order to determine possible roles of cryptochromes on the function of the tomato circadian clock, we compared the changes in the mRNA abundance of photoreceptor genes, in wt, <italic>cry1a-</italic> and <italic>CRY2-OX</italic> plants grown under LL. Our results showed that loss of <italic>CRY1a</italic> as well as over-expression of <italic>CRY2</italic> influenced the circadian transcription profiles of a number of photoreceptor genes, including <italic>CRY1a</italic>, <italic>CRY2</italic>, <italic>PHYA</italic>, <italic>PHYB2</italic> and <italic>PHYF</italic> (##FIG##6##Figure 7##). As already observed in LD, the rhythm in <italic>CRY2</italic> expression was not affected by <italic>CRY2</italic> over-expression, although the transcript levels were 10–15 fold more abundant relative to those observed in wt (##FIG##6##Figure 7A##). Circadian oscillations of <italic>PHYA</italic> appeared to be repressed in <italic>CRY2-OX</italic>, following transfer to LL conditions, with a low amplitude oscillation remaining in <italic>cry1a-</italic> (##FIG##6##Figure 7B##). <italic>CRY2</italic> over-expression dramatically reduced the amplitude of <italic>PHYB2</italic> oscillations (##FIG##6##Figure 7B##). Finally, as for LD experiments, the <italic>PHYF</italic> transcription was dramatically increased under LL conditions in the <italic>cry1a-</italic> genotype (##FIG##6##Figure 7B##).</p> <title>Oscillation of tomato GI and LHC4 mRNAs and the effect of cryptochromes</title> <p>In LD conditions, <italic>GI</italic> transcripts oscillated about 800-fold, with a peak at ZT12, and a trough between ZT20 and ZT0 (##FIG##7##Figure 8A##). The <italic>LHC4</italic> peak occurred 4 h earlier (ZT8) (##FIG##7##Figure 8A##), and the trough 12 h later (ZT20) than <italic>GI</italic> (##FIG##7##Figure 8A##), with an 84-fold difference in transcript abundance. Interestingly, <italic>LHC4</italic> transcripts increased slightly in darkness from ZT20 to ZT24, showing an anticipation of “light-on” that is typical of circadian-regulated, particularly of <italic>LHC</italic> genes (##FIG##7##Figure 8A##) ##REF##11743105##[53]##.</p> <p>Under LL, <italic>GI</italic> and <italic>LHC4</italic> mRNA levels continued to cycle, although with a much reduced amplitude (##FIG##7##Figure 8B##), showing that transcription of these genes is partially controlled by the endogenous circadian clock. LL conditions affected both the amplitude/phase of oscillations as well as the abundance of the mRNAs (##FIG##7##Figure 8B##).</p> <p>In the <italic>cry1a-</italic> genotype, both <italic>GI</italic> and <italic>LHC4</italic> transcripts showed a reduction of the peak at ZT12. In addition, a slight phase shift was evident which anticipated the transcription peak at ZT4 (##FIG##7##Figure 8A##), while no major alterations were observed in the <italic>CRY2-OX</italic> genotype.</p> <p>Under LL conditions, neither <italic>CRY1a</italic> loss nor <italic>CRY2</italic> over-expression have dramatic effect on <italic>GI</italic> and LHC4 expression (##FIG##7##Figure 8B##).</p>
<title>Discussion</title> <title>Diurnal expression of the tomato genome</title> <p>Using the TOM2 microarray, we identified 1016 diurnally regulated genes, corresponding to about 15% of the spots that passed quality controls. Though by using high cutoff (&gt;3x) values we might have excluded a certain number of transcripts cycling with lower amplitude, it is evident that in tomato, like in <italic>Arabidopsis</italic>, diurnal rhythms in gene expression affect a large portion of the transcriptome ##REF##11158533##[36]##. The majority of dirunally regulated genes showed a peak at midday (ZT8), while the other transcription peaks appeared evenly distributed at the other time points, supporting the occurrence of highly coordinated and alternated metabolic processes (see supplementary ##SUPPL##0##Data S1##).</p> <p>Given the cyclic nature of many physiological processes driven by photo- and thermocycles ##REF##18248097##[54]##, it is expected that the majority of transcripts involved in the biosynthesis of mitochondrial and cytosolic proteins peak in the middle of the light phase (ZT8) (##FIG##0##Figure 1A##). This can be attributed to the fact that the biosynthetic processes correlated to photosynthesis and energy metabolism are usually more active in light hours. Similarly, the fact that several transcripts coding for proteins involved in transport, transferase activity and in the transcription control machinery were also abundant at dusk (ZT16) and during the night (ZT20) (##FIG##0##Figure 1B##) indicates that, during the hours of darkness, synthesis of these proteins is still active.</p> <p>Several transcripts with higher levels during daylight (ZT8 to ZT16), grouping in clusters 2 and 4 (##FIG##1##Figure 2##), encode for protein elements involved in photosynthesis and stress response elements. The latter, that include a number of transcription factors – MYB, WRKY, bHLH, salt tolerance proteins, peroxidases, oxygenases and others, could have a major role to adapt tomato plants to day conditions, such as excess of light and higher temperatures.</p> <p>Conversely, several transcripts relatively more abundant in during the dark phase (ZT16 to ZT20), grouped in clusters 2 and 3 (##FIG##1##Figure 2##), are related to biochemical processes occurring in darkness. We found genes involved in nitrogen and sulfur assimilation, as well as key genes involved cell wall loosening, such as xyloglucan endotransglycosylase and expansin. Most of these genes are already known to be down-regulated during daylight in <italic>Arabidopsis</italic>
##REF##11752374##[35]##. Thus, it seems that these processes are diurnally regulated in both <italic>Arabidopsis</italic> and tomato.</p> <p>Plants need protection from the damaging effects of reactive oxygen species generated by the excess of light; in green tissues, carotenoids prevent the chlorophyll-photosensitized formation of highly destructive singlet oxygen by quenching the chlorophyll triplet states, scavenging reactive oxygen species. Furthermore, they have an antenna function and transfer the energy of absorbed light at the singlet excited state level to the chlorophyll system for the execution of photosynthesis ##REF##11884678##[55]##. Structural genes of the pathway appear to be diurnally regulated, although with different phases (##FIG##1##Figure 2##): <italic>ZEAXANTHIN EPOXIDASE</italic> (<italic>ZEP</italic>) is found in cluster 1, <italic>VIOLAXANTHIN DE-EPOXIDASE</italic> (<italic>VDE</italic>) and <italic>PHYTOENE SYNTHASE</italic> (<italic>PSY</italic>) in cluster 2, <italic>β-CAROTENE HYDROXYLASE</italic> (<italic>CHY</italic>) in cluster 3, <italic>DEOXYXYLULOSE 5-PHOSPHATE SYNTHASE</italic> (<italic>DXS</italic>) in cluster 4. Although light regulation of carotenoid gene transcription is a well known phenomenon ##REF##9351247##[56]##, diurnal rhythms in gene expression have been reported, to date, only for <italic>ZEP</italic>
##REF##9808747##[57]##. Our data indicate that these rhythms are widespread in transcripts encoding carotenoid biosynthesis enzymes. This observation remains to be interpreted, in combination with data on the diurnal abundance and activity of the corresponding enzymes.</p> <p>A good example of a possible coordinated response of tomato plants to abiotic stresses is given by the cyclic transcript oscillations of the DREB1A and DREB2 transcription factors (##FIG##1##Figure 2##, clusters 1 and 5). In <italic>Arabidopsis</italic>, <italic>DREB1A</italic> gene and its two homologs are implicated in response to low-temperature stress, in a manner independent of ABA, and its transcripts peak during the presumptive day, whereas expression of the <italic>DREB2A</italic> gene and its single homolog was induced by dehydration ##REF##9707537##[47]##, ##REF##17965713##[58]##. Thus, it is plausible that the observed increment of <italic>DREB</italic> transcripts at ZT8 in tomato plants under LD conditions provide appropriate defense against changing temperature and dehydration occurring during light hours. The expression pattern of <italic>DREB1A</italic> evidenced a consistent increment of its transcripts at dark (ZT20) (##FIG##1##Figure 2##, cluster 5), possibly due to the decrease of ambient temperature at the presumptive dusk (ZT16) (see <xref ref-type="sec" rid="s4">Material and Methods</xref>).</p> <p>Several tomato homologues of the genes involved in the circadian clock feedback-loop in <italic>Arabidopsis</italic>
##UREF##4##[59]## were found to oscillate in a similar phase in tomato: the morning element <italic>LHY</italic> was up-regulated ad dawn (##TAB##0##Table 1##); while <italic>PRR7</italic>, thought to establish a negative loop with <italic>CCA1/LHY</italic>, was more expressed during daylight (ZT8-ZT16) and down-regulated at dawn (##FIG##1##Figure 2##, cluster 4 and ##TAB##0##Table 1##). <italic>ELF4</italic> and <italic>GI</italic> (discussed more in detail below), which are putatively involved in feedback-loops with <italic>CCA1/LHY</italic> and <italic>TOC1/LUX</italic>, respectively ##UREF##2##[33]##, were accordingly more expressed around dusk (##FIG##1##Figure 2##, clusters 6 and 4). These results suggest that the basic molecular machinery of the circadian clock is conserved in higher plants. Furthermore, the fact that a number of other elements, like <italic>FKF1</italic>
##REF##10847687##[11]## and <italic>SPA1</italic>
##REF##16709190##[50]## related to the input/output signalling of the <italic>Arabidopsis</italic> circadian clock, but also involved in flowering and light transduction, showed similar transcript fluctuations in tomato (##FIG##1##Figure 2##, cluster 4 and ##TAB##0##Table 1##) suggests that molecular interactions between the clock core and input/output pathways are also partially conserved. However, it must be considered that our dataset is largely incomplete and does not represent the actual complexity of transcript network interactions described in <italic>Arabidopsis</italic>.</p> <title>Temporal modulation of cryptochrome and phytochrome transcripts</title> <p>Previous experiments in <italic>Arabidopsis</italic> have established a fundamental role of phytochromes and cryptochromes in providing light input to the plant circadian clock ##REF##11181972##[15]##, ##REF##11743105##[53]##. In tomato, as for <italic>Arabidopsis</italic>, we observed a bi-directional regulatory crosstalk between the clock machinery and photoreceptors which allowed the latter to determine significant changes on the temporal transcription pattern of genes under the control of the first.</p> <p>As seen in <italic>Arabidopsis</italic>
##REF##11743105##[53]##, tomato <italic>PHY</italic> and <italic>CRY</italic> genes followed a diurnal rhythm and exhibited maximum expression in the light phase (##FIG##3##Figure 4A-D##). Tomato photoreceptor transcripts, except for <italic>CRY1b</italic>, appeared to be synchronized and peaked during the presumptive afternoon, (##FIG##3##Figure 4A-D##). By contrast, in <italic>Arabidopsis</italic> gene expression trends are different between photo-stable and photo-labile photoreceptors. Indeed, light-stable photoreceptors are highly expressed at the beginning (<italic>PHYC</italic>, <italic>PHYD</italic>, and <italic>PHYE</italic>) or in the first half (<italic>PHYB</italic> and <italic>CRY1</italic>) of the light phase, while photo-labile <italic>PHYA</italic> and <italic>CRY2</italic> reach their maximum transcript abundance close to the end of the light interval. Unfortunately, data on the photo-stability of tomato photoreceptors are not yet available. The massive accumulation, in late afternoon, of most of the tomato photoreceptor transcripts, including <italic>CRY-DASH</italic>
##REF##16876787##[60]##, might reflect the different photoperiodic behaviour of the two species (long-day for <italic>Arabidopsis</italic>, day-neutral for cultivated tomato).</p> <p>The temporal regulation of <italic>CRY1b</italic> expression, whose mRNA was the most abundant among the analyzed photoreceptors, did not show remarkable fluctuations during the day, and was quite different from that of the other <italic>CRY</italic> genes (##FIG##3##Figure 4A##). Despite its high sequence similarity with <italic>CRY1a</italic>, this gene is not yet functionally characterized. The similarities of the expression patterns of <italic>PHYA</italic>, <italic>CRY1a</italic> and <italic>CRY2</italic> genes in both LD and LL (##FIG##3##Figure 4B## and ##FIG##4##Figure 5B##), namely high levels of expression in the second part of the day (ZT8-16) and very low transcript abundances during the night, could be potentially related to overlapping functions and/or cooperation in their physiological roles.</p> <p>
<italic>PHYB1</italic> was the most abundant among phytochrome transcripts, followed by <italic>PHYB2</italic>, <italic>PHYE</italic> and <italic>PHYA</italic>, while <italic>PHYF</italic> is by far the least expressed phytochrome transcript in tomato green tissues (##FIG##3##Figure 4##). In LD, the expression peak of all phytochrome genes was between ZT8-ZT12 with no major phase differences. The amplitude of the oscillations was quite modest, with the sole exception of <italic>PHYB2</italic> that showed a 7-fold difference between through and peak transcript levels (##FIG##3##Figure 4##). These data contrast with a previous report ##REF##9670560##[61]## which evidenced a phase shift of about 10 hours between diurnal transcription rhythms of <italic>PHYB1</italic> and <italic>PHYB2</italic>. However, it must be taken in account that the authors used a quite different experimental set-up, with tomato plants grown in greenhouse and without supplemental illumination. This specific timing of transcript accumulation suggests that photoreceptor-mediated input signalling to the clock machinery may be particularly synchronized in tomato.</p> <p>The fact that under LL all tomato cryptochromes plus <italic>PHYA</italic>, <italic>PHYB2</italic> and <italic>PHYF</italic> seem to keep their oscillations following a period close to 24 hours, though with lower amplitude and minor changes in the phase of the peaks (##FIG##4##Figure 5A-D##), hints that a circadian clock regulates the expression of these photoreceptors, as seen in the <italic>Arabidopsis</italic> closest homologs, <italic>PHYA</italic>, <italic>PHYD</italic> and <italic>PHYF</italic>
##REF##11743105##[53]##. In contrast, <italic>PHYB1</italic> and <italic>PHYE</italic> lose their rhythmicity in LL (##FIG##4##Figure 5A, C##), while the most closely related <italic>Arabidopsis</italic> homologs, <italic>PHYB</italic> and <italic>PHYE</italic> continue to cycle in LL with a peak at the beginning or in the first one-half of the light phase ##REF##11743105##[53]##. The different regulation in the two plants could reflect the different functional organization of the photoreceptor gene families. Differently from <italic>Arabidopsis</italic>, tomato flowering is day-neutral.</p> <p>In LL, early into the presumptive night (ZT20), an increment in the <italic>CRY2</italic>, <italic>PHYB2</italic> and <italic>PHYF</italic> transcript levels with respect to the correspondent LD point was evident (##FIG##4##Figure 5##). This difference may be explained through postulating direct activation by light. The actual transcript levels appear to be then partially restored to the “normal” light/dark oscillation; this is possibly caused by some feedback action mediated by the clock machinery. This hypothetical feedback action is consistent with the model proposed by Tòth and colleagues ##REF##11743105##[53]##, in which the photoreceptors send the “light-on” signal to the clock core, and the core regulates their expression, forming a regulatory loop. This regulatory loop could serve to increase the perception of resetting light signals at the right times, and to neutralize signals from non-predictable environmental cues, which could cause an incorrect resetting of the circadian clock.</p> <title>Influence of cryptochromes on temporal transcription of photoreceptor genes</title> <p>In LD conditions, alterations in cryptochrome gene expression caused a minor increase of cryptochrome transcripts. The peak expression of <italic>CRY1a</italic> was incremented in <italic>cry1a-</italic> plants (##FIG##5##Figure 6A##); this could be the effect of an auto-regulatory feedback mechanism mediated by CRY1a and repressing the transcription of its own gene.</p> <p>The increment of <italic>CRY2</italic> transcripts in <italic>CRY2</italic>-<italic>OX</italic> transgenic tomato is expected. However, the <italic>CRY2</italic> transcript, under the control of the 35S promoter, continues to cycle in both LD and LL (##FIG##5##Figure 6A## and ##FIG##6##Figure 7A##). This strongly suggests that at least part of the <italic>CRY2</italic> oscillations are post-transcriptional. A similar situation has been observed in a transgenic line over-expressing <italic>GI</italic>
##REF##16006578##[62]##. To our knowledge, this is the first time that posttranscriptional (diurnal and circadian) oscillations are described in <italic>CRY2</italic>.</p> <p>Our data demonstrate that cryptochromes regulate phytochrome transcript levels, resulting in changes in transcript abundance, phase and cycling amplitude. Additionally, our data suggest that Cryptochromes 1 and 2 act cooperatively in repressing the transcription of <italic>PHYA</italic> and antagonistically on the transcription of <italic>PHYB2</italic>, which is promoted by CRY1a and repressed by CRY2 (##FIG##5##Figure 6B## and ##FIG##8##Figure 9##). In <italic>Arabidopsis</italic>, there is evidence for a direct interaction between PHYA and CRY1, with PHYA mediating a light-dependent phosphorylation of CRY1 ##REF##9651577##[14]##, and between PHYB and CRY2, with the CRY2 probably suppressing PHYB signaling ##REF##11089975##[63]##. Furthermore, in <italic>Arabidopsis</italic> CRY1 operates as a signal transduction component downstream of PHYA in light input to the clock ##REF##11148293##[6]##. In tomato, an additional level of suppression of PHYB signalling could be represented by the repressive action of CRY1a and CRY2 on <italic>PHYB2</italic> transcript levels (##FIG##5##Figure 6B## and ##FIG##8##Figure 9##). Another interaction is the approximately 3–10 fold increase of <italic>PHYF</italic> transcripts at all time points in plants lacking a functional CRY1a.</p> <p>Under LL conditions, transcriptional oscillations often became more perturbed and were sometimes difficult to interpret. However, the evident arrhythmicity of <italic>PHYA</italic> and <italic>PHYB2</italic> , but not <italic>CRY1a</italic> transcripts caused by the over-expression of CRY2 (##FIG##6##Figure 7B##) is quite intriguing and suggests that this condition specifically disrupts the output signal from the clock to <italic>PHYA and PHYB2</italic> (##FIG##8##figure 9##).</p> <title>Transcription rhythms of <italic>LHC4</italic> and <italic>GI</italic> are under cryptochrome control</title> <p>In tomato wt plants, <italic>GI</italic> and <italic>LHC4</italic> transcripts cycled similarly to their putative counterparts in <italic>Arabidopsis</italic>
##REF##10469647##[21]##, ##UREF##5##[64]## with strong diurnal oscillations (about 800 and 90-fold, respectively). Lack of a functional CRY1a decreased the amplitude of the diurnal oscillation of both <italic>GI</italic> and <italic>LHC4</italic>, meaning that both genes, directly or indirectly, are activated by CRY1a (##FIG##7##Figure 8A## and ##FIG##8##Figure 9##). A recent report showed that <italic>Arabidopsis</italic> mutant <italic>cry1-cry2</italic> plants displayed a severely reduced <italic>GI</italic> response to blue light, while CRY2 had no affect on the diurnal transcription of <italic>GI</italic>
##REF##16775702##[65]##. In agreement with these data, our results demonstrate that CRY1a plays a major role in the activation of tomato <italic>GI</italic> under high fluence white light. If we accept the hypothesis of GI as the “factor Y” in an interlocked feedback loop through light affecting TOC1 expression ##UREF##2##[33]##, we must conclude that CRY1a plays a major role in the input to the tomato circadian clock.</p>
<p>Conceived and designed the experiments: PF LL DWG GG GP. Performed the experiments: PF LL FC GP. Analyzed the data: PF LL FC GP. Contributed reagents/materials/analysis tools: PF LL DWG GP. Wrote the paper: PF LL GG GP.</p> <title>Background</title> <p>Circadian clocks are internal molecular time-keeping mechanisms that provide living organisms with the ability to adjust their growth and physiology and to anticipate diurnal environmental changes. Circadian clocks, without exception, respond to light and, in plants, light is the most potent and best characterized entraining stimulus. The capacity of plants to respond to light is achieved through a number of photo-perceptive proteins including cryptochromes and phytochromes. There is considerable experimental evidence demonstrating the roles of photoreceptors in providing light input to the clock.</p> <title>Methodology</title> <p>In order to identify genes regulated by diurnal and circadian rhythms, and to establish possible functional relations between photoreceptors and the circadian clock in tomato, we monitored the temporal transcription pattern in plants entrained to long-day conditions, either by large scale comparative profiling, or using a focused approach over a number of photosensory and clock-related genes by QRT-PCR. In parallel, focused transcription analyses were performed in <italic>cry1a</italic>- and in <italic>CRY2-OX</italic> tomato genotypes.</p> <title>Conclusions</title> <p>We report a large series of transcript oscillations that shed light on the complex network of interactions among tomato photoreceptors and clock-related genes. Alteration of cryptochrome gene expression induced major changes in the rhythmic oscillations of several other gene transcripts. In particular, over-expression of <italic>CRY2</italic> had an impact not only on day/night fluctuations but also on rhythmicity under constant light conditions. Evidence was found for widespread diurnal oscillations of transcripts encoding specific enzyme classes (e.g. carotenoid biosynthesis enzymes) as well as for post-transcriptional diurnal and circadian regulation of the <italic>CRY2</italic> transcript.</p>
<title>Supporting Information</title>
<p>We are grateful to Leonardo Giliberto for valuable discussions.</p>
<fig id="pone-0002798-g001" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0002798.g001</object-id><label>Figure 1</label><caption><title>The impact of diurnal transcript oscillations on different categories of genes based on Gene Ontology.</title><p>The percentages of diurnally regulated genes classified by cellular compartment (A) and by molecular function (B). Black bars: ZT20/ZT0. White bars: ZT16/ZT0. Grey bars: ZT8/ZT0.</p></caption></fig> <fig id="pone-0002798-g002" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0002798.g002</object-id><label>Figure 2</label><caption><title>Clustering analysis of diurnally regulated genes (P ≤ 0.05).</title><p>Idealized graphs, representing patterns of expression at ZT8/ZT0, ZT16/ZT0 and ZT20/ZT0, are shown.</p></caption></fig> <fig id="pone-0002798-g003" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0002798.g003</object-id><label>Figure 3</label><caption><title>Pearson's correlations between gene expression levels of photoreceptors determined by QRT-PCR and the TOM2 microarray.</title></caption></fig> <fig id="pone-0002798-g004" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0002798.g004</object-id><label>Figure 4</label><caption><title>Diurnal oscillations of Cryptochrome and Phytochrome transcripts analyzed by QRT-PCR in tomato plants grown in LD conditions.</title><p>Results are presented as a proportion of the highest value after normalization with β-actin. Open and closed bars along the horizontal axis represent light and dark periods, respectively. Time points are measured in hours from dawn (zeitgeber Time [ZT]). Data shown are the average of two biological replicates, with error bars representing SEM. Time points of <italic>CRY1a</italic>, <italic>CRY2</italic>, <italic>PHYA</italic>, <italic>PHYB1</italic>, <italic>PHYB2</italic>, <italic>PHYE</italic> and <italic>PHYF</italic> transcripts significantly different from the corresponding ones of the <italic>CRY1b</italic> gene (Student's t test, P ≤ 0.05) are marked with an O. Time points significantly different from the highest transcription value (Student's t test, P ≤ 0.05) are marked with an X.</p></caption></fig> <fig id="pone-0002798-g005" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0002798.g005</object-id><label>Figure 5</label><caption><title>Circadian oscillations of Cryptochrome and Phytochrome transcripts in tomato plants entrained in LD conditions and then transferred to LL.</title><p>Results are presented as a proportion of the highest value after normalization with β-actin. Open and hatched bars along the horizontal axis represent light and subjective night periods, respectively. Time points are measured in hours from dawn (zeitgeber Time [ZT]). Data shown are the average of two biological replicates, with error bars representing SEM. Time points of <italic>CRY1a</italic>, <italic>CRY2</italic>, <italic>PHYA</italic>, <italic>PHYB1</italic>, <italic>PHYB2</italic>, <italic>PHYE</italic> and <italic>PHYF</italic> transcripts significantly different from the corresponding ones of the <italic>CRY1b</italic> gene (Student's t test, P ≤ 0.05) are marked with an O. Time points significantly different from the highest transcription value (Student's t test, P ≤ 0.05) are marked with an X.</p></caption></fig> <fig id="pone-0002798-g006" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0002798.g006</object-id><label>Figure 6</label><caption><title>Effect of <italic>CRY1a</italic> loss-of-function and <italic>CRY2</italic> over-expression on diurnal expression of tomato cryptochrome (A) and phytochrome (B) genes.</title><p>Wt, <italic>cry1a</italic>- and <italic>CRY2-OX</italic> tomato plants were grown under LD conditions. The abundance of the mRNAs were measured by QRT-PCR. Results are presented as a proportion of the highest value after normalization with β-actin. Open and closed bars along the horizontal axis represent light and dark periods, respectively. Time points are measured in hours from dawn (zeitgeber Time [ZT]). An additional panel depicts <italic>CRY2</italic> transcript values in wt and <italic>cry1a</italic>- genotypes to avoid the masking effect of <italic>CRY2-OX</italic> values. Data shown are the average of two biological replicates, with error bars representing SEM. Circles (O) indicate time points of <italic>CRY2-OX</italic> and <italic>cry1a-</italic> genotypes, significantly different from the corresponding ones in wt genotype (Student's t test, P ≤ 0.05). For each genotype X indicate time points significantly different from the highest transcription value (Student's t test, P ≤ 0.05).</p></caption></fig> <fig id="pone-0002798-g007" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0002798.g007</object-id><label>Figure 7</label><caption><title>Effect of <italic>CRY1a</italic> loss-of function and <italic>CRY2</italic> over-expression on circadian expression of tomato cryptochrome (A) and phytochrome (B) genes in LL.</title><p>Wt, <italic>cry1a</italic>- and <italic>CRY2-OX</italic> tomato plants were entrained under LD cycles and then transferred to LL. The abundance of the mRNAs were measured by QRT-PCR. Results are presented as a proportion of the highest value after normalization with β-actin. Open and hatched bars along the horizontal axis represent light and subjective night periods, respectively. Time points are measured in hours from dawn (zeitgeber Time [ZT]). An additional panel depicts <italic>CRY2</italic> transcript values in wt and <italic>cry1a</italic>- genotypes to avoid the masking effect of <italic>CRY2-OX</italic> values. Data shown are the average of two biological replicates, with error bars representing SEM. Circles (O) indicate time points of <italic>CRY2-OX</italic> and <italic>cry1a-</italic> genotypes, significantly different from the corresponding ones in wt genotype (Student's t test, P ≤ 0.05). For each genotype, X indicate time points significantly different from the highest transcription value (Student's t test, P ≤ 0.05).</p></caption></fig> <fig id="pone-0002798-g008" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0002798.g008</object-id><label>Figure 8</label><caption><title>Effect of CRY1a loss and <italic>CRY2</italic> over-expression on light induced transcription of tomato <italic>GI</italic> and <italic>LHC4</italic> genes.</title><p>Wt, <italic>cry1a</italic>- and <italic>CRY2-OX</italic> tomato plants were grown under LD (A) and LL (B) conditions. The abundance of the mRNAs of <italic>GI</italic> and <italic>LHC4</italic> genes were measured by QRT-PCR. Results are presented as a proportion of the highest value after normalization with β-actin. Open, closed and hatched bars along the horizontal axis represent light, dark and subjective night periods, respectively. Time points are measured in hours from dawn (zeitgeber Time [ZT]). Data shown are the average of two biological replicates, with error bars representing SEM. Circles (O) indicate time points of <italic>CRY2-OX</italic> and <italic>cry1a-</italic> genotypes, significantly different from the corresponding ones in wt (Student's t test, P ≤ 0.05). For each genotype, X indicate time points significantly different from the highest transcription value (Student's t test, P ≤ 0.05).</p></caption></fig> <fig id="pone-0002798-g009" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0002798.g009</object-id><label>Figure 9</label><caption><title>A model for genetic interactions among tomato CRY1a and CRY2 proteins and tomato photoreceptor, <italic>LHC4</italic> and <italic>GI</italic> genes, deduced from transcription experiments.</title><p>The arrows represent a stimulatory effect, and the lines terminated with a bar represent an inhibitory effect. Positive and negative effects mediated by the circadian machinery are represented by the clock symbol.</p></caption></fig>
<table-wrap id="pone-0002798-t001" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0002798.t001</object-id><label>Table 1</label><caption><title>Transcripts up-regulated (A) and down-regulated (B) at ZT0</title></caption><alternatives><table frame="hsides" rules="groups"><colgroup span="1"><col align="left" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/></colgroup><thead><tr><td colspan="2" align="left" rowspan="1">zt8/zt0</td><td colspan="2" align="left" rowspan="1">zt16/zt0</td><td colspan="2" align="left" rowspan="1">zt20/zt0</td><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">A</td></tr><tr><td align="left" rowspan="1" colspan="1">Ratio</td><td align="left" rowspan="1" colspan="1">P-value</td><td align="left" rowspan="1" colspan="1">Ratio</td><td align="left" rowspan="1" colspan="1">P-value</td><td align="left" rowspan="1" colspan="1">Ratio</td><td align="left" rowspan="1" colspan="1">P-value</td><td align="left" rowspan="1" colspan="1">Gene_ID</td><td align="left" rowspan="1" colspan="1">Annotation</td></tr></thead><tbody><tr><td align="left" rowspan="1" colspan="1">40,676</td><td align="left" rowspan="1" colspan="1">1,95E-02</td><td align="left" rowspan="1" colspan="1">37,195</td><td align="left" rowspan="1" colspan="1">4,81E-02</td><td align="left" rowspan="1" colspan="1">3,877</td><td align="left" rowspan="1" colspan="1">1,54E-02</td><td align="left" rowspan="1" colspan="1">sgn-U216348</td><td align="left" rowspan="1" colspan="1">APRR7 (APRR1/TOC1 family)</td></tr><tr><td align="left" rowspan="1" colspan="1">24,081</td><td align="left" rowspan="1" colspan="1">4,26E-03</td><td align="left" rowspan="1" colspan="1">46,603</td><td align="left" rowspan="1" colspan="1">1,13E-03</td><td align="left" rowspan="1" colspan="1">6,864</td><td align="left" rowspan="1" colspan="1">2,26E-03</td><td align="left" rowspan="1" colspan="1">sgn-U227341</td><td align="left" rowspan="1" colspan="1">short-chain dehydrogenase/reductase family protein</td></tr><tr><td align="left" rowspan="1" colspan="1">23,893</td><td align="left" rowspan="1" colspan="1">2,95E-03</td><td align="left" rowspan="1" colspan="1">5,277</td><td align="left" rowspan="1" colspan="1">2,44E-05</td><td align="left" rowspan="1" colspan="1">5,734</td><td align="left" rowspan="1" colspan="1">8,85E-04</td><td align="left" rowspan="1" colspan="1">sgn-U243036</td><td align="left" rowspan="1" colspan="1">unknown protein</td></tr><tr><td align="left" rowspan="1" colspan="1">21,624</td><td align="left" rowspan="1" colspan="1">1,73E-02</td><td align="left" rowspan="1" colspan="1">11,585</td><td align="left" rowspan="1" colspan="1">1,59E-02</td><td align="left" rowspan="1" colspan="1">3,264</td><td align="left" rowspan="1" colspan="1">1,18E-03</td><td align="left" rowspan="1" colspan="1">sgn-U214829</td><td align="left" rowspan="1" colspan="1">expressed protein</td></tr><tr><td align="left" rowspan="1" colspan="1">17,494</td><td align="left" rowspan="1" colspan="1">5,50E-04</td><td align="left" rowspan="1" colspan="1">3,399</td><td align="left" rowspan="1" colspan="1">3,04E-05</td><td align="left" rowspan="1" colspan="1">3,174</td><td align="left" rowspan="1" colspan="1">4,03E-03</td><td align="left" rowspan="1" colspan="1">sgn-U221348</td><td align="left" rowspan="1" colspan="1">expressed protein</td></tr><tr><td align="left" rowspan="1" colspan="1">17,069</td><td align="left" rowspan="1" colspan="1">1,99E-02</td><td align="left" rowspan="1" colspan="1">6,311</td><td align="left" rowspan="1" colspan="1">1,08E-02</td><td align="left" rowspan="1" colspan="1">3,050</td><td align="left" rowspan="1" colspan="1">2,09E-02</td><td align="left" rowspan="1" colspan="1">sgn-U231435</td><td align="left" rowspan="1" colspan="1">O-methyltransferase family 2 [Arabidopsis thaliana]</td></tr><tr><td align="left" rowspan="1" colspan="1">16,724</td><td align="left" rowspan="1" colspan="1">1,77E-02</td><td align="left" rowspan="1" colspan="1">8,194</td><td align="left" rowspan="1" colspan="1">8,13E-03</td><td align="left" rowspan="1" colspan="1">3,486</td><td align="left" rowspan="1" colspan="1">6,39E-03</td><td align="left" rowspan="1" colspan="1">sgn-U213589</td><td align="left" rowspan="1" colspan="1">protease inhibitor/seed storage/lipid transfer protein (LTP) family</td></tr><tr><td align="left" rowspan="1" colspan="1">15,621</td><td align="left" rowspan="1" colspan="1">3,66E-02</td><td align="left" rowspan="1" colspan="1">9,834</td><td align="left" rowspan="1" colspan="1">1,67E-02</td><td align="left" rowspan="1" colspan="1">8,262</td><td align="left" rowspan="1" colspan="1">3,62E-03</td><td align="left" rowspan="1" colspan="1">sgn-U214471</td><td align="left" rowspan="1" colspan="1">hydrolase, alpha/beta fold family</td></tr><tr><td align="left" rowspan="1" colspan="1">15,516</td><td align="left" rowspan="1" colspan="1">8,20E-03</td><td align="left" rowspan="1" colspan="1">10,086</td><td align="left" rowspan="1" colspan="1">2,04E-02</td><td align="left" rowspan="1" colspan="1">5,105</td><td align="left" rowspan="1" colspan="1">2,61E-02</td><td align="left" rowspan="1" colspan="1">sgn-U214470</td><td align="left" rowspan="1" colspan="1">hydrolase, alpha/beta fold family</td></tr><tr><td align="left" rowspan="1" colspan="1">15,269</td><td align="left" rowspan="1" colspan="1">1,02E-03</td><td align="left" rowspan="1" colspan="1">10,456</td><td align="left" rowspan="1" colspan="1">1,82E-03</td><td align="left" rowspan="1" colspan="1">5,413</td><td align="left" rowspan="1" colspan="1">6,03E-03</td><td align="left" rowspan="1" colspan="1">sgn-U218302</td><td align="left" rowspan="1" colspan="1">wound-responsive protein -related</td></tr><tr><td align="left" rowspan="1" colspan="1">12,718</td><td align="left" rowspan="1" colspan="1">1,21E-03</td><td align="left" rowspan="1" colspan="1">15,622</td><td align="left" rowspan="1" colspan="1">1,29E-02</td><td align="left" rowspan="1" colspan="1">6,601</td><td align="left" rowspan="1" colspan="1">2,23E-02</td><td align="left" rowspan="1" colspan="1">sgn-U233539</td><td align="left" rowspan="1" colspan="1">unknown protein</td></tr><tr><td align="left" rowspan="1" colspan="1">12,709</td><td align="left" rowspan="1" colspan="1">2,02E-02</td><td align="left" rowspan="1" colspan="1">12,884</td><td align="left" rowspan="1" colspan="1">2,01E-03</td><td align="left" rowspan="1" colspan="1">3,550</td><td align="left" rowspan="1" colspan="1">1,60E-03</td><td align="left" rowspan="1" colspan="1">sgn-U216720</td><td align="left" rowspan="1" colspan="1">cellulose synthase catalytic subunit</td></tr><tr><td align="left" rowspan="1" colspan="1">11,701</td><td align="left" rowspan="1" colspan="1">1,19E-02</td><td align="left" rowspan="1" colspan="1">3,192</td><td align="left" rowspan="1" colspan="1">3,78E-02</td><td align="left" rowspan="1" colspan="1">5,116</td><td align="left" rowspan="1" colspan="1">2,21E-02</td><td align="left" rowspan="1" colspan="1">sgn-U215735</td><td align="left" rowspan="1" colspan="1">heavy-metal-associated domain-containing protein</td></tr><tr><td align="left" rowspan="1" colspan="1">10,576</td><td align="left" rowspan="1" colspan="1">2,96E-02</td><td align="left" rowspan="1" colspan="1">10,126</td><td align="left" rowspan="1" colspan="1">1,65E-02</td><td align="left" rowspan="1" colspan="1">7,662</td><td align="left" rowspan="1" colspan="1">1,66E-03</td><td align="left" rowspan="1" colspan="1">sgn-U226639</td><td align="left" rowspan="1" colspan="1">cysteine protease XBCP3</td></tr><tr><td align="left" rowspan="1" colspan="1">10,292</td><td align="left" rowspan="1" colspan="1">2,30E-03</td><td align="left" rowspan="1" colspan="1">3,273</td><td align="left" rowspan="1" colspan="1">2,51E-02</td><td align="left" rowspan="1" colspan="1">9,745</td><td align="left" rowspan="1" colspan="1">1,96E-02</td><td align="left" rowspan="1" colspan="1">sgn-U222678</td><td align="left" rowspan="1" colspan="1">ABC transporter family protein similar to ABC1 protein</td></tr><tr><td align="left" rowspan="1" colspan="1">10,171</td><td align="left" rowspan="1" colspan="1">3,45E-02</td><td align="left" rowspan="1" colspan="1">3,181</td><td align="left" rowspan="1" colspan="1">8,48E-03</td><td align="left" rowspan="1" colspan="1">3,391</td><td align="left" rowspan="1" colspan="1">2,79E-02</td><td align="left" rowspan="1" colspan="1">sgn-U220461</td><td align="left" rowspan="1" colspan="1">unknown protein</td></tr><tr><td align="left" rowspan="1" colspan="1">9,650</td><td align="left" rowspan="1" colspan="1">1,78E-02</td><td align="left" rowspan="1" colspan="1">3,970</td><td align="left" rowspan="1" colspan="1">3,15E-03</td><td align="left" rowspan="1" colspan="1">6,221</td><td align="left" rowspan="1" colspan="1">2,55E-06</td><td align="left" rowspan="1" colspan="1">sgn-U213637</td><td align="left" rowspan="1" colspan="1">WRKY family transcription factor DNA-binding protein 4 WRKY4</td></tr><tr><td align="left" rowspan="1" colspan="1">9,611</td><td align="left" rowspan="1" colspan="1">1,51E-03</td><td align="left" rowspan="1" colspan="1">4,667</td><td align="left" rowspan="1" colspan="1">2,63E-02</td><td align="left" rowspan="1" colspan="1">3,541</td><td align="left" rowspan="1" colspan="1">1,27E-02</td><td align="left" rowspan="1" colspan="1">sgn-U217986</td><td align="left" rowspan="1" colspan="1">patatin-like latex allergen [Hevea brasiliensis]</td></tr><tr><td align="left" rowspan="1" colspan="1">9,092</td><td align="left" rowspan="1" colspan="1">1,10E-03</td><td align="left" rowspan="1" colspan="1">11,458</td><td align="left" rowspan="1" colspan="1">3,00E-04</td><td align="left" rowspan="1" colspan="1">3,311</td><td align="left" rowspan="1" colspan="1">6,37E-03</td><td align="left" rowspan="1" colspan="1">sgn-U220022</td><td align="left" rowspan="1" colspan="1">CLH2_ARATH Chlorophyllase 2</td></tr><tr><td align="left" rowspan="1" colspan="1">8,844</td><td align="left" rowspan="1" colspan="1">1,46E-02</td><td align="left" rowspan="1" colspan="1">6,183</td><td align="left" rowspan="1" colspan="1">9,89E-03</td><td align="left" rowspan="1" colspan="1">3,503</td><td align="left" rowspan="1" colspan="1">8,60E-03</td><td align="left" rowspan="1" colspan="1">sgn-U213926</td><td align="left" rowspan="1" colspan="1">drought-induced protein SDi-6 -</td></tr><tr><td align="left" rowspan="1" colspan="1">8,427</td><td align="left" rowspan="1" colspan="1">1,26E-02</td><td align="left" rowspan="1" colspan="1">6,856</td><td align="left" rowspan="1" colspan="1">6,09E-04</td><td align="left" rowspan="1" colspan="1">5,731</td><td align="left" rowspan="1" colspan="1">3,95E-03</td><td align="left" rowspan="1" colspan="1">sgn-U212562</td><td align="left" rowspan="1" colspan="1">glutamate decarboxylase 1 (GAD 1)</td></tr><tr><td align="left" rowspan="1" colspan="1">8,292</td><td align="left" rowspan="1" colspan="1">3,03E-02</td><td align="left" rowspan="1" colspan="1">4,683</td><td align="left" rowspan="1" colspan="1">1,71E-02</td><td align="left" rowspan="1" colspan="1">3,999</td><td align="left" rowspan="1" colspan="1">1,13E-02</td><td align="left" rowspan="1" colspan="1">sgn-U225595</td><td align="left" rowspan="1" colspan="1">AUX/IAA family</td></tr><tr><td align="left" rowspan="1" colspan="1">7,781</td><td align="left" rowspan="1" colspan="1">3,55E-03</td><td align="left" rowspan="1" colspan="1">3,635</td><td align="left" rowspan="1" colspan="1">2,40E-03</td><td align="left" rowspan="1" colspan="1">3,167</td><td align="left" rowspan="1" colspan="1">1,91E-04</td><td align="left" rowspan="1" colspan="1">sgn-U213578</td><td align="left" rowspan="1" colspan="1">BEL1-related homeotic protein 11 [Solanum tuberosum]</td></tr><tr><td align="left" rowspan="1" colspan="1">7,479</td><td align="left" rowspan="1" colspan="1">1,53E-02</td><td align="left" rowspan="1" colspan="1">6,450</td><td align="left" rowspan="1" colspan="1">2,52E-02</td><td align="left" rowspan="1" colspan="1">8,916</td><td align="left" rowspan="1" colspan="1">4,75E-03</td><td align="left" rowspan="1" colspan="1">sgn-U222728</td><td align="left" rowspan="1" colspan="1">senescence-associated protein -related</td></tr><tr><td align="left" rowspan="1" colspan="1">7,392</td><td align="left" rowspan="1" colspan="1">1,01E-03</td><td align="left" rowspan="1" colspan="1">3,554</td><td align="left" rowspan="1" colspan="1">6,16E-03</td><td align="left" rowspan="1" colspan="1">3,866</td><td align="left" rowspan="1" colspan="1">6,33E-03</td><td align="left" rowspan="1" colspan="1">sgn-U216076</td><td align="left" rowspan="1" colspan="1">receptor serine/threonine kinase PR5K [Arabidopsis thaliana]</td></tr><tr><td align="left" rowspan="1" colspan="1">7,323</td><td align="left" rowspan="1" colspan="1">1,37E-02</td><td align="left" rowspan="1" colspan="1">3,318</td><td align="left" rowspan="1" colspan="1">2,05E-02</td><td align="left" rowspan="1" colspan="1">6,838</td><td align="left" rowspan="1" colspan="1">9,15E-05</td><td align="left" rowspan="1" colspan="1">sgn-U224875</td><td align="left" rowspan="1" colspan="1">heat shock protein family</td></tr><tr><td align="left" rowspan="1" colspan="1">6,728</td><td align="left" rowspan="1" colspan="1">2,91E-02</td><td align="left" rowspan="1" colspan="1">5,074</td><td align="left" rowspan="1" colspan="1">6,38E-03</td><td align="left" rowspan="1" colspan="1">3,774</td><td align="left" rowspan="1" colspan="1">2,95E-03</td><td align="left" rowspan="1" colspan="1">sgn-U216827</td><td align="left" rowspan="1" colspan="1">cysteine proteinase</td></tr><tr><td align="left" rowspan="1" colspan="1">6,554</td><td align="left" rowspan="1" colspan="1">2,01E-03</td><td align="left" rowspan="1" colspan="1">8,854</td><td align="left" rowspan="1" colspan="1">3,68E-03</td><td align="left" rowspan="1" colspan="1">3,714</td><td align="left" rowspan="1" colspan="1">1,33E-02</td><td align="left" rowspan="1" colspan="1">sgn-U216414</td><td align="left" rowspan="1" colspan="1">unknown protein</td></tr><tr><td align="left" rowspan="1" colspan="1">5,751</td><td align="left" rowspan="1" colspan="1">4,52E-02</td><td align="left" rowspan="1" colspan="1">3,631</td><td align="left" rowspan="1" colspan="1">7,63E-03</td><td align="left" rowspan="1" colspan="1">3,077</td><td align="left" rowspan="1" colspan="1">1,24E-03</td><td align="left" rowspan="1" colspan="1">sgn-U213519</td><td align="left" rowspan="1" colspan="1">dehydration-induced protein (ERD15)</td></tr><tr><td align="left" rowspan="1" colspan="1">5,290</td><td align="left" rowspan="1" colspan="1">4,44E-03</td><td align="left" rowspan="1" colspan="1">7,321</td><td align="left" rowspan="1" colspan="1">2,52E-03</td><td align="left" rowspan="1" colspan="1">3,511</td><td align="left" rowspan="1" colspan="1">2,60E-04</td><td align="left" rowspan="1" colspan="1">sgn-U218536</td><td align="left" rowspan="1" colspan="1">serine/threonine protein kinase</td></tr><tr><td align="left" rowspan="1" colspan="1">5,208</td><td align="left" rowspan="1" colspan="1">4,04E-02</td><td align="left" rowspan="1" colspan="1">3,078</td><td align="left" rowspan="1" colspan="1">6,31E-03</td><td align="left" rowspan="1" colspan="1">3,248</td><td align="left" rowspan="1" colspan="1">2,52E-02</td><td align="left" rowspan="1" colspan="1">sgn-U232570</td><td align="left" rowspan="1" colspan="1">unknown protein</td></tr><tr><td align="left" rowspan="1" colspan="1">5,179</td><td align="left" rowspan="1" colspan="1">6,48E-03</td><td align="left" rowspan="1" colspan="1">3,203</td><td align="left" rowspan="1" colspan="1">4,97E-03</td><td align="left" rowspan="1" colspan="1">4,092</td><td align="left" rowspan="1" colspan="1">4,23E-03</td><td align="left" rowspan="1" colspan="1">sgn-U212706</td><td align="left" rowspan="1" colspan="1">light regulated protein -related</td></tr><tr><td align="left" rowspan="1" colspan="1">4,860</td><td align="left" rowspan="1" colspan="1">8,15E-03</td><td align="left" rowspan="1" colspan="1">4,314</td><td align="left" rowspan="1" colspan="1">1,92E-02</td><td align="left" rowspan="1" colspan="1">4,049</td><td align="left" rowspan="1" colspan="1">1,99E-03</td><td align="left" rowspan="1" colspan="1">sgn-U219226</td><td align="left" rowspan="1" colspan="1">senescence-associated protein</td></tr><tr><td align="left" rowspan="1" colspan="1">3,926</td><td align="left" rowspan="1" colspan="1">3,33E-02</td><td align="left" rowspan="1" colspan="1">4,668</td><td align="left" rowspan="1" colspan="1">3,77E-03</td><td align="left" rowspan="1" colspan="1">3,059</td><td align="left" rowspan="1" colspan="1">9,41E-03</td><td align="left" rowspan="1" colspan="1">sgn-U229252</td><td align="left" rowspan="1" colspan="1">nitrate transporter NRT1-1 [Glycine max]</td></tr><tr><td align="left" rowspan="1" colspan="1">3,864</td><td align="left" rowspan="1" colspan="1">1,81E-02</td><td align="left" rowspan="1" colspan="1">3,788</td><td align="left" rowspan="1" colspan="1">2,82E-02</td><td align="left" rowspan="1" colspan="1">3,152</td><td align="left" rowspan="1" colspan="1">2,06E-02</td><td align="left" rowspan="1" colspan="1">sgn-U214216</td><td align="left" rowspan="1" colspan="1">putative steroid membrane binding protein [Oryza sativa (japonica cultivar-group)]</td></tr><tr><td align="left" rowspan="1" colspan="1">3,697</td><td align="left" rowspan="1" colspan="1">9,00E-03</td><td align="left" rowspan="1" colspan="1">5,706</td><td align="left" rowspan="1" colspan="1">1,16E-02</td><td align="left" rowspan="1" colspan="1">10,814</td><td align="left" rowspan="1" colspan="1">9,84E-04</td><td align="left" rowspan="1" colspan="1">sgn-U215359</td><td align="left" rowspan="1" colspan="1">1-aminocyclopropane-1-carboxylate oxidase homolog (Protein E8)</td></tr><tr><td align="left" rowspan="1" colspan="1">3,071</td><td align="left" rowspan="1" colspan="1">1,25E-02</td><td align="left" rowspan="1" colspan="1">3,300</td><td align="left" rowspan="1" colspan="1">2,47E-02</td><td align="left" rowspan="1" colspan="1">8,554</td><td align="left" rowspan="1" colspan="1">7,00E-03</td><td align="left" rowspan="1" colspan="1">sgn-U217504</td><td align="left" rowspan="1" colspan="1">zinc finger (C3HC4-type RING finger) family protein (RMA1)</td></tr></tbody></table><table frame="hsides" rules="groups"><colgroup span="1"><col align="left" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/></colgroup><thead><tr><td colspan="2" align="left" rowspan="1">zt8/zt0</td><td colspan="2" align="left" rowspan="1">zt16/zt0</td><td colspan="2" align="left" rowspan="1">zt20/zt0</td><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">B</td></tr><tr><td align="left" rowspan="1" colspan="1">Ratio</td><td align="left" rowspan="1" colspan="1">P-value</td><td align="left" rowspan="1" colspan="1">zt16</td><td align="left" rowspan="1" colspan="1">P-value</td><td align="left" rowspan="1" colspan="1">zt20</td><td align="left" rowspan="1" colspan="1">P-value</td><td align="left" rowspan="1" colspan="1">Gene_ID</td><td align="left" rowspan="1" colspan="1">Annotation</td></tr></thead><tbody><tr><td align="left" rowspan="1" colspan="1">0,004</td><td align="left" rowspan="1" colspan="1">1,16E-04</td><td align="left" rowspan="1" colspan="1">0,046</td><td align="left" rowspan="1" colspan="1">8,98E-05</td><td align="left" rowspan="1" colspan="1">0,193</td><td align="left" rowspan="1" colspan="1">5,89E-04</td><td align="left" rowspan="1" colspan="1">sgn-U218628</td><td align="left" rowspan="1" colspan="1">protein kinase family</td></tr><tr><td align="left" rowspan="1" colspan="1">0,001</td><td align="left" rowspan="1" colspan="1">7,75E-05</td><td align="left" rowspan="1" colspan="1">0,072</td><td align="left" rowspan="1" colspan="1">2,18E-03</td><td align="left" rowspan="1" colspan="1">0,251</td><td align="left" rowspan="1" colspan="1">4,60E-06</td><td align="left" rowspan="1" colspan="1">sgn-U227452</td><td align="left" rowspan="1" colspan="1">heat stress transcription factor HSFA9 [Helianthus annuus]</td></tr><tr><td align="left" rowspan="1" colspan="1">0,009</td><td align="left" rowspan="1" colspan="1">7,76E-05</td><td align="left" rowspan="1" colspan="1">0,119</td><td align="left" rowspan="1" colspan="1">1,17E-04</td><td align="left" rowspan="1" colspan="1">0,155</td><td align="left" rowspan="1" colspan="1">2,10E-04</td><td align="left" rowspan="1" colspan="1">sgn-U220755</td><td align="left" rowspan="1" colspan="1">protease inhibitor/seed storage/lipid transfer protein (LTP) family</td></tr><tr><td align="left" rowspan="1" colspan="1">0,026</td><td align="left" rowspan="1" colspan="1">3,27E-07</td><td align="left" rowspan="1" colspan="1">0,009</td><td align="left" rowspan="1" colspan="1">8,04E-09</td><td align="left" rowspan="1" colspan="1">0,021</td><td align="left" rowspan="1" colspan="1">2,25E-08</td><td align="left" rowspan="1" colspan="1">BT012912</td><td align="left" rowspan="1" colspan="1">Lycopersicon esculentum clone 114030R</td></tr><tr><td align="left" rowspan="1" colspan="1">0,044</td><td align="left" rowspan="1" colspan="1">4,58E-04</td><td align="left" rowspan="1" colspan="1">0,094</td><td align="left" rowspan="1" colspan="1">1,45E-02</td><td align="left" rowspan="1" colspan="1">0,018</td><td align="left" rowspan="1" colspan="1">1,52E-05</td><td align="left" rowspan="1" colspan="1">sgn-U237511</td><td align="left" rowspan="1" colspan="1">LHY</td></tr><tr><td align="left" rowspan="1" colspan="1">0,048</td><td align="left" rowspan="1" colspan="1">2,51E-05</td><td align="left" rowspan="1" colspan="1">0,067</td><td align="left" rowspan="1" colspan="1">2,03E-06</td><td align="left" rowspan="1" colspan="1">0,181</td><td align="left" rowspan="1" colspan="1">1,05E-05</td><td align="left" rowspan="1" colspan="1">sgn-U227174</td><td align="left" rowspan="1" colspan="1">pectinesterase family</td></tr><tr><td align="left" rowspan="1" colspan="1">0,054</td><td align="left" rowspan="1" colspan="1">9,85E-07</td><td align="left" rowspan="1" colspan="1">0,142</td><td align="left" rowspan="1" colspan="1">8,72E-06</td><td align="left" rowspan="1" colspan="1">0,299</td><td align="left" rowspan="1" colspan="1">1,38E-04</td><td align="left" rowspan="1" colspan="1">sgn-U226241</td><td align="left" rowspan="1" colspan="1">L-aspartate oxidase family</td></tr><tr><td align="left" rowspan="1" colspan="1">0,082</td><td align="left" rowspan="1" colspan="1">1,55E-07</td><td align="left" rowspan="1" colspan="1">0,099</td><td align="left" rowspan="1" colspan="1">3,47E-06</td><td align="left" rowspan="1" colspan="1">0,010</td><td align="left" rowspan="1" colspan="1">9,01E-08</td><td align="left" rowspan="1" colspan="1">sgn-U217418</td><td align="left" rowspan="1" colspan="1">heat shock protein 81-1 (HSP81-1/heat shock protein 83/HSP83)</td></tr><tr><td align="left" rowspan="1" colspan="1">0,085</td><td align="left" rowspan="1" colspan="1">2,63E-06</td><td align="left" rowspan="1" colspan="1">0,005</td><td align="left" rowspan="1" colspan="1">2,64E-06</td><td align="left" rowspan="1" colspan="1">0,248</td><td align="left" rowspan="1" colspan="1">1,85E-04</td><td align="left" rowspan="1" colspan="1">sgn-U213666</td><td align="left" rowspan="1" colspan="1">expressed protein</td></tr><tr><td align="left" rowspan="1" colspan="1">0,085</td><td align="left" rowspan="1" colspan="1">3,42E-06</td><td align="left" rowspan="1" colspan="1">0,272</td><td align="left" rowspan="1" colspan="1">1,45E-03</td><td align="left" rowspan="1" colspan="1">0,152</td><td align="left" rowspan="1" colspan="1">2,15E-05</td><td align="left" rowspan="1" colspan="1">sgn-U218682</td><td align="left" rowspan="1" colspan="1">abnormal floral organs protein (AFO) (FIL) YABBY1 (YABBY1)</td></tr><tr><td align="left" rowspan="1" colspan="1">0,087</td><td align="left" rowspan="1" colspan="1">2,09E-06</td><td align="left" rowspan="1" colspan="1">0,087</td><td align="left" rowspan="1" colspan="1">4,49E-06</td><td align="left" rowspan="1" colspan="1">0,286</td><td align="left" rowspan="1" colspan="1">3,94E-04</td><td align="left" rowspan="1" colspan="1">sgn-U215829</td><td align="left" rowspan="1" colspan="1">beta-carotene hydroxylase</td></tr><tr><td align="left" rowspan="1" colspan="1">0,090</td><td align="left" rowspan="1" colspan="1">2,96E-06</td><td align="left" rowspan="1" colspan="1">0,177</td><td align="left" rowspan="1" colspan="1">7,57E-05</td><td align="left" rowspan="1" colspan="1">0,109</td><td align="left" rowspan="1" colspan="1">5,24E-06</td><td align="left" rowspan="1" colspan="1">sgn-U222349</td><td align="left" rowspan="1" colspan="1">expressed protein</td></tr><tr><td align="left" rowspan="1" colspan="1">0,097</td><td align="left" rowspan="1" colspan="1">1,03E-06</td><td align="left" rowspan="1" colspan="1">0,163</td><td align="left" rowspan="1" colspan="1">1,16E-05</td><td align="left" rowspan="1" colspan="1">0,256</td><td align="left" rowspan="1" colspan="1">3,41E-05</td><td align="left" rowspan="1" colspan="1">sgn-U222883</td><td align="left" rowspan="1" colspan="1">phytochrome A supressor spa1</td></tr><tr><td align="left" rowspan="1" colspan="1">0,104</td><td align="left" rowspan="1" colspan="1">4,28E-05</td><td align="left" rowspan="1" colspan="1">0,177</td><td align="left" rowspan="1" colspan="1">4,42E-05</td><td align="left" rowspan="1" colspan="1">0,083</td><td align="left" rowspan="1" colspan="1">1,70E-05</td><td align="left" rowspan="1" colspan="1">sgn-U216986</td><td align="left" rowspan="1" colspan="1">DnaJ protein AtJ3</td></tr><tr><td align="left" rowspan="1" colspan="1">0,130</td><td align="left" rowspan="1" colspan="1">1,02E-04</td><td align="left" rowspan="1" colspan="1">0,035</td><td align="left" rowspan="1" colspan="1">1,74E-05</td><td align="left" rowspan="1" colspan="1">0,071</td><td align="left" rowspan="1" colspan="1">1,76E-07</td><td align="left" rowspan="1" colspan="1">sgn-U212621</td><td align="left" rowspan="1" colspan="1">AMP-dependent synthetase and ligase family</td></tr><tr><td align="left" rowspan="1" colspan="1">0,140</td><td align="left" rowspan="1" colspan="1">1,13E-04</td><td align="left" rowspan="1" colspan="1">0,041</td><td align="left" rowspan="1" colspan="1">8,67E-03</td><td align="left" rowspan="1" colspan="1">0,169</td><td align="left" rowspan="1" colspan="1">6,02E-06</td><td align="left" rowspan="1" colspan="1">sgn-U214765</td><td align="left" rowspan="1" colspan="1">copper-binding protein family</td></tr><tr><td align="left" rowspan="1" colspan="1">0,188</td><td align="left" rowspan="1" colspan="1">2,52E-06</td><td align="left" rowspan="1" colspan="1">0,128</td><td align="left" rowspan="1" colspan="1">6,98E-06</td><td align="left" rowspan="1" colspan="1">0,132</td><td align="left" rowspan="1" colspan="1">1,31E-05</td><td align="left" rowspan="1" colspan="1">sgn-U222990</td><td align="left" rowspan="1" colspan="1">expressed protein</td></tr><tr><td align="left" rowspan="1" colspan="1">0,198</td><td align="left" rowspan="1" colspan="1">2,22E-04</td><td align="left" rowspan="1" colspan="1">0,182</td><td align="left" rowspan="1" colspan="1">2,45E-04</td><td align="left" rowspan="1" colspan="1">0,278</td><td align="left" rowspan="1" colspan="1">8,19E-04</td><td align="left" rowspan="1" colspan="1">sgn-U217455</td><td align="left" rowspan="1" colspan="1">linker histone protein -related</td></tr><tr><td align="left" rowspan="1" colspan="1">0,215</td><td align="left" rowspan="1" colspan="1">3,86E-04</td><td align="left" rowspan="1" colspan="1">0,128</td><td align="left" rowspan="1" colspan="1">1,22E-06</td><td align="left" rowspan="1" colspan="1">0,144</td><td align="left" rowspan="1" colspan="1">2,38E-05</td><td align="left" rowspan="1" colspan="1">AY547273</td><td align="left" rowspan="1" colspan="1">early light inducible protein (ELIP)</td></tr><tr><td align="left" rowspan="1" colspan="1">0,216</td><td align="left" rowspan="1" colspan="1">8,70E-06</td><td align="left" rowspan="1" colspan="1">0,223</td><td align="left" rowspan="1" colspan="1">7,48E-05</td><td align="left" rowspan="1" colspan="1">0,285</td><td align="left" rowspan="1" colspan="1">2,17E-06</td><td align="left" rowspan="1" colspan="1">sgn-U222868</td><td align="left" rowspan="1" colspan="1">phytochelatin synthetase</td></tr><tr><td align="left" rowspan="1" colspan="1">0,218</td><td align="left" rowspan="1" colspan="1">1,29E-05</td><td align="left" rowspan="1" colspan="1">0,029</td><td align="left" rowspan="1" colspan="1">1,13E-02</td><td align="left" rowspan="1" colspan="1">0,091</td><td align="left" rowspan="1" colspan="1">4,92E-05</td><td align="left" rowspan="1" colspan="1">sgn-U224520</td><td align="left" rowspan="1" colspan="1">FtsH protease</td></tr><tr><td align="left" rowspan="1" colspan="1">0,240</td><td align="left" rowspan="1" colspan="1">1,21E-04</td><td align="left" rowspan="1" colspan="1">0,111</td><td align="left" rowspan="1" colspan="1">3,91E-02</td><td align="left" rowspan="1" colspan="1">0,244</td><td align="left" rowspan="1" colspan="1">2,01E-04</td><td align="left" rowspan="1" colspan="1">Control-4</td><td align="left" rowspan="1" colspan="1">2-oxoglutarate-dependent dioxygenase</td></tr><tr><td align="left" rowspan="1" colspan="1">0,249</td><td align="left" rowspan="1" colspan="1">2,15E-05</td><td align="left" rowspan="1" colspan="1">0,325</td><td align="left" rowspan="1" colspan="1">1,00E-04</td><td align="left" rowspan="1" colspan="1">0,323</td><td align="left" rowspan="1" colspan="1">6,23E-05</td><td align="left" rowspan="1" colspan="1">sgn-U215717</td><td align="left" rowspan="1" colspan="1">omega-6 fatty acid desaturase</td></tr><tr><td align="left" rowspan="1" colspan="1">0,267</td><td align="left" rowspan="1" colspan="1">2,02E-03</td><td align="left" rowspan="1" colspan="1">0,060</td><td align="left" rowspan="1" colspan="1">1,72E-03</td><td align="left" rowspan="1" colspan="1">0,021</td><td align="left" rowspan="1" colspan="1">6,00E-05</td><td align="left" rowspan="1" colspan="1">sgn-U218114</td><td align="left" rowspan="1" colspan="1">myb family transcription factor [Arabidopsis thaliana]</td></tr><tr><td align="left" rowspan="1" colspan="1">0,270</td><td align="left" rowspan="1" colspan="1">2,49E-04</td><td align="left" rowspan="1" colspan="1">0,130</td><td align="left" rowspan="1" colspan="1">5,06E-06</td><td align="left" rowspan="1" colspan="1">0,132</td><td align="left" rowspan="1" colspan="1">1,18E-05</td><td align="left" rowspan="1" colspan="1">sgn-U221454</td><td align="left" rowspan="1" colspan="1">ammonium transporter</td></tr><tr><td align="left" rowspan="1" colspan="1">0,320</td><td align="left" rowspan="1" colspan="1">2,19E-06</td><td align="left" rowspan="1" colspan="1">0,145</td><td align="left" rowspan="1" colspan="1">1,02E-06</td><td align="left" rowspan="1" colspan="1">0,097</td><td align="left" rowspan="1" colspan="1">2,92E-06</td><td align="left" rowspan="1" colspan="1">sgn-U217678</td><td align="left" rowspan="1" colspan="1">delta 9 acyl-lipid desaturase (ADS1)</td></tr><tr><td align="left" rowspan="1" colspan="1">0,324</td><td align="left" rowspan="1" colspan="1">1,92E-03</td><td align="left" rowspan="1" colspan="1">0,137</td><td align="left" rowspan="1" colspan="1">2,81E-03</td><td align="left" rowspan="1" colspan="1">0,006</td><td align="left" rowspan="1" colspan="1">3,73E-07</td><td align="left" rowspan="1" colspan="1">sgn-U225455</td><td align="left" rowspan="1" colspan="1">CONSTANS-LIKE 1 (COL1)</td></tr></tbody></table></alternatives></table-wrap>
<supplementary-material content-type="local-data" id="pone.0002798.s001"><label>Data S1</label><caption><p>(3.23 MB TXT)</p></caption></supplementary-material>
<fn-group><fn fn-type="COI-statement"><p><bold>Competing Interests: </bold>The authors have declared that no competing interests exist.</p></fn><fn fn-type="financial-disclosure"><p><bold>Funding: </bold>Work supported by the Italian Ministry of Research Project GEPROT DM1570, and by the EU (EUSOL project).</p></fn></fn-group>
<graphic xlink:href="pone.0002798.g001"/> <graphic xlink:href="pone.0002798.g002"/> <graphic id="pone-0002798-t001-1" xlink:href="pone.0002798.t001"/> <graphic xlink:href="pone.0002798.g003"/> <graphic xlink:href="pone.0002798.g004"/> <graphic xlink:href="pone.0002798.g005"/> <graphic xlink:href="pone.0002798.g006"/> <graphic xlink:href="pone.0002798.g007"/> <graphic xlink:href="pone.0002798.g008"/> <graphic xlink:href="pone.0002798.g009"/>
<media xlink:href="pone.0002798.s001.txt"><caption><p>Click here for additional data file.</p></caption></media>
[{"label": ["5"], "element-citation": ["\n"], "surname": ["Folta", "Maruhnich"], "given-names": ["KM", "SA"], "year": ["2007"], "article-title": ["Green light: a signal to slow down or stop."], "source": ["J Exp Bot"], "volume": ["146"], "fpage": ["227"], "lpage": ["228"]}, {"label": ["31"], "element-citation": ["\n"], "surname": ["Locke", "Southern", "Kozma-Bognar", "Hibberd", "Brown"], "given-names": ["JCW", "MM", "L", "V", "PE"], "year": ["2005b"], "article-title": ["Extension of a genetic network model by iterative experimentation and mathematical analysis."], "source": ["Mol Syst Biol"], "volume": ["1"], "fpage": ["13"]}, {"label": ["33"], "element-citation": ["\n"], "surname": ["Locke", "Kozma-Bognar", "Gould", "Feher", "Kevei"], "given-names": ["JC", "L", "PD", "B", "E"], "year": ["2006"], "article-title": ["Experimental validation of a predicted feedback loop in the multi-oscillator clock of "], "italic": ["Arabidopsis thaliana"], "source": ["Mol Syst Biol"], "volume": ["2"], "fpage": ["1"], "lpage": ["6"]}, {"label": ["41"], "element-citation": ["\n"], "surname": ["Perrotta", "Yahoubian", "Nebuloso", "Renzi", "Giuliano"], "given-names": ["G", "G", "E", "L", "G"], "year": ["2001"], "article-title": ["Tomato and barley contain duplicated copies of cryptochrome 1."], "source": ["Plant Cell Environ"], "volume": ["24"], "fpage": ["991"], "lpage": ["997"]}, {"label": ["59"], "element-citation": ["\n"], "surname": ["Hotta", "Gardner", "Hubbard", "Baek", "Dalchau"], "given-names": ["CT", "MJ", "KE", "SJ", "N"], "year": ["2007"], "article-title": ["Modulation of environmental responses of plants by circadian clocks."], "source": ["Plant Cell and Environment"], "volume": ["30"], "fpage": ["333"], "lpage": ["349"]}, {"label": ["64"], "element-citation": ["\n"], "surname": ["Millar"], "given-names": ["AJ"], "year": ["1999"], "article-title": ["Biological clocks in "], "italic": ["Arabidopsis thaliana"], "source": ["New Phytol"], "volume": ["141"], "fpage": ["175"], "lpage": ["97"]}, {"label": ["66"], "element-citation": ["\n"], "surname": ["Sambrook", "Fritsch", "Maniatis"], "given-names": ["J", "EF", "T"], "year": ["1989"], "source": ["Molecular cloning: A Laboratory Manual, 2nd edn"], "publisher-loc": ["Cold Spring Harbor"], "publisher-name": ["Cold Spring Harbor Laboratory Press"]}, {"label": ["67"], "element-citation": ["\n"], "surname": ["Cleveland", "Devlin"], "given-names": ["WS", "SJ"], "year": ["1998"], "article-title": ["Locally weighted regression: an approach to regression analysis by local fitting."], "source": ["J of Am Stat Assoc"], "volume": ["83"], "fpage": ["596"], "lpage": ["610"]}]
PLoS One. 2008 Jul 30; 3(7):e2798